Protein-Ligand Interactions in the STING Binding Site Probed by Rationally Designed Single-Point Mutations: Experiment and Theory

Unraveling the Binding Mode of Cyclic Adenosine-Inosine Monophosphate (cAIMP) to STING through Molecular Dynamics Simulations

The stimulator of interferon genes (STING) plays a significant role in immune defense and protection against tumor proliferation. Many cyclic dinucleotide (CDN) analogues have been reported to regulate its activity, but the dynamic process involved when the ligands activate STING remains unclear. In this work, all-atom molecular dynamics simulations were performed to explore the binding mode between human STING (hSTING) and four cyclic adenosine-inosine monophosphate analogs (cAIMPs), as well as 2',3'-cGMP-AMP (2',3'-cGAMP). The results indicate that these cAIMPs adopt a U-shaped configuration within the binding pocket, forming extensive non-covalent interaction networks with hSTING. These interactions play a significant role in augmenting the binding, particularly in interactions with Tyr167, Arg238, Thr263, and Thr267. Additionally, the presence of hydrophobic interactions between the ligand and the receptor further contributes to the overall stability of the binding. In this work, the conformational changes in hSTING upon binding these cAIMPs were also studied and a significant tendency for hSTING to shift from open to closed state was observed after binding some of the cAIMP ligands.

Fluorinated cGAMP analogs, which act as STING agonists and are not cleavable by poxins: Structural basis of their function

The cGAS-STING pathway is a crucial part of innate immunity; it serves to detect DNA in the cytoplasm and to defend against certain cancers, viruses, and bacteria. We designed and synthesized fluorinated carbocyclic cGAMP analogs, MD1203 and MD1202D (MDs), to enhance their stability and their affinity for STING. These compounds demonstrated exceptional activity against STING. Despite their distinct chemical modifications relative to the canonical cyclic dinucleotides (CDNs), crystallographic analysis revealed a binding mode with STING that was consistent with the canonical CDNs. Importantly, MDs were resistant to cleavage by viral poxin nucleases and MDs-bound poxin adopted an unliganded-like conformation. Moreover, MDs complexed with poxin showed a conformation distinct from cGAMP bound to poxin, closely resembling their conformation when bound to STING. In conclusion, the development of MD1203 and MD1202D showcases their potential as potent STING activators with remarkable stability against poxin-mediated degradation-a crucial characteristic for future development of antivirals.

Vinylphosphonate-based cyclic dinucleotides enhance STING-mediated cancer immunotherapy

Cyclic dinucleotides (CDNs) trigger the cyclic GMP-AMP synthase-stimulator of interferon genes (cGAS-STING) pathway, which plays a key role in cytosolic DNA sensing and thus in immunomodulation against infections, cell damage and cancer. However, cancer immunotherapy trials with CDNs have shown immune activation, but not complete tumor regression. Nevertheless, we designed a novel class of CDNs containing vinylphosphonate based on a STING-affinity screening assay. In vitro, acyloxymethyl phosphate/phosphonate prodrugs of these vinylphosphonate CDNs were up to 1000-fold more potent than the clinical candidate ADU-S100. In vivo, the lead prodrug induced tumor-specific T cell priming and facilitated tumor regression in the 4T1 syngeneic mouse model of breast cancer. Moreover, we solved the crystal structure of this ligand bound to the STING protein. Therefore, our findings not only validate the therapeutic potential of vinylphosphonate CDNs but also open up opportunities for drug development in cancer immunotherapy bridging innate and adaptive immunity.

Recent trends in STING modulators: Structures, mechanisms, and therapeutic potential

The cyclic GMP-AMP synthase stimulator (cGAS)-stimulator of interferon gene (STING) signaling pathway has an integral role in the host immune response through DNA sensing followed by inducing a robust innate immune defense program. STING has become a promising therapeutic target associated with multiple diseases, including various inflammatory diseases, cancer, and infectious diseases, among others. Thus, modulators of STING are regarded as emerging therapeutic agents. Recent progress has been made in STING research, including recently identified STING-mediated regulatory pathways, the development of a new STING modulator, and the new association of STING with disease. In this review, we focus on recent trends in the development of STING modulators, including structures, mechanisms, and clinical application.

Discovery of Novel STING Inhibitors Based on the Structure of the Mouse STING Agonist DMXAA

