Non-B DNA conformations analysis through molecular dynamics simulations

BackgroundNon-B DNA conformations are molecular structures that do not follow the canonical DNA double helix. Mutagenetic instability in nuclear and mitochondrial DNA (mtDNA) genomes has been associated with simple non-B DNA conformations, as hairpins or more complex structures, as G-quadruplexes. One of these structures is Structure A, a cloverleaf-like non-B conformation predicted for a 93-nt (nucleotide) stretch of the mtDNA control region 5′-peripheral domain. Structure A is embedded in a hot spot for the 3′ end of human mtDNA deletions revealing its importance in influencing the mutational instability of the mtDNA genome.MethodsTo better characterize Structure A, we predicted its 3D conformation using state-of-art methods and algorithms. The methodologic workflow consisted in the prediction of non-B conformations using molecular dynamics simulations. The conservation scores of alignments of the Structure A region in humans, primates, and mammals, was also calculated.ResultsOur results show that these computational methods are able to measure the stability of non-B conformations by using the level of base pairing during molecular dynamics. Structure A showed high stability and low flexibility correlated with high conservation scores in mammalian, more specifically in primate lineages.ConclusionsWe showed that 3D non-B conformations can be predicted and characterized by our methodology. This allowed the in-depth analysis of the structure A, and the main results showed the structure remains stable during the simulations.General significanceThe fine-scale atomic molecular determination of this type of non-B conformation opens the way to perform computational molecular studies that can show their involvement in mtDNA cellular mechanisms.     

Functional and structural insights revealed by molecular dynamics simulations of an essential RNA editing ligase in Trypanosoma brucei

RNA editing ligase 1 (TbREL1) is required for the survival of both the insect and bloodstream forms of Trypanosoma brucei, the parasite responsible for the devastating tropical disease African sleeping sickness. The type of RNA editing that TbREL1 is involved in is unique to the trypanosomes, and no close human homolog is known to exist. In addition, the high-resolution crystal structure revealed several unique features of the active site, making this enzyme a promising target for structure-based drug design. In this work, two 20 ns atomistic molecular dynamics (MD) simulations are employed to investigate the dynamics of TbREL1, both with and without the ATP substrate present. The flexibility of the active site, dynamics of conserved residues and crystallized water molecules, and the interactions between TbREL1 and the ATP substrate are investigated and discussed in the context of TbREL1's function. Differences in local and global motion upon ATP binding suggest that two peripheral loops, unique to the trypanosomes, may be involved in interdomain signaling events. Notably, a significant structural rearrangement of the enzyme's active site occurs during the apo simulations, opening an additional cavity adjacent to the ATP binding site that could be exploited in the development of effective inhibitors directed against this protozoan parasite. Finally, ensemble averaged electrostatics calculations over the MD simulations reveal a novel putative RNA binding site, a discovery that has previously eluded scientists. Ultimately, we use the insights gained through the MD simulations to make several predictions and recommendations, which we anticipate will help direct future experimental studies and structure-based drug discovery efforts against this vital enzyme.     

Structural dynamics of precursor and product of the RNA enzyme from the hepatitis delta virus as revealed by molecular dynamics simulations

