1
|
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
|
2
|
Jussupow A, Messias AC, Stehle R, Geerlof A, Solbak SMØ, Paissoni C, Bach A, Sattler M, Camilloni C. The dynamics of linear polyubiquitin. SCIENCE ADVANCES 2020; 6:6/42/eabc3786. [PMID: 33055165 PMCID: PMC7556843 DOI: 10.1126/sciadv.abc3786] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/21/2020] [Accepted: 08/25/2020] [Indexed: 05/17/2023]
Abstract
Polyubiquitin chains are flexible multidomain proteins, whose conformational dynamics enable them to regulate multiple biological pathways. Their dynamic is determined by the linkage between ubiquitins and by the number of ubiquitin units. Characterizing polyubiquitin behavior as a function of their length is hampered because of increasing system size and conformational variability. Here, we introduce a new approach to efficiently integrating small-angle x-ray scattering with simulations allowing us to accurately characterize the dynamics of linear di-, tri-, and tetraubiquitin in the free state as well as of diubiquitin in complex with NEMO, a central regulator in the NF-κB pathway. Our results show that the behavior of the diubiquitin subunits is independent of the presence of additional ubiquitin modules and that the dynamics of polyubiquitins with different lengths follow a simple model. Together with experimental data from multiple biophysical techniques, we then rationalize the 2:1 NEMO:polyubiquitin binding.
Collapse
Affiliation(s)
- Alexander Jussupow
- Department of Chemistry and Institute for Advanced Study, Technical University of Munich, Garching 85747, Germany
| | - Ana C Messias
- Institute of Structural Biology, Helmholtz Zentrum München, Neuherberg 85764, Germany
- Center for Integrated Protein Science Munich at Department of Chemistry, Technical University of Munich, Garching 85747, Germany
| | - Ralf Stehle
- Institute of Structural Biology, Helmholtz Zentrum München, Neuherberg 85764, Germany
- Center for Integrated Protein Science Munich at Department of Chemistry, Technical University of Munich, Garching 85747, Germany
| | - Arie Geerlof
- Institute of Structural Biology, Helmholtz Zentrum München, Neuherberg 85764, Germany
- Center for Integrated Protein Science Munich at Department of Chemistry, Technical University of Munich, Garching 85747, Germany
| | - Sara M Ø Solbak
- Department of Drug Design and Pharmacology, Faculty of Health and Medical Sciences, University of Copenhagen, Universitetsparken 2, DK-2100 Copenhagen, Denmark
| | - Cristina Paissoni
- Dipartimento di Bioscienze, Università degli studi di Milano, 20133 Milano, Italy
| | - Anders Bach
- Department of Drug Design and Pharmacology, Faculty of Health and Medical Sciences, University of Copenhagen, Universitetsparken 2, DK-2100 Copenhagen, Denmark
| | - Michael Sattler
- Institute of Structural Biology, Helmholtz Zentrum München, Neuherberg 85764, Germany.
- Center for Integrated Protein Science Munich at Department of Chemistry, Technical University of Munich, Garching 85747, Germany
| | - Carlo Camilloni
- Department of Chemistry and Institute for Advanced Study, Technical University of Munich, Garching 85747, Germany.
- Dipartimento di Bioscienze, Università degli studi di Milano, 20133 Milano, Italy
| |
Collapse
|
3
|
Miller MS, Lay WK, Li S, Hacker WC, An J, Ren J, Elcock AH. Reparametrization of Protein Force Field Nonbonded Interactions Guided by Osmotic Coefficient Measurements from Molecular Dynamics Simulations. J Chem Theory Comput 2017; 13:1812-1826. [PMID: 28296391 PMCID: PMC5543770 DOI: 10.1021/acs.jctc.6b01059] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
There is a small, but growing, body of literature describing the use of osmotic coefficient measurements to validate and reparametrize simulation force fields. Here we have investigated the ability of five very commonly used force field and water model combinations to reproduce the osmotic coefficients of seven neutral amino acids and five small molecules. The force fields tested include AMBER ff99SB-ILDN, CHARMM36, GROMOS54a7, and OPLS-AA, with the first of these tested in conjunction with the TIP3P and TIP4P-Ew water models. In general, for both the amino acids and the small molecules, the tested force fields produce computed osmotic coefficients that are lower than experiment; this is indicative of excessively favorable solute-solute interactions. The sole exception to this general trend is provided by GROMOS54a7 when applied to amino acids: in this case, the computed osmotic coefficients are consistently too high. Importantly, we show that all of the force fields tested can be made to accurately reproduce the experimental osmotic coefficients of the amino acids when minor modifications-some previously reported by others and some that are new to this study-are made to the van der Waals interactions of the charged terminal groups. Special care is required, however, when simulating Proline with a number of the force fields, and a hydroxyl-group specific modification is required in order to correct Serine and Threonine when simulated with AMBER ff99SB-ILDN. Interestingly, an alternative parametrization of the van der Waals interactions in the latter force field, proposed by the Nerenberg and Head-Gordon groups, is shown to immediately produce osmotic coefficients that are in excellent agreement with experiment. Overall, this study reinforces the idea that osmotic coefficient measurements can be used to identify general shortcomings in commonly used force fields' descriptions of solute-solute interactions and further demonstrates that modifications to van der Waals parameters provide a simple route to optimizing agreement with experiment.
Collapse
Affiliation(s)
- Mark S. Miller
- Department of Biochemistry, University of Iowa, Iowa City, IA 52242
| | - Wesley K. Lay
- Department of Biochemistry, University of Iowa, Iowa City, IA 52242
| | - Shuxiang Li
- Department of Biochemistry, University of Iowa, Iowa City, IA 52242
| | | | - Jiadi An
- Department of Biochemistry, University of Iowa, Iowa City, IA 52242
| | - Jianlan Ren
- Department of Biochemistry, University of Iowa, Iowa City, IA 52242
| | - Adrian H. Elcock
- Department of Biochemistry, University of Iowa, Iowa City, IA 52242
| |
Collapse
|
4
|
Nguyen HT, Pabit SA, Pollack L, Case DA. Extracting water and ion distributions from solution x-ray scattering experiments. J Chem Phys 2017; 144:214105. [PMID: 27276943 DOI: 10.1063/1.4953037] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
Small-angle X-ray scattering measurements can provide valuable information about the solvent environment around biomolecules, but it can be difficult to extract solvent-specific information from observed intensity profiles. Intensities are proportional to the square of scattering amplitudes, which are complex quantities. Amplitudes in the forward direction are real, and the contribution from a solute of known structure (and from the waters it excludes) can be estimated from theory; hence, the amplitude arising from the solvent environment can be computed by difference. We have found that this "square root subtraction scheme" can be extended to non-zero q values, out to 0.1 Å(-1) for the systems considered here, since the phases arising from the solute and from the water environment are nearly identical in this angle range. This allows us to extract aspects of the water and ion distributions (beyond their total numbers), by combining experimental data for the complete system with calculations for the solutes. We use this approach to test molecular dynamics and integral-equation (3D-RISM (three-dimensional reference interaction site model)) models for solvent structure around myoglobin, lysozyme, and a 25 base-pair duplex DNA. Comparisons can be made both in Fourier space and in terms of the distribution of interatomic distances in real space. Generally, computed solvent distributions arising from the MD simulations fit experimental data better than those from 3D-RISM, even though the total small-angle X-ray scattering patterns are very similar; this illustrates the potential power of this sort of analysis to guide the development of computational models.
Collapse
Affiliation(s)
- Hung T Nguyen
- BioMaPS Institute for Quantitative Biology, Rutgers University, Piscataway, New Jersey 08854, USA
| | - Suzette A Pabit
- School of Applied and Engineering Physics, Cornell University, Ithaca, New York 14853, USA
| | - Lois Pollack
- School of Applied and Engineering Physics, Cornell University, Ithaca, New York 14853, USA
| | - David A Case
- Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, USA
| |
Collapse
|
5
|
Matthes D, Gapsys V, Brennecke JT, de Groot BL. An Atomistic View of Amyloidogenic Self-assembly: Structure and Dynamics of Heterogeneous Conformational States in the Pre-nucleation Phase. Sci Rep 2016; 6:33156. [PMID: 27616019 PMCID: PMC5018807 DOI: 10.1038/srep33156] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2016] [Accepted: 08/22/2016] [Indexed: 01/01/2023] Open
Abstract
The formation of well-defined filamentous amyloid structures involves a polydisperse collection of oligomeric states for which relatively little is known in terms of structural organization. Here we use extensive, unbiased explicit solvent molecular dynamics (MD) simulations to investigate the structural and dynamical features of oligomeric aggregates formed by a number of highly amyloidogenic peptides at atomistic resolution on the μs time scale. A consensus approach has been adopted to analyse the simulations in multiple force fields, yielding an in-depth characterization of pre-fibrillar oligomers and their global and local structure properties. A collision cross section analysis revealed structurally heterogeneous aggregate ensembles for the individual oligomeric states that lack a single defined quaternary structure during the pre-nucleation phase. To gain insight into the conformational space sampled in early aggregates, we probed their substructure and found emerging β-sheet subunit layers and a multitude of ordered intermolecular β-structure motifs with growing aggregate size. Among those, anti-parallel out-of-register β-strands compatible with toxic β-barrel oligomers were particularly prevalent already in smaller aggregates and formed prior to ordered fibrillar structure elements. Notably, also distinct fibril-like conformations emerged in the oligomeric state and underscore the notion that pre-nucleated oligomers serve as a critical intermediate step on-pathway to fibrils.
Collapse
Affiliation(s)
- Dirk Matthes
- Computational Biomolecular Dynamics Group, Max Planck Institute for Biophysical Chemistry, Göttingen, Germany
| | - Vytautas Gapsys
- Computational Biomolecular Dynamics Group, Max Planck Institute for Biophysical Chemistry, Göttingen, Germany
| | - Julian T Brennecke
- Computational Biomolecular Dynamics Group, Max Planck Institute for Biophysical Chemistry, Göttingen, Germany
| | - Bert L de Groot
- Computational Biomolecular Dynamics Group, Max Planck Institute for Biophysical Chemistry, Göttingen, Germany
| |
Collapse
|
6
|
Ivanov MV, Talipov MR, Timerghazin QK. Electrostatic point charge fitting as an inverse problem: Revealing the underlying ill-conditioning. J Chem Phys 2015; 143:134102. [DOI: 10.1063/1.4932105] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Affiliation(s)
- Maxim V. Ivanov
- Department of Chemistry, Marquette University, Milwaukee, Wisconsin 53201-1881, USA
| | - Marat R. Talipov
- Department of Chemistry, Marquette University, Milwaukee, Wisconsin 53201-1881, USA
| | - Qadir K. Timerghazin
- Department of Chemistry, Marquette University, Milwaukee, Wisconsin 53201-1881, USA
| |
Collapse
|
7
|
Ivanov MV, Talipov MR, Timerghazin QK. Genetic Algorithm Optimization of Point Charges in Force Field Development: Challenges and Insights. J Phys Chem A 2015; 119:1422-34. [DOI: 10.1021/acs.jpca.5b00218] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Maxim V. Ivanov
- Department of Chemistry, Marquette University,
P.O. Box 1881, Milwaukee, Wisconsin 53201-1881, United States
| | - Marat R. Talipov
- Department of Chemistry, Marquette University,
P.O. Box 1881, Milwaukee, Wisconsin 53201-1881, United States
| | - Qadir K. Timerghazin
- Department of Chemistry, Marquette University,
P.O. Box 1881, Milwaukee, Wisconsin 53201-1881, United States
| |
Collapse
|
8
|
Cerutti DS, Swope WC, Rice J, Case DA. ff14ipq: A Self-Consistent Force Field for Condensed-Phase Simulations of Proteins. J Chem Theory Comput 2014; 10:4515-4534. [PMID: 25328495 PMCID: PMC4196740 DOI: 10.1021/ct500643c] [Citation(s) in RCA: 78] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2014] [Indexed: 01/25/2023]
Abstract
We present the ff14ipq force field, implementing the previously published IPolQ charge set for simulations of complete proteins. Minor modifications to the charge derivation scheme and van der Waals interactions between polar atoms are introduced. Torsion parameters are developed through a generational learning approach, based on gas-phase MP2/cc-pVTZ single-point energies computed of structures optimized by the force field itself rather than the quantum benchmark. In this manner, we sacrifice information about the true quantum minima in order to ensure that the force field maintains optimal agreement with the MP2/cc-pVTZ benchmark for the ensembles it will actually produce in simulations. A means of making the gas-phase torsion parameters compatible with solution-phase IPolQ charges is presented. The ff14ipq model is an alternative to ff99SB and other Amber force fields for protein simulations in programs that accommodate pair-specific Lennard-Jones combining rules. The force field gives strong performance on α-helical and β-sheet oligopeptides as well as globular proteins over microsecond time scale simulations, although it has not yet been tested in conjunction with lipid and nucleic acid models. We show how our choices in parameter development influence the resulting force field and how other choices that may have appeared reasonable would actually have led to poorer results. The tools we developed may also aid in the development of future fixed-charge and even polarizable biomolecular force fields.
Collapse
Affiliation(s)
- David S. Cerutti
- Department
of Chemistry and
Chemical Biology and BioMaPS Institute, Rutgers University, 610 Taylor Road, Piscataway, New Jersey 08854-8066, United States
| | - William C. Swope
- Department
of Chemistry and
Chemical Biology and BioMaPS Institute, Rutgers University, 610 Taylor Road, Piscataway, New Jersey 08854-8066, United States
| | - Julia
E. Rice
- Department
of Chemistry and
Chemical Biology and BioMaPS Institute, Rutgers University, 610 Taylor Road, Piscataway, New Jersey 08854-8066, United States
| | - David A. Case
- Department
of Chemistry and
Chemical Biology and BioMaPS Institute, Rutgers University, 610 Taylor Road, Piscataway, New Jersey 08854-8066, United States
| |
Collapse
|
9
|
Hazel A, Chipot C, Gumbart JC. Thermodynamics of Deca-alanine Folding in Water. J Chem Theory Comput 2014; 10:2836-2844. [PMID: 25061447 PMCID: PMC4095909 DOI: 10.1021/ct5002076] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2014] [Indexed: 11/28/2022]
Abstract
![]()
The determination of the folding
dynamics of polypeptides and proteins
is critical in characterizing their functions in biological systems.
Numerous computational models and methods have been developed for
studying structure formation at the atomic level. Due to its small
size and simple structure, deca-alanine is used as a model system
in molecular dynamics (MD) simulations. The free energy of unfolding
in vacuum has been studied extensively using the end-to-end distance
of the peptide as the reaction coordinate. However, few studies have
been conducted in the presence of explicit solvent. Previous results
show a significant decrease in the free energy of extended conformations
in water, but the α-helical state is still notably favored over
the extended state. Although sufficient in vacuum, we show that end-to-end
distance is incapable of capturing the full complexity of deca-alanine
folding in water. Using α-helical content as a second reaction
coordinate, we deduce a more descriptive free-energy landscape, one
which reveals a second energy minimum in the extended conformations
that is of comparable free energy to the α-helical state. Equilibrium
simulations demonstrate the relative stability of the extended and
α-helical states in water as well as the transition between
the two states. This work reveals both the necessity and challenge
of determining a proper reaction coordinate to fully characterize
a given process.
Collapse
Affiliation(s)
- Anthony Hazel
- School of Physics, Georgia Institute of Technology , Atlanta, Georgia 30332, United States
| | - Christophe Chipot
- Beckman Institute, University of Illinois at Urbana-Champaign , Urbana, Illinois 61801, United States ; Laboratoire Associé International Centre National de la Recherche Scientifique - Unité Mixte de Recherche , N°7565, BP 70239, 54506 Vandœuvre-lès-Nancy, France
| | - James C Gumbart
- School of Physics, Georgia Institute of Technology , Atlanta, Georgia 30332, United States
| |
Collapse
|