INTRODUCTION {#SEC1}
============

Cyclic dinucleotides are second messengers found ubiquitously in prokaryotes and as part of the innate immune system in eukaryotes ([@B1]).  Three cyclic dinucleotides have been identified in bacteria (Figure [1A](#F1){ref-type="fig"})--cyclic-di-GMP (c-di-GMP), cyclic-di-AMP (c-di-AMP), and 3′,5′-cyclic-AMP-GMP (c-AMP-GMP)---and one cyclic dinucleotide in higher order eukaryotes---2′,5′-cyclic-AMP-GMP (2′,5′-cGAMP) ([@B2],[@B3]).  In bacteria, the levels of each of these cyclic dinucleotides fluctuate in response to external stimuli and result in lifestyle changes essential for bacterial survival and host infection.  In general, differing levels of c-di-GMP are associated with biofilm formation and virulence factor production ([@B4]).  c-di-AMP regulates potassium homeostasis, osmolarity, DNA integrity and biofilm formation ([@B9]).  c-AMP-GMP levels are associated with changes in chemotaxis in several species of bacteria and exoelectrogenesis in *Geobacter* ([@B14],[@B15]). 

![ Mechanism of DncV catalysis. (** A**) Chemical structure of three known natural cyclic dinucleotides.  The adenine base is indicated by green and the guanine base is shaded purple. (** B**) Cartoon representation of base recognition and chemical mechanism of c-AMP-GMP 3′-5′-phosphodiester bond based on previous crystal structures ([@B23],[@B24],[@B26]).  Residue, D193, of the catalytic triad is highlighted in red.  Hydrogen bonds are indicated by dashed lines.  Interactions between the triphosphate backbone of the ligand and divalent ion are represented by dashed, orange lines.](gky137fig1){#F1}

Homeostasis of each of these molecules is maintained by distinct proteins in prokaryotic cells.  c-di-GMP, the most well-studied cyclic dinucleotide, is synthesized by diguanylate cyclases containing GGDEF domains ([@B16]) and it is hydrolyzed to either pGpG or GMP by phosphodiesterases containing HD-GYP or EAL domains ([@B6],[@B17]).  c-di-AMP is produced by diadenylate cyclases that specifically cyclize adenosine triphosphate ([@B18]).  These molecules are broken down by phosphodiesterases containing DHH/DHHA or HD domains ([@B19]).  Finally, c-AMP-GMP is synthesized by Hypr GGDEF proteins found in *Geobacter* ([@B20]) or DncV proteins ([@B14]) and broken down by phosphodiesterases containing HD-GYP-domains ([@B21]). 

While diguanylate cyclases typically dimerize to form the necessary active site for nucleotide triphosphate cyclization ([@B22]), DncV contains a single active site ([@B23]). * In vivo*, this unique active site is thought to preferentially recognize the adenine of adenosine triphosphate (ATP) in an anti conformation in the 'acceptor' pocket, while the guanine of guanosine triphosphate (GTP) is recognized in a syn conformation in a 'donor' binding pocket (Figure [1B](#F1){ref-type="fig"}) ([@B23]).  The resulting dinucleotide triphosphate is then released and rebound in the opposite orientation with guanosine bound at the acceptor pocket and adenosine bound at the donor pocket.  As a consequence of this mechanism, both GTP and ATP are recognized by either pocket.  For this reason, the enzyme can also cyclize two ATPs or two GTPs to form c-di-AMP or c-di-GMP, respectively ([@B14],[@B26]). 

Several groups have identified riboswitch and protein receptors that recognize cyclic dinucleotides in bacteria.  For example, there are two different riboswitch classes that recognize c-di-GMP by different mechanisms, the GEMM-I and GEMM-II riboswitches ([@B27]).  The GEMM-Ib riboswitches, which are structurally similar to the GEMM-I riboswitch, exclusively recognize c-AMP-GMP ([@B15],[@B32],[@B33]).  Finally, a structurally distinct riboswitch, the YdaO riboswitch, recognizes c-di-AMP in two different binding sites in a single aptamer ([@B34],[@B35]). 

In order to understand the molecular and phenotypic mechanisms of the cyclic dinucleotide pathway, several groups have synthesized analogs of these second messengers ([@B36]).  In addition to providing information about specific molecular interactions, some studies have found that certain analogs of cyclic dinucleotides bind differentially to RNA and protein receptors ([@B36],[@B38]).  It may be possible, therefore, to selectively target specific cyclic dinucleotide receptors by making such modifications.  Some analogs, such as those containing a phosphorothioate backbone, bind to phosphodiesterases but resist cleavage ([@B37]).  Therefore, cyclic dinucleotide analogs also show promise as tools to study cyclic dinucleotide pathways *in vivo*.  Development of a one-flask synthesis method for producing phosphorothioate modifications has been used to streamline synthesis of many cyclic dinucleotide analogs ([@B39],[@B40]).  However, this process cannot produce all phosphorothioate stereoisomers and 2′-hydroxyl modifications ([@B40]). 

Previous biochemical studies provided valuable insight into molecular interactions between the phosphate backbone and 2′-hydroxyl of c-di-GMP and its cognate riboswitches, GEMM-I and GEMM-II ([@B36],[@B43]).  However, c-di-GMP is a symmetrical molecule and can flip its orientation in the binding pocket.  Therefore, some of the loss in affinity due to modification in these areas could not be definitively attributed to specific hydrogen bond disruptions.  Here, we took advantage of the promiscuity of DncV, and further optimized its catalysis using different divalent ions, to create an expanded library of cyclic dinucleotide analogs modified at the sugar, base and phosphate backbone.  We scaled up production of several c-AMP-GMP and c-di-GMP analogs and used these compounds to biochemically probe phosphate backbone and sugar recognition of these ligands with c-AMP-GMP and c-di-GMP riboswitches. 

MATERIALS AND METHODS {#SEC2}
=====================

Purification of DncV {#SEC2-1}
--------------------

*Vibrio cholerae* DncV was expressed in BL21 cells and purified as previously described ([@B14],[@B24]).  In brief, cells expressing pET-DUET1-DncV were grown to an OD of 0.6 at 37°C and then induced with 1 mM isopropyl β-D-1-thiogalactopyranoside (IPTG) at 18°C for 4 h. The cells were harvested and resuspended in 50 mM KH~2~PO~4~, pH 7.0, 300 mM NaCl and 5 mM β-mercaptoethanol (BME).  Resuspended cells were lysed using a sonicating probe and clarified by centrifugation.  The supernatant was incubated with Nickel-agarose and the protein was eluted with 20, 50 and 250 mM imidazole.  The protein was dialyzed into 50 mM Tris--HCl, 100 mM NaCl, 10% glycerol and 5 mM β-mercaptoethanol and further purified by size exclusion chromatography (SEC) on a Superdex 75 column. 

Synthesis of cyclic dinucleotides {#SEC2-2}
---------------------------------

Nucleotide triphosphate analogs were purchased from Trilink Biotechnologies, Biolog, Sigma-Aldrich or New England Biolabs.  A total of 10 μM DncV was incubated with a total of 2 mM nucleotide triphosphate in 25 mM Tris--HCl, pH 7.5, 100 mM NaCl and 10 mM divalent ion.  For divalent metal screening, 10 mM MgCl~2~ was replaced by 10 mM CaCl~2~, CdCl~2~, CoCl~2~, CuCl~2~, MnCl~2~ or ZnCl~2~. Reactions were incubated as specified at 37°C and heated to 70°C for 10 min to quench the reaction.  Aliquots were separated by High Performance Liquid Chromatography (HPLC) using an Atlantis T3 3 μm, 4.6 × 150 mm column equilibrated in a mobile phase containing 10 mM ammonium acetate, pH 5.0 and acetonitrile as indicated in [Supplementary Table S1](#sup1){ref-type="supplementary-material"}.  The analogs were separated using a gradient of up to 50% acetonitrile with a flow rate of 1.0 ml/min.  Peak area percent was calculated using Agilent Chemstation based on absorbance signal at 254 nm.  A limit of detection was set based on a signal-to-noise ratio of 10 for c-AMP-GMP at 254 nm and limit of quantitation was set based on a signal-to-noise ratio of 3.  Major peaks were identified using electrospray ionization mass spectrometry in a mobile phase of formic acid and acetonitrile. 

Scale-up of cyclic dinucleotide analog synthesis {#SEC2-3}
------------------------------------------------

Analog synthesis used for binding studies were scaled up in 200--500 μl reactions.  Cyclic dinucleotides were purified on an Atlantis T3 3 μm, 9.6 × 300 mm column.  Purified peaks were identified by electrospray mass spectrometry and their purity confirmed to be \>95% using the analytical HPLC method described above.  c-AMP-GMP, c-di-GMP and c-di-AMP analog concentrations were calculated based on extinction coefficients of 25 050 M^−1^cm^−1^ at 256 nm, 23 700 M^−1^cm^−1^ at 254 nm and 27 000 M^−1^cm^−1^ at 259 nm, respectively. 

Synthesis of RNA {#SEC2-4}
----------------

The sequence for GEMM-Ib riboswitch aptamer, Gs 1761, from *Geobacter sulfurreducins* flanked by a T7 promoter at the 5′ end and the restriction enzyme recognition sequence for BsaI on the 3′ end was cloned into pUC19.  The RNA was synthesized *in vitro* by T7 polymerase from an EcoRV linearized plasmid and purified by denaturing gel electrophoresis followed by crush-and-soak.  The GEMM-I riboswitch, VC2 GUAA, and the GEMM-II riboswitch aptamers were synthesized as previously described ([@B30],[@B44]). 

Gel shift assays {#SEC2-5}
----------------

Radiolabeled c-AMP-GMP was prepared by incubating 17.5 μM GTP and α-^32^P-ATP with 10 μM DncV in 25 mM Tris--HCl, pH 7.5, 100 mM NaCl and 10 mM MnCl~2~ for 1 h. Radiolabeled c-di-GMP was prepared by incubating α-^32^P-GTP with the uninhibited mutant of a diguanylate cyclase from *Thermatoga maratime*, tDGC-R158A ([@B45]), for 1 h. Radiolabeled c-AMP-GMP was separated from c-di-AMP and c-di-GMP on a 20% native polyacrylamide gelelectrophoresis (PAGE) gel (100 mM Tris/HEPES, pH 7.5, 0.2 mM EDTA) at 4°C and extracted by soaking the gel slice in MilliQ water. 

To measure direct binding of c-AMP-GMP to the Gs 1761 aptamer, native gel shifts were used as previously described ([@B27]).  Briefly, increasing concentrations of RNA were incubated with trace radiolabeled c-AMP-GMP.  Bound and unbound ligand were separated on a 12% native PAGE gel (100 mM Tris/HEPES, pH 7.5, 10 mM MgCl~2~, 0.2 mM ethylenediaminetetraacetic acid (EDTA)) at 4°C.  Data were analyzed in ImageQuant and fit to a single site-specific binding model in GraphPad Prism. 

A competition gel shift assay previously used for c-di-GMP and c-di-AMP riboswitches ([@B36],[@B46]) was modified to measure the binding affinities of c-AMP-GMP and c-di-GMP analogs.  A total of 20 nM Gs 1761 was folded in the presence of trace radiolabeled c-AMP-GMP and increasing concentrations of cold c-AMP-GMP analog.  The binding reactions were incubated until they had reached equilibrium (at least 19 h) and the bound and unbound cyclic dinucleotide was separated on a native PAGE gel at 4°C.  Data were analyzed as described above and the *K*~d~ of each analog was determined using the following equation:$$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}\begin{equation*}\begin{array}{@{}*{1}{l}@{}} {{\rm{F}} = F_\infty ^C + \frac{{F_0^C}}{{2{{\rm{C}}_{cAG}}}}\left\{ {K_{\rm d}^{cAG} + \frac{{K_{\rm d}^{cAG}\ \times \ {C_{\rm T}}}}{{K_{\rm d}^C}} + {R_{\rm T}} + {C_{cAG\ }}} \right.}\\ {\left. { - {{ \left[ {{{\left( {K_{\rm d}^{cAG} + \ \frac{{K_{\rm d}^{cAG}\ \times \ {C_{\rm T}}}}{{K_{\rm d}^C}} + {R_{\rm T}}\ + \ {C_{cAG\ }}} \right)}^2} - \ \left( {4{R_{\rm T}}{C_{cAG\ }}} \right)\ } \right]}^{1/2}}} \right\}} \end{array}\end{equation*}\end{document}$$where $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}${\rm{F}}$\end{document}$ is the fraction bound, $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$F_\infty ^C$\end{document}$ is the fraction bound at saturating concentrations of cold competitor, $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$F_0^C$\end{document}$ is the fraction bound in the absence of cold competitor, $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$K_{\rm d}^{cAG}$\end{document}$ is the equilibrium dissociation constant of c-AMP-GMP bound to Gs 1761, $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}${{\rm{C}}_{cAG}}$\end{document}$ is the concentration of radiolabeled c-AMP-GMP (estimated to be 0.025 nM as previously calculated ([@B36])), $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}${{{C}}_{\rm T}}$\end{document}$ is the concentration of cold competitor, $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$K_{\rm d}^C$\end{document}$ is the equilibrium dissociation constant of the competitor and $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}${{{R}}_{\rm T}}$\end{document}$ is the total concentration of Gs1761 RNA (20 nM).  Competition assays for c-di-GMP analogs with GEMM-I and GEMM-II aptamers were performed with radiolabeled c-di-GMP as previously described ([@B36]). 

RESULTS AND DISCUSSION {#SEC3}
======================

DncV synthesizes several cyclic dinucleotides in the presence of magnesium {#SEC3-1}
--------------------------------------------------------------------------

DncV is a promiscuous dinucleotide cyclase from *Vibrio cholerae* that preferentially synthesizes c-AMP-GMP *in vivo* ([@B14]).  However, the enzyme can also make c-di-AMP and c-di-GMP *in vitro* when provided only ATP or GTP, respectively ([@B14],[@B26]).  In its preferred conformation, DncV is thought to recognize adenosine at the acceptor site by hydrogen bonding between S259 and the N6 of ATP as depicted in Figure [1](#F1){ref-type="fig"} ([@B23]).  The enzyme recognizes guanosine by hydrogen bonding to N1 and N2 of GTP via D348.  DncV then catalyzes a 3′,5′-phosphodiester bond between the two nucleotide triphosphates via the catalytic triad D131, D133 and D193 ([@B23],[@B24]).  The resulting dinucleotide triphosphate is released and rebound in the opposite orientation with guanosine recognized by S259 and Y179 in the acceptor pocket and adenosine binding to the donor pocket ([@B24]).  The second 3′,5′-linkage is then formed in order to complete the cyclic dinucleotide ([@B23],[@B24]). 

Because the acceptor and donor pockets are promiscuous for both ATP and GTP, we sought to determine if DncV could accept nucleotide triphosphates with base analogs.  We first tested substrates with conservative base modifications, such as inosine triphosphate and 6-thio-guanosine triphosphate (Figure [2A](#F2){ref-type="fig"}).  DncV synthesized cyclic dinucleotides with these modifications using magnesium as a divalent metal (Figure [2C](#F2){ref-type="fig"}).  Because 6-thio-guanosine is capable of chemical crosslinking with proteins, we tried this reaction in the presence and absence of dithiothreitol (DTT) to reverse this reactivity, but found that the cyclic dinucleotide was formed in both cases (data not shown). 

![ Base substitutions to cyclic dinucleotides and divalent ion optimization. (** A**) Base substitutions included on nucleotide triphosphates. (** B**) Divalent ion screen of DncV using GTP as substrate. (** C**) Cyclization of base analog nucleotide triphosphates using different divalent ions.](gky137fig2){#F2}

Under these conditions, DncV could not efficiently synthesize cyclic dinucleotides containing more dramatic modifications.  For example, we incubated DncV with purine triphosphate, which does not contain any exocyclic functional groups, or Z-nucleotide triphosphate (5-amino 4-imidazole carboxamide riboside 5′-triphosphate), the proposed bacterial alarmone and purine metabolite ([@B47]) which doesn't contain a full pyrimidine ring (Figure [2A](#F2){ref-type="fig"}).  In neither case did we observe a significant HPLC peak corresponding to a cyclized product.  Finally, to see if DncV could recognize other bases besides purines, we also tried incubating the pyrimidine nucleotide triphosphates, cytidine triphosphate (CTP) and uridine triphosphate (UTP), with DncV to form cyclic-di-CMP and cyclic-di-UMP, but found that they were also not cyclized efficiently (Figure [2C](#F2){ref-type="fig"}). 

Optimization of divalent metals in DncV increases cyclic dinucleotide production {#SEC3-2}
--------------------------------------------------------------------------------

Replacement of magnesium with other divalent ions is known to change the promiscuity and reaction rate for many nucleic acid enzymes ([@B50]).  This is especially well documented in restriction enzymes that are increasingly promiscuous in the presence of manganese rather than magnesium ([@B54]).  We therefore screened a number divalent metals to determine if any could alter the activity of DncV. We used GTP as a substrate and tested Ca^2+^, Cd^2+^, Co^2+^, Cu^2+^, Mn^2+^, Zn^2+^ for compatibility with DncV. As expected based upon other enzymatic systems, Ca^2+^ inhibited cyclic-di-GMP formation most likely due to its increased atomic radius (Figure [2B](#F2){ref-type="fig"}).  Zn^2+^, Cd^2+^ and Cu^2+^ also did not readily facilitate cyclic dinucleotide formation.  However, the enzyme could cyclize GTP to c-di-GMP in the presence of Mn^2+^ and Co^2+^ (Figure [2B](#F2){ref-type="fig"}). 

To determine if additional cyclic dinucleotides could be prepared in the presence of these alternative metals, we repeated the enzymatic synthesis of cyclic-di-PuMP (from purine triphosphate), cyclic-di-UMP (from uridine triphosphate), cyclic-di-CMP (from cytosine triphosphate) and cyclic-di-ZMP (from z-nucleotide triphosphate) in the presence of either cobalt or manganese.  All of these cyclic dinucleotide analogs, except cyclic-di-ZMP ([Supplementary Figure S1](#sup1){ref-type="supplementary-material"}), could be synthesized within a few hours in the presence of cobalt and the percent peak areas for all analogs were substantially increased when the enzymatic reaction was performed in the presence of manganese (Table [1](#tbl1){ref-type="table"}, Figure [2C](#F2){ref-type="fig"} and [Supplementary Figure S1](#sup1){ref-type="supplementary-material"}).  Because DncV successfully synthesized both purine and pyrimidine cyclic dinucleotide analogs while failing to prepare one with a Z-nucleotide triphosphate, it appears that there simply needs to be a complete pyrimidine ring present for nucleotide recognition in this enzyme. 

###### Synthesis of cyclic dinucleotide base analogs by DncV

  Compound     \% Peak area^a^   Incubation time (h)   Predicted major fragment m/z   Observed major fragment m/z
  ------------ ----------------- --------------------- ------------------------------ -----------------------------
  c-di-GMP     91                1                     689.1                          688.8
  c-di-AMP     98                1                     657.1                          656.9
  c-AMP-GMP    84                1                     673.1                          672.9
  c-di-CMP     63                1                     609.1                          609.0
  c-di-UMP     75                1                     611.3                          610.7
  c-di-IMP     100               1                     659.4                          658.9
  c-di-PuMP    68                2.5                   627.1                          626.8
  c-di-ZMP     \<LOD             N/A                   640.0                          N/A
  c-di-6TGMP   100               1                     722.1                          720.8

^a^Average of two reactions based on A~254~ signal in the presence of MnCl~2~.

LOD: limit of detection (0.1 μM/2 mM substrate). 

N/A: not applicable. 

We also retested the synthesis of cyclic-di-AMP, cyclic-di-GMP, cyclic-di-IMP and cyclic-di-6TGMP to see if the reaction efficiency could be increased by using a different metal.  All cyclic dinucleotide syntheses showed a slight increase in the product peak area when the divalent metal was changed from magnesium to cobalt or manganese, with the exception of cyclic-di-6TGMP (Figure [2C](#F2){ref-type="fig"}).  Interestingly, cobalt inhibited synthesis of cyclic-di-6TGMP while manganese promoted its synthesis. 

DncV synthesizes 2′-modified analogs of cyclic dinucleotides {#SEC3-3}
------------------------------------------------------------

No residues appear to directly interact with the 2′-hydroxyl of nucleotide triphosphates as observed within the crystal structures of DncV with its substrates ([@B23]) (Figure [1](#F1){ref-type="fig"}).  We tested if DncV could synthesize cyclic dinucleotides modified at this position in the presence of MnCl~2~ using nucleotide triphosphate analogs.  DncV reactivity was tested with ATP or GTP analogs that had either a fluorine (2′F), a hydrogen (deoxy, abbreviated with the suffix 'd'), a methoxy (2′OMe) or an amino group (2′NH~2~) (Table [2](#tbl2){ref-type="table"}).  c-AMP-GMP analogs were prepared by adding 1:1 ratios of the ATP and GTP analogs.  c-di-GMP and c-di-AMP analogs were prepared using only the ATP or GTP substrates, respectively. 

###### Synthesis of 2′-hydroxyl analogs of cyclic dinucleotides by DncV

  ---------------------
  ![]( gky137tbl2.jpg)
  ---------------------

Although the DncV does not appear to make direct 2′-OH contacts, the synthesis of 2′-substituted analogs was only possible in the presence of manganese.  In the presence of this divalent metal, the enzyme synthesized 2′-fluoro and 2′-amino analogs of all three cyclic dinucleotides (Table [2](#tbl2){ref-type="table"} and [Supplementary Figure S2](#sup1){ref-type="supplementary-material"}).  However, while nearly all of the 2′-fluoro nucleotide triphosphates were converted to cyclic dinucleotide analogs (c-di-2′F-GMP, c-di-2′F-AMP and c-2′F-AMP-2′F-GMP), much less of the 2′-amino analogs (c-di-2′NH~2~-GMP, c-di-2′NH~2~-AMP and c-2′NH~2~-AMP-2′NH~2~-GMP) were produced from the respective nucleotide triphosphate analogs (Table [2](#tbl2){ref-type="table"} and [Supplementary Figure S2](#sup1){ref-type="supplementary-material"}).  DncV also prepared 2′-amino monosubstitutions of cyclic-di-GMP and cyclic-AMP-GMP within an hour when incubated with 1:1 ratios of modified and unmodified nucleotide triphosphates.  Although the peak area for these monosubstituted products appears to be lower than many disubstitutions, the enzyme simultaneously synthesized other possible cyclic dinucleotides.  For example, DncV converted approximately 90% of substrates into products when incubated with equal amounts of unmodified GTP and modified 2′-NH~2~-GTP.  This included three products: the unsubstituted c-di-GMP cyclic dinucleotide, the disubstituted c-di-NH~2~-GMP analog and the monosubstituted c-2′NH~2~-G-GMP (Table [2](#tbl2){ref-type="table"} and [Supplementary Figure S2](#sup1){ref-type="supplementary-material"}).  The largest peak in the chromatogram from this reaction was c-di-GMP and the smallest peak was c-di-NH~2~-GMP.  This suggests that c-di-GMP is the most favorable product while the disubstituted product is the least favorable. 

There are limitations to the diversity of substitutions at the 2′-position that are accepted by DncV. The 2′-O-methyl was the bulkiest modification tested and no cyclic 2′-O-methylated dinucleotide products were observed even after incubation with the enzyme for \>20 h. In addition, DncV could not cyclize significant quantities of dGTP or dATP to form c-di-dGMP and c-di-dAMP, respectively ([Supplementary Figure S4](#sup1){ref-type="supplementary-material"} and Table [2](#tbl2){ref-type="table"}).  In both cases, large peaks formed earlier than would be expected for the 2′-deoxy cyclic dinucleotide analogs and none of them corresponded to m/z expected for these analogs ([Supplementary Figure S5](#sup1){ref-type="supplementary-material"}).  Both early-eluting peaks had major fragments with an m/z \>800 suggesting that they are likely to be dinucleotide triphosphate intermediates.  This implies the enzyme was able to complete the first, but not the second step of the dimerization reaction.  Interestingly, DncV was able to synthesize c-dAMP-dGMP when provided with equal amounts of dATP and dGTP, though the efficiency of this reaction was low.  The reaction had to be incubated overnight (\>20 h) and even then, only about a quarter of the substrate was converted to product (Table [2](#tbl2){ref-type="table"} and [Supplementary Figure S3](#sup1){ref-type="supplementary-material"}). 

DncV synthesizes phosphorothioate analogs of all three cyclic dinucleotides {#SEC3-4}
---------------------------------------------------------------------------

Phosphorothioate backbone-modified cyclic-di-AMP and cyclic-di-GMP are typically prepared using a one-flask synthesis method ([@B40]).  However, this method only allows synthesis of R~p~R~p~ and R~p~S~p~ phosphorothioate cyclic dinucleotides.  A more complicated synthetic procedure, requiring additional purification steps, is necessary to prepare the S~p~S~p~ phosphorothioate ([@B55]). 

We tested if DncV could convert α-phosphate-modified stereoisomers into the cyclic dinucleotide phosphorothioates.  We incubated different α-thio-GTP analogs with DncV and assigned products based upon the chemically-synthesized standards of c-(R~p~R~p~)-di-G~ps~ and c-(R~p~S~p~)-di-G~ps~ that had been previously reported ([@B36],[@B40]).  When a 1:1 mixture of S~p~-GTP-α-S and R~p~-GTP-α-S was incubated with the DncV enzyme in the presence of MnCl~2,~ a single peak matching the retention time and m/z value of c-(R~p~S~p~)-di-G~ps~ appeared after 1 h (Table [3](#tbl3){ref-type="table"} and [Supplementary Figure S4](#sup1){ref-type="supplementary-material"}).  Because we did not see peaks corresponding to other possible stereoisomers, c-(R~p~R~p~)-di-G~ps~ and c-(S~p~S~p~)-di-G~ps~, we attempted to prepare these individually.  α-S~p~-GTP-α-S was incubated with DncV in the presence of either MnCl~2~ or CoCl~2~. In the presence of CoCl~2~, a significant peak matching the retention time of the chemically synthesized c-(R~p~R~p~)-di-G~ps~ emerged ([Supplementary Figure S4](#sup1){ref-type="supplementary-material"}).  This sample was evaluated by mass spectrometry and had an m/z of 721, consistent with a diphosphorothioate analog of c-di-GMP.  This inversion of stereochemistry is consistent with previous enzymatic synthesis reactions ([@B56],[@B57]).  Finally, R~p~-GTP-α-S was incubated with the enzyme and either CoCl~2~ or MnCl~2~ to prepare c-(S~p~S~p~)-di-G~ps~. After 20 h, a small peak emerged at ∼6.5 min in the presence of MnCl~2~. We identified this peak as c-(S~p~S~p~)-di-G~ps~ because the molecule had an m/z value of 721 like the other two c-di-GMP diphosphorothioates and because the retention time was in between that of c-di-GMP and c-(R~p~S~p~)-di-G~ps~. This same order of elution was previously reported for the three c-di-GMP diphosphorothioate stereoisomers ([@B55]). 

###### Synthesis of phosphorothioate analogs of cyclic dinucleotides by DncV

  ---------------------
  ![]( gky137tbl3.jpg)
  ---------------------

Taken together, it is not unexpected that the c-(R~p~S~p~)-G~ps~ analog is the sole product in a mixture of stereoisomeric starting materials.  c-(R~p~R~p~)-G~ps~ is not readily formed in the presence of MnCl~2~ ([Supplementary Figure S4](#sup1){ref-type="supplementary-material"}).  While it was possible to form the c-(S~p~S~p~)-di-G~ps~ diphosphorothioate in the presence of this divalent ion, the reaction was not efficient and required overnight incubation.  These observations further support the sole production of the c-(R~p~S~p~)-G~ps~ analog, rather than either of the two stereoisomers, when a mixture of S~p~-GTP-α-S and R~p~-GTP-α-S are incubated with the enzyme in the presence of MnCl~2~.

We performed a similar experiment by incubating α-thio-ATP analogs incubated with DncV in the presence of either MnCl~2~ or CoCl~2~ and observed the same patterns for stereochemical inversion as seen for c-di-GMP.  We used the previously characterized diphosphorothioate analogs ([@B46]) c-(R~p~R~p~)-di-A~ps~ and c-(R~p~S~p~)-di-A~ps~ as standards.  When S~p~-ATP-α-S was incubated with DncV in the presence of CoCl~2~, a peak matching the retention time and m/z of c-(R~p~R~p~)-di-A~ps~ was produced.  When a mixture of S~p~-ATP-α-S and R~p~-ATP-α-S was incubated in the presence of MnCl~2~, a single peak matching that of c-(R~p~S~p~)-di-A~ps~ was visible.  This is consistent with the efficient formation of the homologous c-(R~p~S~p~)-G~ps~ analog described above.  When R~p~-ATP-α-S was incubated with DncV in the presence of MnCl~2~, a peak with an m/z value matching a diphosphorothioate analog of c-di-AMP emerged.  This peak had a retention time between c-di-AMP and c-(R~p~S~p~)-di-A~ps~ which is consistent with the elution order of a c-(S~p~S~p~)-di-A~ps~ analog as seen for the c-di-GMP diphosphorothioates ([@B55]). 

As with c-di-GMP phosphorothioates, the divalent metal identity seems to be especially important for production of R~p~R~p~ and S~p~S~p~ phosphorothioates of c-di-GMP and c-di-AMP.  The crystal structures of DncV show divalent ion coordination with the α-pro-R~p~ oxygen of ATP and GTP ([@B23],[@B24]) in the donor pocket ([Supplementary Figure S5](#sup1){ref-type="supplementary-material"}).  Replacement of these oxygens with sulfurs most likely disrupts this coordination.  Both Co^2+^ and Mn^2+^ are thiophilic and should accommodate the phosphorothioate.  However, CoCl~2~ is the preferred divalent ion for production of R~p~R~p~ diphosphorothioates while MnCl~2~ is preferred for S~p~S~p~ diphosphorothioates (Figure [3](#F3){ref-type="fig"} and [Supplementary Figure S4](#sup1){ref-type="supplementary-material"}).  Further structural experiments would be necessary to fully understand the preference of MnCl~2~ or CoCl~2~ in these reactions. 

![ Enzymatic synthesis of phosphorothioate c-di-AMP analogs. (** A**) Chemical structure of c-di-AMP phosphorothioate analogs.  Conformation of the thiol groups is shown in red. (** B**) Chromatogram of chemically synthesized c-di-AMP phosphorothioates. (** C**) Chromatogram of enzymatically synthesized c-(R~p~S~p~)-A~ps~ in the presence of MnCl~2~ from α-R~p~-ATP and α-S~p~-ATP. (** D**) Chromatograms of enzymatically synthesized c-(R~p~R~p~)-A~ps~ from α-S~p~-ATP in the presence of MnCl~2~ (blue trace) and CoCl~2~ (red trace). (** E**) Chromatograms of enzymatically synthesized c-(S~p~S~p~)-A~ps~ from α-R~p~-ATP in the presence of MnCl~2~ (blue trace) and CoCl~2~ (red trace).](gky137fig3){#F3}

There are four possible stereoisomers of c-AMP-GMP phosphorothioates ([Supplementary Figure S6](#sup1){ref-type="supplementary-material"}).  DncV synthesizes all of these analogs when given appropriate substrates (Table [3](#tbl3){ref-type="table"}).  Contrary to the symmetrical cyclic dinucleotides, synthesis of the c-(R~p~R~p~)-A~ps~G~ps~ analog did not require CoCl~2~ as the metal cofactor.  We enzymatically synthesized c-(S~p~R~p~)-A~ps~G~ps~ and c-(R~p~S~p~)-A~ps~G~ps~ and determined that they were different compounds due to the difference in starting material and product peak retention times.  We used R~p~-ATP-α-S + S~p~-GTP-α-S to produce c-(S~p~R~p~)-A~ps~G~ps~ and S~p~-ATP-α-S + R~p~-GTP-α-S to produce c-(R~p~S~p~)-A~ps~G~ps~. These reactions resulted in product peaks eluting at retention times of 6.8 and 7.0 min for c-(S~p~R~p~)-A~ps~G~ps~ and c-(R~p~S~p~)-A~ps~G~ps~, respectively.  Both analogs had the same observed major fragment by mass spectrometry and had relatively good yield of greater than 70% peak area in 1 h as seen in Table [3](#tbl3){ref-type="table"}.  The synthesis of the c-(S~p~S~p~)-A~ps~G~ps~ analog was the least efficient of the four possible analogs, taking up to 2.5 h to convert about half of the starting material into product.  When we scaled up these reactions starting from 400 μg substrate, we were able to obtain 13, 12, 19 and 5% yield following HPLC purification for c-(R~p~R~p~)-A~ps~G~ps~, c-(R~p~S~p~)-A~ps~G~ps~, c-(S~p~R~p~)-A~ps~G~ps~, c-(S~p~S~p~)-A~ps~G~ps~, respectively.  These values are within range of chemical synthesis of c-di-GMP phosphorothioates ([@B40],[@B55]). 

2′-amino analogs of c-AMP-GMP provide insight into cyclic dinucleotide recognition by riboswitches {#SEC3-5}
--------------------------------------------------------------------------------------------------

Previous studies have reported differential binding of 2′-modified cyclic-di-GMP analogs between the GEMM-I and GEMM-II riboswitches ([@B36]) and even among GEMM-I riboswitches from different bacterial species ([@B43]).  The crystal structure of the GEMM-I riboswitch, VC2, shows a non-canonical Watson--Crick/Hoogsteen interaction between a guanine of c-di-GMP, known as the Gα, and G20 of the riboswitch ([@B27]).  The other guanine of c-di-GMP, known as Gβ, binds to C92 of the riboswitch through a Watson--Crick base pair.  This orientation allows for other points of contact between the riboswitch bases and the 2′-hydroxyl of the ligand ([@B27]).  In particular, the 2′-hydroxyl of Gα hydrogen bonds to the phosphate backbone of A47 while the 2′-hydroxyl of Gβ is coordinated to a water molecule (Figure [4A](#F4){ref-type="fig"}).  The GEMM-II riboswitch forms two non-canonical base pairs between both the Gα and Gβ of c-di-GMP, but makes fewer contacts with other functional groups of the ligand ([@B30]).  For this reason, many c-di-GMP analogs show preferential binding to the GEMM-II rather than the GEMM-I riboswitch.  However, the effect of an electropositive group, such as an amino, at the 2′ position of cyclic dinucleotides has not previously been tested.  We synthesized 2′-amino analogs of c-di-GMP and c-AMP-GMP using DncV and tested the effects of this modification for binding to the two classes of c-di-GMP riboswitches (GEMM-I and GEMM-II) and a c-AMP-GMP (GEMM-Ib) riboswitch. 

![ 2′-hydroxyl and phosphate backbone recognition of c-di-GMP and c-AMP-GMP in GEMM-I and GEMM-Ib riboswitches. (** A**) Hydrogen bonding to 2′-hydroxyl and phosphate backbone of c-di-GMP in the GEMM-I riboswitch, VC2, from pdb file 3muh ([@B44]).  The riboswitch is pictured in gray and cyclic-di-GMP is shown in purple. (** B**) Hydrogen bonding to 2′-hydroxyl and phosphate backbone in GEMM-Ib riboswitch, Gs1761, from pdb file 4yaz ([@B33]).  The riboswitch is pictured in gray and cyclic-AMP-GMP is shown in cyan.](gky137fig4){#F4}

We measured the binding of the c-di-2′NH~2~-GMP, which contains two 2′-amino modifications, to the GEMM-I and GEMM-II riboswitch using a competition gel shift assay.  As expected, there was a negligible 1.6-fold decrease in affinity to the GEMM-II riboswitch.  However, this same modification resulted in a significant 1700-fold loss in affinity to the GEMM-I riboswitch (Table [4](#tbl4){ref-type="table"}).  This large effect cannot be explained by disruption of hydrogen bonding or electrostatics.  In its deprotonated state, the 2′-amino should maintain the hydrogen bond donor properties of the 2′-hydroxyl.  If the 2′-amino is protonated, it should stabilize electrostatic interactions with the negatively charged phosphate backbone of A47 in the GEMM-I riboswitch rather than disrupt it.  However, a large binding defect (1800-fold) was also reported for binding of the 2′-deoxy analog, c-di-dGMP, to the GEMM-I riboswitch ([@B36]).  The reported result for the 2′-deoxy analog was attributed to a change in sugar pucker from the 3′-endo conformation seen in c-di-GMP to a 2′-endo sugar pucker upon replacement of the 2′-hydroxyl with a hydrogen in c-di-dGMP ([@B36]).  A 2′-amino group is also thought to alter the ribose sugar pucker of nucleotides from the 3′-endo conformation to the 2′-endo ([@B58]).  Therefore, this conformational change upon introduction of a 2′-amino group is likely responsible for perturbing the binding affinity of the ligand. 

###### Binding affinity of 2′ modified cyclic dinucleotides to riboswitches

  Compound                    Riboswitch          *K* ~d~ (nM)      Fold change   ΔΔG (kcal/mol)
  --------------------------- ------------------- ----------------- ------------- ----------------
  c-di-GMP                    GEMM-I (GUAA VC2)   2.2 ± 0.8^a^      N/A           N/A
  c-2′NH~2~-G-GMP             GEMM-I (GUAA VC2)   68 ± 10           30            2.0
  c-di-2′NH~2~-GMP            GEMM-I (GUAA VC2)   3800 ± 800        1700          4.4
  c-di-GMP                    GEMM-II             2.2 ± 0.2^a,b^    N/A           N/A
  c-di-2′NH~2~-GMP            GEMM-II             3.5 ± 0.4         1.6           0.3
  c-AMP-GMP                   GEMM-Ib (Gs 1761)   1.1 ± 0.04^a,c^   N/A           N/A
  c-2′NH~2~-A-GMP             GEMM-Ib (Gs 1761)   21 ± 2            19            1.3
  c-AMP-2′NH~2~-G             GEMM-Ib (Gs 1761)   9.0 ± 1.8         8.2           0.8
  c-2′NH~2~-AMP-2′NH~2~-GMP   GEMM-Ib (Gs 1761)   2500 ± 230        2300          4.1

^a^Denotes measurements made by direct binding. 

^b^Previously reported in reference ([@B44]). 

^c^Previously reported in reference ([@B30]). 

N/A: not applicable. 

In order to understand the contributions of a single 2′-amino substitution, we also tested binding of c-2′NH~2~-G-GMP to the GEMM-I riboswitch.  This analog was prepared by scaling up the reaction of a 1:1 ratio of GTP and 2′NH~2~-GTP with DncV. Although only the 2′-hydroxyl group of Gα makes direct contact with the riboswitch, this monosubstituted analog still resulted in a 30-fold decrease in binding affinity, corresponding to a 2.0 kcal/mol loss and roughly half of the cost of the disubstituted analog, c-di-2′NH~2~-GMP.  The significant binding defect suggests that a single 2′-amino substitution is enough to alter the sugar pucker and likely disrupts the constraints of the entire ligand in the binding pocket of the riboswitch. 

Because c-di-GMP is a symmetric molecule and can flip its orientation in the binding pocket, the c-2′-NH~2~-G-GMP analog is expected to bind in the orientation that has the lowest energetic effect on the binding affinity.  To eliminate the ambiguity of this possibility, we also tested analog binding to a GEMM-1b riboswitch, Gs 1761 ([Supplementary Figure S7](#sup1){ref-type="supplementary-material"}).  This is a cyclic dinucleotide riboswitch variant that folds in a similar structure to the GEMM-I riboswitch, but specifically recognizes c-AMP-GMP.  This riboswitch has also been crystallized ([@B33]) and the Aα of c-AMP-GMP binds to A14 and Gβ forms a Watson--Crick pair with C75 ([Supplementary Figure S7](#sup1){ref-type="supplementary-material"}).  However, the orientation of c-AMP-GMP allows both 2′-hydroxyls to be in direct contact with the riboswitch (Figure [4B](#F4){ref-type="fig"}).  The 2′-hydroxyl of the Aα hydrogen bonds to the phosphate backbone of A41 similar to interactions in the GEMM-I riboswitch.  In addition, the Gβ 2′-hydroxyl of c-AMP-GMP forms a hydrogen bond with N7 of A11 in the riboswitch.  Because c-AMP-GMP is an asymmetric molecule, it is possible to make specific modifications to the 2′-hydroxyl groups of the α and β sites to individually probe their interactions with the riboswitch. 

We measured the binding affinity of native and modified ligands to the Gs 1761 GEMM-Ib riboswitch using direct binding gel shifts and competition gel shifts, respectively.  In the direct binding assay, c-AMP-GMP bound to the riboswitch with an affinity of 1.1 nM. This is in agreement with the previously-reported, 530 pM, as measured by inline probing ([@B32]).  When we measured binding affinity of c-2′NH~2~-AMP-2′NH~2~-GMP, which has 2′-amino substitutions at both hydroxyls, we found a 2300-fold weaker binding affinity, corresponding to 4.1 kcal/mol (Table [4](#tbl4){ref-type="table"}).  This value is similar to what was observed for c-di-2′NH~2~-GMP binding to the GEMM-I riboswitch (ΔΔG=4.4 kcal/mol) suggesting that the effect of disubstitution is the same in both riboswitches. 

To determine the contribution of 2′-hydroxyl recognition at the two different binding sites, we next tested binding of monosubstituted 2′-amino-modified cyclic-AMP-GMP analogs.  The c-2′NH~2~-A-GMP, which contains a modification of the 2′-hydroxyl at the Aα site, showed a 19-fold change in binding affinity (1.3 kcal/mol), while c-AMP-2′NH~2~-G, which is modified at the 2′-hydroxyl at the Gβ site, showed only an ∼8-fold loss in binding affinity (0.8 kcal/mol).  Even though both 2′-hydroxyls make contact with the riboswitch (Figure [4B](#F4){ref-type="fig"}), these results indicate that the 2′-hydroxyl at the Aα site is slightly more important for recognition than the 2′-hydroxyl at the Gβ site. 

Taken together, these data suggest a difference in ligand recognition of the 2′-hydroxyl not evident from the static crystal structures of the GEMM-I and GEMM-Ib riboswitch.  Because the Gβ 2′-hydroxyl of c-di-GMP does not make direct contact with the riboswitch (Figure [4A](#F4){ref-type="fig"}), it is expected that the single 2′-amino group in the monosubstituted analog is oriented at this site.  However, it still resulted in a 2 kcal/mol energetic loss in affinity, which was attributed to a change in the sugar pucker.  In contrast, modification to the analogous location in c-AMP-GMP resulted in only a 0.8 kcal/mol change.  This suggests that the GEMM-Ib riboswitch is more tolerant of conformational changes at this site. 

Phosphorothioate-modified c-AMP-GMP analogs bind to a c-AMP-GMP riboswitch {#SEC3-6}
--------------------------------------------------------------------------

We next sought to understand the interactions between the phosphate backbone of c-AMP-GMP and the GEMM-Ib riboswitch, Gs1761, using the expanded set of phosphorothioate substituted analogs.  The crystal structure of the GEMM-I riboswitch shows hydrated divalent metals and other water molecules coordinating to at least one of the pro-R~p~ non-bridging oxygens.  However, the crystal structure of the GEMM-Ib bound to c-AMP-GMP does not contain an equivalent set of interactions in the binding site (Figure [4](#F4){ref-type="fig"}). 

Previous studies on the GEMM-I riboswitch reported decreased binding affinities for all phosphorothioate analogs compared to c-di-GMP and attributed this result to disruption of electrostatic interactions in the binding pocket ([@B36]).  However, it is unclear how phosphorothioate modifications are oriented in the binding pocket of the GEMM-I riboswitch due to the symmetry of c-di-GMP.  It is possible to map S~p~ and R~p~ phosphorothioate modifications on c-AMP-GMP within the GEMM-Ib binding pocket.  Given that there unexpected differences between GEMM-I and GEMM-Ib riboswitches when the 2′-hydroxyl at the β site was modified to an amino group, we considered if modification of the phosphate backbone would also reveal differences in conformation or hydrogen bonding between the two riboswitches. 

To elucidate the interactions between the riboswitch and the phosphate backbone, we tested binding of c-AMP-GMP analogs with monosubstitutions of oxygen to sulfur at each non-bridging oxygen of c-AMP-GMP to the GEMM-Ib riboswitch.  Reactions of modified and unmodified nucleotide triphosphates were scaled up to prepare each of the four different possibilities: c-(R~p~)-A~ps~-GMP, c-AMP-(R~p~)-G~ps~, c-(S~p~)-A~ps~-GMP and c-AMP-(S~p~)-G~ps~. c-(R~p~)-A~ps~-GMP and c-(S~p~)-A~ps~-GMP are stereoisomers of sulfur substitution for the non-bridging oxygens at the Aα site.  c-AMP-(R~p~)-G~ps~ and c-AMP-(S~p~)-G~ps~ are modified at the non-bridging oxygens at the Gβ site. 

Neither of the analogs with modifications to the pro-R~p~ oxygens, c-(R~p~)-A~ps~-GMP and c-AMP-(R~p~)-G~ps~, demonstrated large changes in binding affinity (only 1.4- to 1.9-fold changes in *K*~d~).  These values are similar to the *K*~d~ reported for the analogous c-(R~p~)-G~ps~-GMP bound to the GEMM-I construct, which shows only a 4-fold change in *K*~d~ when compared to c-di-GMP ([@B36]).  These data are supported by the crystal structures of GEMM-I and GEMM-Ib that show pro-R~p~ oxygens of the c-di-GMP and c-AMP-GMP ligands pointed toward the solvent rather than making direct contact to their respective riboswitches (Figure [4](#F4){ref-type="fig"}).  We also measured the binding affinity of monophosphorothioate modifications to the pro-S~p~ oxygens.  Both of these analogs showed significant losses in binding affinity when compared to c-AMP-GMP.  The c-AMP-(S~p~)-G~ps~ analog showed a 10-fold loss in *K*~d~ compared to c-AMP-GMP.  c-(S~p~)-A~ps~-GMP showed the largest defect in binding affinity (a 43-fold loss) confirming that this oxygen is especially important for recognition.  This is in agreement with the crystal structure of the GEMM-Ib riboswitch ([@B33]) that shows the pro-S~p~ oxygen of Aα making hydrogen bonds with A41 and A12 (Figure [4B](#F4){ref-type="fig"}). 

We next tested if phosphorothioate modifications were additive or if there might be other factors implicated in recognition of disubstituted analogs as seen for the 2′-amino modifications.  All diphosphorothioate modifications had a significant effect on riboswitch binding.  The c-(R~p~R~p~)-A~ps~G~ps~ analog bound the tightest out of the diphosphorothioates (Table [5](#tbl5){ref-type="table"}), but still had a 29-fold effect on binding affinity.  The c-(R~p~S~p~)-A~ps~G~ps~ and c-(S~p~R~p~)-A~ps~G~ps~ analogs showed even weaker binding when compared to c-AMP-GMP, 140- and 330-fold effects, respectively.  Both of these analogs contain at least one pro-S~p~ oxygen modified to a sulfur and oriented toward the binding pocket.  The pro-S~p~ oxygen at the Aα of c-AMP-GMP shows the most dramatic loss in binding affinity when changed to a sulfur as evidenced by the c-(S~p~)-A~ps~-GMP analog.  This modification is also present in the c-(S~p~R~p~)-A~ps~G~ps~ diphosphorothioate and further supports the importance of this oxygen for ligand recognition. 

###### Binding affinity of c-AMP-GMP and phosphorothioate analogs to Gs1761

  Compound                  *K* ~d~ (nM)                          Fold change   ΔΔG (kcal/mol)
  ------------------------- ------------------------------------- ------------- ----------------
  c-AMP-GMP                 1.1$\documentclass[12pt]{minimal}     N/A           N/A
                            \usepackage{amsmath}                                
                            \usepackage{wasysym}                                
                            \usepackage{amsfonts}                               
                            \usepackage{amssymb}                                
                            \usepackage{amsbsy}                                 
                            \usepackage{upgreek}                                
                            \usepackage{mathrsfs}                               
                            \setlength{\oddsidemargin}{-69pt}                   
                            \begin{document}                                    
                            }{}$ \pm$\end{document}$ 0.04^a^                    
  c-(R~p~)-A~ps~-GMP        2.1 $\documentclass[12pt]{minimal}    1.9           0.4
                            \usepackage{amsmath}                                
                            \usepackage{wasysym}                                
                            \usepackage{amsfonts}                               
                            \usepackage{amssymb}                                
                            \usepackage{amsbsy}                                 
                            \usepackage{upgreek}                                
                            \usepackage{mathrsfs}                               
                            \setlength{\oddsidemargin}{-69pt}                   
                            \begin{document}                                    
                            }{}$ \pm$\end{document}$ 0.1                        
  c-AMP-(R~p~)-G~ps~        1.5 $\documentclass[12pt]{minimal}    1.4           0.2
                            \usepackage{amsmath}                                
                            \usepackage{wasysym}                                
                            \usepackage{amsfonts}                               
                            \usepackage{amssymb}                                
                            \usepackage{amsbsy}                                 
                            \usepackage{upgreek}                                
                            \usepackage{mathrsfs}                               
                            \setlength{\oddsidemargin}{-69pt}                   
                            \begin{document}                                    
                            }{}$ \pm$\end{document}$ 0.2                        
  c-(S~p~)-A~ps~-GMP        48 $\documentclass[12pt]{minimal}     43            2.2
                            \usepackage{amsmath}                                
                            \usepackage{wasysym}                                
                            \usepackage{amsfonts}                               
                            \usepackage{amssymb}                                
                            \usepackage{amsbsy}                                 
                            \usepackage{upgreek}                                
                            \usepackage{mathrsfs}                               
                            \setlength{\oddsidemargin}{-69pt}                   
                            \begin{document}                                    
                            }{}$ \pm$\end{document}$ 5                          
  c-AMP-(S~p~)-G~ps~        11 $\documentclass[12pt]{minimal}     10            1.4
                            \usepackage{amsmath}                                
                            \usepackage{wasysym}                                
                            \usepackage{amsfonts}                               
                            \usepackage{amssymb}                                
                            \usepackage{amsbsy}                                 
                            \usepackage{upgreek}                                
                            \usepackage{mathrsfs}                               
                            \setlength{\oddsidemargin}{-69pt}                   
                            \begin{document}                                    
                            }{}$ \pm$\end{document}$ 2                          
  c-(R~p~R~p~)-A~ps~G~ps~   32 $\documentclass[12pt]{minimal}     29            2.0
                            \usepackage{amsmath}                                
                            \usepackage{wasysym}                                
                            \usepackage{amsfonts}                               
                            \usepackage{amssymb}                                
                            \usepackage{amsbsy}                                 
                            \usepackage{upgreek}                                
                            \usepackage{mathrsfs}                               
                            \setlength{\oddsidemargin}{-69pt}                   
                            \begin{document}                                    
                            }{}$ \pm$\end{document}$ 3                          
  c-(R~p~S~p~)-A~ps~G~ps~   150 $\documentclass[12pt]{minimal}    140           2.9
                            \usepackage{amsmath}                                
                            \usepackage{wasysym}                                
                            \usepackage{amsfonts}                               
                            \usepackage{amssymb}                                
                            \usepackage{amsbsy}                                 
                            \usepackage{upgreek}                                
                            \usepackage{mathrsfs}                               
                            \setlength{\oddsidemargin}{-69pt}                   
                            \begin{document}                                    
                            }{}$ \pm$\end{document}$ 20                         
  c-(S~p~R~p~)-A~ps~G~ps~   360 $\documentclass[12pt]{minimal}    330           3.4
                            \usepackage{amsmath}                                
                            \usepackage{wasysym}                                
                            \usepackage{amsfonts}                               
                            \usepackage{amssymb}                                
                            \usepackage{amsbsy}                                 
                            \usepackage{upgreek}                                
                            \usepackage{mathrsfs}                               
                            \setlength{\oddsidemargin}{-69pt}                   
                            \begin{document}                                    
                            }{}$ \pm$\end{document}$ 60                         
  c-(S~p~S~p~)-A~ps~G~ps~   1200 $\documentclass[12pt]{minimal}   1100          4.1
                            \usepackage{amsmath}                                
                            \usepackage{wasysym}                                
                            \usepackage{amsfonts}                               
                            \usepackage{amssymb}                                
                            \usepackage{amsbsy}                                 
                            \usepackage{upgreek}                                
                            \usepackage{mathrsfs}                               
                            \setlength{\oddsidemargin}{-69pt}                   
                            \begin{document}                                    
                            }{}$ \pm$\end{document}$ 130                        

^a^Denotes measurements made by direct binding. 

N/A: not applicable. 

Because both S~p~ modifications showed weaker binding than the R~p~ phosphorothioate analogs, it is not surprising that the c-(S~p~S~p~)-A~ps~G~ps~ analog bound with the weakest affinity with a 1100-fold effect on binding (Table [5](#tbl5){ref-type="table"}).  Interestingly, c-(S~p~S~p~)-A~ps~G~ps~ shows a nearly additive effect when comparing affinity loss to individual monophosphorothioates.  The combined energy lost with each individual monophosphorothioate, c-(S~p~)-A~ps~-GMP, and c-AMP-(S~p~)-G~ps,~ is 3.6 kcal/mol and the ΔΔG of c-(S~p~S~p~)-A~ps~G~ps~ is 4.1 kcal/mol.  However, all other diphosphorothioates present a discrepancy of at least 1 kcal/mol in comparison to the individual monophosphorothioates.  The largest difference is seen in c-(R~p~R~p~)-A~ps~G~ps~ that has an extra 1.4 kcal/mol loss in binding affinity that cannot be explained by simply adding the ΔΔG values of the individual pro-R~p~ monosubstitutions.  Therefore, it is likely that other interactions and/or alternative ligand conformations beyond the resolution of our experiments are involved in recognition of the pro-R~p~ oxygens in the disubstituted analog. 

CONCLUSION {#SEC4}
==========

The binding partner for the GEMM-I riboswitch was identified as c-di-GMP almost a decade ago ([@B29]).  More recently, a subset of these RNA receptors were shown to selectively bound c-AMP-GMP ([@B15],[@B32]).  While these two riboswitches are globally similar, there are nuances to their molecular recognition that were previously ambiguous.  The enzymatic synthesis method presented in this study allowed us to synthesize novel cyclic dinucleotide analogs to reveal differences in the sugar recognition in these two riboswitches, specifically at the Gβ site, that is not evident from the crystal structures.  The data also suggest a conformational change in the diphosphorothioate cyclic dinucleotide analogs possibly through constraints on the cyclized dinucleotide backbone.  While specific modifications to the phosphate and base of these second messenger molecules were used to probe targeted receptor recognition, still other modifications could be made using enzymatic synthesis in order to identify protein and RNA binding partners in bacterial cyclic dinucleotide pathways. 

Supplementary Material
======================

###### 

Click here for additional data file. 

The authors would like to thank Roger Jones for generously providing chemically synthesized c-di-AMP and c-di-GMP phosphorothioate analogs.  The plasmid containing DncV was a generous gift of John Mekalanos.  The authors are most grateful to David Hiller and Caroline Reiss for providing helpful comments on this manuscript and other members of the Strobel Laboratory for invaluable discussion.  This research made use of the West Campus Analytical Core at Yale University. 

SUPPLEMENTARY DATA {#SEC5}
==================

[Supplementary Data](https://academic.oup.com/nar/article-lookup/doi/10.1093/nar/gky137#supplementary-data) are available at NAR Online. 

FUNDING {#SEC6}
=======

National Institutes of Health \[P01 GM022778 to S.A.S., F32 GM115097 to K.D.L-F\].  Funding for open access charge: National Institutes of Health \[F32 GM115097\]. 

*Conflict of interest statement*.  None declared.