Computational and NMR spectroscopy insights into the conformation of cyclic di-nucleotides

Solvent Dependent Nuclear Magnetic Resonance Molecular Parameters Based on a Polarization Consistent Screened Range Separated Hybrid Density Functional Theory Framework

Nuclear magnetic resonance (NMR) properties of solvated molecules are significantly affected by the solvent. We, therefore, employ a polarization consistent framework that efficiently addresses the solvent polarizing environment effects. Toward this goal a dielectric screened range separated hybrid (SRSH) functional is invoked with a polarizable continuum model (PCM) to properly represent the orbital gap in the condensed phase. We build on the success of range separated hybrid (RSH) functionals to address the erroneous tendency of traditional density functional theory (DFT) to collapse the orbital gap. Recently, the impact of RSH that properly opens up the orbital gap in gas-phase calculations on NMR properties has been assessed. Here, we report the use of SRSH-PCM that produces properly solute orbital gaps in calculating isotropic nuclear magnetic shielding and chemical shift parameters of molecular systems in the condensed phase. We show that in contrast to simpler DFT-PCM approaches, SRSH-PCM successfully follows expected dielectric constant trends.

2',3'-Cyclic GMP-AMP Dinucleotides for STING-Mediated Immune Modulation: Principles, Immunotherapeutic Potential, and Synthesis

The cGAS-STING pathway discovered ten years ago is an important component of the innate immune system. Activation of cGAS-STING triggers downstream signalling, such as TBK1-IRF3, NF-κB and autophagy, which in turn leads to antipathogen responses, durable antitumour immunity or autoimmune diseases. 2',3'-Cyclic GMP-AMP dinucleotides (2',3'-cGAMP), the key second messengers produced by cGAS, play a pivotal role in cGAS-STING signalling by binding and activating STING. Thus, 2',3'-cGAMP has immunotherapeutic potential, which in turn has stimulated research on the design and synthesis of 2',3'-cGAMP analogues for clinical applications over the past ten years. This review presents the discovery, metabolism, and function of 2',3'-cGAMP in the cGAS-STING innate immune signalling axis. The enzymatic and chemical syntheses of 2',3'-cGAMP analogues as STING-targeting therapeutics are also summarized.

Conformational energies and equilibria of cyclic dinucleotides in vacuo and in solution: computational chemistry vs. NMR experiments

Performance of computational methods in modelling cyclic dinucleotides - an important and challenging class of compounds - has been evaluated by two different benchmarks: (1) gas-phase conformational energies and (2) qualitative agreement with NMR observations of the orientation of the χ-dihedral angle in solvent. In gas-phase benchmarks, where CCSD(T) and DLPNO-CCSD(T) methods have been used as the reference, most of the (dispersion corrected) density functional approximations are accurate enough to justify prioritizing computational cost and compatibility with other modelling options as the criterion of choice. NMR experiments of 3'3'-c-di-AMP, 3'3'-c-GAMP, and 3'3'-c-di-GMP show the overall prevalence of the anti-conformation of purine bases, but some population of syn-conformations is observed for guanines. Implicit solvation models combined with quantum-chemical methods struggle to reproduce this behaviour, probably due to a lack of dynamics and explicitly modelled solvent, leading to structures that are too compact. Molecular dynamics simulations overrepresent the syn-conformation of guanine due to the overestimation of an intramolecular hydrogen bond. Our combination of experimental and computational benchmarks provides "error bars" for modelling cyclic dinucleotides in solvent, where such information is generally difficult to obtain, and should help gauge the interpretability of studies dealing with binding of cyclic dinucleotides to important pharmaceutical targets. At the same time, the presented analysis calls for improvement in both implicit solvation models and force-field parameters.

Synthesis and biological evaluation of all possible inosine-mixed cyclic dinucleotides that activate different hSTING variants

Cyclic dinucleotides (CDNs) could activate stimulator of interferon genes (STING) protein to produce type I interferon and other pro-inflammation cytokines in mammalian cells. To explore new types of potentially efficient STING activators targeting all five major hSTING variants (WT, R232H, HAQ, AQ and R293Q), we here reported the synthesis of a total of 19 inosine-containing CDNs based on the combinations of hypoxanthine with four natural bases (A, G, C and U) and three phosphodiester linkage backbones (3'-3', 2'-3', 2'-2'). The IFN-β induction results showed that all of the 2'-3' and 2'-2' CDNs linked by inosine and purine nucleosides favored the stacking interaction with Y167 and R238 residues of hSTING protein, and several CDNs constructed by hypoxanthine and pyrimidine like c[I(2',5')U(2',5')] could also activate all five hSTING variants. The molecular dynamic simulation and the isothermal titration calorimetric (ITC) assay further demonstrated the potential of cAIMP isomers with 2'-5' phosphate to form the hydrogen binding with R232 and R238 residues of hSTING in an entropically driven manner compared to cGAMP isomers. It would be promising to exploit novel inosine-mixed CDNs as activators of hSTING variants in immune therapy.

