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 hSTINGwt protein. Our results revealed that CDN/hSTINGwt 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 hSTINGwt ligand binding cavity, especially in α1-helices, the so-called lid region and α2-tails. The ligand residence time in hSTINGwt 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 hSTINGwt protein. These findings may perhaps guide design of new compounds modulating hSTING activity.Communicated by Ramaswamy H. Sarma.

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.

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

STING protein (stimulator of interferon genes) plays an important role in the innate immune system. A number of potent compounds regulating its activity have been reported, mostly derivatives of cyclic dinucleotides (CDNs), natural STING agonists. Here, we aim to provide complementary information to large-scale "ligand-profiling" studies by probing the importance of STING-CDN protein-ligand interactions on the protein side. We examined in detail six typical CDNs each in complex with 13 rationally devised mutations in STING: S162A, S162T, Y167F, G230A, R232K, R232H, A233L, A233I, R238K, T263A, T263S, R293Q, and G230A/R293Q. The mutations switch on and off various types of protein-ligand interactions: π-π stacking, hydrogen bonding, ionic pairing, and nonpolar contacts. We correlated experimental data obtained by differential scanning fluorimetry, X-ray crystallography, and isothermal titration calorimetry with theoretical calculations. This enabled us to provide a mechanistic interpretation of the differences in the binding of representative CDNs to STING. We observed that the G230A mutation increased the thermal stability of the protein-ligand complex, indicating an increased level of ligand binding, whereas R238K and Y167F led to a complete loss of stabilization (ligand binding). The effects of the other mutations depended on the type of ligand (CDN) and varied, to some extent. A very good correlation (R2 = 0.6) between the experimental binding affinities and interaction energies computed by quantum chemical methods enabled us to explain the effect of the studied mutations in detail and evaluate specific interactions quantitatively. Our work may inspire development of high-affinity ligands against the common STING haplotypes by targeting the key (sometimes non-intuitive) protein-ligand interactions.

Halides with Fifteen Aliphatic C-H···Anion Interaction Sites.

Since the aliphatic C-H···anion interaction is relatively weak, anion binding using hydrophobic aliphatic C-H (Cali-H) groups has generally been considered not possible without the presence of additional binding sites that contain stronger interactions to the anion. Herein, we report X-ray structures of organic crystals that feature a chloride anion bound exclusively by hydrophobic Cali-H groups. An X-ray structure of imidazolium-based scaffolds using Cali-H···A(-) interactions (A(-) = anion) shows that a halide anion is directly interacting with fifteen Cali-H groups (involving eleven hydrogen bonds, two bidentate hydrogen-bond-type binding interactions and two weakly hydrogen-bonding-like binding interactions). Additional supporting interactions and/or other binding sites are not observed. We note that such types of complexes may not be rare since such high numbers of binding sites for an anion are also found in analogous tetraalkylammonium complexes. The Cali-H···A(-) interactions are driven by the formation of a near-spherical dipole layer shell structure around the anion. The alternating layers of electrostatic charge around the anion arise because the repulsions between weakly positively charged H atoms are reduced by the presence of the weakly negatively charged C atoms connected to H atoms.

Intriguing Electrostatic Potential of CO: Negative Bond-ends and Positive Bond-cylindrical-surface.

The strong electronegativity of O dictates that the ground state of singlet CO has positively charged C and negatively charged O, in agreement with ab initio charge analysis, but in disagreement with the dipole direction. Though this unusual phenomenon has been fairly studied, the study of electrostatic potential (EP) for noncovalent interactions of CO is essential for better understanding. Here we illustrate that both C and O atom-ends show negative EP (where the C end gives more negative EP), favoring positively charged species, whereas the cylindrical surface of the CO bond shows positive EP, favoring negatively charged ones. This is demonstrated from the interactions of CO with Na(+), Cl(-), H2O, CO and benzene. It can be explained by the quadrupole driven electrostatic nature of CO (like N2) with very weak dipole moment. The EP is properly described by the tripole model taking into account the electrostatic multipole moments, which has a large negative charge at a certain distance protruded from C, a large positive charge on C, and a small negative charge on O. We also discuss the EP of the first excited triplet CO.

