Guanine quadruplexes (GQs) play crucial roles in various biological processes, and understanding their folding pathways provides insight into their stability, dynamics, and functions. This knowledge aids in designing therapeutic strategies, as GQs are potential targets for anticancer drugs and other therapeutics. Although experimental and theoretical techniques have provided valuable insights into different stages of the GQ folding, the structural complexity of GQs poses significant challenges, and our understanding remains incomplete. This study introduces a novel computational protocol for folding an entire GQ from single-strand conformation to its native state. By combining two complementary enhanced sampling techniques, we were able to model folding pathways, encompassing a diverse range of intermediates. Although our investigation of the GQ free energy surface (FES) is focused solely on the folding of the all-anti parallel GQ topology, this protocol has the potential to be adapted for the folding of systems with more complex folding landscapes.

• Klíčová slova
• DNA quadruplex, computational folding, enhanced sampling, kinetic partitioning mechanism, metadynamics, molecular dynamics, nudged elastic band, pathCV, transition path sampling,
• MeSH
• DNA chemie   MeSH
• G-kvadruplexy * MeSH
• konformace nukleové kyseliny MeSH
• simulace molekulární dynamiky MeSH
• termodynamika MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• DNA MeSH

Guanine quadruplex (GQ) is a noncanonical nucleic acid structure formed by guanine-rich DNA and RNA sequences. Folding of GQs is a complex process, where several aspects remain elusive, despite being important for understanding structure formation and biological functions of GQs. Pulling experiments are a common tool for acquiring insights into the folding landscape of GQs. Herein, we applied a computational pulling strategy─steered molecular dynamics (SMD) simulations─in combination with standard molecular dynamics (MD) simulations to explore the unfolding landscapes of tetrameric parallel GQs. We identified anisotropic properties of elastic conformational changes, unfolding transitions, and GQ mechanical stabilities. Using a special set of structural parameters, we found that the vertical component of pulling force (perpendicular to the average G-quartet plane) plays a significant role in disrupting GQ structures and weakening their mechanical stabilities. We demonstrated that the magnitude of the vertical force component depends on the pulling anchor positions and the number of G-quartets. Typical unfolding transitions for tetrameric parallel GQs involve base unzipping, opening of the G-stem, strand slippage, and rotation to cross-like structures. The unzipping was detected as the first and dominant unfolding event, and it usually started at the 3'-end. Furthermore, results from both SMD and standard MD simulations indicate that partial spiral conformations serve as a transient ensemble during the (un)folding of GQs.

• MeSH
• biomechanika MeSH
• DNA chemie   MeSH
• G-kvadruplexy * MeSH
• mechanické jevy MeSH
• simulace molekulární dynamiky * MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• DNA MeSH

Guanine quadruplexes (GQs) are non-canonical nucleic acid structures involved in many biological processes. GQs formed in single-stranded regions often need to be unwound by cellular machinery, so their mechanochemical properties are important. Here, we performed steered molecular dynamics simulations of human telomeric GQs to study their unfolding. We examined four pulling regimes, including a very slow setup with pulling velocity and force load accessible to high-speed atomic force microscopy. We identified multiple factors affecting the unfolding mechanism, i.e.,: (i) the more the direction of force was perpendicular to the GQ channel axis (determined by GQ topology), the more the base unzipping mechanism happened, (ii) the more parallel the direction of force was, GQ opening and cross-like GQs were more likely to occur, (iii) strand slippage mechanism was possible for GQs with an all-anti pattern in a strand, and (iv) slower pulling velocity led to richer structural dynamics with sampling of more intermediates and partial refolding events. We also identified that a GQ may eventually unfold after a force drop under forces smaller than those that the GQ withstood before the drop. Finally, we found out that different unfolding intermediates could have very similar chain end-to-end distances, which reveals some limitations of structural interpretations of single-molecule spectroscopic data.

• MeSH
• G-kvadruplexy * MeSH
• guanin * chemie   MeSH
• lidé MeSH
• mechanické jevy MeSH
• simulace molekulární dynamiky MeSH
• telomery MeSH
• Check Tag
• lidé MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• guanin * MeSH

Guanine quadruplexes (G4s) are non-canonical nucleic acids structures common in important genomic regions. Parallel-stranded G4 folds are the most abundant, but their folding mechanism is not fully understood. Recent research highlighted that G4 DNA molecules fold via kinetic partitioning mechanism dominated by competition amongst diverse long-living G4 folds. The role of other intermediate species such as parallel G-triplexes and G-hairpins in the folding process has been a matter of debate. Here, we use standard and enhanced-sampling molecular dynamics simulations (total length of ∼0.9 ms) to study these potential folding intermediates. We suggest that parallel G-triplex per se is rather an unstable species that is in local equilibrium with a broad ensemble of triplex-like structures. The equilibrium is shifted to well-structured G-triplex by stacked aromatic ligand and to a lesser extent by flanking duplexes or nucleotides. Next, we study propeller loop formation in GGGAGGGAGGG, GGGAGGG and GGGTTAGGG sequences. We identify multiple folding pathways from different unfolded and misfolded structures leading towards an ensemble of intermediates called cross-like structures (cross-hairpins), thus providing atomistic level of description of the single-molecule folding events. In summary, the parallel G-triplex is a possible, but not mandatory short-living (transitory) intermediate in the folding of parallel-stranded G4.

