1
|
Piotto S, Sessa L, Piotto A, Nardiello AM, Concilio S. Plausible Emergence of Autocatalytic Cycles under Prebiotic Conditions. Life (Basel) 2019; 9:life9020033. [PMID: 30987308 PMCID: PMC6617172 DOI: 10.3390/life9020033] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2019] [Revised: 03/22/2019] [Accepted: 04/02/2019] [Indexed: 11/29/2022] Open
Abstract
The emergence of life in a prebiotic world is an enormous scientific question of paramount philosophical importance. Even when life (in any sense we can define it) can be observed and replicated in the laboratory, it is only an indication of one possible pathway for life emergence, and is by no means be a demonstration of how life really emerged. The best we can hope for is to indicate plausible chemical–physical conditions and mechanisms that might lead to self-organizing and autopoietic systems. Here we present a stochastic simulation, based on chemical reactions already observed in prebiotic environments, that might help in the design of new experiments. We will show how the definition of simple rules for the synthesis of random peptides may lead to the appearance of networks of autocatalytic cycles and the emergence of memory.
Collapse
Affiliation(s)
- Stefano Piotto
- Department of Pharmacy, University of Salerno, Via Giovanni Paolo II, 132, 84084 Fisciano SA, Italy.
| | - Lucia Sessa
- Department of Pharmacy, University of Salerno, Via Giovanni Paolo II, 132, 84084 Fisciano SA, Italy.
| | - Andrea Piotto
- Department of Pharmacy, University of Salerno, Via Giovanni Paolo II, 132, 84084 Fisciano SA, Italy.
| | - Anna Maria Nardiello
- Department of Pharmacy, University of Salerno, Via Giovanni Paolo II, 132, 84084 Fisciano SA, Italy.
| | - Simona Concilio
- Department of Industrial Engineering, University of Salerno, Via Giovanni Paolo II, 132, 84084 Fisciano SA, Italy.
| |
Collapse
|
2
|
König G, Pickard FC, Huang J, Thiel W, MacKerell AD, Brooks BR, York DM. A Comparison of QM/MM Simulations with and without the Drude Oscillator Model Based on Hydration Free Energies of Simple Solutes. Molecules 2018; 23:E2695. [PMID: 30347691 PMCID: PMC6222909 DOI: 10.3390/molecules23102695] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2018] [Revised: 10/15/2018] [Accepted: 10/16/2018] [Indexed: 12/01/2022] Open
Abstract
Maintaining a proper balance between specific intermolecular interactions and non-specific solvent interactions is of critical importance in molecular simulations, especially when predicting binding affinities or reaction rates in the condensed phase. The most rigorous metric for characterizing solvent affinity are solvation free energies, which correspond to a transfer from the gas phase into solution. Due to the drastic change of the electrostatic environment during this process, it is also a stringent test of polarization response in the model. Here, we employ both the CHARMM fixed charge and polarizable force fields to predict hydration free energies of twelve simple solutes. The resulting classical ensembles are then reweighted to obtain QM/MM hydration free energies using a variety of QM methods, including MP2, Hartree⁻Fock, density functional methods (BLYP, B3LYP, M06-2X) and semi-empirical methods (OM2 and AM1 ). Our simulations test the compatibility of quantum-mechanical methods with molecular-mechanical water models and solute Lennard⁻Jones parameters. In all cases, the resulting QM/MM hydration free energies were inferior to purely classical results, with the QM/MM Drude force field predictions being only marginally better than the QM/MM fixed charge results. In addition, the QM/MM results for different quantum methods are highly divergent, with almost inverted trends for polarizable and fixed charge water models. While this does not necessarily imply deficiencies in the QM models themselves, it underscores the need to develop consistent and balanced QM/MM interactions. Both the QM and the MM component of a QM/MM simulation have to match, in order to avoid artifacts due to biased solute⁻solvent interactions. Finally, we discuss strategies to improve the convergence and efficiency of multi-scale free energy simulations by automatically adapting the molecular-mechanics force field to the target quantum method.
Collapse
Affiliation(s)
- Gerhard König
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20892, USA.
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA.
- Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany.
| | - Frank C Pickard
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20892, USA.
| | - Jing Huang
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20892, USA.
- Department of Pharmaceutical Science, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, MD 21201, USA.
- School of Life Sciences, Westlake University, 18 Shilongshan Street, Hangzhou 310024, China.
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany.
| | - Alexander D MacKerell
- Department of Pharmaceutical Science, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, MD 21201, USA.
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20892, USA.
| | - Darrin M York
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA.
| |
Collapse
|
3
|
Zhang H, Jiang Y, Cui Z, Yin C. Force Field Benchmark of Amino Acids. 2. Partition Coefficients between Water and Organic Solvents. J Chem Inf Model 2018; 58:1669-1681. [DOI: 10.1021/acs.jcim.8b00493] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Yang Jiang
- Beijing Key Lab of Bioprocess, College of Life Science and Technology, Beijing University of Chemical Technology, Box 53, 100029 Beijing, China
| | - Ziheng Cui
- Beijing Key Lab of Bioprocess, College of Life Science and Technology, Beijing University of Chemical Technology, Box 53, 100029 Beijing, China
| | - Chunhua Yin
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| |
Collapse
|
4
|
Capping effects on polymorphic Aβ 16-21 amyloids depend on their size: A molecular dynamics simulation study. Biophys Chem 2017; 232:1-11. [PMID: 29046256 DOI: 10.1016/j.bpc.2017.09.003] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2017] [Revised: 09/19/2017] [Accepted: 09/21/2017] [Indexed: 11/24/2022]
Abstract
Understanding Aβ amyloid oligomers associated with neuro-degenerative diseases is needed due to their toxic characteristics and mediation of amyloid fibril growth. Depending on various physiological circumstances such as ionic strength, metal ion, and point-residue mutation, oligomeric amyloids exhibit polymorphic behavior and structural stabilities, i.e. showing different conformation and stabilities. Specifically, experimental and computational researchers have found that the capping modulates the physical and chemical properties of amyloids by preserving electrostatic energy interactions, which is one of the dominant factors for amyloid stability. Still, there is no detailed knowledge for the polymorphic amyloids with reflecting the terminal capping effects. In the present study, we investigated the role of terminal capping (i.e. N-terminal acetylation and C-terminal amidation) on polymorphic Aβ16-21 amyloid oligomer and protofibrils via molecular dynamics (MD) simulations. We found that the capping effects have differently altered the conformation of polymorphic antiparallel-homo and -hetero Aβ16-21 amyloid oligomer, but not Aβ16-21 amyloid protofibrils. However, regardless of polymorphic composition of the amyloids, the capping induces the thermodynamic instabilities of Aβ16-21 amyloid oligomers, but does not show any distinct affect on Aβ16-21 amyloid protofibrils. Specifically, among the molecular mechanic factors, electrostatic energy dominantly contributes the thermodynamic stability of the Aβ16-21 amyloids. We hope that our computation study about the role of the capping effects on the polymorphic amyloids will facilitate additional efforts to enhance degradation of amyloids and to design a selective drug in the future.
Collapse
|
5
|
Drake JA, Harris RC, Pettitt BM. Solvation Thermodynamics of Oligoglycine with Respect to Chain Length and Flexibility. Biophys J 2017; 111:756-767. [PMID: 27558719 DOI: 10.1016/j.bpj.2016.07.013] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2016] [Revised: 06/30/2016] [Accepted: 07/06/2016] [Indexed: 01/24/2023] Open
Abstract
Oligoglycine is a backbone mimic for all proteins and is prevalent in the sequences of intrinsically disordered proteins. We have computed the absolute chemical potential of glycine oligomers at infinite dilution by simulation with the CHARMM36 and Amber ff12SB force fields. We performed a thermodynamic decomposition of the solvation free energy (ΔG(sol)) of Gly2-5 into enthalpic (ΔH(sol)) and entropic (ΔS(sol)) components as well as their van der Waals and electrostatic contributions. Gly2-5 was either constrained to a rigid/extended conformation or allowed to be completely flexible during simulations to assess the effects of flexibility on these thermodynamic quantities. For both rigid and flexible oligoglycine models, the decrease in ΔG(sol) with chain length is enthalpically driven with only weak entropic compensation. However, the apparent rates of decrease of ΔG(sol), ΔH(sol), ΔS(sol), and their elec and vdw components differ for the rigid and flexible models. Thus, we find solvation entropy does not drive aggregation for this system and may not explain the collapse of long oligoglycines. Additionally, both force fields yield very similar thermodynamic scaling relationships with respect to chain length despite both force fields generating different conformational ensembles of various oligoglycine chains.
Collapse
Affiliation(s)
- Justin A Drake
- Department of Biochemistry and Molecular Biology, Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, Galveston, Texas
| | - Robert C Harris
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania
| | - B Montgomery Pettitt
- Department of Biochemistry and Molecular Biology, Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, Galveston, Texas.
| |
Collapse
|
6
|
Tooski HF, Jabbari M, Farajtabar A. Solubility and Preferential Solvation of the Flavonoid Naringenin in Some Aqueous/Organic Solvent Mixtures. J SOLUTION CHEM 2016. [DOI: 10.1007/s10953-016-0526-2] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
7
|
Feizi S, Jabbari M, Farajtabar A. A systematic study on solubility and solvation of bioactive compound chrysin in some water + cosolvent mixtures. J Mol Liq 2016. [DOI: 10.1016/j.molliq.2016.05.019] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
|
8
|
Fu IW, Nguyen HD. Sequence-Dependent Structural Stability of Self-Assembled Cylindrical Nanofibers by Peptide Amphiphiles. Biomacromolecules 2015; 16:2209-19. [DOI: 10.1021/acs.biomac.5b00595] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Iris W. Fu
- Department
of Chemical Engineering
and Materials Science, University of California, Irvine, Irvine, California 92697, United States
| | - Hung D. Nguyen
- Department
of Chemical Engineering
and Materials Science, University of California, Irvine, Irvine, California 92697, United States
| |
Collapse
|
9
|
Porasso RD, Ale NM, Ciocco Aloia F, Masone D, Del Pópolo MG, Ben Altabef A, Gomez-Zavaglia A, Diaz SB, Vila JA. Interaction of glycine, lysine, proline and histidine with dipalmitoylphosphatidylcholine lipid bilayers: a theoretical and experimental study. RSC Adv 2015. [DOI: 10.1039/c5ra03236a] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
The interaction of unblocked glycine, lysine, proline and histidine with a DPPC bilayer was assessed using extensive atomistic molecular dynamics simulations.
Collapse
Affiliation(s)
- Rodolfo D. Porasso
- Instituto de Matemática Aplicada San Luis (IMASL)
- CONICET
- Universidad Nacional de San Luis
- Argentina
| | - Norma M. Ale
- Instituto de Química Física
- Facultad de Bioquímica, Química y Farmacia
- U. N. T
- San Lorenzo 456
- Argentina
| | - Facundo Ciocco Aloia
- CONICET
- Facultad de Ciencias Exactas y Naturales
- Universidad Nacional de Cuyo
- Mendoza
- Argentina
| | - Diego Masone
- CONICET
- Facultad de Ciencias Exactas y Naturales
- Universidad Nacional de Cuyo
- Mendoza
- Argentina
| | - Mario G. Del Pópolo
- CONICET
- Facultad de Ciencias Exactas y Naturales
- Universidad Nacional de Cuyo
- Mendoza
- Argentina
| | - Aida Ben Altabef
- Instituto de Química Física
- Facultad de Bioquímica, Química y Farmacia
- U. N. T
- San Lorenzo 456
- Argentina
| | - Andrea Gomez-Zavaglia
- Center for Research and Development in Food Cryotechnology (CIDCA, CCT-CONICET, La Plata)
- Argentina
| | - Sonia B. Diaz
- Instituto de Química Física
- Facultad de Bioquímica, Química y Farmacia
- U. N. T
- San Lorenzo 456
- Argentina
| | - Jorge A. Vila
- Instituto de Matemática Aplicada San Luis (IMASL)
- CONICET
- Universidad Nacional de San Luis
- Argentina
- Baker Laboratory of Chemistry and Chemical Biology
| |
Collapse
|
10
|
Fu IW, Markegard CB, Nguyen HD. Solvent effects on kinetic mechanisms of self-assembly by peptide amphiphiles via molecular dynamics simulations. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2015; 31:315-24. [PMID: 25488898 DOI: 10.1021/la503399x] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/16/2023]
Abstract
Peptide amphiphiles are known to form a variety of distinctive self-assembled nanostructures (including cylindrical nanofibers in hydrogels) dependent upon the solvent conditions. Using a novel coarse-grained model, large-scale molecular dynamics simulations are performed on a system of 800 peptide amphiphiles (sequence, palmitoyl-Val3Ala3Glu3) to elucidate kinetic mechanisms of molecular assembly as a function of the solvent conditions. The assembly process is found to occur via a multistep process with transient intermediates that ultimately leads to the stabilized nanostructures including open networks of β-sheets, cylindrical nanofibers, and elongated micelles. Different kinetic mechanisms are compared in terms of peptide secondary structures, solvent-accessible surface area, radius of gyration, relative shape anisotropy, intra/intermolecular interactions, and aggregate size dynamics to provide insightful information for the design of functional biomaterials.
Collapse
Affiliation(s)
- Iris W Fu
- Department of Chemical Engineering and Materials Science, University of California, Irvine , Irvine, California 92697, United States
| | | | | |
Collapse
|
11
|
Fu IW, Markegard CB, Chu BK, Nguyen HD. Role of hydrophobicity on self-assembly by peptide amphiphiles via molecular dynamics simulations. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2014; 30:7745-7754. [PMID: 24915982 DOI: 10.1021/la5012988] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Abstract
Using a novel coarse-grained model, large-scale molecular dynamics simulations were performed to examine self-assembly of 800 peptide amphiphiles (sequence palmitoyl-V3A3E3). Under suitable physiological conditions, these molecules readily assemble into nanofibers leading to hydrogel construction as observed in experiments. Our simulations capture this spontaneous self-assembly process, including formation of secondary structure, to identify morphological transitions of distinctive nanostructures. As the hydrophobic interaction is increased, progression from open networks of secondary structures toward closed cylindrical nanostructures containing either β-sheets or random coils are observed. Moreover, temperature effects are also determined to play an important role in regulating formation of secondary structures within those nanostructures. These understandings of the molecular interactions involved and the role of environmental factors on hydrogel formation provide useful insight for development of innovative smart biomaterials for biomedical applications.
Collapse
Affiliation(s)
- Iris W Fu
- Department of Chemical Engineering and Materials Science, University of California-Irvine , Irvine, California 92697, United States
| | | | | | | |
Collapse
|
12
|
Electrochemical analysis of diffusion behavior and nucleation mechanism for Dy(II) and Dy(III) in phosphonium-based ionic liquids. Electrochim Acta 2013. [DOI: 10.1016/j.electacta.2013.09.082] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
13
|
On the transferability of fractional contributions to the hydration free energy of amino acids. Theor Chem Acc 2013. [DOI: 10.1007/s00214-013-1343-y] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
14
|
Han W, Schulten K. Further optimization of a hybrid united-atom and coarse-grained force field for folding simulations: Improved backbone hydration and interactions between charged side chains. J Chem Theory Comput 2012; 8:4413-4424. [PMID: 23204949 PMCID: PMC3507460 DOI: 10.1021/ct300696c] [Citation(s) in RCA: 101] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
PACE, a hybrid force field which couples united-atom protein models with coarse-grained (CG) solvent, has been further optimized, aiming to improve itse ciency for folding simulations. Backbone hydration parameters have been re-optimized based on hydration free energies of polyalanyl peptides through atomistic simulations. Also, atomistic partial charges from all-atom force fields were combined with PACE in order to provide a more realistic description of interactions between charged groups. Using replica exchange molecular dynamics (REMD), ab initio folding using the new PACE has been achieved for seven small proteins (16 - 23 residues) with different structural motifs. Experimental data about folded states, such as their stability at room temperature, melting point and NMR NOE constraints, were also well reproduced. Moreover, a systematic comparison of folding kinetics at room temperature has been made with experiments, through standard MD simulations, showing that the new PACE may speed up the actual folding kinetics 5-10 times. Together with the computational speedup benefited from coarse-graining, the force field provides opportunities to study folding mechanisms. In particular, we used the new PACE to fold a 73-residue protein, 3D, in multiple 10 - 30 μs simulations, to its native states (C(α) RMSD ~ 0.34 nm). Our results suggest the potential applicability of the new PACE for the study of folding and dynamics of proteins.
Collapse
Affiliation(s)
- Wei Han
- Beckman Institute, University of Illinois at Urbana-Champaign, USA
- Center for Biophysics and Computational Biology, University of Illinois at Urbana-Champaign, USA
| | - Klaus Schulten
- Beckman Institute, University of Illinois at Urbana-Champaign, USA
- Center for Biophysics and Computational Biology, University of Illinois at Urbana-Champaign, USA
| |
Collapse
|
15
|
Ren P, Chun J, Thomas DG, Schnieders MJ, Marucho M, Zhang J, Baker NA. Biomolecular electrostatics and solvation: a computational perspective. Q Rev Biophys 2012; 45:427-91. [PMID: 23217364 PMCID: PMC3533255 DOI: 10.1017/s003358351200011x] [Citation(s) in RCA: 135] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
An understanding of molecular interactions is essential for insight into biological systems at the molecular scale. Among the various components of molecular interactions, electrostatics are of special importance because of their long-range nature and their influence on polar or charged molecules, including water, aqueous ions, proteins, nucleic acids, carbohydrates, and membrane lipids. In particular, robust models of electrostatic interactions are essential for understanding the solvation properties of biomolecules and the effects of solvation upon biomolecular folding, binding, enzyme catalysis, and dynamics. Electrostatics, therefore, are of central importance to understanding biomolecular structure and modeling interactions within and among biological molecules. This review discusses the solvation of biomolecules with a computational biophysics view toward describing the phenomenon. While our main focus lies on the computational aspect of the models, we provide an overview of the basic elements of biomolecular solvation (e.g. solvent structure, polarization, ion binding, and non-polar behavior) in order to provide a background to understand the different types of solvation models.
Collapse
Affiliation(s)
- Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin
| | | | | | | | - Marcelo Marucho
- Department of Physics and Astronomy, The University of Texas at San Antonio
| | - Jiajing Zhang
- Department of Biomedical Engineering, The University of Texas at Austin
| | - Nathan A. Baker
- To whom correspondence should be addressed. Pacific Northwest National Laboratory, PO Box 999, MSID K7-29, Richland, WA 99352. Phone: +1-509-375-3997,
| |
Collapse
|
16
|
König G, Boresch S. Non-Boltzmann sampling and Bennett's acceptance ratio method: How to profit from bending the rules. J Comput Chem 2010; 32:1082-90. [DOI: 10.1002/jcc.21687] [Citation(s) in RCA: 60] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2010] [Revised: 08/25/2010] [Accepted: 09/05/2010] [Indexed: 11/07/2022]
|
17
|
Hu CY, Kokubo H, Lynch GC, Bolen DW, Pettitt BM. Backbone additivity in the transfer model of protein solvation. Protein Sci 2010; 19:1011-22. [PMID: 20306490 DOI: 10.1002/pro.378] [Citation(s) in RCA: 65] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
The transfer model implying additivity of the peptide backbone free energy of transfer is computationally tested. Molecular dynamics simulations are used to determine the extent of change in transfer free energy (DeltaG(tr)) with increase in chain length of oligoglycine with capped end groups. Solvation free energies of oligoglycine models of varying lengths in pure water and in the osmolyte solutions, 2M urea and 2M trimethylamine N-oxide (TMAO), were calculated from simulations of all atom models, and DeltaG(tr) values for peptide backbone transfer from water to the osmolyte solutions were determined. The results show that the transfer free energies change linearly with increasing chain length, demonstrating the principle of additivity, and provide values in reasonable agreement with experiment. The peptide backbone transfer free energy contributions arise from van der Waals interactions in the case of transfer to urea, but from electrostatics on transfer to TMAO solution. The simulations used here allow for the calculation of the solvation and transfer free energy of longer oligoglycine models to be evaluated than is currently possible through experiment. The peptide backbone unit computed transfer free energy of -54 cal/mol/M compares quite favorably with -43 cal/mol/M determined experimentally.
Collapse
Affiliation(s)
- Char Y Hu
- Graduate Program in Structural and Computational Biology and Molecular Biophysics, Baylor College of Medicine, Houston, Texas 77030, USA
| | | | | | | | | |
Collapse
|
18
|
König G, Boresch S. Hydration free energies of amino acids: why side chain analog data are not enough. J Phys Chem B 2009; 113:8967-74. [PMID: 19507836 DOI: 10.1021/jp902638y] [Citation(s) in RCA: 53] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Using molecular dynamics based free energy simulations, we computed relative solvation free energies for pairs of N-acetyl-methylamide amino acids (Ala-Ser, Val-Thr, Phe-Tyr, Val-Ala, Thr-Ser, Phe-Ala, and Tyr-Ser) and compared the results with the relative solvation free energies of the corresponding pairs of side chain analogs. We observed differences in (relative) solvent affinity DeltaDeltaDeltaA between amino acids and side chain analogs of up to 66% or, in absolute numbers, 4.9 kcal/mol (Ala-Ser). To rationalize these findings, we estimated separately contributions from what we refer to as solvent exclusion and self-solvation. While the former accounts for the reduction in solute-solvent interactions as one part of the solute occludes other parts of the solute, the latter turned out to be the determining contribution for small polar amino acids and could be shown to arise from interactions between the polar backbone and the polar functional group of the respective side chain in the gas phase. Consequently, the solvent affinity of small polar amino acids depends strongly on the backbone conformation. Our results indicate that the still widely used group additivity-solvent exclusion assumption to estimate solvation free energies for large(r) molecules (such as peptides and proteins) from model compound data (such as side chain analogs) is insufficient. To illustrate practical consequences, we compare the explicit solvent results with those of implicit solvent models. While approaches based on the generalized Born model give results in (mostly) good agreement with explicit solvent, approaches relying (primarily) on the group additivity-solvent exclusion assumption fail to reproduce DeltaDeltaDeltaA. Finally, we briefly discuss the implications of our results for hydrophobicity scales.
Collapse
Affiliation(s)
- Gerhard König
- Department of Computional Biological Chemistry, Faculty of Chemistry, University of Vienna, Vienna, Austria
| | | |
Collapse
|
19
|
Leontidis E, Aroti A. Liquid Expanded Monolayers of Lipids As Model Systems to Understand the Anionic Hofmeister Series: 2. Ion Partitioning Is Mostly a Matter of Size. J Phys Chem B 2009; 113:1460-7. [DOI: 10.1021/jp809444n] [Citation(s) in RCA: 47] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- E. Leontidis
- Department of Chemistry, University of Cyprus, Nicosia 1678, Cyprus
| | - A. Aroti
- Department of Chemistry, University of Cyprus, Nicosia 1678, Cyprus
| |
Collapse
|
20
|
Partitioning of amino-acid analogues in a five-slab membrane model. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2008; 1778:2234-43. [DOI: 10.1016/j.bbamem.2008.06.014] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/19/2007] [Revised: 05/15/2008] [Accepted: 06/17/2008] [Indexed: 11/20/2022]
|
21
|
Benbrahim N, Rahmouni A, Ruiz-López MF. A theoretical study of medium effects on the structure of the glycine analogue aminomethylphosphonic acid. Phys Chem Chem Phys 2008; 10:5624-32. [DOI: 10.1039/b804694h] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
22
|
|
23
|
Kuwajima S, Kikuchi H, Fukuda M. Molecular-dynamics evaluation of fluid-phase equilibrium properties by a novel free-energy perturbation approach: Application to gas solubility and vapor pressure of liquid hexane. J Chem Phys 2006; 124:124111. [PMID: 16599666 DOI: 10.1063/1.2178321] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A novel free-energy perturbation method is developed for the computation of the free energy of transferring a molecule between fluid phases. The methodology consists in drawing a free-energy profile of the target molecule moving across a binary-phase structure built in the computer. The novelty of the method lies in the difference of the definition of the free-energy profile from the common definition. As an important element of the method, the process of making a correction to the transfer free energy with respect to the cutoff of intermolecular forces is elucidated. In order to examine the performance of the method in the application to fluid-phase equilibrium properties, molecular-dynamics computations are carried out for the evaluation of gas solubility and vapor pressure of liquid n-hexane at 298.15 K. The gas species treated are methane, ethane, propane, and n-butane, with the gas solubility expressed as Henry's constant. It is shown that the method works fine and calculated results are generally in good agreement with experiments. It is found that the cutoff correction is strikingly large, constituting a dominant part of the calculated transfer free energy at the cutoff of 8 A.
Collapse
Affiliation(s)
- Satoru Kuwajima
- NanoSimulation Associates, 825-1 Amado-cho, Villa DE 201, Hanamigawa-ku, Chiba-shi, Chiba 262-0043, Japan.
| | | | | |
Collapse
|