The stimulator-of-interferon-gene (STING) protein is involved in innate immunity. The drug DMXAA (5,6-dimethylxanthenone-4-acetic acid) proved to be a potent murine-STING (mSTING) agonist but had little effect on human-STING (hSTING). In this paper, we draw upon the comparison of different crystal structures and protein-ligand interaction relationships analysis to venture the hypothesis that the drug design of DMXAA variants has the potential to convert STING agonists to inhibitors. Based on our previous discovery of two DMXAA analogs, 3 and 4 (both could bind to STING), we structurally optimized them and synthesized new derivatives, respectively. In binding assays, we found compounds 11 and 27 to represent STING binders that were superior to the original structures and discussed the structure-activity relationships. All target compounds were inactive in cellular assays for the screening of STING agonistic activity. Gratifyingly, we identified 11 and 27 as STING inhibitors with micromolar activity in both hSTING and mSTING pathways. In addition, 11 and 27 inhibited the induction of interferon and inflammatory cytokines activated by 2'3'-cGAMP without apparent cytotoxicity. These findings break the rigid thinking that DMXAA provides the structural basis specifically for STING agonists and open up more possibilities for developing novel STING agonists or inhibitors.

Molecular dynamics simulations provide structural insight into binding of cyclic dinucleotides to human STING protein

Human stimulator of interferon genes (hSTING) is a signaling adaptor protein that triggers innate immune system by response to cytosolic DNA and second messenger cyclic dinucleotides (CDNs). Natural CDNs contain purine nucleobase with different phosphodiester linkage types (3'-3', 2'-2' or mixed 2'-3'-linkages) and exhibit different binding affinity towards hSTING, ranging from micromolar to nanomolar. High-affinity CDNs are considered as suitable candidates for treatment of chronic hepatitis B and cancer. We have used molecular dynamics simulations to investigate dynamical aspects of binding of natural CDNs (specifically, 2'-2'-cGAMP, 2'-3'-cGAMP, 3'-3'-cGAMP, 3'-3'-c-di-AMP, and 3'-3'-c-di-GMP) with hSTING^wt protein. Our results revealed that CDN/hSTING^wt interactions are controlled by the balance between fluctuations (conformational changes) in the CDN ligand and the protein dynamics. Binding of different CDNs induces different degrees of conformational/dynamics changes in hSTING^wt ligand binding cavity, especially in α₁-helices, the so-called lid region and α₂-tails. The ligand residence time in hSTING^wt protein pocket depends on different contribution of R232 and R238 residues interacting with oxygen atoms of phosphodiester groups in ligand, water distribution around interacting charged centers (in protein residues and ligand) and structural stability of closed conformation state of hSTING^wt protein. These findings may perhaps guide design of new compounds modulating hSTING activity.Communicated by Ramaswamy H. Sarma.

Design, Synthesis, and Biochemical and Biological Evaluation of Novel 7-Deazapurine Cyclic Dinucleotide Analogues as STING Receptor Agonists

Cyclic dinucleotides (CDNs) are second messengers that activate stimulator of interferon genes (STING). The cGAS-STING pathway plays a promising role in cancer immunotherapy. Here, we describe the synthesis of CDNs containing 7-substituted 7-deazapurine moiety. We used mouse cyclic GMP-AMP synthase and bacterial dinucleotide synthases for the enzymatic synthesis of CDNs. Alternatively, 7-(het)aryl 7-deazapurine CDNs were prepared by Suzuki-Miyaura cross-couplings. New CDNs were tested in biochemical and cell-based assays for their affinity to human STING. Eight CDNs showed better activity than 2'3'-cGAMP, the natural ligand of STING. The effect on cytokine and chemokine induction was also evaluated. The best activities were observed for CDNs bearing large aromatic substituents that point above the CDN molecule. We solved four X-ray structures of complexes of new CDNs with human STING. We observed π-π stacking interactions between the aromatic substituents and Tyr240 that are involved in the stabilization of CDN-STING complexes.

African Swine Fever Virus EP364R and C129R Target Cyclic GMP-AMP To Inhibit the cGAS-STING Signaling Pathway

