Do PDB bruto a uma trajetória analisável — sete etapas, cada uma com um widget pra arrastar, girar e avançar enquanto lê. O sistema-exemplo é o mesmo do início ao fim: a cruzaína, principal cisteíno-protease de Trypanosoma cruzi, em complexo com o inibidor não-covalente B95 (PDB 3KKU2).
Uma estrutura cristalográfica como a 3KKU é uma fotografia: uma única conformação, obtida sob condições de cristalização distantes das fisiológicas, sem nenhuma informação sobre como a proteína realmente se move. A dinâmica molecular (MD) resolve isso numericamente — integra as equações de movimento de Newton para cada átomo do sistema segundo um campo de força que aproxima as interações físicas relevantes (ligações, ângulos, torções, van der Waals, eletrostática), permitindo observar como proteína, ligante e solvente se comportam ao longo do tempo.
Este material segue o pipeline padrão de uma simulação de complexo proteína–ligante no GROMACS1, da estrutura crua até uma trajetória analisável, usando como exemplo a cruzaína em complexo com o inibidor B95 (PDB 3KKU2).
Dois detalhes técnicos aparecem em praticamente toda etapa e valem uma explicação só: as ligações envolvendo hidrogênio são mantidas rígidas pelo algoritmo LINCS3, o que permite um passo de integração de 2 fs em vez de ~0,5 fs; e a eletrostática de longo alcance é tratada por Particle Mesh Ewald4, necessário porque a caixa é periódica e a lei de Coulomb não pode simplesmente ser truncada sem introduzir artefatos.
Cada seção abaixo tem um widget interativo — arraste, gire, avance os controles. Os números e as trajetórias específicos são esquemáticos (nada aqui é uma simulação real rodando no seu navegador), mas a física por trás de cada um foi escolhida pra ser consistente com o que você veria numa simulação de verdade desse sistema.
Antes de rodar pdb2gmx, alguém precisa decidir o que fica e o que sai — não é um
processo no tempo, é uma sequência de decisões: resolver ambiguidades da cristalografia,
separar o que é biológico do que é artefato de cristalização, e computar os hidrogênios que a
difração de raios-X normalmente não resolve. O GROMACS1
automatiza boa parte disso — as decisões de "o que fica" continuam sendo de quem monta o sistema.
Com a topologia pronta, o sistema ainda está no vácuo — sem solvente, sem contra-íons, sem
volume definido. O editconf estabelece uma caixa periódica ao redor do soluto, o
solvate a preenche com água (aqui, o modelo SPC/E5),
e o genion substitui algumas moléculas de água por íons — primeiro para neutralizar
a carga líquida do sistema, depois, opcionalmente, para aproximar a força iônica fisiológica (~0,15 M).
Colocar milhares de moléculas de água ao redor de uma proteína quase sempre gera alguns contatos artificialmente próximos — geometricamente aceitáveis para quem gerou a caixa, mas fisicamente catastróficos: forças de repulsão de centenas de milhares de kJ/mol/nm que, integradas diretamente, explodiriam a simulação nos primeiros passos. A minimização por steepest descent caminha na direção oposta ao gradiente da energia potencial, em passos discretos, até um mínimo local — não o mínimo global da paisagem conformacional, apenas o vale mais próximo de onde o sistema começou.
esquemático — paisagem de energia ilustrativa; ordens de grandeza reaisDepois da minimização, o sistema tem posições razoáveis mas nenhuma velocidade real — está, para todos os efeitos, a 0 K. O acoplamento a um termostato aquece o sistema até a temperatura alvo e mantém a distribuição de velocidades consistente com o ensemble canônico (NVT: partículas, volume e temperatura constantes). Aqui usamos o termostato V-rescale6, que reamostra a energia cinética estocasticamente e evita as patologias de esquemas mais simples, como o termostato de Berendsen.
estrutura real (PDB 3KKU) · dinâmica da água ilustrativa
Com a temperatura estável, falta o volume: a caixa que o editconf definiu foi um
chute razoável, não o volume de equilíbrio real do sistema solvatado naquela temperatura e
pressão. A equilibração NPT libera o volume e acopla um barostato — aqui,
C-rescale7 — que ajusta o tamanho da caixa até a
pressão média convergir para o valor alvo (tipicamente 1 bar), o que na prática significa a
densidade da água convergindo para ~1 g/cm³.
Estrutura real (bruta e processada), cartoon e sticks renderizados com 3Dmol.js. Água, caixa, íons, paisagem de energia e todos os números de temperatura/pressão/densidade/RMSD são esquemáticos — construídos pra ter a forma e a ordem de grandeza corretas, não uma trajetória real.