Structural Analysis of (p)ppGpp Reveals Its Versatile Binding Pattern for Diverse Types of Target Proteins

(p)ppGpp, highly phosphorylated guanosine, are global regulatory nucleotides that modulate several biochemical events in bacterial physiology ranging from core central dogma to various metabolic pathways. Conventionally, (p)ppGpp collectively refers to two nucleotides, ppGpp, and pppGpp in the literature. Initially, (p)ppGpp has been discovered as a transcription regulatory molecule as it binds to RNA polymerase and regulates transcriptional gene regulation. During the past decade, several other target proteins of (p)ppGpp have been discovered and as of now, more than 30 proteins have been reported to be regulated by the binding of these two signaling nucleotides. The regulation of diverse biochemical activities by (p)ppGpp requires fine-tuned molecular interactions with various classes of proteins so that it can moderate varied functions. Here we report a structural dynamics of (p)ppGpp in the unbound state using well-defined computational tools and its interactions with target proteins to understand the differential regulation by (p)ppGpp at the molecular level. We carried out replica exchange molecular dynamics simulation studies to enhance sampling of conformations during (p)ppGpp simulation. The detailed comparative analysis of torsion angle conformation of ribose sugar of unbound (p)ppGpp and bound states of (p)ppGpp was carried out. The structural dynamics shows that two linear phosphate chains provide plasticity to (p)ppGpp nucleotides for the binding to diverse proteins. Moreover, the intermolecular interactions between (p)ppGpp and target proteins were characterized through various physicochemical parameters including, hydrogen bonds, van der Waal’s interactions, aromatic stacking, and side chains of interacting residues of proteins. Surprisingly, we observed that interactions of (p)ppGpp to target protein have a consensus binding pattern for a particular functional class of enzymes. For example, the binding of (p)ppGpp to RNA polymerase is significantly different from the binding of (p)ppGpp to the proteins involved in the ribosome biogenesis pathway. Whereas, (p)ppGpp binding to enzymes involved in nucleotide metabolism facilitates the functional regulation through oligomerization. Analysis of these datasets revealed that guanine base-specific contacts are key determinants to discriminate functional class of protein. Altogether, our studies provide significant information to understand the differential interaction pattern of (p)ppGpp to its target and this information may be useful to design antibacterial compounds based on (p)ppGpp analogs.

Highly Efficient Cyclic Dinucleotide Based Artificial Metalloribozymes for Enantioselective Friedel-Crafts Reactions in Water

The diverse secondary structures of nucleic acids are emerging as attractive chiral scaffolds to construct artificial metalloenzymes (ArMs) for enantioselective catalysis. DNA-based ArMs containing duplex and G-quadruplex scaffolds have been widely investigated, yet RNA-based ArMs are scarce. Here we report that a cyclic dinucleotide of c-di-AMP and Cu²⁺ ions assemble into an artificial metalloribozyme (c-di-AMP⋅Cu²⁺ ) that enables catalysis of enantioselective Friedel-Crafts reactions in aqueous media with high reactivity and excellent enantioselectivity of up to 97 % ee. The assembly of c-di-AMP⋅Cu²⁺ gives rise to a 20-fold rate acceleration compared to Cu²⁺ ions. Based on various biophysical techniques and density function theory (DFT) calculations, a fine coordination structure of c-di-AMP⋅Cu²⁺ metalloribozyme is suggested in which two c-di-AMP form a dimer scaffold and the Cu²⁺ ion is located in the center of an adenine-adenine plane through binding to two N7 nitrogen atoms and one phosphate oxygen atom.

Solution Structures of a G-Quadruplex Bound to Linear- and Cyclic-Dinucleotides

Cyclic dinucleotides have emerged as important secondary messengers and cell signaling molecules that regulate several cell responses. A guanine-deficit G-quadruplex structure formation by a sequence containing (4n - 1) guanines, n denoting the number of G-tetrad layers, was previously reported. Here, a (4n - 1) G-quadruplex structure is shown to be capable of binding guanine-containing dinucleotides in micromolar affinity. The guanine base of the dinucleotides interacts with a vacant G-triad, forming four additional Hoogsteen hydrogen bonds to complete a G-tetrad. Solution structures of two complexes, both comprised of a (4n - 1) G-quadruplex structure, one bound to a linear dinucleotide (d(AG)) and the other to a cyclic dinucleotide (cGAMP), are solved using NMR spectroscopy. The latter suggests sufficiently strong interaction between the guanine base of the dinucleotide and the vacant G-triad, which acts as an anchor point of binding. The binding interfaces from the two solution structures provide useful information for specific ligand design. The results also infer that other guanine-containing metabolites of a similar size have the capability of binding G-quadruplexes, potentially affecting the expression of the metabolites and functionality of the bound G-quadruplexes.