1
|
Krevert CS, Gunkel L, Haese C, Hunger J. Ion-specific binding of cations to the carboxylate and of anions to the amide of alanylalanine. Commun Chem 2022; 5:173. [PMID: 36697920 PMCID: PMC9814750 DOI: 10.1038/s42004-022-00789-y] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2022] [Accepted: 11/30/2022] [Indexed: 12/24/2022] Open
Abstract
Studies of ion-specific effects on oligopeptides have aided our understanding of Hofmeister effects on proteins, yet the use of different model peptides and different experimental sensitivities have led to conflicting conclusions. To resolve these controversies, we study a small model peptide, L-Alanyl-L-alanine (2Ala), carrying all fundamental chemical protein motifs: C-terminus, amide bond, and N-terminus. We elucidate the effect of GdmCl, LiCl, KCl, KI, and KSCN by combining dielectric relaxation, nuclear magnetic resonance (1H-NMR), and (two-dimensional) infrared spectroscopy. Our dielectric results show that all ions reduce the rotational mobility of 2Ala, yet the magnitude of the reduction is larger for denaturing cations than for anions. The NMR chemical shifts of the amide group are particularly sensitive to denaturing anions, indicative of anion-amide interactions. Infrared experiments reveal that LiCl alters the spectral homogeneity and dynamics of the carboxylate, but not the amide group. Interaction of LiCl with the negatively charged pole of 2Ala, the COO- group, can explain the marked cationic effect on dipolar rotation, while interaction of anions between the poles, at the amide, only weakly perturbs dipolar dynamics. As such, our results provide a unifying view on ions' preferential interaction sites at 2Ala and help rationalize Hofmeister effects on proteins.
Collapse
Affiliation(s)
- Carola Sophie Krevert
- grid.419547.a0000 0001 1010 1663Department of Molecular Spectroscopy, Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Lucas Gunkel
- grid.419547.a0000 0001 1010 1663Department of Molecular Spectroscopy, Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Constantin Haese
- grid.419547.a0000 0001 1010 1663Department of Molecular Electronics, Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Johannes Hunger
- grid.419547.a0000 0001 1010 1663Department of Molecular Spectroscopy, Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| |
Collapse
|
2
|
Niemöller H, Blasius J, Hollóczki O, Kirchner B. How do alternative amino acids behave in water? A comparative ab initio molecular dynamics study of solvated α-amino acids and α-amino amidines. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120282] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
3
|
Barone V, Alessandrini S, Biczysko M, Cheeseman JR, Clary DC, McCoy AB, DiRisio RJ, Neese F, Melosso M, Puzzarini C. Computational molecular spectroscopy. ACTA ACUST UNITED AC 2021. [DOI: 10.1038/s43586-021-00034-1] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
4
|
Burley KH, Gill SC, Lim NM, Mobley DL. Enhancing Side Chain Rotamer Sampling Using Nonequilibrium Candidate Monte Carlo. J Chem Theory Comput 2019; 15:1848-1862. [PMID: 30677291 DOI: 10.1021/acs.jctc.8b01018] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Molecular simulations are a valuable tool for studying biomolecular motions and thermodynamics. However, such motions can be slow compared to simulation time scales, yet critical. Specifically, adequate sampling of side chain motions in protein binding pockets is crucial for obtaining accurate estimates of ligand binding free energies from molecular simulations. The time scale of side chain rotamer flips can range from a few ps to several hundred ns or longer, particularly in crowded environments like the interior of proteins. Here, we apply a mixed nonequilibrium candidate Monte Carlo (NCMC)/molecular dynamics (MD) method to enhance sampling of side chain rotamers. The NCMC portion of our method applies a switching protocol wherein the steric and electrostatic interactions between target side chain atoms and the surrounding environment are cycled off and then back on during the course of a move proposal. Between NCMC move proposals, simulation of the system continues via traditional molecular dynamics. Here, we first validate this approach on a simple, solvated valine-alanine dipeptide system and then apply it to a well-studied model ligand binding site in T4 lysozyme L99A. We compute the rate of rotamer transitions for a valine side chain using our approach and compare it to that of traditional molecular dynamics simulations. Here, we show that our NCMC/MD method substantially enhances side chain sampling, especially in systems where the torsional barrier to rotation is high (≥10 kcal/mol). These barriers can be intrinsic torsional barriers or steric barriers imposed by the environment. Overall, this may provide a promising strategy to selectively improve side chain sampling in molecular simulations.
Collapse
Affiliation(s)
- Kalistyn H Burley
- Department of Pharmaceutical Sciences , University of California, Irvine , Irvine , California 92697 , United States
| | - Samuel C Gill
- Department of Chemistry , University of California, Irvine , Irvine , California 92617 , United States
| | - Nathan M Lim
- Department of Pharmaceutical Sciences , University of California, Irvine , Irvine , California 92697 , United States
| | - David L Mobley
- Department of Pharmaceutical Sciences , University of California, Irvine , Irvine , California 92697 , United States.,Department of Chemistry , University of California, Irvine , Irvine , California 92617 , United States
| |
Collapse
|
5
|
Mironov V, Alexeev Y, Mulligan VK, Fedorov DG. A systematic study of minima in alanine dipeptide. J Comput Chem 2018; 40:297-309. [DOI: 10.1002/jcc.25589] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2018] [Revised: 07/12/2018] [Accepted: 08/07/2018] [Indexed: 12/11/2022]
Affiliation(s)
- Vladimir Mironov
- Department of Chemistry Lomonosov Moscow State University Leninskie Gory 1/3, Moscow 119991 Russia
| | - Yuri Alexeev
- Argonne National Laboratory Computational Science Division Argonne Illinois 60439
| | - Vikram Khipple Mulligan
- Department of Biochemistry University of Washington, Institute for Protein Design Seattle Washington 98195
| | - Dmitri G. Fedorov
- CD‐FMat National Institute of Advanced Industrial Science and Technology Central 2, Umezono 1‐1‐1, Tsukuba 305‐8568 Japan
| |
Collapse
|
6
|
Cormanich RA, Bühl M, Rittner R. Understanding the conformational behaviour of Ac-Ala-NHMe in different media. A joint NMR and DFT study. Org Biomol Chem 2015. [PMID: 26219244 DOI: 10.1039/c5ob01296a] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
The conformational behaviour of Ac-Ala-NHMe was studied in the gas-phase and in solution by theoretical calculations (B3LYP-D3/aug-cc-pVDZ level) and experimental (1)H NMR. The conformational preferences of this compound were shown to result from a complex interplay between the strengths of possible intramolecular hydrogen bonds, steric interactions, hyperconjugation, entropy effects and the overall dipole moments. The Ac-Ala-N(Me)2 derivative was studied in addition, to design a system akin to Ac-Ala-NHMe, but with disrupted intramolecular hydrogen bonds involving the -NHMe group, mimicking the effect of polar protic solvents.
Collapse
Affiliation(s)
- Rodrigo A Cormanich
- Chemistry Institute, University of Campinas, Campinas, SP 13083-970, Brazil.
| | | | | |
Collapse
|
7
|
Lanza G, Chiacchio MA. Ab Initio MP2 and Density Functional Theory Computational Study of AcAlaNH2Peptide Hydration: A Bottom-Up Approach. Chemphyschem 2014; 15:2785-93. [DOI: 10.1002/cphc.201402222] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2014] [Indexed: 01/15/2023]
|
8
|
Lindner B, Yi Z, Prinz JH, Smith JC, Noé F. Dynamic neutron scattering from conformational dynamics. I. Theory and Markov models. J Chem Phys 2013; 139:175101. [DOI: 10.1063/1.4824070] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
|
9
|
Beck JP, Gaigeot MP, Lisy JM. Anharmonic vibrations of N-H in Cl(-)(N-methylacetamide)1(H2O)(0-2)Ar2 cluster ions. Combined IRPD experiments and BOMD simulations. Phys Chem Chem Phys 2013; 15:16736-45. [PMID: 23986352 DOI: 10.1039/c3cp52418c] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Infrared Predissociation (IRPD) spectra of Cl(-)(NMA)1(H2O)0-2Ar2 combined with Born-Oppenheimer Molecular Dynamics (BOMD) IR spectra have been acquired, providing the structure and dynamics of these systems. We show that the chloride ion is bound to the hydrogen of the amide N-H group, forming a strong ionic hydrogen bond, weakening the N-H stretch, and shifting it to lower frequency. The presence of water molecules enhances the ionic hydrogen bond by binding to the amide carbonyl oxygen of NMA and shifts the N-H stretch further to lower frequency. The BOMD IR spectra can recapture all, but about 100 cm(-1), of the 600 to 700 cm(-1) shifts due to the strong N-H stretch anharmonicities observed in experiments. This residual error was found to be due to the lack of zero point energy in the classical treatment of motion in the BOMD method.
Collapse
Affiliation(s)
- Jordan P Beck
- Concordia University Wisconsin, 12800 N. Lakeshore Drive, Mequon, Wisconsin 53097, USA
| | | | | |
Collapse
|
10
|
Jiang F, Han W, Wu YD. The intrinsic conformational features of amino acids from a protein coil library and their applications in force field development. Phys Chem Chem Phys 2013; 15:3413-28. [PMID: 23385383 DOI: 10.1039/c2cp43633g] [Citation(s) in RCA: 52] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The local conformational (φ, ψ, χ) preferences of amino acid residues remain an active research area, which are important for the development of protein force fields. In this perspective article, we first summarize spectroscopic studies of alanine-based short peptides in aqueous solution. While most studies indicate a preference for the P(II) conformation in the unfolded state over α and β conformations, significant variations are also observed. A statistical analysis from various coil libraries of high-resolution protein structures is then summarized, which gives a more coherent view of the local conformational features. The φ, ψ, χ distributions of the 20 amino acids have been obtained from a protein coil library, considering both backbone and side-chain conformational preferences. The intrinsic side-chain χ(1) rotamer preference and χ(1)-dependent Ramachandran plot can be generally understood by combining the interaction of the side-chain Cγ/Oγ atom with two neighboring backbone peptide groups. Current all-atom force fields such as AMBER ff99sb-ILDN, ff03 and OPLS-AA/L do not reproduce these distributions well. A method has been developed by combining the φ, ψ plot of alanine with the influence of side-chain χ(1) rotamers to derive the local conformational features of various amino acids. It has been further applied to improve the OPLS-AA force field. The modified force field (OPLS-AA/C) reproduces experimental (3)J coupling constants for various short peptides quite well. It also better reproduces the temperature-dependence of the helix-coil transition for alanine-based peptides. The new force field can fold a series of peptides and proteins with various secondary structures to their experimental structures. MD simulations of several globular proteins using the improved force field give significantly less deviation (RMSD) to experimental structures. The results indicate that the local conformational features from coil libraries are valuable for the development of balanced protein force fields.
Collapse
Affiliation(s)
- Fan Jiang
- Laboratory of Computational Chemistry and Drug Design, Laboratory of Chemical Genomics, Peking University Shenzhen Graduate School, Shenzhen 518055, China
| | | | | |
Collapse
|
11
|
Lanza G, Chiacchio MA. Comprehensive and Accurate Ab Initio Energy Surface of Simple Alanine Peptides. Chemphyschem 2013; 14:3284-93. [DOI: 10.1002/cphc.201300445] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2013] [Indexed: 11/07/2022]
Affiliation(s)
- Giuseppe Lanza
- Dipartimento di Scienze del Farmaco, Università di Catania, Viale A. Doria 6, 95125 Catania (Italy)
| | - Maria A. Chiacchio
- Dipartimento di Scienze del Farmaco, Università di Catania, Viale A. Doria 6, 95125 Catania (Italy)
| |
Collapse
|
12
|
Li Y, Gao Y, Zhang X, Wang X, Mou L, Duan L, He X, Mei Y, Zhang JZH. A coupled two-dimensional main chain torsional potential for protein dynamics: generation and implementation. J Mol Model 2013; 19:3647-57. [PMID: 23765039 DOI: 10.1007/s00894-013-1879-8] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2013] [Accepted: 05/01/2013] [Indexed: 11/29/2022]
Abstract
Main chain torsions of alanine dipeptide are parameterized into coupled 2-dimensional Fourier expansions based on quantum mechanical (QM) calculations at M06 2X/aug-cc-pvtz//HF/6-31G** level. Solvation effect is considered by employing polarizable continuum model. Utilization of the M06 2X functional leads to precise potential energy surface that is comparable to or even better than MP2 level, but with much less computational demand. Parameterization of the 2D expansions is against the full main chain torsion space instead of just a few low energy conformations. This procedure is similar to that for the development of AMBER03 force field, except unique weighting factor was assigned to all the grid points. To avoid inconsistency between quantum mechanical calculations and molecular modeling, the model peptide is further optimized at molecular mechanics level with main chain dihedral angles fixed before the calculation of the conformational energy on molecular mechanical level at each grid point, during which generalized Born model is employed. Difference in solvation models at quantum mechanics and molecular mechanics levels makes this parameterization procedure less straightforward. All force field parameters other than main chain torsions are taken from existing AMBER force field. With this new main chain torsion terms, we have studied the main chain dihedral distributions of ALA dipeptide and pentapeptide in aqueous solution. The results demonstrate that 2D main chain torsion is effective in delineating the energy variation associated with rotations along main chain dihedrals. This work is an implication for the necessity of more accurate description of main chain torsions in the future development of ab initio force field and it also raises a challenge to the development of quantum mechanical methods, especially the quantum mechanical solvation models.
Collapse
Affiliation(s)
- Yongxiu Li
- Center for Laser and Computational Biophysics, State Key Laboratory of Precision Spectroscopy and Department of Physics and Institute of Theoretical and Computational Science, East China Normal University, Shanghai, 200062, China
| | | | | | | | | | | | | | | | | |
Collapse
|
13
|
Simon A, Spiegelman F. Water clusters adsorbed on polycyclic aromatic hydrocarbons: Energetics and conformational dynamics. J Chem Phys 2013; 138:194309. [DOI: 10.1063/1.4805015] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
14
|
Abstract
Transitions between metastable conformations of a dipeptide are investigated using classical molecular dynamics simulation with explicit water molecules. The distribution of the surrounding water at different moments before the transitions and the dynamical correlations of water with the peptide's configurational motions indicate that the water molecules represent an integral part of the molecular system during the conformational changes, in contrast with the metastable periods when water and peptide dynamics are essentially decoupled.
Collapse
Affiliation(s)
- Dmitry Nerukh
- †Nonlinearity and Complexity Research Group, Aston University, Birmingham, B4 7ET, United Kingdom
| | - Sergey Karabasov
- ‡School of Engineering and Materials Science, Queen Mary University of London, Mile End Road, London E1 4NS, United Kingdom
| |
Collapse
|
15
|
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
|
16
|
Okamoto T, Ishikawa T, Koyano Y, Yamamoto N, Kuwata K, Nagaoka M. A Minimal Implementation of the AMBER-PAICS Interface for Ab Initio FMO-QM/MM-MD Simulation. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2013. [DOI: 10.1246/bcsj.20120216] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Takuya Okamoto
- Graduate School of Information Science, Nagoya University
| | - Takeshi Ishikawa
- Division of Prion Research, Center for Emerging Infectious Disease, Gifu University
| | | | | | - Kazuo Kuwata
- Division of Prion Research, Center for Emerging Infectious Disease, Gifu University
| | | |
Collapse
|
17
|
Mirkin NG, Krimm S. Water interaction differences determine the relative energetic stability of the polyproline II conformation of the alanine dipeptide in aqueous environments. Biopolymers 2012; 97:789-94. [DOI: 10.1002/bip.22064] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
|
18
|
Han C, Wang J. Influence of an Unnatural Amino Acid Side Chain on the Conformational Dynamics of Peptides. Chemphyschem 2012; 13:1522-34. [DOI: 10.1002/cphc.201100995] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2011] [Indexed: 11/09/2022]
|
19
|
Jin L, Gerber RB. Migration and chemical reaction of H+ in protonated β-galactose. Phys Chem Chem Phys 2012; 14:13522-6. [DOI: 10.1039/c2cp41766a] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
20
|
Ingrosso F, Monard G, Hamdi Farag M, Bastida A, Ruiz-López MF. Importance of Polarization and Charge Transfer Effects to Model the Infrared Spectra of Peptides in Solution. J Chem Theory Comput 2011; 7:1840-9. [DOI: 10.1021/ct2000588] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Francesca Ingrosso
- Equipe de Chimie et Biochimie Théoriques, UMR 7565 SRSMC, CNRS-Nancy Université, BP 70239 Vandœuvre-lès Nancy, France
| | - Gérald Monard
- Equipe de Chimie et Biochimie Théoriques, UMR 7565 SRSMC, CNRS-Nancy Université, BP 70239 Vandœuvre-lès Nancy, France
| | - Marwa Hamdi Farag
- Departamento de Química Física, Universidad de Murcia, 30100 Murcia, Spain
| | - Adolfo Bastida
- Departamento de Química Física, Universidad de Murcia, 30100 Murcia, Spain
| | - Manuel F. Ruiz-López
- Equipe de Chimie et Biochimie Théoriques, UMR 7565 SRSMC, CNRS-Nancy Université, BP 70239 Vandœuvre-lès Nancy, France
| |
Collapse
|
21
|
Biswal HS, Loquais Y, Tardivel B, Gloaguen E, Mons M. Isolated monohydrates of a model peptide chain: effect of a first water molecule on the secondary structure of a capped phenylalanine. J Am Chem Soc 2011; 133:3931-42. [PMID: 21361380 DOI: 10.1021/ja108643p] [Citation(s) in RCA: 61] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The formation of monohydrates of capped phenylalanine model peptides, CH(3)-CO-Phe-NH(2) and CH(3)-CO-Phe-NH-CH(3), in a supersonic expansion has been investigated using laser spectroscopy and quantum chemistry methods. Conformational distributions of the monohydrates have been revealed by IR/UV double-resonance spectroscopy and their structures assigned by comparison with DFT-D calculations. A careful analysis of the final hydrate distribution together with a detailed theoretical investigation of the potential energy surface of the monohydrates demonstrates that solvation occurs from the conformational distribution of the isolated peptide monomers. The distribution of the monohydrates appears to be strongly dependent on both the initial monomer conformation (extended or folded backbone) and the solvation site initially occupied by the water molecule. The solvation processes taking place during the cooling can be categorized as follows: (a) solvation without significant structural changes of the peptide, (b) solvation inducing significant distortions of the backbone but retaining the secondary structure, and (c) solvation triggering backbone isomerizations, leading to a modification of the peptide secondary structure. It is observed that solvation by a single water molecule can fold a β-strand into a γ-turn structure (type c) or induce a significant opening of a γ-turn characterized by an elongated C(7) hydrogen bond (type b). These structural changes can be considered as a first step toward the polyproline II condensed-phase structure, illustrating the role played by the very first water molecule in the solvation process.
Collapse
Affiliation(s)
- Himansu S Biswal
- Laboratoire Francis Perrin, CEA/DSM/IRAMIS/SPAM-CNRS URA 2453, CEA/Saclay, 91191 Gif-sur-Yvette, France
| | | | | | | | | |
Collapse
|
22
|
Gaigeot MP, Besley NA, Hirst JD. Modeling the infrared and circular dichroism spectroscopy of a bridged cyclic diamide. J Phys Chem B 2011; 115:5526-35. [PMID: 21344909 DOI: 10.1021/jp111140f] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Density functional theory based molecular dynamics simulations are used to study the structure, infrared (IR) spectroscopy, circular dichroism (CD) spectroscopy, and coupling between the amide I vibrations of a bridged cyclic diamide in the gas phase and in aqueous solution. IR spectra computed via the dipole moment time correlation function show a large red-shift of 30 cm(-1) in the amide I vibration in solution compared to the gas phase, and are in good agreement with experiment. Conformationally averaged CD spectra computed using the CIS(D) method are highly sensitive to the structures used, and structures sampled in the aqueous phase simulation are required to obtain qualitatively correct CD spectra. Analysis of the coupling between the amide I modes shows that in the aqueous phase there is an increased localization of the vibrations on the individual peptide groups and a reduction in the mode coupling parameter compared to the gas phase. Overall, the results illustrate the significance of incorporating molecular dynamics in the simulation of IR and CD spectra.
Collapse
Affiliation(s)
- Marie-Pierre Gaigeot
- Université d'Evry Val d'Essonne, LAMBE UMR8587 Laboratoire Analyse et Modélisation pour la Biologie et l'Environment, Blvd F. Mitterrand, Bât. Maupertuis, 91025 Evry, France
| | | | | |
Collapse
|
23
|
Liu Z, Ensing B, Moore PB. Quantitative Assessment of Force Fields on Both Low-Energy Conformational Basins and Transition-State Regions of the (ϕ-ψ) Space. J Chem Theory Comput 2010; 7:402-19. [PMID: 26596162 DOI: 10.1021/ct100395n] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
The free energy surfaces (FESs) of alanine dipeptide are studied to illustrate a new strategy to assess the performance of classical molecular mechanics force field on the full range of the (ϕ-ψ) conformational space. The FES is obtained from metadynamics simulations with five commonly used force fields and from ab initio density functional theory calculations in both gas phase and aqueous solution. The FESs obtained at the B3LYP/6-311+G(2d,p)//B3LYP/6-31G(d,p) level of theory are validated by comparison with previously reported MP2 and LMP2 results as well as with experimentally obtained probability distribution between the C5-β (or β-PPII) and αR states. A quantitative assessment is made for each force field in three conformational basins, LeRI (C5-β-C7eq), LeRII (β2-αR), and LeRIII(αL-C7ax-αD) as well as three transition-state regions linking the above conformational basins. The performance of each force field is evaluated in terms of the average free energy of each region in comparison with that of the ab initio results. We quantify how well a force field FES matches the ab initio FES through the calculation of the standard deviation of a free energy difference map between the two FESs. The results indicate that the performance varies largely from region to region or from force field to force field. Although not one force field is able to outperform all others in all conformational areas, the OPLSAA/L force field gives the best performance overall, followed by OPLSAA and AMBER03. For the three top performers, the average free energies differ from the corresponding ab initio values from within the error range (<0.4 kcal/mol) to ∼1.5 kcal/mol for the low-energy regions and up to ∼2.0 kcal/mol for the transition-state regions. The strategy presented and the results obtained here should be useful for improving the parametrization of force fields targeting both accuracy in the energies of conformers and the transition-state barriers.
Collapse
Affiliation(s)
- Zhiwei Liu
- West Center for Computational Chemistry and Drug Design, Department of Chemistry & Biochemistry, University of the Sciences in Philadelphia, 600 South 43rd Street, Philadelphia, Pennsylvania 19104, United States and Van 't Hoff Institute for Molecular Sciences, University of Amsterdam, Nieuwe Achtergracht 166, 1018 WV Amsterdam, The Netherlands
| | - Bernd Ensing
- West Center for Computational Chemistry and Drug Design, Department of Chemistry & Biochemistry, University of the Sciences in Philadelphia, 600 South 43rd Street, Philadelphia, Pennsylvania 19104, United States and Van 't Hoff Institute for Molecular Sciences, University of Amsterdam, Nieuwe Achtergracht 166, 1018 WV Amsterdam, The Netherlands
| | - Preston B Moore
- West Center for Computational Chemistry and Drug Design, Department of Chemistry & Biochemistry, University of the Sciences in Philadelphia, 600 South 43rd Street, Philadelphia, Pennsylvania 19104, United States and Van 't Hoff Institute for Molecular Sciences, University of Amsterdam, Nieuwe Achtergracht 166, 1018 WV Amsterdam, The Netherlands
| |
Collapse
|
24
|
Gaigeot MP, Cimas A, Seydou M, Kim JY, Lee S, Schermann JP. Proton Transfer from the Inactive Gas-Phase Nicotine Structure to the Bioactive Aqueous-Phase Structure. J Am Chem Soc 2010; 132:18067-77. [DOI: 10.1021/ja103759v] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Marie-Pierre Gaigeot
- Laboratoire Analyse et Modélisation pour la Biologie et l’Environnement, LAMBE UMR8587 CNRS, Université d’Evry val d’Essonne, boulevard F. Mitterrand, Bat. Maupertuis, 91025 Evry Cedex, France, Institut Universitaire de France IUF, 103 boulevard St Michel, 75005 Paris, France, Subatech, Ecole des Mines de Nantes, 4 rue Alfred Kastler, 44307 Nantes, France, Department of Applied Chemistry, Kyunghee University, Kyungki 446-701, South Korea, Department of Biophysics and Biochemical Chemistry, WCU, Seoul
| | - Alvaro Cimas
- Laboratoire Analyse et Modélisation pour la Biologie et l’Environnement, LAMBE UMR8587 CNRS, Université d’Evry val d’Essonne, boulevard F. Mitterrand, Bat. Maupertuis, 91025 Evry Cedex, France, Institut Universitaire de France IUF, 103 boulevard St Michel, 75005 Paris, France, Subatech, Ecole des Mines de Nantes, 4 rue Alfred Kastler, 44307 Nantes, France, Department of Applied Chemistry, Kyunghee University, Kyungki 446-701, South Korea, Department of Biophysics and Biochemical Chemistry, WCU, Seoul
| | - Mahamadou Seydou
- Laboratoire Analyse et Modélisation pour la Biologie et l’Environnement, LAMBE UMR8587 CNRS, Université d’Evry val d’Essonne, boulevard F. Mitterrand, Bat. Maupertuis, 91025 Evry Cedex, France, Institut Universitaire de France IUF, 103 boulevard St Michel, 75005 Paris, France, Subatech, Ecole des Mines de Nantes, 4 rue Alfred Kastler, 44307 Nantes, France, Department of Applied Chemistry, Kyunghee University, Kyungki 446-701, South Korea, Department of Biophysics and Biochemical Chemistry, WCU, Seoul
| | - Ju-Young Kim
- Laboratoire Analyse et Modélisation pour la Biologie et l’Environnement, LAMBE UMR8587 CNRS, Université d’Evry val d’Essonne, boulevard F. Mitterrand, Bat. Maupertuis, 91025 Evry Cedex, France, Institut Universitaire de France IUF, 103 boulevard St Michel, 75005 Paris, France, Subatech, Ecole des Mines de Nantes, 4 rue Alfred Kastler, 44307 Nantes, France, Department of Applied Chemistry, Kyunghee University, Kyungki 446-701, South Korea, Department of Biophysics and Biochemical Chemistry, WCU, Seoul
| | - Sungyul Lee
- Laboratoire Analyse et Modélisation pour la Biologie et l’Environnement, LAMBE UMR8587 CNRS, Université d’Evry val d’Essonne, boulevard F. Mitterrand, Bat. Maupertuis, 91025 Evry Cedex, France, Institut Universitaire de France IUF, 103 boulevard St Michel, 75005 Paris, France, Subatech, Ecole des Mines de Nantes, 4 rue Alfred Kastler, 44307 Nantes, France, Department of Applied Chemistry, Kyunghee University, Kyungki 446-701, South Korea, Department of Biophysics and Biochemical Chemistry, WCU, Seoul
| | - Jean-Pierre Schermann
- Laboratoire Analyse et Modélisation pour la Biologie et l’Environnement, LAMBE UMR8587 CNRS, Université d’Evry val d’Essonne, boulevard F. Mitterrand, Bat. Maupertuis, 91025 Evry Cedex, France, Institut Universitaire de France IUF, 103 boulevard St Michel, 75005 Paris, France, Subatech, Ecole des Mines de Nantes, 4 rue Alfred Kastler, 44307 Nantes, France, Department of Applied Chemistry, Kyunghee University, Kyungki 446-701, South Korea, Department of Biophysics and Biochemical Chemistry, WCU, Seoul
| |
Collapse
|
25
|
Oh KI, Kim W, Joo C, Yoo DG, Han H, Hwang GS, Cho M. Azido Gauche Effect on the Backbone Conformation of β-Azidoalanine Peptides. J Phys Chem B 2010; 114:13021-9. [DOI: 10.1021/jp107359m] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Kwang-Im Oh
- Department of Chemistry, Korea University, Seoul 136-701, Korea, and Korea Basic Science Institute, Seoul 136-713, Korea
| | - Woosung Kim
- Department of Chemistry, Korea University, Seoul 136-701, Korea, and Korea Basic Science Institute, Seoul 136-713, Korea
| | - Cheonik Joo
- Department of Chemistry, Korea University, Seoul 136-701, Korea, and Korea Basic Science Institute, Seoul 136-713, Korea
| | - Dong-Geun Yoo
- Department of Chemistry, Korea University, Seoul 136-701, Korea, and Korea Basic Science Institute, Seoul 136-713, Korea
| | - Hogyu Han
- Department of Chemistry, Korea University, Seoul 136-701, Korea, and Korea Basic Science Institute, Seoul 136-713, Korea
| | - Geum-Sook Hwang
- Department of Chemistry, Korea University, Seoul 136-701, Korea, and Korea Basic Science Institute, Seoul 136-713, Korea
| | - Minhaeng Cho
- Department of Chemistry, Korea University, Seoul 136-701, Korea, and Korea Basic Science Institute, Seoul 136-713, Korea
| |
Collapse
|
26
|
Gaigeot MP. Infrared spectroscopy of the alanine dipeptide analog in liquid water with DFT-MD. Direct evidence for P(II)/beta conformations. Phys Chem Chem Phys 2010; 12:10198-209. [PMID: 20539891 DOI: 10.1039/c003485a] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Following our previous work [J. Phys. Chem. B. Lett., 2009, 113, 10059], DFT-based molecular dynamics (DFTMD) simulations of 2-Ala peptide (i.e. Ac-Ala-NHMe dialanine peptide analog with methyl group caps at the extremities) immersed in liquid water at room temperature are reported. Our goal here is the theoretical calculation of the infrared spectrum of aqueous 2-Ala, in order to provide a definitive understanding of the average conformation adopted by this peptide in the liquid phase, taking into account solute and solvent at the same theoretical level of representation. We find that the experimental Amide I-II band predominantly results from a mixture of partially unfolded P(II) and unfolded beta conformational equilibrium of aqueous 2-Ala at room temperature.
Collapse
Affiliation(s)
- Marie-Pierre Gaigeot
- Université d'Evry val d'Essonne, LAMBE UMR8587 Laboratoire Analyse et Modélisation pour la Biologie et l'Environnement, Blvd F. Mitterrand, Bât. Maupertuis, 91025 Evry, France.
| |
Collapse
|
27
|
Gaigeot MP. Theoretical spectroscopy of floppy peptides at room temperature. A DFTMD perspective: gas and aqueous phase. Phys Chem Chem Phys 2010; 12:3336-59. [PMID: 20336243 DOI: 10.1039/b924048a] [Citation(s) in RCA: 135] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Theoretical spectroscopy is mandatory for a precise understanding and assignment of experimental spectra recorded at finite temperature. We review here room temperature DFT-based molecular dynamics simulations for the purpose of interpreting finite temperature infrared spectra of peptides of increasing size and complexity, in terms of temperature-dependent conformational dynamics and flexibility, and vibrational anharmonicities (potential energy surface anharmonicities, vibrational mode couplings and dipole anharmonicities). We take examples from our research projects in order to illustrate the main key-points and strengths of dynamical spectra modeling in that context. The calculations are presented in relation to room temperature gas phase IR-MPD experiments and room temperature liquid phase IR absorption experiments. These illustrations of floppy polypeptides have been chosen in order to convey the following ideas: temperature-dependent spectra modeling is pivotal for a precise understanding of gas phase spectra recorded at room temperature, including conformational dynamics and vibrational anharmonicities; harmonic spectroscopy (as commonly performed in the literature) can be misleading and even erroneous for a proper interpretation of spectra recorded at finite temperature; taking into account vibrational anharmonicities is pivotal for a proper interplay between theory and experiments; amide I-III bands are not necessarily the most relevant fingerprints for unraveling the local structures of peptides and more complex systems; liquid phase simulations have unraveled relationships between the zwitterionic properties of the peptide bonds and infrared signatures. The review presents a state-of-the-art account of the domain and offers perspectives and new developments for future still more challenging applications.
Collapse
Affiliation(s)
- Marie-Pierre Gaigeot
- Université d'Evry val d'Essonne, LAMBE UMR8587 Laboratoire Analyse et Modélisation pour la Biologie et l'Environnement, Blvd F. Mitterrand, Bat Maupertuis, 91025 Evry, France.
| |
Collapse
|
28
|
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]
|