The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA enzyme involved in the replication of a human pathogen, the hepatitis delta virus. Recent crystal structures of the precursor and product of self-cleavage, together with detailed kinetic analyses, have led to hypotheses on the catalytic strategies employed by the HDV ribozyme. We report molecular dynamics (MD) simulations (approximately 120 ns total simulation time) to test the plausibility that specific conformational rearrangements are involved in catalysis. Site-specific self-cleavage requires cytidine in position 75 (C75). A precursor simulation with unprotonated C75 reveals a rather weak dynamic binding of C75 in the catalytic pocket with spontaneous, transient formation of a H-bond between U-1(O2') and C75(N3). This H-bond would be required for C75 to act as the general base. Upon protonation in the precursor, C75H+ has a tendency to move towards its product location and establish a firm H-bonding network within the catalytic pocket. However, a C75H+(N3)-G1(O5') H-bond, which would be expected if C75 acted as a general acid catalyst, is not observed on the present simulation timescale. The adjacent loop L3 is relatively dynamic and may serve as a flexible structural element, possibly gated by the closing U20.G25 base-pair, to facilitate a conformational switch induced by a protonated C75H+. L3 also controls the electrostatic environment of the catalytic core, which in turn may modulate C75 base strength and metal ion binding. We find that a distant RNA tertiary interaction involving a protonated cytidine (C41) becomes unstable when left unprotonated, leading to disruptive conformational rearrangements adjacent to the catalytic core. A Na ion temporarily compensates for the loss of the protonated hydrogen bond, which is strikingly consistent with the experimentally observed synergy between low pH and high Na+ concentrations in mediating residual self-cleavage of the HDV ribozyme in the absence of divalents.     

Assembly of lipoprotein particles revealed by coarse-grained molecular dynamics simulations

High-density lipoproteins (HDL) function as cholesterol transporters, facilitating the removal of excess cholesterol from the body. Due to the heterogeneity of native HDL particles (both in size and shape), the details on how these protein-lipid particles form and the structure they assume in their lipid-associated states are not well characterized. We report here a study of the self-assembly of discoidal HDL particles using coarse-grained (CG) molecular dynamics. The microsecond simulations reveal the self-assembly of HDL particles from disordered protein-lipid complexes to form structures containing many of the features of the generally accepted double-belt model for discoidal HDL particles. HDL assembly is found to proceed in two broad steps, aggregation of proteins and lipids driven by the hydrophobic effect which occurs on a approximately 1 micros time scale, followed by the optimization of the protein structure driven by increasingly specific protein-protein interactions.     

Structure and interactions of RecA: plasticity revealed by molecular dynamics simulations

Eleven independent simulations, each involving three consecutive molecules in the RecA filament, carried out on the protein from Mycobacterium tuberculosis, Mycobacterium smegmatis and Escherichia coli and their Adenosine triphosphate (ATP) complexes, provide valuable information which is complementary to that obtained from crystal structures, in addition to confirming the robust common structural framework within which RecA molecules from different eubacteria function. Functionally important loops, which are largely disordered in crystal structures, appear to adopt in each simulation subsets of conformations from larger ensembles. The simulations indicate the possibility of additional interactions involving the P-loop which remains largely invariant. The phosphate tail of the ATP is firmly anchored on the loop while the nucleoside moiety exhibits substantial structural variability. The most important consequence of ATP binding is the movement of the ‘switch’ residue. The relevant simulations indicate the feasibility of a second nucleotide binding site, but the pathway between adjacent molecules in the filament involving the two nucleotide binding sites appears to be possible only in the mycobacterial proteins.     

Collective dynamics of EcoRI-DNA complex by elastic network model and molecular dynamics simulations

Anisotropic network model (ANM) is used to analyze the collective motions of restriction enzyme EcoRI in free form and in complex with DNA. For comparison, three independent molecular dynamics (MD) simulations, each of 1.5 ns duration, are also performed for the EcoRI-DNA complex in explicit water. Although high mobility (equilibrium fluctuations) of inner and outer loops that surround the DNA is consistent in both methods and experiments, MD runs sample different conformational subspaces from which reliable collective dynamics cannot be extracted. However, ANM employed on different conformations from MD simulations indicates very similar collective motions. The stems of the inner loops are quite immobile even in the free enzyme and form a large, almost fixed, pocket for DNA binding. As a result, the residues that make specific and non-specific interactions with the DNA exhibit very low fluctuations in the free enzyme. The vibrational entropy difference between the EcoRI complex and free protein + unkinked DNA is positive (favorable), which may partially counteract the unfavorable enthalpy difference of DNA kink formation. Dynamic domains in EcoRI complex and cross-correlations between residue fluctuations indicate possible means of communication between the distal active sites.     

Protective mechanisms of α,α-trehalose revealed by molecular dynamics simulations

