Hairpins participating in folding of human telomeric sequence quadruplexes studied by standard and T-REMD simulations

DNA G-hairpins are potential key structures participating in folding of human telomeric guanine quadruplexes (GQ). We examined their properties by standard MD simulations starting from the folded state and long T-REMD starting from the unfolded state, accumulating ∼130 μs of atomistic simulations. Antiparallel G-hairpins should spontaneously form in all stages of the folding to support lateral and diagonal loops, with sub-μs scale rearrangements between them. We found no clear predisposition for direct folding into specific GQ topologies with specific syn/anti patterns. Our key prediction stemming from the T-REMD is that an ideal unfolded ensemble of the full GQ sequence populates all 4096 syn/anti combinations of its four G-stretches. The simulations can propose idealized folding pathways but we explain that such few-state pathways may be misleading. In the context of the available experimental data, the simulations strongly suggest that the GQ folding could be best understood by the kinetic partitioning mechanism with a set of deep competing minima on the folding landscape, with only a small fraction of molecules directly folding to the native fold. The landscape should further include non-specific collapse processes where the molecules move via diffusion and consecutive random rare transitions, which could, e.g. structure the propeller loops.     

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.     

Structural characterization of quadruplex DNA with in-cell EPR approaches

Guanosine-rich DNA sequences have the potential to adopt four-stranded conformations termed quadruplexes. The chromosomes of higher organisms are capped by so-called telomeres that are composed of repeats of the sequence TTAGGG. Up to 200 nucleotides of the G-rich strand form an overhang that is suspected to fold into intramolecular G-quadruplexes. Since induction of quadruplexes at the telomeres results in anti-proliferative effects, the intracellular structure of G-quadruplexes is of high interest as an anti-cancer drug target. Here we give a perspective on the elucidation of DNA sequence folds by electron paramagnetic resonance (EPR) distance measurements. The technique complements X-ray crystallography and NMR spectroscopy, as it can be applied in noncrystalline states, is not intrinsically limited by the size of the bio-macromolecular complex, and is able to analyze flexible structures or coexisting DNA conformation.     

The high-resolution crystal structure of a parallel intermolecular DNA G-4 quadruplex/drug complex employing syn glycosyl linkages

We have determined the X-ray structure of the complex between the DNA quadruplex d(5′-GGGG-3′)₄ and daunomycin, as a potential model for studying drug–telomere interactions. The structure was solved at 1.08 Å by direct methods in space group I4. The asymmetric unit comprises a linear arrangement of one d(GGGG) strand, four daunomycin molecules, a second d(GGGG) strand facing in the opposite direction to the first, and Na and Mg cations. The crystallographic 4-fold axis generates the biological unit, which is a 12-layered structure comprising two sets of four guanine layers, with four layers each of four daunomycins stacked between the 5′ faces of the two quadruplexes. The daunomycin layers fall into two groups which are novel in their mode of self assembly. The only contacts between daunomycin molecules within any one of these layers are van der Waals interactions, however there is substantial π–π stacking between successive daunomycin layers and also with adjacent guanine layers. The structure differs significantly from all other parallel d(TGGGGT)₄ quadruplexes in that the 5′ guanine adopts the unusual syn glycosyl linkage, refuting the widespread belief that such conformations should all be anti. In contrast to the related d(TGGGGT)/daunomycin complex, there are no ligand–quadruplex groove insertion interactions.     

Formation pathways of a guanine-quadruplex DNA revealed by molecular dynamics and thermodynamic analysis of the substates

