1
|
Bowles J, Jähnigen S, Agostini F, Vuilleumier R, Zehnacker A, Calvo F, Clavaguéra C. Vibrational Circular Dichroism Spectroscopy with a Classical Polarizable Force Field: Alanine in the Gas and Condensed Phases. Chemphyschem 2024; 25:e202300982. [PMID: 38318765 DOI: 10.1002/cphc.202300982] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2023] [Revised: 02/05/2024] [Accepted: 02/06/2024] [Indexed: 02/07/2024]
Abstract
Polarizable force fields are an essential component for the chemically accurate modeling of complex molecular systems with a significant degree of fluxionality, beyond harmonic or perturbative approximations. In this contribution we examine the performance of such an approach for the vibrational spectroscopy of the alanine amino acid, in the gas and condensed phases, from the Fourier transform of appropriate time correlation functions generated along molecular dynamics (MD) trajectories. While the infrared (IR) spectrum only requires the electric dipole moment, the vibrational circular dichroism (VCD) spectrum further requires knowledge of the magnetic dipole moment, for which we provide relevant expressions to be used with polarizable force fields. The AMOEBA force field was employed here to model alanine in the neutral and zwitterionic isolated forms, solvated by water or nitrogen, and as a crystal. Within this framework, comparison of the electric and magnetic dipole moments to those obtained with nuclear velocity perturbation theory based on density-functional theory for the same MD trajectories are found to agree well with one another. The statistical convergence of the IR and VCD spectra is examined and found to be more demanding in the latter case. Comparisons with experimental frequencies are also provided for the condensed phases.
Collapse
Affiliation(s)
- Jessica Bowles
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405, Orsay, France
| | - Sascha Jähnigen
- PASTEUR Laboratory, Département de Chimie, Ecole Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005, Paris, France
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405, Orsay, France
| | - Rodolphe Vuilleumier
- PASTEUR Laboratory, Département de Chimie, Ecole Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005, Paris, France
| | - Anne Zehnacker
- Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d'Orsay UMR8214, 91405, Orsay, France
| | - Florent Calvo
- Université Grenoble Alpes, CNRS, LIPhy, 38000, Grenoble, France
| | - Carine Clavaguéra
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405, Orsay, France
| |
Collapse
|
2
|
Stylianakis I, Zervos N, Lii JH, Pantazis DA, Kolocouris A. Conformational energies of reference organic molecules: benchmarking of common efficient computational methods against coupled cluster theory. J Comput Aided Mol Des 2023; 37:607-656. [PMID: 37597063 PMCID: PMC10618395 DOI: 10.1007/s10822-023-00513-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2023] [Accepted: 06/03/2023] [Indexed: 08/21/2023]
Abstract
We selected 145 reference organic molecules that include model fragments used in computer-aided drug design. We calculated 158 conformational energies and barriers using force fields, with wide applicability in commercial and free softwares and extensive application on the calculation of conformational energies of organic molecules, e.g. the UFF and DREIDING force fields, the Allinger's force fields MM3-96, MM3-00, MM4-8, the MM2-91 clones MMX and MM+, the MMFF94 force field, MM4, ab initio Hartree-Fock (HF) theory with different basis sets, the standard density functional theory B3LYP, the second-order post-HF MP2 theory and the Domain-based Local Pair Natural Orbital Coupled Cluster DLPNO-CCSD(T) theory, with the latter used for accurate reference values. The data set of the organic molecules includes hydrocarbons, haloalkanes, conjugated compounds, and oxygen-, nitrogen-, phosphorus- and sulphur-containing compounds. We reviewed in detail the conformational aspects of these model organic molecules providing the current understanding of the steric and electronic factors that determine the stability of low energy conformers and the literature including previous experimental observations and calculated findings. While progress on the computer hardware allows the calculations of thousands of conformations for later use in drug design projects, this study is an update from previous classical studies that used, as reference values, experimental ones using a variety of methods and different environments. The lowest mean error against the DLPNO-CCSD(T) reference was calculated for MP2 (0.35 kcal mol-1), followed by B3LYP (0.69 kcal mol-1) and the HF theories (0.81-1.0 kcal mol-1). As regards the force fields, the lowest errors were observed for the Allinger's force fields MM3-00 (1.28 kcal mol-1), ΜΜ3-96 (1.40 kcal mol-1) and the Halgren's MMFF94 force field (1.30 kcal mol-1) and then for the MM2-91 clones MMX (1.77 kcal mol-1) and MM+ (2.01 kcal mol-1) and MM4 (2.05 kcal mol-1). The DREIDING (3.63 kcal mol-1) and UFF (3.77 kcal mol-1) force fields have the lowest performance. These model organic molecules we used are often present as fragments in drug-like molecules. The values calculated using DLPNO-CCSD(T) make up a valuable data set for further comparisons and for improved force field parameterization.
Collapse
Affiliation(s)
- Ioannis Stylianakis
- Department of Medicinal Chemistry, Faculty of Pharmacy, National and Kapodistrian University of Athens, Panepistimioupolis Zografou, 15771, Athens, Greece
| | - Nikolaos Zervos
- Department of Medicinal Chemistry, Faculty of Pharmacy, National and Kapodistrian University of Athens, Panepistimioupolis Zografou, 15771, Athens, Greece
| | - Jenn-Huei Lii
- Department of Chemistry, National Changhua University of Education, Changhua City, Taiwan
| | - Dimitrios A Pantazis
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470, Mülheim an der Ruhr, Germany
| | - Antonios Kolocouris
- Department of Medicinal Chemistry, Faculty of Pharmacy, National and Kapodistrian University of Athens, Panepistimioupolis Zografou, 15771, Athens, Greece.
- Laboratory of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department of Pharmacy, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou, 15771, Athens, Greece.
| |
Collapse
|
3
|
Sun Z, He Q, Gong Z, Kalhor P, Huai Z, Liu Z. A General Picture of Cucurbit[8]uril Host–Guest Binding: Recalibrating Bonded Interactions. Molecules 2023; 28:molecules28073124. [PMID: 37049887 PMCID: PMC10095826 DOI: 10.3390/molecules28073124] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2023] [Revised: 03/15/2023] [Accepted: 03/27/2023] [Indexed: 04/03/2023] Open
Abstract
Atomic-level understanding of the dynamic feature of host–guest interactions remains a central challenge in supramolecular chemistry. The remarkable guest binding behavior of the Cucurbiturils family of supramolecular containers makes them promising drug carriers. Among Cucurbit[n]urils, Cucurbit[8]uril (CB8) has an intermediate portal size and cavity volume. It can exploit almost all host–guest recognition motifs formed by this host family. In our previous work, an extensive computational investigation of the binding of seven commonly abused and structurally diverse drugs to the CB8 host was performed, and a general dynamic binding picture of CB8-guest interactions was obtained. Further, two widely used fixed-charge models for drug-like molecules were investigated and compared in great detail, aiming at providing guidelines in choosing an appropriate charge scheme in host-guest modelling. Iterative refitting of atomic charges leads to improved binding thermodynamics and the best root-mean-squared deviation from the experimental reference is 2.6 kcal/mol. In this work, we focus on a thorough evaluation of the remaining parts of classical force fields, i.e., the bonded interactions. The widely used general Amber force fields are assessed and refitted with generalized force-matching to improve the intra-molecular conformational preference, and thus the description of inter-molecular host–guest interactions. The interaction pattern and binding thermodynamics show a significant dependence on the modelling parameters. The refitted system-specific parameter set improves the consistency of the modelling results and the experimental reference significantly. Finally, combining the previous charge-scheme comparison and the current force-field refitting, we provide general guidelines for the theoretical modelling of host–guest binding.
Collapse
|
4
|
Chakraborty D, Banerjee A, Wales DJ. Side-Chain Polarity Modulates the Intrinsic Conformational Landscape of Model Dipeptides. J Phys Chem B 2021; 125:5809-5822. [PMID: 34037392 PMCID: PMC8279551 DOI: 10.1021/acs.jpcb.1c02412] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
The
intrinsic conformational preferences of small peptides may
provide additional insight into the thermodynamics and kinetics of
protein folding. In this study, we explore the underlying energy landscapes
of two model peptides, namely, Ac-Ala-NH2 and Ac-Ser-NH2, using geometry-optimization-based tools developed within
the context of energy landscape theory. We analyze not only how side-chain
polarity influences the structural preferences of the dipeptides,
but also other emergent properties of the landscape, including heat
capacity profiles, and kinetics of conformational rearrangements.
The contrasting topographies of the free energy landscape agree with
recent results from Fourier transform microwave spectroscopy experiments,
where Ac-Ala-NH2 was found to exist as a mixture of two
conformers, while Ac-Ser-NH2 remained structurally locked,
despite exhibiting an apparently rich conformational landscape.
Collapse
Affiliation(s)
- Debayan Chakraborty
- Department of Chemistry, The University of Texas at Austin, 24th Street Stop A5300, Austin, Texas 78712, United States
| | - Atreyee Banerjee
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom.,Max Planck Institute for Polymer Research, 55128 Mainz, Germany
| | - David J Wales
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
5
|
Lewis-Atwell T, Townsend PA, Grayson MN. Comparisons of different force fields in conformational analysis and searching of organic molecules: A review. Tetrahedron 2021. [DOI: 10.1016/j.tet.2020.131865] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
|
6
|
Otero-de-la-Roza A, DiLabio GA. Transferable Atom-Centered Potentials for the Correction of Basis Set Incompleteness Errors in Density-Functional Theory. J Chem Theory Comput 2017. [DOI: 10.1021/acs.jctc.7b00300] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- A. Otero-de-la-Roza
- Department
of Chemistry and ‡Faculty of Management, University of British Columbia, Okanagan, 3247
University Way, Kelowna, British Columbia, Canada V1V 1V7
| | - Gino A. DiLabio
- Department
of Chemistry and ‡Faculty of Management, University of British Columbia, Okanagan, 3247
University Way, Kelowna, British Columbia, Canada V1V 1V7
| |
Collapse
|
7
|
Balaji GA, Nagendra HG, Balaji VN, Rao SN. Experimental conformational energy maps of proteins and peptides. Proteins 2017; 85:979-1001. [PMID: 28168743 DOI: 10.1002/prot.25266] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2016] [Revised: 01/26/2017] [Accepted: 01/30/2017] [Indexed: 01/26/2023]
Abstract
We have presented an extensive analysis of the peptide backbone dihedral angles in the PDB structures and computed experimental Ramachandran plots for their distributions seen under a various constraints on X-ray resolution, representativeness at different sequence identity percentages, and hydrogen bonding distances. These experimental distributions have been converted into isoenergy contour plots using the approach employed previously by F. M. Pohl. This has led to the identification of energetically favored minima in the Ramachandran (ϕ, ψ) plots in which global minima are predominantly observed either in the right-handed α-helical or the polyproline II regions. Further, we have identified low energy pathways for transitions between various minima in the (ϕ,ψ) plots. We have compared and presented the experimental plots with published theoretical plots obtained from both molecular mechanics and quantum mechanical approaches. In addition, we have developed and employed a root mean square deviation (RMSD) metric for isoenergy contours in various ranges, as a measure (in kcal.mol-1 ) to compare any two plots and determine the extent of correlation and similarity between their isoenergy contours. In general, we observe a greater degree of compatibility with experimental plots for energy maps obtained from molecular mechanics methods compared to most quantum mechanical methods. The experimental energy plots we have investigated could be helpful in refining protein structures obtained from X-ray, NMR, and electron microscopy and in refining force field parameters to enable simulations of peptide and protein structures that have higher degree of consistency with experiments. Proteins 2017; 85:979-1001. © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Govardhan A Balaji
- Department of Biotechnology, Sir M Visvesvaraya Institute of Technology, Bangalore, 562157, India
| | - H G Nagendra
- Department of Biotechnology, Sir M Visvesvaraya Institute of Technology, Bangalore, 562157, India
| | - Vitukudi N Balaji
- Department of Biotechnology, Sir M Visvesvaraya Institute of Technology, Bangalore, 562157, India
| | - Shashidhar N Rao
- Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey, 08552
| |
Collapse
|
8
|
Baggioli A, Cavallotti CA, Famulari A. Exploring short intramolecular interactions in alkylaromatic substrates. Phys Chem Chem Phys 2016; 18:29616-29628. [PMID: 27753437 DOI: 10.1039/c6cp03323g] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
From proteins and peptides to semiconducting polymers, aliphatic chains on aromatic groups are recurring motifs in macromolecules from very diverse application fields. Fields in which molecular folding and packing determine the macroscopic physical properties that make such advanced materials appealing in the first place. Within each macromolecule, the intrinsic structure of each unit defines how it interacts with its neighbours, ultimately opening up or denying certain backbone conformations. This eventually also determines how macromolecules interact with each other. This account deals specifically with the conformational problem of many common alkylaromatic units, examining the features of an intramolecular interaction involving a side chain with as few as three methylene groups. A set of 23 model compounds featuring an intramolecular interaction between an aliphatic X-H (X = C, N, O, and S) bond and an aromatic ring was considered. Quantitative computational analysis was made possible, thanks to complete basis set extrapolated CCSD(T) calculations and NCI topological analysis, the latter of which revealed an elaborate network of dispersive and steric interactions leading to somewhat unintuitive and unexpected results, such as the higher energetic stability of certain twisted conformational isomers over those with extended side chains. Vicinal covalent effects from polarizing groups and various heteroatoms, along with the occurrence of non-dispersive phenomena, were also investigated. The conclusions drawn from the investigation include a comprehensive set of guidelines intended to aid in the prediction of the most stable conformation for this class of building blocks. Our findings affect a variety of different research fields, including the tailoring of functional materials for organic electronics and photovoltaics, with insights into a rational treatment of conformational disorder, and the study of protein- and peptide-folding preferences, putting an emphasis on peculiar interactions between the backbone and aromatic residues.
Collapse
Affiliation(s)
- Alberto Baggioli
- Dipartimento di Chimica, Materiali e Ingegneria Chimica "G. Natta", Politecnico di Milano, via Mancinelli 7, I-20131 Milano, Italy.
| | - Carlo A Cavallotti
- Dipartimento di Chimica, Materiali e Ingegneria Chimica "G. Natta", Politecnico di Milano, via Mancinelli 7, I-20131 Milano, Italy.
| | - Antonino Famulari
- Dipartimento di Chimica, Materiali e Ingegneria Chimica "G. Natta", Politecnico di Milano, via Mancinelli 7, I-20131 Milano, Italy.
| |
Collapse
|
9
|
Jangra H, Haindl MH, Achrainer F, Hioe J, Gschwind RM, Zipse H. Conformational Preferences in Small Peptide Models: The Relevance ofcis/trans-Conformations. Chemistry 2016; 22:13328-35. [DOI: 10.1002/chem.201601828] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2016] [Indexed: 11/12/2022]
Affiliation(s)
- Harish Jangra
- Department of Chemistry; LMU München; Butenandstrasse 5-14 81377 München Germany
| | - Michael H. Haindl
- Institut für Organische Chemie; Universität Regensburg; 93053 Regensburg Germany
| | - Florian Achrainer
- Department of Chemistry; LMU München; Butenandstrasse 5-14 81377 München Germany
| | - Johnny Hioe
- Institut für Organische Chemie; Universität Regensburg; 93053 Regensburg Germany
| | - Ruth M. Gschwind
- Institut für Organische Chemie; Universität Regensburg; 93053 Regensburg Germany
| | - Hendrik Zipse
- Department of Chemistry; LMU München; Butenandstrasse 5-14 81377 München Germany
| |
Collapse
|
10
|
Assessing the performance of popular QM methods for calculation of conformational energies of trialanine. Chem Phys Lett 2016. [DOI: 10.1016/j.cplett.2016.04.037] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
11
|
Kaminský J, Jensen F. Conformational Interconversions of Amino Acid Derivatives. J Chem Theory Comput 2016; 12:694-705. [PMID: 26691979 DOI: 10.1021/acs.jctc.5b00911] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Exhaustive conformational interconversions including transition structure analyses of N-acetyl-l-glycine-N-methylamide as well as its alanine, serine, and cysteine analogues have been investigated at the MP2/6-31G** level, yielding a total of 142 transition states. Improved estimates of relative energies were obtained by separately extrapolating the Hartree-Fock and MP2 energies to the basis set limit and adding the difference between CCSD(T) and MP2 results with the cc-pVDZ basis set to the extrapolated MP2 results. The performance of eight empirical force fields (AMBER94, AMBER14SB, MM2, MM3, MMFFs, CHARMM22_CMAP, OPLS_2005, and AMOEBAPRO13) in reproducing ab initio energies of transition states was tested. Our results indicate that commonly used class I force fields employing a fixed partial charge model for the electrostatic interaction provide mean errors in the ∼10 kJ/mol range for energies of conformational transition states for amino acid conformers. Modern reparametrized versions, such as CHARMM22_CMAP, and polarizable force fields, such as AMOEBAPRO13, have slightly lower mean errors, but maximal errors are still in the 35 kJ/mol range. There are differences between the force fields in their ability for reproducing conformational transitions classified according to backbone/side-chain or regions in the Ramachandran angles, but the data set is likely too small to draw any general conclusions. Errors in conformational interconversion barriers by ∼10 kJ/mol suggest that the commonly used force field may bias certain types of transitions by several orders of magnitude in rate and thus lead to incorrect dynamics in simulations. It is therefore suggested that information for conformational transition states should be included in parametrizations of new force fields.
Collapse
Affiliation(s)
- Jakub Kaminský
- Institute of Organic Chemistry and Biochemistry, Flemingovo nám. 2, 166 10 Prague, Czech Republic
| | - Frank Jensen
- Department of Chemistry, Aarhus University , Langelandsgade 140, DK-8000 Aarhus C, Denmark
| |
Collapse
|
12
|
Jensen F. An Atomic Counterpoise Method for Estimating Inter- and Intramolecular Basis Set Superposition Errors. J Chem Theory Comput 2015; 6:100-6. [PMID: 26614323 DOI: 10.1021/ct900436f] [Citation(s) in RCA: 59] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
An atomic counterpoise method is proposed to calculate estimates of inter- and intramolecular basis set superposition errors. The method estimates the basis set superposition error as a sum of atomic contributions and can be applied for both independent particle and electron correlation models. It is shown that the atomic counterpoise method provides results very similar to the molecular counterpoise method for intermolecular basis set superposition errors at both the HF and MP2 levels of theory with a sequence of increasingly larger basis sets. The advantage of the atomic counterpoise method is that it can be applied with equal ease to estimate intramolecular basis set superposition errors, for which few other methods exist. The atomic counterpoise method is computationally quite efficient, requiring typically double the amount of computer time as required for calculating the uncorrected energy.
Collapse
Affiliation(s)
- Frank Jensen
- Department of Chemistry, University of Aarhus Langelandsgade 140, DK-8000 Aarhus, Denmark
| |
Collapse
|
13
|
Riffet V, Bouchoux G, Frison G. Microhydration of Protonated Nα-Acetylhistidine: A Theoretical Approach. J Phys Chem B 2015; 119:11527-39. [PMID: 26252717 DOI: 10.1021/acs.jpcb.5b05581] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Extensive exploration of the potential energy surfaces of protonated Nα-acetylhistidine hydrated by 0-3 molecules of water was performed. The methodology combined hierarchical and genealogical (Darwin family tree) approaches using polarizable AMOEBA force field and M06 functional. It is demonstrated that this mixed approach allows recovering a larger number of conformers than the number recovered by using any one of the two methods alone. Hydration enthalpies of protonated Nα-acetylhistidine and of model compounds have been computed using higher theoretical methods, up to the G4MP2 procedure. Excellent agreement with experiment is observed for successive hydration of methylamonium and imidazolium cations using MP2/6-311++G(2d,2p)//M06/6-311++G(d,p) and G4MP2 methods, thereby validating the theory levels used for hydrated protonated Nα-acetylhistidine. It is found that the first hydration enthalpy of protonated Nα-acetylhistidine is ca. 10 kJ mol(-1) lower than that of imidazolium, a result explained by the local environment of the positively charged imidazolium moiety.
Collapse
Affiliation(s)
- Vanessa Riffet
- Laboratoire de Chimie Moléculaire, Ecole Polytechnique and CNRS , 91128 Palaiseau cedex, France
| | - Guy Bouchoux
- Laboratoire de Chimie Moléculaire, Ecole Polytechnique and CNRS , 91128 Palaiseau cedex, France
| | - Gilles Frison
- Laboratoire de Chimie Moléculaire, Ecole Polytechnique and CNRS , 91128 Palaiseau cedex, France
| |
Collapse
|
14
|
Antila HS, Salonen E. On combining Thole's induced point dipole model with fixed charge distributions in molecular mechanics force fields. J Comput Chem 2015; 36:739-50. [PMID: 25753482 DOI: 10.1002/jcc.23850] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2014] [Revised: 12/19/2014] [Accepted: 01/11/2015] [Indexed: 11/12/2022]
Abstract
The Thole induced point dipole model is combined with three different point charge fitting methods, Merz-Kollman (MK), charges from electrostatic potentials using a grid (CHELPG), and restrained electrostatic potential (RESP), and two multipole algorithms, distributed multipole analysis (DMA) and Gaussian multipole model (GMM), which can be used to describe the electrostatic potential (ESP) around molecules in molecular mechanics force fields. This is done to study how the different methods perform when intramolecular polarizability contributions are self-consistently removed from the fitting done in the force field parametrization. It is demonstrated that the polarizable versions of the partial charge models provide a good compromise between accuracy and computational efficiency in describing the ESP of small organic molecules undergoing conformational changes. For the point charge models, the inclusion of polarizability reduced the the average root mean square error of ESP over the test set by 4-10%.
Collapse
|
15
|
Effect of the external electric field on selected tripeptides. Amino Acids 2015; 47:1399-408. [DOI: 10.1007/s00726-015-1971-8] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2014] [Accepted: 03/23/2015] [Indexed: 10/23/2022]
|
16
|
Topham CM, Smith JC. Tri-peptide reference structures for the calculation of relative solvent accessible surface area in protein amino acid residues. Comput Biol Chem 2014; 54:33-43. [PMID: 25544680 DOI: 10.1016/j.compbiolchem.2014.11.007] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2014] [Revised: 11/23/2014] [Accepted: 11/30/2014] [Indexed: 10/24/2022]
Abstract
Relative amino acid residue solvent accessibility values allow the quantitative comparison of atomic solvent-accessible surface areas in different residue types and physical environments in proteins and in protein structural alignments. Geometry-optimised tri-peptide structures in extended solvent-exposed reference conformations have been obtained for 43 amino acid residue types at a high level of quantum chemical theory. Significant increases in side-chain solvent accessibility, offset by reductions in main-chain atom solvent exposure, were observed for standard residue types in partially geometry-optimised structures when compared to non-minimised models built from identical sets of proper dihedral angles abstracted from the literature. Optimisation of proper dihedral angles led most notably to marked increases of up to 54% in proline main-chain atom solvent accessibility compared to literature values. Similar effects were observed for fully-optimised tri-peptides in implicit solvent. The relief of internal strain energy was associated with systematic variation in N, C(α) and C(β) atom solvent accessibility across all standard residue types. The results underline the importance of optimisation of 'hard' degrees of freedom (bond lengths and valence bond angles) and improper dihedral angle values from force field or other context-independent reference values, and impact on the use of standardised fixed internal co-ordinate geometry in sampling approaches to the determination of absolute values of protein amino acid residue solvent accessibility. Quantum chemical methods provide a useful and accurate alternative to molecular mechanics methods to perform energy minimisation of peptides containing non-standard (chemically modified) amino acid residues frequently present in experimental protein structure data sets, for which force field parameters may not be available. Reference tri-peptide atomic co-ordinate sets including hydrogen atoms are made freely available.
Collapse
Affiliation(s)
- Christopher M Topham
- Molecular Forces Consulting, 40 Rue Boyssonne, Toulouse 31400, France; Computational Molecular Biophysics, IWR der Universität Heidelberg, Im Neuenheimer Feld 368, Heidelberg D-69120, Germany; University of Tennessee/Oak Ridge National Laboratory, Center for Molecular Biophysics, P.O. Box 2008, Oak Ridge, TN 37831-6309, USA; Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, M407 Walters Life Sciences, 1414 Cumberland Avenue, Knoxville, TN 37996, USA.
| | - Jeremy C Smith
- Computational Molecular Biophysics, IWR der Universität Heidelberg, Im Neuenheimer Feld 368, Heidelberg D-69120, Germany; University of Tennessee/Oak Ridge National Laboratory, Center for Molecular Biophysics, P.O. Box 2008, Oak Ridge, TN 37831-6309, USA; Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, M407 Walters Life Sciences, 1414 Cumberland Avenue, Knoxville, TN 37996, USA
| |
Collapse
|
17
|
Gregori B, Guidoni L, Chiavarino B, Scuderi D, Nicol E, Frison G, Fornarini S, Crestoni ME. Vibrational Signatures of S-Nitrosoglutathione as Gaseous, Protonated Species. J Phys Chem B 2014; 118:12371-82. [DOI: 10.1021/jp5072742] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
| | - Leonardo Guidoni
- Dipartimento
di Scienza Fisiche e Chimiche, Università degli Studi dell’Aquila, Via Vetoio 2, Coppito, L’Aquila I-64100, Italy
| | | | - Debora Scuderi
- Laboratoire
de Chimie Physique, UMR8000 CNRS, Faculté des Sciences, Université Paris-Sud, Batiment 350, 91405 Orsay Cedex, France
| | - Edith Nicol
- Laboratoire
de Chimie Moléculaire, Ecole Polytechnique and CNRS, 91128 Palaiseau Cedex, France
| | - Gilles Frison
- Laboratoire
de Chimie Moléculaire, Ecole Polytechnique and CNRS, 91128 Palaiseau Cedex, France
| | | | | |
Collapse
|
18
|
Xue Y, Yuwen T, Zhu F, Skrynnikov NR. Role of electrostatic interactions in binding of peptides and intrinsically disordered proteins to their folded targets. 1. NMR and MD characterization of the complex between the c-Crk N-SH3 domain and the peptide Sos. Biochemistry 2014; 53:6473-95. [PMID: 25207671 DOI: 10.1021/bi500904f] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
Intrinsically disordered proteins (IDPs) often rely on electrostatic interactions to bind their structured targets. To obtain insight into the mechanism of formation of the electrostatic encounter complex, we investigated the binding of the peptide Sos (PPPVPPRRRR), which serves as a minimal model for an IDP, to the c-Crk N-terminal SH3 domain. Initially, we measured ¹⁵N relaxation rates at two magnetic field strengths and determined the binding shifts for the complex of Sos with wild-type SH3. We have also recorded a 3 μs molecular dynamics (MD) trajectory of this complex using the Amber ff99SB*-ILDN force field. The comparison of the experimental and simulated data shows that MD simulation consistently overestimates the strength of salt bridge interactions at the binding interface. The series of simulations using other advanced force fields also failed to produce any satisfactory results. To address this issue, we have devised an empirical correction to the Amber ff99SB*-ILDN force field whereby the Lennard-Jones equilibrium distance for the nitrogen-oxygen pair across the Arg-to-Asp and Arg-to-Glu salt bridges has been increased by 3%. Implementing this correction resulted in a good agreement between the simulations and the experiment. Adjusting the strength of salt bridge interactions removed a certain amount of strain contained in the original MD model, thus improving the binding of the hydrophobic N-terminal portion of the peptide. The arginine-rich C-terminal portion of the peptide, freed from the effect of the overstabilized salt bridges, was found to interconvert more rapidly between its multiple conformational states. The modified MD protocol has also been successfully used to simulate the entire binding process. In doing so, the peptide was initially placed high above the protein surface. It then arrived at the correct bound pose within ∼2 Å of the crystallographic coordinates. This simulation allowed us to analyze the details of the dynamic binding intermediate, i.e., the electrostatic encounter complex. However, an experimental characterization of this transient, weakly populated state remains out of reach. To overcome this problem, we designed the double mutant of c-Crk N-SH3 in which mutations Y186L and W169F abrogate tight Sos binding and shift the equilibrium toward the intermediate state resembling the electrostatic encounter complex. The results of the combined NMR and MD study of this engineered system will be reported in the next part of this paper.
Collapse
Affiliation(s)
- Yi Xue
- Department of Chemistry, Purdue University , West Lafayette, Indiana 47907, United States
| | | | | | | |
Collapse
|
19
|
Yuan Y, Mills MJL, Popelier PLA, Jensen F. Comprehensive analysis of energy minima of the 20 natural amino acids. J Phys Chem A 2014; 118:7876-91. [PMID: 25084473 DOI: 10.1021/jp503460m] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Energy minima of the 20 natural amino acids (capped by a peptide bond at both the N and the C termini, CH3-C(═O)-N(H)-(H)Cα(R)-C(═O)-N(H)-CH3), were obtained by ab initio geometry optimization. Starting with a large number of minima, quickly generated by MarvinView, geometry optimization at the HF/6-31G(d,p) level of theory reduced the number of minima, followed by further optimization at the B3LYP/apc-1 and MP2/cc-pVDZ levels, which caused some minima to disappear and some stable minima to migrate on the Ramachandran map. There is a relation between the number of minima and the size and the flexibility of the side chain. The energy minima of the 20 amino acids are mainly located in the regions of βL, γL, δL, and αL of the Ramachandran map. Multipole moments of atoms occurring in the fragment [-NH-Cα-C(═O)-] common to all 20 amino acids were calculated at the three levels of theory mentioned above. The near parallelism in behavior of these moments between levels of theory is beneficial toward estimating moments with the more expensive B3LYP and MP2 methods from data calculated with the cheaper HF method. Finally, we explored the transferability of properties between different amino acids: the bond length and angles of the common fragment [-NH-Cα(HαCβ)-C'(═O)-] in all amino acids except Gly and Pro. All bond lengths are highly transferable between different amino acids, and the standard deviations are small.
Collapse
Affiliation(s)
- Yongna Yuan
- Manchester Institute of Biotechnology (MIB) , 131 Princess Street, Manchester M1 7DN, Great Britain and
| | | | | | | |
Collapse
|
20
|
Dong C, Yong-Zhi L, Zhi-Chao W, Bo L. Performance of four different force fields for simulations of dipeptide conformations: GlyGly, GlyGly-, GlyGly · Cl-, GlyGly · Na+ and GlyGly · (H2O)2. J Mol Model 2014; 20:2279. [PMID: 24863532 DOI: 10.1007/s00894-014-2279-4] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2013] [Accepted: 04/25/2014] [Indexed: 11/30/2022]
Abstract
Monte-Carlo conformational searches with four kinds of force fields (AMBER94, MM3*, MMFFs, and OPLS-2005) were performed on glycylglycine (GlyGly), deprotonated glycylglycine (GlyGly(-)), glycylglycine chloride anion complex (GlyGly · Cl(-)), glycylglycine sodium cation complex (GlyGly · Na(+)) and glycylglycine dihydrate [GlyGly · (H2O)2]. Combined with Hartree-Fock (HF) and second-order Møller-Plesset (MP2) optimizations, conformations within an energy of 20 kJ mol(-1) were predicted. After MP2 calculations, the geometries and relative energies of the predicted structures were the same regardless of the force field used. Therefore, the performance of different force fields reflects mainly the conformational search process. For GlyGly, there was practically no difference among the four force fields. Due to the complex hydrogen bonding network when involving water, the total number of resulting conformers for GlyGly · (H2O)2 increased drastically. Moreover, the MMFFs force field fared best in finding the global minimum compared to the remaining three force fields. In describing hydrogen bonded and inter-molecular complexes, we recommend application of the MMFFs and AMBER94 force fields. Furthermore, the MMFFs and OPLS-2005 force fields have a good description of electrostatic interactions. This work will contribute to helping the reader make an optimal choice of force field, taking into account the latter's strengths and limitations.
Collapse
Affiliation(s)
- Chen Dong
- Institute of Photo-Biophysics, Physics and Electronics Department, Henan University, 475004, Kaifeng, China
| | | | | | | |
Collapse
|
21
|
Kang YK, Park HS. Assessment of CCSD(T), MP2, DFT-D, CBS-QB3, and G4(MP2) methods for conformational study of alanine and proline dipeptides. Chem Phys Lett 2014. [DOI: 10.1016/j.cplett.2014.03.067] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
22
|
Cisneros GA, Karttunen M, Ren P, Sagui C. Classical electrostatics for biomolecular simulations. Chem Rev 2014; 114:779-814. [PMID: 23981057 PMCID: PMC3947274 DOI: 10.1021/cr300461d] [Citation(s) in RCA: 192] [Impact Index Per Article: 19.2] [Reference Citation Analysis] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
23
|
|
24
|
Walker M, Sen A, Harvey AJ, Dessent CE. Complexation of anions to gas-phase amino acids: Conformation is critical in determining if the global minimum is canonical or zwitterionic. Chem Phys Lett 2013. [DOI: 10.1016/j.cplett.2013.09.074] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
25
|
Hoffgaard F, Heil J, Kast SM. Three-Dimensional RISM Integral Equation Theory for Polarizable Solute Models. J Chem Theory Comput 2013; 9:4718-26. [PMID: 26583390 DOI: 10.1021/ct400699q] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
Modeling solute polarizability is a key ingredient for improving the description of solvation phenomena. In recent years, polarizable molecular mechanics force fields have emerged that circumvent the limitations of classical fixed charge force fields by the ability to adapt their electrostatic potential distribution to a polarizing environment. Solvation phenomena are characterized by the solute's excess chemical potential, which can be computed by expensive fully atomistic free energy simulations. The alternative is to employ an implicit solvent model, which poses a challenge to the formulation of the solute-solvent interaction term within a polarizable framework. Here, we adapt the three-dimensional reference interaction site model (3D RISM) integral equation theory as a solvent model, which analytically yields the chemical potential, to the polarizable AMOEBA force field using an embedding cluster (EC-RISM) strategy. The methodology is analogous to our earlier approach to the coupling of a quantum-chemical solute description with a classical 3D RISM solvent. We describe the conceptual physical and algorithmic basis as well as the performance for several benchmark cases as a proof of principle. The results consistently show reasonable agreement between AMOEBA and quantum-chemical free energies in solution in general and allow for separate assessment of energetic and solvation-related contributions. We find that, depending on the parametrization, AMOEBA reproduces the chemical potential in better agreement with reference quantum-chemical calculations than the intramolecular energies, which suggests possible routes toward systematic improvement of polarizable force fields.
Collapse
Affiliation(s)
- Franziska Hoffgaard
- Physikalische Chemie III, TU Dortmund , Otto-Hahn-Str. 6, 44227 Dortmund, Germany
| | - Jochen Heil
- Physikalische Chemie III, TU Dortmund , Otto-Hahn-Str. 6, 44227 Dortmund, Germany
| | - Stefan M Kast
- Physikalische Chemie III, TU Dortmund , Otto-Hahn-Str. 6, 44227 Dortmund, Germany
| |
Collapse
|
26
|
Aliev AE, Kulke M, Khaneja HS, Chudasama V, Sheppard TD, Lanigan RM. Motional timescale predictions by molecular dynamics simulations: case study using proline and hydroxyproline sidechain dynamics. Proteins 2013; 82:195-215. [PMID: 23818175 PMCID: PMC4282583 DOI: 10.1002/prot.24350] [Citation(s) in RCA: 159] [Impact Index Per Article: 14.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2013] [Revised: 06/01/2013] [Accepted: 06/06/2013] [Indexed: 01/08/2023]
Abstract
We propose a new approach for force field optimizations which aims at reproducing dynamics characteristics using biomolecular MD simulations, in addition to improved prediction of motionally averaged structural properties available from experiment. As the source of experimental data for dynamics fittings, we use 13C NMR spin-lattice relaxation times T1 of backbone and sidechain carbons, which allow to determine correlation times of both overall molecular and intramolecular motions. For structural fittings, we use motionally averaged experimental values of NMR J couplings. The proline residue and its derivative 4-hydroxyproline with relatively simple cyclic structure and sidechain dynamics were chosen for the assessment of the new approach in this work. Initially, grid search and simplexed MD simulations identified large number of parameter sets which fit equally well experimental J couplings. Using the Arrhenius-type relationship between the force constant and the correlation time, the available MD data for a series of parameter sets were analyzed to predict the value of the force constant that best reproduces experimental timescale of the sidechain dynamics. Verification of the new force-field (termed as AMBER99SB-ILDNP) against NMR J couplings and correlation times showed consistent and significant improvements compared to the original force field in reproducing both structural and dynamics properties. The results suggest that matching experimental timescales of motions together with motionally averaged characteristics is the valid approach for force field parameter optimization. Such a comprehensive approach is not restricted to cyclic residues and can be extended to other amino acid residues, as well as to the backbone. Proteins 2014; 82:195–215. © 2013 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Abil E Aliev
- Department of Chemistry, University College London, London, WC1H 0AJ, United Kingdom
| | | | | | | | | | | |
Collapse
|
27
|
Abstract
Molecular dynamics simulations of biomolecules have matured into powerful tools of structural biology. In addition to the commonly used empirical force field potentials, quantum mechanical descriptions are gaining popularity for structure optimization and dynamic simulations of peptides and proteins. In this chapter, we introduce methodological developments such as the QM/MM framework and linear-scaling QM that make efficient calculations on large biomolecules possible. We identify the most common scenarios in which quantum descriptions of peptides and proteins are employed, such as structural refinement, force field development, treatment of unusual residues, and predicting spectroscopic and exited state properties. The benefits and shortcomings of QM potentials, in comparison to classical force fields, are discussed, with special emphasis on the sampling problems of protein conformational space. Finally, recent examples of QM/MM calculations in light-sensitive membrane proteins illustrate typical applications of the reviewed methods.
Collapse
Affiliation(s)
- Thomas Steinbrecher
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany
| | | |
Collapse
|
28
|
Semrouni D, Clavaguéra C, Ohanessian G, Parks JH. Relationship between Conformational Dynamics and Electron Transfer in a Desolvated Peptide. Part I. Structures. J Phys Chem B 2013; 117:1746-55. [DOI: 10.1021/jp3078375] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Affiliation(s)
- David Semrouni
- Laboratoire des Mécanismes Réactionnels,
Department of Chemistry, Ecole Polytechnique, CNRS, 91128 Palaiseau Cedex, France
| | - Carine Clavaguéra
- Laboratoire des Mécanismes Réactionnels,
Department of Chemistry, Ecole Polytechnique, CNRS, 91128 Palaiseau Cedex, France
| | - Gilles Ohanessian
- Laboratoire des Mécanismes Réactionnels,
Department of Chemistry, Ecole Polytechnique, CNRS, 91128 Palaiseau Cedex, France
| | - Joel H. Parks
- Rowland Institute at Harvard, 100 Edwin H. Land Boulevard, Cambridge,
Massachusetts 02142, United
States
| |
Collapse
|
29
|
Monti S, Corozzi A, Fristrup P, Joshi KL, Shin YK, Oelschlaeger P, van Duin ACT, Barone V. Exploring the conformational and reactive dynamics of biomolecules in solution using an extended version of the glycine reactive force field. Phys Chem Chem Phys 2013; 15:15062-77. [DOI: 10.1039/c3cp51931g] [Citation(s) in RCA: 87] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
30
|
Abstract
This chapter provides an overview of the most common methods for including an explicit description of electronic polarization in molecular mechanics force fields: the induced point dipole, shell, and fluctuating charge models. The importance of including polarization effects in biomolecular simulations is discussed, and some of the most important achievements in the development of polarizable biomolecular force fields to date are highlighted.
Collapse
Affiliation(s)
- Hanne S Antila
- Department of Chemistry, Aalto University, Espoo, Finland
| | | |
Collapse
|
31
|
Zhang J, Yang W, Piquemal JP, Ren P. Modeling Structural Coordination and Ligand Binding in Zinc Proteins with a Polarizable Potential. J Chem Theory Comput 2012; 8:1314-1324. [PMID: 22754403 PMCID: PMC3383645 DOI: 10.1021/ct200812y] [Citation(s) in RCA: 92] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
As the second most abundant cation in human body, zinc is vital for the structures and functions of many proteins. Zinc-containing matrix metalloproteinases (MMPs) have been widely investigated as potential drug targets in a range of diseases ranging from cardiovascular disorders to cancers. However, it remains a challenge in theoretical studies to treat zinc in proteins with classical mechanics. In this study, we examined Zn(2+) coordination with organic compounds and protein side chains using a polarizable atomic multipole based electrostatic model. We find that polarization effect plays a determining role in Zn(2+) coordination geometry in both matrix metalloproteinase (MMP) complexes and in zinc-finger proteins. In addition, the relative binding free energies of selected inhibitors binding with MMP13 have been estimated and compared with experimental results. While not directly interacting with the small molecule inhibitors, the permanent and polarizing field of Zn(2+) exerts a strong influence on the relative affinities of the ligands. The simulation results also reveal the polarization effect on binding is ligand dependent and thus difficult to be incorporated into fixed-charge models implicitly.
Collapse
Affiliation(s)
- Jiajing Zhang
- Department of Biomedical Engineering, The University of Texas at Austin, TX 78712
| | - Wei Yang
- The Institute of Molecular Biophysics, Florida State University, Tallahassee, FL 32306
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, FL 32306
| | - Jean-Philip Piquemal
- UPMC Univ. Paris 06, UMR 7616, Laboratoire de Chimie Théorique, case courrier 137, 4 place Jussieu, F-75005, Paris, France
- CNRS, UMR 7616, Laboratoire de Chimie Théorique, case courrier 137, 4 place Jussieu, F-75005, Paris, France
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, TX 78712
| |
Collapse
|
32
|
Atwood RE, Urban JJ. Conformations of the Glycine Tripeptide Analog Ac-Gly-Gly-NHMe: A Computational Study Including Aqueous Solvation Effects. J Phys Chem A 2012; 116:1396-408. [DOI: 10.1021/jp206152d] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Affiliation(s)
- Rex E. Atwood
- Chemistry Department, United States Naval Academy, 572 Holloway Road, Annapolis,
Maryland 21402, United States
| | - Joseph J. Urban
- Chemistry Department, United States Naval Academy, 572 Holloway Road, Annapolis,
Maryland 21402, United States
| |
Collapse
|
33
|
Aliev AE, Mia ZA, Khaneja HS, King FD. Structures in Solutions from Joint Experimental-Computational Analysis: Applications to Cyclic Molecules and Studies of Noncovalent Interactions. J Phys Chem A 2012; 116:1093-109. [DOI: 10.1021/jp211083f] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Affiliation(s)
- Abil E. Aliev
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Zakirin A. Mia
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Harmeet S. Khaneja
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Frank D. King
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| |
Collapse
|
34
|
Kim CK, Park BH, Lee HW, Kim CK. Comprehensive Studies on the Free Energies of Solvation and Conformers of Glycine: A Theoretical Study. B KOREAN CHEM SOC 2011. [DOI: 10.5012/bkcs.2011.32.6.1985] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
35
|
Affiliation(s)
- Julie Grouleff
- Department of Chemistry, Aarhus University, DK-8000 Aarhus, Denmark
| | - Frank Jensen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus, Denmark
| |
Collapse
|
36
|
Hioe J, Savasci G, Brand H, Zipse H. The stability of Cα peptide radicals: why glycyl radical enzymes? Chemistry 2011; 17:3781-9. [PMID: 21341321 DOI: 10.1002/chem.201002620] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2010] [Indexed: 11/06/2022]
Abstract
The conformational space of dipeptide models derived from glycine, alanine, phenylalanine, proline, tyrosine, and cysteine has been searched extensively and compared with the corresponding C(α) dipeptide radicals at the G3(MP2)-RAD level of theory. The results indicate that the (least-substituted) glycine dipeptide radical is the thermochemically most stable of these species. Analysis of the structural parameters indicates that this is due to repulsive interactions between the C(α) substituents and peptide units in the radical. A comparison of the conformational preferences of dipeptide radicals and their closed-shell parents also indicates that radical stability is a strongly conformation-dependent property.
Collapse
Affiliation(s)
- Johnny Hioe
- Department of Chemistry, LMU München, Butenandtstrasse 5-13, 81377 München, Germany
| | | | | | | |
Collapse
|
37
|
Aliev AE, Courtier-Murias D. Experimental verification of force fields for molecular dynamics simulations using Gly-Pro-Gly-Gly. J Phys Chem B 2011; 114:12358-75. [PMID: 20825228 DOI: 10.1021/jp101581h] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Experimental NMR verification of MD simulations using 12 different force fields (AMBER, CHARMM, GROMOS, and OPLS-AA) and 5 different water models has been undertaken to identify reliable MD protocols for structure and dynamics elucidations of small open chain peptides containing Gly and Pro. A conformationally flexible tetrapeptide Gly-Pro-Gly-Gly was selected for NMR (3)J-coupling, chemical shift, and internuclear distance measurements, followed by their calculations using 2 μs long MD simulations in water. In addition, Ramachandran population maps for Pro-2 and Gly-3 residues of GPGG obtained from MD simulations were used for detailed comparisons with similar maps from the protein data bank (PDB) for large number of Gly and Pro residues in proteins. The MD simulations revealed strong dependence of the populations and geometries of preferred backbone and side chain conformations, as well as the time scales of the peptide torsional transitions on the force field used. On the basis of the analysis of the measured and calculated data, AMBER99SB is identified as the most reliable force field for reproducing NMR measured parameters, which are dependent on the peptide backbone and the Pro side chain geometries and dynamics. Ramachandran maps showing the dependence of conformational populations as a function of backbone ϕ/ψ angles for Pro-2 and Gly-3 residues of GPGG from MD simulations using AMBER99SB, AMBER03, and CHARMM were found to resemble similar maps for Gly and Pro residues from the PDB survey. Three force fields (AMBER99, AMBER99ϕ, and AMBER94) showed the least satisfactory agreement with both the solution NMR and the PDB survey data. The poor performance of these force fields is attributed to their propensity to overstabilize helical peptide backbone conformations at the Pro-2 and Gly-3 residues. On the basis of the similarity of the MD and PDB Ramachandran plots, the following sequence of transitions is suggested for the Gly backbone conformation: α(L) ⇆ β(PR) ⇆ β(S) ⇆ β(P) ⇆ α, where backbone secondary structures α(L) and α are associated with helices and turns, β(P) and β(PR) correspond to the left- and right-handed polyproline II structures and β(S) denotes the fully stretched backbone conformation. Compared to the force field dependence, less significant, but noteworthy, variations in the populations of the peptide backbone conformations were observed. For different solvent models considered, a correlation was noted between the number of torsional transitions in GPGG and the water self-diffusion coefficient on using TIP3P, TIP4P, and TIP5P models. In addition to MD results, we also report DFT derived Karplus relationships for Gly and Pro residues using B972 and B3LYP functionals.
Collapse
Affiliation(s)
- Abil E Aliev
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, UK.
| | | |
Collapse
|
38
|
Riffet V, Frison G, Bouchoux G. Acid–base thermochemistry of gaseous oxygen and sulfur substituted amino acids (Ser, Thr, Cys, Met). Phys Chem Chem Phys 2011; 13:18561-80. [DOI: 10.1039/c1cp22206f] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
39
|
Hawe GI, Popelier PL. A water potential based on multipole moments trained by machine learning — Reducing maximum energy errors. CAN J CHEM 2010. [DOI: 10.1139/v10-075] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
A potential that strives to represent the Coulomb interaction realistically must include polarization. In our approach, three decisions were made to accomplish this: (i) define an atom according to quantum chemical topology (QCT), (ii) express the interaction between atoms via their multipole moments, and (iii) use machine learning to capture the response of an atomic multipole moment to a change in this atom’s environment. This approach avoids explicit (distributed) polarizabilities and eliminates the problem of polarization catastrophe. Previously, we showed ( Phys. Chem. Chem. Phys. 2009, 11, 6365 ) that a machine learning method called kriging predicted atomic multipole moments more accurately than competing machine learning methods. This was established for the atoms of a central water molecule in water clusters, from the dimer to the hexamer. The prediction errors in all multipole moments were collectively assessed by errors in total interaction energy, for thousands of clusters configurations. Here, we target the maximum errors, with an eye on reducing the worst predictions that the potential may return. We demonstrate proof-of-principle for the water dimer using local kriging.
Collapse
Affiliation(s)
- Glenn I. Hawe
- Manchester Interdisciplinary Biocentre (MIB), 131 Princess Street, Manchester M1 7DN, UK
- School of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL, UK
| | - Paul L.A. Popelier
- Manchester Interdisciplinary Biocentre (MIB), 131 Princess Street, Manchester M1 7DN, UK
- School of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL, UK
| |
Collapse
|
40
|
Taylor CJ, Nix MGD, Dessent CEH. Noncovalent Interactions in the Gas-Phase Conformers of Anionic Iduronate (methyl 2-O-sulfo-α-L-iduronate): Variation of Subconformer versus Ring Conformer Energetics for a Prototypical Anionic Monosaccharide Studied Using Computational Methods. J Phys Chem A 2010; 114:11153-60. [DOI: 10.1021/jp102657t] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Affiliation(s)
| | - Michael G. D. Nix
- Department of Chemistry, University of York, Heslington, York, YO10 5DD, U.K
| | | |
Collapse
|
41
|
Dong L, Qin S, Su Z, Yang H, Hu C. Computational investigation on the mechanism and the stereoselectivity of Morita-Baylis-Hillman reaction and the effect of the bifunctional catalyst N-methylprolinol. Org Biomol Chem 2010; 8:3985-91. [PMID: 20623055 DOI: 10.1039/c004932h] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The mechanism of the Morita-Baylis-Hillman (MBH) reaction between formaldehyde and methyl vinyl ketone (MVK) catalyzed by N-methylprolinol was investigated using density functional theory (DFT) method. The overall reaction includes two steps: C-C bond formation and hydrogen migration. In the presence of water, the hydrogen migration occurs via a six-membered ring transition state and the corresponding energy barrier decreases dramatically, and therefore the RDS is the C-C bond formation step. The calculations indicate that the C-C bond formation step controls the stereochemistry of the reaction. In this step, the hydrogen bonding induces the direction of the attack of enamine to aldehyde from the -OH group side of N-methylprolinol. The energy-favored transition states are mainly stabilized by hydrogen bonding, while the chirality of the products is affected by the hydrogen bonding and the steric hindrance. The calculations correctly reproduce the major product in (R)-configuration, which is consistent with the experimental observation.
Collapse
Affiliation(s)
- Liang Dong
- Key Laboratory of Green Chemistry and Technology, Ministry of Education, College of Chemistry, Sichuan University, Chengdu, Sichuan 610064, China
| | | | | | | | | |
Collapse
|
42
|
Kang YK, Byun BJ. Assessment of density functionals with long-range and/or empirical dispersion corrections for conformational energy calculations of peptides. J Comput Chem 2010; 31:2915-23. [DOI: 10.1002/jcc.21587] [Citation(s) in RCA: 44] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
43
|
Jiang J, Wu Y, Wang ZX, Wu C. Assessing the Performance of Popular Quantum Mechanics and Molecular Mechanics Methods and Revealing the Sequence-Dependent Energetic Features Using 100 Tetrapeptide Models. J Chem Theory Comput 2010. [DOI: 10.1021/ct100008q] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Jinliang Jiang
- College of Chemistry and Chemical Engineering, Graduate University of Chinese Academy of Sciences, Beijing, 100049, China, and Department of Chemistry and Biochemistry, University of California—Santa Barbara, Santa Barbara, California 93106
| | - Yanbo Wu
- College of Chemistry and Chemical Engineering, Graduate University of Chinese Academy of Sciences, Beijing, 100049, China, and Department of Chemistry and Biochemistry, University of California—Santa Barbara, Santa Barbara, California 93106
| | - Zhi-Xiang Wang
- College of Chemistry and Chemical Engineering, Graduate University of Chinese Academy of Sciences, Beijing, 100049, China, and Department of Chemistry and Biochemistry, University of California—Santa Barbara, Santa Barbara, California 93106
| | - Chun Wu
- College of Chemistry and Chemical Engineering, Graduate University of Chinese Academy of Sciences, Beijing, 100049, China, and Department of Chemistry and Biochemistry, University of California—Santa Barbara, Santa Barbara, California 93106
| |
Collapse
|
44
|
Kaminský J, Raich I, Tomčáková K, Bouř P. Conformational behavior of simple furanosides studied by optical rotation. J Comput Chem 2010; 31:2213-24. [DOI: 10.1002/jcc.21511] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
45
|
Zhao DX, Liu C, Wang FF, Yu CY, Gong LD, Liu SB, Yang ZZ. Development of a Polarizable Force Field Using Multiple Fluctuating Charges per Atom. J Chem Theory Comput 2010; 6:795-804. [DOI: 10.1021/ct9006647] [Citation(s) in RCA: 74] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Dong-Xia Zhao
- Chemistry and Chemical Engineering Faculty, Liaoning Normal University, Dalian, 116029, China and Research Computing Center, University of North Carolina, Chapel Hill, North Carolina 27599-3420
| | - Cui Liu
- Chemistry and Chemical Engineering Faculty, Liaoning Normal University, Dalian, 116029, China and Research Computing Center, University of North Carolina, Chapel Hill, North Carolina 27599-3420
| | - Fang-Fang Wang
- Chemistry and Chemical Engineering Faculty, Liaoning Normal University, Dalian, 116029, China and Research Computing Center, University of North Carolina, Chapel Hill, North Carolina 27599-3420
| | - Chun-Yang Yu
- Chemistry and Chemical Engineering Faculty, Liaoning Normal University, Dalian, 116029, China and Research Computing Center, University of North Carolina, Chapel Hill, North Carolina 27599-3420
| | - Li-Dong Gong
- Chemistry and Chemical Engineering Faculty, Liaoning Normal University, Dalian, 116029, China and Research Computing Center, University of North Carolina, Chapel Hill, North Carolina 27599-3420
| | - Shu-Bin Liu
- Chemistry and Chemical Engineering Faculty, Liaoning Normal University, Dalian, 116029, China and Research Computing Center, University of North Carolina, Chapel Hill, North Carolina 27599-3420
| | - Zhong-Zhi Yang
- Chemistry and Chemical Engineering Faculty, Liaoning Normal University, Dalian, 116029, China and Research Computing Center, University of North Carolina, Chapel Hill, North Carolina 27599-3420
| |
Collapse
|
46
|
Mata RA. Application of high level wavefunction methods in quantum mechanics/molecular mechanics hybrid schemes. Phys Chem Chem Phys 2010; 12:5041-52. [DOI: 10.1039/b918608e] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
47
|
Semrouni D, Ohanessian G, Clavaguéra C. Structural, energetic and dynamical properties of sodiated oligoglycines: relevance of a polarizable force field. Phys Chem Chem Phys 2010; 12:3450-62. [DOI: 10.1039/b924317h] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
|
48
|
McKinney BE, Urban JJ. Fluoroolefins as Peptide Mimetics. 2. A Computational Study of the Conformational Ramifications of Peptide Bond Replacement. J Phys Chem A 2009; 114:1123-33. [DOI: 10.1021/jp9094535] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Affiliation(s)
- Brian E. McKinney
- Chemistry Department, United States Naval Academy, 572 Holloway Road, Annapolis, Maryland 21402
| | - Joseph J. Urban
- Chemistry Department, United States Naval Academy, 572 Holloway Road, Annapolis, Maryland 21402
| |
Collapse
|
49
|
Iwaoka M, Kimura N, Yosida D, Minezaki T. The SAAP force field: development of the single amino acid potentials for 20 proteinogenic amino acids and Monte Carlo molecular simulation for short peptides. J Comput Chem 2009; 30:2039-55. [PMID: 19140140 DOI: 10.1002/jcc.21196] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Molecular simulation by using force field parameters has been widely applied in the fields of peptide and protein research for various purposes. We recently proposed a new all-atom protein force field, called the SAAP force field, which utilizes single amino acid potentials (SAAPs) as the fundamental elements. In this article, whole sets of the SAAP force field parameters in vacuo, in ether, and in water have been developed by ab initio calculation for all 20 proteinogenic amino acids and applied to Monte Carlo molecular simulation for two short peptides. The side-chain separation approximation method was employed to obtain the SAAP parameters for the amino acids with a long side chain. Monte Carlo simulation for Met-enkephalin (CHO-Tyr-Gly-Gly-Phe-Met-NH2) by using the SAAP force field revealed that the conformation in vacuo is mainly controlled by strong electrostatic interactions between the amino acid residues, while the SAAPs and the interamino acid Lennard-Jones potentials are predominant in water. In ether, the conformation would be determined by the combination of the three components. On the other hand, the SAAP simulation for chignolin (H-Gly-Tyr-Asp-Pro-Glu-Thr-Gly-Thr-Trp-Gly-OH) reasonably reproduced a native-like beta-hairpin structure in water although the C-terminal and side-chain conformations were different from the native ones. It was suggested that the SAAP force field is a useful tool for analyzing conformations of polypeptides in terms of intrinsic conformational propensities of the single amino acid units.
Collapse
Affiliation(s)
- Michio Iwaoka
- Department of Chemistry, School of Science, Tokai University, Kitakaname, Hiratsuka-shi, Kanagawa 259-1292, Japan.
| | | | | | | |
Collapse
|
50
|
Handley CM, Popelier PLA. Dynamically Polarizable Water Potential Based on Multipole Moments Trained by Machine Learning. J Chem Theory Comput 2009; 5:1474-89. [DOI: 10.1021/ct800468h] [Citation(s) in RCA: 51] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Chris M. Handley
- Manchester Interdisciplinary Biocentre (MIB), 131 Princess Street, Manchester M1 7DN, Great Britain and School of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL, Great Britain
| | - Paul L. A. Popelier
- Manchester Interdisciplinary Biocentre (MIB), 131 Princess Street, Manchester M1 7DN, Great Britain and School of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL, Great Britain
| |
Collapse
|