Molecular dynamics simulations were carried out to model aqueous solution with different concentration of α,α-trehalose, one kind of non-reducing sugars possessing outstanding freeze-drying protective effect on biological system. The dihedral angles of the intraglycosidic linkage in trehalose were measured to estimate its structure rigidity. The dynamics and hydrogen bonding properties were studied by calculating the self-diffusion coefficient of trehalose and the distributions and lifetimes of various types of H-bonds in the solution. Through analysing the results as well as comparison with another common sugar sucrose, the freeze-drying protective mechanism of trehalose was explained at molecule level. First, trehalose is able to maintain the local structure around it as a frame due to its relatively rigid conformation. Second, the addition of trehalose restrains the water molecules from rearrangement as a result of low mobility, thus reduces the probability of freezing; trehalose has lower diffusion coefficient than water and bigger thermal diffusivity, which are favourable for vitrification. Third, the formation of H-bonds between trehalose and water and between trehalose molecules is the essence of the protective effect. Trehalose does not work via strengthening the H-bonds formed between water molecules (W–W H-bonds), instead of which it breaks the potential tetrahedral pattern of W–W H-bonds, thus suppresses the tendency of ice formation. It was also found that trehalose realises protective action better at higher concentration as far as this study is concerned.     

Backbone conformations and side chain flexibility of two somatostatin mimics investigated by molecular dynamics simulations

Molecular dynamics simulations with two designed somatostatin mimics, SOM230 and SMS 201‐995, were performed in explicit water for a total aggregated time of 208 ns. Analysis of the runs with SOM230 revealed the presence of two clusters of conformations. Strikingly, the two sampled conformers correspond to the two main X‐ray structures in the asymmetric unit of SMS 201‐995. Structural comparison between the residues of SOM230 and SMS 201‐995 provides an explanation for the high binding affinity of SOM230 to four of five somatostatin receptors. Similarly, cluster analysis of the simulations with SMS 201‐995 shows that the backbone of the peptide interconverts between its two main crystallographic conformers. The conformations of SMS 201‐995 sampled in the two clusters violated two different sets of NOE distance constraints in agreement with a previous NMR study. Differences in side chain fluctuations between SOM230 and SMS 201‐995 observed in the simulations may contribute to the relatively higher binding affinity of SOM230 to most somatostatin receptors. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

DNA tension-modulated translocation and loop extrusion by SMC complexes revealed by molecular dynamics simulations

Structural Maintenance of Chromosomes (SMC) complexes play essential roles in genome organization across all domains of life. To determine how the activities of these large (≈50 nm) complexes are controlled by ATP binding and hydrolysis, we developed a molecular dynamics model that accounts for conformational motions of the SMC and DNA. The model combines DNA loop capture with an ATP-induced ‘power stroke’ to translocate the SMC complex along DNA. This process is sensitive to DNA tension: at low tension (0.1 pN), the model makes loop-capture steps of average 60 nm and up to 200 nm along DNA (larger than the complex itself), while at higher tension, a distinct inchworm-like translocation mode appears. By tethering DNA to an experimentally-observed additional binding site (‘safety belt’), the model SMC complex can perform loop extrusion (LE). The dependence of LE on DNA tension is distinct for fixed DNA tension vs. fixed DNA end points: LE reversal occurs above 0.5 pN for fixed tension, while LE stalling without reversal occurs at about 2 pN for fixed end points. Our model matches recent experimental results for condensin and cohesin, and makes testable predictions for how specific structural variations affect SMC function.     

Dynamics of the base of ribosomal A-site finger revealed by molecular dynamics simulations and Cryo-EM

