Conformational Plasticity of Parallel G-Quadruplex─Implications on Duplex-Quadruplex Motifs.

DNA G-quadruplexes are essential motifs in molecular biology performing a wide range of functions enabled by their unique and diverse structures. In this study, we focus on the conformational plasticity of the most abundant and biologically relevant parallel G-quadruplex topology. A multipronged approach of structure survey, solution-state NMR spectroscopy, and molecular dynamics simulations unravels subtle yet essential features of the parallel G-quadruplex topology. Stark differences in flexibility are observed for the nucleotides depending upon their positioning in the tetrad planes that are intricately correlated with the conformational sampling of the propeller loop. Importantly, the terminal nucleotides in the 5'-end versus the 3'-end of the parallel quadruplex display differential dynamics that manifests their ability to accommodate a duplex on either end of the G-quadruplex. The conformational plasticity characterized in this study provides essential cues toward biomolecular processes such as small molecular binding, intermolecular quadruplex stacking, and implications on how a duplex influences the structure of a neighboring quadruplex.

Why are G-quadruplexes good at preventing protein aggregation?

Maintaining a healthy protein folding environment is essential for cellular function. Recently, we found that nucleic acids, G-quadruplexes in particular, are potent chaperones for preventing protein aggregation. With the aid of structure-function and NMR analyses of two G-quadruplex forming sequences, PARP-I and LTR-III, we uncovered several contributing factors that affect G-quadruplexes in preventing protein aggregation. Notably, three factors emerged as vital in determining holdase activity of G-quadruplexes: their structural topology, G-quadruplex accessibility and dynamics, and oligomerization state. These factors together appear to largely dictate whether a G-quadruplex is able to prevent partially misfolded proteins from aggregating. Understanding the physical traits that govern the ability of G-quadruplexes to modulate protein aggregation will help elucidate their possible roles in neurodegenerative disease.

5-Formylcytosine weakens the G-C pair and imparts local conformational fluctuations to DNA duplexes.

DNA epigenetic modifications such as 5-methyl (5mC), 5-hydroxymethyl (5hmC), 5-formyl (5fC) and 5-carboxyl (5caC) cytosine have unique and specific biological roles. Crystallographic studies of 5mC containing duplexes were conducted in the A-, B- or the intermediate E-DNA polymorphic forms. 5fC-modified duplexes initially observed in the disputed F-DNA architecture were subsequently crystallized in the A-form, suggesting that epigenetic modifications enable DNA sequences to adopt diverse conformational states that plausibly contribute to their function. Solution-state studies of these modifications were found in the B-DNA form, with marked differences in the conformational flexibility of 5fC containing duplexes in comparison to C/5mC containing duplexes, compromising the DNA duplex's stability. Herein, we systematically evaluate sensitive and commonly inaccessible NMR parameters to map the subtle differences between C, 5mC, and their oxidized (5hmC/5fC) counterparts. We observe that 15N/1H chemical shifts effectively report on the weakening of 5fC-G Watson-Crick base-pair H-bonding, extending the instability beyond any achievable within the sequence-specific changes in DNA. Triple 5fC containing sequences propagate the destabilization farther from the site of modifications, explaining reduced duplex stability upon multiple modifications. Additionally, scalar and residual dipolar coupling measurements unravel local sugar pucker fluctuations. One-bond 13C-1H scalar coupling measurements point towards a significant deviation away from the anticipated C2'-endo pucker for the 5fC modified nucleotide. Structural models obtained employing 13C-1H residual dipolar couplings and inter-proton distances corroborate the sugar pucker's deviation for 5fC modified DNA duplexes. The changes in the sugar pucker equilibria remain local to the 5fC modified nucleotide sans additive/long-range effects arising from multiple contiguous modifications. These observations highlight the impact of a major groove modification that alters the physical properties of DNA duplex without disturbing the Watson-Crick face. The changes observed in our studies for the 5fC containing DNA contrast with the perturbations induced by damage/lesion highlight the varied conformational preferences that modified nucleobases impart to the DNA duplex. As sequence-specific DNA transactions are rooted in the base-pair stability and pucker deviations, the observed structural perturbations for 5fC-modified DNA potentially play critical functional roles, such as protein-DNA recognition and interactions.