African swine fever virus (ASFV) is a highly pathogenic swine DNA virus with high mortality that causes African swine fever (ASF) in domestic pigs and wild boars. For efficient viral infection, ASFV has developed complex strategies to evade key components of antiviral innate immune responses. However, the immune escape mechanism of ASFV remains unclear. Upon ASFV infection, cyclic GMP-AMP (2′,3′-cGAMP) synthase (cGAS), a cytosolic DNA sensor, recognizes ASFV DNA and synthesizes the second messenger 2′,3′-cGAMP, which triggers interferon (IFN) production to interfere with viral replication. In this study, we demonstrated a novel immune evasion mechanism of ASFV EP364R and C129R, which blocks cellular cyclic 2′,3′-cGAMP-mediated antiviral responses. ASFV EP364R and C129R with nuclease homology inhibit IFN-mediated responses by specifically interacting with 2′,3′-cGAMP and exerting their phosphodiesterase (PDE) activity to cleave 2′,3′-cGAMP. Particularly notable is that ASFV EP364R had a region of homology with the stimulator of interferon genes (STING) protein containing a 2′,3′-cGAMP-binding motif and point mutations in the Y76S and N78A amino acids of EP364R that impaired interaction with 2′,3′-cGAMP and restored subsequent antiviral responses. These results highlight a critical role for ASFV EP364R and C129R in the inhibition of IFN responses and could be used to develop ASFV live attenuated vaccines.

IMPORTANCE African swine fever (ASF) is a highly contagious hemorrhagic disease in domestic pigs and wild boars caused by African swine fever virus (ASFV). ASF is a deadly epidemic disease in the global pig industry, but no drugs or vaccines are available. Understanding the pathogenesis of ASFV is essential to developing an effective live attenuated ASFV vaccine, and investigating the immune evasion mechanisms of ASFV is crucial to improve the understanding of its pathogenesis. In this study, for the first time, we identified the EP364R and C129R, uncharacterized proteins that inhibit type I interferon signaling. ASFV EP364R and C129R specifically interacted with 2′,3′-cGAMP, the mammalian second messenger, and exerted phosphodiesterase activity to cleave 2′,3′-cGAMP. In this study, we discovered a novel mechanism by which ASFV inhibits IFN-mediated antiviral responses, and our findings can guide the understanding of ASFV pathogenesis and the development of live attenuated ASFV vaccines.

How Fragile We Are: Influence of Stimulator of Interferon Genes (STING) Variants on Pathogen Recognition and Immune Response Efficiency

The stimulator of interferon genes (STING) protein is a cornerstone of the human immune response. Its activation by cGAMP in the presence of cytosolic DNA stimulates the production of type I interferons and inflammatory cytokines. In the human population, several STING variants exist and exhibit dramatic differences in their activity, impacting the efficiency of the host defense against infections. Understanding the molecular mechanisms of these variants opens perspectives for personalized medicine treatments against diseases such as viral infections, cancers, or autoinflammatory diseases. Through microsecond-scale molecular modeling simulations, contact analyses, and machine learning techniques, we reveal the dynamic behavior of four STING variants (wild type, G230A, R293Q, and G230A/R293Q) and rationalize the variability of efficiency observed experimentally. Our results show that the decrease in STING activity is linked to a stiffening of key structural elements of the binding cavity together with changes in the interaction patterns within the protein.

Discovery of isonucleotidic CDNs as potent STING agonists with immunomodulatory potential

Stimulator of interferon genes (STING) is an adaptor protein of the cGAS-STING signaling pathway involved in the sensing of cytosolic DNA. It functions as a receptor for cyclic dinucleotides (CDNs) and, upon their binding, mediates cytokine expression and host immunity. Besides naturally occurring CDNs, various synthetic CDNs, such as ADU-S100, have been reported to effectively activate STING and are being evaluated in clinical trials for the treatment of cancer. Here, we describe the preparation of a unique new class of STING agonists: isonucleotidic cyclic dinucleotides and the synthesis of their prodrugs. The presented CDNs stimulate STING with comparable efficiency to ADU-S100, whereas their prodrugs demonstrate activity up to four orders of magnitude better due to the improved cellular uptake. The compounds are very potent inducers of inflammatory cytokines by peripheral blood mononuclear cells (PBMCs). We also report the X-ray crystal structure of the lead inhibitor bound to the wild-type (WT) STING.

Predicting Effects of Site-Directed Mutagenesis on Enzyme Kinetics by QM/MM and QM Calculations: A Case of Glutamate Carboxypeptidase II

