1
|
Genchev GZ, Kobayashi M, Kobayashi T, Lu H. Molecular dynamics provides new insights into the mechanism of calcium signal transduction and interdomain interactions in cardiac troponin. FEBS Open Bio 2021; 11:1841-1853. [PMID: 33085832 PMCID: PMC8255835 DOI: 10.1002/2211-5463.13009] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2020] [Revised: 10/05/2020] [Accepted: 10/17/2020] [Indexed: 12/16/2022] Open
Abstract
Understanding the regulation of cardiac muscle contraction at a molecular level is crucial for the development of therapeutics for heart conditions. Despite the availability of atomic structures of the protein components of cardiac muscle thin filaments, detailed insights into their dynamics and response to calcium are yet to be fully depicted. In this study, we used molecular dynamics simulations of the core domains of the cardiac muscle protein troponin to characterize the equilibrium dynamics of its calcium-bound and calcium-free forms, with a focus on elements of cardiac muscle contraction activation and deactivation, that is, calcium binding to the cardiac troponin Ca2+ -binding subunit (TnC) and the release of the switch region of the troponin inhibitory subunit (TnI) from TnC. The process of calcium binding to the TnC binding site is described as a three-step process commencing with calcium capture by the binding site residues, followed by cooperative residue interplay bringing the calcium ion to the binding site, and finally, calcium-water exchange. Furthermore, we uncovered a set of TnC-TnI interdomain interactions that are critical for TnC N-lobe hydrophobic pocket dynamics. Absence of these interactions allows the closure of the TnC N-lobe hydrophobic pocket while the TnI switch region remains expelled, whereas if the interactions are maintained, the hydrophobic pocket remains open. Modification of these interactions may fine-tune the ability of the TnC N-lobe hydrophobic pocket to close or remain open, modulate cardiac contractility and present potential therapy-relevant targets.
Collapse
Affiliation(s)
- Georgi Z Genchev
- Center for Biomedical Informatics, Shanghai Children's Hospital, Shanghai, China.,SJTU-Yale Joint Center for Biostatistics, Shanghai Jiao Tong University, Shanghai, China.,Bulgarian Institute for Genomics and Precision Medicine, Sofia, Bulgaria.,Bioinformatics Program, Department of Bioengineering, University of Illinois at Chicago, Chicago, IL, USA
| | - Minae Kobayashi
- Department of Physiology and Biophysics and Center for Cardiovascular Research, University of Illinois at Chicago College of Medicine, Chicago, IL, USA
| | - Tomoyoshi Kobayashi
- Department of Physiology and Biophysics and Center for Cardiovascular Research, University of Illinois at Chicago College of Medicine, Chicago, IL, USA
| | - Hui Lu
- Center for Biomedical Informatics, Shanghai Children's Hospital, Shanghai, China.,SJTU-Yale Joint Center for Biostatistics, Shanghai Jiao Tong University, Shanghai, China.,Department of Bioinformatics and Biostatistics, Shanghai Jiao Tong University, Shanghai, China
| |
Collapse
|
2
|
Aboelkassem Y, McCabe KJ, Huber GA, Regnier M, McCammon JA, McCulloch AD. A Stochastic Multiscale Model of Cardiac Thin Filament Activation Using Brownian-Langevin Dynamics. Biophys J 2019; 117:2255-2272. [PMID: 31547973 PMCID: PMC6990154 DOI: 10.1016/j.bpj.2019.08.003] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2019] [Revised: 07/31/2019] [Accepted: 08/02/2019] [Indexed: 11/16/2022] Open
Abstract
We use Brownian-Langevin dynamics principles to derive a coarse-graining multiscale myofilament model that can describe the thin-filament activation process during contraction. The model links atomistic molecular simulations of protein-protein interactions in the thin-filament regulatory unit to sarcomere-level activation dynamics. We first calculate the molecular interaction energy between tropomyosin and actin surface using Brownian dynamics simulations. This energy profile is then generalized to account for the observed tropomyosin transitions between its regulatory stable states. The generalized energy landscape then served as a basis for developing a filament-scale model using Langevin dynamics. This integrated analysis, spanning molecular to thin-filament scales, is capable of tracking the events of the tropomyosin conformational changes as it moves over the actin surface. The tropomyosin coil with flexible overlap regions between adjacent tropomyosins is represented in the model as a system of coupled stochastic ordinary differential equations. The proposed multiscale approach provides a more detailed molecular connection between tropomyosin dynamics, the trompomyosin-actin interaction-energy landscape, and the generated force by the sarcomere.
Collapse
Affiliation(s)
- Yasser Aboelkassem
- Department of Bioengineering, University of California San Diego, La Jolla, California.
| | - Kimberly J McCabe
- Department of Bioengineering, University of California San Diego, La Jolla, California
| | - Gary A Huber
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California
| | - Michael Regnier
- Department of Bioengineering, University of Washington, Seattle, Washington
| | - J Andrew McCammon
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California
| | - Andrew D McCulloch
- Department of Bioengineering, University of California San Diego, La Jolla, California
| |
Collapse
|
3
|
Bowman JD, Lindert S. Computational Studies of Cardiac and Skeletal Troponin. Front Mol Biosci 2019; 6:68. [PMID: 31448287 PMCID: PMC6696891 DOI: 10.3389/fmolb.2019.00068] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2019] [Accepted: 07/25/2019] [Indexed: 12/16/2022] Open
Abstract
Troponin is a key regulatory protein in muscle contraction, consisting of three subunits troponin C (TnC), troponin I (TnI), and troponin T (TnT). Calcium association to TnC initiates contraction by causing a series of dynamic and conformational changes that allow the switch peptide of TnI to bind and subsequently cross bridges to form between the thin and thick filament of the sarcomere. Owing to its pivotal role in contraction regulation, troponin has been the focus of numerous computational studies over the last decade. These studies elegantly supplemented a large volume of experimental work and focused on the structure, dynamics and function of the whole troponin complex, individual subunits, and even on segments of the thin filament. Molecular dynamics, Brownian dynamics, and free energy simulations have been used to elucidate the conformational dynamics and underlying free energy landscape of troponin, calcium, and switch peptide binding, as well as the effect of disease mutations, small molecules and post-translational modifications such as phosphorylation. Frequently, simulations have been used to confirm or explain experimental observations. Computer-aided drug discovery tools have been employed to identify novel potential calcium sensitizing agents binding to the TnC-TnI interface. Finally, Markov modeling has contributed to simulating contraction within the sarcomere on the mesoscale. Here we are reviewing and classifying the existing computational work on troponin and its subunits, outline current gaps in simulations elucidating troponin's role in contraction and suggest potential future developments in the field.
Collapse
Affiliation(s)
- Jacob D Bowman
- Department of Chemistry and Biochemistry, Ohio State University, Columbus, OH, United States
| | - Steffen Lindert
- Department of Chemistry and Biochemistry, Ohio State University, Columbus, OH, United States
| |
Collapse
|
4
|
Bowman JD, Coldren WH, Lindert S. Mechanism of Cardiac Troponin C Calcium Sensitivity Modulation by Small Molecules Illuminated by Umbrella Sampling Simulations. J Chem Inf Model 2019; 59:2964-2972. [PMID: 31141358 DOI: 10.1021/acs.jcim.9b00256] [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/14/2022]
Abstract
Cardiac troponin C (cTnC) binds intracellular calcium and subsequently cardiac troponin I (cTnI), initiating cardiac muscle contraction. Due to its role in contraction, cTnC has been a therapeutic target in the search for small molecules to treat conditions that interfere with normal muscle contraction like the heritable cardiomyopathies. Structural studies have shown the binding location of small molecules such as bepridil, dfbp-o, 3-methyldiphenylamine (DPA), and W7 to be a hydrophobic pocket in the regulatory domain of cTnC (cNTnC) but have not shown the influence of these small molecules on the energetics of opening this domain. Here we describe an application of an umbrella sampling method used to elucidate the impact these calcium sensitivity modulators have on the free energy of cNTnC hydrophobic patch opening. We found that all these molecules lowered the free energy of opening in the absence of the cTnI, with bepridil facilitating the least endergonic transformation. In the presence of cTnI, however, we saw a stabilization of the open configuration due to DPA and dfbp-o binding, and a destabilization of the open configuration imparted by bepridil and W7. Predicted poor binding molecule NSC34337 left the hydrophobic patch in under 3 ns in conventional MD simulations suggesting that only hydrophobic patch binders stabilized the open conformation. In conclusion, this study presents a novel approach to study the impact of small molecules on hydrophobic patch opening through umbrella sampling, and it proposes mechanisms for calcium sensitivity modulation.
Collapse
Affiliation(s)
- Jacob D Bowman
- Department of Chemistry and Biochemistry , Ohio State University , 2114 Newman & Wolfrom Laboratory, 100 West 18th Avenue , Columbus , Ohio 43210 , United States
| | - William H Coldren
- Department of Chemistry and Biochemistry , Ohio State University , 2114 Newman & Wolfrom Laboratory, 100 West 18th Avenue , Columbus , Ohio 43210 , United States
| | - Steffen Lindert
- Department of Chemistry and Biochemistry , Ohio State University , 2114 Newman & Wolfrom Laboratory, 100 West 18th Avenue , Columbus , Ohio 43210 , United States
| |
Collapse
|
5
|
Bowman JD, Lindert S. Molecular Dynamics and Umbrella Sampling Simulations Elucidate Differences in Troponin C Isoform and Mutant Hydrophobic Patch Exposure. J Phys Chem B 2018; 122:7874-7883. [PMID: 30070845 DOI: 10.1021/acs.jpcb.8b05435] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Troponin C (TnC) facilitates muscle contraction through calcium-binding within its N-terminal region (NTnC). As previously observed using molecular dynamics (MD) simulations, this calcium-binding event leads to an increase in the dynamics of helices lining a hydrophobic patch on TnC. Simulation times of multiple microseconds were required to even see a partial opening of the hydrophobic patch, limiting the ability to thoroughly and quantitatively investigate these rare events. Here we describe the application of umbrella sampling to probe the TnC hydrophobic patch opening in a more targeted and quantitative fashion. Umbrella sampling was utilized to investigate the differences in the free energy of opening between cardiac (cTnC) and fast skeletal TnC (sTnC). We found that, in agreement with previous reports, holo (calcium-bound) sTnC had a lower free energy of opening compared with holo cTnC. Additionally, differences in the free energy of opening of hypertrophic (HCM) and dilated cardiomyopathy (DCM) cTnC systems were investigated. MD simulations and umbrella sampling revealed a lower free energy of opening for the HCM mutations A8V and A31S, as well as the calcium-sensitizing mutation L48Q. The DCM mutations, Y5H, Q50R, and E59D/D75Y, all exhibited a higher free energy of opening. An umbrella sampling simulation of cTnI-bound holo cTnC exhibited the lowest free energy in the open configuration, in agreement with experimental data. In conclusion, this study presents a novel and successful protocol for applying umbrella sampling simulations to quantitatively study the molecular basis of muscle contraction and proposes a mechanism by which HCM and DCM-associated mutations influence contraction.
Collapse
Affiliation(s)
- Jacob D Bowman
- Department of Chemistry and Biochemistry , Ohio State University , Columbus , Ohio 43210 , United States
| | - Steffen Lindert
- Department of Chemistry and Biochemistry , Ohio State University , Columbus , Ohio 43210 , United States
| |
Collapse
|
6
|
Karami M, Jalali C, Mirzaie S. Combined virtual screening, MMPBSA, molecular docking and dynamics studies against deadly anthrax: An in silico effort to inhibit Bacillus anthracis nucleoside hydrolase. J Theor Biol 2017; 420:180-189. [DOI: 10.1016/j.jtbi.2017.03.010] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2016] [Revised: 03/03/2017] [Accepted: 03/10/2017] [Indexed: 10/20/2022]
|
7
|
Dewan S, McCabe KJ, Regnier M, McCulloch AD. Insights and Challenges of Multi-Scale Modeling of Sarcomere Mechanics in cTn and Tm DCM Mutants-Genotype to Cellular Phenotype. Front Physiol 2017; 8:151. [PMID: 28352236 PMCID: PMC5348544 DOI: 10.3389/fphys.2017.00151] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2016] [Accepted: 02/24/2017] [Indexed: 01/18/2023] Open
Abstract
Dilated Cardiomyopathy (DCM) is a leading cause of sudden cardiac death characterized by impaired pump function and dilatation of cardiac ventricles. In this review we discuss various in silico approaches to elucidating the mechanisms of genetic mutations leading to DCM. The approaches covered in this review focus on bridging the spatial and temporal gaps that exist between molecular and cellular processes. Mutations in sarcomeric regulatory thin filament proteins such as the troponin complex (cTn) and Tropomyosin (Tm) have been associated with DCM. Despite the experimentally-observed myofilament measures of contractility in the case of these mutations, the mechanisms by which the underlying molecular changes and protein interactions scale up to organ failure by these mutations remains elusive. The review highlights multi-scale modeling approaches and their applicability to study the effects of sarcomeric gene mutations in-silico. We discuss some of the insights that can be gained from computational models of cardiac biomechanics when scaling from molecular states to cellular level.
Collapse
Affiliation(s)
- Sukriti Dewan
- Departments of Bioengineering and Medicine, University of California San Diego, La Jolla, CA, USA
| | - Kimberly J McCabe
- Departments of Bioengineering and Medicine, University of California San Diego, La Jolla, CA, USA
| | - Michael Regnier
- Departments of Bioengineering and Medicine, University of Washington Seattle, WA, USA
| | - Andrew D McCulloch
- Departments of Bioengineering and Medicine, University of California San Diego, La Jolla, CA, USA
| |
Collapse
|
8
|
Dewan S, McCabe KJ, Regnier M, McCulloch AD, Lindert S. Molecular Effects of cTnC DCM Mutations on Calcium Sensitivity and Myofilament Activation-An Integrated Multiscale Modeling Study. J Phys Chem B 2016; 120:8264-75. [PMID: 27133568 PMCID: PMC5001916 DOI: 10.1021/acs.jpcb.6b01950] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Mutations in cardiac troponin C (D75Y, E59D, and G159D), a key regulatory protein of myofilament contraction, have been associated with dilated cardiomyopathy (DCM). Despite reports of altered myofilament function in these mutants, the underlying molecular alterations caused by these mutations remain elusive. Here we investigate in silico the intramolecular mechanisms by which these mutations affect myofilament contraction. On the basis of the location of cardiac troponin C (cTnC) mutations, we tested the hypothesis that intramolecular effects can explain the altered myofilament calcium sensitivity of force development for D75Y and E59D cTnC, whereas altered cardiac troponin C-troponin I (cTnC-cTnI) interaction contributes to the reported contractile effects of the G159D mutation. We employed a multiscale approach combining molecular dynamics (MD) and Brownian dynamics (BD) simulations to estimate cTnC calcium association and hydrophobic patch opening. We then integrated these parameters into a Markov model of myofilament activation to compute the steady-state force-pCa relationship. The analysis showed that myofilament calcium sensitivity with D75Y and E59D can be explained by changes in calcium binding affinity of cTnC and the rate of hydrophobic patch opening, if a partial cTnC interhelical opening angle (110°) is sufficient for cTnI switch peptide association to cTnC. In contrast, interactions between cTnC and cTnI within the cardiac troponin complex must also be accounted for to explain contractile alterations due to G159D. In conclusion, this is the first multiscale in silico study to elucidate how direct molecular effects of genetic mutations in cTnC translate to altered myofilament contractile function.
Collapse
Affiliation(s)
- Sukriti Dewan
- Department of Bioengineering, University of California at San Diego, La Jolla, CA, 92093
| | - Kimberly J. McCabe
- Department of Bioengineering, University of California at San Diego, La Jolla, CA, 92093
| | - Michael Regnier
- Dept. of Bioengineering, University of Washington, Seattle, WA 98195
- Center for Cardiovascular Biology, University of Washington, Seattle, WA 98109
| | - Andrew D. McCulloch
- Department of Bioengineering, University of California at San Diego, La Jolla, CA, 92093
| | - Steffen Lindert
- Department of Chemistry & Biochemistry, Ohio State University, Columbus, OH, 43210
| |
Collapse
|
9
|
Cheng Y, Regnier M. Cardiac troponin structure-function and the influence of hypertrophic cardiomyopathy associated mutations on modulation of contractility. Arch Biochem Biophys 2016; 601:11-21. [PMID: 26851561 PMCID: PMC4899195 DOI: 10.1016/j.abb.2016.02.004] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2015] [Revised: 01/30/2016] [Accepted: 02/02/2016] [Indexed: 11/29/2022]
Abstract
Cardiac troponin (cTn) acts as a pivotal regulator of muscle contraction and relaxation and is composed of three distinct subunits (cTnC: a highly conserved Ca(2+) binding subunit, cTnI: an actomyosin ATPase inhibitory subunit, and cTnT: a tropomyosin binding subunit). In this mini-review, we briefly summarize the structure-function relationship of cTn and its subunits, its modulation by PKA-mediated phosphorylation of cTnI, and what is known about how these properties are altered by hypertrophic cardiomyopathy (HCM) associated mutations of cTnI. This includes recent work using computational modeling approaches to understand the atomic-based structural level basis of disease-associated mutations. We propose a viewpoint that it is alteration of cTnC-cTnI interaction (rather than the Ca(2+) binding properties of cTn) per se that disrupt the ability of PKA-mediated phosphorylation at cTnI Ser-23/24 to alter contraction and relaxation in at least some HCM-associated mutations. The combination of state of the art biophysical approaches can provide new insight on the structure-function mechanisms of contractile dysfunction resulting cTnI mutations and exciting new avenues for the diagnosis, prevention, and even treatment of heart diseases.
Collapse
Affiliation(s)
- Yuanhua Cheng
- University of Washington, Department of Bioengineering, Seattle, WA, USA
| | - Michael Regnier
- University of Washington, Department of Bioengineering, Seattle, WA, USA.
| |
Collapse
|
10
|
Cheng Y, Lindert S, Kekenes-Huskey P, Rao VS, Solaro RJ, Rosevear PR, Amaro R, McCulloch AD, McCammon JA, Regnier M. Computational studies of the effect of the S23D/S24D troponin I mutation on cardiac troponin structural dynamics. Biophys J 2014; 107:1675-85. [PMID: 25296321 PMCID: PMC4190606 DOI: 10.1016/j.bpj.2014.08.008] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2014] [Revised: 07/25/2014] [Accepted: 08/11/2014] [Indexed: 11/27/2022] Open
Abstract
During β-adrenergic stimulation, cardiac troponin I (cTnI) is phosphorylated by protein kinase A (PKA) at sites S23/S24, located at the N-terminus of cTnI. This phosphorylation has been shown to decrease KCa and pCa50, and weaken the cTnC-cTnI (C-I) interaction. We recently reported that phosphorylation results in an increase in the rate of early, slow phase of relaxation (kREL,slow) and a decrease in its duration (tREL,slow), which speeds up the overall relaxation. However, as the N-terminus of cTnI (residues 1-40) has not been resolved in the whole cardiac troponin (cTn) structure, little is known about the molecular-level behavior within the whole cTn complex upon phosphorylation of the S23/S24 residues of cTnI that results in these changes in function. In this study, we built up the cTn complex structure (including residues cTnC 1-161, cTnI 1-172, and cTnT 236-285) with the N-terminus of cTnI. We performed molecular-dynamics (MD) simulations to elucidate the structural basis of PKA phosphorylation-induced changes in cTn structure and Ca(2+) binding. We found that introducing two phosphomimic mutations into sites S23/S24 had no significant effect on the coordinating residues of Ca(2+) binding site II. However, the overall fluctuation of cTn was increased and the C-I interaction was altered relative to the wild-type model. The most significant changes involved interactions with the N-terminus of cTnI. Interestingly, the phosphomimic mutations led to the formation of intrasubunit interactions between the N-terminus and the inhibitory peptide of cTnI. This may result in altered interactions with cTnC and could explain the increased rate and decreased duration of slow-phase relaxation seen in myofibrils.
Collapse
Affiliation(s)
- Yuanhua Cheng
- Department of Bioengineering, University of Washington, Seattle, Washington; National Biomedical Computational Resource, University of California, San Diego, La Jolla, California
| | - Steffen Lindert
- National Biomedical Computational Resource, University of California, San Diego, La Jolla, California; Department of Pharmacology, University of California, San Diego, La Jolla, California
| | - Peter Kekenes-Huskey
- National Biomedical Computational Resource, University of California, San Diego, La Jolla, California; Department of Pharmacology, University of California, San Diego, La Jolla, California
| | - Vijay S Rao
- Department of Bioengineering, University of Washington, Seattle, Washington
| | - R John Solaro
- Department of Physiology and Biophysics, College of Medicine, University of Illinois at Chicago, Chicago, Illinois
| | - Paul R Rosevear
- Department of Molecular Genetics, Biochemistry, and Microbiology, University of Cincinnati, Cincinnati, Ohio
| | - Rommie Amaro
- National Biomedical Computational Resource, University of California, San Diego, La Jolla, California
| | - Andrew D McCulloch
- National Biomedical Computational Resource, University of California, San Diego, La Jolla, California; Department of Bioengineering, University of California, San Diego, La Jolla, California
| | - J Andrew McCammon
- National Biomedical Computational Resource, University of California, San Diego, La Jolla, California; Department of Pharmacology, University of California, San Diego, La Jolla, California
| | - Michael Regnier
- Department of Bioengineering, University of Washington, Seattle, Washington.
| |
Collapse
|
11
|
Sevrieva I, Knowles AC, Kampourakis T, Sun YB. Regulatory domain of troponin moves dynamically during activation of cardiac muscle. J Mol Cell Cardiol 2014; 75:181-7. [PMID: 25101951 PMCID: PMC4169182 DOI: 10.1016/j.yjmcc.2014.07.015] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 05/01/2014] [Revised: 07/09/2014] [Accepted: 07/26/2014] [Indexed: 11/24/2022]
Abstract
Heart muscle is activated by Ca2+ to generate force and shortening, and the signaling pathway involves allosteric mechanisms in the thin filament. Knowledge about the structure-function relationship among proteins in the thin filament is critical in understanding the physiology and pathology of the cardiac function, but remains obscure. We investigate the conformation of the cardiac troponin (Tn) on the thin filament and its response to Ca2+ activation and propose a molecular mechanism for the regulation of cardiac muscle contraction by Tn based uniquely on information from in situ protein domain orientation. Polarized fluorescence from bifunctional rhodamine is used to determine the orientation of the major component of Tn core domain on the thin filaments of cardiac muscle. We show that the C-terminal lobe of TnC (CTnC) does not move during activation, suggesting that CTnC, together with the coiled coil formed by the TnI and TnT chains (IT arm), acts as a scaffold that holds N-terminal lobe of TnC (NTnC) and the actin binding regions of troponin I. The NTnC, on the other hand, exhibits multiple orientations during both diastole and systole. By combining the in situ orientation data with published in vitro measurements of intermolecular distances, we construct a model for the in situ structure of the thin filament. The conformational dynamics of NTnC plays an important role in the regulation of cardiac muscle contraction by moving the C-terminal region of TnI from its actin-binding inhibitory location and enhancing the movement of tropomyosin away from its inhibitory position. In situ conformational changes of troponin in myocardium were investigated. A model for the cardiac thin filament was constructed based on the in situ data. The IT arm of cardiac troponin acts as a scaffold that holds the regulatory domain. The regulatory domain of cardiac troponin moves dynamically during activation. The dynamics of regulatory domain is important in cardiac muscle regulation.
Collapse
Affiliation(s)
- Ivanka Sevrieva
- Randall Division of Cell and Molecular Biophysics and British Heart Foundation Centre of Research Excellence, King's College London, London, SE1 1UL, UK
| | - Andrea C Knowles
- Randall Division of Cell and Molecular Biophysics and British Heart Foundation Centre of Research Excellence, King's College London, London, SE1 1UL, UK
| | - Thomas Kampourakis
- Randall Division of Cell and Molecular Biophysics and British Heart Foundation Centre of Research Excellence, King's College London, London, SE1 1UL, UK
| | - Yin-Biao Sun
- Randall Division of Cell and Molecular Biophysics and British Heart Foundation Centre of Research Excellence, King's College London, London, SE1 1UL, UK.
| |
Collapse
|
12
|
Jayasundar JJ, Xing J, Robinson JM, Cheung HC, Dong WJ. Molecular dynamics simulations of the cardiac troponin complex performed with FRET distances as restraints. PLoS One 2014; 9:e87135. [PMID: 24558365 PMCID: PMC3928104 DOI: 10.1371/journal.pone.0087135] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2013] [Accepted: 12/18/2013] [Indexed: 11/22/2022] Open
Abstract
Cardiac troponin (cTn) is the Ca2+-sensitive molecular switch that controls cardiac muscle activation and relaxation. However, the molecular detail of the switching mechanism and how the Ca2+ signal received at cardiac troponin C (cTnC) is communicated to cardiac troponin I (cTnI) are still elusive. To unravel the structural details of troponin switching, we performed ensemble Förster resonance energy transfer (FRET) measurements and molecular dynamic (MD) simulations of the cardiac troponin core domain complex. The distance distributions of forty five inter-residue pairs were obtained under Ca2+-free and saturating Ca2+ conditions from time-resolved FRET measurements. These distances were incorporated as restraints during the MD simulations of the cardiac troponin core domain. Compared to the Ca2+-saturated structure, the absence of regulatory Ca2+ perturbed the cTnC N-domain hydrophobic pocket which assumed a closed conformation. This event partially unfolded the cTnI regulatory region/switch. The absence of Ca2+, induced flexibility to the D/E linker and the cTnI inhibitory region, and rotated the cTnC N-domain with respect to rest of the troponin core domain. In the presence of saturating Ca2+ the above said phenomenon were absent. We postulate that the secondary structure perturbations experienced by the cTnI regulatory region held within the cTnC N-domain hydrophobic pocket, coupled with the rotation of the cTnC N-domain would control the cTnI mobile domain interaction with actin. Concomitantly the rotation of the cTnC N-domain and perturbation of the D/E linker rigidity would control the cTnI inhibitory region interaction with actin to effect muscle relaxation.
Collapse
Affiliation(s)
- Jayant James Jayasundar
- Voiland School of Chemical Engineering and Bioengineering and The Department of Integrated Physiology and Neuroscience, Washington State University, Pullman, Washington, United States of America
| | - Jun Xing
- Voiland School of Chemical Engineering and Bioengineering and The Department of Integrated Physiology and Neuroscience, Washington State University, Pullman, Washington, United States of America
| | - John M. Robinson
- Department of Chemistry and Biochemistry, South Dakota State University, Brookings, South Dakota, United States of America
| | - Herbert C. Cheung
- The Department of Biochemistry and Molecular Genetics, University of Alabama at Birmingham, Birmingham, Alabama, United States of America
| | - Wen-Ji Dong
- Voiland School of Chemical Engineering and Bioengineering and The Department of Integrated Physiology and Neuroscience, Washington State University, Pullman, Washington, United States of America
- * E-mail:
| |
Collapse
|
13
|
Genchev GZ, Kobayashi T, Lu H. Calcium induced regulation of skeletal troponin--computational insights from molecular dynamics simulations. PLoS One 2013; 8:e58313. [PMID: 23554884 PMCID: PMC3598806 DOI: 10.1371/journal.pone.0058313] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2012] [Accepted: 02/01/2013] [Indexed: 01/11/2023] Open
Abstract
The interaction between calcium and the regulatory site(s) of striated muscle regulatory protein troponin switches on and off muscle contraction. In skeletal troponin binding of calcium to sites I and II of the TnC subunit results in a set of structural changes in the troponin complex, displaces tropomyosin along the actin filament and allows myosin-actin interaction to produce mechanical force. In this study, we used molecular dynamics simulations to characterize the calcium dependent dynamics of the fast skeletal troponin molecule and its TnC subunit in the calcium saturated and depleted states. We focused on the N-lobe and on describing the atomic level events that take place subsequent to removal of the calcium ion from the regulatory sites I and II. A main structural event - a closure of the A/B helix hydrophobic pocket results from the integrated effect of the following conformational changes: the breakage of H-bond interactions between the backbone nitrogen atoms of the residues at positions 2, 9 and sidechain oxygen atoms of the residue at position 12 (N2-OE12/N9-OE12) in sites I and II; expansion of sites I and II and increased site II N-terminal end-segment flexibility; strengthening of the β-sheet scaffold; and the subsequent re-packing of the N-lobe hydrophobic residues. Additionally, the calcium release allows the N-lobe to rotate relative to the rest of the Tn molecule. Based on the findings presented herein we propose a novel model of skeletal thin filament regulation.
Collapse
Affiliation(s)
- Georgi Z. Genchev
- Bioinformatics Program, Department of Bioengineering, University of Illinois at Chicago, Chicago, Illinois, United States of America
| | - Tomoyoshi Kobayashi
- Department of Physiology and Biophysics and Center for Cardiovascular Research, University of Illinois at Chicago, Chicago, Illinois, United States of America
- * E-mail: (HL); (TK)
| | - Hui Lu
- Bioinformatics Program, Department of Bioengineering, University of Illinois at Chicago, Chicago, Illinois, United States of America
- Shanghai Institute of Medical Genetics, Children’s Hospital of Shanghai, Shanghai, China
- Key Lab of Embryo Molecular Biology, Ministry of Health, Shanghai, China
- Shanghai Lab of Embryo and Reproduction Engineering, Shanghai, China
- * E-mail: (HL); (TK)
| |
Collapse
|
14
|
Skowronsky RA, Schroeter M, Baxley T, Li Y, Chalovich JM, Spuches AM. Thermodynamics and molecular dynamics simulations of calcium binding to the regulatory site of human cardiac troponin C: evidence for communication with the structural calcium binding sites. J Biol Inorg Chem 2012; 18:49-58. [PMID: 23111626 DOI: 10.1007/s00775-012-0948-2] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2012] [Accepted: 09/30/2012] [Indexed: 01/10/2023]
Abstract
Human cardiac troponin C (HcTnC), a member of the EF hand family of proteins, is a calcium sensor responsible for initiating contraction of the myocardium. Ca(2+) binding to the regulatory domain induces a slight change in HcTnC conformation which modifies subsequent interactions in the troponin-tropomyosin-actin complex. Herein, we report a calorimetric study of Ca(2+) binding to HcTnC. Isotherms obtained at 25 °C (10 mM 2-morpholinoethanesulfonic acid, 50 mM KCl, pH 7.0) provided thermodynamic parameters for Ca(2+) binding to both the high-affinity and the low-affinity domain of HcTnC. Ca(2+) binding to the N-domain was shown to be endothermic in 2-morpholinoethanesulfonic acid buffer and allowed us to extract the thermodynamics of Ca(2+) binding to the regulatory domain. This pattern stems from changes that occur at the Ca(2+) site rather than structural changes of the protein. Molecular dynamics simulations performed on apo and calcium-bound HcTnC(1-89) support this claim. The values of the Gibbs free energy for Ca(2+) binding to the N-domain in the full-length protein and to the isolated domain (HcTnC(1-89)) are similar; however, differences in the entropic and enthalpic contributions to the free energy provide supporting evidence for the cooperativity of the C-domain and the N-domain. Thermograms obtained at two additional temperatures (10 and 37 °C) revealed interesting trends in the enthalpies and entropies of binding for both thermodynamic events. This allowed the determination of the change in heat capacity (∆C(p)) from a plot of ∆H verses temperature and may provide evidence for positive cooperativity of Ca(2+) binding to the C-domain.
Collapse
Affiliation(s)
- Rachel A Skowronsky
- Department of Chemistry, East Carolina University, 300 Science and Technology, Greenville, NC 27858, USA
| | | | | | | | | | | |
Collapse
|
15
|
John S, Thangapandian S, Lee KW. Potential human cholesterol esterase inhibitor design: benefits from the molecular dynamics simulations and pharmacophore modeling studies. J Biomol Struct Dyn 2012; 29:921-36. [PMID: 22292952 DOI: 10.1080/07391102.2012.10507419] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Abstract
Human pancreatic cholesterol esterase (hCEase) is one of the lipases found to involve in the digestion of large and broad spectrum of substrates including triglycerides, phospholipids, cholesteryl esters, etc. The presence of bile salts is found to be very important for the activation of hCEase. Molecular dynamic simulations were performed for the apoform and bile salt complexed form of hCEase using the co-ordinates of two bile salts from bovine CEase. The stability of the systems throughout the simulation time was checked and two representative structures from the highly populated regions were selected using cluster analysis. These two representative structures were used in pharmacophore model generation. The generated pharmacophore models were validated and used in database screening. The screened hits were refined for their drug-like properties based on Lipinski's rule of five and ADMET properties. The drug-like compounds were further refined by molecular docking simulation using GOLD program based on the GOLD fitness score, mode of binding, and molecular interactions with the active site amino acids. Finally, three hits of novel scaffolds were selected as potential leads to be used in novel and potent hCEase inhibitor design. The stability of binding modes and molecular interactions of these final hits were re-assured by molecular dynamics simulations.
Collapse
Affiliation(s)
- Shalini John
- Division of Applied Life Science_(BK21 Program), Systems and Synthetic Agrobiotech Center (SSAC) Gyeongsang National University (GNU), 501 Jinju-daero, Gazha-dong, Jinju 660-701, Republic of Korea
| | | | | |
Collapse
|
16
|
Matos KS, da Cunha EF, da Silva Gonçalves A, Wilter A, Kuča K, França TC, Ramalho TC. First principles calculations of thermodynamics and kinetic parameters and molecular dynamics simulations of acetylcholinesterase reactivators: can mouse data provide new insights into humans? J Biomol Struct Dyn 2012; 30:546-58. [DOI: 10.1080/07391102.2012.687521] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
|
17
|
Wang ZJ, Si YX, Oh S, Yang JM, Yin SJ, Park YD, Lee J, Qian GY. The effect of fucoidan on tyrosinase: computational molecular dynamics integrating inhibition kinetics. J Biomol Struct Dyn 2012; 30:460-73. [PMID: 22694253 DOI: 10.1080/07391102.2012.682211] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Fucoidan is a complex sulfated polysaccharide extracted from brown seaweed and has a wide variety of biological activities. In this study, we investigated the inhibitory effect of fucoidan on tyrosinase via a combination of inhibition kinetics and computational simulations. Fucoidan reversibly inhibited tyrosinase in a mixed-type manner. Time-interval kinetics showed that the inhibition was processed as first order with biphasic processes. For further insight, we simulated dockings with various sizes of molecular models (monomer to decamer) of fucoidan and showed that the best binding energy change results were obtained from the pentamer (-1.89 kcal/mol) and the hexamer (-1.97 kcal/mol) models of AutoDock Vina. The molecular dynamics simulation confirmed the binding mechanisms between tyrosinase and fucoidan and suggested that fucoidan mostly interacts with several residues including copper ions located in the active site. Our study suggests that fucoidan might be a potential natural antipigment agent.
Collapse
Affiliation(s)
- Zhi-Jiang Wang
- College of Biological and Environmental Sciences, Zhejiang Wanli University, Ningbo, 315100, P.R. China
| | | | | | | | | | | | | | | |
Collapse
|
18
|
Chalovich JM. Disease causing mutations of troponin alter regulated actin state distributions. J Muscle Res Cell Motil 2012; 33:493-9. [PMID: 22678497 DOI: 10.1007/s10974-012-9305-x] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2012] [Accepted: 05/24/2012] [Indexed: 11/29/2022]
Abstract
Striated muscle contraction is regulated primarily through the action of tropomyosin and troponin that are bound to actin. Activation requires Ca(2+) binding to troponin and/or binding of high affinity myosin complexes to actin. Mutations within components of the regulatory complex may lead to familial cardiomyopathies and myopathies. In several cases examined, either physiological or pathological changes in troponin alter the distribution among states of actin-tropomyosin-troponin that differ in their abilities to stimulate myosin ATPase activity. These observations open possibilities for managing disorders of the troponin complex. Furthermore, analyses of mutant forms of troponin give insights into the regulation of striated muscle contraction.
Collapse
Affiliation(s)
- Joseph M Chalovich
- Department of Biochemistry and Molecular Biology, Brody School of Medicine at East Carolina University, Greenville, NC, USA.
| |
Collapse
|
19
|
Franklin A, Baxley T, Kobayashi T, Chalovich J. The C-terminus of troponin T is essential for maintaining the inactive state of regulated actin. Biophys J 2012; 102:2536-44. [PMID: 22713569 PMCID: PMC3368147 DOI: 10.1016/j.bpj.2012.04.037] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2011] [Revised: 03/23/2012] [Accepted: 04/24/2012] [Indexed: 01/11/2023] Open
Abstract
Striated muscle contraction is regulated by the actin binding proteins tropomyosin and troponin. Defects in these proteins lead to myopathies and cardiomyopathies. Deletion of the 14 C-terminal residues of cardiac troponin T leads to hypertrophic cardiomyopathy. We showed earlier that regulated actin containing Δ14 TnT was more readily activated than wild-type regulated actin. We suggested that the equilibria among the inactive (blocked), intermediate (closed or calcium), and active (open or myosin) states was shifted to the active state. We now show that, in addition, such regulated actin filaments cannot enter the inactive or blocked state. Regulated actin containing Δ14 TnT had ATPase activities in the absence of Ca2+ that were higher than wild-type filaments but far below the fully active rate. The rapid dissociation of S1-ATP from regulated actin filaments containing Δ14 TnT and acrylodan-labeled tropomyosin did not show the fluorescence increase characteristic of moving to the inactive state. Replacing wild-type TnI with S45E TnI, that favors the inactive state, did not restore the fluorescence change. We conclude that TnT has a previously unrecognized role in forming the inactive state of regulated actin.
Collapse
Affiliation(s)
- Andrew J. Franklin
- Department of Biochemistry and Molecular Biology, Brody School of Medicine, East Carolina University, Greenville, North Carolina
| | - Tamatha Baxley
- Department of Biochemistry and Molecular Biology, Brody School of Medicine, East Carolina University, Greenville, North Carolina
| | - Tomoyoshi Kobayashi
- Department of Physiology and Biophysics and Center for Cardiovascular Research, College of Medicine, University of Illinois, Chicago, Illinois
| | - Joseph M. Chalovich
- Department of Biochemistry and Molecular Biology, Brody School of Medicine, East Carolina University, Greenville, North Carolina
| |
Collapse
|
20
|
Majumder R, Roy S, Thakur AR. Analysis of Delta–Notch interaction by molecular modeling and molecular dynamic simulation studies. J Biomol Struct Dyn 2012; 30:13-29. [DOI: 10.1080/07391102.2012.674184] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
21
|
Serikov R, Petyuk V, Vorobijev Y, Koval V, Fedorova O, Vlassov V, Zenkova M. Mechanism of antisense oligonucleotide interaction with natural RNAs. J Biomol Struct Dyn 2011; 29:27-50. [PMID: 21696224 DOI: 10.1080/073911011010524987] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
Oligonucleotides find several numbers of applications: as diagnostic probes, RT and PCR primers and antisense agents due to their ability of forming specific interactions with complementary nucleotide sequences within nucleic acids. These interactions are strongly affected by accessibility of the target sequence in the RNA structure. In the present work the mechanism of invasion of RNA structure by oligonucleotide was investigated using a model system: yeast tRNA(Phe) and oligonucleotides complementary to the 3'-part of this molecule. Kinetics of interaction of oligonucleotides with in vitro transcript of yeast tRNAPhe was studied using stopped-flow technique with fluorescence quenching detection, 5'-DABCYL labeled oligonucleotide was hybridized with 3'-fluorescein labeled tRNA(Phe). The results evidence for a four-step invasion process of the oligonucleotide-RNA complex formation. The process is initiated by formation of transition complexes with nucleotides in the T-loop and ACCA sequence. This complex formation is followed by RNA unfolding and formation of an extended heteroduplex with the oligonucleotide via strand displacement process. Computer modeling of oligonucleotide-tRNA(Phe) interaction revealed potential factors that could favor transition complexes formation and confirmed the proposed mechanism, showing the oligonucleotide to be a molecular "wedge". Our data evidence that oligonucleotide invasion into structured RNA is initiated by loop-single strand interactions, similar to the initial step of the antisense RNA-RNA interactions. The obtained results can be used for choosing efficient oligonucleotide probes.
Collapse
Affiliation(s)
- R Serikov
- Institute of Chemical Biology and Fundamental Medicine SB RAS, 8 Lavrentiev Ave., 630090 Novosibirsk, Russian Federation
| | | | | | | | | | | | | |
Collapse
|
22
|
Semighini EP, Resende JA, de Andrade P, Morais PAB, Carvalho I, Taft CA, Silva CHTP. Using computer-aided drug design and medicinal chemistry strategies in the fight against diabetes. J Biomol Struct Dyn 2011; 28:787-96. [PMID: 21294589 DOI: 10.1080/07391102.2011.10508606] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
The aim of this work is to present a simple, practical and efficient protocol for drug design, in particular Diabetes, which includes selection of the illness, good choice of a target as well as a bioactive ligand and then usage of various computer aided drug design and medicinal chemistry tools to design novel potential drug candidates in different diseases. We have selected the validated target dipeptidyl peptidase IV (DPP-IV), whose inhibition contributes to reduce glucose levels in type 2 diabetes patients. The most active inhibitor with complex X-ray structure reported was initially extracted from the BindingDB database. By using molecular modification strategies widely used in medicinal chemistry, besides current state-of-the-art tools in drug design (including flexible docking, virtual screening, molecular interaction fields, molecular dynamics, ADME and toxicity predictions), we have proposed 4 novel potential DPP-IV inhibitors with drug properties for Diabetes control, which have been supported and validated by all the computational tools used herewith.
Collapse
Affiliation(s)
- Evandro P Semighini
- Departamento de Ciencias Farmaceuticas, Faculdade de Ciencias Farmaceuticas de Ribeirao Preto, Universidade de Sao Paulo, Av. do Cafe, s/n, Monte Alegre, 14040-903, Ribeirao Preto-SP, Brazil
| | | | | | | | | | | | | |
Collapse
|
23
|
Xu X, Su J, Chen W, Wang C. Thermal stability and unfolding pathways of Sso7d and its mutant F31A: insight from molecular dynamics simulation. J Biomol Struct Dyn 2011; 28:717-27. [PMID: 21294584 DOI: 10.1080/07391102.2011.10508601] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
The thermo-stability and unfolding behaviors of a small hyperthermophilic protein Sso7d as well as its single-point mutation F31A are studied by molecular dynamics simulation at temperatures of 300 K, 371 K and 500 K. Simulations at 300 K show that the F31A mutant displays a much larger flexibility than the wild type, which implies that the mutation obviously decreases the protein's stability. In the simulations at 371 K, although larger fluctuations were observed, both of these two maintain their stable conformations. High temperature simulations at 500 K suggest that the unfolding of these two proteins evolves along different pathways. For the wild-type protein, the C-terminal alpha-helix is melted at the early unfolding stage, whereas it is destroyed much later in the unfolding process of the F31A mutant. The results also show that the mutant unfolds much faster than its parent protein. The deeply buried aromatic cluster in the F31A mutant dissociates quickly relative to the wild-type protein at high temperature. Besides, it is found that the triple-stranded antiparallel β-sheet in the wild-type protein plays an important role in maintaining the stability of the entire structure.
Collapse
Affiliation(s)
- Xianjin Xu
- College of Life Science and Bioengineering, Beijing University of Technology, Beijing 100124, China
| | | | | | | |
Collapse
|
24
|
Majumder R, Roy S, Thakur AR. Molecular Modeling and Molecular Dynamics Simulation Studies of Delta-Notch Complex. J Biomol Struct Dyn 2011; 29:297-310. [DOI: 10.1080/07391102.2011.10507386] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
|
25
|
Behmard E, Abdolmaleki P, Asadabadi EB, Jahandideh S. Prevalent Mutations of Human Prion Protein: A Molecular Modeling and Molecular Dynamics Study. J Biomol Struct Dyn 2011; 29:379-89. [DOI: 10.1080/07391102.2011.10507392] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
|
26
|
Oliveira AA, Rennó MN, de Matos CAS, Bertuzzi MD, Ramalho TC, Fraga CA, França TCC. Molecular Modeling Studies ofYersinia pestisDihydrofolate Reductase. J Biomol Struct Dyn 2011; 29:351-67. [DOI: 10.1080/07391102.2011.10507390] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
|
27
|
Jani V, Sonavane UB, Joshi R. Microsecond scale replica exchange molecular dynamic simulation of villin headpiece: an insight into the folding landscape. J Biomol Struct Dyn 2011; 28:845-60. [PMID: 21469746 DOI: 10.1080/07391102.2011.10508612] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Abstract
Reaching the experimental time scale of millisecond is a grand challenge for protein folding simulations. The development of advanced Molecular Dynamics techniques like Replica Exchange Molecular Dynamics (REMD) makes it possible to reach these experimental timescales. In this study, an attempt has been made to reach the multi microsecond simulation time scale by carrying out folding simulations on a three helix bundle protein, Villin, by combining REMD and Amber United Atom model. Twenty replicas having different temperatures ranging from 295 K to 390 K were simulated for 1.5 µs each. The lowest Root Mean Square Deviation (RMSD) structure of 2.5 Å was obtained with respect to native structure (PDB code 1VII), with all the helices formed. The folding population landscapes were built using segment-wise RMSD and Principal Components as reaction coordinates. These analyses suggest the two-stage folding for Villin. The combination of REMD and Amber United Atom model may be useful to understand the folding mechanism of various fast folding proteins.
Collapse
Affiliation(s)
- Vinod Jani
- Bioinformatics Team, Scientific and Engineering Computing Group, Centre for Development of Advanced Computing, Pune University Campus, Ganeshkhind, Pune - 411007, Maharashtra, India
| | | | | |
Collapse
|
28
|
Su JG, Xu XJ, Li CH, Chen WZ, Wang CX. An Analysis of the Influence of Protein Intrinsic Dynamical Properties on its Thermal Unfolding Behavior. J Biomol Struct Dyn 2011; 29:105-21. [DOI: 10.1080/07391102.2011.10507377] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
|
29
|
Purohit R, Rajendran V, Sethumadhavan R. Studies on Adaptability of Binding Residues Flap Region of TMC-114 Resistance HIV-1 Protease Mutants. J Biomol Struct Dyn 2011; 29:137-52. [DOI: 10.1080/07391102.2011.10507379] [Citation(s) in RCA: 53] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
|
30
|
Lin CH, Chang TT, Sun MF, Chen HY, Tsai FJ, Chang KL, Fisher M, Chen CYC. Potent inhibitor design against H1N1 swine influenza: structure-based and molecular dynamics analysis for M2 inhibitors from traditional Chinese medicine database. J Biomol Struct Dyn 2011; 28:471-82. [PMID: 21142218 DOI: 10.1080/07391102.2011.10508589] [Citation(s) in RCA: 50] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
The rapid spread of influenza virus subtype H1N1 poses a great threat to million lives worldwide. To search for new anti-influenza compounds, we performed molecular docking and molecular dynamics simulation to identify potential traditional Chinese medicine (TCM) constituents that could block influenza M2 channel activity. Quinic acid, genipin, syringic acid, cucurbitine, fagarine, and methyl isoferulate all have extremely well docking results as compared to control amantadine. Further de novo drug design suggests that derivatives of genipin and methyl isoferulate could have enhanced binding affinity towards M2 channel. Selected molecular dynamics simulations of M2-derivative complexes show stable hydrogen bond interactions between the derivatives and M2 residues, Ser10 and Ala9. To our best knowledge, this is the first study on the anti-viral activity of the above listed TCM compounds.
Collapse
Affiliation(s)
- Chia-Hui Lin
- Department of Chinese Medicine, China Medical University Hospital, Taiwan
| | | | | | | | | | | | | | | |
Collapse
|
31
|
Zhou ZL, Zhao JH, Liu HL, Wu JW, Liu KT, Chuang CK, Tsai WB, Ho Y. The Possible Structural Models for Polyglutamine Aggregation: A Molecular Dynamics Simulations Study. J Biomol Struct Dyn 2011; 28:743-58. [DOI: 10.1080/07391102.2011.10508603] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
32
|
Olmez EO, Alakent B. Alpha7 Helix Plays an Important Role in the Conformational Stability of PTP1B. J Biomol Struct Dyn 2011; 28:675-93. [DOI: 10.1080/07391102.2011.10508599] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
|
33
|
Roy S, Thakur AR. Two models of Smad4 and Hoxa9 Complex are Proposed: Structural and Interactional Perspective. J Biomol Struct Dyn 2011; 28:729-42. [DOI: 10.1080/07391102.2011.10508602] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
|
34
|
Zhao Y, Gong Z, Xiao Y. Improvements of the Hierarchical Approach for Predicting RNA Tertiary Structure. J Biomol Struct Dyn 2011; 28:815-26. [DOI: 10.1080/07391102.2011.10508609] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
35
|
Varughese JF, Baxley T, Chalovich JM, Li Y. A computational and experimental approach to investigate bepridil binding with cardiac troponin. J Phys Chem B 2011; 115:2392-400. [PMID: 21332124 DOI: 10.1021/jp1094504] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Cardiac troponin is a Ca(2+)-dependent switch for the contraction in heart muscle and a potential target for drugs in the therapy of heart failure. Bepridil is a drug that binds to troponin and increases calcium sensitivity of muscle contraction. Because bepridil has been well studied, it is a good model for analysis by computational and experimental methods. Molecular dynamics (MD) simulations were performed on troponin complexes of different sizes in the presence and absence of bepridil bound within the hydrophobic pocket at the N-terminal domain of troponin C. About 100 ns of simulation trajectory data were generated, which were analyzed using cross-correlation analyses and MMPBSA and MMGBSA techniques. The results indicated that bepridil binding within the hydrophobic pocket of cardiac TnC decreases the interaction of TnC with TnI at both the N-domain of TnC and the C-domain of TnC, and decreases the correlations of motions among the segments of the troponin subunits. The estimated calcium-binding affinities using MMPBSA showed that bepridil has a sensitizing effect for the isolated system of TnC, but loses this effect for the complex. Our experimental measurements of calcium dissociation rates were consistent with that prediction. We also observed that while bepridil enhanced the troponin-tropomyosin-actin-activated ATPase activity of myosin S1 at low calcium concentrations it was slightly inhibitory at high calcium concentrations. Bepridil increases the ATPase activity and force generation in muscle fibers, but its effects appear to depend on the concentration of calcium.
Collapse
Affiliation(s)
- Jayson F Varughese
- Department of Chemistry, Brody School of Medicine, East Carolina University, Greenville, North Carolina 27858, USA
| | | | | | | |
Collapse
|
36
|
Guimarães AP, Oliveira AA, da Cunha EFF, Ramalho TC, França TCC. Design of New Chemotherapeutics Against the Deadly Anthrax Disease. Docking and Molecular Dynamics studies of Inhibitors Containing Pyrrolidine and Riboamidrazone Rings on Nucleoside Hydrolase fromBacillus anthracis. J Biomol Struct Dyn 2011; 28:455-69. [DOI: 10.1080/07391102.2011.10508588] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
37
|
Chakrabarti B, Bairagya HR, Mallik P, Mukhopadhyay BP, Bera AK. An Insight to Conserved Water Molecular Dynamics of Catalytic and Structural Zn+2ions in Matrix Metalloproteinase 13 of Human. J Biomol Struct Dyn 2011; 28:503-16. [DOI: 10.1080/07391102.2011.10508591] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|