The formation of a cation-stabilized guanine quadruplex (G-DNA) stem is an exceptionally slow process involving complex kinetics that has not yet been characterized at atomic resolution. Here, we investigate the formation of a parallel stranded G-DNA stem consisting of four strands of d(GGGG) using molecular dynamics simulations with explicit inclusion of counterions and solvent. Due to the limitations imposed by the nanosecond timescale of the simulations, rather than watching for the spontaneous formation of G-DNA, our approach probes the stability of possible supramolecular intermediates (including two-, three-, and four-stranded assemblies with out-of-register base pairing between guanines) on the formation pathway. The simulations suggest that "cross-like" two-stranded assemblies may serve as nucleation centers in the initial formation of parallel stranded G-DNA quadruplexes, proceeding through a series of rearrangements involving trapping of cations, association of additional strands, and progressive slippage of strands toward the full stem. To supplement the analysis, approximate free energies of the models are obtained with explicit consideration of the integral cations. The approach applied here serves as a prototype for qualitatively investigating other G-DNA molecules using molecular dynamics simulation and free-energy analysis.     

Effects of 8-methylguanine on structure, stability and kinetics of formation of tetramolecular quadruplexes

Tetramolecular G-quadruplexes result from the association of four guanine-rich strands. Modification of the backbone strand or the guanine bases of the oligonucleotide may improve stability or introduce new functionalities. In this regard, the 8 position of a guanosine is particularly suitable for introduction of modifications since as it is positioned in the groove of the quadruplex structure. Modifications at this position should not interfere with structural assembly as would changes at Watson-Crick and Hoogsteen sites. In this study, we investigated the effect of an 8-methyl-2′-deoxyguanosine residue (M) on the structure and stability of tetramolecular parallel G-quadruplexes. In some cases, the presence of this residue resulted in the formation of unusual quadruplex structures containing all-syn tetrads. Furthermore, the modified nucleoside M at the 5′-end of the sequence accelerated quadruplex formation by 15-fold or more relative to the unmodified oligonucleotide, which makes this nucleobase an attractive replacement for guanine in the context of tetramolecular parallel quadruplexes.     

A new anti conformation for N-(deoxyguanosin-8-yl)-2-acetylaminofluorene (AAF-dG) allows Watson-Crick pairing in the Sulfolobus solfataricus P2 DNA polymerase IV (Dpo4)

Primer extension studies have shown that the Y-family DNA polymerase IV (Dpo4) from Sulfolobus solfataricus P2 can preferentially insert C opposite N-(deoxyguanosin-8-yl)-2-acetylaminofluorene (AAF-dG) [F. Boudsocq, S. Iwai, F. Hanaoka and R. Woodgate (2001) Nucleic Acids Res., 29, 4607–4616]. Our goal is to elucidate on a structural level how AAF-dG can be harbored in the Dpo4 active site opposite an incoming dCTP, using molecular modeling and molecular dynamics simulations, since AAF-dG prefers the syn glycosidic torsion. Both anti and syn conformations of the templating AAF-dG in a Dpo4 ternary complex were investigated. All four dNTPs were studied. We found that an anti glycosidic torsion with C1′-exo deoxyribose conformation allows AAF-dG to be Watson–Crick hydrogen-bonded with dCTP with modest polymerase perturbation, but other nucleotides are more distorting. The AAF is situated in the Dpo4 major groove open pocket with fluorenyl rings 3′- and acetyl 5′-directed along the modified strand, irrespective of dNTP. With AAF-dG syn, the fluorenyl rings are in the small minor groove pocket and the active site region is highly distorted. The anti-AAF-dG conformation with C1′-exo sugar pucker can explain the preferential incorporation of dC by Dpo4. Possible relevance of our new major groove structure for AAF-dG to other polymerases, lesion repair and solution conformations are discussed.     

Dynameomics: A Consensus View of the Protein Unfolding/Folding Transition State Ensemble across a Diverse Set of Protein Folds

