Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2015 Feb 15;36(5):348-54.
doi: 10.1002/jcc.23804. Epub 2014 Dec 8.

pmx: Automated protein structure and topology generation for alchemical perturbations

Vytautas Gapsys  1 , Servaas Michielssens, Daniel Seeliger, Bert L de Groot
Affiliations

pmx: Automated protein structure and topology generation for alchemical perturbations

Vytautas Gapsys et al. J Comput Chem. .

Abstract

Computational protein design requires methods to accurately estimate free energy changes in protein stability or binding upon an amino acid mutation. From the different approaches available, molecular dynamics-based alchemical free energy calculations are unique in their accuracy and solid theoretical basis. The challenge in using these methods lies in the need to generate hybrid structures and topologies representing two physical states of a system. A custom made hybrid topology may prove useful for a particular mutation of interest, however, a high throughput mutation analysis calls for a more general approach. In this work, we present an automated procedure to generate hybrid structures and topologies for the amino acid mutations in all commonly used force fields. The described software is compatible with the Gromacs simulation package. The mutation libraries are readily supported for five force fields, namely Amber99SB, Amber99SB*-ILDN, OPLS-AA/L, Charmm22*, and Charmm36.

Keywords: alchemy; free energy calculations; molecular dynamics; mutations; thermostability.

PubMed Disclaimer

Figures

Figure 1
Figure 1
An overall workflow of the pmx hybrid topology and structure generation. The force field specific mutation libraries are created by generate_hybrid_residue.py using the amino acid topologies, bonded, and nonbonded parameters as defined in the force field of interest. By default pmx provides mutation libraries for five commonly used force fields. A hybrid mutated structure is generated by mutate.py. The Gromacs standard tool pdb2gmx uses the hybrid structure to create its topology, which subsequently is processed by generate_hybrid_topology.py by adding the required parameters for the A and B states.
Figure 2
Figure 2
A schematic illustration of the two main pmx data structures used in the hybrid structure and topology generation. The top row describes the classes containing protein structure information, whereas the bottom rows describe the topological information. A detailed description can be found in the main text.
Figure 3
Figure 3
Thermodynamic cycles used for validation. A) The double system in a single simulation box setup. Two tripeptides with the hybrid structures, V2F and F2V, are placed in the same box. During an alchemical transition, valine is morphed into phenylalanine and phenylalanine into valine. The overall state of the system remains unchanged, therefore, the calculated ΔΔG must be equal to 0 kJ/mol. B: results of the ΔΔG calculation over a closed thermodynamic cycle for nine amino acid mutation pairs for five force fields.
Figure 4
Figure 4
Illustration of the input/output structures and an output topology of the pmx scripts. A) A capped Gly-Val-Gly tripeptide as a starting structure for the mutation. B) Structure of the tripeptide with the hybrid V2F amino acid. Valine (dark gray) is in the physical state A, whereas phenylalanine (light gray) is in the state B. C) An excerpt from the hybrid V2F topology.
See this image and copyright information in PMC

Similar articles

See all similar articles

Cited by

See all "Cited by" articles

References

    1. Pantazes R, Grisewood M, Maranas C. Curr. Opin. Struct. Biol. 2011;21:467. - PubMed
    1. Schreiber G, Fleishman S. Curr. Opin. Struct. Biol. 2013;23:903. - PubMed
    1. Khan S, Vihinen M. Hum. Mutat. 2010;31:675. - PubMed
    1. Capriotti E, Fariselli P, Casadio R. Nucl. Acids Res. 2005;33:W306. - PMC - PubMed
    1. Capriotti E, Fariselli P, Rossi I, Casadio R. BMC Bioinf. 2008;9:S6. - PMC - PubMed

LinkOut - more resources