Streamlining effects of extra telomeric repeat on telomeric DNA folding revealed by fluorescence-force spectroscopy

Single-Macromolecule Studies of Eukaryotic Genomic Maintenance

Genomes are self-organized and self-maintained as long, complex macromolecules of chromatin. The inherent heterogeneity, stochasticity, phase separation, and chromatin dynamics of genome operation make it challenging to study genomes using ensemble methods. Various single-molecule force-, fluorescent-, and sequencing-based techniques rooted in different disciplines have been developed to fill critical gaps in the capabilities of bulk measurements, each providing unique, otherwise inaccessible, insights into the structure and maintenance of the genome. Capable of capturing molecular-level details about the organization, conformational changes, and packaging of genetic material, as well as processive and stochastic movements of maintenance factors, a single-molecule toolbox provides an excellent opportunity for collaborative research to understand how genetic material functions in health and malfunctions in disease. In this review, we discuss novel insights brought to genomic sciences by single-molecule techniques and their potential to continue to revolutionize the field-one molecule at a time.

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.

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.

DNA-RNA hybrid G-quadruplex tends to form near the 3' end of telomere overhang

Telomeric repeat-containing RNA (TERRA) has been suggested to participate in telomere maintenance. TERRA consisting of UUAGGG repeats is capable of forming an intermolecular G-quadruplex (GQ) with single-stranded TTAGGG-repeat DNA in the telomere 3' overhang. To explore the structural features and potential functions of this DNA-RNA hybrid GQ (HGQ), we used single-molecule FRET to study the folding patterns of DNA with four to seven telomeric tandem repeats annealed with a short RNA consisting of two or five telomeric repeats. Our data highlight that RNA prefers to form DNA-RNA HGQ near the 3' end of telomeric DNA. Furthermore, the unfolding of secondary structures by a complementary C-rich sequence was observed for DNA GQ but not for DNA-RNA HGQ, which demonstrated the enhanced stability of the telomere 3' end via hybridization with RNA. These conformational and physical properties of telomeric DNA-RNA HGQ suggest that TERRA might limit access to the 3' end of the telomeric DNA overhang, which is known to be critical for the interaction with telomerase and other telomere-associated proteins.

Vectorial folding of telomere overhang promotes higher accessibility

Human telomere overhang composed of tandem repeats of TTAGGG folds into G-quadruplex (G4). Unlike in an experimental setting in the test tube in which the entire length is allowed to fold at once, inside the cell, the overhang is expected to fold as it is synthesized directionally (5' to 3') and released segmentally by a specialized enzyme, the telomerase. To mimic such vectorial G4 folding process, we employed a superhelicase, Rep-X which can unwind DNA to release the TTAGGG repeats in 5' to 3' direction. We demonstrate that the folded conformation achieved by the refolding of full sequence is significantly different from that of the vectorial folding for two to eight TTAGGG repeats. Strikingly, the vectorially folded state leads to a remarkably higher accessibility to complementary C-rich strand and the telomere binding protein POT1, reflecting a less stably folded state resulting from the vectorial folding. Importantly, our study points to an inherent difference between the co-polymerizing and post-polymerized folding of telomere overhang that can impact telomere architecture and downstream processes.

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.

Multi-start Evolutionary Nonlinear OpTimizeR (MENOTR): A hybrid parameter optimization toolbox

