High Mechanical Stability and Slow Unfolding Rates Are Prevalent in Parallel-Stranded DNA G-Quadruplexes

Unfolding rates of 1:1 and 2:1 complex of CX-5461 and c-MYC promoter G-quadruplexes revealed by single-molecule force spectroscopy

CX-5461, also known as pidnarulex, is a strong G4 stabilizer and has received FDA fast-track designation for BRCA1- and BRCA2- mutated cancers. However, quantitative measurements of the unfolding rates of CX-5461-G4 complexes which are important for the regulation function of G4s, remain lacking. Here, we employ single-molecule magnetic tweezers to measure the unfolding force distributions of c-MYC G4s in the presence of different concentrations of CX-5461. The unfolding force distributions exhibit three discrete levels of unfolding force peaks, corresponding to three binding modes. In combination with a fluorescent quenching assay and molecular docking to previously reported ligand-c-MYC G4 structure, we assigned the ~69 pN peak corresponding to the 1:1 (ligand:G4) complex where CX-5461 binds at the G4's 5'-end. The ~84 pN peak is attributed to the 2:1 complex where CX-5461 occupies both the 5' and 3'. Furthermore, using the Bell-Arrhenius model to fit the unfolding force distributions, we determined the zero-force unfolding rates of 1:1, and 2:1 complexes to be (2.4 ± 0.9) × 10-8 s-1 and (1.4 ± 1.0) × 10-9 s-1 respectively. These findings provide valuable insights for the development of G4-targeted ligands to combat c-MYC-driven cancers.

Mechanical Stability and Unfolding Pathways of Parallel Tetrameric G-Quadruplexes Probed by Pulling Simulations

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.

Distinguishing G-Quadruplexes Stabilizer and Chaperone for c-MYC Promoter G-Quadruplexes through Single-Molecule Manipulation

G-quadruplex (G4) selective stabilizing ligands can regulate c-MYC gene expression, but the kinetic basis remains unclear. Determining the effects of ligands on c-MYC promoter G4s' folding/unfolding kinetics is challenging due to the polymorphic nature of G4s and the high energy barrier to unfold c-MYC promoter G4s. Here, we used single-molecule magnetic tweezers to manipulate a duplex hairpin containing a c-MYC promoter sequence to mimic the transiently denatured duplex during transcription. We measured the effects of six commonly used G4s binding ligands on the competition between quadruplex and duplex structures, as well as the folding/unfolding kinetics of G4s. Our results revealed two distinct roles for G4s selective stabilization: CX-5461 is mainly acting as c-MYC G4s stabilizer, reducing the unfolding rate (ku) of c-MYC G4s, whereas PDS and 360A also act as G4s chaperone, accelerating the folding rates (kf) of c-MYC G4s. qRT-PCR results obtained from CA46 and Raji cell lines demonstrated that G4s stabilizing ligands can downregulate c-MYC expression, while G4s stabilizer CX-5461 exhibited the strongest c-MYC gene suppression. These results shed light on the potential of manipulating G4s' folding/unfolding kinetics by ligands for precise regulation of promoter G4-associated biological activities.

A TeZla micromixer for interrogating the early and broad folding landscape of G-quadruplex via multistage velocity descending

Biomacromolecular folding kinetics involves fast folding events and broad timescales. Current techniques face limitations in either the required time resolution or the observation window. In this study, we developed the TeZla micromixer, integrating Tesla and Zigzag microstructures with a multistage velocity descending strategy. TeZla achieves a significant short mixing dead time (40 µs) and a wide time window covering four orders of magnitude (up to 300 ms). Using this unique micromixer, we explored the folding landscape of c-Myc G4 and its noncanonical-G4 derivatives with different loop lengths or G-vacancy sites. Our findings revealed that c-Myc can bypass folding intermediates and directly adopt a G4 structure in the cation-deficient buffer. Moreover, we found that the loop length and specific G-vacancy site could affect the folding pathway and significantly slow down the folding rates. These results were also cross-validated with real-time NMR and circular dichroism. In conclusion, TeZla represents a versatile tool for studying biomolecular folding kinetics, and our findings may ultimately contribute to the design of drugs targeting G4 structures.

Complexity of Guanine Quadruplex Unfolding Pathways Revealed by Atomistic Pulling Simulations

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.

An overview of single-molecule techniques and applications in the study of nucleic acid structure and function

