1
|
Varvdekar B, Prabhakant A, Krishnan M. Response of Terahertz Protein Vibrations to Ligand Binding: Calmodulin-Peptide Complexes as a Case Study. J Chem Inf Model 2022; 62:1669-1679. [PMID: 35312312 DOI: 10.1021/acs.jcim.1c01344] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Terahertz vibrations are sensitive reporters of the structure and interactions of proteins. Ligand binding alters the nature and distribution of these collective vibrations. The ligand-induced changes in the terahertz protein vibrations contribute to the binding entropy and to the overall thermodynamic stability of the resultant protein-ligand complexes. Here, we have examined the response of the low-frequency (below 6 terahertz) collective vibrations of the calcium-loaded calmodulin (CaM) to binding to five different ligands, both in the presence and absence of water, using normal-mode analysis and molecular dynamics simulations. A comparison of the vibrational spectra of hydrated and dry systems reveals that protein-solvent interactions stiffen the terahertz protein vibrations and that these solvent-coupled collective vibrations contribute significantly to the hydration-sensitive variation in the vibrational entropy of CaM. In the absence of water, the low-frequency vibrations of CaM are stiffened by ligand binding. On the contrary, the number and the cumulative vibrational entropy of low-frequency vibrational modes (ω < 200 cm-1) of the hydrated CaM are increased noticeably after binding to the peptides, indicating binding-induced softening of collective vibrations of the protein. Although the calculated and experimental binding affinities of the chosen complexes correlated reasonably well, no systematic correlation was observed between the protein vibrational entropy and the binding affinity. The results underscored the importance of the interplay of protein-ligand and solvent interactions in modulating the low-frequency vibrations of proteins.
Collapse
Affiliation(s)
- Bhagyesh Varvdekar
- Center for Computational Natural Sciences and Bioinformatics (CCNSB), International Institute of Information Technology, Gachibowli, Hyderabad 500032, India
| | - Akshay Prabhakant
- Center for Computational Natural Sciences and Bioinformatics (CCNSB), International Institute of Information Technology, Gachibowli, Hyderabad 500032, India
| | - Marimuthu Krishnan
- Center for Computational Natural Sciences and Bioinformatics (CCNSB), International Institute of Information Technology, Gachibowli, Hyderabad 500032, India
| |
Collapse
|
2
|
Gumbart JC, Ulmschneider MB, Hazel A, White SH, Ulmschneider JP. Computed Free Energies of Peptide Insertion into Bilayers are Independent of Computational Method. J Membr Biol 2018; 251:345-356. [PMID: 29520628 PMCID: PMC6030508 DOI: 10.1007/s00232-018-0026-y] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2017] [Accepted: 02/27/2018] [Indexed: 01/15/2023]
Abstract
We show that the free energy of inserting hydrophobic peptides into lipid bilayer membranes from surface-aligned to transmembrane inserted states can be reliably calculated using atomistic models. We use two entirely different computational methods: high temperature spontaneous peptide insertion calculations as well as umbrella sampling potential-of-mean-force (PMF) calculations, both yielding the same energetic profiles. The insertion free energies were calculated using two different protein and lipid force fields (OPLS protein/united-atom lipids and CHARMM36 protein/all-atom lipids) and found to be independent of the simulation parameters. In addition, the free energy of insertion is found to be independent of temperature for both force fields. However, we find major difference in the partitioning kinetics between OPLS and CHARMM36, likely due to the difference in roughness of the underlying free energy surfaces. Our results demonstrate not only a reliable method to calculate insertion free energies for peptides, but also represent a rare case where equilibrium simulations and PMF calculations can be directly compared.
Collapse
Affiliation(s)
| | | | | | - Stephen H White
- Department of Physiology & Biophysics, University of California at Irvine, Irvine, CA, USA
| | - Jakob P Ulmschneider
- Department of Physics and the Institute of Natural Sciences, Shanghai Jiao Tong University, Shanghai, 200240, China.
| |
Collapse
|
3
|
Romero E, Ladani ST, Hamelberg D, Gadda G. Solvent-Slaved Motions in the Hydride Tunneling Reaction Catalyzed by Human Glycolate Oxidase. ACS Catal 2016. [DOI: 10.1021/acscatal.5b02889] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Elvira Romero
- Department of Chemistry, ¶Department of Biology, ∥Center for Biotechnology
and Drug
Design, and #Center
for Diagnostics and Therapeutics, Georgia State University, Atlanta, Georgia 30302-3965, United States
| | - Safieh Tork Ladani
- Department of Chemistry, ¶Department of Biology, ∥Center for Biotechnology
and Drug
Design, and #Center
for Diagnostics and Therapeutics, Georgia State University, Atlanta, Georgia 30302-3965, United States
| | - Donald Hamelberg
- Department of Chemistry, ¶Department of Biology, ∥Center for Biotechnology
and Drug
Design, and #Center
for Diagnostics and Therapeutics, Georgia State University, Atlanta, Georgia 30302-3965, United States
| | - Giovanni Gadda
- Department of Chemistry, ¶Department of Biology, ∥Center for Biotechnology
and Drug
Design, and #Center
for Diagnostics and Therapeutics, Georgia State University, Atlanta, Georgia 30302-3965, United States
| |
Collapse
|
4
|
Anandakrishnan R, Drozdetski A, Walker RC, Onufriev AV. Speed of conformational change: comparing explicit and implicit solvent molecular dynamics simulations. Biophys J 2016; 108:1153-64. [PMID: 25762327 DOI: 10.1016/j.bpj.2014.12.047] [Citation(s) in RCA: 130] [Impact Index Per Article: 16.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2014] [Revised: 12/23/2014] [Accepted: 12/29/2014] [Indexed: 11/24/2022] Open
Abstract
Adequate sampling of conformation space remains challenging in atomistic simulations, especially if the solvent is treated explicitly. Implicit-solvent simulations can speed up conformational sampling significantly. We compare the speed of conformational sampling between two commonly used methods of each class: the explicit-solvent particle mesh Ewald (PME) with TIP3P water model and a popular generalized Born (GB) implicit-solvent model, as implemented in the AMBER package. We systematically investigate small (dihedral angle flips in a protein), large (nucleosome tail collapse and DNA unwrapping), and mixed (folding of a miniprotein) conformational changes, with nominal simulation times ranging from nanoseconds to microseconds depending on system size. The speedups in conformational sampling for GB relative to PME simulations, are highly system- and problem-dependent. Where the simulation temperatures for PME and GB are the same, the corresponding speedups are approximately onefold (small conformational changes), between ∼1- and ∼100-fold (large changes), and approximately sevenfold (mixed case). The effects of temperature on speedup and free-energy landscapes, which may differ substantially between the solvent models, are discussed in detail for the case of miniprotein folding. In addition to speeding up conformational sampling, due to algorithmic differences, the implicit solvent model can be computationally faster for small systems or slower for large systems, depending on the number of solute and solvent atoms. For the conformational changes considered here, the combined speedups are approximately twofold, ∼1- to 60-fold, and ∼50-fold, respectively, in the low solvent viscosity regime afforded by the implicit solvent. For all the systems studied, 1) conformational sampling speedup increases as Langevin collision frequency (effective viscosity) decreases; and 2) conformational sampling speedup is mainly due to reduction in solvent viscosity rather than possible differences in free-energy landscapes between the solvent models.
Collapse
Affiliation(s)
| | | | - Ross C Walker
- San Diego Supercomputer Center and Department of Chemistry and Biochemistry, University of California, San Diego, San Diego, California
| | - Alexey V Onufriev
- Department of Computer Science, Virginia Tech, Blacksburg, Virginia; Department of Physics, Virginia Tech, Blacksburg, Virginia.
| |
Collapse
|
5
|
Barman A, Smitherman C, Souffrant M, Gadda G, Hamelberg D. Conserved Hydration Sites in Pin1 Reveal a Distinctive Water Recognition Motif in Proteins. J Chem Inf Model 2015; 56:139-47. [DOI: 10.1021/acs.jcim.5b00560] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Arghya Barman
- Departments
of Chemistry and ‡Biology and the §Centers for Diagnostics and Therapeutics and ∥Biotechnology
and Drug Design, Georgia State University, Atlanta, Georgia 30302-3965, United States
| | - Crystal Smitherman
- Departments
of Chemistry and ‡Biology and the §Centers for Diagnostics and Therapeutics and ∥Biotechnology
and Drug Design, Georgia State University, Atlanta, Georgia 30302-3965, United States
| | - Michael Souffrant
- Departments
of Chemistry and ‡Biology and the §Centers for Diagnostics and Therapeutics and ∥Biotechnology
and Drug Design, Georgia State University, Atlanta, Georgia 30302-3965, United States
| | - Giovanni Gadda
- Departments
of Chemistry and ‡Biology and the §Centers for Diagnostics and Therapeutics and ∥Biotechnology
and Drug Design, Georgia State University, Atlanta, Georgia 30302-3965, United States
| | - Donald Hamelberg
- Departments
of Chemistry and ‡Biology and the §Centers for Diagnostics and Therapeutics and ∥Biotechnology
and Drug Design, Georgia State University, Atlanta, Georgia 30302-3965, United States
| |
Collapse
|
6
|
Moritsugu K, Kidera A, Smith JC. Solvent Friction Effects Propagate over the Entire Protein Molecule through Low-Frequency Collective Modes. J Phys Chem B 2014; 118:8559-65. [DOI: 10.1021/jp503956m] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Kei Moritsugu
- Graduate
School of Medical Life Science, Yokohama City University, 1-7-29
Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Akinori Kidera
- Graduate
School of Medical Life Science, Yokohama City University, 1-7-29
Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Jeremy C. Smith
- Center
for Molecular Biophysics, University of Tennessee/Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| |
Collapse
|
7
|
Hayes RL, Noel JK, Whitford PC, Mohanty U, Sanbonmatsu KY, Onuchic JN. Reduced model captures Mg(2+)-RNA interaction free energy of riboswitches. Biophys J 2014; 106:1508-19. [PMID: 24703312 PMCID: PMC3976530 DOI: 10.1016/j.bpj.2014.01.042] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2013] [Revised: 01/22/2014] [Accepted: 01/29/2014] [Indexed: 12/17/2022] Open
Abstract
The stability of RNA tertiary structures depends heavily on Mg(2+). The Mg(2+)-RNA interaction free energy that stabilizes an RNA structure can be computed experimentally through fluorescence-based assays that measure Γ2+, the number of excess Mg(2+) associated with an RNA molecule. Previous explicit-solvent simulations predict that the majority of excess Mg(2+) ions interact closely and strongly with the RNA, unlike monovalent ions such as K(+), suggesting that an explicit treatment of Mg(2+) is important for capturing RNA dynamics. Here we present a reduced model that accurately reproduces the thermodynamics of Mg(2+)-RNA interactions. This model is able to characterize long-timescale RNA dynamics coupled to Mg(2+) through the explicit representation of Mg(2+) ions. KCl is described by Debye-Hückel screening and a Manning condensation parameter, which represents condensed K(+) and models its competition with condensed Mg(2+). The model contains one fitted parameter, the number of condensed K(+) ions in the absence of Mg(2+). Values of Γ2+ computed from molecular dynamics simulations using the model show excellent agreement with both experimental data on the adenine riboswitch and previous explicit-solvent simulations of the SAM-I riboswitch. This agreement confirms the thermodynamic accuracy of the model via the direct relation of Γ2+ to the Mg(2+)-RNA interaction free energy, and provides further support for the predictions from explicit-solvent calculations. This reduced model will be useful for future studies of the interplay between Mg(2+) and RNA dynamics.
Collapse
Affiliation(s)
- Ryan L Hayes
- Center for Theoretical Biological Physics and Department of Physics and Astronomy, Rice University, Houston, Texas
| | - Jeffrey K Noel
- Center for Theoretical Biological Physics and Department of Physics and Astronomy, Rice University, Houston, Texas
| | - Paul C Whitford
- Department of Physics, Northeastern University, Boston, Massachusetts
| | - Udayan Mohanty
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts
| | - Karissa Y Sanbonmatsu
- Theoretical Biology and Biophysics, Theoretical Division, Los Alamos National Labs, Los Alamos, New Mexico.
| | - José N Onuchic
- Center for Theoretical Biological Physics and Department of Physics and Astronomy, Rice University, Houston, Texas; Department of Chemistry and Department of Biochemistry and Cell Biology, Rice University, Houston, Texas.
| |
Collapse
|
8
|
Normal-mode-analysis-guided investigation of crucial intersubunit contacts in the cAMP-dependent gating in HCN channels. Biophys J 2012; 103:19-28. [PMID: 22828328 DOI: 10.1016/j.bpj.2012.05.030] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2011] [Revised: 05/21/2012] [Accepted: 05/22/2012] [Indexed: 11/23/2022] Open
Abstract
Protein structures define a complex network of atomic interactions in three dimensions. Direct visualization of the structure and analysis of the interaction potential energy are not straightforward approaches to pinpoint the atomic contacts that are crucial for protein function. We used the tetrameric hyperpolarization-activated cAMP-regulated (HCN) channel as a model system to study the intersubunit contacts in cAMP-dependent gating. To obtain a systematic survey of the contacts between each pair of residues, we used normal-mode analysis, a computational approach for studying protein dynamics, and constructed the covariance matrix for C-α atoms. The significant contacts revealed by covariance analysis were further investigated by means of mutagenesis and functional assays. Among the mutant channels that show phenotypes different from those of the wild-type, we focused on two mutant channels that express opposite changes in cAMP-dependent gating. Subsequent biochemical assays on isolated C-terminal fragments, including the cAMP binding domain, revealed only minimal effects on cAMP binding, suggesting the necessity of interpreting the cAMP-dependent allosteric regulation at the whole-channel level. For this purpose, we applied the patch-clamp fluorometry technique and observed correlated changes in the dynamic, state-dependent cAMP binding in the mutant channels. This study not only provides further understanding of the intersubunit contacts in allosteric coupling in the HCN channel, it also illustrates an effective strategy for delineating important atomic contacts within a structure.
Collapse
|
9
|
Ladani ST, Hamelberg D. Entropic and Surprisingly Small Intramolecular Polarization Effects in the Mechanism of Cyclophilin A. J Phys Chem B 2012; 116:10771-8. [DOI: 10.1021/jp305917c] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Safieh Tork Ladani
- Department of Chemistry and Center
for Biotechnology
and Drug Design, Georgia State University, Atlanta, Georgia 30302-4098, United States
| | - Donald Hamelberg
- Department of Chemistry and Center
for Biotechnology
and Drug Design, Georgia State University, Atlanta, Georgia 30302-4098, United States
| |
Collapse
|
10
|
Nellas RB, Glover MM, Hamelberg D, Shen T. High-pressure effect on the dynamics of solvated peptides. J Chem Phys 2012; 136:145103. [DOI: 10.1063/1.3700183] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
11
|
Voltz K, Trylska J, Calimet N, Smith JC, Langowski J. Unwrapping of nucleosomal DNA ends: a multiscale molecular dynamics study. Biophys J 2012; 102:849-58. [PMID: 22385856 DOI: 10.1016/j.bpj.2011.11.4028] [Citation(s) in RCA: 59] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2011] [Revised: 11/10/2011] [Accepted: 11/28/2011] [Indexed: 11/25/2022] Open
Abstract
To permit access to DNA-binding proteins involved in the control and expression of the genome, the nucleosome undergoes structural remodeling including unwrapping of nucleosomal DNA segments from the nucleosome core. Here we examine the mechanism of DNA dissociation from the nucleosome using microsecond timescale coarse-grained molecular dynamics simulations. The simulations exhibit short-lived, reversible DNA detachments from the nucleosome and long-lived DNA detachments not reversible on the timescale of the simulation. During the short-lived DNA detachments, 9 bp dissociate at one extremity of the nucleosome core and the H3 tail occupies the space freed by the detached DNA. The long-lived DNA detachments are characterized by structural rearrangements of the H3 tail including the formation of a turn-like structure at the base of the tail that sterically impedes the rewrapping of DNA on the nucleosome surface. Removal of the H3 tails causes the long-lived detachments to disappear. The physical consistency of the CG long-lived open state was verified by mapping a CG structure representative of this state back to atomic resolution and performing molecular dynamics as well as by comparing conformation-dependent free energies. Our results suggest that the H3 tail may stabilize the nucleosome in the open state during the initial stages of the nucleosome remodeling process.
Collapse
Affiliation(s)
- Karine Voltz
- Biophysics of Macromolecules, German Cancer Research Center, Heidelberg, Germany
| | | | | | | | | |
Collapse
|
12
|
Xue Y, Skrynnikov NR. Motion of a disordered polypeptide chain as studied by paramagnetic relaxation enhancements, 15N relaxation, and molecular dynamics simulations: how fast is segmental diffusion in denatured ubiquitin? J Am Chem Soc 2011; 133:14614-28. [PMID: 21819149 DOI: 10.1021/ja201605c] [Citation(s) in RCA: 48] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Molecular dynamics (MD) simulations have been widely used to analyze dynamic conformational equilibria of folded proteins, especially in relation to NMR observables. However, this approach found little use in the studies of disordered proteins, where the sampling of vast conformational space presents a serious problem. In this paper, we demonstrate that the latest advances in computation technology make it possible to overcome this limitation. The experimentally validated (calibrated) MD models allow for new insights into structure/dynamics of disordered proteins. As a test system, we have chosen denatured ubiquitin in solution with 8 M urea at pH 2. High-temperature MD simulations in implicit solvent have been carried out for the wild-type ubiquitin as well as MTSL-tagged Q2C, D32C, and R74C mutants. To recalibrate the MD data (500 K) in relation to the experimental conditions (278 K, 8 M urea), the time axes of the MD trajectories were rescaled. The scaling factor was adjusted such as to maximize the agreement between the simulated and experimental (15)N relaxation rates. The resulting effective length of the trajectories, 311 μs, ensures good convergence properties of the MD model. The constructed MD model was validated against the array of experimental data, including additional (15)N relaxation parameters, multiple sets of paramagnetic relaxation enhancements (PREs), and the radius of gyration. In each case, a near-quantitative agreement has been obtained, suggesting that the model is successful. Of note, the MD-based approach rigorously predicts the quantities that are inherently dynamic, i.e., dependent on the motional correlation times. This cannot be accomplished, other than in empirical fashion, on the basis of static structural models (conformational ensembles). The MD model was further used to investigate the relative translational motion of the MTSL label and the individual H(N) atoms. The derived segmental diffusion coefficients proved to be nearly uniform along the peptide chain, averaging to D = 0.49-0.55 × 10(-6) cm(2)/s. This result was verified by direct analysis of the experimental PRE data using the recently proposed Ullman-Podkorytov model. In this model, MTSL and H(N) moieties are treated as two tethered spheres undergoing mutual diffusion in a harmonic potential. The fitting of the experimental data involving D as a single adjustable parameter leads to D = 0.45 × 10(-6) cm(2)/s, in good agreement with the MD-based analyses. This result can be compared with the range of estimates obtained from the resonance energy transfer experiments, D = 0.2-6.0 × 10(-6) cm(2)/s.
Collapse
Affiliation(s)
- Yi Xue
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907-2084, USA
| | | |
Collapse
|
13
|
Wang Y, Markwick PRL, de Oliveira CAF, McCammon JA. Enhanced Lipid Diffusion and Mixing in Accelerated Molecular Dynamics. J Chem Theory Comput 2011; 7:3199-3207. [PMID: 22003320 PMCID: PMC3191728 DOI: 10.1021/ct200430c] [Citation(s) in RCA: 56] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2011] [Indexed: 11/28/2022]
Abstract
Accelerated molecular dynamics (aMD) is an enhanced sampling technique that expedites conformational space sampling by reducing the barriers separating various low-energy states of a system. Here, we present the first application of the aMD method on lipid membranes. Altogether, ∼1.5 μs simulations were performed on three systems: a pure POPC bilayer, a pure DMPC bilayer, and a mixed POPC:DMPC bilayer. Overall, the aMD simulations are found to produce significant speedup in trans–gauche isomerization and lipid lateral diffusion versus those in conventional MD (cMD) simulations. Further comparison of a 70-ns aMD run and a 300-ns cMD run of the mixed POPC:DMPC bilayer shows that the two simulations yield similar lipid mixing behaviors, with aMD generating a 2–3-fold speedup compared to cMD. Our results demonstrate that the aMD method is an efficient approach for the study of bilayer structural and dynamic properties. On the basis of simulations of the three bilayer systems, we also discuss the impact of aMD parameters on various lipid properties, which can be used as a guideline for future aMD simulations of membrane systems.
Collapse
Affiliation(s)
- Yi Wang
- Center for Theoretical Biological Physics, Howard Hughes Medical Institute, Department of Chemistry and Biochemistry, Department of Pharmacology, University of California, San Diego , La Jolla, California 92093, United States
| | | | | | | |
Collapse
|
14
|
Role of hydration in collagen recognition by bacterial adhesins. Biophys J 2011; 100:2253-61. [PMID: 21539794 DOI: 10.1016/j.bpj.2011.03.033] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2010] [Revised: 02/22/2011] [Accepted: 03/09/2011] [Indexed: 11/23/2022] Open
Abstract
Protein-protein recognition regulates the vast majority of physiological or pathological processes. We investigated the role of hydration in collagen recognition by bacterial adhesin CNA by means of first principle molecular-dynamics samplings. Our characterization of the hydration properties of the isolated partners highlights dewetting-prone areas on the surface of CNA that closely match the key regions involved in hydrophobic intermolecular interactions upon complex formation, suggesting that the hydration state of the ligand-free CNA predisposes the protein to the collagen recognition. Moreover, hydration maps of the CNA-collagen complex reveal the presence of a number of structured water molecules that mediate intermolecular interactions at the interface between the two proteins. These hydration sites feature long residence times, significant binding free energies, and a geometrical distribution that closely resembles the hydration pattern of the isolated collagen triple helix. These findings are striking evidence that CNA recognizes the collagen triple helix as a hydrated molecule. For this structural motif, the exposure of several unsatisfied backbone carbonyl groups results in a strong interplay with the solvent, which is shown to also play a role in collagen recognition.
Collapse
|
15
|
Peters MH. Langevin dynamics for the transport of flexible biological macromolecules in confined geometries. J Chem Phys 2011; 134:025105. [DOI: 10.1063/1.3525381] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
16
|
Markwick PRL, McCammon JA. Studying functional dynamics in bio-molecules using accelerated molecular dynamics. Phys Chem Chem Phys 2011; 13:20053-65. [DOI: 10.1039/c1cp22100k] [Citation(s) in RCA: 139] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
17
|
Xin Y, Hamelberg D. Deciphering the role of glucosamine-6-phosphate in the riboswitch action of glmS ribozyme. RNA (NEW YORK, N.Y.) 2010; 16:2455-63. [PMID: 20971809 PMCID: PMC2995406 DOI: 10.1261/rna.2334110] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/22/2010] [Accepted: 09/16/2010] [Indexed: 05/23/2023]
Abstract
The GlmS ribozyme is believed to exploit a general acid-base catalytic mechanism in the presence of glucosamine-6-phosphate (GlcN6P) to accelerate self-cleavage by approximately six orders of magnitude. The general acid and general base are not known, and the role of the GlcN6P cofactor is even less well understood. The amine group of GlcN6P has the ability to either accept or donate a proton and could therefore potentially act as an acid or a base. In order to decipher the role of GlcN6P in the self-cleavage of glmS, we have determined the preferred protonation state of the amine group in the wild-type and an inactive G40A mutant using molecular dynamics simulations and free energy calculations. Here we show that, upon binding of GlcN6P to wild-type glmS, the pK(a) of the amine moiety is altered by the active site environment, decreasing by about 2.2 from a solution pK(a) of about 8.2. On the other hand, we show that the pK(a) of the amine group slightly increases to about 8.4 upon binding to the G40A inactive mutant of glmS. These results suggest that GlcN6P acts as a general acid in the self-cleavage of glmS. Upon binding to glmS, GlcN6P can easily release a proton to the 5'-oxygen of G1 during self-cleavage of the backbone phosphodiester bond. However, in the G40A inactive mutant of glmS, the results suggest that the ability of GlcN6P to easily release its proton is diminished, in addition to the possible lack of G40 as an effective base.
Collapse
Affiliation(s)
- Yao Xin
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302-4098, USA
| | | |
Collapse
|
18
|
Arora K, Brooks Iii CL. Functionally important conformations of the Met20 loop in dihydrofolate reductase are populated by rapid thermal fluctuations. J Am Chem Soc 2010; 131:5642-7. [PMID: 19323547 DOI: 10.1021/ja9000135] [Citation(s) in RCA: 47] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Conformational changes in enzymes are well recognized to play an important role in the organization of the reactive groups for efficient catalysis. This study reveals atomic and energetic details of the conformational change process that precedes the catalytic reaction of the enzyme dihydrofolate reductase. The computed free energy profile provides insights into the ligand binding mechanism and a quantitative estimate of barrier heights separating different conformational states along the pathway. Studies show that the ternary complex comprised of NADPH cofactor and substrate dihydrofolate undergoes transitions between a closed state and an occluded state via an intermediate "open" conformation. During these transitions the largest conformational change occurs in the Met20 loop of DHFR and is accompanied by the motion of the cofactor into and out of the binding pocket. When the cofactor is out of the binding pocket, the enzyme frequently samples open and occluded conformations with a small (approximately 5 k(B)T) free energy barrier between the two states. However, when the cofactor is in the binding pocket, the closed conformation is thermodynamically most favored. The determination of a profile characterizing the position-dependent diffusion of the Met20 loop allowed us to apply reaction rate theory and deduce the kinetics of loop motions based on the computed free energy landscape.
Collapse
Affiliation(s)
- Karunesh Arora
- Department of Chemistry and Biophysics Program, University of Michigan, 930 North University Avenue, Ann Arbor, Michigan 48109, USA
| | | |
Collapse
|
19
|
Xin Y, Doshi U, Hamelberg D. Examining the limits of time reweighting and Kramers' rate theory to obtain correct kinetics from accelerated molecular dynamics. J Chem Phys 2010; 132:224101. [PMID: 20550384 DOI: 10.1063/1.3432761] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
Accelerated molecular dynamics simulations are routinely being used to recover the correct canonical probability distributions corresponding to the original potential energy landscape of biomolecular systems. However, the limits of time reweighting, based on transition state theory, in obtaining true kinetic rates from accelerated molecular dynamics for biomolecular systems are less obvious. Here, we investigate this issue by studying the kinetics of cis-trans isomerization of peptidic omega bond by accelerated molecular dynamics. We find that time reweighting is valid for obtaining true kinetics when the original potential is not altered at the transition state regions, as expected. When the original potential landscape is modified such that the applied boost potential alters the transition state regions, time reweighting fails to reproduce correct kinetics and the reweighted rate is much slower than the true rate. By adopting the overdamped limit of Kramers' rate theory, we are successful in recovering correct kinetics irrespective of whether or not the transition state regions are modified. Furthermore, we tested the validity of the acceleration weight factor from the path integral formalism for obtaining the correct kinetics of cis-trans isomerization. It was found that this formulation of the weight factor is not suitable for long time scale processes such as cis-trans isomerization with high energy barriers.
Collapse
Affiliation(s)
- Yao Xin
- Department of Chemistry and The Center for Biotechnology and Drug Design, Georgia State University, Atlanta, Georgia 30302-4098, USA
| | | | | |
Collapse
|
20
|
Johnson Q, Doshi U, Shen T, Hamelberg D. Water’s Contribution to the Energetic Roughness from Peptide Dynamics. J Chem Theory Comput 2010; 6:2591-7. [DOI: 10.1021/ct100183s] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Quentin Johnson
- Department of Chemistry and the Center for Biotechnology and Drug Design, Georgia State University, Atlanta, Georgia 30302-4098, Department of Biochemistry, Cellular & Molecular Biology, University of Tennessee, Knoxville, Tennessee 37996, and Center for Molecular Biophysics, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37830
| | - Urmi Doshi
- Department of Chemistry and the Center for Biotechnology and Drug Design, Georgia State University, Atlanta, Georgia 30302-4098, Department of Biochemistry, Cellular & Molecular Biology, University of Tennessee, Knoxville, Tennessee 37996, and Center for Molecular Biophysics, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37830
| | - Tongye Shen
- Department of Chemistry and the Center for Biotechnology and Drug Design, Georgia State University, Atlanta, Georgia 30302-4098, Department of Biochemistry, Cellular & Molecular Biology, University of Tennessee, Knoxville, Tennessee 37996, and Center for Molecular Biophysics, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37830
| | - Donald Hamelberg
- Department of Chemistry and the Center for Biotechnology and Drug Design, Georgia State University, Atlanta, Georgia 30302-4098, Department of Biochemistry, Cellular & Molecular Biology, University of Tennessee, Knoxville, Tennessee 37996, and Center for Molecular Biophysics, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37830
| |
Collapse
|
21
|
Ding Y, Mamonov AB, Zuckerman DM. Efficient equilibrium sampling of all-atom peptides using library-based Monte Carlo. J Phys Chem B 2010; 114:5870-7. [PMID: 20380366 PMCID: PMC2882875 DOI: 10.1021/jp910112d] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We applied our previously developed library-based Monte Carlo (LBMC) to equilibrium sampling of several implicitly solvated all-atom peptides. LBMC can perform equilibrium sampling of molecules using precalculated statistical libraries of molecular-fragment configurations and energies. For this study, we employed residue-based fragments distributed according to the Boltzmann factor of the optimized potential for liquid simulations all-atom (OPLS-AA) forcefield describing the individual fragments. Two solvent models were employed: a simple uniform dielectric and the generalized Born/surface area (GBSA) model. The efficiency of LBMC was compared to standard Langevin dynamics (LD) using three different statistical tools. The statistical analyses indicate that LBMC is more than 100 times faster than LD not only for the simple solvent model but also for GBSA.
Collapse
Affiliation(s)
- Ying Ding
- Department of Computational Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania 15213
| | - Artem B. Mamonov
- Department of Computational Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania 15213
| | - Daniel M. Zuckerman
- Department of Computational Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania 15213
| |
Collapse
|
22
|
Yeh IC, Wallqvist A. Structure and dynamics of end-to-end loop formation of the penta-peptide Cys-Ala-Gly-Gln-Trp in implicit solvents. J Phys Chem B 2009; 113:12382-90. [PMID: 19685925 DOI: 10.1021/jp904064z] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
To investigate the effects of implicit solvents on peptide structure and dynamics, we performed extensive molecular dynamics simulations on the penta-peptide Cys-Ala-Gly-Gln-Trp. Two different implicit solvent models based on the CHARMM22 all-atom force field were used. Structural properties of the peptide such as distributions of end-to-end distances and dihedral angles obtained from molecular dynamics simulations with implicit solvent models were in a good agreement with those obtained from a previous explicit solvent simulation using the same force field. Representative structures observed in explicit solvent were sampled by implicit solvent models but with different relative probabilities. However, we observed significant differences in dynamical properties in explicit and implicit solvent models when we used traditional methods for the temperature control, such as Nose-Hoover or Berendsen thermostats. The explicitly solvated peptide displayed the slowest dynamics in both end-to-end contact formation and intrinsic diffusive motion of end-to-end distances. A closer agreement between implicit and explicit solvated peptide dynamics was observed when Langevin dynamics with a friction coefficient of 10 ps(-1) was used to maintain the temperature of the systems.
Collapse
Affiliation(s)
- In-Chul Yeh
- Biotechnology High Performance Computing Software Applications Institute, Telemedicine and Advanced Technology Research Center, U.S. Army Medical Research and Materiel Command, Fort Detrick, Maryland 21702-5012, USA.
| | | |
Collapse
|
23
|
Hamelberg D, McCammon JA. Mechanistic insight into the role of transition-state stabilization in cyclophilin A. J Am Chem Soc 2009; 131:147-52. [PMID: 19128175 PMCID: PMC2651649 DOI: 10.1021/ja806146g] [Citation(s) in RCA: 47] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Peptidyl prolyl cis-trans isomerases (PPIases) are ubiquitous enzymes in biology that catalyze the cis-trans isomerization of the proline imide peptide bond in many cell signaling pathways. The local change of the isomeric state of the prolyl peptide bond acts as a switching mechanism in altering the conformation of proteins. A complete understanding of the mechanism of PPIases is still lacking, and current experimental techniques have not been able to provide a detailed atomistic picture. Here we have carried out several accelerated molecular dynamics simulations with explicit solvent, and we have provided a detailed description of cis-trans isomerization of the free and cyclophilin A-catalyzed process. We show that the catalytic mechanism of cyclophilin is due mainly to the stabilization and preferential binding of the transition state that is achieved by a favorable hydrogen bond interaction with a backbone NH group. We also show that the substrate in the transition state interacts more favorably with the enzyme than the cis isomer, which in turn interacts more favorably than the trans isomer. The stability of the enzyme-substrate complex is directly correlated with the interaction the substrate makes with a highly conserved arginine residue. Finally, we show that catalysis is achieved through the rotation of the carbonyl oxygen on the N-terminal of the prolyl peptide bond in a predominately unidirectional fashion.
Collapse
Affiliation(s)
- Donald Hamelberg
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302-4098, USA
| | | |
Collapse
|
24
|
Samsonov S, Teyra J, Pisabarro MT. A molecular dynamics approach to study the importance of solvent in protein interactions. Proteins 2009; 73:515-25. [PMID: 18452208 DOI: 10.1002/prot.22076] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Water constitutes the cellular environment for biomolecules to interact. Solvent is important for protein folding and stability, and it is also known to actively participate in many catalytic processes in the cell. However, solvent is often ignored in molecular recognition and not taken into account in protein-protein interaction studies and rational design. Previously we developed SCOWLP, a database and its web application (http://www.scowlp.org), to perform studies on the contribution of solvent to protein interface definition in all protein complexes of the PDB. We introduced the concept of wet spots, interfacial residues interacting only through one water molecule, which were shown to considerably enrich protein interface descriptions. Analysis of interfacial solvent in a nonredundant dataset of protein complexes suggested the importance of including interfacial water molecules in protein interaction studies. In this work we use a molecular dynamics approach to gain deeper insights into solvent contribution to protein interfaces. We characterize the dynamic and energetic properties of water-mediated protein interactions by comparing different interfacial interaction types (direct, dual and wet spot) at residue and solvent level. For this purpose, we perform an analysis of 17 representative complexes from two protein families of different interface nature. Energetically wet spots are quantitatively comparable to other residues in interfaces, and their mobility is shown to be lower than protein surface residues. The residence time of water molecules in wet spots sites is higher than of those on the surface of the protein. In terms of free energy, though wet-spots-forming water molecules are very heterogeneous, their contribution to the free energy of complex formation is considerable. We find that water molecules can play an important role in interaction conservation in protein interfaces by allowing sequence variability in the corresponding binding partner, and we discuss the important implications of our observations related to the use of the correlated mutations concept in protein interactions studies. The results obtained in this work help to deepen our understanding of the physico-chemical nature underlying protein-protein interactions and strengthen the idea of using the wet spots concept to qualitatively improve the accuracy of folding, docking and rational design algorithms.
Collapse
Affiliation(s)
- Sergey Samsonov
- Structural Bioinformatics, BIOTEC TU Dresden, Tatzberg 47-51, 01307 Dresden, Germany
| | | | | |
Collapse
|
25
|
Yonezawa Y, Nakata K, Sakakura K, Takada T, Nakamura H. Intra- and Intermolecular Interaction Inducing Pyramidalization on Both Sides of a Proline Dipeptide during Isomerization: An Ab Initio QM/MM Molecular Dynamics Simulation Study in Explicit Water. J Am Chem Soc 2009; 131:4535-40. [DOI: 10.1021/ja807814x] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Yasushige Yonezawa
- Institute for Protein Research, Osaka University, 3-2 Yamada-oka, Suita, Osaka 565-0871, Japan, CREST, Japan Science and Technology Agency, HPC Marketing Promotion Division, NEC Corporation, 1-10, Nisshin-Cho, Fuchu, Tokyo 183-8501, Japan, and RIKEN, Research Program for Computational Science, 2-1-1 Marunouchi, Chiyodaku, Tokyo 100-0005, Japan
| | - Kazuto Nakata
- Institute for Protein Research, Osaka University, 3-2 Yamada-oka, Suita, Osaka 565-0871, Japan, CREST, Japan Science and Technology Agency, HPC Marketing Promotion Division, NEC Corporation, 1-10, Nisshin-Cho, Fuchu, Tokyo 183-8501, Japan, and RIKEN, Research Program for Computational Science, 2-1-1 Marunouchi, Chiyodaku, Tokyo 100-0005, Japan
| | - Kota Sakakura
- Institute for Protein Research, Osaka University, 3-2 Yamada-oka, Suita, Osaka 565-0871, Japan, CREST, Japan Science and Technology Agency, HPC Marketing Promotion Division, NEC Corporation, 1-10, Nisshin-Cho, Fuchu, Tokyo 183-8501, Japan, and RIKEN, Research Program for Computational Science, 2-1-1 Marunouchi, Chiyodaku, Tokyo 100-0005, Japan
| | - Toshikazu Takada
- Institute for Protein Research, Osaka University, 3-2 Yamada-oka, Suita, Osaka 565-0871, Japan, CREST, Japan Science and Technology Agency, HPC Marketing Promotion Division, NEC Corporation, 1-10, Nisshin-Cho, Fuchu, Tokyo 183-8501, Japan, and RIKEN, Research Program for Computational Science, 2-1-1 Marunouchi, Chiyodaku, Tokyo 100-0005, Japan
| | - Haruki Nakamura
- Institute for Protein Research, Osaka University, 3-2 Yamada-oka, Suita, Osaka 565-0871, Japan, CREST, Japan Science and Technology Agency, HPC Marketing Promotion Division, NEC Corporation, 1-10, Nisshin-Cho, Fuchu, Tokyo 183-8501, Japan, and RIKEN, Research Program for Computational Science, 2-1-1 Marunouchi, Chiyodaku, Tokyo 100-0005, Japan
| |
Collapse
|
26
|
Shen T, Hamelberg D. A statistical analysis of the precision of reweighting-based simulations. J Chem Phys 2008; 129:034103. [PMID: 18647012 DOI: 10.1063/1.2944250] [Citation(s) in RCA: 98] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
Various advanced simulation techniques, which are used to sample the statistical ensemble of systems with complex Hamiltonians, such as those displayed in condensed matters and biomolecular systems, rely heavily on successfully reweighting the sampled configurations. The sampled points of a system from an elevated thermal environment or on a modified Hamiltonian are reused with different statistical weights to evaluate its properties at the initial desired temperature or of the original Hamiltonian. Often, the decrease of accuracy induced by this procedure is ignored and the final results can be far from what is expected. We have addressed the reasons behind such a phenomenon and have provided a quantitative method to estimate the number of sampled points required in the crucial step of reweighting of these advanced simulation methods. We also provided examples from temperature histogram reweighting and accelerated molecular dynamics reweighting to illustrate this idea, which can be generalized to the dynamic reweighting as well. The study shows that this analysis may provide a priori guidance for the strategy of setting up the parameters of advanced simulations before a lengthy one is carried out. The method can therefore provide insights for optimizing the parameters for high accuracy simulations with finite amount of computational resources.
Collapse
Affiliation(s)
- Tongye Shen
- Theoretical Biology & Biophysics Group, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA.
| | | |
Collapse
|
27
|
Abstract
Normal mode analysis (NMA) has received much attention as a direct approach to extract the collective motions of macromolecules. However, the stringent requirement of computational resources by classical all-atom NMA limits the size of the macromolecules to which the method is normally applied. We implemented a novel coarse-grained normal mode approach based on partitioning the all-atom Hessian matrix into relevant and nonrelevant parts. It is interesting to note that, using classical all-atom NMA results as a reference, we found that this method generates more accurate results than do other coarse-grained approaches, including elastic network model and block normal mode approaches. Moreover, this new method is effective in incorporating the energetic contributions from the nonrelevant atoms, including surface water molecules, into the coarse-grained protein motions. The importance of such improvements is demonstrated by the effect of surface water to shift vibrational modes to higher frequencies and by an increase in overlap of the coarse-grained eigenvector space (the motion directions) with that obtained from molecular dynamics simulations of solvated protein in a water box. These results not only confirm the quality of our method but also point out the importance of incorporating surface structural water in studying protein dynamics.
Collapse
|
28
|
de Oliveira CAF, Hamelberg D, McCammon JA. Estimating kinetic rates from accelerated molecular dynamics simulations: alanine dipeptide in explicit solvent as a case study. J Chem Phys 2008; 127:175105. [PMID: 17994855 DOI: 10.1063/1.2794763] [Citation(s) in RCA: 52] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Molecular dynamics (MD) simulation is the standard computational technique used to obtain information on the time evolution of the conformations of proteins and many other molecular systems. However, for most biological systems of interest, the time scale for slow conformational transitions is still inaccessible to standard MD simulations. Several sampling methods have been proposed to address this issue, including the accelerated molecular dynamics method. In this work, we study the extent of sampling of the phi/psi space of alanine dipeptide in explicit water using accelerated molecular dynamics and present a framework to recover the correct kinetic rate constant for the helix to beta-strand transition. We show that the accelerated MD can drastically enhance the sampling of the phi/psi conformational phase space when compared to normal MD. In addition, the free energy density plots of the phi/psi space show that all minima regions are accurately sampled and the canonical distribution is recovered. Moreover, the kinetic rate constant for the helix to beta-strand transition is accurately estimated from these simulations by relating the diffusion coefficient to the local energetic roughness of the energy landscape. Surprisingly, even for such a low barrier transition, it is difficult to obtain enough transitions to accurately estimate the rate constant when one uses normal MD.
Collapse
Affiliation(s)
- César Augusto F de Oliveira
- Howard Hughes Medical Institute, Center for Theoretical Biological Physics, Department of Chemistry and Biochemistry, and Department of Pharmacology, University of California at San Diego, La Jolla, California 92093, USA.
| | | | | |
Collapse
|
29
|
Hamelberg D, de Oliveira CAF, McCammon JA. Sampling of slow diffusive conformational transitions with accelerated molecular dynamics. J Chem Phys 2007; 127:155102. [PMID: 17949218 DOI: 10.1063/1.2789432] [Citation(s) in RCA: 197] [Impact Index Per Article: 11.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Abstract
Slow diffusive conformational transitions play key functional roles in biomolecular systems. Our ability to sample these motions with molecular dynamics simulation in explicit solvent is limited by the slow diffusion of the solvent molecules around the biomolecules. Previously, we proposed an accelerated molecular dynamics method that has been shown to efficiently sample the torsional degrees of freedom of biomolecules beyond the millisecond timescale. However, in our previous approach, large-amplitude displacements of biomolecules are still slowed by the diffusion of the solvent. Here we present a unified approach of efficiently sampling both the torsional degrees of freedom and the diffusive motions concurrently. We show that this approach samples the configuration space more efficiently than normal molecular dynamics and that ensemble averages converge faster to the correct values.
Collapse
Affiliation(s)
- Donald Hamelberg
- Howard Hughes Medical Institute, Center for Theoretical Biological Physics, Department of Chemistry and Biochemistry, University of California at San Diego, La Jolla, California 92093-0365, USA.
| | | | | |
Collapse
|
30
|
Murarka RK, Liwo A, Scheraga HA. Separation of time scale and coupling in the motion governed by the coarse-grained and fine degrees of freedom in a polypeptide backbone. J Chem Phys 2007; 127:155103. [DOI: 10.1063/1.2784200] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
|
31
|
Xing J. Nonequilibrium dynamic mechanism for allosteric effect. PHYSICAL REVIEW LETTERS 2007; 99:168103. [PMID: 17995300 DOI: 10.1103/physrevlett.99.168103] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/16/2007] [Indexed: 05/25/2023]
Abstract
Allosteric regulation is often viewed as thermodynamic in nature. However, protein internal motions during an enzymatic reaction cycle can slow the hoping processes over numerous potential barriers. We propose that regulating molecules may function by modifying the nonequilibrium protein dynamics. The theory predicts that an enzyme under the new mechanism has a different temperature dependence, waiting time distribution of the turnover cycle, and dynamic fluctuation patterns with and without an effector. Experimental tests of the theory are proposed.
Collapse
Affiliation(s)
- Jianhua Xing
- Chemistry, Materials and Life Sciences Directorate, University of California, Livermore, California 94550, USA.
| |
Collapse
|
32
|
Abstract
The influence of solvent on the slow internal dynamics of proteins is studied by comparing molecular dynamics simulations of solvated and unsolvated lysozyme. The dynamical trajectories are projected onto the protein's normal modes in order to obtain a separate analysis for each of the associated time scales. The results show that solvent effects are important for the slowest motions (below approximately 1 ps(-1)) but negligible for faster motions. The damping effects seen in the latter show that the principal source of friction in protein dynamics is not the solvent, but the protein itself.
Collapse
Affiliation(s)
- Konrad Hinsen
- Centre de Biophysique Moléculaire, CNRS UPR 4301, Rue Charles Sadron, 45071 Orléans Cedex 2, France.
| | | |
Collapse
|
33
|
Hamelberg D, Shen T, McCammon JA. A proposed signaling motif for nuclear import in mRNA processing via the formation of arginine claw. Proc Natl Acad Sci U S A 2007; 104:14947-51. [PMID: 17823247 PMCID: PMC1968059 DOI: 10.1073/pnas.0703151104] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
Phosphorylation of proteins by kinases is the most commonly studied class of posttranslational modification, yet its structural consequences are not well understood. The human SR (serine-arginine) protein ASF/SF2 relies on the processive phosphorylation of the serine residues of eight consecutive arginine-serine (RS) dipeptide repeats at the C terminus by SRPK1 before it can be transported into the nucleus. This SR protein plays critical roles in spliceosome assembly, pre-mRNA splicing, and mRNA export, and the phosphorylation process of the RS repeats has been extensively studied experimentally. However, knowledge of the conformational changes associated with the phosphorylation of this simple sequence and how it triggers the importation of the SR protein is lacking. Here, we have carried out extensive molecular dynamics simulations to show that phosphorylation of the eight RS repeats significantly alters the peptide's conformation and leads to the formation of very stable structures that are likely to be involved in the recognition, binding, and transport of the SR protein. Specifically, we found an unusual symmetry-broken phase of conformations of the repetitive and quasi-symmetric phosphorylated peptide sequence. One of the main characteristics of these conformations is the exposed phosphate groups on the periphery, which possibly could serve as the recognition platform for the transport protein transportin-SR2.
Collapse
Affiliation(s)
- Donald Hamelberg
- Howard Hughes Medical Institute, Center for Theoretical Biological Physics, Department of Chemistry and Biochemistry, Department of Pharmacology, University of California at San Diego, La Jolla, CA 92093-0365, USA.
| | | | | |
Collapse
|