Bisindolylmaleimide Ligands Stabilize c-MYC G-Quadruplex DNA Structure and Downregulate Gene Expression.

G-Quadruplex (G4) structures play a pivotal role in diverse biological functions, including essential processes, such as telomere maintenance and gene regulation. G4 structures formed in functional regions of genomes are actively pursued toward therapeutics and are targeted by small-molecule ligands that alter their structure and/or stability. Herein, we report the synthesis of bisindolylmaleimide-based (BIM) ligands, which preferentially stabilize parallel G4 structures of c-MYC and c-KIT oncogenes over the telomeric h-RAS1 G4 and duplex DNAs. The preferential stabilization of parallel G4s with BIM ligands is further validated by the DNA polymerase stop assay, where stop products were only observed for templates containing the c-MYC G4 sequence. Nuclear magnetic resonance (NMR) titration studies indicate that the lead ligand BIM-Pr1 forms a 2:1 complex with c-MYC G4 DNA with a KD of 38 ± 5 µM. The BIM ligand stacks at the 5' and 3' quartets, with molecular modeling and dynamics studies supporting the proposed binding mode. The ligand is cytotoxic to HeLa cells and downregulates c-MYC gene expression. Collectively, the results present bisindolylmaleimide scaffolds as novel and powerful G4 targeting agents.

NIX-mediated mitophagy regulate metabolic reprogramming in phagocytic cells during mycobacterial infection.

RNASeq analysis of PBMCs from treatment naïve TB patients and healthy controls revealed that M. tuberculosis (Mtb) infection dysregulates several metabolic pathways and upregulates BNIP3L/NIX receptor mediated mitophagy. Analysis of publicly available transcriptomic data from the NCBI-GEO database indicated that M. bovis (BCG) infection also induces similar rewiring of metabolic and mitophagy pathways. Mtb chronic infection and BCG in-vitro infection both downregulated oxidative phosphorylation and upregulated glycolysis and mitophagy; therefore, we used non-pathogenic mycobacterial species BCG as a model for Mtb infection to gain molecular insights and outcomes of this phenomenon. BCG infection in PBMCs and THP-1 macrophages induce mitophagy and glycolysis, leading to differentiation of naïve macrophage to M1 phenotype. Glucose consumption and lactate production were quantified by NMR, while the mitochondrial mass assessment was performed by mitotracker red uptake assay. Infected macrophages predominantly exhibit M1-phenotype, which is indicated by an increase in M1 specific cytokines (IL-6, TNF-α, and IL-1ß) and increased NOS2/ARG1, CD86/CD206 ratio. NIX knockdown abrogates this upregulation of glycolysis, mitophagy, and secretion of pro-inflammatory cytokines in BCG infected cells, indicating that mycobacterial infection-induced immunometabolic changes are executed via NIX mediated mitophagy and are essential for macrophage differentiation and resolution of infection.

Characterization of DNA G-quadruplex Topologies with NMR Chemical Shifts.

G-quadruplexes are nucleic acid motifs formed by stacking of guanosine-tetrad pseudoplanes. They perform varied biological roles, and their distinctive structural features enable diverse applications. High-resolution structural characterization of G-quadruplexes is often time-consuming and expensive, calling for effective methods. Herein, we develop NMR chemical shifts and machine learning-based methodology that allows direct, rapid, and reliable analysis of canonical three-plane DNA G-quadruplexes sans isotopic enrichment. We show, for the first time, that each unique topology enforces a specific distribution of glycosidic torsion angles. Newly acquired carbon chemical shifts are exquisite probes for the dihedral angle distribution and provide immediate and unambiguous backbone topology assignment. The support vector machine learning methodology aids resonance assignment by providing plane indices for tetrad-forming guanosines. We further demonstrate the robustness by successful application of the methodology to a sequence that folds in two dissimilar topologies under different ionic conditions, providing its first atomic-level characterization.