Helix 38 (H38) of the large ribosomal subunit, with a length of 110 Å, reaches the small subunit through intersubunit bridge B1a. Previous cryo-EM studies revealed that the tip of H38 moves by more than 10 Å from the non-ratcheted to the ratcheted state of the ribosome while mutational studies implicated a key role of flexible H38 in attenuation of translocation and in dynamical signaling between ribosomal functional centers. We investigate a region including the elbow-shaped kink-turn (Kt-38) in the Haloarcula marismortui archaeal ribosome, and equivalently positioned elbows in three eubacterial species, located at the H38 base. We performed explicit solvent molecular dynamics simulations on the H38 elbows in all four species. They are formed by at first sight unrelated sequences resulting in diverse base interactions but built with the same overall topology, as shown by X-ray crystallography. The elbows display similar fluctuations and intrinsic flexibilities in simulations indicating that the eubacterial H38 elbows are structural and dynamical analogs of archaeal Kt-38. We suggest that this structural element plays a pivotal role in the large motions of H38 and may act as fulcrum for the abovementioned tip motion. The directional flexibility inferred from simulations correlates well with the cryo-EM results.     

Recognition of aminoacyl-tRNA: a common molecular mechanism revealed by cryo-EM

The accuracy of ribosomal translation is achieved by an initial selection and a proofreading step, mediated by EF-Tu, which forms a ternary complex with aminoacyl(aa)-tRNA. To study the binding modes of different aa-tRNAs, we compared cryo-EM maps of the kirromycin-stalled ribosome bound with ternary complexes containing Phe-tRNA^Phe, Trp-tRNA^Trp, or Leu-tRNA^LeuI. The three maps suggest a common binding manner of cognate aa-tRNAs in their specific binding with both the ribosome and EF-Tu. All three aa-tRNAs have the same ‘loaded spring' conformation with a kink and twist between the D-stem and anticodon stem. The three complexes are similarly integrated in an interaction network, extending from the anticodon loop through h44 and protein S12 to the EF-Tu-binding CCA end of aa-tRNA, proposed to signal cognate codon–anticodon interaction to the GTPase centre and tune the accuracy of aa-tRNA selection.     

The gating mechanism of the bacterial mechanosensitive channel MscL revealed by molecular dynamics simulations

One of the ultimate goals of the study on mechanosensitive (MS) channels is to understand the biophysical mechanisms of how the MS channel protein senses forces and how the sensed force induces channel gating. The bacterial MS channel MscL is an ideal subject to reach this goal owing to its resolved 3D protein structure in the closed state on the atomic scale and large amounts of electrophysiological data on its gating kinetics. However, the structural basis of the dynamic process from the closed to open states in MscL is not fully understood. In this study, we performed molecular dynamics (MD) simulations on the initial process of MscL opening in response to a tension increase in the lipid bilayer. To identify the tension-sensing site(s) in the channel protein, we calculated interaction energy between membrane lipids and candidate amino acids (AAs) facing the lipids. We found that Phe78 has a conspicuous interaction with the lipids, suggesting that Phe78 is the primary tension sensor of MscL. Increased membrane tension by membrane stretch dragged radially the inner (TM1) and outer (TM2) helices of MscL at Phe78, and the force was transmitted to the pentagon-shaped gate that is formed by the crossing of the neighboring TM1 helices in the inner leaflet of the bilayer. The radial dragging force induced radial sliding of the crossing portions, leading to a gate expansion. Calculated energy for this expansion is comparable to an experimentally estimated energy difference between the closed and the first subconductance state, suggesting that our model simulates the initial step toward the full opening of MscL. The model also successfully mimicked the behaviors of a gain of function mutant (G22N) and a loss of function mutant (F78N), strongly supporting that our MD model did simulate some essential biophysical aspects of the mechano-gating in MscL.     

Membrane-binding mechanism of the EEA1 FYVE domain revealed by multi-scale molecular dynamics simulations

