Molecular dynamics studies on troponin (TnI-TnT-TnC) complexes: insight into the regulation of muscle contraction

Molecular dynamics provides new insights into the mechanism of calcium signal transduction and interdomain interactions in cardiac troponin

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 Ca²⁺ -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.

Computational Studies of Cardiac and Skeletal Troponin

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.

A Stochastic Multiscale Model of Cardiac Thin Filament Activation Using Brownian-Langevin Dynamics

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.

Mechanism of Cardiac Troponin C Calcium Sensitivity Modulation by Small Molecules Illuminated by Umbrella Sampling Simulations

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.

Molecular Dynamics and Umbrella Sampling Simulations Elucidate Differences in Troponin C Isoform and Mutant Hydrophobic Patch Exposure

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.

Insights and Challenges of Multi-Scale Modeling of Sarcomere Mechanics in cTn and Tm DCM Mutants-Genotype to Cellular Phenotype

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.

Molecular Effects of cTnC DCM Mutations on Calcium Sensitivity and Myofilament Activation-An Integrated Multiscale Modeling Study

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.

Cardiac troponin structure-function and the influence of hypertrophic cardiomyopathy associated mutations on modulation of contractility

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.

Computational studies of the effect of the S23D/S24D troponin I mutation on cardiac troponin structural dynamics

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.

Regulatory domain of troponin moves dynamically during activation of cardiac muscle

Heart muscle is activated by Ca²⁺ 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 Ca²⁺ 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.

Molecular dynamics simulations of the cardiac troponin complex performed with FRET distances as restraints

Cardiac troponin (cTn) is the Ca²⁺-sensitive molecular switch that controls cardiac muscle activation and relaxation. However, the molecular detail of the switching mechanism and how the Ca²⁺ 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 Ca²⁺-free and saturating Ca²⁺ 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 Ca²⁺-saturated structure, the absence of regulatory Ca²⁺ perturbed the cTnC N-domain hydrophobic pocket which assumed a closed conformation. This event partially unfolded the cTnI regulatory region/switch. The absence of Ca²⁺, 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 Ca²⁺ 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.

Calcium induced regulation of skeletal troponin--computational insights from molecular dynamics simulations

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 (N²-OE¹²/N⁹-OE¹²) 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.

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

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.

Potential human cholesterol esterase inhibitor design: benefits from the molecular dynamics simulations and pharmacophore modeling studies

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.

The effect of fucoidan on tyrosinase: computational molecular dynamics integrating inhibition kinetics

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.

Disease causing mutations of troponin alter regulated actin state distributions

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.

The C-terminus of troponin T is essential for maintaining the inactive state of regulated actin

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.

Mechanism of antisense oligonucleotide interaction with natural RNAs

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.

Using computer-aided drug design and medicinal chemistry strategies in the fight against diabetes

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.

Thermal stability and unfolding pathways of Sso7d and its mutant F31A: insight from molecular dynamics simulation

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.

Microsecond scale replica exchange molecular dynamic simulation of villin headpiece: an insight into the folding landscape

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.

Potent inhibitor design against H1N1 swine influenza: structure-based and molecular dynamics analysis for M2 inhibitors from traditional Chinese medicine database

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.

A computational and experimental approach to investigate bepridil binding with cardiac troponin

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.