Molecular dynamics of DNA quadruplex molecules containing inosine, 6-thioguanine and 6-thiopurine

The ability of the four-stranded guanine (G)-DNA motif to incorporate nonstandard guanine analogue bases 6-oxopurine (inosine, I), 6-thioguanine (tG), and 6-thiopurine (tI) has been investigated using large-scale molecular dynamics simulations. The simulations suggest that a G-DNA stem can incorporate inosines without any marked effect on its structure and dynamics. The all-inosine quadruplex stem d(IIII)(4) shows identical dynamical properties as d(GGGG)(4) on the nanosecond time scale, with both molecular assemblies being stabilized by monovalent cations residing in the channel of the stem. However, simulations carried out in the absence of these cations show dramatic differences in the behavior of d(GGGG)(4) and d(IIII)(4). Whereas vacant d(GGGG)(4) shows large fluctuations but does not disintegrate, vacant d(IIII)(4) is completely disrupted within the first nanosecond. This is a consequence of the lack of the H-bonds involving the N2 amino group that is not present in inosine. This indicates that formation of the inosine quadruplex could involve entirely different intermediate structures than formation of the guanosine quadruplex, and early association of cations in this process appears to be inevitable. In the simulations, the incorporation of 6-thioguanine and 6-thiopurine sharply destabilizes four-stranded G-DNA structures, in close agreement with experimental data. The main reason is the size of the thiogroup leading to considerable steric conflicts and expelling the cations out of the channel of the quadruplex stem. The G-DNA stem can accommodate a single thioguanine base with minor perturbations. Incorporation of a thioguanine quartet layer is associated with a large destabilization of the G-DNA stem whereas the all-thioguanine quadruplex immediately collapses.     

Structural dynamics of possible late-stage intermediates in folding of quadruplex DNA studied by molecular simulations

Explicit solvent molecular dynamics simulations have been used to complement preceding experimental and computational studies of folding of guanine quadruplexes (G-DNA). We initiate early stages of unfolding of several G-DNAs by simulating them under no-salt conditions and then try to fold them back using standard excess salt simulations. There is a significant difference between G-DNAs with all-anti parallel stranded stems and those with stems containing mixtures of syn and anti guanosines. The most natural rearrangement for all-anti stems is a vertical mutual slippage of the strands. This leads to stems with reduced numbers of tetrads during unfolding and a reduction of strand slippage during refolding. The presence of syn nucleotides prevents mutual strand slippage; therefore, the antiparallel and hybrid quadruplexes initiate unfolding via separation of the individual strands. The simulations confirm the capability of G-DNA molecules to adopt numerous stable locally and globally misfolded structures. The key point for a proper individual folding attempt appears to be correct prior distribution of syn and anti nucleotides in all four G-strands. The results suggest that at the level of individual molecules, G-DNA folding is an extremely multi-pathway process that is slowed by numerous misfolding arrangements stabilized on highly variable timescales.     

Molecular dynamics simulations of G-DNA and perspectives on the simulation of nucleic acid structures

The article reviews the application of biomolecular simulation methods to understand the structure, dynamics and interactions of nucleic acids with a focus on explicit solvent molecular dynamics simulations of guanine quadruplex (G-DNA and G-RNA) molecules. While primarily dealing with these exciting and highly relevant four-stranded systems, where recent and past simulations have provided several interesting results and novel insight into G-DNA structure, the review provides some general perspectives on the applicability of the simulation techniques to nucleic acids.     

Cation binding to 15-TBA quadruplex DNA is a multiple-pathway cation-dependent process

A combination of explicit solvent molecular dynamics simulation (30 simulations reaching 4 μs in total), hybrid quantum mechanics/molecular mechanics approach and isothermal titration calorimetry was used to investigate the atomistic picture of ion binding to 15-mer thrombin-binding quadruplex DNA (G-DNA) aptamer. Binding of ions to G-DNA is complex multiple pathway process, which is strongly affected by the type of the cation. The individual ion-binding events are substantially modulated by the connecting loops of the aptamer, which play several roles. They stabilize the molecule during time periods when the bound ions are not present, they modulate the route of the ion into the stem and they also stabilize the internal ions by closing the gates through which the ions enter the quadruplex. Using our extensive simulations, we for the first time observed full spontaneous exchange of internal cation between quadruplex molecule and bulk solvent at atomistic resolution. The simulation suggests that expulsion of the internally bound ion is correlated with initial binding of the incoming ion. The incoming ion then readily replaces the bound ion while minimizing any destabilization of the solute molecule during the exchange.     