Early Endosomal Antigen 1 (EEA1) is a key protein in endosomal trafficking and is implicated in both autoimmune and neurological diseases. The C-terminal FYVE domain of EEA1 binds endosomal membranes, which contain phosphatidylinositol-3-phosphate (PI(3)P). Although it is known that FYVE binds PI(3)P specifically, it has not previously been described of how FYVE attaches and binds to endosomal membranes. In this study, we employed both coarse-grained (CG) and atomistic (AT) molecular dynamics (MD) simulations to determine how FYVE binds to PI(3)P-containing membranes. CG-MD showed that the dominant membrane binding mode resembles the crystal structure of EEA1 FYVE domain in complex with inositol-1,3-diphospate (PDB ID 1JOC). FYVE, which is a homodimer, binds the membrane via a hinge mechanism, where the C-terminus of one monomer first attaches to the membrane, followed by the C-terminus of the other monomer. The estimated total binding energy is ~70 kJ/mol, of which 50–60 kJ/mol stems from specific PI(3)P-interactions. By AT-MD, we could partition the binding mode into two types: (i) adhesion by electrostatic FYVE-PI(3)P interaction, and (ii) insertion of amphipathic loops. The AT simulations also demonstrated flexibility within the FYVE homodimer between the C-terminal heads and coiled-coil stem. This leads to a dynamic model whereby the 200 nm long coiled coil attached to the FYVE domain dimer can amplify local hinge-bending motions such that the Rab5-binding domain at the other end of the coiled coil can explore an area of 0.1 μm² in the search for a second endosome with which to interact.     

Solution conformations of a trimannoside from nuclear magnetic resonance and molecular dynamics simulations

N-linked oligosaccharides often act as ligands for receptor proteins in a variety of cell recognition processes. Knowledge of the solution conformations, as well as protein-bound conformations, of these oligosaccharides is required to understand these important interactions. In this paper we present a model for the solution conformations sampled by a simple trimannoside, methyl 3, 6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranoside, which contains two of the most commonly found glycosidic linkages in N-linked oligosaccharides. This model was derived from simulated annealing protocols incorporating distance restraints extracted from NOESY spectra along with torsional restraints computed from three-bond (1)H-(13)C coupling constants measured across the glycosidic bonds. The model was refined in light of unrestrained molecular dynamics simulations conducted in the presence of solvent water. The resulting model depicts a molecule undergoing conformational averaging in solution, adopting four major and two minor conformations. The four major conformations arise from a pair of two-state transitions, one each at the alpha(1-->3) and alpha(1-->6) linkages, whereas the minor conformations result from an additional transition of the alpha(1-->6) linkage. Our data also suggest that the alpha(1-->3) transition is fast and changes the molecular shape slightly, whereas the alpha(1-->6) is much slower and alters the molecular shape dramatically.     

Mechanism of CDK5 activation revealed by steered molecular dynamics simulations and energy calculations

In the current work, CDK5/p25 complexes were pulled apart by applying external forces with steered molecular dynamics (SMD) simulations. The crucial interactions between the kinase and the activation protein were investigated and the SMD simulations showed that several activation-relevant motifs of CDK5 leave p25 in sequence during the pulling and lead to an apo-CDK2 like CDK5 structure after separation. Based on systematic examination of hydrogen bond breaking and classical MD/molecular mechanics-generalized Born/surface area) (MM-GBSA) calculations, a CDK5 activation mechanism by p25 is suggested. This is the first step towards the systemic development of CDK inhibitors and the mechanism proposed could lead to a better understanding of the protein–protein recognition characteristics between the kinase and its activator.     

Structured pathway across the transition state for peptide folding revealed by molecular dynamics simulations

Small globular proteins and peptides commonly exhibit two-state folding kinetics in which the rate limiting step of folding is the surmounting of a single free energy barrier at the transition state (TS) separating the folded and the unfolded states. An intriguing question is whether the polypeptide chain reaches, and leaves, the TS by completely random fluctuations, or whether there is a directed, stepwise process. Here, the folding TS of a 15-residue β-hairpin peptide, Peptide 1, is characterized using independent 2.5 μs-long unbiased atomistic molecular dynamics (MD) simulations (a total of 15 μs). The trajectories were started from fully unfolded structures. Multiple (spontaneous) folding events to the NMR-derived conformation are observed, allowing both structural and dynamical characterization of the folding TS. A common loop-like topology is observed in all the TS structures with native end-to-end and turn contacts, while the central segments of the strands are not in contact. Non-native sidechain contacts are present in the TS between the only tryptophan (W11) and the turn region (P7-G9). Prior to the TS the turn is found to be already locked by the W11 sidechain, while the ends are apart. Once the ends have also come into contact, the TS is reached. Finally, along the reactive folding paths the cooperative loss of the W11 non-native contacts and the formation of the central inter-strand native contacts lead to the peptide rapidly proceeding from the TS to the native state. The present results indicate a directed stepwise process to folding the peptide.     

