We have carried out a series of extended unbiased molecular dynamics (MD) simulations (up to 10 μs long, ∼162 μs in total) complemented by replica-exchange with the collective variable tempering (RECT) approach for several human telomeric DNA G-quadruplex (GQ) topologies with TTA propeller loops. We used different AMBER DNA force-field variants and also processed simulations by Markov State Model (MSM) analysis. The slow conformational transitions in the propeller loops took place on a scale of a few μs, emphasizing the need for long simulations in studies of GQ dynamics. The propeller loops sampled similar ensembles for all GQ topologies and for all force-field dihedral-potential variants. The outcomes of standard and RECT simulations were consistent and captured similar spectrum of loop conformations. However, the most common crystallographic loop conformation was very unstable with all force-field versions. Although the loss of canonical γ-trans state of the first propeller loop nucleotide could be related to the indispensable bsc0 α/γ dihedral potential, even supporting this particular dihedral by a bias was insufficient to populate the experimentally dominant loop conformation. In conclusion, while our simulations were capable of providing a reasonable albeit not converged sampling of the TTA propeller loop conformational space, the force-field description still remained far from satisfactory.

• MeSH
• DNA chemistry   genetics   metabolism   MeSH
• G-Quadruplexes * MeSH
• Humans MeSH
• Base Sequence MeSH
• Cluster Analysis MeSH
• Molecular Dynamics Simulation * MeSH
• Telomere genetics   MeSH
• Water chemistry   MeSH
• Check Tag
• Humans MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• DNA MeSH
• Water MeSH

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.

• MeSH
• DNA chemistry   MeSH
• G-Quadruplexes * MeSH
• Cations chemistry   MeSH
• Humans MeSH
• Oxytricha genetics   MeSH
• Molecular Dynamics Simulation * MeSH
• Telomere chemistry   MeSH
• Check Tag
• Humans MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• DNA MeSH
• Cations MeSH

Abasic (AP) lesions are the most frequent type of damages occurring in cellular DNA. Here we describe the conformational effects of AP sites substituted for 2'-deoxyadenosine in the first (ap7), second (ap13) or third (ap19) loop of the quadruplex formed in K(+) by the human telomere DNA 5'-d[AG3(TTAG3)3]. CD spectra and electrophoresis reveal that the presence of AP sites does not hinder the formation of intramolecular quadruplexes. NMR spectra show that the structural heterogeneity is substantially reduced in ap7 and ap19 as compared to that in the wild-type. These two (ap7 and ap19) sequences are shown to adopt the hybrid-1 and hybrid-2 quadruplex topology, respectively, with AP site located in a propeller-like loop. All three studied sequences transform easily into parallel quadruplex in dehydrating ethanol solution. Thus, the AP site in any loop region facilitates the formation of the propeller loop. Substitution of all adenines by AP sites stabilizes the parallel quadruplex even in the absence of ethanol. Whereas guanines are the major determinants of quadruplex stability, the presence or absence of loop adenines substantially influences quadruplex folding. The naturally occurring adenine-lacking sites in the human telomere DNA can change the quadruplex topology in vivo with potentially vital biological consequences.

• MeSH
• Adenine chemistry   MeSH
• Potassium chemistry   MeSH
• G-Quadruplexes * MeSH
• Guanine chemistry   MeSH
• Humans MeSH
• Models, Molecular MeSH
• Nuclear Magnetic Resonance, Biomolecular MeSH
• DNA Damage * MeSH
• Telomere chemistry   MeSH
• Check Tag
• Humans MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• Adenine MeSH
• Potassium MeSH
• Guanine MeSH

Four-stranded DNA structures were structurally characterized in vitro by NMR, X-ray and Circular Dichroism spectroscopy in detail. Among the different types of quadruplexes (i-Motifs, minor groove quadruplexes, G-quadruplexes, etc.), the best described are G-quadruplexes which are featured by Hoogsteen base-paring. Sequences with the potential to form quadruplexes are widely present in genome of all organisms. They are found often in repetitive sequences such as telomeric ones, and also in promoter regions and 5' non-coding sequences. Recently, many proteins with binding affinity to G-quadruplexes have been identified. One of the initially portrayed G-rich regions, the human telomeric sequence (TTAGGG)n, is recognized by many proteins which can modulate telomerase activity. Sequences with the potential to form G-quadruplexes are often located in promoter regions of various oncogenes. The NHE III1 region of the c-MYC promoter has been shown to interact with nucleolin protein as well as other G-quadruplex-binding proteins. A number of G-rich sequences are also present in promoter region of estrogen receptor alpha. In addition to DNA quadruplexes, RNA quadruplexes, which are critical in translational regulation, have also been predicted and observed. For example, the RNA quadruplex formation in telomere-repeat-containing RNA is involved in interaction with TRF2 (telomere repeat binding factor 2) and plays key role in telomere regulation. All these fundamental examples suggest the importance of quadruplex structures in cell processes and their understanding may provide better insight into aging and disease development.

• MeSH
• DNA-Binding Proteins chemistry   metabolism   MeSH
• DNA chemistry   metabolism   MeSH
• G-Quadruplexes * MeSH
• Nucleic Acid Conformation MeSH
• Humans MeSH
• Promoter Regions, Genetic MeSH
• RNA chemistry   metabolism   MeSH
• Aging MeSH
• Telomere MeSH
• Check Tag
• Humans MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Review MeSH
• Names of Substances
• DNA-Binding Proteins MeSH
• DNA MeSH
• RNA MeSH