Guiding the folding of G-quadruplexes through loop residue interactions

Structural elucidation of HIV-1 G-quadruplexes in a cellular environment and their ligand binding using responsive 19F-labeled nucleoside probes

Understanding the structure and recognition of highly conserved regulatory segments of the integrated viral DNA genome that forms unique topologies can greatly aid in devising novel therapeutic strategies to counter chronic infections. In this study, we configured a probe system using highly environment-sensitive nucleoside analogs, 5-fluoro-2'-deoxyuridine (FdU) and 5-fluorobenzofuran-2'-deoxyuridine (FBFdU), to investigate the structural polymorphism of HIV-1 long terminal repeat (LTR) G-quadruplexes (GQs) by fluorescence and 19F NMR. FdU and FBFdU, serving as hairpin and GQ sensors, produced distinct spectral signatures for different GQ topologies adopted by LTR G-rich oligonucleotides. Importantly, systematic 19F NMR analysis in Xenopus laevis oocytes gave unprecedented information on the structure adopted by the LTR G-rich region in the cellular environment. The results indicate that it forms a unique GQ-hairpin hybrid architecture, a potent hotspot for selective targeting. Furthermore, structural models generated using MD simulations provided insights on how the probe system senses different GQs. Using the responsiveness of the probes and Taq DNA polymerase stop assay, we monitored GQ- and hairpin-specific ligand interactions and their synergistic inhibitory effect on the replication process. Our findings suggest that targeting GQ and hairpin motifs simultaneously using bimodal ligands could be a new strategy to selectively block the viral replication.

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.

Solution-State Structure of a Long-Loop G-Quadruplex Formed Within Promoters of Plasmodium falciparum B var Genes

We report the high-resolution NMR solution-state structure of an intramolecular G-quadruplex with a diagonal loop of ten nucleotides. The G-quadruplex is formed by a 34-nt DNA sequence, d[CAG3T2A2G3TATA2CT3AG4T2AG3T2], named UpsB-Q-1. This sequence is found within promoters of the var genes of Plasmodium falciparum, which play a key role in malaria pathogenesis and evasion of the immune system. The [3+1]-hybrid G-quadruplex formed under physiologically relevant conditions exhibits a unique equilibrium between two structures, both stabilized by base stacking and non-canonical hydrogen bonding. Unique equilibrium of the two closely related 3D structures originates from a North-South repuckering of deoxyribose moiety of residue T27 in the lateral loop. Besides the 12 guanines involved in three G-quartets, most residues in loop regions are involved in interactions at both G-quartet-loop interfaces.

Loop nucleobases-dependent folding of G-quadruplex in normal and cancer cell-mimicking KCl microenvironments

In this study, the folding of G-quadruplex (G4) from the telomeric DNA sequences having loop nucleobases of different chemical natures, numbers, and arrangements in 10 mM and 100 mM KCl salt conditions mimicking the cancerous and normal KCl salt microenvironments have been investigated. The data suggest that the structure and stability of the G4 are highly dependent on the KCl salt concentration. In general, the conformational flexibility of the folded G4 is higher in KCl salt relevant to cancer than in the normal case for any loop arrangements with the same number of nucleobases. The stability of the G4 decreases with the increase in the number of loop nucleobases for both salt conditions. However, the decrease in the stability of G4 having adenine in the loop region is significantly higher than the case of thymine, particularly more prominent in the KCl salt relevant to the cancer. The topology of the folded G4 and its stability also depend delicately on the permutation of the nucleobases in the loop and the salt concentrations for a particular sequence. The findings indicate that the structure and stability of G4 are noticeably different in KCl salt relevant to physiological and cancer conditions.

Structural Differences at Quadruplex-Duplex Interfaces Enable Ligand-Induced Topological Transitions