Fast lipid disorientation at the onset of membrane fusion revealed by molecular dynamics simulations.

Membrane fusion is a key event in vesicular trafficking in every cell, and many fusion-related proteins have been identified. However, how the actual fusion event occurs has not been elucidated. By using molecular dynamics simulations we found that when even a small region of two membranes is closely apposed such that only a limited number of water molecules remain in the apposed area (e.g., by a fusogenic protein and thermal membrane fluctuations), dramatic lipid disorientation results within 100 ps-2 ns, which might initiate membrane fusion. Up to 12% of phospholipid molecules in the apposing layers had their alkyl chains outside the hydrophobic region, lying almost parallel to the membrane surface or protruding out of the bilayer by 2 ns after two membranes were closely apposed.     

Structural differences between allelic variants of the ovine prion protein revealed by molecular dynamics simulations

We have modeled ovine prion protein (residues 119-233) based on NMR structures of PrP from other mammalian species. Modeling of the C-terminal domain of ovine PrP predicts three helices: helix-1 (residues 147-155), flanked by two short beta-strands; helix-2 (residues 176-197), and helix-3 (residues 203-229). Molecular dynamics simulations on this model of ovine PrP have determined structural differences between allelic variants. At neutral pH, limited root mean-squared (RMS) fluctuations were seen in the region of helix-1; between beta-strand-2 and residue 171, and the loop connecting helix-2 and helix-3. At low pH, these RMS fluctuations increased and showed allelic variation. The extent of RMS fluctuation between beta-strand 2 and residue 171 was ARR > ARQ > VRQ. This order was reversed for the loop region connecting helix-2 and helix-3. Although all three variants have the potential to display an extended helix at the C-terminal region of helix-1, the major influence of the VRQ allele was to restrict the conformations of the Asn162 and Arg139 side-chains. Variations observed in the simulations in the vicinity of helix-1 correlated with reactivity of C-terminal specific anti-PrP monoclonal antibodies with peripheral blood cells from scrapie-susceptible and -resistant genotypes of sheep: cells from VRQ homozygous sheep showed uniform reactivity, while cells from ARQ and ARR homozygous sheep showed variable binding. Our data show that molecular dynamics simulations can be used to determine structural differences between allelic variants of ovine PrP. The binding of anti-PrP monoclonal antibodies to ovine blood cells may validate these structural predictions.     

Locally accessible conformations of proteins: multiple molecular dynamics simulations of crambin.

Multiple molecular dynamics (MD) simulations of crambin with different initial atomic velocities are used to sample conformations in the vicinity of the native structure. Individual trajectories of length up to 5 ns sample only a fraction of the conformational distribution generated by ten independent 120 ps trajectories at 300 K. The backbone atom conformational space distribution is analyzed using principal components analysis (PCA). Four different major conformational regions are found. In general, a trajectory samples only one region and few transitions between the regions are observed. Consequently, the averages of structural and dynamic properties over the ten trajectories differ significantly from those obtained from individual trajectories. The nature of the conformational sampling has important consequences for the utilization of MD simulations for a wide range of problems, such as comparisons with X-ray or NMR data. The overall average structure is significantly closer to the X-ray structure than any of the individual trajectory average structures. The high frequency (less than 10 ps) atomic fluctuations from the ten trajectories tend to be similar, but the lower frequency (100 ps) motions are different. To improve conformational sampling in molecular dynamics simulations of proteins, as in nucleic acids, multiple trajectories with different initial conditions should be used rather than a single long trajectory.