1
|
Kunzmann P, Krumbach JH, Saponaro A, Moroni A, Thiel G, Hamacher K. Anisotropic Network Analysis of Open/Closed HCN4 Channel Advocates Asymmetric Subunit Cooperativity in cAMP Modulation of Gating. J Chem Inf Model 2024; 64:4727-4738. [PMID: 38830626 PMCID: PMC11203669 DOI: 10.1021/acs.jcim.4c00360] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2024] [Revised: 05/07/2024] [Accepted: 05/10/2024] [Indexed: 06/05/2024]
Abstract
Hyperpolarization-activated cyclic nucleotide-modulated (HCN) channels are opened in an allosteric manner by membrane hyperpolarization and cyclic nucleotides such as cAMP. Because of conflicting reports from experimental studies on whether cAMP binding to the four available binding sites in the channel tetramer operates cooperatively in gating, we employ here a computational approach as a promising route to examine ligand-induced conformational changes after binding to individual sites. By combining an elastic network model (ENM) with linear response theory (LRT) for modeling the apo-holo transition of the cyclic nucleotide-binding domain (CNBD) in HCN channels, we observe a distinct pattern of cooperativity matching the "positive-negative-positive" cooperativity reported from functional studies. This cooperativity pattern is highly conserved among HCN subtypes (HCN4, HCN1), but only to a lesser extent visible in structurally related channels, which are only gated by voltage (KAT1) or cyclic nucleotides (TAX4). This suggests an inherent cooperativity between subunits in HCN channels as part of a ligand-triggered gating mechanism in these channels.
Collapse
Affiliation(s)
- Patrick Kunzmann
- Department
of Biology, Computational Biology & Simulation, TU Darmstadt, Schnittspahnstrasse 10, 64287 Darmstadt, Germany
| | - Jan H. Krumbach
- Department
of Biology, Computational Biology & Simulation, TU Darmstadt, Schnittspahnstrasse 10, 64287 Darmstadt, Germany
| | - Andrea Saponaro
- Department
of Pharmacology and Biomolecular Sciences, University of Milano, via Balzaretti 9, 20133 Milano, Italy
| | - Anna Moroni
- Department
of Biosciences, Ion Channel Biophysics, University of Milan, via Celoria 26, 20133 Milan, Italy
| | - Gerhard Thiel
- Department
of Biology, Membrane Biophysics, TU Darmstadt, Schnittspahnstrasse 10, 64287 Darmstadt, Germany
- Centre
for Synthetic Biology, TU Darmstadt, Schnittspahnstrasse 10, 64287 Darmstadt, Germany
| | - Kay Hamacher
- Department
of Biology, Computational Biology & Simulation, TU Darmstadt, Schnittspahnstrasse 10, 64287 Darmstadt, Germany
- Centre
for Synthetic Biology, TU Darmstadt, Schnittspahnstrasse 10, 64287 Darmstadt, Germany
| |
Collapse
|
2
|
Tenenbaum A. Energy condensation and dipole alignment in protein dynamics. Phys Rev E 2024; 109:044401. [PMID: 38755822 DOI: 10.1103/physreve.109.044401] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2023] [Accepted: 02/26/2024] [Indexed: 05/18/2024]
Abstract
The possibility that distant biomolecules in a cell interact via electromagnetic (e.m.) radiation was proposed many years ago to explain the high rate of encounters of partners in some enzymatic reactions. The results of two recent experiments designed to test the propensity of protein bovine serum albumin (BSA) to interact via e.m. radiation with other proteins were interpreted in a theoretical framework based on three main assumptions: (i) in order to experience this kind of interaction the protein must be in an out-of-equilibrium state; (ii) in this state there is a condensation of energy in low-frequency vibrational modes; and (iii) the hydration layers of water around the protein sustain the energy condensation. In the present paper we present the results of molecular dynamics simulations of BSA in four states: at equilibrium and out-of-equilibrium in water, and at room and high temperature in vacuum. By comparing physical properties of the system in the four states, our simulations provide a qualitative and quantitative assessment of the three assumptions on which the theoretical framework is based. Our results confirm the assumptions of the theoretical model showing energy condensation at low frequency and electretlike alignment between the protein's and the water's dipoles; they also allow a quantitative estimate of the contribution of the out-of-equilibrium state and of the water to the observed behavior of the protein. In particular, it has been found that in the out-of-equilibrium state the amplitude of the oscillation of the protein's dipole moment greatly increases, thereby enhancing a possible absorption or emission of e.m. radiation. The analysis of BSA's dynamics outlined in the present paper provides a procedure for checking the propensity of a biomolecule to interact via e.m. radiation with its biochemical partners.
Collapse
Affiliation(s)
- Alexander Tenenbaum
- Physics Department, Sapienza University, Piazzale Aldo Moro 5, 00185 Roma, Italy
| |
Collapse
|
3
|
Gogovi GK, Silayi S, Shehu A. Computing the Structural Dynamics of RVFV L Protein Domain in Aqueous Glycerol Solutions. Biomolecules 2021; 11:biom11101427. [PMID: 34680060 PMCID: PMC8533350 DOI: 10.3390/biom11101427] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2021] [Revised: 09/24/2021] [Accepted: 09/25/2021] [Indexed: 11/23/2022] Open
Abstract
Many biological and biotechnological processes are controlled by protein–protein and protein–solvent interactions. In order to understand, predict, and optimize such processes, it is important to understand how solvents affect protein structure during protein–solvent interactions. In this study, all-atom molecular dynamics are used to investigate the structural dynamics and energetic properties of a C-terminal domain of the Rift Valley Fever Virus L protein solvated in glycerol and aqueous glycerol solutions in different concentrations by molecular weight. The Generalized Amber Force Field is modified by including restrained electrostatic potential atomic charges for the glycerol molecules. The peptide is considered in detail by monitoring properties like the root-mean-squared deviation, root-mean-squared fluctuation, radius of gyration, hydrodynamic radius, end-to-end distance, solvent-accessible surface area, intra-potential energy, and solvent–peptide interaction energies for hundreds of nanoseconds. Secondary structure analysis is also performed to examine the extent of conformational drift for the individual helices and sheets. We predict that the peptide helices and sheets are maintained only when the modeling strategy considers the solvent with lower glycerol concentration. We also find that the solvent-peptide becomes more cohesive with decreasing glycerol concentrations. The density and radial distribution function of glycerol solvent calculated when modeled with the modified atomic charges show a very good agreement with experimental results and other simulations at 298.15K.
Collapse
Affiliation(s)
- Gideon K. Gogovi
- Department of Mathematics and Statistics, University of Houston—Downtown, Houston, TX 77054, USA
- Correspondence:
| | - Swabir Silayi
- Office of Research Computing, George Mason University, Fairfax, VA 22030, USA;
| | - Amarda Shehu
- Department of Computer Science, George Mason University, Fairfax, VA 22030, USA;
- Department of Bioengineering, George Mason University, Fairfax, VA 22030, USA
- School of Systems Biology, George Mason University, Fairfax, VA 22030, USA
- Center for Advancing Human-Machine Partnerships, George Mason University, Fairfax, VA 22030, USA
| |
Collapse
|
4
|
Nakagawa H, Saio T, Nagao M, Inoue R, Sugiyama M, Ajito S, Tominaga T, Kawakita Y. Conformational dynamics of a multidomain protein by neutron scattering and computational analysis. Biophys J 2021; 120:3341-3354. [PMID: 34242590 PMCID: PMC8391080 DOI: 10.1016/j.bpj.2021.07.001] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2021] [Revised: 06/07/2021] [Accepted: 07/01/2021] [Indexed: 11/25/2022] Open
Abstract
The flexible conformations of a multidomain protein are responsible for its biological functions. Although MurD, a 47-kDa protein that consists of three domains, sequentially changes its domain conformation from an open form to a closed form through a semiclosed form in its enzymatic reaction, the domain dynamics in each conformation remains unclear. In this study, we verify the conformational dynamics of MurD in the corresponding three states (apo and ATP- and inhibitor-bound states) with a combination of small-angle x-ray and neutron scattering (SAXS and SANS), dynamic light scattering (DLS), neutron backscattering (NBS), neutron spin echo (NSE) spectroscopy, and molecular dynamics (MD) simulations. Applying principal component analysis of the MD trajectories, twisting and open-closed domain modes are identified as the major collective coordinates. The deviations of the experimental SAXS profiles from the theoretical calculations based on the known crystal structures become smaller in the ATP-bound state than in the apo state, and a further decrease is evident upon inhibitor binding. These results suggest that domain motions of the protein are suppressed step by step of each ligand binding. The DLS and NBS data yield collective and self-translational diffusion constants, respectively, and we used them to extract collective domain motions in nanometer and nanosecond scales from the NSE data. In the apo state, MurD shows both twisting and open-closed domain modes, whereas an ATP binding suppresses twisting domain motions, and a further reduction of open-closed mode is seen in the inhibitor-binding state. These observations are consistent with the structure modifications measured by the small-angle scattering as well as the MD simulations. Such changes in the domain dynamics associated with the sequential enzymatic reactions should be related to the affinity and reaction efficiency with a ligand that binds specifically to each reaction state.
Collapse
Affiliation(s)
- Hiroshi Nakagawa
- Materials Sciences Research Center, Japan Atomic Energy Agency, Tokai, Ibaraki, Japan; 2 J-PARC Center, Japan Atomic Energy Agency, Tokai, Ibaraki, Japan.
| | - Tomohide Saio
- Division of Molecular Life Science, Institute of Advanced Medical Sciences, Tokushima University, Tokushima, Japan
| | - Michihiro Nagao
- NIST Centre for Neutron Research, National Institute of Standards and Technology, Gaithersburg, Maryland; Department of Materials Science and Engineering, University of Maryland, College Park, Maryland; Department of Physics and Astronomy, University of Delaware, Newark, Delaware
| | - Rintaro Inoue
- Institute for Integrative Radiation and Nuclear Science, Kyoto University, Kumatori, Sennan-gun, Osaka, Japan
| | - Masaaki Sugiyama
- Institute for Integrative Radiation and Nuclear Science, Kyoto University, Kumatori, Sennan-gun, Osaka, Japan
| | - Satoshi Ajito
- Materials Sciences Research Center, Japan Atomic Energy Agency, Tokai, Ibaraki, Japan
| | - Taiki Tominaga
- Neutron Science and Technology Center, CROSS, Tokai, Ibaraki, Japan
| | | |
Collapse
|
5
|
Bose Majumdar A, Kim IJ, Na H. Effect of solvent on protein structure and dynamics. Phys Biol 2020; 17:036006. [DOI: 10.1088/1478-3975/ab74b3] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
6
|
Leherte L, Vercauteren DP. Reduced Point Charge Models of Proteins: Effect of Protein–Water Interactions in Molecular Dynamics Simulations of Ubiquitin Systems. J Phys Chem B 2017; 121:9771-9784. [DOI: 10.1021/acs.jpcb.7b06355] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Laurence Leherte
- Laboratoire de Physico-Chimie Informatique, Unité de Chimie Physique Théorique et Structurale, Department of Chemistry, Namur Medicine & Drug Innovation Center (NAMEDIC), Namur Institute of Structured Matter (NISM), University of Namur, Rue de Bruxelles 61, B-5000 Namur, Belgium
| | - Daniel P. Vercauteren
- Laboratoire de Physico-Chimie Informatique, Unité de Chimie Physique Théorique et Structurale, Department of Chemistry, Namur Medicine & Drug Innovation Center (NAMEDIC), Namur Institute of Structured Matter (NISM), University of Namur, Rue de Bruxelles 61, B-5000 Namur, Belgium
| |
Collapse
|
7
|
Thirumuruganandham SP, Gómez EA, Lakshmanan S, Hamblin MR. Terahertz Frequency Spectroscopy to Determine Cold Shock Protein Stability upon Solvation and Evaporation - A Molecular Dynamics Study. IEEE TRANSACTIONS ON TERAHERTZ SCIENCE AND TECHNOLOGY 2017; 7:131-143. [PMID: 30881732 PMCID: PMC6419770 DOI: 10.1109/tthz.2016.2637380] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/10/2023]
Abstract
Infrared (IR) and Terahertz (THz) spectroscopy simulations were carried out using CHARMM35b2 to determine protein stability. The stabilities of three bacterial cold shock proteins (Csps) originating from mesophiles, thermophiles and hyper- thermophiles respectively were investigated in this study. The three different Csps were investigated by Normal-Mode analysis and Molecular Dynamics simulation of THz spectra using the Hessian matrix for solvated systems, interpreted in the harmonic approximation at optimum near-melting temperatures of each homologue, by incorporating differences in the hydrous and anhydrous states of the Csps. The results show slight variations in the large scale protein motion. However, the IR spectra of Csps observed at the low frequency saddle surface region, clearly distinguishes the thermophilic and mesophilic proteins based on their stability. Further studies on protein stability employing low-frequency collective modes have the potential to reveal functionally important conformational changes that are biologically significant.
Collapse
Affiliation(s)
| | - Edgar A Gómez
- Programa de Física, Universidad del Quindío, Armenia, Colombia
| | - Shanmugamurthy Lakshmanan
- Department of Dermatology, Harvard Medical School, Boston, MA 02114, USA
- Wellman Center for Photomedicine, Massachusetts General Hospital, Boston, MA 02114, USA
| | - Michael R Hamblin
- Department of Dermatology, Harvard Medical School, Boston, MA 02114, USA
- Wellman Center for Photomedicine, Massachusetts General Hospital, Boston, MA 02114, USA
| |
Collapse
|
8
|
Biswas S, Mukherjee SK, Chowdhury PK. Crowder-Induced Rigidity in a Multidomain Protein: Insights from Solvation. J Phys Chem B 2016; 120:12501-12510. [DOI: 10.1021/acs.jpcb.6b10478] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Saikat Biswas
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| | - Sanjib Kumar Mukherjee
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| | - Pramit Kumar Chowdhury
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| |
Collapse
|
9
|
Mott AJ, Rez P. Calculation of the infrared spectra of proteins. EUROPEAN BIOPHYSICS JOURNAL: EBJ 2014; 44:103-12. [PMID: 25538002 DOI: 10.1007/s00249-014-1005-6] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/02/2014] [Revised: 10/28/2014] [Accepted: 12/05/2014] [Indexed: 11/24/2022]
Abstract
The CHARMM22 force field with associated partial charges is used to calculate the infrared spectra of a number of small proteins and some larger biothreat proteins. The calculated high-frequency region, from about 2,500 to 3,500 cm(-1), is dominated by stretching modes of hydrogen bonded to other atoms, and is very similar in all proteins. There is a peak at 3,430 cm(-1) whose intensity is predicted by these calculations to be a direct measure of arginine content. The calculated low-frequency THz region, up to 300 cm(-1), is also very similar in all the proteins and just reflects the vibrational density of states in agreement with experimental results. Calculations show that the intermediate-frequency region between 500 and 1,200 cm(-1) shows the greatest difference between individual proteins and is also the least affected by water absorption. However, to match experimental measurements in the amide region, it was necessary to reduce the hydrogen partial charges.
Collapse
Affiliation(s)
- Adam J Mott
- Department of Physics, Arizona State University, Tempe, AZ, 85287, USA
| | | |
Collapse
|
10
|
Fuglebakk E, Tiwari SP, Reuter N. Comparing the intrinsic dynamics of multiple protein structures using elastic network models. Biochim Biophys Acta Gen Subj 2014; 1850:911-922. [PMID: 25267310 DOI: 10.1016/j.bbagen.2014.09.021] [Citation(s) in RCA: 61] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2014] [Revised: 09/15/2014] [Accepted: 09/16/2014] [Indexed: 12/15/2022]
Abstract
BACKGROUND Elastic network models (ENMs) are based on the simple idea that a protein can be described as a set of particles connected by springs, which can then be used to describe its intrinsic flexibility using, for example, normal mode analysis. Since the introduction of the first ENM by Monique Tirion in 1996, several variants using coarser protein models have been proposed and their reliability for the description of protein intrinsic dynamics has been widely demonstrated. Lately an increasing number of studies have focused on the meaning of slow dynamics for protein function and its potential conservation through evolution. This leads naturally to comparisons of the intrinsic dynamics of multiple protein structures with varying levels of similarity. SCOPE OF REVIEW We describe computational strategies for calculating and comparing intrinsic dynamics of multiple proteins using elastic network models, as well as a selection of examples from the recent literature. MAJOR CONCLUSIONS The increasing interest for comparing dynamics across protein structures with various levels of similarity, has led to the establishment and validation of reliable computational strategies using ENMs. Comparing dynamics has been shown to be a viable way for gaining greater understanding for the mechanisms employed by proteins for their function. Choices of ENM parameters, structure alignment or similarity measures will likely influence the interpretation of the comparative analysis of protein motion. GENERAL SIGNIFICANCE Understanding the relation between protein function and dynamics is relevant to the fundamental understanding of protein structure-dynamics-function relationship. This article is part of a Special Issue entitled Recent developments of molecular dynamics.
Collapse
Affiliation(s)
- Edvin Fuglebakk
- Department of Molecular Biology, University of Bergen, Pb. 7803, N-5020 Bergen, Norway; Computational Biology Unit, Department of Informatics, University of Bergen, Pb. 7803, N-5020 Bergen, Norway.
| | - Sandhya P Tiwari
- Department of Molecular Biology, University of Bergen, Pb. 7803, N-5020 Bergen, Norway; Computational Biology Unit, Department of Informatics, University of Bergen, Pb. 7803, N-5020 Bergen, Norway.
| | - Nathalie Reuter
- Department of Molecular Biology, University of Bergen, Pb. 7803, N-5020 Bergen, Norway; Computational Biology Unit, Department of Informatics, University of Bergen, Pb. 7803, N-5020 Bergen, Norway.
| |
Collapse
|
11
|
Seyler SL, Beckstein O. Sampling large conformational transitions: adenylate kinase as a testing ground. MOLECULAR SIMULATION 2014. [DOI: 10.1080/08927022.2014.919497] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
12
|
Conformational response to solvent interaction and temperature of a protein (Histone h3.1) by a multi-grained monte carlo simulation. PLoS One 2013; 8:e76069. [PMID: 24204592 PMCID: PMC3799992 DOI: 10.1371/journal.pone.0076069] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2013] [Accepted: 08/19/2013] [Indexed: 12/01/2022] Open
Abstract
Interaction with the solvent plays a critical role in modulating the structure and dynamics of a protein. Because of the heterogeneity of the interaction strength, it is difficult to identify multi-scale structural response. Using a coarse-grained Monte Carlo approach, we study the structure and dynamics of a protein (H3.1) in effective solvent media. The structural response is examined as a function of the solvent-residue interaction strength (based on hydropathy index) in a range of temperatures (spanning low to high) involving a knowledge-based (Miyazawa-Jernigan(MJ)) residue-residue interaction. The protein relaxes rapidly from an initial random configuration into a quasi-static structure at low temperatures while it continues to diffuse at high temperatures with fluctuating conformation. The radius of gyration (Rg) of the protein responds non-monotonically to solvent interaction, i.e., on increasing the residue-solvent interaction strength (fs), the increase in Rg (fs≤fsc) is followed by decay (fs≥fsc) with a maximum at a characteristic value (fsc) of the interaction. Raising the temperature leads to wider spread of the distribution of the radius of gyration with higher magnitude of fsc. The effect of solvent on the multi-scale (λ: residue to Rg) structures of the protein is examined by analyzing the structure factor (S(q),|q| = 2π/λ is the wave vector of wavelength, λ) in detail. Random-coil to globular transition with temperature of unsolvated protein (H3.1) is dramatically altered by the solvent at low temperature while a systematic change in structure and scale is observed on increasing the temperature. The interaction energy profile of the residues is not sufficient to predict its mobility in the solvent. Fine-grain representation of protein with two-node and three-node residue enhances the structural resolution; results of the fine-grained simulations are consistent with the finding described above of the coarse-grained description with one-node residue.
Collapse
|
13
|
Roy U, Woods AG, Sokolowska I, Darie CC. Structural Evaluation and Analyses of Tumor Differentiation Factor. Protein J 2013; 32:512-8. [DOI: 10.1007/s10930-013-9510-5] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
14
|
Weinkam P, Sali A. Mapping polymerization and allostery of hemoglobin S using point mutations. J Phys Chem B 2013; 117:13058-68. [PMID: 23957820 DOI: 10.1021/jp4025156] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Hemoglobin is a complex system that undergoes conformational changes in response to oxygen, allosteric effectors, mutations, and environmental changes. Here, we study allostery and polymerization of hemoglobin and its variants by application of two previously described methods: (i) AllosMod for simulating allostery dynamics given two allosterically related input structures and (ii) a machine-learning method for dynamics- and structure-based prediction of the mutation impact on allostery (Weinkam et al. J. Mol. Biol. 2013, 425, 647-661), now applicable to systems with multiple coupled binding sites, such as hemoglobin. First, we predict the relative stabilities of substates and microstates of hemoglobin, which are determined primarily by entropy within our model. Next, we predict the impact of 866 annotated mutations on hemoglobin's oxygen binding equilibrium. We then discuss a subset of 30 mutations that occur in the presence of the sickle cell mutation and whose effects on polymerization have been measured. Seven of these HbS mutations occur in three predicted druggable binding pockets that might be exploited to directly inhibit polymerization; one of these binding pockets is not apparent in the crystal structure, but only in structures generated by AllosMod. For the 30 mutations, we predict that mutation-induced conformational changes within a single tetramer tend not to significantly impact polymerization; instead, these mutations more likely impact polymerization by directly perturbing a polymerization interface. Finally, our analysis of allostery allows us to hypothesize why hemoglobin evolved to have multiple subunits and a persistent low frequency sickle cell mutation.
Collapse
Affiliation(s)
- Patrick Weinkam
- Department of Bioengineering and Therapeutic Sciences, ‡Department of Pharmaceutical Chemistry, and California Institute for Quantitative Biosciences (QB3), University of California, San Francisco , San Francisco, California 94158, United States
| | | |
Collapse
|
15
|
Assessing the effect of loop mutations in the folding space of β2-microglobulin with molecular dynamics simulations. Int J Mol Sci 2013; 14:17256-78. [PMID: 23975166 PMCID: PMC3794727 DOI: 10.3390/ijms140917256] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2013] [Revised: 07/27/2013] [Accepted: 07/30/2013] [Indexed: 12/15/2022] Open
Abstract
We use molecular dynamics simulations of a full atomistic Gō model to explore the impact of selected DE-loop mutations (D59P and W60C) on the folding space of protein human β2-microglobulin (Hβ2m), the causing agent of dialysis-related amyloidosis, a conformational disorder characterized by the deposition of insoluble amyloid fibrils in the osteoarticular system. Our simulations replicate the effect of mutations on the thermal stability that is observed in experiments in vitro. Furthermore, they predict the population of a partially folded state, with 60% of native internal free energy, which is akin to a molten globule. In the intermediate state, the solvent accessible surface area increases up to 40 times relative to the native state in 38% of the hydrophobic core residues, indicating that the identified species has aggregation potential. The intermediate state preserves the disulfide bond established between residue Cys25 and residue Cys80, which helps maintain the integrity of the core region, and is characterized by having two unstructured termini. The movements of the termini dominate the essential modes of the intermediate state, and exhibit the largest displacements in the D59P mutant, which is the most aggregation prone variant. PROPKA predictions of pKa suggest that the population of the intermediate state may be enhanced at acidic pH explaining the larger amyloidogenic potential observed in vitro at low pH for the WT protein and mutant forms.
Collapse
|
16
|
Micheletti C. Comparing proteins by their internal dynamics: exploring structure-function relationships beyond static structural alignments. Phys Life Rev 2012. [PMID: 23199577 DOI: 10.1016/j.plrev.2012.10.009] [Citation(s) in RCA: 70] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
The growing interest for comparing protein internal dynamics owes much to the realisation that protein function can be accompanied or assisted by structural fluctuations and conformational changes. Analogously to the case of functional structural elements, those aspects of protein flexibility and dynamics that are functionally oriented should be subject to evolutionary conservation. Accordingly, dynamics-based protein comparisons or alignments could be used to detect protein relationships that are more elusive to sequence and structural alignments. Here we provide an account of the progress that has been made in recent years towards developing and applying general methods for comparing proteins in terms of their internal dynamics and advance the understanding of the structure-function relationship.
Collapse
Affiliation(s)
- Cristian Micheletti
- Scuola Internazionale Superiore di Studi Avanzati, via Bonomea 265, Trieste, Italy.
| |
Collapse
|
17
|
Hong L, Cheng X, Glass DC, Smith JC. Surface hydration amplifies single-well protein atom diffusion propagating into the macromolecular core. PHYSICAL REVIEW LETTERS 2012; 108:238102. [PMID: 23003993 DOI: 10.1103/physrevlett.108.238102] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/09/2012] [Indexed: 06/01/2023]
Abstract
The effect of surface hydration water on internal protein motion is of fundamental interest in molecular biophysics. Here, by decomposing the picosecond to nanosecond atomic motion in molecular dynamics simulations of lysozyme at different hydration levels into three components--localized single-well diffusion, methyl group rotation, and nonmethyl jumps--we show that the effect of surface hydration is mainly to increase the volume of the localized single-well diffusion. These diffusive motions are coupled in such a way that the hydration effect propagates from the protein surface into the dry core.
Collapse
Affiliation(s)
- Liang Hong
- UT/ORNL Center for Molecular Biophysics, Oak Ridge National Laboratory, P.O. Box 2008, Oak Ridge, Tennessee 37831, USA
| | | | | | | |
Collapse
|
18
|
Matyushov DV. Nanosecond Stokes Shift Dynamics, Dynamical Transition, and Gigantic Reorganization Energy of Hydrated Heme Proteins. J Phys Chem B 2011; 115:10715-24. [DOI: 10.1021/jp200409z] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- Dmitry V. Matyushov
- Center for Biological Physics, Arizona State University, P.O. Box 871504, Tempe, Arizona 85287-1504, United States
| |
Collapse
|
19
|
Matyushov DV, Morozov AY. Electrostatics of the protein-water interface and the dynamical transition in proteins. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2011; 84:011908. [PMID: 21867214 DOI: 10.1103/physreve.84.011908] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/25/2010] [Indexed: 05/31/2023]
Abstract
Atomic displacements of hydrated proteins are dominated by phonon vibrations at low temperatures and by dissipative large-amplitude motions at high temperatures. A crossover between the two regimes is known as a dynamical transition. Recent experiments indicate a connection between the dynamical transition and the dielectric response of the hydrated protein. We analyze two mechanisms of the coupling between the protein atomic motions and the protein-water interface. The first mechanism considers viscoelastic changes in the global shape of the protein plasticized by its coupling to the hydration shell. The second mechanism involves modulations of the local motions of partial charges inside the protein by electrostatic fluctuations. The model is used to analyze mean-square displacements of iron of metmyoglobin reported by Mössbauer spectroscopy. We show that high displacement of heme iron at physiological temperatures is dominated by electrostatic fluctuations. Two onsets, one arising from the viscoelastic response and the second from electrostatic fluctuations, are seen in the temperature dependence of the mean-square displacements when the corresponding relaxation times enter the instrumental resolution window.
Collapse
Affiliation(s)
- Dmitry V Matyushov
- Center for Biological Physics, Arizona State University, PO Box 871604, Tempe, AZ 85287-1604, USA.
| | | |
Collapse
|
20
|
Guo J, Budarz T, Ward JM, Prohofsky EW. Dynamical transition in proteins and non-Gaussian behavior of low-frequency modes in self-consistent normal mode analysis. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2010; 82:041917. [PMID: 21230323 DOI: 10.1103/physreve.82.041917] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/26/2010] [Revised: 07/22/2010] [Indexed: 05/30/2023]
Abstract
Self-consistent normal mode analysis (SCNMA) is applied to heme c type cytochrome f to study temperature-dependent protein motion. Classical normal mode analysis assumes harmonic behavior and the protein mean-square displacement has a linear dependence on temperature. This is only consistent with low-temperature experimental results. To connect the protein vibrational motions between low and physiological temperatures, we have incorporated a fitted set of anharmonic potentials into SCNMA. In addition, quantum harmonic-oscillator theory has been used to calculate the displacement distribution for individual vibrational modes. We find that the modes involving soft bonds exhibit significant non-Gaussian dynamics at physiological temperature, which suggests that it may be the cause of the non-Gaussian behavior of the protein motions probed by elastic incoherent neutron scattering. The combined theory displays a dynamical transition caused by the softening of few "torsional" modes in the low-frequency regime ( <50 cm(-1) or <6 meV or >0.6 ps). These modes change from Gaussian to a classical distribution upon heating. Our theory provides an alternative way to understand the microscopic origin of the protein dynamical transition.
Collapse
Affiliation(s)
- Jianguang Guo
- Department of Physics, Purdue University, West Lafayette, Indiana 47907, USA
| | | | | | | |
Collapse
|
21
|
Magazù S, Migliardo F, Affouard F, Descamps M, Telling MTF. Study of the relaxational and vibrational dynamics of bioprotectant glass-forming mixtures by neutron scattering and molecular dynamics simulation. J Chem Phys 2010. [DOI: 10.1063/1.3407428] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
|
22
|
Bahar I, Lezon TR, Bakan A, Shrivastava IH. Normal mode analysis of biomolecular structures: functional mechanisms of membrane proteins. Chem Rev 2010; 110:1463-97. [PMID: 19785456 PMCID: PMC2836427 DOI: 10.1021/cr900095e] [Citation(s) in RCA: 377] [Impact Index Per Article: 26.9] [Reference Citation Analysis] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Affiliation(s)
- Ivet Bahar
- Department of Computational Biology, School of Medicine, University of Pittsburgh, 3064 BST3, 3501 Fifth Avenue, Pittsburgh, Pennsylvania 15213, USA.
| | | | | | | |
Collapse
|
23
|
Feng Y, Yang L, Kloczkowski A, Jernigan RL. The energy profiles of atomic conformational transition intermediates of adenylate kinase. Proteins 2010; 77:551-8. [PMID: 19507242 DOI: 10.1002/prot.22467] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
The elastic network interpolation (ENI) (Kim et al., Biophys J 2002;83:1620-1630) is a computationally efficient and physically realistic method to generate conformational transition intermediates between two forms of a given protein. However it can be asked whether these calculated conformations provide good representatives for these intermediates. In this study, we use ENI to generate conformational transition intermediates between the open form and the closed form of adenylate kinase (AK). Based on C(alpha)-only intermediates, we construct atomic intermediates by grafting all the atoms of known AK structures onto the C(alpha) atoms and then perform CHARMM energy minimization to remove steric conflicts and optimize these intermediate structures. We compare the energy profiles for all intermediates from both the CHARMM force-field and from knowledge-based energy functions. We find that the CHARMM energies can successfully capture the two energy minima representing the open AK and closed AK forms, while the energies computed from the knowledge-based energy functions can detect the local energy minimum representing the closed AK form and show some general features of the transition pathway with a somewhat similar energy profile as the CHARMM energies. The combinatorial extension structural alignment (Shindyalov et al., 1998;11:739-747) and the k-means clustering algorithm are then used to show that known PDB structures closely resemble computed intermediates along the transition pathway.
Collapse
Affiliation(s)
- Yaping Feng
- Department of Biochemistry, Biophysics, and Molecular Biology, Iowa State University Ames, Iowa 50011-0320, USA
| | | | | | | |
Collapse
|
24
|
Matsunaga Y, Fuchigami S, Kidera A. Multivariate frequency domain analysis of protein dynamics. J Chem Phys 2009; 130:124104. [DOI: 10.1063/1.3090812] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
25
|
|
26
|
Thirumuruganandham SP, Urbassek HM. Low-frequency vibrational modes and infrared absorbance of red, blue and green opsin. J Mol Model 2009; 15:959-69. [DOI: 10.1007/s00894-008-0446-1] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2008] [Accepted: 12/23/2008] [Indexed: 11/25/2022]
|