1
|
Marien J, Prévost C, Sacquin-Mora S. nP-Collabs: Investigating Counterion-Mediated Bridges in the Multiply Phosphorylated Tau-R2 Repeat. J Chem Inf Model 2024; 64:6570-6582. [PMID: 39092904 DOI: 10.1021/acs.jcim.4c00742] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/04/2024]
Abstract
Tau is an intrinsically disordered (IDP) microtubule-associated protein (MAP) that plays a key part in microtubule assembly and organization. The function of tau can be regulated by multiple phosphorylation sites. These post-translational modifications are known to decrease the binding affinity of tau for microtubules, and abnormal tau phosphorylation patterns are involved in Alzheimer's disease. Using all-atom molecular dynamics simulations, we compared the conformational landscapes explored by the tau R2 repeat domain (which comprises a strong tubulin binding site) in its native state and with multiple phosphorylations on the S285, S289, and S293 residues, with four different standard force field (FF)/water model combinations. We find that the different parameters used for the phosphate groups (which can be more or less flexible) in these FFs and the specific interactions between bulk cations and water lead to the formation of a specific type of counterion bridge, termed nP-collab (for nphosphate collaboration, with n being an integer), where counterions form stable structures binding with two or three phosphate groups simultaneously. The resulting effect of nP-collabs on the tau-R2 conformational space differs when using sodium or potassium cations and is likely to impact the peptide overall dynamics and how this MAP interacts with tubulins. We also investigated the effect of phosphoresidue spacing and ionic concentration by modeling polyalanine peptides containing two phosphoserines located one-six residues apart. Three new metrics specifically tailored for IDPs (proteic Menger curvature, local curvature, and local flexibility) were introduced, which allow us to fully characterize the impact of nP-collabs on the dynamics of disordered peptides at the residue level.
Collapse
Affiliation(s)
- Jules Marien
- Laboratoire de Biochimie Théorique, Université Paris-Cité, CNRS, 13 Rue Pierre et Marie Curie, 75005 Paris, France
| | - Chantal Prévost
- Laboratoire de Biochimie Théorique, Université Paris-Cité, CNRS, 13 Rue Pierre et Marie Curie, 75005 Paris, France
| | - Sophie Sacquin-Mora
- Laboratoire de Biochimie Théorique, Université Paris-Cité, CNRS, 13 Rue Pierre et Marie Curie, 75005 Paris, France
| |
Collapse
|
2
|
Falbo E, Lavecchia A. HessFit: A Toolkit to Derive Automated Force Fields from Quantum Mechanical Information. J Chem Inf Model 2024; 64:5634-5645. [PMID: 38897917 DOI: 10.1021/acs.jcim.4c00540] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/21/2024]
Abstract
In this study, we introduce a novel approach to enhance the accuracy of molecular dynamics simulations by refining the force fields (FFs) through a combination of transferable parameters and molecule-specific characteristics derived from quantum mechanical (QM) calculations. Traditional FFs often prioritize generality over precision, leading to limitations in the accuracy of accurately capturing intra- and intermolecular interactions. To address this, we present an open-source toolkit, called HessFit, designed to integrate QM-derived bonded parameters and atomic charges into existing FFs. In combination with bond, angle, torsional, and nonbonded parameters derivation, HessFit can easily extract multiple barrier terms of dihedrals from QM Hessian and gradient or return all terms through a fitting procedure scheme of QM potential energy surface. We showcase the effectiveness of HessFit through comprehensive evaluations of vibrational properties across a diverse set of small molecules and demonstrate that experimental results support its ability in predicting thermodynamic properties of organic molecules compared to previous state-of-the-art approaches. We further explore its application to Zn2+ metalloprotein models, which are hard systems to treat with automatic approaches. Our results demonstrate that HessFit parameters compete with GAFF2 and OPLS parameters to describing small organic molecules, and its feasibility is also comparable to current FFs used to modeling nonstandard residues in Zn proteins for molecular dynamics simulations. The effectiveness of the HessFit protocol makes it a valuable tool for deriving or extending force field parameters for novel compounds in several molecular modeling applications.
Collapse
Affiliation(s)
- Emanuele Falbo
- Department of Pharmacy, Drug Discovery Lab, University of Naples Federico II, Naples 80131, Italy
| | - Antonio Lavecchia
- Department of Pharmacy, Drug Discovery Lab, University of Naples Federico II, Naples 80131, Italy
| |
Collapse
|
3
|
Adhikari S, Mondal J. Machine Learning Subtle Conformational Change due to Phosphorylation in Intrinsically Disordered Proteins. J Phys Chem B 2023; 127:9433-9449. [PMID: 37905972 DOI: 10.1021/acs.jpcb.3c05136] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2023]
Abstract
Phosphorylation of intrinsically disordered proteins/regions (IDPs/IDRs) has a profound effect in biological functions such as cell signaling, protein folding or unfolding, and long-range allosteric effects. However, here we focus on two IDPs, namely 83-residue IDR transcription factor Ash1 and 92-residue long N-terminal region of CDK inhibitor Sic1 protein, found in Saccharomyces cerevisiae, for which experimental measurements of average conformational properties, namely, radius of gyration and structure factor, indicate negligible changes upon phosphorylation. Here, we show that a judicious dissection of conformational ensemble via combination of unsupervised machine learning and extensive molecular dynamics (MD) trajectories can highlight key differences and similarities among the phosphorylated and wild-type IDP. In particular, we develop Markov state model (MSM) using the latent-space dimensions of an autoencoder, trained using multi-microsecond long MD simulation trajectories. Examination of structural changes among the states, prior to and upon phosphorylation, captured several similarities and differences in their backbone contact maps, secondary structure, and torsion angles. Hydrogen bonding analysis revealed that phosphorylation not only increases the number of hydrogen bonds but also switches the pattern of hydrogen bonding between the backbone and side chain atoms with the phosphorylated residues. We also observe that although phosphorylation introduces salt bridges, there is a loss of the cation-π interaction. Phosphorylation also improved the probability for long-range hydrophobic contacts and also enhanced interaction with water molecules and improved the local structure of water as evident from the geometric order parameters. The observations on these machine-learnt states gave important insights, as it would otherwise be difficult to determine experimentally which is important, if we were to understand the role of phosphorylation of IDPs in their biological functions.
Collapse
|
4
|
Wych DC, Wall ME. Molecular-dynamics simulations of macromolecular diffraction, part I: Preparation of protein crystal simulations. Methods Enzymol 2023; 688:87-114. [PMID: 37748833 DOI: 10.1016/bs.mie.2023.06.008] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/27/2023]
Abstract
Molecular-dynamics (MD) simulations of protein crystals enable the prediction of structural and dynamical features of both the protein and the solvent components of macromolecular crystals, which can be validated against diffraction data from X-ray crystallographic experiments. The simulations have been useful for studying and predicting both Bragg and diffuse scattering in protein crystallography; however, the preparation is not yet automated and includes choices and tradeoffs that can impact the results. Here we examine some of the intricacies and consequences of the choices involved in setting up MD simulations of protein crystals for the study of diffraction data, and provide a recipe for preparing the simulations, packaged in an accompanying Jupyter notebook. This article and the accompanying notebook are intended to serve as practical resources for researchers wishing to put these models to work.
Collapse
Affiliation(s)
- David C Wych
- Computer, Computational and Statistical Sciences Division, Los Alamos, NM, United States; Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, NM, United States
| | - Michael E Wall
- Computer, Computational and Statistical Sciences Division, Los Alamos, NM, United States.
| |
Collapse
|
5
|
Song G, Zhong B, Zhang B, Rehman AU, Chen HF. Phosphorylation Modification Force Field FB18CMAP Improving Conformation Sampling of Phosphoproteins. J Chem Inf Model 2023; 63:1602-1614. [PMID: 36800279 DOI: 10.1021/acs.jcim.3c00112] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/18/2023]
Abstract
Phosphorylation of proteins plays an important regulatory role at almost all levels of cellular organization. Molecular dynamics (MD) simulation is a promising tool to reveal the mechanism of how phosphorylation regulates many key biological processes at the atomistic level. MD simulation accuracy depends on force field precision, while the current force fields for phospho-amino acids have resulted in notable inconsistency with experimental data. Here, a new force field parameter (named FB18CMAP) is generated by fitting against quantum mechanics (QM) energy in aqueous solution with φ/ψ dihedral potential-energy surfaces optimized using CMAP parameters. MD simulations of phosphorylated dipeptides, intrinsically disordered proteins (IDPs), and ordered (folded) proteins show that FB18CMAP can mimic NMR observables and structural characteristics of phosphorylated dipeptides and proteins more accurately than the FB18 force field. These findings suggest that FB18CMAP performs well in both the simulation of ordered and disordered states of phosphorylated proteins.
Collapse
Affiliation(s)
- Ge Song
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic and Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Bozitao Zhong
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic and Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Bo Zhang
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic and Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Ashfaq Ur Rehman
- Departments of Molecular Biology and Biochemistry, University of California, Irvine, California 92697, United States
| | - Hai-Feng Chen
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic and Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China.,Shanghai Center for Bioinformation Technology, Shanghai 200240, China
| |
Collapse
|
6
|
Weigle AT, Feng J, Shukla D. Thirty years of molecular dynamics simulations on posttranslational modifications of proteins. Phys Chem Chem Phys 2022; 24:26371-26397. [PMID: 36285789 PMCID: PMC9704509 DOI: 10.1039/d2cp02883b] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/06/2023]
Abstract
Posttranslational modifications (PTMs) are an integral component to how cells respond to perturbation. While experimental advances have enabled improved PTM identification capabilities, the same throughput for characterizing how structural changes caused by PTMs equate to altered physiological function has not been maintained. In this Perspective, we cover the history of computational modeling and molecular dynamics simulations which have characterized the structural implications of PTMs. We distinguish results from different molecular dynamics studies based upon the timescales simulated and analysis approaches used for PTM characterization. Lastly, we offer insights into how opportunities for modern research efforts on in silico PTM characterization may proceed given current state-of-the-art computing capabilities and methodological advancements.
Collapse
Affiliation(s)
- Austin T Weigle
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Jiangyan Feng
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Diwakar Shukla
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
- Department of Bioengineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
- Department of Plant Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA.
| |
Collapse
|
7
|
Stoppelman JP, Ng TT, Nerenberg PS, Wang LP. Correction to Development and Validation of AMBER-FB15-Compatible Force Field Parameters for Phosphorylated Amino Acids. J Phys Chem B 2022; 126:8596. [PMID: 36223551 DOI: 10.1021/acs.jpcb.2c06820] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
8
|
Balanced Force Field ff03CMAP Improving the Dynamics Conformation Sampling of Phosphorylation Site. Int J Mol Sci 2022; 23:ijms231911285. [PMID: 36232586 PMCID: PMC9569523 DOI: 10.3390/ijms231911285] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2022] [Revised: 09/08/2022] [Accepted: 09/11/2022] [Indexed: 11/30/2022] Open
Abstract
Phosphorylation plays a key role in plant biology, such as the accumulation of plant cells to form the observed proteome. Statistical analysis found that many phosphorylation sites are located in disordered regions. However, current force fields are mainly trained for structural proteins, which might not have the capacity to perfectly capture the dynamic conformation of the phosphorylated proteins. Therefore, we evaluated the performance of ff03CMAP, a balanced force field between structural and disordered proteins, for the sampling of the phosphorylated proteins. The test results of 11 different phosphorylated systems, including dipeptides, disordered proteins, folded proteins, and their complex, indicate that the ff03CMAP force field can better sample the conformations of phosphorylation sites for disordered proteins and disordered regions than ff03. For the solvent model, the results strongly suggest that the ff03CMAP force field with the TIP4PD water model is the best combination for the conformer sampling. Additional tests of CHARMM36m and FB18 force fields on two phosphorylated systems suggest that the overall performance of ff03CMAP is similar to that of FB18 and better than that of CHARMM36m. These results can help other researchers to choose suitable force field and solvent models to investigate the dynamic properties of phosphorylation proteins.
Collapse
|
9
|
Yang D, Gronenborn AM, Chong LT. Development and Validation of Fluorinated, Aromatic Amino Acid Parameters for Use with the AMBER ff15ipq Protein Force Field. J Phys Chem A 2022; 126:2286-2297. [PMID: 35352936 PMCID: PMC9014858 DOI: 10.1021/acs.jpca.2c00255] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2022] [Revised: 03/15/2022] [Indexed: 12/27/2022]
Abstract
We developed force field parameters for fluorinated, aromatic amino acids enabling molecular dynamics (MD) simulations of fluorinated proteins. These parameters are tailored to the AMBER ff15ipq protein force field and enable the modeling of 4, 5, 6, and 7F-tryptophan, 3F- and 3,5F-tyrosine, and 4F- or 4-CF3-phenylalanine. The parameters include 181 unique atomic charges derived using the implicitly polarized charge (IPolQ) scheme in the presence of SPC/Eb explicit water molecules and 9 unique bond, angle, or torsion terms. Our simulations of benchmark peptides and proteins maintain expected conformational propensities on the μs time scale. In addition, we have developed an open-source Python program to calculate fluorine relaxation rates from MD simulations. The extracted relaxation rates from protein simulations are in good agreement with experimental values determined by 19F NMR. Collectively, our results illustrate the power and robustness of the IPolQ lineage of force fields for modeling the structure and dynamics of fluorine-containing proteins at the atomic level.
Collapse
Affiliation(s)
- Darian
T. Yang
- Molecular
Biophysics and Structural Biology Graduate Program, University of Pittsburgh and Carnegie Mellon University, Pittsburgh, Pennsylvania 15260, United States
- Department
of Structural Biology, University of Pittsburgh
School of Medicine, Pittsburgh, Pennsylvania 15260, United States
- Department
of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Angela M. Gronenborn
- Department
of Structural Biology, University of Pittsburgh
School of Medicine, Pittsburgh, Pennsylvania 15260, United States
- Department
of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Lillian T. Chong
- Department
of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| |
Collapse
|