On the Use of the Discrete Constant pH Molecular Dynamics to Describe the Conformational Space of Peptides
Open Access
- 29 December 2020
- Vol. 13 (1), 99
- https://doi.org/10.3390/polym13010099
Abstract
Solvent pH is an important property that defines the protonation state of the amino acids and, therefore, modulates the interactions and the conformational space of the biochemical systems. Generally, this thermodynamic variable is poorly considered in Molecular Dynamics (MD) simulations. Fortunately, this lack has been overcome by means of the Constant pH Molecular Dynamics (CPHMD) methods in the recent decades. Several studies have reported promising results from these approaches that include pH in simulations but focus on the prediction of the effective pKa of the amino acids. In this work, we want to shed some light on the CPHMD method and its implementation in the AMBER suitcase from a conformational point of view. To achieve this goal, we performed CPHMD and conventional MD (CMD) simulations of six protonatable amino acids in a blocked tripeptide structure to compare the conformational sampling and energy distributions of both methods. The results reveal strengths and weaknesses of the CPHMD method in the implementation of AMBER18 version. The change of the protonation state according to the chemical environment is presumably an improvement in the accuracy of the simulations. However, the simulations of the deprotonated forms are not consistent, which is related to an inaccurate assignment of the partial charges of the backbone atoms in the CPHMD residues. Therefore, we recommend the CPHMD methods of AMBER program but pointing out the need to compare structural properties with experimental data to bring reliability to the conformational sampling of the simulations.Keywords
Funding Information
- Generalitat de Catalunya (2017SGR1033)
This publication has 59 references indexed in Scilit:
- GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkitBioinformatics, 2013
- pH replica‐exchange method based on discrete protonation statesProteins-Structure Function and Bioinformatics, 2011
- Measuring the successes and deficiencies of constant pH molecular dynamics: A blind prediction studyProteins-Structure Function and Bioinformatics, 2011
- Constant pH Molecular Dynamics in Explicit Solvent with λ-DynamicsJournal of Chemical Theory and Computation, 2011
- Net charge per residue modulates conformational ensembles of intrinsically disordered proteinsProceedings of the National Academy of Sciences, 2010
- Constant pH Replica Exchange Molecular Dynamics in Biomolecules Using a Discrete Protonation ModelJournal of Chemical Theory and Computation, 2010
- Coupling Constant pH Molecular Dynamics with Accelerated Molecular DynamicsJournal of Chemical Theory and Computation, 2010
- CHARMM: The biomolecular simulation programJournal of Computational Chemistry, 2009
- Exploring atomistic details of pH-dependent peptide foldingProceedings of the National Academy of Sciences, 2006
- Replica-exchange molecular dynamics method for protein foldingChemical Physics Letters, 1999