The stability of cylindrin β‐barrel amyloid oligomer models—A molecular dynamics study

Small‐soluble amyloid oligomers are believed to play a significant role in the pathology of amyloid diseases. Recently, the atomic structure of a toxic oligomer formed by an 11 residue and its tandem repeat was found to have an out‐off register antiparallel β‐strands in the shape of a β‐barrel. In the present article we investigate the effect of mutations in the hydrophobic cores on the structure and dynamic of the β‐barrels using all atom multiple molecular dynamics simulations with an explicit solvent. Extending previous experiments with molecular dynamics simulations we systematically test how stability and formation of cylindrin depends on the interplay between hydrophobicity and steric effects of the core residues. We find that strong hydrophobic interactions between geometrically fitting residues keep the strands (both in register and out‐off‐register interface) in close proximity, which in turn stabilizes the side‐chain and main‐chain hydrogen bonds, and the salt bridges on the outer surface along the weak out‐of‐register interface. Our simulations also indicate presence of water molecules in the hydrophobic interior of the cylindrin β‐barrel.Proteins 2013. © 2013 Wiley Periodicals, Inc.     

Mechanisms of the self-assembly of EAK16-family peptides into fibrillar and globular structures: molecular dynamics simulations from nano- to micro-seconds

The self-assembly of EAK16-family peptides in a bulk solution was studied using a combination of all-atom and coarse-grained molecular dynamics simulations. In addition, specified concentrations of EAK16 peptides were induced to form fibrillary or globular assemblies in vitro. The results show that the combination of all-atom molecular dynamics simulations on the single- and double-chain levels and coarse-grained simulations on the many-chain level predicts the experimental observations reasonably well. At neutral pH conditions, EAK16-I and EAK16-II assemble into fibrillary structures, whereas EAK16-IV aggregates into globular assemblies. Mechanisms of the formation of fibrillar and globular assemblies are described using the simulation results.     

Refolding molecular dynamics simulations of small- and middle-sized proteins in an explicit solvent

In order to elucidate the protein folding problem, we performed molecular dynamics simulations for small- and middle-sized two unfolding and six refolding proteins in an explicit solvent. Histidine-containing phosphocarrier protein and small designed protein were chosen for the simulations. We found that the protein folding process of these proteins was divided into three phases: an α -helix formation phase, a packing phase and a β -sheet formation phase. In the α -helix formation phase, an α -helix was developed from a β -turn structure through a 3₁₀-helix state. In the packing phase, proteins became compact, and tertiary structures (α / α or pre- β / β packing) were formed. Formation of a hydrophobic nucleus occurred concomitant with the α -helix formation and packing phase. Finally, in the β -sheet formation phase, a β -sheet was developed owing to the sequential formation of hydrogen bonds between two neighbouring strands, just like a "closing zipper".     

Stability of Aβ‐fibril fragments in the presence of fatty acids

We consider the effect of lauric acid on the stability of various fibril‐like assemblies of Aβ peptides. For this purpose, we have performed molecular dynamics simulations of these assemblies either in complex with lauric acid or without presence of the ligand. While we do not observe a stabilizing effect on Aβ₄₀‐fibrils, we find that addition of lauric acid strengthens the stability of fibrils built from the triple‐stranded S‐shaped Aβ₄₂‐peptides considered to be more toxic. Or results may help to understand how the specifics of the brain‐environment modulate amyloid formation and propagation.     

Association of DNA quadruplexes through G:C:G:C tetrads. Solution structure of d(GCGGTGGAT)

