1
|
Chialvo AA. Preferential Solvation Phenomena as Solute-Induced Structure-Making/Breaking Processes: Linking Thermodynamic Preferential Interaction Parameters to Fundamental Structure Making/Breaking Functions. J Phys Chem B 2024; 128:5228-5245. [PMID: 38754065 DOI: 10.1021/acs.jpcb.4c00385] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/18/2024]
Abstract
In this work, we identify the explicit macroscopic-to-microscopic rigorous links between existing thermodynamic preferential interaction parameters Γ Q α Q β ( χ i ) and microstructural descriptors based on total correlation function integrals, leading to their unambiguous characterization in terms of fundamental structure making/breaking functions S α β . First, we provide the statistics mechanical framework to identify a universal molecular-based signature for the preferential solvation P S phenomenon involving solutes at infinite dilution in mixed-solvent environments and discuss its fundamental properties. Then, we characterize the S α β functions relevant to the P S process, identify the microscopic markers for the existing preferential interaction parameters Γ Q α Q β ( χ i ) in terms of the S α β functions, and test their compliance with a pair of essential microstructural constraints linked to the properties of the universal P S signature. Moreover, we illustrate the analysis by probing the behavior of a representative ternary system comprising the solubility of methane in aqueous 1,4-dioxane mixed-solvent environments under ambient conditions. Finally, we discuss some relevant issues surrounding the statistical mechanical (microstructural) interpretation of the thermodynamic (macroscopic) preferential interaction parameters, review some pitfalls in their evaluation from molecular simulation, and provide an outlook.
Collapse
Affiliation(s)
- Ariel A Chialvo
- Retired Scientist, Knoxville, Tennessee 37922-3108, United States
| |
Collapse
|
2
|
Sreenivasan S, Heffren P, Suh K, Rodnin MV, Kosa E, Fenton AW, Ladokhin AS, Smith PE, Fontes JD, Swint‐Kruse L. The intrinsically disordered transcriptional activation domain of CIITA is functionally tuneable by single substitutions: An exception or a new paradigm? Protein Sci 2024; 33:e4863. [PMID: 38073129 PMCID: PMC10806935 DOI: 10.1002/pro.4863] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2023] [Revised: 12/04/2023] [Accepted: 12/07/2023] [Indexed: 01/27/2024]
Abstract
During protein evolution, some amino acid substitutions modulate protein function ("tuneability"). In most proteins, the tuneable range is wide and can be sampled by a set of protein variants that each contains multiple amino acid substitutions. In other proteins, the full tuneable range can be accessed by a set of variants that each contains a single substitution. Indeed, in some globular proteins, the full tuneable range can be accessed by the set of site-saturating substitutions at an individual "rheostat" position. However, in proteins with intrinsically disordered regions (IDRs), most functional studies-which would also detect tuneability-used multiple substitutions or small deletions. In disordered transcriptional activation domains (ADs), studies with multiple substitutions led to the "acidic exposure" model, which does not anticipate the existence of rheostat positions. In the few studies that did assess effects of single substitutions on AD function, results were mixed: the ADs of two full-length transcription factors did not show tuneability, whereas a fragment of a third AD was tuneable by single substitutions. In this study, we tested tuneability in the AD of full-length human class II transactivator (CIITA). Sequence analyses and experiments showed that CIITA's AD is an IDR. Functional assays of singly-substituted AD variants showed that CIITA's function was highly tuneable, with outcomes not predicted by the acidic exposure model. Four tested positions showed rheostat behavior for transcriptional activation. Thus, tuneability of different IDRs can vary widely. Future studies are needed to illuminate the biophysical features that govern whether an IDR is tuneable by single substitutions.
Collapse
Affiliation(s)
- Shwetha Sreenivasan
- Department of Biochemistry and Molecular BiologyUniversity of Kansas Medical CenterKansas CityKansasUSA
| | - Paul Heffren
- Department of Biochemistry and Molecular BiologyUniversity of Kansas Medical CenterKansas CityKansasUSA
- Present address:
Department of BiosciencesKansas City UniversityKansas CityMissouriUSA
| | - Kyung‐Shin Suh
- Department of ChemistryKansas State UniversityManhattanKansasUSA
| | - Mykola V. Rodnin
- Department of Biochemistry and Molecular BiologyUniversity of Kansas Medical CenterKansas CityKansasUSA
| | - Edina Kosa
- Department of Biochemistry and Molecular BiologyUniversity of Kansas Medical CenterKansas CityKansasUSA
| | - Aron W. Fenton
- Department of Biochemistry and Molecular BiologyUniversity of Kansas Medical CenterKansas CityKansasUSA
| | - Alexey S. Ladokhin
- Department of Biochemistry and Molecular BiologyUniversity of Kansas Medical CenterKansas CityKansasUSA
| | - Paul E. Smith
- Department of ChemistryKansas State UniversityManhattanKansasUSA
| | - Joseph D. Fontes
- Department of Biochemistry and Molecular BiologyUniversity of Kansas Medical CenterKansas CityKansasUSA
| | - Liskin Swint‐Kruse
- Department of Biochemistry and Molecular BiologyUniversity of Kansas Medical CenterKansas CityKansasUSA
| |
Collapse
|
3
|
Kariyawasam NL, Ploetz EA, Swint-Kruse L, Smith PE. Simulated pressure changes in LacI suggest a link between hydration and functional conformational changes. Biophys Chem 2024; 304:107126. [PMID: 37924711 PMCID: PMC10842697 DOI: 10.1016/j.bpc.2023.107126] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2023] [Revised: 10/18/2023] [Accepted: 10/25/2023] [Indexed: 11/06/2023]
Abstract
The functions of many proteins are associated with interconversions among conformational substates. However, these substates can be difficult to measure experimentally, and determining contributions from hydration changes can be especially difficult. Here, we assessed the use of pressure perturbations to sample the substates accessible to the Escherichia coli lactose repressor protein (LacI) in various liganded forms. In the presence of DNA, the regulatory domain of LacI adopts an Open conformation that, in the absence of DNA, changes to a Closed conformation. Increasing the simulation pressure prevented the transition from an Open to a Closed conformation, in a similar manner to the binding of DNA and anti-inducer, ONPF. The results suggest the hydration of specific residues play a significant role in determining the population of different LacI substates and that simulating pressure perturbation could be useful for assessing the role of hydration changes that accompany functionally-relevant amino acid substitutions.
Collapse
Affiliation(s)
- Nilusha L Kariyawasam
- Department of Chemistry, 213 CBC Building, 1212 Mid-Campus Dr. North, Kansas State University, Manhattan, KS 66506, USA
| | - Elizabeth A Ploetz
- Department of Chemistry, 213 CBC Building, 1212 Mid-Campus Dr. North, Kansas State University, Manhattan, KS 66506, USA
| | - Liskin Swint-Kruse
- Department of Biochemistry and Molecular Biology, MSN 3030, 3901 Rainbow Blvd, The University of Kansas Medical Center, Kansas City, KS 66160, USA
| | - Paul E Smith
- Department of Chemistry, 213 CBC Building, 1212 Mid-Campus Dr. North, Kansas State University, Manhattan, KS 66506, USA.
| |
Collapse
|
4
|
Zhang P, Yang W. Toward a general neural network force field for protein simulations: Refining the intramolecular interaction in protein. J Chem Phys 2023; 159:024118. [PMID: 37431910 PMCID: PMC10481389 DOI: 10.1063/5.0142280] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2023] [Accepted: 06/22/2023] [Indexed: 07/12/2023] Open
Abstract
Molecular dynamics (MD) is an extremely powerful, highly effective, and widely used approach to understanding the nature of chemical processes in atomic details for proteins. The accuracy of results from MD simulations is highly dependent on force fields. Currently, molecular mechanical (MM) force fields are mainly utilized in MD simulations because of their low computational cost. Quantum mechanical (QM) calculation has high accuracy, but it is exceedingly time consuming for protein simulations. Machine learning (ML) provides the capability for generating accurate potential at the QM level without increasing much computational effort for specific systems that can be studied at the QM level. However, the construction of general machine learned force fields, needed for broad applications and large and complex systems, is still challenging. Here, general and transferable neural network (NN) force fields based on CHARMM force fields, named CHARMM-NN, are constructed for proteins by training NN models on 27 fragments partitioned from the residue-based systematic molecular fragmentation (rSMF) method. The NN for each fragment is based on atom types and uses new input features that are similar to MM inputs, including bonds, angles, dihedrals, and non-bonded terms, which enhance the compatibility of CHARMM-NN to MM MD and enable the implementation of CHARMM-NN force fields in different MD programs. While the main part of the energy of the protein is based on rSMF and NN, the nonbonded interactions between the fragments and with water are taken from the CHARMM force field through mechanical embedding. The validations of the method for dipeptides on geometric data, relative potential energies, and structural reorganization energies demonstrate that the CHARMM-NN local minima on the potential energy surface are very accurate approximations to QM, showing the success of CHARMM-NN for bonded interactions. However, the MD simulations on peptides and proteins indicate that more accurate methods to represent protein-water interactions in fragments and non-bonded interactions between fragments should be considered in the future improvement of CHARMM-NN, which can increase the accuracy of approximation beyond the current mechanical embedding QM/MM level.
Collapse
Affiliation(s)
- Pan Zhang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| |
Collapse
|
5
|
Panigrahy S, Sahu R, Reddy SK, Nayar D. Structure, energetics and dynamics in crowded amino acid solutions: a molecular dynamics study. Phys Chem Chem Phys 2023; 25:5430-5442. [PMID: 36744506 DOI: 10.1039/d2cp04238j] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
A comprehensive understanding of crowding effects on biomolecular processes necessitates investigating the bulk thermodynamic and kinetic properties of the solutions with an accurate molecular representation of the crowded milieu. Recent studies have reparameterized the non-bonded dispersion interaction of solutes to precisely model intermolecular interactions, which would circumvent artificial aggregation as shown by the original force-fields. However, the performance of this reparameterization is yet to be assessed for concentrated crowded solutions in terms of investigating the hydration shell structure, energetics and dynamics. In this study, we perform molecular dynamics simulations of crowded aqueous solutions of five zwitterionic neutral amino acids (Gly, Ala, Thr, Pro, and Ser), mimicking the molecular crowding environment, using a modified AMBER ff99SB-ILDN force-field. We systematically examine and show that the reproducibility of the osmotic coefficients, density, viscosity and self-diffusivity of amino acids improves using the modified force-field in crowded concentrations. The modified force-field also improves the structuring of the solute solvation shells, solute interaction energy and convergence of tails of radial distribution functions, indicating reduction in the artificial aggregation. Our results also indicate that the hydrogen bonding network of water weakens and water molecules anomalously diffuse at small time scales in the crowded solutions. These results underscore the significance of examining the solution properties and anomalous hydration behaviour of water in crowded solutions, which have implications in shaping the structure and dynamics of biomolecules. The findings also illustrate the improvement in predicting bulk solution properties using the modified force-field, thereby providing an approach towards accurate modeling of crowded molecular solutions.
Collapse
Affiliation(s)
- Sibasankar Panigrahy
- Department of Materials Science and Engineering, Indian Institute of Technology Delhi, New Delhi 110016, India.
| | - Rahul Sahu
- Center for Computational and Data Sciences, Indian Institute of Technology Kharagpur, West Bengal 721302, India
| | - Sandeep K Reddy
- Center for Computational and Data Sciences, Indian Institute of Technology Kharagpur, West Bengal 721302, India
| | - Divya Nayar
- Department of Materials Science and Engineering, Indian Institute of Technology Delhi, New Delhi 110016, India.
| |
Collapse
|
6
|
McIvor JAP, Larsen DS, Mercadante D. Simulating Polyproline II-Helix-Rich Peptides with the Latest Kirkwood-Buff Force Field: A Direct Comparison with AMBER and CHARMM. J Phys Chem B 2022; 126:7833-7846. [PMID: 36125334 DOI: 10.1021/acs.jpcb.2c03837] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We simulated the dynamics of a set of peptides characterized by ensembles rich in PPII-helical content, to assess the ability of the most recent Kirkwood-Buff force field (KBFF20) to sample this conformational peculiarity. KBFF has been previously shown to capably reproduce experimental dimensions of disordered proteins, while being limited in confidently sampling structured proteins. Further development of the force field bridged this gap. It is however still unknown what are the main differences between KBFF and AMBER/CHARMM force fields. A direct comparison is now possible as both AMBER/CHARMM force fields have been used to sample peptides rich in PPII-helical content. We found that KBFF20 samples' PPII-helical content qualitatively matches both AMBER and CHARMM force fields, with the main difference being the KBFF ability to populate the αR region of the Ramachandran plot in the set of simulated peptides. Overall, KBFF20 is a well-balanced force field, able to sample the dynamics of both structured and unstructured proteins.
Collapse
Affiliation(s)
- Jordan A P McIvor
- School of Chemical Sciences, The University of Auckland, 23 Symonds Street, Auckland 1010, New Zealand
| | - Danaé S Larsen
- School of Chemical Sciences, The University of Auckland, 23 Symonds Street, Auckland 1010, New Zealand
| | - Davide Mercadante
- School of Chemical Sciences, The University of Auckland, 23 Symonds Street, Auckland 1010, New Zealand
| |
Collapse
|
7
|
Ganguly P, Bubák D, Polák J, Fagan P, Dračínský M, van der Vegt NFA, Heyda J, Shea JE. Cosolvent Exclusion Drives Protein Stability in Trimethylamine N-Oxide and Betaine Solutions. J Phys Chem Lett 2022; 13:7980-7986. [PMID: 35984361 DOI: 10.1021/acs.jpclett.2c01692] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Using a combination of molecular dynamics simulation, dialysis experiments, and electronic circular dichroism measurements, we studied the solvation thermodynamics of proteins in two osmolyte solutions, trimethylamine N-oxide (TMAO) and betaine. We showed that existing force fields are unable to capture the solvation properties of the proteins lysozyme and ribonuclease T1 and that the inaccurate parametrization of protein-osmolyte interactions in these force fields promoted an unphysical strong thermal denaturation of the trpcage protein. We developed a novel force field for betaine (the KBB force field) which reproduces the experimental solution Kirkwood-Buff integrals and density. We further introduced appropriate scaling to protein-osmolyte interactions in both the betaine and TMAO force fields which led to successful reproduction of experimental protein-osmolyte preferential binding coefficients for lysozyme and ribonuclease T1 and prevention of the unphysical denaturation of trpcage in osmolyte solutions. Correct parametrization of protein-TMAO interactions also led to the stabilization of the collapsed conformations of a disordered elastin-like peptide, while the uncorrected parameters destabilized the collapsed structures. Our results establish that the thermodynamic stability of proteins in both betaine and TMAO solutions is governed by osmolyte exclusion from proteins.
Collapse
Affiliation(s)
- Pritam Ganguly
- Department of Chemistry and Biochemistry, University of California at Santa Barbara, Santa Barbara, California93106, United States
| | - Dominik Bubák
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 166 28Prague 6, Czech Republic
| | - Jakub Polák
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 166 28Prague 6, Czech Republic
| | - Patrik Fagan
- Department of Analytical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 166 28Prague 6, Czech Republic
| | - Martin Dračínský
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nám. 2, 160 00Prague, Czech Republic
| | - Nico F A van der Vegt
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, Alarich-Weiss-Straße 10, Darmstadt64287, Germany
| | - Jan Heyda
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 166 28Prague 6, Czech Republic
| | - Joan-Emma Shea
- Department of Chemistry and Biochemistry, University of California at Santa Barbara, Santa Barbara, California93106, United States
- Department of Physics, University of California at Santa Barbara, Santa Barbara, California93106, United States
| |
Collapse
|
8
|
Venetsanos F, Anogiannakis SD, Theodorou DN. Mixing Thermodynamics and Flory–Huggins Interaction Parameter of Polyethylene Oxide/Polyethylene Oligomeric Blends from Kirkwood–Buff Theory and Molecular Simulations. Macromolecules 2022. [DOI: 10.1021/acs.macromol.2c00642] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Fotis Venetsanos
- School of Chemical Engineering, National Technical University of Athens, 9 Heroon Polytechniou Street, Athens 15780, Greece
| | - Stefanos D. Anogiannakis
- School of Chemical Engineering, National Technical University of Athens, 9 Heroon Polytechniou Street, Athens 15780, Greece
| | - Doros N. Theodorou
- School of Chemical Engineering, National Technical University of Athens, 9 Heroon Polytechniou Street, Athens 15780, Greece
| |
Collapse
|
9
|
Boothroyd S, Madin OC, Mobley DL, Wang LP, Chodera JD, Shirts MR. Improving Force Field Accuracy by Training against Condensed-Phase Mixture Properties. J Chem Theory Comput 2022; 18:3577-3592. [PMID: 35533269 DOI: 10.1021/acs.jctc.1c01268] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Developing a sufficiently accurate classical force field representation of molecules is key to realizing the full potential of molecular simulations as a route to gaining a fundamental insight into a broad spectrum of chemical and biological phenomena. This is only possible, however, if the many complex interactions between molecules of different species in the system are accurately captured by the model. Historically, the intermolecular van der Waals (vdW) interactions have primarily been trained against densities and enthalpies of vaporization of pure (single-component) systems, with occasional usage of hydration free energies. In this study, we demonstrate how including physical property data of binary mixtures can better inform these parameters, encoding more information about the underlying physics of the system in complex chemical mixtures. To demonstrate this, we retrain a select number of Lennard-Jones parameters describing the vdW interactions of the OpenFF 1.0.0 (Parsley) fixed charge force field against training sets composed of densities and enthalpies of mixing for binary liquid mixtures as well as densities and enthalpies of vaporization of pure liquid systems and assess the performance of each of these combinations. We show that retraining against the mixture data improves the force field's ability to reproduce mixture properties, including solvation free energies, correcting some systematic errors that exist when training vdW interactions against properties of pure systems only.
Collapse
Affiliation(s)
- Simon Boothroyd
- Boothroyd Scientific Consulting Ltd., 71-75 Shelton Street, London WC2H 9JQ, Greater London, U.K
| | - Owen C Madin
- Department of Chemical and Biological Engineering, University of Colorado Boulder, Boulder, Colorado 80309, United States
| | - David L Mobley
- Department of Chemistry, University of California, Irvine, California 92617, United States.,Department of Pharmaceutical Sciences, University of California, Irvine, California 92617, United States
| | - Lee-Ping Wang
- Department of Chemistry, University of California, Davis, California 95616, United States
| | - John D Chodera
- Computational & Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, New York 10065, United States
| | - Michael R Shirts
- Department of Chemical and Biological Engineering, University of Colorado Boulder, Boulder, Colorado 80309, United States
| |
Collapse
|
10
|
Heo L, Sugita Y, Feig M. Protein assembly and crowding simulations. Curr Opin Struct Biol 2022; 73:102340. [PMID: 35219215 PMCID: PMC8957576 DOI: 10.1016/j.sbi.2022.102340] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2021] [Revised: 01/07/2022] [Accepted: 01/18/2022] [Indexed: 11/17/2022]
Abstract
Proteins encounter frequent molecular interactions in biological environments. Computer simulations have become an increasingly important tool in providing mechanistic insights into how such interactions in vivo relate to their biological function. The review here focuses on simulations describing protein assembly and molecular crowding effects as two important aspects that are distinguished mainly by how specific and long-lived protein contacts are. On the topic of crowding, recent simulations have provided new insights into how crowding affects protein folding and stability, modulates enzyme activity, and affects diffusive properties. Recent studies of assembly processes focus on assembly pathways, especially for virus capsids, amyloid aggregation pathways, and the role of multivalent interactions leading to phase separation. Also, discussed are technical challenges in achieving increasingly realistic simulations of interactions in cellular environments.
Collapse
Affiliation(s)
- Lim Heo
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI 48824, USA. https://twitter.com/huhlim
| | - Yuji Sugita
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan; Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, 6-7-1 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan; Computational Biophysics Research Team, RIKEN Center for Computational Science, 6-7-1 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.
| | - Michael Feig
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI 48824, USA.
| |
Collapse
|
11
|
Rivas G, Minton A. Influence of Nonspecific Interactions on Protein Associations: Implications for Biochemistry In Vivo. Annu Rev Biochem 2022; 91:321-351. [PMID: 35287477 DOI: 10.1146/annurev-biochem-040320-104151] [Citation(s) in RCA: 17] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
The cellular interior is composed of a variety of microenvironments defined by distinct local compositions and composition-dependent intermolecular interactions. We review the various types of nonspecific interactions between proteins and between proteins and other macromolecules and supramolecular structures that influence the state of association and functional properties of a given protein existing within a particular microenvironment at a particular point in time. The present state of knowledge is summarized, and suggestions for fruitful directions of research are offered. Expected final online publication date for the Annual Review of Biochemistry, Volume 91 is June 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.
Collapse
Affiliation(s)
- Germán Rivas
- Department of Structural and Chemical Biology, Centro de Investigaciones Biológicas Margarita Salas, Consejo Superior de Investigaciones Científicas, Madrid, Spain;
| | - Allen Minton
- Laboratory of Biochemistry and Genetics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland, USA;
| |
Collapse
|
12
|
Biomolecular Simulations with the Three-Dimensional Reference Interaction Site Model with the Kovalenko-Hirata Closure Molecular Solvation Theory. Int J Mol Sci 2021; 22:ijms22105061. [PMID: 34064655 PMCID: PMC8151972 DOI: 10.3390/ijms22105061] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2021] [Revised: 05/08/2021] [Accepted: 05/10/2021] [Indexed: 11/17/2022] Open
Abstract
The statistical mechanics-based 3-dimensional reference interaction site model with the Kovalenko-Hirata closure (3D-RISM-KH) molecular solvation theory has proven to be an essential part of a multiscale modeling framework, covering a vast region of molecular simulation techniques. The successful application ranges from the small molecule solvation energy to the bulk phase behavior of polymers, macromolecules, etc. The 3D-RISM-KH successfully predicts and explains the molecular mechanisms of self-assembly and aggregation of proteins and peptides related to neurodegeneration, protein-ligand binding, and structure-function related solvation properties. Upon coupling the 3D-RISM-KH theory with a novel multiple time-step molecular dynamic (MD) of the solute biomolecule stabilized by the optimized isokinetic Nosé-Hoover chain thermostat driven by effective solvation forces obtained from 3D-RISM-KH and extrapolated forward by generalized solvation force extrapolation (GSFE), gigantic outer time-steps up to picoseconds to accurately calculate equilibrium properties were obtained in this new quasidynamics protocol. The multiscale OIN/GSFE/3D-RISM-KH algorithm was implemented in the Amber package and well documented for fully flexible model of alanine dipeptide, miniprotein 1L2Y, and protein G in aqueous solution, with a solvent sampling rate ~150 times faster than a standard MD simulation in explicit water. Further acceleration in computation can be achieved by modifying the extent of solvation layers considered in the calculation, as well as by modifying existing closure relations. This enhanced simulation technique has proven applications in protein-ligand binding energy calculations, ligand/solvent binding site prediction, molecular solvation energy calculations, etc. Applications of the RISM-KH theory in molecular simulation are discussed in this work.
Collapse
|
13
|
Ploetz EA, Karunaweera S, Smith PE. Kirkwood-Buff-Derived Force Field for Peptides and Proteins: Applications of KBFF20. J Chem Theory Comput 2021; 17:2991-3009. [PMID: 33878264 DOI: 10.1021/acs.jctc.1c00076] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Here, we perform structural, thermodynamic, and kinetics tests of the Kirkwood-Buff-derived force field, KBFF20, for peptides and proteins developed in the previous article. The physical/structural tests measure the ability of KBFF20 to capture the experimental J-couplings for small peptides, to keep globular monomeric and oligomeric proteins folded, and to produce the experimentally relevant expanded conformational ensembles of intrinsically disordered proteins. The thermodynamic-based tests probe KBFF20's ability to quantify the preferential interactions of sodium chloride around native β-lactoglobulin and urea around native lysozyme, to reproduce the melting curves for small helix- and sheet-based peptides, and to fold the small proteins Trp-cage and Villin. The kinetics-based tests quantify how well KBFF20 can match the experimental contact formation rates of small, repeat-sequence peptides of variable lengths and the rotational diffusion coefficients of globular proteins. The results suggest that KBFF20 is naturally able to reproduce properties of both folded and disordered proteins, which we attribute to the use of the Kirkwood-Buff theory as the foundation of the force field's development. However, we show that KBFF20 tends to lose some well-defined secondary structural elements and increases the percentage of coil regions, indicating that the perfect balance of all interactions remains elusive. Nevertheless, we argue that KBFF20 is an improvement over recently modified force fields that require ad hoc interventions to prevent the collapse of intrinsically disordered proteins.
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
| | - Paul E Smith
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| |
Collapse
|