Quadruplex-duplex (QD) junctions, which represent unique structural motifs of both biological and technological significance, have been shown to constitute high-affinity binding sites for various ligands. A QD hybrid construct based on a human telomeric sequence, which harbors a duplex stem-loop in place of a short lateral loop, is structurally characterized by NMR. It folds into two major species with a (3+1) hybrid and a chair-type (2+2) antiparallel quadruplex domain coexisting in a K+ buffer solution. The antiparallel species is stabilized by an unusual capping structure involving a thymine and protonated adenine base AH+ of the lateral loop facing the hairpin duplex to form a T·AH+ ·G·C quartet with the interfacial G·C base pair at neutral pH. Addition and binding of Phen-DC3 to the QD hybrid mixture by its partial intercalation at corresponding QD junctions leads to a topological transition with exclusive formation of the (3+1) hybrid fold. In agreement with the available experimental data, such an unprecedented discrimination of QD junctions by a ligand can be rationalized following an induced fit mechanism.

Stabilization of G-Quadruplex Structures of the SARS-CoV-2 Genome by TMPyP4, BRACO19, and PhenDC3

The G-quadruplex is one of the non-canonical structures formed by nucleic acids, which can be formed by guanine-rich sequences. They became the focus of much research when they were found in several oncogene promoter regions and also in the telomeres. Later on, they were discovered in viruses as well. Various ligands have been developed in order to stabilize DNA G-quadruplexes, which were believed to have an anti-cancer or antiviral effect. We investigated three of these ligands, and whether they can also affect the stability of the G-quadruplex-forming sequences of the RNA genome of SARS-CoV-2. All three investigated oligonucleotides showed the G-quadruplex form. We characterized their stability and measured their thermodynamic parameters using the Förster resonance energy transfer method. The addition of the ligands caused an increase in the unfolding temperature, but this effect was smaller compared to that found earlier in the case of G-quadruplexes of the hepatitis B virus, which has a DNA genome.

DNA Quadruplex Structure with a Unique Cation Dependency

DNA quadruplex structures provide an additional layer of regulatory control in genome maintenance and gene expression and are widely used in nanotechnology. We report the discovery of an unprecedented tetrastranded structure formed from a native G-rich DNA sequence originating from the telomeric region of Caenorhabditis elegans. The structure is defined by multiple properties that distinguish it from all other known DNA quadruplexes. Most notably, the formation of a stable so-called KNa-quadruplex (KNaQ) requires concurrent coordination of K+ and Na+ ions at two distinct binding sites. This structure provides novel insight into G-rich DNA folding under ionic conditions relevant to eukaryotic cell physiology and the structural evolution of telomeric DNA. It highlights the differences between the structural organization of human and nematode telomeric DNA, which should be considered when using C. elegans as a model in telomere biology, particularly in drug screening applications. Additionally, the absence/presence of KNaQ motifs in the host/parasite introduces an intriguing possibility of exploiting the KNaQ fold as a plausible antiparasitic drug target. The structure's unique shape and ion dependency and the possibility of controlling its folding by using low-molecular-weight ligands can be used for the design or discovery of novel recognition DNA elements and sensors.

Impact of loop length and duplex extensions on the design of hybrid-type G-quadruplexes

A G-rich core sequence G3-TCA-G3-T1,2-G3-T1,2-G3 can be designed to fold into a parallel or into two different (3+1) hybrid-type G-quadruplexes, among them an elusive topology with one lateral followed by two propeller loops. Favored folds can be rationalized based on the number of intervening thymidines and on additional complementary flanking sequences.

Showcasing Different G-Quadruplex Folds of a G-Rich Sequence: Between Rule-Based Prediction and Butterfly Effect