The structure formed by the DNA sequence d(GCGGTGGAT) in a 100 mM Na(+) solution has been determined using molecular dynamics calculations constrained by distance and dihedral restraints derived from NMR experiments performed at isotopic natural abundance. The sequence folds into a dimer of dimers. Each symmetry-related half contains two parallel stranded G:G:G:G tetrads flanked by an A:A mismatch and by four-stranded G:C:G:C tetrads. Each of the two juxtaposed G:C:G:C tetrads is composed of alternating antiparallel strands from the two halves of the dimer. For each single strand, a thymine intersperses a double chain reversal connecting the juxtaposed G:G:G:G tetrads. This architecture has potential implications in genetic recombination. It suggests a pathway for oligomerization involving association of quadruplex entities through GpC steps.     

Monovalent cation induced structural transitions in telomeric DNAs: G-DNA folding intermediates

Telomeric DNA consists of G- and C-rich strands that are always polarized such that the G-rich strand extends past the 3' end of the duplex to form a 12-16-base overhang. These overhanging strands can self-associate in vitro to form intramolecular structures that have several unusual physical properties and at least one common feature, the presence of non-Watson-Crick G.G base pairs. The term "G-DNA" was coined for this class of structures (Cech, 1988). On the basis of gel electrophoresis, imino proton NMR, and circular dichroism (CD) results, we find that changing the counterions from sodium to potassium (in 20 mM phosphate buffers) specifically induces conformational transitions in the G-rich telomeric DNA from Tetrahymena, d(T2G4)4 (TET4), which results in a change from the intramolecular species to an apparent multistranded structure, accompanied by an increase in the melting temperature of the base pairs of greater than 25 degrees, as monitored by loss of the imino proton NMR signals. NMR semiselective spin-lattice relaxation rate measurements and HPLC size-exclusion chromatography studies show that in 20 mM potassium phosphate (pH 7) buffer (KP) TET4 is approximately twice the length of the form obtained in 20 mM sodium phosphate (pH 7) buffer (NaP) and that mixtures of Na+ and K+ produce mixtures of the two forms whose populations depend on the ratio of the cations. Since K+ and NH4+ are known to stabilize a parallel-stranded quadruplex structure of poly[r(I)4], we infer that the multistranded structure is a quadruplex. Our results indicate that specific differences in ionic interactions can result in a switch in telomeric DNAs between intramolecular hairpin-like or quadruplex-containing species and intermolecular quadruplex structures, all of which involve G.G base pairing interactions. We propose a model in which duplex or hairpin forms of G-DNA are folding intermediates in the formation of either 1-, 2-, or 4-stranded quadruplex structures. In this model monovalent cations stabilize the duplex and quadruplex forms via two distinct mechanisms, counterion condensation and octahedral coordination to the carbonyl groups in stacked planar guanine "quartet" base assemblies. Substituting one of the guanosine residues in each of the repeats of the Tetrahymena sequence to give the human telomeric DNA, d(T2AG3)4, results in less effective K(+)-dependent stabilization. Thus, the ion-dependent stabilization is attenuated by altering the sequence. Upon addition of the Watson-Crick (WC) complementary strand, only the Na(+)-stabilized structure dissociates quickly to form a WC double helix.(ABSTRACT TRUNCATED AT 400 WORDS)     

Molecular dynamics simulations of Guanine quadruplex loops: advances and force field limitations

A computational analysis of d(GGGGTTTTGGGG)(2) guanine quadruplexes containing either lateral or diagonal four-thymidine loops was carried out using molecular dynamics (MD) simulations in explicit solvent, locally enhanced sampling (LES) simulations, systematic conformational search, and free energy molecular-mechanics, Poisson Boltzmann, surface area (MM-PBSA) calculations with explicit inclusion of structural monovalent cations. The study provides, within the approximations of the applied all-atom additive force field, a qualitatively complete analysis of the available loop conformational space. The results are independent of the starting structures. Major conformational transitions not seen in conventional MD simulations are observed when LES is applied. The favored LES structures consistently provide lower free energies (as estimated by molecular-mechanics, Poisson Boltzmann, surface area) than other structures. Unfortunately, the predicted optimal structure for the diagonal loop arrangement differs substantially from the atomic resolution experiments. This result is attributed to force field deficiencies, such as the potential misbalance between solute-cation and solvent-cation terms. The MD simulations are unable to maintain the stable coordination of the monovalent cations inside the diagonal loops as reported in a recent x-ray study. The optimal diagonal and lateral loop arrangements appear to be close in energy although a proper inclusion of the loop monovalent cations could stabilize the diagonal architecture.     