Nucleic acids are an indispensable component in all known life forms. The biological processes are regulated by Nucleic acids, which associate to form special high-order structures. since the high-level structures of nucleic acids are related to gene expression in cancer cells or viruses, it is very likely to become a potential drug target. Traditional biochemical methods are limited to distinguish the conformational distribution and dynamic transition process of single nucleic acid structure. The ligands based on the intermediate and transition states between different conformations are not designed by traditional biochemical methods. The single-molecule techniques enable real-time observation of the individual nucleic acid behavior due to its high resolution. Here, we introduce the application of single-molecule techniques in the study of small molecules to recognize nucleic acid structures, such as single-molecule FRET, magnetic tweezers, optical tweezers and atomic force microscopy. At the same time, we also introduce the specific advantages of single-molecule technology compared with traditional biochemical methods and some problems arisen in current research.

Characterization of G-Quadruplexes Folding/Unfolding Dynamics and Interactions with Proteins from Single-Molecule Force Spectroscopy

G-quadruplexes (G4s) are stable secondary nucleic acid structures that play crucial roles in many fundamental biological processes. The folding/unfolding dynamics of G4 structures are associated with the replication and transcription regulation functions of G4s. However, many DNA G4 sequences can adopt a variety of topologies and have complex folding/unfolding dynamics. Determining the dynamics of G4s and their regulation by proteins remains challenging due to the coexistence of multiple structures in a heterogeneous sample. Here, in this mini-review, we introduce the application of single-molecule force-spectroscopy methods, such as magnetic tweezers, optical tweezers, and atomic force microscopy, to characterize the polymorphism and folding/unfolding dynamics of G4s. We also briefly introduce recent studies using single-molecule force spectroscopy to study the molecular mechanisms of G4-interacting proteins.

Mechanical diversity and folding intermediates of parallel-stranded G-quadruplexes with a bulge

A significant number of sequences in the human genome form noncanonical G-quadruplexes (G4s) with bulges or a guanine vacancy. Here, we systematically characterized the mechanical stability of parallel-stranded G4s with a one to seven nucleotides bulge at various positions. Our results show that G4-forming sequences with a bulge form multiple conformations, including fully-folded G4 with high mechanical stability (unfolding forces > 40 pN), partially-folded intermediates (unfolding forces < 40 pN). The folding probability and folded populations strongly depend on the positions and lengths of the bulge. By combining a single-molecule unfolding assay, dimethyl sulfate (DMS) footprinting, and a guanine-peptide conjugate that selectively stabilizes guanine-vacancy-bearing G-quadruplexes (GVBQs), we identified that GVBQs are the major intermediates of G4s with a bulge near the 5′ or 3′ ends. The existence of multiple structures may induce different regulatory functions in many biological processes. This study also demonstrates a new strategy for selectively stabilizing the intermediates of bulged G4s to modulate their functions.

Folding Dynamics of DNA G-Quadruplexes Probed by Millisecond Temperature Jump Circular Dichroism

G-quadruplexes play important roles in cellular regulatory functions, but despite significant experimental and theoretical efforts, their folding mechanisms remain poorly understood. In this context, we developed a T-jump experiment to access the thermal denaturation and renaturation dynamics of short intramolecular G-quadruplexes in vitro, on the time scale of a few hundred milliseconds. With this new setup, we compared the thermal denaturation and renaturation kinetics of three antiparallel topologies made of the human telomeric sequences d[(5'-GGG(TTAGGG)₃-3']/Na⁺ and d[5'-AGGG(TTAGGG)₃-3']/Na⁺ and the thrombin-binding aptamer sequence d[5'-GGTTGGTGTGGTTGG-3']/K⁺, with those of the parallel topology made of the human CEB25 minisatellite d[5'-AAGGGTGGGTGTAAGTGTGGGTGGGT-3']/Na⁺. In all cases, exponential kinetics of the order of several hundred milliseconds were observed. Measurements performed for different initial temperatures revealed distinct denaturation and renaturation dynamics, ruling out a simple two-state mechanism. The parallel topology, in which all guanines adopt an anti conformation, displays much slower dynamics than antiparallel topologies associated with very low activation barriers. This behavior can be explained by the constrained conformational space due to the presence of the single-base propeller loops that likely hinders the movement of the coiled DNA strand and reduces the contribution of the entropy during the renaturation process at high temperatures.