• MeSH
• DNA chemie   genetika   metabolismus   MeSH
• G-kvadruplexy * MeSH
• guanin chemie   metabolismus   MeSH
• jednovláknová DNA chemie   genetika   metabolismus   MeSH
• kinetika MeSH
• konformace nukleové kyseliny * MeSH
• lidé MeSH
• sekvence nukleotidů MeSH
• simulace molekulární dynamiky * MeSH
• zvířata MeSH
• Check Tag
• lidé MeSH
• zvířata MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• DNA MeSH
• guanin MeSH
• jednovláknová DNA MeSH

We have carried out an extended set of standard and enhanced-sampling MD simulations (for a cumulative simulation time of 620 μs) with the aim to study folding landscapes of the rGGGUUAGGG and rGGGAGGG parallel G-hairpins (PH) with propeller loop. We identify folding and unfolding pathways of the PH, which is bridged with the unfolded state via an ensemble of cross-like structures (CS) possessing mutually tilted or perpendicular G-strands interacting via guanine-guanine H-bonding. The oligonucleotides reach the PH conformation from the unfolded state via a conformational diffusion through the folding landscape, i.e. as a series of rearrangements of the H-bond interactions starting from compacted anti-parallel hairpin-like structures. Although isolated PHs do not appear to be thermodynamically stable we suggest that CS and PH-types of structures are sufficiently populated during RNA guanine quadruplex (GQ) folding within the context of complete GQ-forming sequences. These structures may participate in compact coil-like ensembles that involve all four G-strands and already some bound ions. Such ensembles can then rearrange into the fully folded parallel GQs via conformational diffusion. We propose that the basic atomistic folding mechanism of propeller loops suggested in this work may be common for their formation in RNA and DNA GQs.

• MeSH
• G-kvadruplexy * MeSH
• guanin chemie   metabolismus   MeSH
• kinetika MeSH
• RNA chemie   metabolismus   MeSH
• sbalování RNA * MeSH
• sekvence nukleotidů MeSH
• simulace molekulární dynamiky MeSH
• termodynamika MeSH
• vodíková vazba MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• guanin MeSH
• RNA MeSH

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 chemie   genetika   metabolismus   MeSH
• G-kvadruplexy * MeSH
• lidé MeSH
• sekvence nukleotidů MeSH
• shluková analýza MeSH
• simulace molekulární dynamiky * MeSH
• telomery genetika   MeSH
• voda chemie   MeSH
• Check Tag
• lidé MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• DNA MeSH
• voda MeSH

Guanine-rich oligonucleotides can form a unique G-quadruplex (GQ) structure with stacking units of four guanine bases organized in a plane through Hoogsteen bonding. GQ structures have been detected in vivo and shown to exert their roles in maintaining genome integrity and regulating gene expression. Understanding GQ conformation is important for understanding its inherent biological role and for devising strategies to control and manipulate functions based on targeting GQ. Although a number of biophysical methods have been used to investigate structure and dynamics of GQs, our understanding is far from complete. As such, this work explores the use of the site-directed spin labeling technique, complemented by molecular dynamics simulations, for investigating GQ conformations. A nucleotide-independent nitroxide label (R5), which has been previously applied for probing conformations of noncoding RNA and DNA duplexes, is attached to multiple sites in a 22-nucleotide DNA strand derived from the human telomeric sequence (hTel-22) that is known to form GQ. The R5 labels are shown to minimally impact GQ folding, and inter-R5 distances measured using double electron-electron resonance spectroscopy are shown to adequately distinguish the different topological conformations of hTel-22 and report variations in their occupancies in response to changes of the environment variables such as salt, crowding agent, and small molecule ligand. The work demonstrates that the R5 label is able to probe GQ conformation and establishes the base for using R5 to study more complex sequences, such as those that may potentially form multimeric GQs in long telomeric repeats.