Parameter optimization or "data fitting" is a computational process that identifies a set of parameter values that best describe an experimental data set. Parameter optimization is commonly carried out using a computer program utilizing a non-linear least squares (NLLS) algorithm. These algorithms work by continuously refining a user supplied initial guess resulting in a systematic increase in the goodness of fit. A well-understood problem with this class of algorithms is that in the case of models with correlated parameters the optimized output parameters are initial guess dependent. This dependency can potentially introduce user bias into the resultant analysis. While many optimization programs exist, few address this dilemma. Here we present a data analysis tool, MENOTR, that is capable of overcoming the initial guess dependence in parameter optimization. Several case studies with published experimental data are presented to demonstrate the capabilities of this tool. The results presented here demonstrate how to effectively overcome the initial guess dependence of NLLS leading to greater confidence that the resultant optimized parameters are the best possible set of parameters to describe an experimental data set. While the optimization strategies implemented within MENOTR are not entirely novel, the application of these strategies to optimize parameters in kinetic and thermodynamic biochemical models is uncommon. MENOTR was designed to require minimal modification to accommodate a new model making it immediately accessible to researchers with a limited programming background. We anticipate that this toolbox can be used in a wide variety of data analysis applications. Prototype versions of this toolbox have been used in a number of published investigations already, as well as ongoing work with chemical-quenched flow, stopped-flow, and molecular tweezers data sets. STATEMENT OF SIGNIFICANCE: Non-linear least squares (NLLS) is a common form of parameter optimization in biochemistry kinetic and thermodynamic investigations These algorithms are used to fit experimental data sets and report corresponding parameter values. The algorithms are fast and able to provide good quality solutions for models involving few parameters. However, initial guess dependence is a well-known drawback of this optimization strategy that can introduce user bias. An alternative method of parameter optimization are genetic algorithms (GA). Genetic algorithms do not have an initial guess dependence but are slow at arriving at the best set of fit parameters. Here, we present MENOTR, a parameter optimization toolbox utilizing a hybrid GA/NLLS algorithm. The toolbox maximizes the strength of each strategy while minimizing the inherent drawbacks.

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.

Nanomechanics of G-quadruplexes within the promoter of the KIT oncogene

G-quadruplexes (G4s) are tetrahelical DNA structures stabilized by four guanines paired via Hoogsteen hydrogen bonds into quartets. While their presence within eukaryotic DNA is known to play a key role in regulatory processes, their functional mechanisms are still under investigation. In the present work, we analysed the nanomechanical properties of three G4s present within the promoter of the KIT proto-oncogene from a single-molecule point of view through the use of magnetic tweezers (MTs). The study of DNA extension fluctuations under negative supercoiling allowed us to identify a characteristic fingerprint of G4 folding. We further analysed the energetic contribution of G4 to the double-strand denaturation process in the presence of negative supercoiling, and we observed a reduction in the energy required for strands separation.

Transcriptional regulation of telomeric repeat-containing RNA by acridine derivatives

Telomere is a specialized DNA-protein complex that plays an important role in maintaining chromosomal integrity. Shelterin is a protein complex formed by six different proteins, with telomeric repeat factors 1 (TRF1) and 2 (TRF2) binding to double-strand telomeric DNA. Telomeric DNA consists of complementary G-rich and C-rich repeats, which could form G-quadruplex and intercalated motif (i-motif), respectively, during cell cycle. Its G-rich transcription product, telomeric repeat-containing RNA (TERRA), is essential for telomere stability and heterochromatin formation. After extensive screening, we found that acridine derivative 2c and acridine dimer DI26 could selectively interact with TRF1 and telomeric i-motif, respectively. Compound 2c blocked the binding of TRF1 with telomeric duplex DNA, resulting in up-regulation of TERRA. Accumulated TERRA could bind with TRF1 at its allosteric site and further destabilize its binding with telomeric DNA. In contrast, DI26 could destabilize telomeric i-motif, resulting in down-regulation of TERRA. Both compounds exhibited anti-tumour activity for A549 cells, but induced different DNA damage pathways. Compound 2c significantly suppressed tumour growth in A549 xenograft mouse model. The function of telomeric i-motif structure was first studied with a selective binding ligand, which could play an important role in regulating TERRA transcription. Our results showed that appropriate level of TERRA transcript could be important for stability of telomere, and acridine derivatives could be further developed as anti-cancer agents targeting telomere. This research increased understanding for biological roles of telomeric i-motif, TRF1 and TERRA, as potential anti-cancer drug targets.

Unraveling protein's structural dynamics: from configurational dynamics to ensemble switching guides functional mesoscale assemblies