Do coinage metal anions interact with substituted benzene derivatives?

The nature of the anion-π interaction has been investigated by carrying out ab initio calculations of the complexes of coinage metal anions (Au(-), Ag(-), and Cu(-)) with different kinds of π-systems. The binding energies indicate that gold anion has the highest and copper anion has the lowest affinity for interactions with π-systems. Different aspects of the anion-π interaction in these systems have been investigated, including charge-transfer effects (using the Merz-Kollman method), "atoms-in-molecules" (AIM) topological parameters, and interaction energies (using energy decomposition analysis, EDA). Our results indicated that, for most M(-)···π interactions, the electrostatic term provides the dominant contribution, whereas polarization, charge transfer, and dispersion effects contribute less than 25 % of the interaction. We believe that the present results should lead to a greater understanding of the basis for anion-π interactions of coinage metal anions.

Comparison of gas phase intrinsic properties of cytosine and thymine nucleobases with their O-alkyl adducts: different hydrogen bonding preferences for thymine versus O-alkyl thymine.

In recent years, there has been increasing interest in damaged DNA and RNA nucleobases. These damaged nucleobases can cause DNA mutation, resulting in various diseases such as cancer. Alkylating agents are mutagenic and carcinogenic in a variety of prokaryotic and eukaryotic organisms. The present study employs density functional theory (DFT/B3LYP) with the 6-311++G(d,p) basis set to investigate the effect of chemical damage in O-alkyl pyrimidines such as O(4)-methylthymine, O(2)-methylcytosine and O(2)-methylthymine. We compared the intrinsic properties, such as proton affinities, gas phase acidities, equilibrium tautomerization and nucleobase pair's hydrogen bonding properties, of these molecules with those in the normal nucleobases thymine and cytosine. The results are of interest for chemical reasons and also possibly for biological purposes since biological media can be quite non-polar. Furthermore, we found that N1-H of O(4)-methylthymine is less acidic than N1-H of thymine, suggesting that alkyl DNA glycosylase enzyme cannot discriminate this damaged nucleobase from a normal thymine nucleobase. This result indicates that the conjugated base anion of O(4)-methylthymine would be a worse leaving group and O(4)-methylthymine is repaired in genome by demethylation rather than enzyme-catalyzed excision at N1.

Interactions of glutathione tripeptide with gold cluster: influence of intramolecular hydrogen bond on complexation behavior.

Understanding the nature of the interaction between metal nanoparticles and biomolecules has been important in the development and design of sensors. In this paper, structural, electronic, and bonding properties of the neutral and anionic forms of glutathione tripeptide (GSH) complexes with a Au(3) cluster were studied using the DFT-B3LYP with 6-31+G**-LANL2DZ mixed basis set. Binding of glutathione with the gold cluster is governed by two different kinds of interactions: Au-X (X = N, O, and S) anchoring bond and Au···H-X nonconventional hydrogen bonding. The influence of the intramolecular hydrogen bonding of glutathione on the interaction of this peptide with the gold cluster has been investigated. To gain insight on the role of intramolecular hydrogen bonding on Au-GSH interaction, we compared interaction energies of Au-GSH complexes with those of cystein and glycine components. Our results demonstrated that, in spite of the ability of cystein to form highly stable metal-sulfide interaction, complexation behavior of glutathione is governed by its intramolecular backbone hydrogen bonding. The quantum theory of atom in molecule (QTAIM) and natural bond orbital analysis (NBO) have also been applied to interpret the nature of interactions in Au-GSH complexes. Finally, conformational flexibility of glutathione during complexation with the Au(3) cluster was investigated by means of monitoring Ramachandran angles.