1
|
Are all-atom any better than united-atom force fields for the description of liquid properties of alkanes? 2. A systematic study considering different chain lengths. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.118829] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
2
|
Zhao L, Damodaran S. Hofmeister Order of Anions on Protein Stability Originates from Lifshitz-van der Waals Dispersion Interaction with the Protein Phase. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2019; 35:12993-13002. [PMID: 31512478 DOI: 10.1021/acs.langmuir.9b00486] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
The mechanism underpinning the Hofmeister order of anions on protein stability and other physical and biological processes has been a mystery since its discovery in 1888. In the present study, we investigated electrostatic and Lifshitz-van der Waals (L-vdW) dispersion (electrodynamic) interactions between Hofmeister salts and four monomeric globular proteins. It is shown that structure-stabilizing salts exerted positive L-vdW pressure, whereas structure-destabilizing salts exerted negative L-vdW pressure on proteins. The relative order of the L-vdW pressure followed the Hofmeister series and it overshadowed the electrostatic pressure at high salt concentrations. The net change in the thermal denaturation temperature (ΔTd) of proteins in 0.8 M Hofmeister salt solutions followed a linear relationship (r2 > 0.8) with the net electrodynamic pressure regardless of the physicochemical differences between proteins. This study also revealed that segregation of anions into structure stabilizers and destabilizers depended on the dielectric susceptibility of the anion in the ultraviolet region: ions having absorbtion spectrum in the ultraviolet region (e.g., Cl-, Br-, I-, and SCN-) exerted a negative electrodynamic pressure, whereas those with absorbtion spectrum only in the infrared region, for example, SO42-, exerted a positive electrodynamic pressure. The lack of ultraviolet absorption of SO42- ions was because of quenching of ultraviolet radiation by water at below 170 nm.
Collapse
Affiliation(s)
- Lei Zhao
- Department of Food Science, College of Agriculture and Life Sciences , University of Wisconsin-Madison , Madison 53706 , Wisconsin , United States
| | - Srinivasan Damodaran
- Department of Food Science, College of Agriculture and Life Sciences , University of Wisconsin-Madison , Madison 53706 , Wisconsin , United States
| |
Collapse
|
3
|
Zhao C, Caplan DA, Noskov SY. Evaluations of the Absolute and Relative Free Energies for Antidepressant Binding to the Amino Acid Membrane Transporter LeuT with Free Energy Simulations. J Chem Theory Comput 2015; 6:1900-14. [PMID: 26615849 DOI: 10.1021/ct9006597] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The binding of ligands to protein receptors with high affinity and specificity is central to many cellular processes. The quest for the development of computational models capable of accurately evaluating binding affinity remains one of the main goals of modern computational biophysics. In this work, free energy perturbation/molecular dynamics simulations were used to evaluate absolute and relative binding affinity for three different antidepressants to a sodium-dependent membrane transporter, LeuT, a bacterial homologue of human serotonin and dopamine transporters. Dysfunction of these membrane transporters in mammals has been implicated in multiple diseases of the nervous system, including bipolar disorder and depression. Furthermore, these proteins are key targets for antidepressants including fluoxetine (aka Prozac) and tricyclic antidepressants known to block transport activity. In addition to being clinically relevant, this system, where multiple crystal structures are readily available, represents an ideal testing ground for methods used to study the molecular mechanisms of ligand binding to membrane proteins. We discuss possible pitfalls and different levels of approximation required to evaluate binding affinity, such as the dependence of the computed affinities on the strength of constraints and the sensitivity of the computed affinities to the particular partial charges derived from restrained electrostatic potential fitting of quantum mechanics electrostatic potential. Finally, we compare the effects of different constraint schemes on the absolute and relative binding affinities obtained from free energy simulations.
Collapse
Affiliation(s)
- Chunfeng Zhao
- Institute for Biocomplexity and Informatics and Department of Biological Sciences, University of Calgary, 2500 University Drive, BI558, Calgary, AB, Canada T2N 1N4 and Molecular Structure and Function, Hospital for Sick Children and Department of Biochemistry, University of Toronto, Ontario, Canada
| | - David A Caplan
- Institute for Biocomplexity and Informatics and Department of Biological Sciences, University of Calgary, 2500 University Drive, BI558, Calgary, AB, Canada T2N 1N4 and Molecular Structure and Function, Hospital for Sick Children and Department of Biochemistry, University of Toronto, Ontario, Canada
| | - Sergei Yu Noskov
- Institute for Biocomplexity and Informatics and Department of Biological Sciences, University of Calgary, 2500 University Drive, BI558, Calgary, AB, Canada T2N 1N4 and Molecular Structure and Function, Hospital for Sick Children and Department of Biochemistry, University of Toronto, Ontario, Canada
| |
Collapse
|
4
|
Matsumura Y, Sato H. An integral equation theory for solvation effects on the molecular structural fluctuation. J Chem Phys 2015; 143:014104. [DOI: 10.1063/1.4923038] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yoshihiro Matsumura
- Department of Molecular Engineering, Kyoto University, Nishikyo-ku, Kyoto 615-8510, Japan
| | - Hirofumi Sato
- Department of Molecular Engineering, Kyoto University, Nishikyo-ku, Kyoto 615-8510, Japan
- Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Nishikyo-ku, Kyoto 615-8510, Japan
| |
Collapse
|
5
|
Sonti R, Rai R, Ragothama S, Balaram P. NMR Analysis of Cross Strand Aromatic Interactions in an 8 Residue Hairpin and a 14 Residue Three Stranded β-Sheet Peptide. J Phys Chem B 2012; 116:14207-15. [DOI: 10.1021/jp3034769] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Rajesh Sonti
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore, India 560012
- NMR Research Centre, Indian Institute of Science, Bangalore, India 560012
| | - Rajkishor Rai
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore, India 560012
| | | | - Padmanabhan Balaram
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore, India 560012
| |
Collapse
|
6
|
Brooks B, Brooks C, MacKerell A, Nilsson L, Petrella R, Roux B, Won Y, Archontis G, Bartels C, Boresch S, Caflisch A, Caves L, Cui Q, Dinner A, Feig M, Fischer S, Gao J, Hodoscek M, Im W, Kuczera K, Lazaridis T, Ma J, Ovchinnikov V, Paci E, Pastor R, Post C, Pu J, Schaefer M, Tidor B, Venable RM, Woodcock HL, Wu X, Yang W, York D, Karplus M. CHARMM: the biomolecular simulation program. J Comput Chem 2009; 30:1545-614. [PMID: 19444816 PMCID: PMC2810661 DOI: 10.1002/jcc.21287] [Citation(s) in RCA: 6096] [Impact Index Per Article: 406.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
CHARMM (Chemistry at HARvard Molecular Mechanics) is a highly versatile and widely used molecular simulation program. It has been developed over the last three decades with a primary focus on molecules of biological interest, including proteins, peptides, lipids, nucleic acids, carbohydrates, and small molecule ligands, as they occur in solution, crystals, and membrane environments. For the study of such systems, the program provides a large suite of computational tools that include numerous conformational and path sampling methods, free energy estimators, molecular minimization, dynamics, and analysis techniques, and model-building capabilities. The CHARMM program is applicable to problems involving a much broader class of many-particle systems. Calculations with CHARMM can be performed using a number of different energy functions and models, from mixed quantum mechanical-molecular mechanical force fields, to all-atom classical potential energy functions with explicit solvent and various boundary conditions, to implicit solvent and membrane models. The program has been ported to numerous platforms in both serial and parallel architectures. This article provides an overview of the program as it exists today with an emphasis on developments since the publication of the original CHARMM article in 1983.
Collapse
Affiliation(s)
- B.R. Brooks
- Laboratory of Computational Biology, National Heart, Lung, and
Blood Institute, National Institutes of Health, Bethesda, MD 20892
| | - C.L. Brooks
- Departments of Chemistry & Biophysics, University of
Michigan, Ann Arbor, MI 48109
| | - A.D. MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy,
University of Maryland, Baltimore, MD, 21201
| | - L. Nilsson
- Karolinska Institutet, Department of Biosciences and Nutrition,
SE-141 57, Huddinge, Sweden
| | - R.J. Petrella
- Department of Chemistry and Chemical Biology, Harvard University,
Cambridge, MA 02138
- Department of Medicine, Harvard Medical School, Boston, MA
02115
| | - B. Roux
- Department of Biochemistry and Molecular Biology, University of
Chicago, Gordon Center for Integrative Science, Chicago, IL 60637
| | - Y. Won
- Department of Chemistry, Hanyang University, Seoul
133–792 Korea
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | - M. Karplus
- Department of Chemistry and Chemical Biology, Harvard University,
Cambridge, MA 02138
- Laboratoire de Chimie Biophysique, ISIS, Université de
Strasbourg, 67000 Strasbourg France
| |
Collapse
|
7
|
Kasezawa K, Kato M. Effect of Pressure on Conformational Equilibria of 1-Chloropropane and 1-Bromopropane in Water and Organic Solvents: A Raman Spectroscopic Study. J Phys Chem B 2009; 113:8607-12. [DOI: 10.1021/jp900073p] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
- Kunihiro Kasezawa
- Graduate School of Science and Engineering, and Department of Pharmacy, College of Pharmaceutical Sciences, Ritsumeikan University, 1-1-1 Nojihigashi, Kusatsu, Shiga 525-8577, Japan
| | - Minoru Kato
- Graduate School of Science and Engineering, and Department of Pharmacy, College of Pharmaceutical Sciences, Ritsumeikan University, 1-1-1 Nojihigashi, Kusatsu, Shiga 525-8577, Japan
| |
Collapse
|
8
|
Morozov AN, Lin SH. Thermodynamics of a conformational change using a random walk in energy-reaction coordinate space: Application to methane dimer hydrophobic interactions. J Chem Phys 2009; 130:074903. [PMID: 19239312 DOI: 10.1063/1.3077658] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A random walk sampling algorithm allows the extraction of the density of states distribution in energy-reaction coordinate space. As a result, the temperature dependences of thermodynamic quantities such as relative energy, entropy, and heat capacity can be calculated using first-principles statistical mechanics. The strategies for optimal convergence of the algorithm and control of its accuracy are proposed. We show that the saturation of the error [Q. Yan and J. J. de Pablo, Phys. Rev. Lett. 90, 035701 (2003); E. Belardinelli and V. D. Pereyra, J. Chem. Phys. 127, 184105 (2007)] is due to the use of histogram flatness as a criterion of convergence. An application of the algorithm to methane dimer hydrophobic interactions is presented. We obtained a quantitatively accurate energy-entropy decomposition of the methane dimer cavity potential. The presented results confirm the previous results, and they provide new information regarding the thermodynamics of hydrophobic interactions. We show that the finite-difference approximation, which is widely used in molecular dynamic simulations for the energy-entropy decomposition of a free energy potential, can lead to a significant error.
Collapse
Affiliation(s)
- A N Morozov
- National Chiao Tung University, 1001 Ta Hsuen Road, Hsinchu, Taiwan Republic of China.
| | | |
Collapse
|
9
|
Calculation of the Free Energy and the Entropy of Macromolecular Systems by Computer Simulation. REVIEWS IN COMPUTATIONAL CHEMISTRY 2007. [DOI: 10.1002/9780470125892.ch1] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register]
|
10
|
Depaepe J, Ryckaert J, Bellemans A. Kinetics of the geometric isomerization of cyclohexene in a stochastic bath. Mol Phys 2006. [DOI: 10.1080/00268979300101021] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
11
|
Sun L, Siepmann JI, Schure MR. Conformation and Solvation Structure for an Isolated n-Octadecane Chain in Water, Methanol, and Their Mixtures. J Phys Chem B 2006; 110:10519-25. [PMID: 16722762 DOI: 10.1021/jp0602631] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Configurational-bias Monte Carlo simulations in the isobaric-isothermal ensemble (T = 323 K and p = 10 atm) were carried out to probe structural properties of an isolated n-octadecane chain solvated in water, methanol, water-rich, or methanol-rich mixtures and, for comparison, of an isolated chain in the gas phase and for neat liquid n-octadecane. The united-atom version of the TraPPE (transferable potentials for phase equilibria) force field was used to represent n-octadecane and methanol and the TIP-4P model was used for water. In all six environments, broad conformational distributions are observed and the n-octadecane chains are found to predominantly adopt extended, but not all-trans conformations. In addition, a small fraction of more collapsed conformations in which the chain ends approach each other is observed for aqueous hydration, the water-rich solvent mixture and the gas phase, but the simulation data do not support a simple two-state picture with folded and unfolded basins of attraction. For chains in these three "poor" solvent environments, the dihedral angles near the center of the chain show an enhancement of the gauche population. The ensemble of water-solvated chains with end-to-end contacts is preferentially found in a U-shaped conformation rather than a more globular state. An analysis of the local solvation structures in the water-methanol mixtures shows, as expected, an enrichment of the methyl group of methanol near the methylene and methyl segments of the n-octadecane chain. Interestingly, these local bead fractions are enhanced by factors of 2.5 and 1.5 for methyl and methylene segments reflecting the more hydrophobic nature of the former segments.
Collapse
Affiliation(s)
- Li Sun
- Department of Chemistry, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455, USA
| | | | | |
Collapse
|
12
|
Kato M, Tsuchiya H, Taniguchi Y. Raman Study of Pressure and Temperature Effects on the Conformational Equilibrium of Diethyl Disulfide in Solutions. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2005. [DOI: 10.1246/bcsj.78.1411] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/13/2022]
|
13
|
San Román-Zimbrón ML, Costas ME, Acevedo-Chávez R. Neutral hypoxanthine in aqueous solution: quantum chemical and Monte-Carlo studies. ACTA ACUST UNITED AC 2004. [DOI: 10.1016/j.theochem.2004.07.044] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
14
|
Imai T, Hirata F. Partial molar volume and compressibility of a molecule with internal degrees of freedom. J Chem Phys 2003. [DOI: 10.1063/1.1600437] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023] Open
|
15
|
|
16
|
|
17
|
Cui Q, Smith VH. Solvation Structure, Thermodynamics, and Molecular Conformational Equilibria for n-Butane in Water Analyzed by Reference Interaction Site Model Theory Using an All-Atom Solute Model. J Phys Chem B 2002. [DOI: 10.1021/jp020191n] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Affiliation(s)
- Qizhi Cui
- Department of Chemistry, Queen's University, Kingston, Ontario, Canada K7L 3N6
| | - Vedene H. Smith
- Department of Chemistry, Queen's University, Kingston, Ontario, Canada K7L 3N6
| |
Collapse
|
18
|
Abstract
This paper reviews the molecular theory of hydrophobic effects relevant to biomolecular structure and assembly in aqueous solution. Recent progress has resulted in simple, validated molecular statistical thermodynamic theories and clarification of confusing theories of decades ago. Current work is resolving effects of wider variations of thermodynamic state, e.g., pressure denaturation of soluble proteins, and more exotic questions such as effects of surface chemistry in treating stability of macromolecular structures in aqueous solution.
Collapse
Affiliation(s)
- Lawrence R Pratt
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA.
| |
Collapse
|
19
|
Wescott JT, Fisher LR, Hanna S. Use of thermodynamic integration to calculate the hydration free energies of n-alkanes. J Chem Phys 2002. [DOI: 10.1063/1.1431588] [Citation(s) in RCA: 42] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
20
|
Zuckerman DM, Woolf TB. Transition events in butane simulations: Similarities across models. J Chem Phys 2002. [DOI: 10.1063/1.1433501] [Citation(s) in RCA: 50] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
21
|
Cui Q, Smith VH. Analysis of solvation structure and thermodynamics of ethane and propane in water by reference interaction site model theory using all-atom models. J Chem Phys 2001. [DOI: 10.1063/1.1384421] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
|
22
|
Affiliation(s)
- Robert Günther
- Institut für Biochemie, Fakultät für Biowissenschaften, Pharmazie und Psychologie, Universität Leipzig, Talstr. 33, D-04103 Leipzig, Germany
| | - Hans-Jörg Hofmann
- Institut für Biochemie, Fakultät für Biowissenschaften, Pharmazie und Psychologie, Universität Leipzig, Talstr. 33, D-04103 Leipzig, Germany
| | - Krzysztof Kuczera
- Department of Chemistry and Department of Molecular Biosciences, University of Kansas, 2010 Mallot Hall, Lawrence, Kansas 66045, U.S.A
| |
Collapse
|
23
|
Ashbaugh HS, Garde S, Hummer G, Kaler EW, Paulaitis ME. Conformational equilibria of alkanes in aqueous solution: relationship to water structure near hydrophobic solutes. Biophys J 1999; 77:645-54. [PMID: 10423414 PMCID: PMC1300360 DOI: 10.1016/s0006-3495(99)76920-1] [Citation(s) in RCA: 41] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022] Open
Abstract
Conformational free energies of butane, pentane, and hexane in water are calculated from molecular simulations with explicit waters and from a simple molecular theory in which the local hydration structure is estimated based on a proximity approximation. This proximity approximation uses only the two nearest carbon atoms on the alkane to predict the local water density at a given point in space. Conformational free energies of hydration are subsequently calculated using a free energy perturbation method. Quantitative agreement is found between the free energies obtained from simulations and theory. Moreover, free energy calculations using this proximity approximation are approximately four orders of magnitude faster than those based on explicit water simulations. Our results demonstrate the accuracy and utility of the proximity approximation for predicting water structure as the basis for a quantitative description of n-alkane conformational equilibria in water. In addition, the proximity approximation provides a molecular foundation for extending predictions of water structure and hydration thermodynamic properties of simple hydrophobic solutes to larger clusters or assemblies of hydrophobic solutes.
Collapse
Affiliation(s)
- H S Ashbaugh
- Department of Chemical Engineering and Center for Molecular and Engineering Thermodynamics, University of Delaware, Newark, Delaware 19716, USA
| | | | | | | | | |
Collapse
|
24
|
Kato M, Abe I, Taniguchi Y. Raman study of the trans–gauche conformational equilibrium of 1,2-dichloroethane in water: Experimental evidence for the hydrophobic effect. J Chem Phys 1999. [DOI: 10.1063/1.479136] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
25
|
Abstract
A Gaussian solvent-exclusion model for the solvation free energy is developed. It is based on theoretical considerations and parametrized with experimental data. When combined with the CHARMM 19 polar hydrogen energy function, it provides an effective energy function (EEF1) for proteins in solution. The solvation model assumes that the solvation free energy of a protein molecule is a sum of group contributions, which are determined from values for small model compounds. For charged groups, the self-energy contribution is accounted for primarily by the exclusion model. Ionic side-chains are neutralized, and a distance-dependent dielectric constant is used to approximate the charge-charge interactions in solution. The resulting EEF1 is subjected to a number of tests. Molecular dynamics simulations at room temperature of several proteins in their native conformation are performed, and stable trajectories are obtained. The deviations from the experimental structures are similar to those observed in explicit water simulations. The calculated enthalpy of unfolding of a polyalanine helix is found to be in good agreement with experimental data. Results reported elsewhere show that EEF1 clearly distinguishes correctly from incorrectly folded proteins, both in static energy evaluations and in molecular dynamics simulations and that unfolding pathways obtained by high-temperature molecular dynamics simulations agree with those obtained by explicit water simulations. Thus, this energy function appears to provide a realistic first approximation to the effective energy hypersurface of proteins.
Collapse
Affiliation(s)
- T Lazaridis
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA
| | | |
Collapse
|
26
|
Madan B, Sharp K. Changes in water structure induced by a hydrophobic solute probed by simulation of the water hydrogen bond angle and radial distribution functions. Biophys Chem 1999; 78:33-41. [PMID: 10343383 DOI: 10.1016/s0301-4622(98)00227-0] [Citation(s) in RCA: 42] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
Abstract
In order to better characterize changes in water structure induced by a hydrophobic solute the oxygen-oxygen and hydrogen-hydrogen radial distribution functions (goo(r), ghh(r)) and the hydrogen bond angle distribution function p(theta) for water molecules in the first hydration shell of the tetramethyl ammonium (TMA) cation were computed using Monte Carlo simulations. goo(r) and ghh(r) were corrected for the effect of solute volume exclusion on the local solvent density so that intrinsic structural changes independent of local solvent density variations could be detected. Comparison of ghh(r) of TMA's first hydration shell water with ghh(r) for bulk water shows subtle but clear evidence of structure formation induced by the ion. These changes in ghh(r) are very similar to those seen experimentally for larger tetra-alkyl ammonium ions in previous neutron diffraction experiments. Larger changes in p(theta) in the first hydration shell of TMA were seen. Comparison of changes in p(theta) with changes in goo(r) and ghh(r) show that the angle distribution function provides the most sensitive way to analyze water structure changes associated with hydrophobic solvation.
Collapse
Affiliation(s)
- B Madan
- Johnson Research Foundation, Department of Biochemistry and Biophysics, University of Pennsylvania, Philadelphia 19104-6059, USA
| | | |
Collapse
|
27
|
Levy RM, Gallicchio E. Computer simulations with explicit solvent: recent progress in the thermodynamic decomposition of free energies and in modeling electrostatic effects. Annu Rev Phys Chem 1999; 49:531-67. [PMID: 9933909 DOI: 10.1146/annurev.physchem.49.1.531] [Citation(s) in RCA: 203] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
This review focuses on recent progress in two areas in which computer simulations with explicit solvent are being applied: the thermodynamic decomposition of free energies, and modeling electrostatic effects. The computationally intensive nature of these simulations has been an obstacle to the systematic study of many problems in solvation thermodynamics, such as the decomposition of solvation and ligand binding free energies into component enthalpies and entropies. With the revolution in computer power continuing, these problems are ripe for study but require the judicious choice of algorithms and approximations. We provide a critical evaluation of several numerical approaches to the thermodynamic decomposition of free energies and summarize applications in the current literature. Progress in computer simulations with explicit solvent of charge perturbations in biomolecules was slow in the early 1990s because of the widespread use of truncated Coulomb potentials in these simulations, among other factors. Development of the sophisticated technology described in this review to handle the long-range electrostatic interactions has increased the predictive power of these simulations to the point where comparisons between explicit and continuum solvent models can reveal differences that have their true physical origin in the inherent molecularity of the surrounding medium.
Collapse
Affiliation(s)
- R M Levy
- Department of Chemistry, Rutgers, State University of New Jersey, Piscataway 08855-0939, USA.
| | | |
Collapse
|
28
|
|
29
|
Ashbaugh HS, Kaler EW, Paulaitis ME. Hydration and conformational equilibria of simple hydrophobic and amphiphilic solutes. Biophys J 1998; 75:755-68. [PMID: 9675177 PMCID: PMC1299750 DOI: 10.1016/s0006-3495(98)77565-4] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023] Open
Abstract
We consider whether the continuum model of hydration optimized to reproduce vacuum-to-water transfer free energies simultaneously describes the hydration free energy contributions to conformational equilibria of the same solutes in water. To this end, transfer and conformational free energies of idealized hydrophobic and amphiphilic solutes in water are calculated from explicit water simulations and compared to continuum model predictions. As benchmark hydrophobic solutes, we examine the hydration of linear alkanes from methane through hexane. Amphiphilic solutes were created by adding a charge of +/-1e to a terminal methyl group of butane. We find that phenomenological continuum parameters fit to transfer free energies are significantly different from those fit to conformational free energies of our model solutes. This difference is attributed to continuum model parameters that depend on solute conformation in water, and leads to effective values for the free energy/surface area coefficient and Born radii that best describe conformational equilibrium. In light of these results, we believe that continuum models of hydration optimized to fit transfer free energies do not accurately capture the balance between hydrophobic and electrostatic contributions that determines the solute conformational state in aqueous solution.
Collapse
Affiliation(s)
- H S Ashbaugh
- Department of Chemical Engineering and Center for Molecular and Engineering Thermodynamics, University of Delaware, Newark, Delaware 19716, USA.
| | | | | |
Collapse
|
30
|
Lazaridis T. Inhomogeneous Fluid Approach to Solvation Thermodynamics. 2. Applications to Simple Fluids. J Phys Chem B 1998. [DOI: 10.1021/jp972358w] [Citation(s) in RCA: 136] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Themis Lazaridis
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138
| |
Collapse
|
31
|
Hermone A, Kuczera K. Free-energy simulations of the retinal cis --> trans isomerization in bacteriorhodopsin. Biochemistry 1998; 37:2843-53. [PMID: 9485435 DOI: 10.1021/bi9717789] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Free-energy profiles for ground-state cis --> trans isomerization of retinal in vacuum, in solution, and in the protein bacteriorhodopsin are calculated using free-energy simulations. The free-energy barriers in the protein were 9 kcal/mol for ionized Asp85 and 14 kcal/mol for neutral Asp85, significantly lower than those found in solution (18 kcal/mol) or vacuum (19 kcal/mol). Therefore, bacteriorhodopsin can be said to act as a catalyst in the isomerization. The barrier in the protein is due mainly to stabilization of the transition state through favorable nonbonded interactions with the protein part of the system, with internal strain and interactions with solvent playing minor roles. The protonated Asp85 simulation models the behavior of the system in the N --> O transition. Our calculated 14 kcal/mol barrier and 4-ms relaxation time for this process are in excellent agreement with experimentally measured values of 12 kcal/mol and 5 ms, respectively. The ionized Asp85 simulation models two hypothetical processes: the N --> O transition with a proton removed from Asp85 and the initial BR568 --> L transition on the ground-state energy surface. The cis-trans isomerization barrier in this system is 9 kcal/mol, the lowest of all the studied cases. The presence of the charged carboxylate group in the ionized Asp85 system leads to strong stabilization of the transition state by interactions with the surroundings and changes the distance between Asp85 and the Schiff base proton compared to the corresponding distance in the neutral Asp85 system. This suggests that the protonation of Asp85 plays an important role in regulating access to the Schiff base proton. For both Asp85 ionization states the calculated cis-trans free-energy difference was close to 0, indicating that the protein can accommodate both retinal isomers equally well. The computed negligible difference between the N and O free-energy levels is in accord with experimental data.
Collapse
Affiliation(s)
- A Hermone
- Department of Chemistry, University of Kansas, 2010 Malott Hall, Lawrence, Kansas 66045, USA
| | | |
Collapse
|
32
|
SanRomán-Zimbrón ML, Ortega-Blake I. On the molecular basis of hydrophobicity: A Monte Carlo study of propionic acid hydration. J Chem Phys 1997. [DOI: 10.1063/1.474676] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
33
|
Wang Y, Kuczera K. Multidimensional Conformational Free Energy Surface Exploration: Helical States of Alanand AibnPeptides. J Phys Chem B 1997. [DOI: 10.1021/jp964027+] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
34
|
Engkvist O, Karlström G. A Monte Carlo simulation study of the temperature dependence for the conformation distribution of 1,2-dimethoxyethane in water. J Chem Phys 1997. [DOI: 10.1063/1.473152] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
35
|
|
36
|
Engkvist O, Karlström G. A method to calculate the probability distribution for systems with large energy barriers. Chem Phys 1996. [DOI: 10.1016/s0301-0104(96)00247-9] [Citation(s) in RCA: 26] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
37
|
Kuczera K. Free energy simulations of axial contacts in sickle-cell hemoglobin. Biopolymers 1996; 39:221-42. [PMID: 8679951 DOI: 10.1002/(sici)1097-0282(199608)39:2%3c221::aid-bip10%3e3.0.co;2-h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
Molecular dynamics simulations have been used to investigate the thermodynamic stability of axial contacts in sickle-cell hemoglobin (HbS). Free energy changes were evaluated for the point mutation beta 121 Glu --> Gln in the axial contact region of HbS crystals. The calculations predict a free energy change of-3.6 kcal/mol per contact for the mutation, which is in qualitative agreement with experimental observations of aggravated sickling found in the double mutant Hb D Los Angeles (beta 6 Glu --> Val. beta 121 Glu --> Gln) relative to HbS (beta 6 Glu --> Val). The beta 121 Glu is sequestered in a salt link with beta 17 Lys located on the same polypeptide chain, making the Glu interactions with its surroundings similar in aggregates and individual hemoglobins. Due to this cancellation of the large electrostatic Glu contributions, the weak nonspecific interactions between the Gln and the neighboring polypeptide chain are the main contributing factor to the enhanced aggregation of Hb D Los Angeles relative to HbS. Together with the previous study of the lateral contact [K. Kuczera et al. (1990) Proceedings of the National Academy of Science USA, Vol. 87, pp, 8481-8485], the present results provide a more complete picture of the forces driving the sickling aggregation. A comparison of different treatments of internal flexibility in free energy simulations and analysis of rate of convergence of the different calculated properties has also been performed.
Collapse
Affiliation(s)
- K Kuczera
- Department of Chemistry, University of Kansas, Lawrence 66045, USA
| |
Collapse
|
38
|
Ashbaugh HS, Paulaitis ME. Entropy of Hydrophobic Hydration: Extension to Hydrophobic Chains. ACTA ACUST UNITED AC 1996. [DOI: 10.1021/jp952387b] [Citation(s) in RCA: 86] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Henry S. Ashbaugh
- Center for Molecular and Engineering Thermodynamics, Department of Chemical Engineering, University of Delaware, Newark, Delaware 19716
| | - Michael E. Paulaitis
- Center for Molecular and Engineering Thermodynamics, Department of Chemical Engineering, University of Delaware, Newark, Delaware 19716
| |
Collapse
|
39
|
Polson JM, Burnell EE. Conformational equilibrium and orientational ordering:1H‐nuclear magnetic resonance of butane in a nematic liquid crystal. J Chem Phys 1995. [DOI: 10.1063/1.470367] [Citation(s) in RCA: 34] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
40
|
Pellegrini M, Doniach S. Modeling solvation contributions to conformational free energy changes of biomolecules using a potential of mean force expansion. J Chem Phys 1995. [DOI: 10.1063/1.470503] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
41
|
Beglov D, Roux B. Finite representation of an infinite bulk system: Solvent boundary potential for computer simulations. J Chem Phys 1994. [DOI: 10.1063/1.466711] [Citation(s) in RCA: 793] [Impact Index Per Article: 26.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
42
|
|
43
|
Is the water-induced potential of mean torsion of n-butane transferable to longer n-alkanes? Chem Phys Lett 1993. [DOI: 10.1016/0009-2614(93)89316-a] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
44
|
Depaepe JM, Ryckaert JP, Paci E, Ciccotti G. Sampling of molecular conformations by molecular dynamics techniques. Mol Phys 1993. [DOI: 10.1080/00268979300101411] [Citation(s) in RCA: 20] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
45
|
Smith DE, Haymet ADJ. Free energy, entropy, and internal energy of hydrophobic interactions: Computer simulations. J Chem Phys 1993. [DOI: 10.1063/1.464809] [Citation(s) in RCA: 236] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
46
|
Rashin AA. Aspects of protein energetics and dynamics. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 1993; 60:73-200. [PMID: 8362069 DOI: 10.1016/0079-6107(93)90017-e] [Citation(s) in RCA: 61] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Affiliation(s)
- A A Rashin
- Biosym Technologies Inc, Parsippany, NJ 07054
| |
Collapse
|
47
|
Dunn WJ, Nagy PI. Relative logP and solution structure for small organic solutes in the chloroform/water system using monte carlo methods. J Comput Chem 1992. [DOI: 10.1002/jcc.540130409] [Citation(s) in RCA: 18] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
48
|
Tobias DJ, Sneddon SF, Brooks CL. Reverse turns in blocked dipeptides are intrinsically unstable in water. J Mol Biol 1990; 216:783-96. [PMID: 2258940 DOI: 10.1016/0022-2836(90)90399-7] [Citation(s) in RCA: 71] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
We have carried out molecular dynamics simulations to study the conformational equilibria of two blocked dipeptides, Ac-Ala-Ala-NHMe and trans-Ac-Pro-Ala-NHMe, in water (Ac, amino-terminal blocking group COCH3; NHMe, carboxy-terminal blocking group NHCH3). Using specialized sampling techniques we computed free-energy surfaces as functions of a conformation co-ordinate that corresponds to hydrogen-bonded reverse turns at small values and to extended conformations at large values. The free-energy difference between hydrogen-bonded reverse turn conformations and extended conformations, determined from the equilibrium constants for reverse turn unfolding, is approximately -5 kcal/mole for Ac-Ala-Ala-NHMe, and -10 kcal/mole for Ac-Pro-Ala-NHMe. These results demonstrate that reverse turns in blocked dipeptides are intrinsically unstable in water. That is, in the absence of strongly stabilizing sequence-specific inter-residue interactions involving side-chains and/or charged terminal groups, the extended conformations of small peptides are highly favored in solution. By thermodynamically decomposing the free-energy differences, we found that the peptide-water entropy is the primary reason for the exceptional stability of the extended conformations of both peptides, and that the differences between the two peptides are primarily due to differences in the peptide-water interactions. In addition, we assessed the "proline effect" on the conformational equilibria by comparing the differences in configurational entropies between the reverse turn and extended conformations of the two peptides. As expected, the extended conformation of the Pro-Ala peptide is destabilized by reduced configurational entropy, but the effect is negligible in the blocked dipeptides. Finally, we compared our results with the results of several other experimental studies to identify some of the specific interactions that may be responsible for stabilizing reverse turns in small peptides in solution.
Collapse
Affiliation(s)
- D J Tobias
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, PA 15213
| | | | | |
Collapse
|