Evidence regarding protein structure and function manifest the imperative role that dynamics play in proteins, underlining reconsideration of the unanimated sequence-to-structure-to-function paradigm. Structural dynamics portray a heterogeneous energy landscape described by conformational ensembles where each structural representation can be responsible for unique functions or enable macromolecular assemblies. Using the human p27/Cdk2/Cyclin A ternary complex as an example, we highlight the vital role of intramolecular and intermolecular dynamics for target recognition, binding, and inhibition as a critical modulator of cell division. Rapidly sampling configurations is critical for the population of different conformational ensembles encoding functional roles. To garner this knowledge, we present how the integration of (sub)ensemble and single-molecule fluorescence spectroscopy with molecular dynamic simulations can characterize structural dynamics linking the heterogeneous ensembles to function. The incorporation of dynamics into the sequence-to-structure-to-function paradigm promises to assist in tackling various challenges, including understanding the formation and regulation of mesoscale assemblies inside cells.

Antiviral Activity of the G-Quadruplex Ligand TMPyP4 against Herpes Simplex Virus-1

The herpes simplex virus 1 (HSV-1) genome is extremely rich in guanine tracts that fold into G-quadruplexes (G4s), nucleic acid secondary structures implicated in key biological functions. Viral G4s were visualized in HSV-1 infected cells, with massive virus cycle-dependent G4-formation peaking during viral DNA replication. Small molecules that specifically interact with G4s have been shown to inhibit HSV-1 DNA replication. We here investigated the antiviral activity of TMPyP4, a porphyrin known to interact with G4s. The analogue TMPyP2, with lower G4 affinity, was used as control. We showed by biophysical analysis that TMPyP4 interacts with HSV-1 G4s, and inhibits polymerase progression in vitro; in infected cells, it displayed good antiviral activity which, however, was independent of inhibition of virus DNA replication or entry. At low TMPyP4 concentration, the virus released by the cells was almost null, while inside the cell virus amounts were at control levels. TEM analysis showed that virus particles were trapped inside cytoplasmatic vesicles, which could not be ascribed to autophagy, as proven by RT-qPCR, western blot, and immunofluorescence analysis. Our data indicate a unique mechanism of action of TMPyP4 against HSV-1, and suggest the unprecedented involvement of currently unknown G4s in viral or antiviral cellular defense pathways.

Using Magnetic Tweezers to Unravel the Mechanism of the G-quadruplex Binding and Unwinding Activities of DHX36 Helicase

Single-molecule manipulation methods are useful techniques to probe the interactions of proteins and nucleic acid structures. Here, we describe the magnetic tweezers-based single-molecule investigation of the binding of helicases to G-quadruplex structures and their ATP-dependent unwinding activity, using DHX36 (also known as RHAU and G4R1) helicase and a DNA G-quadruplex structure for an example. We specifically emphasize on the principle and method to probe the interactions between DHX36 and the DNA G-quadruplex in different intermediate states during an ATPase cycle of DHX36, based on detecting the DHX36-induced changes in the lifetime of the DNA G-quadruplex under tension. The principle of the measurement can be broadly extended to the studies of other DNA or RNA G-quadruplex helicases.

Real-time monitoring of single ZTP riboswitches reveals a complex and kinetically controlled decision landscape

RNAs begin to fold and function during transcription. Riboswitches undergo cotranscriptional switching in the context of transcription elongation, RNA folding, and ligand binding. To investigate how these processes jointly modulate the function of the folate stress-sensing Fusobacterium ulcerans ZTP riboswitch, we apply a single-molecule vectorial folding (VF) assay in which an engineered superhelicase Rep-X sequentially releases fluorescently labeled riboswitch RNA from a heteroduplex in a 5′-to-3′ direction, at ~60 nt s⁻¹ [comparable to the speed of bacterial RNA polymerase (RNAP)]. We demonstrate that the ZTP riboswitch is kinetically controlled and that its activation is favored by slower unwinding, strategic pausing between but not before key folding elements, or a weakened transcription terminator. Real-time single-molecule monitoring captures folding riboswitches in multiple states, including an intermediate responsible for delayed terminator formation. These results show how individual nascent RNAs occupy distinct channels within the folding landscape that controls the fate of the riboswitch.

Many RNAs become functional before their synthesis completes. Here the authors employ a single-molecule vectorial folding assay mimicking RNA transcription and show that the ZTP riboswitch is kinetically controlled and activated by slower unwinding and strategic pausing.