In better understanding the interactions of G-quadruplexes in a cellular or noncellular environment, a reliable sequence-based prediction of their three-dimensional fold would be extremely useful, yet is often limited by their remarkable structural diversity. A G-rich sequence related to a promoter sequence of the PDGFR-β nuclease hypersensitivity element (NHE) comprises a G3-G3-G2-G4-G3 pattern of five G-runs with two to four G residues. Although the predominant formation of three-layered canonical G-quadruplexes with uninterrupted G-columns can be expected, minimal base substitutions in a non-G-tract domain were shown to guide folding into either a basket-type antiparallel quadruplex, a parallel-stranded quadruplex with an interrupted G-column, a quadruplex with a V-shaped loop, or a (3+1) hybrid quadruplex. A 3D NMR structure for each of the different folds was determined. Supported by thermodynamic profiling on additional sequence variants, formed topologies were rationalized by the identification and assessment of specific critical interactions of loop and overhang residues, giving valuable insights into their contribution to favor a particular conformer. The variability of such tertiary interactions, together with only small differences in quadruplex free energies, emphasizes current limits for a reliable sequence-dependent prediction of favored topologies from sequences with multiple irregularly positioned G-tracts.

Terminal Alkyne-Modified DNA Aptamers with Enhanced Protein Binding Affinities

Nucleic acid-based receptors, known as aptamers, are relatively fast to discover and manufacture but lack the diverse functional groups of protein receptors (e.g., antibodies). The binding properties of DNA aptamers can be enhanced by attaching abiotic functional groups; for example, aromatic groups such as naphthalene slow dissociation from proteins. Although the terminal alkyne is a π-electron-rich functional group that has been used in small molecule drugs to enhance binding to proteins through noncovalent interactions, it remains unexplored for enhancing DNA aptamer binding affinity. Here, we demonstrate the utility of the terminal alkyne for improving the binding of DNA to proteins. We prepared a library of 256 terminal-alkyne-bearing variants of HD22, a DNA aptamer that binds the protein thrombin with nanomolar affinity. After a one-step thrombin-binding selection, a high-affinity aptamer containing two alkynes was discovered, exhibiting 3.2-fold tighter thrombin binding than the corresponding unmodified sequence. The tighter binding was attributable to a slower rate of dissociation from thrombin (5.2-fold slower than HD22). Molecular dynamics simulations with enhanced sampling by Replica Exchange with Solute Tempering (REST2) suggest that the π-electron-rich alkyne interacts with an asparagine side chain N-H group on thrombin, forming a noncovalent interaction that stabilizes the aptamer-protein interface. Overall, this work represents the first case of terminal alkynes enhancing the binding properties of an aptamer and underscores the utility of the terminal alkyne as an atom economical π-electron-rich functional group to enhance binding affinity with minimal steric perturbation.

Understanding the Interactions of Guanine Quadruplexes with Peptides as Novel Strategies for Diagnosis or Tuning Biological Functions

Guanine quadruplexes (G4s) are nucleic acid structures exhibiting a complex structural behavior and exerting crucial biological functions in both cells and viruses. The specific interactions of peptides with G4s, as well as an understanding of the factors driving the specific recognition are important for the rational design of both therapeutic and diagnostic agents. In this review, we examine the most important studies dealing with the interactions between G4s and peptides, highlighting the strengths and limitations of current analytic approaches. We also show how the combined use of high-level molecular simulation techniques and experimental spectroscopy is the best avenue to design specifically tuned and selective peptides, thus leading to the control of important biological functions.

Structural and Functional Classification of G-Quadruplex Families within the Human Genome

G-quadruplexes (G4s) are short secondary DNA structures located throughout genomic DNA and transcribed RNA. Although G4 structures have been shown to form in vivo, no current search tools that examine these structures based on previously identified G-quadruplexes and filter them based on similar sequence, structure, and thermodynamic properties are known to exist. We present a framework for clustering G-quadruplex sequences into families using the CD-HIT, MeShClust, and DNACLUST methods along with a combination of Starcode and BLAST. Utilizing this framework to filter and annotate clusters, 95 families of G-quadruplex sequences were identified within the human genome. Profiles for each family were created using hidden Markov models to allow for the identification of additional family members and generate homology probability scores. The thermodynamic folding energy properties, functional annotation of genes associated with the sequences, scores from different prediction algorithms, and transcription factor binding motifs within a family were used to annotate and compare the diversity within and across clusters. The resulting set of G-quadruplex families can be used to further understand how different regions of the genome are regulated by factors targeting specific structures common to members of a specific cluster.