Publisher Correction: High-performance virtual screening by targeting a high-resolution RNA dynamic ensemble.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

A metabolomics footprint approach to understanding the benefits of synbiotics in functional foods and dietary therapeutics for health, communicable and non-communicable diseases.

Gut microbiota have been shown to affect various cellular and host response elements such as immunological, neurological, energy, storage, etc. In recent years, this has led to rapid expansion in dietary products containing probiotics, prebiotics and combination thereof in synbiotics. While benefits of consuming functional foods derived from probiotics strains have been demonstrated for various metabolites, a detailed analysis of the biochemical footprints and their benefits remain under-studied. Herein, using a combination of NMR metabolomics, microbial techniques and cell-culture assays, we have characterized metabolite profiles of probiotic viz. Lactobacillus delbruekii ATCC 9649, Lactobacillus casei ATCC 335, Lactobacillus plantarum NRC 716 and Bacillus coagulans ATCC 12425 cultures in fermented milk. We identified predominance of sugars, small chain fatty acids, organic acids and branched chain amino acids from natural abundance 13C NMR studies. Additionally, we identified myriad metabolites and their respective pathways using 1H NMR spectroscopy. Based on our findings, synbiotic fermented dairy products were customized with co-cultures and complemented with pro- and pre- biotics. Furthermore, we demonstrate epithelial cell interaction and anti-microbial activity of L. plantarum based ferment against a range of bacterial pathogens highlighting possible biochemical mechanisms for anti-microbial activity, quorum sensing, gut colonization and other beneficial factors that may be crucial. Furthermore, we propose plausible explanation against non-communicable diseases such as tumor-inhibitory, anti-proliferative and pro-apoptotic effects which has direct implications for dietary therapeutics.

Characterizing Watson-Crick versus Hoogsteen Base Pairing in a DNA-Protein Complex Using Nuclear Magnetic Resonance and Site-Specifically 13C- and 15N-Labeled DNA.

A( syn)-T and G( syn)-C+ Hoogsteen base pairs in protein-bound DNA duplexes can be difficult to resolve by X-ray crystallography due to ambiguous electron density and by nuclear magnetic resonance (NMR) spectroscopy due to poor chemical shift dispersion and size limitations with solution-state NMR spectroscopy. Here we describe an NMR strategy for characterizing Hoogsteen base pairs in protein-DNA complexes, which relies on site-specifically incorporating 13C- and 15N-labeled nucleotides into DNA duplexes for unambiguous resonance assignment and to improve spectral resolution. The approach was used to resolve the conformation of an A-T base pair in a crystal structure of an ∼43 kDa complex between a 34 bp duplex DNA and the integration host factor (IHF) protein. In the crystal structure (Protein Data Bank entry 1IHF ), this base pair adopts an unusual Hoogsteen conformation with a distorted sugar backbone that is accommodated by a nearby nick used to aid in crystallization. The NMR chemical shifts and interproton nuclear Overhauser effects indicate that this base pair predominantly adopts a Watson-Crick conformation in the intact DNA-IHF complex under solution conditions. Consistent with these NMR findings, substitution of 7-deazaadenine at this base pair resulted in only a small (∼2-fold) decrease in the IHF-DNA binding affinity. The NMR strategy provides a new approach for resolving crystallographic ambiguity and more generally for studying the structure and dynamics of protein-DNA complexes in solution.

High-performance virtual screening by targeting a high-resolution RNA dynamic ensemble.

Dynamic ensembles hold great promise in advancing RNA-targeted drug discovery. Here we subjected the transactivation response element (TAR) RNA from human immunodeficiency virus type-1 to experimental high-throughput screening against ~100,000 drug-like small molecules. Results were augmented with 170 known TAR-binding molecules and used to generate sublibraries optimized for evaluating enrichment when virtually screening a dynamic ensemble of TAR determined by combining NMR spectroscopy data and molecular dynamics simulations. Ensemble-based virtual screening scores molecules with an area under the receiver operator characteristic curve of ~0.85-0.94 and with ~40-75% of all hits falling within the top 2% of scored molecules. The enrichment decreased significantly for ensembles generated from the same molecular dynamics simulations without input NMR data and for other control ensembles. The results demonstrate that experimentally determined RNA ensembles can significantly enrich libraries with true hits and that the degree of enrichment is dependent on the accuracy of the ensemble.