The Dynameomics project aims to simulate a representative sample of all globular protein metafolds under both native and unfolding conditions. We have identified protein unfolding transition state (TS) ensembles from multiple molecular dynamics simulations of high-temperature unfolding in 183 structurally distinct proteins. These data can be used to study individual proteins and individual protein metafolds and to mine for TS structural features common across all proteins. Separating the TS structures into four different fold classes (all proteins, all-α, all-β, and mixed α/β and α + β) resulted in no significant difference in the overall protein properties. The residues with the most contacts in the native state lost the most contacts in the TS ensemble. On average, residues beginning in an α-helix maintained more structure in the TS ensemble than did residues starting in β-strands or any other conformation. The metafolds studied here represent 67% of all known protein structures, and this is, to our knowledge, the largest, most comprehensive study of the protein folding/unfolding TS ensemble to date. One might have expected broad distributions in the average global properties of the TS relative to the native state, indicating variability in the amount of structure present in the TS. Instead, the average global properties converged with low standard deviations across metafolds, suggesting that there are general rules governing the structure and properties of the TS.     

Multistranded DNA structures.

DNA oligonucleotides can form multistranded helices through either the folding of a single strand or the association of two, three or four strands of DNA. Structures of several new DNA triplexes, G-quartet DNA quadruplexes and I-motif DNA quadruplexes have been reported recently. These structures provide new insights into helix stability and folding, loop conformations and cation interactions.     

Strand directionality affects cation binding and movement within tetramolecular G-quadruplexes

Nuclear magnetic resonance study of G-quadruplex structures formed by d(TG₃T) and its modified analogs containing a 5′-5′ or 3′-3′ inversion of polarity sites, namely d(3′TG5′-5′G₂T3′), d(3′T5′-5′G₃T3′) and d(5′TG3′-3′G₂T5’) demonstrates formation of G-quadruplex structures with tetrameric topology and distinct cation-binding preferences. All oligonucleotides are able to form quadruplex structures with two binding sites, although the modified oligonucleotides also form, in variable amounts, quadruplex structures with only one bound cation. The inter-quartet cavities at the inversion of polarity sites bind ammonium ions less tightly than a naturally occurring 5′-3′ backbone. Exchange of ¹⁵ ions between G-quadruplex and bulk solution is faster at the 3′-end in comparison to the 5′-end. In addition to strand directionality, cation movement is influenced by formation of an all-syn G-quartet. Formation of such quartet has been observed also for the parent d(TG₃T) that besides the canonical quadruplex with only all-anti G-quartets, forms a tetramolecular parallel quadruplex containing one all-syn G-quartet, never observed before in unmodified quadruplex structures.     

On the importance of orientation in general base catalysis by carboxylate

The carboxylate group serves as a general-base catalyst in numerous chemical and enzymatic reactions. The importance of the direction in which the proton is transferred to the carboxylate is discussed. syn (on the same side of the CO bond as the forming CO) protonation is estimated to be 10⁴-fold more favorable than anti protonation. To date, in the models studied where intramolecular reactions involve carboxylate, only anti protonation can occur. In these intramolecular models the catalytic efficiency of a carboxylate group may be underestimated due to this inability to achieve an optimal orientation for protonation; whereas for enzymatic reactions involving carboxylate side chains, structural studies support mechanisms involving syn protonation.     

(Un)Folding Mechanisms of the FBP28 WW Domain in Explicit Solvent Revealed by Multiple Rare Event Simulation Methods

We report a numerical study of the (un)folding routes of the truncated FBP28 WW domain at ambient conditions using a combination of four advanced rare event molecular simulation techniques. We explore the free energy landscape of the native state, the unfolded state, and possible intermediates, with replica exchange molecular dynamics. Subsequent application of bias-exchange metadynamics yields three tentative unfolding pathways at room temperature. Using these paths to initiate a transition path sampling simulation reveals the existence of two major folding routes, differing in the formation order of the two main hairpins, and in hydrophobic side-chain interactions. Having established that the hairpin strand separation distances can act as reasonable reaction coordinates, we employ metadynamics to compute the unfolding barriers and find that the barrier with the lowest free energy corresponds with the most likely pathway found by transition path sampling. The unfolding barrier at 300 K is ∼17 k_BT ≈ 42 kJ/mol, in agreement with the experimental unfolding rate constant. This work shows that combining several powerful simulation techniques provides a more complete understanding of the kinetic mechanism of protein folding.     