Conformational behavior of nucleotide-sugar in solution: molecular dynamics and NMR study of solvated uridine diphosphate-glucose in the presence of monovalent cations

The nucleotide-sugars are metabolites of primary importance in the biosynthesis of polysaccharides and glycoconjugates since they serve as sugar donors in the reactions of glycosyltransferases, enzymes that displays a high specificity for both donors and acceptors. In order to determine the conformational behavior of uridinediphosphoglucose in dilute aqueous solution that includes a physiologically relevant concentration of salt, parallel NMR and molecular modeling investigations have been conducted. Nine molecular dynamics trajectories of 3 ns each were calculated in presence of explicit water and monovalent cations with the use of the AMBER force field with recently developed energy parameters for nucleotide-sugars (P. Petrova, J. Koca, and A. Imberty, Journal of American Chemical Society, 1999, vol. 121, pp. 5535-5547). Theoretical nuclear Overhauser effect data were calculated from these simulations using a model-free approach that takes into account internal motions. Comparison of theoretical and experimental data gives excellent agreement for the region surrounding the glucose-phosphate linkage including the pyrophosphate linkage itself. Less satisfactory agreement is obtained for the ribose ring and the base orientations. On the whole, both NMR and molecular dynamics simulations predict the molecule to be flexible, and to visit a large number of conformations while maintaining an extended overall shape.     

Molecular dynamics simulations of the unfolding of beta(2)-microglobulin and its variants

In this study, we examined the unfolding processes of native beta(2)-microglobulin and two related variants, one with an N-terminal hexapeptide deletion DeltaN6 and another with Lys57-Asp58 cleavage, by high-temperature molecular dynamics simulations. Three simulation models were used, molecular dynamics (MD) simulations with explicit water solvation, MD simulations with the CHARMM EEF1 force field and Langevin dynamics with the CHARMM EEF1 force field. Our simulations reproduce many of the experimentally observed structural changes. The most striking agreement is in the beta-strands to alpha-helix transition. In our simulations, strands beta(3), beta(4) and beta(5) consistently change to alpha-helix, whereas beta(8) changes to an alpha-helix only briefly. Through comparisons of the conformational behavior of the native, the DeltaN6 and the Lys-cut beta(2)-m, using the three simulation methods, we identified the consensus conformational changes that differentiate between the native beta(2)-m and its two variants. We found that the main effect of the removal of the N-terminal hexapeptide is to increase the separation between strands beta(2) and beta(6) and to facilitate the beta to alpha transition. On the other hand, the lysine cleavage only increases the flexibility of strand beta(5) and does not affect the interactions between strands beta(2) and beta(6). These conformational changes may relate to polymerization tendencies of these variants.     

Molecular dynamics simulations on the oligomer-formation process of the GNNQQNY peptide from yeast prion protein Sup35

Oligomeric intermediates are possible cytotoxic species in diseases associated with amyloid deposits. Understanding the early steps of fibril formation at atomic details may provide useful information for the rational therapeutic design. In this study, using the heptapeptide GNNQQNY from the yeast prion-like protein Sup35 as a model system, for which a detailed atomic structure of the fibril formed has been determined by x-ray microcrystallography, we investigated its oligomer-formation process from monomer to tetramer at the atomistic level by means of a molecular dynamics simulation with explicit water. Although the number of simulations was limited, the qualitative statistical data gave some interesting results, which indicated that the oligomer formation might start from antiparallel beta-sheet-like dimers. When a new single peptide strand was added to the preformed dimers to form trimers and then tetramers, the transition time from disorder aggregates to regular ones for the parallel alignment was found to be obviously much less than for the antiparallel one. Moreover, the parallel pattern also statistically stayed longer, providing more chances for oligomer extending, although the number of parallel stack events was almost equal to antiparallel ones. Therefore, our simulations showed that new strands might prefer to extend in a parallel arrangement to form oligomers, which agrees with the microcrystal structure of the amyloid fibril formed by this peptide. In addition, analysis of the pi-pi stacking of aromatic residues showed that this type of interaction did not play an important role in giving directionality for beta-strand alignment but played a great influence on stabilizing the structures formed in the oligomer-formation process.     