Atomic structures of excited state A-T Hoogsteen base pairs in duplex DNA by combining NMR relaxation dispersion, mutagenesis, and chemical shift calculations.

NMR relaxation dispersion studies indicate that in canonical duplex DNA, Watson-Crick base pairs (bps) exist in dynamic equilibrium with short-lived low abundance excited state Hoogsteen bps. N1-methylated adenine (m1A) and guanine (m1G) are naturally occurring forms of damage that stabilize Hoogsteen bps in duplex DNA. NMR dynamic ensembles of DNA duplexes with m1A-T Hoogsteen bps reveal significant changes in sugar pucker and backbone angles in and around the Hoogsteen bp, as well as kinking of the duplex towards the major groove. Whether these structural changes also occur upon forming excited state Hoogsteen bps in unmodified duplexes remains to be established because prior relaxation dispersion probes provided limited information regarding the sugar-backbone conformation. Here, we demonstrate measurements of C3' and C4' spin relaxation in the rotating frame (R1ρ) in uniformly 13C/15N labeled DNA as sensitive probes of the sugar-backbone conformation in DNA excited states. The chemical shifts, combined with structure-based predictions using an automated fragmentation quantum mechanics/molecular mechanics method, show that the dynamic ensemble of DNA duplexes containing m1A-T Hoogsteen bps accurately model the excited state Hoogsteen conformation in two different sequence contexts. Formation of excited state A-T Hoogsteen bps is accompanied by changes in sugar-backbone conformation that allow the flipped syn adenine to form hydrogen-bonds with its partner thymine and this in turn results in overall kinking of the DNA toward the major groove. Results support the assignment of Hoogsteen bps as the excited state observed in canonical duplex DNA, provide an atomic view of DNA dynamics linked to formation of Hoogsteen bps, and lay the groundwork for a potentially general strategy for solving structures of nucleic acid excited states.

Dynamic basis for dG•dT misincorporation via tautomerization and ionization.

Tautomeric and anionic Watson-Crick-like mismatches have important roles in replication and translation errors through mechanisms that are not fully understood. Here, using NMR relaxation dispersion, we resolve a sequence-dependent kinetic network connecting G•T/U wobbles with three distinct Watson-Crick mismatches: two rapidly exchanging tautomeric species (Genol•T/UG•Tenol/Uenol; population less than 0.4%) and one anionic species (G•T-/U-; population around 0.001% at neutral pH). The sequence-dependent tautomerization or ionization step was inserted into a minimal kinetic mechanism for correct incorporation during replication after the initial binding of the nucleotide, leading to accurate predictions of the probability of dG•dT misincorporation across different polymerases and pH conditions and for a chemically modified nucleotide, and providing mechanisms for sequence-dependent misincorporation. Our results indicate that the energetic penalty for tautomerization and/or ionization accounts for an approximately 10-2 to 10-3-fold discrimination against misincorporation, which proceeds primarily via tautomeric dGenol•dT and dG•dTenol, with contributions from anionic dG•dT- dominant at pH 8.4 and above or for some mutagenic nucleotides.

Amiloride as a new RNA-binding scaffold with activity against HIV-1 TAR.

Diversification of RNA-targeted scaffolds offers great promise in the search for selective ligands of therapeutically relevant RNA such as HIV-1 TAR. We herein report the establishment of amiloride as a novel RNA-binding scaffold along with synthetic routes for combinatorial C(5)- and C(6)-diversification. Iterative modifications at the C(5)- and C(6)- positions yielded derivative 24, which demonstrated a 100-fold increase in activity over the parent dimethylamiloride in peptide displacement assays. NMR chemical shift mapping was performed using the 2D SOFAST- [1H-13C] HMQC NMR method, which allowed for facile and rapid evaluation of binding modes for all library members. Cheminformatic analysis revealed distinct differences between selective and non-selective ligands. In this study, we evolved dimethylamiloride from a weak TAR ligand to one of the tightest binding selective TAR ligands reported to date through a novel combination of synthetic methods and analytical techniques. We expect these methods to allow for rapid library expansion and tuning of the amiloride scaffold for a range of RNA targets and for SOFAST NMR to allow unprecedented evaluation of small molecule:RNA interactions.