Quantum and molecular mechanics (QM/MM) and QM-only (cluster model) modeling techniques represent the two workhorses in mechanistic understanding of enzyme catalysis. One of the stringent tests for QM/MM and/or QM approaches is to provide quantitative answers to real-world biochemical questions, such as the effect of single-point mutations on enzyme kinetics. This translates into predicting the relative activation energies to 1-2 kcal·mol^-1 accuracy; such predictions can be used for the rational design of novel enzyme variants with desired/improved characteristics. Herein, we employ glutamate carboxypeptidase II (GCPII), a dizinc metallopeptidase, also known as the prostate specific membrane antigen, as a model system. The structure and activity of this major cancer antigen have been thoroughly studied, both experimentally and computationally, which makes it an ideal model system for method development. Its reaction mechanism is quite well understood: the reaction coordinate comprises a "tetrahedral intermediate" and two transition states and experimental activation Gibbs free energy of ∼17.5 kcal·mol^-1 can be inferred for the known k_cat ≈ 1 s^-1. We correlate experimental kinetic data (including the E424H variant, newly characterized in this work) for various GCPII mutants (k_cat = 8.6 × 10^-5 s^-1 to 2.7 s^-1) with the energy profiles calculated by QM/MM and QM-only (cluster model) approaches. We show that the near-quantitative agreement between the experimental values and the calculated activation energies (ΔH^⧧) can be obtained and recommend the combination of the two protocols: QM/MM optimized structures and cluster model (QM) energetics. The trend in relative activation energies is mostly independent of the QM method (DFT functional) used. Last but not least, a satisfactory correlation between experimental and theoretical data allows us to provide qualitative and fairly simple explanations of the observed kinetic effects which are thus based on a rigorous footing.

Enzymatic Synthesis of 3'-5', 3'-5' Cyclic Dinucleotides, Their Binding Properties to the Stimulator of Interferon Genes Adaptor Protein, and Structure/Activity Correlations

The 3'-5', 3'-5' cyclic dinucleotides (3'3'CDNs) are bacterial second messengers that can also bind to the stimulator of interferon genes (STING) adaptor protein in vertebrates and activate the host innate immunity. Here, we profiled the substrate specificity of four bacterial dinucleotide synthases from Vibrio cholerae (DncV), Bacillus thuringiensis (btDisA), Escherichia coli (dgcZ), and Thermotoga maritima (tDGC) using a library of 33 nucleoside-5'-triphosphate analogues and then employed these enzymes to synthesize 24 3'3'CDNs. The STING affinity of CDNs was evaluated in cell-based and biochemical assays, and their ability to induce cytokines was determined by employing human peripheral blood mononuclear cells. Interestingly, the prepared heterodimeric 3'3'CDNs bound to the STING much better than their homodimeric counterparts and showed similar or better potency than bacterial 3'3'CDNs. We also rationalized the experimental findings by in-depth STING-CDN structure-activity correlations by dissecting computed interaction free energies into a set of well-defined and intuitive terms. To this aim, we employed state-of-the-art methods of computational chemistry, such as quantum mechanics/molecular mechanics (QM/MM) calculations, and complemented the computed results with the {STING:3'3'c-di-ara-AMP} X-ray crystallographic structure. QM/MM identified three outliers (mostly homodimers) for which we have no clear explanation of their impaired binding with respect to their heterodimeric counterparts, whereas the R² = 0.7 correlation between the computed ΔG'_{int_rel} and experimental ΔT_m's for the remaining ligands has been very encouraging.

Synthesis and Biological Evaluation of Phosphoester and Phosphorothioate Prodrugs of STING Agonist 3',3'-c-Di(2'F,2'dAMP)

Cyclic dinucleotides (CDNs) are second messengers that bind to the stimulator of interferon genes (STING) and trigger the expression of type I interferons and proinflammatory cytokines. Here we evaluate the activity of 3',3'-c-di(2'F,2'dAMP) and its phosphorothioate analogues against five STING allelic forms in reporter-cell-based assays and rationalize our findings with X-ray crystallography and quantum mechanics/molecular mechanics calculations. We show that the presence of fluorine in the 2' position of 3',3'-c-di(2'F,2'dAMP) improves its activity not only against the wild type (WT) but also against REF and Q STING. Additionally, we describe the synthesis of the acyloxymethyl and isopropyloxycarbonyl phosphoester prodrugs of CDNs. Masking the negative charges of the CDNs results in an up to a 1000-fold improvement of the activities of the prodrugs relative to those of their parent CDNs. Finally, the uptake and intracellular cleavage of pivaloyloxymethyl prodrugs to the parent CDN is rapid, reaching a peak intracellular concentration within 2 h.

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.