Guia Prático - Seção 3: Construção e Equilibração do Modelo de Membrana

3.1. Racional para o Modelo POPC:Colesterol (70:30)

A escolha de um modelo de membrana apropriado é um balanço entre a relevância biofísica e a tratabilidade computacional. O sistema POPC:Colesterol na proporção molar de 70:30 é um padrão-ouro para estudos de permeação, pois representa uma abstração justificada da complexa membrana intestinal.

Justificativa Biofísica

• POPC (1-palmitoil-2-oleoil-sn-glicero-3-fosfocolina): Como fosfolipídio PC, é o componente majoritário das membranas de mamíferos, e sua cadeia insaturada (oleoil) garante a fluidez necessária da membrana na temperatura fisiológica.
• Colesterol (30 mol%): Esta concentração induz a formação da fase líquido-ordenada (Lo). Esta fase é fisicamente distinta da fase líquido-desordenada (Ld) do POPC puro. A fase Lo é caracterizada por um alto ordenamento das cadeias lipídicas (similar a um gel) mas com alta mobilidade lateral (similar a um líquido). Este estado é mais denso, mais espesso e, crucialmente, **menos permeável** a pequenas moléculas, representando muito melhor a função de barreira da membrana plasmática real.

---

3.2. Construção da Bicamada com CHARMM-GUI

A construção manual de uma bicamada lipídica é um processo demorado e propenso a erros. A ferramenta web CHARMM-GUI Membrane Builder é a solução recomendada e padrão na área para automatizar e padronizar esta tarefa, gerando sistemas prontos para uso em diversos pacotes de DM, incluindo AMBER.

Protocolo Passo a Passo no CHARMM-GUI

1. Navegação: Vá para CHARMM-GUI e navegue até Input Generator → Membrane Builder → Bilayer Builder.
2. Tipo de Sistema: Selecione "Membrane Only System".
3. Componentes Lipídicos: Em "Heterogeneous Lipid", adicione os lípidos **POPC** e **CHOL** (colesterol). Para um sistema de aproximadamente 128 lípidos por folheto, especifique **90 moléculas de POPC e 38 de CHOL** tanto no folheto superior (Upperleaflet) quanto no inferior (Lowerleaflet), resultando na proporção ~70:30.
4. Hidratação e Íons: Defina uma hidratação de 40 águas por lípido e adicione íons (ex: KCL ou NACL) para uma concentração fisiológica de **0.15 M**. Escolha o método "Monte-Carlo" para posicionamento dos íons, que garante uma distribuição mais homogênea.
5. Geração para AMBER: Prossiga pelas etapas. No passo final ("Generate Equilibration and Dynamics Input"), **selecione AMBER** como o pacote de simulação alvo e certifique-se de que a opção de conversão de campo de força está ativa.
6. Download: Faça o download do arquivo .tgz completo gerado pelo projeto.

---

3.3. Protocolo de Equilibração Multifásico para AMBER

Dentro do arquivo baixado, no diretório amber/, você encontrará um README e uma série de scripts de input (`.in`). Estes arquivos compõem um protocolo de equilibração validado. Não é necessário criar scripts do zero; adapte e execute os que foram fornecidos. O princípio é a relaxação gradual para evitar artefatos.

Tabela Resumo do Protocolo de Equilibração Padrão

Estágio	Ensemble	Duração Típica	Propósito Principal	Restrições nos Lípidos
1. Minimização	-	~5000 passos	Remover contatos estéricos desfavoráveis.	Fortes
2. Aquecimento	NVT	200 ps	Aquecer gradualmente o sistema até 310 K em volume constante.	Fortes
3. Pressurização	NPT	~500 ps	Ajustar a pressão do sistema de forma controlada.	Fortes → Médias
4. Equilibração Curta	NPT	5 ns	Relaxar as restrições nos lípidos gradualmente.	Médias → Fracas
5. Equilibração Final	NPT	50-100 ns	Permitir a convergência total das propriedades coletivas da membrana.	Nenhuma

---

3.4. Validação do Equilíbrio: Métricas e Valores de Referência

A simulação só está equilibrada quando as propriedades macroscópicas do sistema param de exibir um desvio sistemático (*drift*) e flutuam em torno de um valor médio estável. Utilize o módulo cpptraj do AmberTools para analisar a trajetória da equilibração final (Estágio 5).

Critérios de Validação Essenciais

• Área por Lípido (APL): É a métrica mais importante. Para uma mistura POPC:Colesterol 70:30 a 310K, o valor esperado é de ~48-52 Ų. Se o valor da sua simulação estabilizar nesta faixa, é um forte indicativo de sucesso. (Comando cpptraj: areapermol).
• Espessura da Bicamada (d_P-P): A distância entre os picos de densidade dos átomos de fósforo. Para este sistema, um valor entre 38-42 Å é esperado. (Comando cpptraj: density).
• Parâmetros de Ordem (S_CD): Devem mostrar o platô característico do efeito de ordenamento do colesterol, com valores na região central da cadeia de ~0.40. (Comando cpptraj: lipidorder).

Conclusão da Fase

Apenas após a confirmação de que estas três propriedades atingiram um platô estável e estão em concordância com os valores de referência da literatura, o seu modelo de membrana pode ser considerado **fisicamente realista e pronto para receber a molécula CNFD**, iniciando os estudos de permeação.