Insights into Watson-Crick/Hoogsteen breathing dynamics and damage repair from the solution structure and dynamic ensemble of DNA duplexes containing m1A.

In the canonical DNA double helix, Watson-Crick (WC) base pairs (bps) exist in dynamic equilibrium with sparsely populated (∼0.02-0.4%) and short-lived (lifetimes ∼0.2-2.5 ms) Hoogsteen (HG) bps. To gain insights into transient HG bps, we used solution-state nuclear magnetic resonance spectroscopy, including measurements of residual dipolar couplings and molecular dynamics simulations, to examine how a single HG bp trapped using the N1-methylated adenine (m1A) lesion affects the structural and dynamic properties of two duplexes. The solution structure and dynamic ensembles of the duplexes reveals that in both cases, m1A forms a m1A•T HG bp, which is accompanied by local and global structural and dynamic perturbations in the double helix. These include a bias toward the BI backbone conformation; sugar repuckering, major-groove directed kinking (∼9°); and local melting of neighboring WC bps. These results provide atomic insights into WC/HG breathing dynamics in unmodified DNA duplexes as well as identify structural and dynamic signatures that could play roles in m1A recognition and repair.

m(1)A and m(1)G disrupt A-RNA structure through the intrinsic instability of Hoogsteen base pairs.

The B-DNA double helix can dynamically accommodate G-C and A-T base pairs in either Watson-Crick or Hoogsteen configurations. Here, we show that G-C(+) (in which + indicates protonation) and A-U Hoogsteen base pairs are strongly disfavored in A-RNA. As a result,N(1)-methyladenosine and N(1)-methylguanosine, which occur in DNA as a form of alkylation damage and in RNA as post-transcriptional modifications, have dramatically different consequences. Whereas they create G-C(+) and A-T Hoogsteen base pairs in duplex DNA, thereby maintaining the structural integrity of the double helix, they block base-pairing and induce local duplex melting in RNA. These observations provide a mechanism for disrupting RNA structure through post-transcriptional modifications. The different propensities to form Hoogsteen base pairs in B-DNA and A-RNA may help cells meet the opposing requirements of maintaining genome stability, on the one hand, and of dynamically modulating the structure of the epitranscriptome, on the other.

Characterizing RNA Excited States Using NMR Relaxation Dispersion.

Changes in RNA secondary structure play fundamental roles in the cellular functions of a growing number of noncoding RNAs. This chapter describes NMR-based approaches for characterizing microsecond-to-millisecond changes in RNA secondary structure that are directed toward short-lived and low-populated species often referred to as "excited states." Compared to larger scale changes in RNA secondary structure, transitions toward excited states do not require assistance from chaperones, are often orders of magnitude faster, and are localized to a small number of nearby base pairs in and around noncanonical motifs. Here, we describe a procedure for characterizing RNA excited states using off-resonance R1ρ NMR relaxation dispersion utilizing low-to-high spin-lock fields (25-3000 Hz). R1ρ NMR relaxation dispersion experiments are used to measure carbon and nitrogen chemical shifts in base and sugar moieties of the excited state. The chemical shift data are then interpreted with the aid of secondary structure prediction to infer potential excited states that feature alternative secondary structures. Candidate structures are then tested by using mutations, single-atom substitutions, or by changing physiochemical conditions, such as pH and temperature, to either stabilize or destabilize the candidate excited state. The resulting chemical shifts of the mutants or under different physiochemical conditions are then compared to those of the ground and excited states. Application is illustrated with a focus on the transactivation response element from the human immune deficiency virus type 1, which exists in dynamic equilibrium with at least two distinct excited states.

Discrete stacking of aromatic oligoamide macrocycles.