• MeSH
• G-kvadruplexy * MeSH
• konformace nukleové kyseliny MeSH
• konformace proteinů MeSH
• lidé MeSH
• oligonukleotidy chemie   MeSH
• oxid dusný chemie   MeSH
• Check Tag
• lidé MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Research Support, N.I.H., Extramural MeSH
• Research Support, U.S. Gov't, Non-P.H.S. MeSH
• Názvy látek
• oligonukleotidy MeSH
• oxid dusný 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 chemie   MeSH
• G-kvadruplexy * MeSH
• kationty chemie   MeSH
• lidé MeSH
• Oxytricha genetika   MeSH
• simulace molekulární dynamiky * MeSH
• telomery chemie   MeSH
• Check Tag
• lidé MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• DNA MeSH
• kationty MeSH

The 22-mer c-kit promoter sequence folds into a parallel-stranded quadruplex with a unique structure, which has been elucidated by crystallographic and NMR methods and shows a high degree of structural conservation. We have carried out a series of extended (up to 10 μs long, ∼50 μs in total) molecular dynamics simulations to explore conformational stability and loop dynamics of this quadruplex. Unfolding no-salt simulations are consistent with a multi-pathway model of quadruplex folding and identify the single-nucleotide propeller loops as the most fragile part of the quadruplex. Thus, formation of propeller loops represents a peculiar atomistic aspect of quadruplex folding. Unbiased simulations reveal μs-scale transitions in the loops, which emphasizes the need for extended simulations in studies of quadruplex loops. We identify ion binding in the loops which may contribute to quadruplex stability. The long lateral-propeller loop is internally very stable but extensively fluctuates as a rigid entity. It creates a size-adaptable cleft between the loop and the stem, which can facilitate ligand binding. The stability gain by forming the internal network of GA base pairs and stacks of this loop may be dictating which of the many possible quadruplex topologies is observed in the ground state by this promoter quadruplex.

• MeSH
• denaturace nukleových kyselin MeSH
• draslík chemie   MeSH
• G-kvadruplexy * MeSH
• kationty MeSH
• párování bází MeSH
• promotorové oblasti (genetika) * MeSH
• protoonkogenní proteiny c-kit genetika   MeSH
• simulace molekulární dynamiky MeSH
• sodík chemie   MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• draslík MeSH
• kationty MeSH
• protoonkogenní proteiny c-kit MeSH
• sodík MeSH

We report here that a tetra-substituted naphthalene-diimide derivative (MM41) has significant in vivo anti-tumour activity against the MIA PaCa-2 pancreatic cancer xenograft model. IV administration with a twice-weekly 15 mg/kg dose produces ca 80% tumour growth decrease in a group of tumour-bearing animals. Two animals survived tumour-free after 279 days. High levels of MM41 are rapidly transported into cell nuclei and were found to accumulate in the tumour. MM41 is a quadruplex-interactive compound which binds strongly to the quadruplexes encoded in the promoter sequences of the BCL-2 and k-RAS genes, both of which are dis-regulated in many human pancreatic cancers. Levels of BCL-2 were reduced by ca 40% in tumours from MM41-treated animals relative to controls, consistent with BCL-2 being a target for MM41. Molecular modelling suggests that MM41 binds to a BCL-2 quadruplex in a manner resembling that previously observed in co-crystal structures with human telomeric quadruplexes. This supports the concept that MM41 (and by implication other quadruplex-targeting small molecules) can bind to quadruplex-forming promoter regions in a number of genes and down-regulate their transcription. We suggest that quadruplexes within those master genes that are up-regulated drivers for particular cancers, may be selective targets for compounds such as MM41.

• MeSH
• antitumorózní látky chemie   farmakologie   MeSH
• exprese genu MeSH
• G-kvadruplexy MeSH
• genetická transkripce MeSH
• imidy chemie   farmakologie   MeSH
• injekce intravenózní MeSH
• lidé MeSH
• myši nahé MeSH
• myši MeSH
• nádorové buněčné linie MeSH
• nádory slinivky břišní farmakoterapie   genetika   metabolismus   patologie   MeSH
• naftaleny chemie   farmakologie   MeSH
• promotorové oblasti (genetika) * MeSH
• protoonkogenní proteiny c-bcl-2 antagonisté a inhibitory   chemie   genetika   metabolismus   MeSH
• protoonkogenní proteiny p21(ras) genetika   metabolismus   MeSH
• rozvrh dávkování léků MeSH
• simulace molekulární dynamiky MeSH
• tumor burden účinky léků   MeSH
• xenogenní modely - testy antitumorózní aktivity MeSH
• zvířata MeSH
• Check Tag
• lidé MeSH
• myši MeSH
• ženské pohlaví MeSH
• zvířata MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• antitumorózní látky MeSH
• Bcl2 protein, mouse MeSH Prohlížeč
• Hras protein, mouse MeSH Prohlížeč
• imidy MeSH
• naftaleny MeSH
• naphthalenediimide MeSH Prohlížeč
• protoonkogenní proteiny c-bcl-2 MeSH
• protoonkogenní proteiny p21(ras) MeSH