Atomistic Picture for the Folding Pathway of a Hybrid-1 Type Human Telomeric DNA G-quadruplex

In this work we studied the folding process of the hybrid-1 type human telomeric DNA G-quadruplex with solvent and ions explicitly modeled. Enabled by the powerful bias-exchange metadynamics and large-scale conventional molecular dynamic simulations, the free energy landscape of this G-DNA was obtained for the first time and four folding intermediates were identified, including a triplex and a basically formed quadruplex. The simulations also provided atomistic pictures for the structures and cation binding patterns of the intermediates. The results showed that the structure formation and cation binding are cooperative and mutually supporting each other. The syn/anti reorientation dynamics of the intermediates was also investigated. It was found that the nucleotides usually take correct syn/anti configurations when they form native and stable hydrogen bonds with the others, while fluctuating between two configurations when they do not. Misfolded intermediates with wrong syn/anti configurations were observed in the early intermediates but not in the later ones. Based on the simulations, we also discussed the roles of the non-native interactions. Besides, the formation process of the parallel conformation in the first two G-repeats and the associated reversal loop were studied. Based on the above results, we proposed a folding pathway for the hybrid-1 type G-quadruplex with atomistic details, which is new and more complete compared with previous ones. The knowledge gained for this type of G-DNA may provide a general insight for the folding of the other G-quadruplexes.     

Multiscale modeling of lipids and lipid bilayers

A multiscale modeling approach is applied for simulations of lipids and lipid assemblies on mesoscale. First, molecular dynamics simulation of initially disordered system of lipid molecules in water within all-atomic model was carried out. On the next stage, structural data obtained from the molecular dynamics (MD) simulation were used to build a coarse-grained (ten sites) lipid model, with effective interaction potentials computed by the inverse Monte Carlo method. Finally, several simulations of the coarse-grained model on longer length- and time-scale were performed, both within Monte Carlo and molecular dynamics simulations: a periodical sample of lipid molecules ordered in bilayer, a free sheet of such bilayer without periodic boundary conditions, formation of vesicle from a plain membrane, process of self-assembly of lipids randomly dispersed in volume. It was shown that the coarse-grained model, developed exclusively from all-atomic simulation data, reproduces well all the basic features of lipids in water solution.     

Induced β-Barrel Formation of the Alzheimer's Aβ25-35 Oligomers on Carbon Nanotube Surfaces: Implication for Amyloid Fibril Inhibition

Recent experimental studies show that carbon nanotubes impact the aggregation process of proteins associated with neurodegenerative diseases. However, the details of molecular interactions between proteins and carbon nanotubes are still not well understood. In this study, we investigate the initial adsorption features and dynamics of the Alzheimer's amyloid-β peptide spanning residues 25-35 (Aβ25-35) on a single-walled carbon nanotube (SWNT) surface using fully atomic molecular dynamics simulations (MD) in explicit solvent. The initial configurations of the Aβ25-35 peptides consist of two preformed bilayer β-sheets, each with four or five β-strands in parallel or mixed antiparallel-parallel orientations. Our simulations show, for what we believe is the first time, that two disjointed Aβ25-35 β-sheets with mixed antiparallel-parallel strands can assemble into β-barrels wrapping the SWNT. In contrast, both simulations of Aβ25-35 without SWNT, and simulations of SWNT−Aβ25-35 with purely parallel β-strands, lead to disordered aggregates. We find that Aβ25-35 β-barrel formation involves at least two steps: i), curving of the Aβ25-35 β-sheets as a result of strong hydrophobic interactions with carbon nanotube concomitantly with dehydration of the SWNT-peptide interface; and ii), intersheet backbone hydrogen bond formation with fluctuating intrasheet hydrogen bonds. Detailed analysis of the conversion shows that β-barrel formation on SWNT surface results from the interplay of dehydration and peptide-SWNT/peptide-peptide interactions. Implications of our results on amyloid fibril inhibition are discussed.