Unlike the precise structural control typical of closed assemblies, curbing the stacking of disc- and ring-shaped molecules is quite challenging. Here we report the discrete stacking of rigid aromatic oligoamide macrocycles 1. With increasing concentration, the aggregation of 1 quickly plateaus, forming a discrete oligomer, as suggested by 1D (1)H, 2D nuclear Overhauser effect, and diffusion-ordered NMR spectroscopy. Quantum-chemical calculations indicate that the tetramer of 1 is the most stable among oligomeric stacks. X-ray crystallography revealed a tetrameric stack containing identical molecules adopting two different conformations. With a defined length and an inner pore capable of accommodating distinctly different guests, the tetramers of 1 densely pack into 2D layers. Besides being a rare system of conformation-regulated supramolecular oligomerization, the discrete stacks of 1, along with their higher-order assemblies, may offer new nanotechnological applications.

New insights into Hoogsteen base pairs in DNA duplexes from a structure-based survey.

Hoogsteen (HG) base pairs (bps) provide an alternative pairing geometry to Watson-Crick (WC) bps and can play unique functional roles in duplex DNA. Here, we use structural features unique to HG bps (syn purine base, HG hydrogen bonds and constricted C1'-C1' distance across the bp) to search for HG bps in X-ray structures of DNA duplexes in the Protein Data Bank. The survey identifies 106 A•T and 34 G•C HG bps in DNA duplexes, many of which are undocumented in the literature. It also uncovers HG-like bps with syn purines lacking HG hydrogen bonds or constricted C1'-C1' distances that are analogous to conformations that have been proposed to populate the WC-to-HG transition pathway. The survey reveals HG preferences similar to those observed for transient HG bps in solution by nuclear magnetic resonance, including stronger preferences for A•T versus G•C bps, TA versus GG steps, and also suggests enrichment at terminal ends with a preference for 5'-purine. HG bps induce small local perturbations in neighboring bps and, surprisingly, a small but significant degree of DNA bending (∼14°) directed toward the major groove. The survey provides insights into the preferences and structural consequences of HG bps in duplex DNA.

Visualizing transient Watson-Crick-like mispairs in DNA and RNA duplexes.

Rare tautomeric and anionic nucleobases are believed to have fundamental biological roles, but their prevalence and functional importance has remained elusive because they exist transiently, in low abundance, and involve subtle movements of protons that are difficult to visualize. Using NMR relaxation dispersion, we show here that wobble dG•dT and rG•rU mispairs in DNA and RNA duplexes exist in dynamic equilibrium with short-lived, low-populated Watson-Crick-like mispairs that are stabilized by rare enolic or anionic bases. These mispairs can evade Watson-Crick fidelity checkpoints and form with probabilities (10(-3) to 10(-5)) that strongly imply a universal role in replication and translation errors. Our results indicate that rare tautomeric and anionic bases are widespread in nucleic acids, expanding their structural and functional complexity beyond that attainable with canonical bases.

Development and application of aromatic [(13)C, (1)H] SOFAST-HMQC NMR experiment for nucleic acids.

Higher sensitivity of NMR spectrometers and novel isotopic labeling schemes have ushered the development of rapid data acquisition methodologies, improving the time resolution with which NMR data can be acquired. For nucleic acids, longitudinal relaxation optimization in conjunction with Ernst angle excitation (SOFAST-HMQC) for imino protons, in addition to rendering rapid pulsing, has been demonstrated to yield significant improvements in sensitivity per unit time. Extending such methodology to other spins offers a viable prospect to measure additional chemical shifts, thereby broadening their utilization for various applications. Here, we introduce the 2D [(13)C, (1)H] aromatic SOFAST-HMQC that results in overall sensitivity gain of 1.4- to 1.7-fold relative to the conventional HMQC and can also be extended to yield long-range heteronuclear chemical shifts such as the adenine imino nitrogens N1, N3, N7 and N9. The applications of these experiments range from monitoring real-time biochemical processes, drug/ligand screening, and to collecting data at very low sample concentration and/or in cases where isotopic enrichment cannot be achieved.