1
|
Kalmer TL, Ancajas CMF, Cohen CI, McDaniel JM, Oyedele AS, Thirman HL, Walker AS. Statistical Coupling Analysis Predicts Correlated Motions in Dihydrofolate Reductase. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.06.18.599103. [PMID: 38948820 PMCID: PMC11213021 DOI: 10.1101/2024.06.18.599103] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/02/2024]
Abstract
The role of dynamics in enzymatic function is a highly debated topic. Dihydrofolate reductase (DHFR), due to its universality and the depth with which it has been studied, is a model system in this debate. Myriad previous works have identified networks of residues in positions near to and remote from the active site that are involved in dynamics and others that are important for catalysis. For example, specific mutations on the Met20 loop in E. coli DHFR (N23PP/S148A) are known to disrupt millisecond-timescale motions and reduce catalytic activity. However, how and if networks of dynamically coupled residues influence the evolution of DHFR is still an unanswered question. In this study, we first identify, by statistical coupling analysis and molecular dynamic simulations, a network of coevolving residues, which possess increased correlated motions. We then go on to show that allosteric communication in this network is selectively knocked down in N23PP/S148A mutant E. coli DHFR. Finally, we identify two sites in the human DHFR sector which may accommodate the Met20 loop double proline mutation while preserving dynamics. These findings strongly implicate protein dynamics as a driving force for evolution.
Collapse
Affiliation(s)
- Thomas L. Kalmer
- Department of Chemistry, Vanderbilt University Nashville, TN, USA
| | | | - Cameron I. Cohen
- Department of Biological Sciences, Vanderbilt University, Nashville, TN, USA
- Center for Structural Biology, Vanderbilt University, Nashville, TN, USA
| | - Jade M. McDaniel
- Department of Biological Sciences, Vanderbilt University, Nashville, TN, USA
| | | | - Hannah L. Thirman
- Department of Cell and Developmental Biology, Vanderbilt University, Nashville, TN, USA
- Department of Pathology, Microbiology and Immunology, Vanderbilt University Medical Center, Nashville, TN, USA
- Vanderbilt Center for Immunobiology, Vanderbilt University Medical Center, Nashville, TN, USA
- Chemical & Physical Biology Program, Vanderbilt University, Nashville, TN, USA
| | - Allison S. Walker
- Department of Chemistry, Vanderbilt University Nashville, TN, USA
- Department of Biological Sciences, Vanderbilt University, Nashville, TN, USA
- Department of Pathology, Microbiology and Immunology, Vanderbilt University Medical Center, Nashville, TN, USA
- Evolutionary Studies Initiative, Vanderbilt University, Nashville, TN, USA
| |
Collapse
|
2
|
Li M, Bao Y, Xu R, Zhang X, La H, Guo J. Mechanism of enhanced sensitivity of mutated β-adrenergic-like octopamine receptor to amitraz in honeybee Apis mellifera: An insight from MD simulations. PEST MANAGEMENT SCIENCE 2022; 78:5423-5431. [PMID: 36057136 DOI: 10.1002/ps.7164] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/02/2022] [Revised: 08/29/2022] [Accepted: 09/03/2022] [Indexed: 06/15/2023]
Abstract
BACKGROUND Amitraz is one of the critical acaricides/insecticides for effective control of pest infestation of Varroa destructor mite, a devastating parasite of Apis mellifera, because of its low toxicity to honeybees. Previous assays verified that a typical G protein-coupled receptor, β-adrenergic-like octopamine receptor (Octβ2R), is the unique target of amitraz, but the honeybee Octβ2R resists to amitraz. However, the underlying molecular mechanism of the enhanced sensitivity or toxicity of amitraz to mutated honeybee Octβ2RE208V/I335T/I350V is not fully understood. Here, molecular dynamics simulations are employed to explore the implied mechanism of the enhanced sensitivity to amitraz in mutant honeybee Octβ2R. RESULTS We found that amitraz binding stabilized the structure of Octβ2R, particularly the intracellular loop 3 associated with the Octβ2R signaling. Then, it was further demonstrated that both mutations and ligand binding resulted in a more rigid and compact amitraz binding site, as well as the outward movement of the transmembrane helix 6, which was a prerequisite for G protein coupling and activation. Moreover, mutations were found to promote the binding between Octβ2R and amitraz. Finally, community analysis illuminated that mutations and amitraz strengthened the residue-residue communication within the transmembrane domain, which might facilitate the allosteric signal propagation and activation of Octβ2R. CONCLUSION Our results unveiled structural determinants of improved sensitivity in the Octβ2R-amitraz complex and may contribute to further structure-based drug design for safer and less toxic selective insecticides. © 2022 Society of Chemical Industry.
Collapse
Affiliation(s)
- Mengrong Li
- College of Life Sciences, Nanjing Agricultural University, Nanjing, China
| | - Yiqiong Bao
- College of Life Sciences, Nanjing Agricultural University, Nanjing, China
| | - Ran Xu
- College of Life Sciences, Nanjing Agricultural University, Nanjing, China
| | - Xiaoxiao Zhang
- College of Life Sciences, Nanjing Agricultural University, Nanjing, China
| | - Honggui La
- College of Life Sciences, Nanjing Agricultural University, Nanjing, China
| | - Jingjing Guo
- College of Life Sciences, Nanjing Agricultural University, Nanjing, China
- Engineering Research Centre of Applied Technology on Machine Translation and Artificial Intelligence, Centre in Artificial Intelligence Driven Drug Discovery, Faculty of Applied Science, Macao Polytechnic University, Macao, China
| |
Collapse
|
3
|
Ali AAAI, Gulzar A, Wolf S, Stock G. Nonequilibrium Modeling of the Elementary Step in PDZ3 Allosteric Communication. J Phys Chem Lett 2022; 13:9862-9868. [PMID: 36251493 DOI: 10.1021/acs.jpclett.2c02821] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
While allostery is of paramount importance for protein signaling and regulation, the underlying dynamical process of allosteric communication is not well understood. The PDZ3 domain represents a prime example of an allosteric single-domain protein, as it features a well-established long-range coupling between the C-terminal α3-helix and ligand binding. In an intriguing experiment, Hamm and co-workers employed photoswitching of the α3-helix to initiate a conformational change of PDZ3 that propagates from the C-terminus to the bound ligand within 200 ns. Performing extensive nonequilibrium molecular dynamics simulations, the modeling of the experiment reproduces the measured time scales and reveals a detailed picture of the allosteric communication in PDZ3. In particular, a correlation analysis identifies a network of contacts connecting the α3-helix and the core of the protein, which move in a concerted manner. Representing a one-step process and involving direct α3-ligand contacts, this cooperative transition is considered as the elementary step in the propagation of conformational change.
Collapse
Affiliation(s)
- Ahmed A A I Ali
- Biomolecular Dynamics, Institute of Physics, University of Freiburg, 79104Freiburg, Germany
| | - Adnan Gulzar
- Biomolecular Dynamics, Institute of Physics, University of Freiburg, 79104Freiburg, Germany
| | - Steffen Wolf
- Biomolecular Dynamics, Institute of Physics, University of Freiburg, 79104Freiburg, Germany
| | - Gerhard Stock
- Biomolecular Dynamics, Institute of Physics, University of Freiburg, 79104Freiburg, Germany
| |
Collapse
|
4
|
Haliloglu T, Hacisuleyman A, Erman B. Prediction of Allosteric Communication Pathways in Proteins. Bioinformatics 2022; 38:3590-3599. [PMID: 35674396 DOI: 10.1093/bioinformatics/btac380] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2021] [Revised: 04/12/2022] [Accepted: 06/01/2022] [Indexed: 12/24/2022] Open
Abstract
MOTIVATION Allostery in proteins is an essential phenomenon in biological processes. In this paper, we present a computational model to predict paths of maximum information transfer between active and allosteric sites. In this information theoretic study, we use mutual information as the measure of information transfer, where transition probability of information from one residue to its contacting neighbors is proportional to the magnitude of mutual information between the two residues. Starting from a given residue and using a Hidden Markov Model, we successively determine the neighboring residues that eventually lead to a path of optimum information transfer. The Gaussian approximation of mutual information between residue pairs is adopted. The limits of validity of this approximation are discussed in terms of a nonlinear theory of mutual information and its reduction to the Gaussian form. RESULTS Predictions of the model are tested on six widely studied cases, CheY Bacterial Chemotaxis, B-cell Lymphoma extra-large Bcl-xL, Human proline isomerase cyclophilin A (CypA), Dihydrofolate reductase DHFR, HRas GTPase, and Caspase-1. The communication transmission rendering the propagation of local fluctuations from the active sites throughout the structure in multiple paths correlate well with the known experimental data. Distinct paths originating from the active site may likely represent a multi functionality such as involving more than one allosteric site and/or preexistence of some other functional states. Our model is computationally fast and simple, and can give allosteric communication pathways, which are crucial for the understanding and control of protein functionality. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Turkan Haliloglu
- Polymer Research Center and Chemical Engineering Department, Bogazici University, 34342, Turkey
| | - Aysima Hacisuleyman
- Institute of Bioengineering, Swiss Federal Institute of Technology (EPFL), 1015, Switzerland
| | - Burak Erman
- Chemical and Biological Engineering, Koc University, 34450, Turkey
| |
Collapse
|
5
|
Krucinska J, Lombardo MN, Erlandsen H, Estrada A, Si D, Viswanathan K, Wright DL. Structure-guided functional studies of plasmid-encoded dihydrofolate reductases reveal a common mechanism of trimethoprim resistance in Gram-negative pathogens. Commun Biol 2022; 5:459. [PMID: 35562546 PMCID: PMC9106665 DOI: 10.1038/s42003-022-03384-y] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2021] [Accepted: 04/20/2022] [Indexed: 11/20/2022] Open
Abstract
Two plasmid-encoded dihydrofolate reductase (DHFR) isoforms, DfrA1 and DfrA5, that give rise to high levels of resistance in Gram-negative bacteria were structurally and biochemically characterized to reveal the mechanism of TMP resistance and to support phylogenic groupings for drug development against antibiotic resistant pathogens. Preliminary screening of novel antifolates revealed related chemotypes that showed high levels of inhibitory potency against Escherichia coli chromosomal DHFR (EcDHFR), DfrA1, and DfrA5. Kinetics and biophysical analysis, coupled with crystal structures of trimethoprim bound to EcDHFR, DfrA1 and DfrA5, and two propargyl-linked antifolates (PLA) complexed with EcDHFR, DfrA1 and DfrA5, were determined to define structural features of the substrate binding pocket and guide synthesis of pan-DHFR inhibitors. Critical residue variations in two of the most clinically prevalent DHFR isoforms are identified as a common structural element in trimethoprim-resistant DHFR which impose changes on enzyme catalysis and ligand-cofactor cooperativity.
Collapse
Affiliation(s)
- Jolanta Krucinska
- Department of Pharmaceutical Sciences, University of Connecticut, 69N. Eagleville Rd., Storrs, CT, 06269, USA
| | - Michael N Lombardo
- Department of Pharmaceutical Sciences, University of Connecticut, 69N. Eagleville Rd., Storrs, CT, 06269, USA
| | - Heidi Erlandsen
- Center for Open Research Resources & Equipment (COR2E), University of Connecticut, 91N. Eagleville Rd., Storrs, CT, 06269, USA
| | - Alexavier Estrada
- Department of Pharmaceutical Sciences, University of Connecticut, 69N. Eagleville Rd., Storrs, CT, 06269, USA
| | - Debjani Si
- Department of Pharmaceutical Sciences, University of Connecticut, 69N. Eagleville Rd., Storrs, CT, 06269, USA
| | - Kishore Viswanathan
- Department of Pharmaceutical Sciences, University of Connecticut, 69N. Eagleville Rd., Storrs, CT, 06269, USA
| | - Dennis L Wright
- Department of Pharmaceutical Sciences, University of Connecticut, 69N. Eagleville Rd., Storrs, CT, 06269, USA.
| |
Collapse
|
6
|
Hasnain S, Mugnai ML, Thirumalai D. Effects of Gold Nanoparticles on the Stepping Trajectories of Kinesin. J Phys Chem B 2021; 125:10432-10444. [PMID: 34499499 DOI: 10.1021/acs.jpcb.1c02218] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
A substantial increase in the temporal resolution of the stepping of dimeric molecular motors is possible by tracking the position of a large gold nanoparticle (GNP) attached to a labeled site on one of the heads. This technique was employed to measure the stepping trajectories of conventional kinesin (Kin1) using the time-dependent position of the GNP as a proxy. The trajectories revealed that the detached head always passes to the right of the head that is tightly bound to the microtubule (MT) during a step. In interpreting the results of such experiments, it is assumed that the GNP does not significantly alter the diffusive motion of the detached head. We used coarse-grained simulations of a system consisting of the MT-Kin1 complex with and without attached GNP to investigate how the stepping trajectories are affected. The two significant findings are: (1) The GNP does not faithfully track the position of the stepping head, and (2) the rightward bias is typically exaggerated by the GNP. Both these findings depend on the precise residue position to which the GNP is attached. Surprisingly, the stepping trajectories of kinesin are not significantly affected if, in addition to the GNP, a 1 μm diameter cargo is attached to the coiled coil. Our simulations suggest the effects of the large probe have to be considered when inferring the stepping mechanisms using GNP tracking experiments.
Collapse
Affiliation(s)
- Sabeeha Hasnain
- Department of Chemistry, The University of Texas at Austin, Austin 78712, Texas, United States
| | - Mauro L Mugnai
- Department of Chemistry, The University of Texas at Austin, Austin 78712, Texas, United States
| | - D Thirumalai
- Department of Chemistry, The University of Texas at Austin, Austin 78712, Texas, United States
| |
Collapse
|
7
|
Civera M, Moroni E, Sorrentino L, Vasile F, Sattin S. Chemical and Biophysical Approaches to Allosteric Modulation. European J Org Chem 2021. [DOI: 10.1002/ejoc.202100506] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Monica Civera
- Department of Chemistry Università degli Studi di Milano via C. Golgi, 19 20133 Milan Italy
| | - Elisabetta Moroni
- Istituto di Scienze e Tecnologie Chimiche Giulio Natta, SCITEC Via Mario Bianco 9 20131 Milan Italy
| | - Luca Sorrentino
- Department of Chemistry Università degli Studi di Milano via C. Golgi, 19 20133 Milan Italy
| | - Francesca Vasile
- Department of Chemistry Università degli Studi di Milano via C. Golgi, 19 20133 Milan Italy
| | - Sara Sattin
- Department of Chemistry Università degli Studi di Milano via C. Golgi, 19 20133 Milan Italy
| |
Collapse
|
8
|
Huang BC, Yang LW. Molecular dynamics simulations and linear response theories jointly describe biphasic responses of myoglobin relaxation and reveal evolutionarily conserved frequent communicators. Biophys Physicobiol 2020; 16:473-484. [PMID: 31984199 PMCID: PMC6975898 DOI: 10.2142/biophysico.16.0_473] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2019] [Accepted: 09/20/2019] [Indexed: 12/01/2022] Open
Abstract
In this study, we provide a time-dependent mechanical model, taking advantage of molecular dynamics simulations, quasiharmonic analysis of molecular dynamics trajectories, and time-dependent linear response theories to describe vibrational energy redistribution within the protein matrix. The theoretical description explained the observed biphasic responses of specific residues in myoglobin to CO-photolysis and photoexcitation on heme. The fast responses were found to be triggered by impulsive forces and propagated mainly by principal modes <40 cm−1. The predicted fast responses for individual atoms were then used to study signal propagation within the protein matrix and signals were found to propagate ~8 times faster across helices (4076 m/s) than within the helices, suggesting the importance of tertiary packing in the sensitivity of proteins to external perturbations. We further developed a method to integrate multiple intramolecular signal pathways and discover frequent “communicators”. These communicators were found to be evolutionarily conserved including those distant from the heme.
Collapse
Affiliation(s)
- Bang-Chieh Huang
- Institute of Bioinformatics and Structural Biology, National Tsing-Hua University, Hsinchu 30013, Taiwan
| | - Lee-Wei Yang
- Institute of Bioinformatics and Structural Biology, National Tsing-Hua University, Hsinchu 30013, Taiwan.,Bioinformatics Program, Taiwan International Graduate Program, Institute of Information Sciences, Academia Sinica, Taipei 11529, Taiwan.,Department of Life Science, National Tsing Hua University, Hsinchu 30013, Taiwan.,Physics Division, National Center for Theoretical Sciences, Hsinchu 30013, Taiwan
| |
Collapse
|
9
|
Thirumalai D, Hyeon C, Zhuravlev PI, Lorimer GH. Symmetry, Rigidity, and Allosteric Signaling: From Monomeric Proteins to Molecular Machines. Chem Rev 2019; 119:6788-6821. [DOI: 10.1021/acs.chemrev.8b00760] [Citation(s) in RCA: 51] [Impact Index Per Article: 10.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- D. Thirumalai
- Department of Chemistry, The University of Texas, Austin, Texas 78712, United States
| | - Changbong Hyeon
- Korea Institute for Advanced Study, Seoul 02455, Republic of Korea
| | - Pavel I. Zhuravlev
- Biophysics Program, Institute for Physical Science and Technology and Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742, United States
| | - George H. Lorimer
- Biophysics Program, Institute for Physical Science and Technology and Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742, United States
| |
Collapse
|
10
|
Nagai T, Tama F, Miyashita O. Cryo-Cooling Effect on DHFR Crystal Studied by Replica-Exchange Molecular Dynamics Simulations. Biophys J 2019; 116:395-405. [PMID: 30638963 DOI: 10.1016/j.bpj.2018.11.3139] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2018] [Revised: 11/12/2018] [Accepted: 11/26/2018] [Indexed: 11/25/2022] Open
Abstract
Cryo-cooling is routinely performed before x-ray diffraction image collection to reduce the damage to crystals due to ionizing radiation. It has been suggested that although backbone structures are usually very similar between room temperature and cryo-temperature, cryo-cooling may hamper biologically relevant dynamics. In this study, the crystal of Escherichia coli dihydrofolate reductase is studied with replica-exchange molecular dynamics simulation, and the results are compared with the crystal structure determined at cryo-temperature and room temperature with the time-averaged ensemble method. Although temperature dependence of unit cell compaction and root mean-square fluctuation of Cα is found in accord with experiment, it is found that the protein structure at low temperature can be more heterogeneous than the ensemble of structures reported by using the time-averaged ensemble method, encouraging further development of the time-averaged ensemble method and indicating that data should be examined carefully to avoid overinterpretation of one average structure.
Collapse
Affiliation(s)
- Tetsuro Nagai
- Department of Physics, Graduate School of Science, Nagoya University, Nagoya, Aichi, Japan
| | - Florence Tama
- Department of Physics, Graduate School of Science, Nagoya University, Nagoya, Aichi, Japan; ITbM, Nagoya University, Nagoya, Aichi, Japan; RIKEN Center for Computational Science, Kobe, Hyogo, Japan
| | | |
Collapse
|
11
|
Yu M, Chen Y, Wang ZL, Liu Z. Fluctuation correlations as major determinants of structure- and dynamics-driven allosteric effects. Phys Chem Chem Phys 2019; 21:5200-5214. [DOI: 10.1039/c8cp07859a] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Both structure- and dynamics-driven allosteric effects are determined by the correlation of distance fluctuations in proteins.
Collapse
Affiliation(s)
- Miao Yu
- College of Chemistry and Molecular Engineering
- Peking University
- Beijing 100871
- China
| | - Yixin Chen
- College of Chemistry and Molecular Engineering
- Peking University
- Beijing 100871
- China
| | - Zi-Le Wang
- Department of Physics
- Tsinghua University
- Beijing 100084
- China
| | - Zhirong Liu
- College of Chemistry and Molecular Engineering
- Peking University
- Beijing 100871
- China
- Center for Quantitative Biology
| |
Collapse
|
12
|
Srinivasan B, Tonddast-Navaei S, Roy A, Zhou H, Skolnick J. Chemical space of Escherichia coli dihydrofolate reductase inhibitors: New approaches for discovering novel drugs for old bugs. Med Res Rev 2018; 39:684-705. [PMID: 30192413 DOI: 10.1002/med.21538] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2018] [Revised: 07/16/2018] [Accepted: 08/09/2018] [Indexed: 12/15/2022]
Abstract
Escherichia coli Dihydrofolate reductase is an important enzyme that is essential for the survival of the Gram-negative microorganism. Inhibitors designed against this enzyme have demonstrated application as antibiotics. However, either because of poor bioavailability of the small-molecules resulting from their inability to cross the double membrane in Gram-negative bacteria or because the microorganism develops resistance to the antibiotics by mutating the DHFR target, discovery of new antibiotics against the enzyme is mandatory to overcome drug-resistance. This review summarizes the field of DHFR inhibition with special focus on recent efforts to effectively interface computational and experimental efforts to discover novel classes of inhibitors that target allosteric and active-sites in drug-resistant variants of EcDHFR.
Collapse
Affiliation(s)
- Bharath Srinivasan
- Center for the Study of Systems Biology, School of Biology, Georgia Institute of Technology, Atlanta, Georgia
| | - Sam Tonddast-Navaei
- Center for the Study of Systems Biology, School of Biology, Georgia Institute of Technology, Atlanta, Georgia
| | - Ambrish Roy
- Center for the Study of Systems Biology, School of Biology, Georgia Institute of Technology, Atlanta, Georgia
| | - Hongyi Zhou
- Center for the Study of Systems Biology, School of Biology, Georgia Institute of Technology, Atlanta, Georgia
| | - Jeffrey Skolnick
- Center for the Study of Systems Biology, School of Biology, Georgia Institute of Technology, Atlanta, Georgia
| |
Collapse
|
13
|
Chen J, Thirumalai D. Interface Residues That Drive Allosteric Transitions Also Control the Assembly of l-Lactate Dehydrogenase. J Phys Chem B 2018; 122:11195-11205. [PMID: 30102042 DOI: 10.1021/acs.jpcb.8b06430] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
The allosteric enzyme, l-lactate dehydrogenase (LDH), is activated by fructose 1,6-metaphosphate (FBP) to reduce pyruvate to lactate. The molecular details of the FBP-driven transition from the low affinity T state to the high affinity R state in LDH, a tetramer composed of identical subunits, are not known. The dynamics of the T → R allosteric transition, investigated using Brownian dynamics (BD) simulations of the self-organized polymer (SOP) model, revealed that coordinated rotations of the subunits drive the T → R transition. We used the structural perturbation method (SPM), which requires only the static structure, to identify the allostery wiring diagram (AWD), a network of residues that transmits signals across the tetramer, as LDH undergoes the T → R transition. Interestingly, the residues that play a major role in the dynamics, which are predominantly localized at the interfaces, coincide with the AWD identified using the SPM. Although the allosteric pathways are highly heterogeneous, on the basis of our simulations, we surmise that predominantly the conformational changes in the T → R transition start from the region near the active site, comprised of helix αC, helix α1/2G, helix α3G, and helix α2F, and proceed to other structural units, thus completing the global motion. Brownian dynamics simulations of the tetramer assembly, triggered by a temperature quench from the fully disrupted conformations, show that the bottleneck for assembly is the formation of the correct orientational registry between the subunits, requiring contacts between the interface residues. Surprisingly, these residues are part of the AWD, which was identified using the SPM. Taken together, our results show that LDH, and perhaps other multidomain proteins, may have evolved to stabilize distinct states of allosteric enzymes using precisely the same AWD that also controls the functionally relevant allosteric transitions.
Collapse
Affiliation(s)
- Jie Chen
- Biophysics Program, Institute for Physical Science and Technology , University of Maryland , College Park , Maryland 20742 , United States
| | - D Thirumalai
- Department of Chemistry , The University of Texas , Austin , Texas 78712 , United States
| |
Collapse
|
14
|
Mhashal AR, Pshetitsky Y, Eitan R, Cheatum CM, Kohen A, Major DT. Effect of Asp122 Mutation on the Hydride Transfer in E. coli DHFR Demonstrates the Goldilocks of Enzyme Flexibility. J Phys Chem B 2018; 122:8006-8017. [PMID: 30040418 DOI: 10.1021/acs.jpcb.8b05556] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
Dihydrofolate reductase (DHFR) catalyzes the reduction of dihydrofolate (DHF) to tetrahydrofolate (THF) in the presence of NADPH. The key hydride transfer step in the reaction is facilitated by a combination of enzyme active site preorganization and correlated protein motions in the Michaelis-Menten (E:NADPH:DHF) complex. The present theoretical study employs mutagenesis to examine the relation between structural and functional properties of the enzyme. We mutate Asp122 in Escherichia coli DHFR, which is a conserved amino acid in the DHFR family. The consequent effect of the mutation on enzyme catalysis is examined from an energetic, structural and short-time dynamic perspective. Our investigations suggest that the structural and short-time dynamic perturbations caused by Asp122X mutations (X = Asn, Ser, Ala) are along the reaction coordinate and lower the rate of hydride transfer. Importantly, analysis of the correlated and principle component motions in the enzyme suggest that the mutation alters the coupled motions that are present in the wild-type enzyme. In the case of D122N and D122S, the mutations inhibit coupled motion, whereas in the case of D122A, the mutation enhances coupled motion, although all mutations result in similar rate reduction. These results emphasize a Goldilocks principle of enzyme flexibility, that is, enzymes should neither be too rigid nor too flexible.
Collapse
Affiliation(s)
- Anil R Mhashal
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry , Bar-Ilan University , Ramat-Gan 52900 , Israel
| | - Yaron Pshetitsky
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry , Bar-Ilan University , Ramat-Gan 52900 , Israel
| | - Reuven Eitan
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry , Bar-Ilan University , Ramat-Gan 52900 , Israel
| | - Christopher M Cheatum
- Department of Chemistry , University of Iowa , Iowa City , Iowa 52242 , United States
| | - Amnon Kohen
- Department of Chemistry , University of Iowa , Iowa City , Iowa 52242 , United States
| | - Dan Thomas Major
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry , Bar-Ilan University , Ramat-Gan 52900 , Israel
| |
Collapse
|
15
|
Parsing the roles of neck-linker docking and tethered head diffusion in the stepping dynamics of kinesin. Proc Natl Acad Sci U S A 2017; 114:E9838-E9845. [PMID: 29087307 DOI: 10.1073/pnas.1706014114] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
Kinesin walks processively on microtubules (MTs) in an asymmetric hand-over-hand manner consuming one ATP molecule per 16-nm step. The individual contributions due to docking of the approximately 13-residue neck linker to the leading head (deemed to be the power stroke) and diffusion of the trailing head (TH) that contributes in propelling the motor by 16 nm have not been quantified. We use molecular simulations by creating a coarse-grained model of the MT-kinesin complex, which reproduces the measured stall force as well as the force required to dislodge the motor head from the MT, to show that nearly three-quarters of the step occurs by bidirectional stochastic motion of the TH. However, docking of the neck linker to the leading head constrains the extent of diffusion and minimizes the probability that kinesin takes side steps, implying that both the events are necessary in the motility of kinesin and for the maintenance of processivity. Surprisingly, we find that during a single step, the TH stochastically hops multiple times between the geometrically accessible neighboring sites on the MT before forming a stable interaction with the target binding site with correct orientation between the motor head and the [Formula: see text] tubulin dimer.
Collapse
|
16
|
Abstract
Allosteric transition, defined as conformational changes induced by ligand binding, is one of the fundamental properties of proteins. Allostery has been observed and characterized in many proteins, and has been recently utilized to control protein function via regulation of protein activity. Here, we review the physical and evolutionary origin of protein allostery, as well as its importance to protein regulation, drug discovery, and biological processes in living systems. We describe recently developed approaches to identify allosteric pathways, connected sets of pairwise interactions that are responsible for propagation of conformational change from the ligand-binding site to a distal functional site. We then present experimental and computational protein engineering approaches for control of protein function by modulation of allosteric sites. As an example of application of these approaches, we describe a synergistic computational and experimental approach to rescue the cystic-fibrosis-associated protein cystic fibrosis transmembrane conductance regulator, which upon deletion of a single residue misfolds and causes disease. This example demonstrates the power of allosteric manipulation in proteins to both elucidate mechanisms of molecular function and to develop therapeutic strategies that rescue those functions. Allosteric control of proteins provides a tool to shine a light on the complex cascades of cellular processes and facilitate unprecedented interrogation of biological systems.
Collapse
Affiliation(s)
- Nikolay V Dokholyan
- Department of Biochemistry and Biophysics, University of North Carolina , Chapel Hill, North Carolina 27599, United States
| |
Collapse
|
17
|
Papaleo E, Saladino G, Lambrughi M, Lindorff-Larsen K, Gervasio FL, Nussinov R. The Role of Protein Loops and Linkers in Conformational Dynamics and Allostery. Chem Rev 2016; 116:6391-423. [DOI: 10.1021/acs.chemrev.5b00623] [Citation(s) in RCA: 239] [Impact Index Per Article: 29.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Elena Papaleo
- Computational
Biology Laboratory, Unit of Statistics, Bioinformatics and Registry, Danish Cancer Society Research Center, Strandboulevarden 49, 2100 Copenhagen, Denmark
- Structural
Biology and NMR Laboratory, Department of Biology, University of Copenhagen, 2200 Copenhagen, Denmark
| | - Giorgio Saladino
- Department
of Chemistry, University College London, London WC1E 6BT, United Kingdom
| | - Matteo Lambrughi
- Department
of Biotechnology and Biosciences, University of Milano-Bicocca, Piazza
della Scienza 2, 20126 Milan, Italy
| | - Kresten Lindorff-Larsen
- Structural
Biology and NMR Laboratory, Department of Biology, University of Copenhagen, 2200 Copenhagen, Denmark
| | | | - Ruth Nussinov
- Cancer
and Inflammation Program, Leidos Biomedical Research, Inc., Frederick
National Laboratory for Cancer Research, National Cancer Institute Frederick, Frederick, Maryland 21702, United States
- Sackler Institute
of Molecular Medicine, Department of Human Genetics and Molecular
Medicine Sackler School of Medicine, Tel Aviv University, Tel Aviv 69978, Israel
| |
Collapse
|
18
|
Effects of Non-Natural Amino Acid Incorporation into the Enzyme Core Region on Enzyme Structure and Function. Int J Mol Sci 2015; 16:22735-53. [PMID: 26402667 PMCID: PMC4613333 DOI: 10.3390/ijms160922735] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2015] [Revised: 09/02/2015] [Accepted: 09/14/2015] [Indexed: 11/17/2022] Open
Abstract
Techniques to incorporate non-natural amino acids (NNAAs) have enabled biosynthesis of proteins containing new building blocks with unique structures, chemistry, and reactivity that are not found in natural amino acids. It is crucial to understand how incorporation of NNAAs affects protein function because NNAA incorporation may perturb critical function of a target protein. This study investigates how the site-specific incorporation of NNAAs affects catalytic properties of an enzyme. A NNAA with a hydrophobic and bulky sidechain, 3-(2-naphthyl)-alanine (2Nal), was site-specifically incorporated at six different positions in the hydrophobic core of a model enzyme, murine dihydrofolate reductase (mDHFR). The mDHFR variants with a greater change in van der Waals volume upon 2Nal incorporation exhibited a greater reduction in the catalytic efficiency. Similarly, the steric incompatibility calculated using RosettaDesign, a protein stability calculation program, correlated with the changes in the catalytic efficiency.
Collapse
|
19
|
Xiong L, Liu Z. Molecular dynamics study on folding and allostery in RfaH. Proteins 2015; 83:1582-92. [DOI: 10.1002/prot.24839] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2015] [Revised: 05/18/2015] [Accepted: 05/22/2015] [Indexed: 12/18/2022]
Affiliation(s)
- Liqin Xiong
- Department of Physics; Beijing Normal University; Beijing 100875 China
| | - Zhenxing Liu
- Department of Physics; Beijing Normal University; Beijing 100875 China
| |
Collapse
|
20
|
Both protein dynamics and ligand concentration can shift the binding mechanism between conformational selection and induced fit. Proc Natl Acad Sci U S A 2014; 111:10197-202. [PMID: 24982141 DOI: 10.1073/pnas.1407545111] [Citation(s) in RCA: 70] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
This study aimed to shed light on the long debate over whether conformational selection (CS) or induced fit (IF) is the governing mechanism for protein-ligand binding. The main difference between the two scenarios is whether the conformational transition of the protein from the unbound form to the bound form occurs before or after encountering the ligand. Here we introduce the IF fraction (i.e., the fraction of binding events achieved via IF), to quantify the binding mechanism. Using simulations of a model protein-ligand system, we demonstrate that both the rate of the conformational transition and the concentration of ligand molecules can affect the IF fraction. CS dominates at slow conformational transition and low ligand concentration. An increase in either quantity results in a higher IF fraction. Despite the many-body nature of the system and the involvement of multiple, disparate types of dynamics (i.e., ligand diffusion, protein conformational transition, and binding reaction), the overall binding kinetics over wide ranges of parameters can be fit to a single exponential, with the apparent rate constant exhibiting a linear dependence on ligand concentration. The present study may guide future kinetics experiments and dynamics simulations in determining the IF fraction.
Collapse
|
21
|
Iakhiaev MA, Iakhiaev AV. Mapping the intramolecular signal propagation pathways in protein using Bayesian change point analysis of atomic motions. Comput Biol Chem 2013; 47:89-95. [PMID: 24025705 DOI: 10.1016/j.compbiolchem.2013.08.004] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2013] [Revised: 05/27/2013] [Accepted: 08/09/2013] [Indexed: 11/29/2022]
Abstract
We propose to use change points of atomic positions in the molecular dynamics trajectory as indicators of the propagating signals in protein. We designate these changes as signals because they can propagate within the molecule in the form of "perturbation wave", transmit energy or information between different parts of protein, and serve as allosteric signals. We found that change points can distinguish between thermal fluctuations of atoms (noise) and signals in a protein despite the differences in the motility of amino acid residues. Clustering of the spatially close residues that were experiencing change points close in time, allowed us to map pathways of signal propagation in a protein at the atomic level of resolution. We propose a potential mechanism for the origin of the signal and its propagation that relies on the autonomic coherence resonance in atomic fluctuations. According to this mechanism, random synchronization of fluctuations of neighboring atoms results in a resonance, which increases amplitude of vibration of these atoms. This increase can be transmitted to the atoms colliding with the resonant atoms, leading to the propagating signal. The wavelet-based coherence analysis of the inter-atomic distances between carbon-alpha atoms and surrounding atoms for the residue pairs that belong to the same communication pathway allowed us to find time periods with temporarily locked phases, confirming the occurrence of conditions for resonance. Analysis of the mapped pathways demonstrated that they form a network that connects different regions of the protein.
Collapse
|
22
|
Abstract
The spatial structure of the cell is highly organized at all levels: from small complexes and assemblies, to local nano- and microclusters, to global, micrometer scales across and between cells. We suggest that this multiscale spatial cell organization also organizes signaling and coordinates cellular behavior. We propose a new view of the spatial structure of cell signaling systems. This new view describes cell signaling in terms of dynamic allosteric interactions within and among distinct, spatially organized transient clusters. The clusters vary over time and space and are on length scales from nanometers to micrometers. When considered across these length scales, primary factors in the spatial organization are cell membrane domains and the actin cytoskeleton, both also highly dynamic. A key challenge is to understand the interplay across these multiple scales, link it to the physicochemical basis of the conformational behavior of single molecules and ultimately relate it to cellular function. Overall, our premise is that at these scales, cell signaling should be thought of not primarily as a sequence of diffusion-controlled molecular collisions, but instead transient, allostery-driven cluster re-forming interactions.
Collapse
Affiliation(s)
- Ruth Nussinov
- Basic Research Program, SAIC-Frederick, Inc., Cancer and Inflammation Program, National Cancer Institute, Frederick, MD 21702
- Sackler Inst. of Molecular Medicine, Department of Human Genetics and Molecular Medicine, Sackler School of Medicine, Tel Aviv University, Tel Aviv 69978, Israel
| |
Collapse
|
23
|
NMR spectroscopy on domain dynamics in biomacromolecules. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2013; 112:58-117. [DOI: 10.1016/j.pbiomolbio.2013.05.001] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/03/2013] [Revised: 05/06/2013] [Accepted: 05/07/2013] [Indexed: 12/22/2022]
|
24
|
Boehr DD, Schnell JR, McElheny D, Bae SH, Duggan BM, Benkovic SJ, Dyson HJ, Wright PE. A distal mutation perturbs dynamic amino acid networks in dihydrofolate reductase. Biochemistry 2013; 52:4605-19. [PMID: 23758161 DOI: 10.1021/bi400563c] [Citation(s) in RCA: 70] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
Correlated networks of amino acids have been proposed to play a fundamental role in allostery and enzyme catalysis. These networks of amino acids can be traced from surface-exposed residues all the way into the active site, and disruption of these networks can decrease enzyme activity. Substitution of the distal Gly121 residue in Escherichia coli dihydrofolate reductase results in an up to 200-fold decrease in the hydride transfer rate despite the fact that the residue is located 15 Å from the active-site center. In this study, nuclear magnetic resonance relaxation experiments are used to demonstrate that dynamics on the picosecond to nanosecond and microsecond to millisecond time scales are changed significantly in the G121V mutant of dihydrofolate reductase. In particular, picosecond to nanosecond time scale dynamics are decreased in the FG loop (containing the mutated residue at position 121) and the neighboring active-site loop (the Met20 loop) in the mutant compared to those of the wild-type enzyme, suggesting that these loops are dynamically coupled. Changes in methyl order parameters reveal a pathway by which dynamic perturbations can be propagated more than 25 Å across the protein from the site of mutation. All of the enzyme complexes, including the model Michaelis complex with folate and nicotinamide adenine dinucleotide phosphate bound, assume an occluded ground-state conformation, and we do not observe sampling of a higher-energy closed conformation by (15)N R2 relaxation dispersion experiments. This is highly significant, because it is only in the closed conformation that the cofactor and substrate reactive centers are positioned for reaction. The mutation also impairs microsecond to millisecond time scale fluctuations that have been implicated in the release of product from the wild-type enzyme. Our results are consistent with an important role for Gly121 in controlling protein dynamics critical for enzyme function and further validate the dynamic energy landscape hypothesis of enzyme catalysis.
Collapse
Affiliation(s)
- David D Boehr
- Department of Integrative Structural and Computational Biology and the Skaggs Institute for Chemical Biology, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, CA 92037, USA
| | - Jason R Schnell
- Department of Integrative Structural and Computational Biology and the Skaggs Institute for Chemical Biology, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, CA 92037, USA
| | - Dan McElheny
- Department of Integrative Structural and Computational Biology and the Skaggs Institute for Chemical Biology, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, CA 92037, USA
| | - Sung-Hun Bae
- Department of Integrative Structural and Computational Biology and the Skaggs Institute for Chemical Biology, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, CA 92037, USA
| | - Brendan M Duggan
- Department of Integrative Structural and Computational Biology and the Skaggs Institute for Chemical Biology, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, CA 92037, USA
| | - Stephen J Benkovic
- The Pennsylvania State University, Department of Chemistry, 414 Wartik Laboratory, University Park, PA 16802, USA
| | - H Jane Dyson
- Department of Integrative Structural and Computational Biology and the Skaggs Institute for Chemical Biology, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, CA 92037, USA
| | - Peter E Wright
- Department of Integrative Structural and Computational Biology and the Skaggs Institute for Chemical Biology, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, CA 92037, USA
| |
Collapse
|
25
|
Terada TP, Kimura T, Sasai M. Entropic mechanism of allosteric communication in conformational transitions of dihydrofolate reductase. J Phys Chem B 2013; 117:12864-77. [PMID: 23705773 DOI: 10.1021/jp402071m] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
The mechanism of allosteric conformational transitions of Escherichia coli dihydrofolate reductase (DHFR) is investigated theoretically by applying a newly developed coarse-grained model. Functional forms of interaction potentials in the model depend on the local structural environments around those interactions to represent the many-residue effects due to atomic packing in each local region, and hence, this model is called "the chameleon model". The chameleon model consistently describes the free-energy landscape of two conformational transitions in the catalytic cycle of DHFR, which we call conformational transition 1 (CT1) and conformational transition 2 (CT2); CT1 is accompanied by the hydride transfer reaction, and CT2 is accompanied by the product ligand release. The transition state of CT1 is entropically stabilized by the disordering of loops at the peripheral regions of the protein, which enhances the positively correlated fluctuations at the center part of the protein, showing that the allosteric communication between distant regions through the central region is intrinsically associated with the entropic stabilization of the transition state. The transition state of CT2 is entropically stabilized through the mechanism that enhances the breathing motion of two domains, showing that the difference in the distribution of interactions brings about the difference in the transition mechanism between CT1 and CT2. The chameleon model opens a way to consistently describe the dynamical energy landscape of enzymatic reactions.
Collapse
Affiliation(s)
- Tomoki P Terada
- Department of Computational Science and Engineering and ‡Department of Applied Physics, Nagoya University , Nagoya 464-8603, Japan
| | | | | |
Collapse
|
26
|
Abstract
Allosteric propagation results in communication between distinct sites in the protein structure; it also encodes specific effects on cellular pathways, and in this way it shapes cellular response. One example of long-range effects is binding of morphogens to cell surface receptors, which initiates a cascade of protein interactions that leads to genome activation and specific cellular action. Allosteric propagation results from combinations of multiple factors, takes place through dynamic shifts of conformational ensembles, and affects the equilibria of macromolecular interactions. Here, we (a) emphasize the well-known yet still underappreciated role of allostery in conveying explicit signals across large multimolecular assemblies and distances to specify cellular action; (b) stress the need for quantitation of the allosteric effects; and finally, (c) propose that each specific combination of allosteric effectors along the pathway spells a distinct function. The challenges are colossal; the inspiring reward will be predicting function, misfunction, and outcomes of drug regimes.
Collapse
Affiliation(s)
- Ruth Nussinov
- Basic Research Program, SAIC-Frederick, Inc., Center for Cancer Research Nanobiology Program, National Cancer Institute, Frederick, Maryland 21702, USA.
| | | | | |
Collapse
|
27
|
Arora K, Brooks CL. Multiple intermediates, diverse conformations, and cooperative conformational changes underlie the catalytic hydride transfer reaction of dihydrofolate reductase. Top Curr Chem (Cham) 2013; 337:165-87. [PMID: 23420416 PMCID: PMC4394636 DOI: 10.1007/128_2012_408] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
It has become increasingly clear that protein motions play an essential role in enzyme catalysis. However, exactly how these motions are related to an enzyme's chemical step is still intensely debated. This chapter examines the possible role of protein motions that display a hierarchy of timescales in enzyme catalysis. The linkage between protein motions and catalysis is investigated in the context of a model enzyme, E. coli dihydrofolate reductase (DHFR), that catalyzes the hydride transfer reaction in the conversion of dihydrofolate to tetrahydrofolate. The results of extensive computer simulations probing the protein motions that are manifest during different steps along the turnover cycle of DHFR are summarized. Evidence is presented that the protein motions modulate the catalytic efficacy of DHFR by generating a conformational ensemble conducive to the hydride transfer. The alteration of the equilibrium conformational ensemble rather than any protein dynamical effects is found to be sufficient to explain the rate-diminishing effects of mutation on the kinetics of the enzyme. These data support the view that the protein motions facilitate catalysis by establishing reaction competent conformations of the enzyme, but they do not directly couple to the chemical reaction itself. These findings have broad implications for our understanding of enzyme mechanisms and the design of novel protein catalysts.
Collapse
Affiliation(s)
- Karunesh Arora
- Department of Chemistry and Biophysics Program, University of Michigan, Ann Arbor, MI 48109
| | - Charles L. Brooks
- Department of Chemistry and Biophysics Program, University of Michigan, Ann Arbor, MI 48109
| |
Collapse
|
28
|
Chen J, Thirumalai D. Helices 2 and 3 are the initiation sites in the PrP(C) → PrP(SC) transition. Biochemistry 2012; 52:310-9. [PMID: 23256626 DOI: 10.1021/bi3005472] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
It is established that prion protein is the sole causative agent in a number of diseases in humans and animals. However, the nature of conformational changes that the normal cellular form, PrP(C), undergoes in its conversion to a self-replicating state is still not fully understood. The ordered C-terminus of PrP(C) proteins has three helices (H1-H3). Here, we use statistical coupling analysis (SCA) to infer covariations at various locations using a family of evolutionarily related sequences and the response of mouse and human PrP(C)s to mechanical force to decipher the initiation sites for the transition from PrP(C) to an aggregation-prone PrP* state. Sequence-based SCA predicts that the clustered residues in nonmammals are localized in the stable core (near H1) of PrP(C), whereas in mammalian PrP(C), they are localized in frustrated helices H2 and H3 where most of the pathogenic mutations are found. Force-extension curves and free energy profiles as a function of extension of mouse and human PrP(C) in the absence of a disulfide (SS) bond between residues Cys179 and Cys214, generated by applying mechanical force to the ends of the molecule, show a sequence of unfolding events starting first with rupture of H2 and H3. This is followed by disruption of structure in two strands. Helix H1, stabilized by three salt bridges, resists substantial force before unfolding. Force extension profiles and the dynamics of rupture of tertiary contacts also show that even in the presence of an SS bond the instabilities in most of H3 and parts of H2 still determine the propensity to form the PrP* state. In mouse PrP(C) with an SS bond, there are ∼10 residues that retain their order even at high forces. Both SCA and single-molecule force simulations show that in the conversion from PrP(C) to PrP(SC) major conformational changes occur (at least initially) in H2 and H3, which because of their sequence compositions are frustrated in the helical state. Implications of our findings for the structural model for the scrapie form of PrP(C) are discussed.
Collapse
Affiliation(s)
- Jie Chen
- Biophysics Program, Institute for Physical Science and Technology, and ‡Department of Chemistry and Biochemistry, University of Maryland , College Park, Maryland 20742, USA
| | | |
Collapse
|
29
|
Weikl TR, Boehr DD. Conformational selection and induced changes along the catalytic cycle of Escherichia coli
dihydrofolate reductase. Proteins 2012; 80:2369-83. [DOI: 10.1002/prot.24123] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2012] [Revised: 05/02/2012] [Accepted: 05/09/2012] [Indexed: 01/21/2023]
|
30
|
Long D, Brüschweiler R. Structural and Entropic Allosteric Signal Transduction Strength via Correlated Motions. J Phys Chem Lett 2012; 3:1722-1726. [PMID: 26285736 DOI: 10.1021/jz300488e] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Allosteric signal transduction in biomacromolecules can play an essential role in their function. Internal motional correlations in proteins provide a possible communication mechanism, but the quantitative relationship between statistical correlations and allostery is unknown. Quantitative relationships between internal motional correlations and the efficiency of propagation of allosteric structural and entropic effects are introduced and validated against conformational ensembles obtained from molecular dynamics simulations. This framework can explain a range of phenomena, such as the occurrence of an allosteric entropy change in the absence of any noticeable structural change.
Collapse
Affiliation(s)
- Dong Long
- Department of Chemistry and Biochemistry and National High Magnetic Field Laboratory, Florida State University, Tallahassee, Florida 32306, United States
| | - Rafael Brüschweiler
- Department of Chemistry and Biochemistry and National High Magnetic Field Laboratory, Florida State University, Tallahassee, Florida 32306, United States
| |
Collapse
|
31
|
Hot spots for allosteric regulation on protein surfaces. Cell 2012; 147:1564-75. [PMID: 22196731 DOI: 10.1016/j.cell.2011.10.049] [Citation(s) in RCA: 263] [Impact Index Per Article: 21.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2011] [Revised: 08/10/2011] [Accepted: 10/19/2011] [Indexed: 11/22/2022]
Abstract
Recent work indicates a general architecture for proteins in which sparse networks of physically contiguous and coevolving amino acids underlie basic aspects of structure and function. These networks, termed sectors, are spatially organized such that active sites are linked to many surface sites distributed throughout the structure. Using the metabolic enzyme dihydrofolate reductase as a model system, we show that: (1) the sector is strongly correlated to a network of residues undergoing millisecond conformational fluctuations associated with enzyme catalysis, and (2) sector-connected surface sites are statistically preferred locations for the emergence of allosteric control in vivo. Thus, sectors represent an evolutionarily conserved "wiring" mechanism that can enable perturbations at specific surface positions to rapidly initiate conformational control over protein function. These findings suggest that sectors enable the evolution of intermolecular communication and regulation.
Collapse
|
32
|
Lee Y, Jeong LS, Choi S, Hyeon C. Link between allosteric signal transduction and functional dynamics in a multisubunit enzyme: S-adenosylhomocysteine hydrolase. J Am Chem Soc 2011; 133:19807-15. [PMID: 22023331 DOI: 10.1021/ja2066175] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
S-adenosylhomocysteine hydrolase (SAHH), a cellular enzyme that plays a key role in methylation reactions including those required for maturation of viral mRNA, is an important drug target in the discovery of antiviral agents. While targeting the active site is a straightforward strategy of enzyme inhibition, evidence of allosteric modulation of active site in many enzymes underscores the molecular origin of signal transduction. Information of co-evolving sequences in SAHH family and the key residues for functional dynamics that can be identified using native topology of the enzyme provide glimpses into how the allosteric signaling network, dispersed over the molecular structure, coordinates intra- and intersubunit conformational dynamics. To study the link between the allosteric communication and functional dynamics of SAHHs, we performed Brownian dynamics simulations by building a coarse-grained model based on the holo and ligand-bound structures. The simulations of ligand-induced transition revealed that the signal of intrasubunit closure dynamics is transmitted to form intersubunit contacts, which in turn invoke a precise alignment of active site, followed by the dimer-dimer rotation that compacts the whole tetrameric structure. Further analyses of SAHH dynamics associated with ligand binding provided evidence of both induced fit and population shift mechanisms and also showed that the transition-state ensemble is akin to the ligand-bound state. Besides the formation of enzyme-ligand contacts at the active site, the allosteric couplings from the residues distal to the active site are vital to the enzymatic function.
Collapse
Affiliation(s)
- Yoonji Lee
- College of Pharmacy, Division of Life and Pharmaceutical Sciences and National Core Research Center for Cell Signaling and Drug Discovery Research, Ewha Womans University, Seoul 120-750, Republic of Korea
| | | | | | | |
Collapse
|
33
|
Dixit A, Verkhivker GM. Computational modeling of allosteric communication reveals organizing principles of mutation-induced signaling in ABL and EGFR kinases. PLoS Comput Biol 2011; 7:e1002179. [PMID: 21998569 PMCID: PMC3188506 DOI: 10.1371/journal.pcbi.1002179] [Citation(s) in RCA: 63] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2011] [Accepted: 07/16/2011] [Indexed: 12/15/2022] Open
Abstract
The emerging structural information about allosteric kinase complexes and the growing number of allosteric inhibitors call for a systematic strategy to delineate and classify mechanisms of allosteric regulation and long-range communication that control kinase activity. In this work, we have investigated mechanistic aspects of long-range communications in ABL and EGFR kinases based on the results of multiscale simulations of regulatory complexes and computational modeling of signal propagation in proteins. These approaches have been systematically employed to elucidate organizing molecular principles of allosteric signaling in the ABL and EGFR multi-domain regulatory complexes and analyze allosteric signatures of the gate-keeper cancer mutations. We have presented evidence that mechanisms of allosteric activation may have universally evolved in the ABL and EGFR regulatory complexes as a product of a functional cross-talk between the organizing αF-helix and conformationally adaptive αI-helix and αC-helix. These structural elements form a dynamic network of efficiently communicated clusters that may control the long-range interdomain coupling and allosteric activation. The results of this study have unveiled a unifying effect of the gate-keeper cancer mutations as catalysts of kinase activation, leading to the enhanced long-range communication among allosterically coupled segments and stabilization of the active kinase form. The results of this study can reconcile recent experimental studies of allosteric inhibition and long-range cooperativity between binding sites in protein kinases. The presented study offers a novel molecular insight into mechanistic aspects of allosteric kinase signaling and provides a quantitative picture of activation mechanisms in protein kinases at the atomic level. Despite recent progress in computational and experimental studies of dynamic regulation in protein kinases, a mechanistic understanding of long-range communication and mechanisms of mutation-induced signaling controlling kinase activity remains largely qualitative. In this study, we have performed a systematic modeling and analysis of allosteric activation in ABL and EGFR kinases at the increasing level of complexity - from catalytic domain to multi-domain regulatory complexes. The results of this study have revealed organizing structural and mechanistic principles of allosteric signaling in protein kinases. Although activation mechanisms in ABL and EGFR kinases have evolved through acquisition of structurally different regulatory complexes, we have found that long-range interdomain communication between common functional segments (αF-helix and αC-helix) may be important for allosteric activation. The results of study have revealed molecular signatures of activating cancer mutations and have shed the light on general mechanistic aspects of mutation-induced signaling in protein kinases. An advanced understanding and further characterization of molecular signatures of kinase mutations may aid in a better rationalization of mutational effects on clinical outcomes and facilitate molecular-based therapeutic strategies to combat kinase mutation-dependent tumorigenesis.
Collapse
Affiliation(s)
- Anshuman Dixit
- Department of Pharmaceutical Chemistry, School of Pharmacy, The University of Kansas, Lawrence, Kansas, United States of America
| | - Gennady M. Verkhivker
- Department of Pharmaceutical Chemistry, School of Pharmacy, The University of Kansas, Lawrence, Kansas, United States of America
- Department of Pharmacology, University of California San Diego, La Jolla, California, United States of America
- * E-mail:
| |
Collapse
|
34
|
Dixit A, Verkhivker GM. The energy landscape analysis of cancer mutations in protein kinases. PLoS One 2011; 6:e26071. [PMID: 21998754 PMCID: PMC3188581 DOI: 10.1371/journal.pone.0026071] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2011] [Accepted: 09/19/2011] [Indexed: 11/18/2022] Open
Abstract
The growing interest in quantifying the molecular basis of protein kinase activation and allosteric regulation by cancer mutations has fueled computational studies of allosteric signaling in protein kinases. In the present study, we combined computer simulations and the energy landscape analysis of protein kinases to characterize the interplay between oncogenic mutations and locally frustrated sites as important catalysts of allostetric kinase activation. While structurally rigid kinase core constitutes a minimally frustrated hub of the catalytic domain, locally frustrated residue clusters, whose interaction networks are not energetically optimized, are prone to dynamic modulation and could enable allosteric conformational transitions. The results of this study have shown that the energy landscape effect of oncogenic mutations may be allosteric eliciting global changes in the spatial distribution of highly frustrated residues. We have found that mutation-induced allosteric signaling may involve a dynamic coupling between structurally rigid (minimally frustrated) and plastic (locally frustrated) clusters of residues. The presented study has demonstrated that activation cancer mutations may affect the thermodynamic equilibrium between kinase states by allosterically altering the distribution of locally frustrated sites and increasing the local frustration in the inactive form, while eliminating locally frustrated sites and restoring structural rigidity of the active form. The energy landsape analysis of protein kinases and the proposed role of locally frustrated sites in activation mechanisms may have useful implications for bioinformatics-based screening and detection of functional sites critical for allosteric regulation in complex biomolecular systems.
Collapse
Affiliation(s)
- Anshuman Dixit
- Department of Pharmaceutical Chemistry, School of Pharmacy, The University of Kansas, Lawrence, Kansas, United States of America
| | - Gennady M. Verkhivker
- School of Computational Sciences and Crean School of Health and Life Sciences, Schmid College of Science and Technology, Chapman University, Orange, California, United States of America
- Department of Pharmacology, University of California, San Diego, La Jolla, California, United States of America
- * E-mail:
| |
Collapse
|
35
|
Duan L, Zhmurov A, Barsegov V, Dima RI. Exploring the mechanical stability of the C2 domains in human synaptotagmin 1. J Phys Chem B 2011; 115:10133-46. [PMID: 21776988 DOI: 10.1021/jp2025945] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Human synaptotagmin 1 (Syt1) plays a crucial role in the bending of the membrane during neurotransmitter release at the synapse. Hence, resolving the structural details of Syt1 that underlie its biological function is fundamental for providing mechanistic insights into the nature of the synaptic response. We explored the unfolding micromechanics of Syt1 by analyzing the free energy landscape of the whole molecule and its C2A and C2B domains. We employed a self-organized polymer (SOP) model of a protein chain to carry out pulling simulations, accelerated on graphics processing units (GPUs), under experimental force loads. To resolve the atomic-level details, we complemented the SOP model simulations with atomistic simulations. On the basis of the results obtained, we hypothesize that (1) isolated single domains C2A and C2B present similar mechanical resistance against an applied pulling force but unfold following different kinetic pathways and that (2) C2B is more mechanically resistant in the C2AB complex due to stabilizing interactions with other domains. These findings correlate well with recent atomic force microscopy (AFM) studies on the Syt1 molecule, in which the increase in the unfolding force for C2B was detected when this domain was joined with C2A. Our results also suggest that the linkers (I27 domains) used in the experimental setup can modulate the mechanical behavior of this synaptic protein complex and alter not only the critical force for unfolding but also the unfolding pathways for the C2 domains. Interestingly, the presence of the C2A-C2B domain interface in the C2AB complex confers mechanical stability to either of the C2 domains. Our findings provide new insights into the relative conformational variability of the C2 domains, which we believe to be modulated, to a large extent, by intermolecular coupling with other proteins.
Collapse
Affiliation(s)
- Li Duan
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221, USA
| | | | | | | |
Collapse
|
36
|
Kamerlin SCL, Vicatos S, Dryga A, Warshel A. Coarse-grained (multiscale) simulations in studies of biophysical and chemical systems. Annu Rev Phys Chem 2011; 62:41-64. [PMID: 21034218 DOI: 10.1146/annurev-physchem-032210-103335] [Citation(s) in RCA: 145] [Impact Index Per Article: 11.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Recent years have witnessed an explosion in computational power, leading to attempts to model ever more complex systems. Nevertheless, there remain cases for which the use of brute-force computer simulations is clearly not the solution. In such cases, great benefit can be obtained from the use of physically sound simplifications. The introduction of such coarse graining can be traced back to the early usage of a simplified model in studies of proteins. Since then, the field has progressed tremendously. In this review, we cover both key developments in the field and potential future directions. Additionally, particular emphasis is given to two general approaches, namely the renormalization and reference potential approaches, which allow one to move back and forth between the coarse-grained (CG) and full models, as these approaches provide the foundation for CG modeling of complex systems.
Collapse
Affiliation(s)
- Shina C L Kamerlin
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | | | | | |
Collapse
|
37
|
Bhattacharyya M, Vishveshwara S. Probing the allosteric mechanism in pyrrolysyl-tRNA synthetase using energy-weighted network formalism. Biochemistry 2011; 50:6225-36. [PMID: 21650159 DOI: 10.1021/bi200306u] [Citation(s) in RCA: 45] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Pyrrolysyl-tRNA synthetase (PylRS) is an atypical enzyme responsible for charging tRNA(Pyl) with pyrrolysine, despite lacking precise tRNA anticodon recognition. This dimeric protein exhibits allosteric regulation of function, like any other tRNA synthetases. In this study we examine the paths of allosteric communication at the atomic level, through energy-weighted networks of Desulfitobacterium hafniense PylRS (DhPylRS) and its complexes with tRNA(Pyl) and activated pyrrolysine. We performed molecular dynamics simulations of the structures of these complexes to obtain an ensemble conformation-population perspective. Weighted graph parameters relevant to identifying key players and ties in the context of social networks such as edge/node betweenness, closeness index, and the concept of funneling are explored in identifying key residues and interactions leading to shortest paths of communication in the structure networks of DhPylRS. Further, the changes in the status of important residues and connections and the costs of communication due to ligand induced perturbations are evaluated. The optimal, suboptimal, and preexisting paths are also investigated. Many of these parameters have exhibited an enhanced asymmetry between the two subunits of the dimeric protein, especially in the pretransfer complex, leading us to conclude that encoding of function goes beyond the sequence/structure of proteins. The local and global perturbations mediated by appropriate ligands and their influence on the equilibrium ensemble of conformations also have a significant role to play in the functioning of proteins. Taking a comprehensive view of these observations, we propose that the origin of many functional aspects (allostery and half-sites reactivity in the case of DhPylRS) lies in subtle rearrangements of interactions and dynamics at a global level.
Collapse
|
38
|
Armenta-Medina D, Pérez-Rueda E, Segovia L. Identification of functional motions in the adenylate kinase (ADK) protein family by computational hybrid approaches. Proteins 2011; 79:1662-71. [PMID: 21365689 DOI: 10.1002/prot.22995] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2010] [Revised: 12/03/2010] [Accepted: 12/07/2010] [Indexed: 02/02/2023]
Abstract
Based on integrative computational hybrid approaches that combined statistical coupling analysis (SCA), molecular dynamics (MD), and normal mode analysis (NMA), evolutionarily coupled residues involved in functionally relevant motion in the adenylate kinase protein family were identified. The hybrids identified four top-ranking site pairs that belong to a conserved hydrogen bond network that is involved in the enzyme's flexibility. A second group of top-ranking site pairs was identified in critical regions for functional dynamics, such as those related to enzymatic turnover. The high consistency of the results obtained by SCA with NMA (SCA.NMA) and by SCA.MD hybrid analyses suggests that suitable replacement of the matrix of cross-correlation analysis of atomic fluctuations (derived by using NMA) with those based on MD contributes to the identification of such sites by means of a fast computational calculation. The analysis presented here strongly supports the hypothesis that evolutionary forces, such as coevolution at the sequence level, have promoted functional dynamic properties of the adenylate kinase protein family. Finally, these hybrid approaches can be used to identify, at the residue level, protein motion coordination patterns not previously observed, such as in hinge regions.
Collapse
Affiliation(s)
- Dagoberto Armenta-Medina
- Departamento de Ingeniería Celular y Biocatálisis, Instituto de Biotecnología, Universidad Nacional Autónoma de México, Cuernavaca, Morelos, México.
| | | | | |
Collapse
|
39
|
Elber R. Simulations of allosteric transitions. Curr Opin Struct Biol 2011; 21:167-72. [PMID: 21333527 DOI: 10.1016/j.sbi.2011.01.012] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2010] [Revised: 01/25/2011] [Accepted: 01/27/2011] [Indexed: 11/30/2022]
Abstract
Allosteric transitions are one of the subtler mechanisms used by nature to fine tune protein activity. Effector binding to a specific site on the protein surface induces significant activity change, and initiates a conformational transition that frequently includes domain motions and is very large. From a theoretical and biophysical perspective two problems are particularly intriguing. The first is the way in which a launching signal, which is spatially confined and includes only a few interacting atoms, is propagated to a large-scale conformational transition we frequently see in allosteric transitions. Hence, there is the question of how a small perturbation is magnified to yield motions of thousands of atoms. The second puzzle is of focus, coherence, and efficiency. The impact of the binding of the effector is sometimes extended over tens of angstroms. How the signal is transmitted and kept significant over such large distances in the 'noisy' molecular environment is another major direction of investigation. In the present review we examined different theoretical and computational attempts to solve the questions.
Collapse
Affiliation(s)
- Ron Elber
- Department of Chemistry and Biochemistry, Institute of Computational Engineering and Sciences, 1 University Station, ICES, C0200, The University of Texas at Austin, Austin, TX 78712, USA.
| |
Collapse
|
40
|
Brut M, Estève A, Landa G, Renvez G, Djafari Rouhani M, Vaisset M. Atomic Scale Determination of Enzyme Flexibility and Active Site Stability through Static Modes: Case of Dihydrofolate Reductase. J Phys Chem B 2011; 115:1616-22. [DOI: 10.1021/jp109874z] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Marie Brut
- LAAS-CNRS, Université de Toulouse, 7 avenue du Colonel Roche, F-31077 Toulouse, France
- Université de Toulouse, UPS, INSA, INP, ISAE, LAAS, F-31077 Toulouse, France
- Department of Structural Biology, Stanford University School of Medicine, Stanford, California, United States
| | - Alain Estève
- LAAS-CNRS, Université de Toulouse, 7 avenue du Colonel Roche, F-31077 Toulouse, France
- Université de Toulouse, UPS, INSA, INP, ISAE, LAAS, F-31077 Toulouse, France
| | - Georges Landa
- LAAS-CNRS, Université de Toulouse, 7 avenue du Colonel Roche, F-31077 Toulouse, France
- Université de Toulouse, UPS, INSA, INP, ISAE, LAAS, F-31077 Toulouse, France
| | - Guillaume Renvez
- LAAS-CNRS, Université de Toulouse, 7 avenue du Colonel Roche, F-31077 Toulouse, France
- Université de Toulouse, UPS, INSA, INP, ISAE, LAAS, F-31077 Toulouse, France
| | - Mehdi Djafari Rouhani
- LAAS-CNRS, Université de Toulouse, 7 avenue du Colonel Roche, F-31077 Toulouse, France
- Université de Toulouse, UPS, INSA, INP, ISAE, LAAS, F-31077 Toulouse, France
| | - Marc Vaisset
- LAAS-CNRS, Université de Toulouse, 7 avenue du Colonel Roche, F-31077 Toulouse, France
- Université de Toulouse, UPS, INSA, INP, ISAE, LAAS, F-31077 Toulouse, France
| |
Collapse
|
41
|
Zhmurov A, Dima RI, Kholodov Y, Barsegov V. Sop-GPU: accelerating biomolecular simulations in the centisecond timescale using graphics processors. Proteins 2011; 78:2984-99. [PMID: 20715052 DOI: 10.1002/prot.22824] [Citation(s) in RCA: 55] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Theoretical exploration of fundamental biological processes involving the forced unraveling of multimeric proteins, the sliding motion in protein fibers and the mechanical deformation of biomolecular assemblies under physiological force loads is challenging even for distributed computing systems. Using a C(α)-based coarse-grained self organized polymer (SOP) model, we implemented the Langevin simulations of proteins on graphics processing units (SOP-GPU program). We assessed the computational performance of an end-to-end application of the program, where all the steps of the algorithm are running on a GPU, by profiling the simulation time and memory usage for a number of test systems. The ∼90-fold computational speedup on a GPU, compared with an optimized central processing unit program, enabled us to follow the dynamics in the centisecond timescale, and to obtain the force-extension profiles using experimental pulling speeds (v(f) = 1-10 μm/s) employed in atomic force microscopy and in optical tweezers-based dynamic force spectroscopy. We found that the mechanical molecular response critically depends on the conditions of force application and that the kinetics and pathways for unfolding change drastically even upon a modest 10-fold increase in v(f). This implies that, to resolve accurately the free energy landscape and to relate the results of single-molecule experiments in vitro and in silico, molecular simulations should be carried out under the experimentally relevant force loads. This can be accomplished in reasonable wall-clock time for biomolecules of size as large as 10(5) residues using the SOP-GPU package.
Collapse
Affiliation(s)
- A Zhmurov
- Department of Chemistry, University of Massachusetts, Lowell, Massachusetts 01854, USA
| | | | | | | |
Collapse
|
42
|
Protein stability, flexibility and function. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2010; 1814:969-76. [PMID: 21094283 DOI: 10.1016/j.bbapap.2010.11.005] [Citation(s) in RCA: 161] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/29/2010] [Revised: 09/24/2010] [Accepted: 11/11/2010] [Indexed: 02/08/2023]
Abstract
Proteins rely on flexibility to respond to environmental changes, ligand binding and chemical modifications. Potentially, a perturbation that changes the flexibility of a protein may interfere with its function. Millions of mutations have been performed on thousands of proteins in quests for a delineation of the molecular details of their function. Several of these mutations interfered with the binding of a specific ligand with a concomitant effect on the stability of the protein scaffold. It has been ambiguous and not straightforward to recognize if any relationships exist between the stability of a protein and the affinity for its ligand. In this review, we present examples of proteins where changes in stability results in changes in affinity and of proteins where stability and affinity are uncorrelated. We discuss the possibility for a relationship between stability and binding. From the data presented is it clear that there are specific sites (flexibility hotspots) in proteins that are important for both binding and stability. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.
Collapse
|
43
|
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
|
44
|
Lambry JC, Beaumont E, Tarus B, Blanchard-Desce M, Slama-Schwok A. Selective probing of a NADPH site controlled light-induced enzymatic catalysis. J Mol Recognit 2010; 23:379-88. [PMID: 20029835 DOI: 10.1002/jmr.1009] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
Achieving molecular recognition of NADPH binding sites is a compelling strategy to control many redox biological processes. The NADPH sites recognize the ubiquitous NADPH cofactor via highly conserved binding interactions, despite differences in the regulation of the hydride transfer in redox active proteins. We recently developed a photoactive NADPH substitute, called nanotrigger NT synchronizing the initiation of enzymatic catalysis of the endothelial NO-synthase (eNOS) with a laser pulse. Spatial and temporal control of enzymatic activity by such a designed light-driven activator would benefit from achieving molecular selectivity, i.e. activation of a single NADPH-mediated enzyme.In this work, we probe the ability of NT to discriminate between two NADPH sites with light. The selected NADPH sites belong to dihydrofolate reductase dihydrofolate reductase enzyme (DHFR) and endothelial NO-synthase (eNOS). Ultrafast kinetics showed that NT could not activate DHFR catalysis with a laser pulse in contrast with the observed trigger of eNOS catalysis leading to NO formation. Homology modelling, molecular dynamics simulations showed that NT discriminated between the two NADPH sites by different donor to acceptor distances and by local steric effects hindering light activation of DHFR catalysis. The data suggested that the narrow NADPH site required a tight fit of the nanotrigger at a suitable distance/angle to the electron acceptor for a specific activation of the catalysis. The ability of the nanotrigger to activate eNOS combined with a low reactivity in unfavourable NADPH sites makes NT a highly promising tool for targeting eNOS in endothelial cells with a laser pulse.
Collapse
Affiliation(s)
- Jean-Christophe Lambry
- Unité INSERM 696, Laboratory for Optics & Biosciences, Ecole Polytechnique, Palaiseau, France
| | | | | | | | | |
Collapse
|
45
|
Tehver R, Thirumalai D. Rigor to post-rigor transition in myosin V: link between the dynamics and the supporting architecture. Structure 2010; 18:471-81. [PMID: 20399184 DOI: 10.1016/j.str.2010.01.019] [Citation(s) in RCA: 47] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2009] [Revised: 01/01/2010] [Accepted: 01/05/2010] [Indexed: 10/19/2022]
Abstract
The detachment kinetics from actin upon ATP binding is a key step in the reaction cycle of myosin V. We show that a network of residues, constituting the allostery wiring diagram (AWD), that trigger the rigor (R) to post-rigor (PR) transition, span key structural elements from the ATP and actin-binding regions. Several of the residues are in the 33 residue helix (H18), P loop, and switch I. Brownian dynamics simulations show that a hierarchy of kinetically controlled local structural changes leads to the opening of the "cleft" region, resulting in the detachment of the motor domain from actin. Movements in switch I and P loop facilitate changes in the rest of the motor domain, in particular the rotation of H18, whose stiffness within the motor domain is crucial in the R --> PR transition. The finding that residues in the AWD also drive the kinetics of the R --> PR transition shows how the myosin architecture regulates the allosteric movements during the reaction cycle.
Collapse
Affiliation(s)
- Riina Tehver
- Biophysics Program, Institute for Physical Science and Technology, University of Maryland, College Park, MD 20742, USA
| | | |
Collapse
|
46
|
Brüschweiler S, Schanda P, Kloiber K, Brutscher B, Kontaxis G, Konrat R, Tollinger M. Direct observation of the dynamic process underlying allosteric signal transmission. J Am Chem Soc 2010; 131:3063-8. [PMID: 19203263 DOI: 10.1021/ja809947w] [Citation(s) in RCA: 102] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Allosteric regulation is an effective mechanism of control in biological processes. In allosteric proteins a signal originating at one site in the molecule is communicated through the protein structure to trigger a specific response at a remote site. Using NMR relaxation dispersion techniques we directly observe the dynamic process through which the KIX domain of CREB binding protein communicates allosteric information between binding sites. KIX mediates cooperativity between pairs of transcription factors through binding to two distinct interaction surfaces in an allosteric manner. We show that binding the activation domain of the mixed lineage leukemia (MLL) transcription factor to KIX induces a redistribution of the relative populations of KIX conformations toward a high-energy state in which the allosterically activated second binding site is already preformed, consistent with the Monod-Wyman-Changeux (WMC) model of allostery. The structural rearrangement process that links the two conformers and by which allosteric information is communicated occurs with a time constant of 3 ms at 27 degrees C. Our dynamic NMR data reveal that an evolutionarily conserved network of hydrophobic amino acids constitutes the pathway through which information is transmitted.
Collapse
Affiliation(s)
- Sven Brüschweiler
- Department of Computational & Structural Biology, Max F. Perutz Laboratories, Campus Vienna Biocenter 5, A-1030 Vienna, Austria
| | | | | | | | | | | | | |
Collapse
|
47
|
Promoter melting triggered by bacterial RNA polymerase occurs in three steps. Proc Natl Acad Sci U S A 2010; 107:12523-8. [PMID: 20615963 DOI: 10.1073/pnas.1003533107] [Citation(s) in RCA: 51] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Abstract
RNA synthesis, carried out by DNA-dependent RNA polymerase (RNAP) in a process called transcription, involves several stages. In bacteria, transcription initiation starts with promoter recognition and binding of RNAP holoenzyme, resulting in the formation of the closed (R.P(c)) RNAP-promoter DNA complex. Subsequently, a transition to the open R.P(o) complex occurs, characterized by separation of the promoter DNA strands in an approximately 12 base-pair region to form the transcription bubble. Using coarse-grained self-organized polymer models of Thermus aquatics RNAP holoenzyme and promoter DNA complexes, we performed Brownian dynamics simulations of the R.P(c) --> R.P(o) transition. In the fast trajectories, unwinding of the promoter DNA begins by local melting around the -10 element, which is followed by sequential unzipping of DNA till the +2 site. The R.P(c) --> R.P(o) transition occurs in three steps. In step I, dsDNA melts and the nontemplate strand makes stable interactions with RNAP. In step II, DNA scrunches into RNA polymerase and the downstream base pairs sequentially open to form the transcription bubble, which results in strain build up. Subsequently, downstream dsDNA bending relieves the strain as R.P(o) forms. Entry of the dsDNA into the active-site channel of RNAP requires widening of the channel, which occurs by a swing mechanism involving transient movements of a subdomain of the beta subunit caused by steric repulsion with the DNA template strand. If premature local melting away from the -10 element occurs first then the transcription bubble formation is slow involving reformation of the opened base pairs and subsequent sequential unzipping as in the fast trajectories.
Collapse
|
48
|
Joshi H, Momin F, Haines KE, Dima RI. Exploring the contribution of collective motions to the dynamics of forced-unfolding in tubulin. Biophys J 2010; 98:657-66. [PMID: 20159162 DOI: 10.1016/j.bpj.2009.10.043] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2009] [Revised: 10/05/2009] [Accepted: 10/27/2009] [Indexed: 12/28/2022] Open
Abstract
Decomposition of the intrinsic dynamics of proteins into collective motions among distant regions of the protein structure provides a physically appealing approach that couples the dynamics of the system with its functional role. The cellular functions of microtubules (an essential component of the cytoskeleton in all eukaryotic cells) depend on their dynamic instability, which is altered by various factors among which applied forces are central. To shed light on the coupling between forces and the dynamic instability of microtubules, we focus on the investigation of the response of the microtubule subunits (tubulin) to applied forces. We address this point by adapting an approach designed to survey correlations for the equilibrium dynamics of proteins to the case of correlations for proteins forced-dynamics. The resulting collective motions in tubulin have a number of functional implications, such as the identification of long-range couplings with a role in blocking the dynamic instability of microtubules. A fundamental implication of our study for the life of a cell is that, to increase the likelihood of unraveling of large cytoskeletal filaments under physiological forces, molecular motors must use a combination of pulling and torsion rather than just pulling.
Collapse
Affiliation(s)
- Harshad Joshi
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio, USA
| | | | | | | |
Collapse
|
49
|
Mauldin RV, Lee AL. Nuclear magnetic resonance study of the role of M42 in the solution dynamics of Escherichia coli dihydrofolate reductase. Biochemistry 2010; 49:1606-15. [PMID: 20073522 DOI: 10.1021/bi901798g] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
It is widely recognized that key positions throughout a protein's structure contribute unequally to function. In light of recent studies that suggest protein dynamics are required for function, a number of these residues may serve to promote motions required for ligand binding and catalysis. In this nuclear magnetic resonance (NMR) study, the conformational dynamics of the dihydrofolate reductase (DHFR) mutant M42W, in the presence of methotrexate and NADPH, are characterized and compared to those of the wild-type enzyme. M42 is distal to the active site, yet the M42W substitution regulates catalysis and ligand affinity and is therefore analogous to an allosteric modulator of DHFR function. To gain understanding of how this mutation regulates activity, we employ a "pandynamic" strategy by measuring conformational fluctuations of backbone amide and side-chain methyl groups on multiple time scales. Changes in pico- to nanosecond dynamics indicate that the mutational effects are propagated throughout a network of interacting residues within DHFR, consistent with a role for M42 as a dynamic communication hub. On the micro- to millisecond time scale, mutation increases the rate of switching in the catalytic core. Mutation also introduces switching in the adenosine binding subdomain that occurs at a higher frequency than in the catalytic core and which correlates with the rate of product release for M42W-DHFR. Finally, a structurally inferred analysis of side-chain dynamics suggests that the M42W mutation dampens motional contributions from nonlocal sources. These data show that the M42W mutation alters the dynamics of DHFR and are consistent with theoretical analysis that suggests this mutation disrupts motion that promotes catalysis.
Collapse
Affiliation(s)
- Randall V Mauldin
- Department of Biochemistry and Biophysics, School of Medicine, University of North Carolina, Chapel Hill, North Carolina 27599-7568, USA
| | | |
Collapse
|
50
|
Millisecond timescale fluctuations in dihydrofolate reductase are exquisitely sensitive to the bound ligands. Proc Natl Acad Sci U S A 2010; 107:1373-8. [PMID: 20080605 DOI: 10.1073/pnas.0914163107] [Citation(s) in RCA: 118] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Enzyme catalysis can be described as progress over a multi-dimensional energy landscape where ensembles of interconverting conformational substates channel the enzyme through its catalytic cycle. We applied NMR relaxation dispersion to investigate the role of bound ligands in modulating the dynamics and energy landscape of Escherichia coli dihydrofolate reductase to obtain insights into the mechanism by which the enzyme efficiently samples functional conformations as it traverses its reaction pathway. Although the structural differences between the occluded substrate binary complexes and product ternary complexes are very small, there are substantial differences in protein dynamics. Backbone fluctuations on the micros-ms timescale in the cofactor binding cleft are similar for the substrate and product binary complexes, but fluctuations on this timescale in the active site loops are observed only for complexes with substrate or substrate analog and are not observed for the binary product complex. The dynamics in the substrate and product binary complexes are governed by quite different kinetic and thermodynamic parameters. Analogous dynamic differences in the E:THF:NADPH and E:THF:NADP(+) product ternary complexes are difficult to rationalize from ground-state structures. For both of these complexes, the nicotinamide ring resides outside the active site pocket in the ground state. However, they differ in the structure, energetics, and dynamics of accessible higher energy substates where the nicotinamide ring transiently occupies the active site. Overall, our results suggest that dynamics in dihydrofolate reductase are exquisitely "tuned" for every intermediate in the catalytic cycle; structural fluctuations efficiently channel the enzyme through functionally relevant conformational space.
Collapse
|