Conformations of cyclic 3',5'-nucleotides. Effect of the base on the synanti conformer distribution

The furanose and the phosphate rings of cyclic 3′,5′-nucleotides are locked in the ₄T³ and chair conformations respectively. The only variable which shows major conformational flexibility in these molecules is the rotation about the glycosyl bond which describes the orientation of the base relative to the sugar-phosphate bicyclic system. The glycosyl torsion angle has been analyzed for cyclic nucleotides with different purine and pyrimidine bases by use of conformational energy calculations. The results indicate that all the pyrimidine bases, U, T and C show a very strong energetic preference for the anti range of conformations. The calculations predict that among cyclic 3′,5′-purine nucleotides cyclic GMP and cyclic IMP favor the syn conformation to the anti by 95:5 and 70:30 respectively, while cyclic AMP shows a preference for the anti conformation to syn by 70:30. Thus the purines show a greater probability for the syn conformation than the pyrimidines in cyclic 3′,5′-nucleotides.     

A detailed unfolding pathway of a (beta/alpha)8-barrel protein as studied by molecular dynamics simulations

The (beta/alpha)(8)-barrel is the most common protein fold. Similar structural properties for folding intermediates of (beta/alpha)(8)-barrel proteins involved in tryptophan biosynthesis have been reported in a number of experimental studies; these intermediates have the last two beta-strands and three alpha-helices partially unfolded, with other regions of the polypeptide chain native-like in conformation. To investigate the detailed folding/unfolding pathways of these (beta/alpha)(8)-barrel proteins, temperature-induced unfolding simulations of N-(5'-phosphoribosyl)anthranilate isomerase from Escherichia coli were carried out using a special-purpose parallel computer system. Unfolding simulations at five different temperatures showed a sequential unfolding pathway comprised of several events. Early events in unfolding involved disruption of the last two strands and three helices, producing an intermediate ensemble similar to those detected in experimental studies. Then, denaturation of the first two betaalpha units and separation of the sixth strand from the fifth took place independently. The remaining central betaalphabetaalphabeta module persisted the longest during all simulations, suggesting an important role for this module as the incipient folding scaffold. Our simulations also predicted the presence of a nucleation site, onto which several hydrophobic residues condensed forming the foundation for the central betaalphabetaalphabeta module.     

Conformationally rigid nucleoside probes help understand the role of sugar pucker and nucleobase orientation in the thrombin-binding aptamer

Modified thrombin-binding aptamers carrying 2′-deoxyguanine (dG) residues with locked North- or South-bicyclo[3.1.0]hexane pseudosugars were synthesized. Individual 2′-deoxyguanosines at positions dG5, dG10, dG14 and dG15 of the aptamer were replaced by these analogues where the North/anti and South/syn conformational states were confined. It was found that the global structure of the DNA aptamer was, for the most part, very accommodating. The substitution at positions 5, 10 and 14 with a locked South/syn-dG nucleoside produced aptamers with the same stability and global structure as the innate, unmodified one. Replacing position 15 with the same South/syn-dG nucleoside induced a strong destabilization of the aptamer, while the antipodal North/anti-dG nucleoside was less destabilizing. Remarkably, the insertion of a North/anti-dG nucleoside at position 14, where both pseudosugar conformation and glycosyl torsion angle are opposite with respect to the native structure, led to the complete disruption of the G-tetraplex structure as detected by NMR and confirmed by extensive molecular dynamics simulations. We conclude that conformationally locked bicyclo[3.1.0]hexane nucleosides appear to be excellent tools for studying the role of key conformational parameters that are critical for the formation of a stable, antiparallel G-tetrad DNA structures.     

Influence of local sequence context on damaged base conformation in human DNA polymerase ι: molecular dynamics studies of nucleotide incorporation opposite a benzo[a]pyrene-derived adenine lesion

Human DNA polymerase ι is a lesion bypass polymerase of the Y family, capable of incorporating nucleotides opposite a variety of lesions in both near error-free and error-prone bypass. With undamaged templating purines polymerase ι normally favors Hoogsteen base pairing. Polymerase ι can incorporate nucleotides opposite a benzo[a]pyrene-derived adenine lesion (dA*); while mainly error-free, the identity of misincorporated bases is influenced by local sequence context. We performed molecular modeling and molecular dynamics simulations to elucidate the structural basis for lesion bypass. Our results suggest that hydrogen bonds between the benzo[a]pyrenyl moiety and nearby bases limit the movement of the templating base to maintain the anti glycosidic bond conformation in the binary complex in a 5′-CAGA*TT-3′ sequence. This facilitates correct incorporation of dT via a Watson−Crick pair. In a 5′-TTTA*GA-3′ sequence the lesion does not form these hydrogen bonds, permitting dA* to rotate around the glycosidic bond to syn and incorporate dT via a Hoogsteen pair. With syn dA*, there is also an opportunity for increased misincorporation of dGTP. These results expand our understanding of the versatility and flexibility of polymerase ι and its lesion bypass functions in humans.     

Molecular Origin of the Stability Difference in Four Shark IgNAR Constant Domains

Improving the stability of antibodies for manufacture and shelf life is one of the main focuses of antibody engineering. One stabilization strategy is to perform specific mutations in human antibodies based on highly stable antibodies in other species. To identify the key residues for mutagenesis, it is necessary to understand the roles of these residues in stabilizing the antibody. Here, we use molecular dynamics simulations to study the molecular origin of the four shark immunoglobulin new antigen receptors constant domains (C1–C4). According to the unfolding pathways and the conformational free energy surfaces in 8 M urea at 380 K, the C2 domain is the most stable, followed by C4, C1, and C3, which agrees with the experimental findings. The C1 and C3 domains follow a common unfolding pathway in which the unfolding starts from the edge strands, particularly strand g, and then gradually progresses to the inner strands. Detailed structural analysis of the C2 domain reveals a “sandwich-like” R339-E322-R341 salt-bridge cluster on strand g, which grants ultrahigh stability to the C2 domain. We further design two sets of mutations by mutating E322 to alanine or setting all atomic charges in E322 to zero to break the salt-bridge cluster in the C2 domain, which confirms the importance of the salt bridges in stability. In the C4 domain, the D80-K104 salt bridge on strand g also strengthens the stability. On the other hand, in the C1 and C3 domains, there is no salt bridge on strand g. In addition to the salt bridges, the overall hydrophobicity score of the hydrophobic core is also positively correlated with the domain stability. Our findings provide a detailed microscopic picture of the molecular origin of the four shark immunoglobulin new antigen receptors constant domains that not only explains the differences in their structural stability but also provides important insights into future antibody design.     

The beginning and the end: flanking nucleotides induce a parallel G-quadruplex topology

Genomic sequences susceptible to form G-quadruplexes (G4s) are always flanked by other nucleotides, but G4 formation in vitro is generally studied with short synthetic DNA or RNA oligonucleotides, for which bases adjacent to the G4 core are often omitted. Herein, we systematically studied the effects of flanking nucleotides on structural polymorphism of 371 different oligodeoxynucleotides that adopt intramolecular G4 structures. We found out that the addition of nucleotides favors the formation of a parallel fold, defined as the ‘flanking effect’ in this work. This ‘flanking effect’ was more pronounced when nucleotides were added at the 5′-end, and depended on loop arrangement. NMR experiments and molecular dynamics simulations revealed that flanking sequences at the 5′-end abolish a strong syn-specific hydrogen bond commonly found in non-parallel conformations, thus favoring a parallel topology. These analyses pave a new way for more accurate prediction of DNA G4 folding in a physiological context.