1
|
Chen J, Qiu Z, Huang J. Structure and Dynamics of Confined Water Inside Diphenylalanine Peptide Nanotubes. ACS OMEGA 2023; 8:42936-42950. [PMID: 38024738 PMCID: PMC10652825 DOI: 10.1021/acsomega.3c06071] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/16/2023] [Revised: 09/22/2023] [Accepted: 10/19/2023] [Indexed: 12/01/2023]
Abstract
Diphenylalanine (FF) peptides exhibit a unique ability to self-assemble into nanotubes with confined water molecules playing pivotal roles in their structure and function. This study investigates the structure and dynamics of diphenylalanine peptide nanotubes (FFPNTs) using all-atom molecular dynamics (MD) and grand canonical Monte Carlo combined with MD (GCMC/MD) simulations with both the CHARMM additive and Drude polarizable force fields. The occupancy and dynamics of confined water molecules were also examined. It was found that less than 2 confined water molecules per FF help stabilize the FFPNTs on the x-y plane. Analyses of the kinetics of confined water molecules revealed distinctive transport behaviors for bound and free water, and their respective diffusion coefficients were compared. Our results validate the importance of polarizable force field models in studying peptide nanotubes and provide insights into our understanding of nanoconfined water.
Collapse
Affiliation(s)
- Jinfeng Chen
- College
of Life Sciences, Zhejiang University, Hangzhou, Zhejiang 310027, China
- Key
Laboratory of Structural Biology of Zhejiang Province, School of Life
Sciences, Westlake University, Hangzhou, Zhejiang 310024, China
- Westlake
AI Therapeutics Lab, Westlake Laboratory
of Life Sciences and Biomedicine, Hangzhou, Zhejiang 310024, China
| | - Zongyang Qiu
- Key
Laboratory of Structural Biology of Zhejiang Province, School of Life
Sciences, Westlake University, Hangzhou, Zhejiang 310024, China
- Westlake
AI Therapeutics Lab, Westlake Laboratory
of Life Sciences and Biomedicine, Hangzhou, Zhejiang 310024, China
| | - Jing Huang
- Key
Laboratory of Structural Biology of Zhejiang Province, School of Life
Sciences, Westlake University, Hangzhou, Zhejiang 310024, China
- Westlake
AI Therapeutics Lab, Westlake Laboratory
of Life Sciences and Biomedicine, Hangzhou, Zhejiang 310024, China
| |
Collapse
|
2
|
Wang A, Peng X, Li Y, Zhang D, Zhang Z, Li G. Quality of force fields and sampling methods in simulating pepX peptides: a case study for intrinsically disordered proteins. Phys Chem Chem Phys 2021; 23:2430-2437. [PMID: 33459730 DOI: 10.1039/d0cp05484d] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
Intrinsically disordered proteins (IDPs) are a group of proteins that lack well-defined structures under native conditions and carry out crucial physiological functions in various biochemical pathways. Due to the heterogeneous nature of IDPs, molecular dynamics simulations have been extensively adopted to investigate the conformational ensembles and dynamic properties of these proteins. However, their accuracy remains limited by the development of force fields and sampling algorithms. Here, we evaluated the quality of both force fields and enhanced sampling algorithms based on five short pepX peptides. Our results show that the more extended conformational ensembles sampled by the AMOEBA polarizable force field present a higher ability to reproduce experimental NMR observables than AMBER and CHARMM classical force fields. Moreover, a better agreement with experiments is achieved in the simulation of IaMD (integrated accelerated molecular dynamics) than in aMD (accelerated molecular dynamics). The results together indicate that the combination of AMOEBA force field and IaMD enhanced sampling might be a better choice for simulating IDPs. This work may provide important clues for developments and applications of force fields and enhanced sampling methods in future simulations of IDPs.
Collapse
Affiliation(s)
- Anhui Wang
- State Key Laboratory of Fine Chemicals, School of Chemistry, Dalian University of Technology, Dalian 116024, China.
| | | | | | | | | | | |
Collapse
|
3
|
Ploetz EA, Karunaweera S, Bentenitis N, Chen F, Dai S, Gee MB, Jiao Y, Kang M, Kariyawasam NL, Naleem N, Weerasinghe S, Smith PE. Kirkwood-Buff-Derived Force Field for Peptides and Proteins: Philosophy and Development of KBFF20. J Chem Theory Comput 2021; 17:2964-2990. [PMID: 33878263 DOI: 10.1021/acs.jctc.1c00075] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
A new classical nonpolarizable force field, KBFF20, for the simulation of peptides and proteins is presented. The force field relies heavily on the use of Kirkwood-Buff theory to provide a comparison of simulated and experimental Kirkwood-Buff integrals for solutes containing the functional groups common in proteins, thus ensuring intermolecular interactions that provide a good balance between the peptide-peptide, peptide-solvent, and solvent-solvent distributions observed in solution mixtures. In this way, it differs significantly from other biomolecular force fields. Further development and testing of the intermolecular potentials are presented here. Subsequently, rotational potentials for the ϕ/ψ and χ dihedral degrees of freedom are obtained by analysis of the Protein Data Bank, followed by small modifications to provide a reasonable balance between simulated and observed α and β percentages for small peptides. This, the first of two articles, describes in detail the philosophy and development behind KBFF20.
Collapse
Affiliation(s)
- Elizabeth A Ploetz
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Sadish Karunaweera
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Nikolaos Bentenitis
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Feng Chen
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Shu Dai
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Moon B Gee
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Yuanfang Jiao
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Myungshim Kang
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Nilusha L Kariyawasam
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Nawavi Naleem
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | | | - Paul E Smith
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| |
Collapse
|
4
|
Kulkarni A, Bortz M, Küfer KH, Kohns M, Hasse H. Multicriteria Optimization of Molecular Models of Water Using a Reduced Units Approach. J Chem Theory Comput 2020; 16:5127-5138. [PMID: 32609517 DOI: 10.1021/acs.jctc.0c00301] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Multicriteria optimization (MCO) is used to parametrize molecular models of water. The set of the best possible compromises between different objectives, the Pareto set, is determined. Calculating Pareto sets for optimization problems involving molecular simulations is computationally expensive. Therefore, we use a novel, highly efficient method, which is based on the fact that numerical results from molecular simulations can be interpreted as dimensionless numbers. Hence, they carry information on an entire class of models in physical units. This approach was applied here for the MCO of water models of the "one-center Lennard-Jones + point charge" type, in which the objectives were the quality of the description of the vapor pressure, liquid density, and enthalpy of vaporization. The results were compared to models from the literature. Significant improvements were observed. The new optimization method for the development of molecular models is efficient, robust, and broadly applicable.
Collapse
Affiliation(s)
- Aditya Kulkarni
- Laboratory of Engineering Thermodynamics (LTD), Technische Universität Kaiserslautern (TUK), 67663 Kaiserslautern, Germany
| | - Michael Bortz
- Fraunhofer Institute for Industrial Mathematics (ITWM), 67663 Kaiserslautern, Germany
| | - Karl-Heinz Küfer
- Fraunhofer Institute for Industrial Mathematics (ITWM), 67663 Kaiserslautern, Germany
| | - Maximilian Kohns
- Laboratory of Engineering Thermodynamics (LTD), Technische Universität Kaiserslautern (TUK), 67663 Kaiserslautern, Germany
| | - Hans Hasse
- Laboratory of Engineering Thermodynamics (LTD), Technische Universität Kaiserslautern (TUK), 67663 Kaiserslautern, Germany
| |
Collapse
|
5
|
Wang A, Zhang Z, Li G. Higher Accuracy Achieved in the Simulations of Protein Structure Refinement, Protein Folding, and Intrinsically Disordered Proteins Using Polarizable Force Fields. J Phys Chem Lett 2018; 9:7110-7116. [PMID: 30514082 DOI: 10.1021/acs.jpclett.8b03471] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
The accuracy of molecular mechanics force fields is of vital importance in biomolecular simulations. However, the admittedly more accurate polarizable force fields were recently reported to be less able to reproduce the experimental properties in comparison to additive force fields in some cases. Here, we perform long-time-scale molecular dynamics simulations to systematically evaluate the effect of explicit electronic polarization in polarizable force fields. The results show that the inclusion of electrostatic polarization effect in polarizable force fields can improve their accuracies in protein structure refinement and generate conformational ensembles more approximate to experiments for intrinsically disordered proteins. In contrast, it is difficult for polarizable force fields to approach the native structure, let alone to predict the native state when it is unknown a priori in the real protein structure predictions. We speculate that these effects might be attributed to the preference of protein-water interactions in polarizable force fields.
Collapse
Affiliation(s)
- Anhui Wang
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics , Dalian Institute of Chemical Physics, Chinese Academy of Sciences , Dalian 116023 , China
- State Key Laboratory of Fine Chemicals, School of Chemistry , Dalian University of Technology , Dalian 116024 , China
| | - Zhichao Zhang
- State Key Laboratory of Fine Chemicals, School of Chemistry , Dalian University of Technology , Dalian 116024 , China
| | - Guohui Li
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics , Dalian Institute of Chemical Physics, Chinese Academy of Sciences , Dalian 116023 , China
| |
Collapse
|
6
|
Long GE, Dhakal P, Redeker BN, Paluch AS. Using limiting activity coefficients to efficiently evaluate the ability of fixed-charge force fields to model miscible water plus cosolvent mixtures. MOLECULAR SIMULATION 2018. [DOI: 10.1080/08927022.2018.1531399] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Garrett E. Long
- Department of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, OH, USA
| | - Pratik Dhakal
- Department of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, OH, USA
| | - Bryce N. Redeker
- Department of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, OH, USA
| | - Andrew S. Paluch
- Department of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, OH, USA
| |
Collapse
|
7
|
Milzetti J, Nayar D, van der Vegt NFA. Convergence of Kirkwood–Buff Integrals of Ideal and Nonideal Aqueous Solutions Using Molecular Dynamics Simulations. J Phys Chem B 2018; 122:5515-5526. [DOI: 10.1021/acs.jpcb.7b11831] [Citation(s) in RCA: 53] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jasmin Milzetti
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Center of Smart Interfaces, Technische Universität Darmstadt, Alarich-Weiss-Straße 10, 64287, Darmstadt, Germany
| | - Divya Nayar
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Center of Smart Interfaces, Technische Universität Darmstadt, Alarich-Weiss-Straße 10, 64287, Darmstadt, Germany
| | - Nico F. A. van der Vegt
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Center of Smart Interfaces, Technische Universität Darmstadt, Alarich-Weiss-Straße 10, 64287, Darmstadt, Germany
| |
Collapse
|
8
|
Martínez L, Shimizu S. Molecular Interpretation of Preferential Interactions in Protein Solvation: A Solvent-Shell Perspective by Means of Minimum-Distance Distribution Functions. J Chem Theory Comput 2017; 13:6358-6372. [DOI: 10.1021/acs.jctc.7b00599] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Leandro Martínez
- Institute of Chemistry and Center for Computational Engineering & Science, University of Campinas, Campinas, São Paulo 13083-970, Brazil
| | - Seishi Shimizu
- York
Structural Biology Laboratory, Department of Chemistry, University of York, Heslington, York YO10
5DD, U.K
| |
Collapse
|
9
|
Markthaler D, Zeman J, Baz J, Smiatek J, Hansen N. Validation of Trimethylamine-N-oxide (TMAO) Force Fields Based on Thermophysical Properties of Aqueous TMAO Solutions. J Phys Chem B 2017; 121:10674-10688. [DOI: 10.1021/acs.jpcb.7b07774] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Daniel Markthaler
- Institute of Thermodynamics and Thermal Process Engineering and ‡Institute for Computational
Physics, University of Stuttgart, D-70569 Stuttgart, Germany
| | - Johannes Zeman
- Institute of Thermodynamics and Thermal Process Engineering and ‡Institute for Computational
Physics, University of Stuttgart, D-70569 Stuttgart, Germany
| | - Jörg Baz
- Institute of Thermodynamics and Thermal Process Engineering and ‡Institute for Computational
Physics, University of Stuttgart, D-70569 Stuttgart, Germany
| | - Jens Smiatek
- Institute of Thermodynamics and Thermal Process Engineering and ‡Institute for Computational
Physics, University of Stuttgart, D-70569 Stuttgart, Germany
| | - Niels Hansen
- Institute of Thermodynamics and Thermal Process Engineering and ‡Institute for Computational
Physics, University of Stuttgart, D-70569 Stuttgart, Germany
| |
Collapse
|
10
|
Ion-induced alterations of the local hydration environment elucidate Hofmeister effect in a simple classical model of Trp-cage miniprotein. J Mol Model 2017; 23:298. [PMID: 28956172 DOI: 10.1007/s00894-017-3471-0] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2017] [Accepted: 09/10/2017] [Indexed: 10/18/2022]
Abstract
Protein stability is known to be influenced by the presence of Hofmeister active ions in the solution. In addition to direct ion-protein interactions, this influence manifests through the local alterations of the interfacial water structure induced by the anions and cations present in this region. In our earlier works it was pointed out that the effects of Hofmeister active salts on the stability of Trp-cage miniprotein can be modeled qualitatively using non-polarizable force fields. These simulations reproduced the structure-stabilization and structure-destabilization effects of selected kosmotropic and chaotropic salts, respectively. In the present study we use the same model system to elucidate atomic processes behind the chaotropic destabilization and kosmotropic stabilization of the miniprotein. We focus on changes of the local hydration environment of the miniprotein upon addition of NaClO4 and NaF salts to the solution. The process is separated into two parts. In the first, 'promotion' phase, the protein structure is fixed, and the local hydration properties induced by the simultaneous presence of protein and ions are investigated, with a special focus on the interaction of Hofmeister active anions with the charged and polar sites. In the second, 'rearrangement' phase we follow changes of the hydration of ions and the protein, accompanying the conformational relaxation of the protein. We identify significant factors of an enthalpic and entropic nature behind the ion-induced free energy changes of the protein-water system, and also propose a possible atomic mechanism consistent with the Collins's rule, for the chaotropic destabilization and kosmotropic stabilization of protein conformation.
Collapse
|
11
|
Shimizu S, Smith PE. How Osmolytes Counteract Pressure Denaturation on a Molecular Scale. Chemphyschem 2017; 18:2243-2249. [PMID: 28678423 PMCID: PMC5626881 DOI: 10.1002/cphc.201700503] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2017] [Revised: 05/30/2017] [Indexed: 11/08/2022]
Abstract
Life in the deep sea exposes enzymes to high hydrostatic pressure, which decreases their stability. For survival, deep sea organisms tend to accumulate various osmolytes, most notably trimethylamine N-oxide used by fish, to counteract pressure denaturation. However, exactly how these osmolytes work remains unclear. Here, a rigorous statistical thermodynamics approach is used to clarify the mechanism of osmoprotection. It is shown that the weak, nonspecific, and dynamic interactions of water and osmolytes with proteins can be characterized only statistically, and that the competition between protein-osmolyte and protein-water interactions is crucial in determining conformational stability. Osmoprotection is driven by a stronger exclusion of osmolytes from the denatured protein than from the native conformation, and water distribution has no significant effect on these changes at low osmolyte concentrations.
Collapse
Affiliation(s)
- Seishi Shimizu
- York Structural Biology Laboratory, Department of Chemistry, University of York, Heslington, York, YO10 5DD, UK
| | - Paul E Smith
- Department of Chemistry, Kansas State University, 213 CBC Building, Manhattan, Kansas, 66506-0401, USA
| |
Collapse
|
12
|
Ploetz EA, Smith PE. Simulated pressure denaturation thermodynamics of ubiquitin. Biophys Chem 2017; 231:135-145. [PMID: 28576277 DOI: 10.1016/j.bpc.2017.04.006] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2016] [Revised: 04/14/2017] [Accepted: 04/17/2017] [Indexed: 01/09/2023]
Abstract
Simulations of protein thermodynamics are generally difficult to perform and provide limited information. It is desirable to increase the degree of detail provided by simulation and thereby the potential insight into the thermodynamic properties of proteins. In this study, we outline how to analyze simulation trajectories to decompose conformation-specific, parameter free, thermodynamically defined protein volumes into residue-based contributions. The total volumes are obtained using established methods from Fluctuation Solution Theory, while the volume decomposition is new and is performed using a simple proximity method. Native and fully extended ubiquitin are used as the test conformations. Changes in the protein volumes are then followed as a function of pressure, allowing for conformation-specific protein compressibility values to also be obtained. Residue volume and compressibility values indicate significant contributions to protein denaturation thermodynamics from nonpolar and coil residues, together with a general negative compressibility exhibited by acidic residues.
Collapse
Affiliation(s)
- Elizabeth A Ploetz
- Department of Chemistry, 213 CBC Building, 1212 Mid Campus Dr. North, Kansas State University, Manhattan, KS 66506-0401, United States
| | - Paul E Smith
- Department of Chemistry, 213 CBC Building, 1212 Mid Campus Dr. North, Kansas State University, Manhattan, KS 66506-0401, United States.
| |
Collapse
|
13
|
Fingerhut BP, Costard R, Elsaesser T. Predominance of short range Coulomb forces in phosphate-water interactions—a theoretical analysis. J Chem Phys 2016. [DOI: 10.1063/1.4962755] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Affiliation(s)
- Benjamin P. Fingerhut
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, D-12489 Berlin, Germany
| | - Rene Costard
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, D-12489 Berlin, Germany
| | - Thomas Elsaesser
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, D-12489 Berlin, Germany
| |
Collapse
|
14
|
Chiba S, Furuta T, Shimizu S. Kirkwood–Buff Integrals for Aqueous Urea Solutions Based upon the Quantum Chemical Electrostatic Potential and Interaction Energies. J Phys Chem B 2016; 120:7714-23. [DOI: 10.1021/acs.jpcb.6b05611] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Affiliation(s)
- Shuntaro Chiba
- Education
Academy of Computational Life Sciences, Tokyo Institute of Technology, J3-141 4259 Nagatsuta-cho, Midori-ku, Yokohama 226-8501, Japan
| | - Tadaomi Furuta
- School
of Life Science and Technology, Tokyo Institute of Technology, B-62 4259
Nagatsuta-cho, Midori-ku, Yokohama 226-8501, Japan
| | - Seishi Shimizu
- York
Structural Biology Laboratory, Department of Chemistry, University of York, Heslington, York YO10 5DD, United Kingdom
| |
Collapse
|