# Structure of a Diguanylate Cyclase from Thermotoga maritima: Insights into Activation, Feedback Inhibition and Thermostability

Angeline Deepthi, Chong Wai Liew, Zhao-Xun Liang, Kunchithapadam Swaminathan, Julien Lescar

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Angeline Deepthi, Chong Wai Liew, Zhao-Xun Liang, Kunchithapadam Swaminathan, Julien Lescar. Structure of a Diguanylate Cyclase from Thermotoga maritima: Insights into Activation, Feedback Inhibition and Thermostability. PLoS ONE, 2014, 9 (10), pp.e110912. 10.1371/journal.pone.0110912. hal-01345084

HAL Id: hal-01345084

https://hal.sorbonne-universite.fr/hal-01345084v1

Submitted on 13 Jul 2016

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# Structure of a Diguanylate Cyclase from *Thermotoga maritima*: Insights into Activation, Feedback Inhibition and Thermostability

Angeline Deepthi¹, Chong Wai Liew², Zhao-Xun Liang², Kunchithapadam Swaminathan¹*,
Julien Lescar²,³*

¹Department of Biological Sciences, National University of Singapore, Singapore, Singapore, ²School of Biological Sciences, Nanyang Technological University, Singapore, Singapore, ³Centre d'Immunologie et des Maladies Infectieuses, Centre Hospitalier Universitaire Pitié-Salpêtrière Faculté de Médecine Pierre et Marie Curie, Paris, France

## Abstract

Large-scale production of bis-3'-5'-cyclic-di-GMP (c-di-GMP) would facilitate biological studies of numerous bacterial signaling pathways and phenotypes controlled by this second messenger molecule, such as virulence and biofilm formation. C-di-GMP constitutes also a potentially interesting molecule as a vaccine adjuvant. Even though chemical synthesis of c-di-GMP can be done, the yields are incompatible with mass-production. tDGC, a stand-alone diguanylate cyclase (DGC or GGDEF domain) from *Thermotoga maritima*, enables the robust enzymatic production of large quantities of c-di-GMP. To understand the structural correlates of tDGC thermostability, its catalytic mechanism and feedback inhibition, we determined structures of an active-like dimeric conformation with both active (A) sites facing each other and of an inactive dimeric conformation, locked by c-di-GMP bound at the inhibitory (I) site. We also report the structure of a single mutant of tDGC, with the R158A mutation at the I-site, abolishing product inhibition and unproductive dimerization. A comparison with structurally characterized DGC homologues from mesophiles reveals the presence of a higher number of salt bridges in the hyperthermophile enzyme tDGC. Denaturation experiments of mutants disrupting in turn each of the salt bridges unique to tDGC identified three salt-bridges critical to confer thermostability.

**Citation:** Deepthi A, Liew CW, Liang Z-X, Swaminathan K, Lescar J (2014) Structure of a Diguanylate Cyclase from *Thermotoga maritima*: Insights into Activation, Feedback Inhibition and Thermostability. PLoS ONE 9(10): e110912. doi:10.1371/journal.pone.0110912

**Editor:** Emily Parker, University of Cantebury, New Zealand

**Received** July 19, 2014; **Accepted** September 24, 2014; **Published** October 31, 2014

**Copyright:** © 2014 Deepthi et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability:** The authors confirm that all data underlying the findings are fully available without restriction. All PDB files are available from the RCSB database (accession numbers 4URG, 4URQ, 4URS).

**Funding:** Supported by Competitive Research Programme Grant CRP2008 from the National Research Foundation (to the J. L. laboratory). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

* Email: dbsks@nus.edu.sg (KS); julien@ntu.edu.sg (JL)

## Introduction

Bis-(3'-5') cyclic diguanosine monophosphate (c-di-GMP) is a second messenger that was initially discovered as an allosteric activator of cellulose synthase in *Gluconobacter xylinus* [1]. Binding of c-di-GMP to various effector or receptor proteins regulates the expression of genes controlling the motility, virulence and biofilm formation of several bacterial species [2–4]. Cellular levels of c-di-GMP are tightly regulated by the opposing activities of diguanylate cyclases (DGCs) [5–7] that synthesize c-di-GMP from two GTP molecules and phosphodiesterases (PDEs) [8,9] that hydrolyze c-di-GMP into 5'-phosphoguanlyl- (3'-5') - guanosine (pGpG) and GMP. A large number of genes bearing the PDE or DGC sequence signature motifs were identified in bacterial genomes, indicating a prominent role of c-di-GMP signaling pathways in regulating their community behavior [10].

Catalytically active DGCs bear the GG(D/E)EF consensus amino-acid sequence at their active site (A-site) and a conserved RxD motif at an inhibitory site (I-site) [11]. The I-site binds to the product and mediates feedback inhibition. The A- and I-sites are located at opposite ends of the GGDEF domain [11].

Crystal structures of several proteins containing DGC domains have been reported, such as PleD [5] from *Caulobacter crescentus*, DgcZ [12] from *Escherichia coli* K-12 and WspR [6,13] from *Pseudomonas aeruginosa* and *Pseudomonas syringae*. Isolated GGDEF domains structurally characterized include YfiN from *Pseudomonas aeruginosa* [14] and Xcc4471 from *Xanthomonas campestris* [15]. C-di-GMP synthesis from two separate GTP molecules requires dimerization of two active GGDEF domains, with the loops bearing the GGDEF motifs of each monomer brought in close proximity to initiate the cyclization reaction. Several proteins containing a DGC catalytic domain also possess an N-terminal domain such as PAS, GAF, CheY and HAMP, that control their diguanylate cyclase activity in response to upstream activating signals, often as part of two-component signaling systems [16]. Thus, the N-terminal domain of PleD plays a crucial role in orienting the two GGDEF domains in either a productive dimeric conformation for c-di-GMP synthesis or in a non-productive arrangement of the two half active sites, to prevent catalysis [7]. A similar mechanism is observed in DgcZ from *E. coli* where binding of Zinc to the N-terminal CZB domain induces an inactive dimeric arrangement of the GGDEF domains [12]. In the case of
WspR, inhibition proceeds through the formation of different oligomeric states mediated by a helical stalk domain⁴. Thus, two mechanisms can chiefly account for the feedback inhibition of DGCs mediated by the c-di-GMP product: allostery [11] and domain immobilization [7]. According to the allosteric model, binding of c-di-GMP to the enzyme I-site triggers structural changes across the protein, leading to a rearrangement of the A-site residues into an inactive conformation. In contrast, inhibition could proceed from the bridging action of c-di-GMP, which, when bound to the I-site, locks two GGDEF domains in a catalytically incompetent state, where the two A sites of the dimer are maintained apart from each other. This was observed in PleD, where an inactive dimer is stabilized by a bound c-di-GMP molecule on both monomers at the primary and secondary I-sites [17]. A similar inactive dimer was also observed for a GGDEF domain from *Marinobacter aquaeoli* [18].

Large-scale production of c-di-GMP would facilitate biological studies of several bacterial signaling pathways and phenotypes controlled by this molecule. Moreover, the immunomodulatory properties of c-di-GMP open up the possibility of using it as a vaccine adjuvant [19,20]. Chemical synthesis of c-di-GMP was achieved, but yields are not economical [21]. Recently, Rao *et al.* [22] cloned and expressed a GGDEF domain from *Thermotoga maritima*, hereafter named tDGC, that enables robust enzymatic production of large quantities of c-di-GMP. Enzyme assays showed that tDGC is active even after days of incubation at 30°C but with very low product production, irrespective of substrate concentration, due to feedback inhibition. However, a mutation (R158A) at the I-site of tDGC significantly alleviates production inhibition, enabling synthesis of several hundred milligrams of c-di-GMP with just 10 mg of the mutant enzyme [22].

In order to understand the structural determinants of tDGC thermostability and to investigate the catalytic mechanism and feedback inhibition of this biotechnologically important enzyme, we crystallized tDGC (i) in an active-like conformation with two A-sites in proximity, and (ii) in an inactive conformation, with c-di-GMP bound to the I-site. A comparison between these two structures enables us to examine potential conformational changes propagated between the I- and A- site upon c-diGMP binding to the I-site. We also report the structure of the R158A mutant of tDGC. Structural comparison with DGC homologues from mesophilic bacteria reveals the presence of a higher number of salt bridges in tDGC. Using circular dichroism (CD), three key salt bridges responsible for the enhanced thermostability of tDGC were identified.

## Materials and Methods

### Protein expression and purification

The gene encoding the complete (248 amino-acids) wild type protein TM1788 (Uniprot: Q9X2A8) from *Thermotoga maritima*, with codons optimized for bacterial expression, was purchased from Genscript. The gene corresponding to residues 82–248 (comprising the GGDEF domain, referred to as “tDGC protein”) was cloned into the pET28b vector between Nde1 and Xho1 restriction sites). The R158A feedback inhibition null mutant was purchased from Genscript. All proteins containing six Histidine residues at their N-terminus were overexpressed in *E. coli* BL21 (DE3) cells by the addition of 0.8 mM IPTG at 37°C for 4 hours. The tDGC protein was initially purified by immobilized metal affinity chromatography (IMAC) using Ni-NTA agarose (Qiagen) and further purified by ion exchange using a HiTrap column (GE Lifesciences) and size exclusion on a Superdex 75 column (GE

Lifesciences). The RocR PDE was expressed and purified as described [23]. Nucleotide-free tDGC protein was obtained by treatment with RocR in the presence of 5 mM MgCl₂. The mixture of RocR, tDGC and pGpG was separated by gel filtration. The identity of the eluted pGpG was confirmed by mass spectrometry.

### Protein crystallization

Purified tDGC was concentrated to 10 mg/ml and its homogeneity assessed by dynamic light scattering on a DynaPro instrument. Initial crystallization screening was performed by sitting-drop vapor diffusion at 18°C, using commercial crystallization screens from Hampton Research (Crystal Screens 1 and 2, PEG/Ion, Index and SaltRx). Intelli96–3 wells sitting drop plates were used to test three precipitant:protein ratios (1:1, 1:2 and 2:1) using a Phoenix crystallization robot (Art Robbins Instruments). Crystals of the inactive dimer were obtained with 0.1 M monobasic sodium phosphate monohydrate, 0.1 M monobasic potassium phosphate, 0.1 M MES monohydrate (pH 6.5), 2 M NaCl. The R158A mutant was crystallized with 2 M NaCl, 10% (w/v) PEG 6000. The active-like dimer was crystallized with 0.15 M lithium sulfate, 0.1 M citric acid (pH 4.0), 10% (w/v) PEG 6000. The protein used to obtain the latter crystal form was also the c-di-GMP-locked inactive dimer. However, upon mixing with precipitant, the crystallization drops turn milky white, indicating protein precipitation. Upon overnight incubation at 18°C, the crystallization drops become clear again and crystal formation is observed after 4–5 days. Thus, the active-like dimer is likely to originate from a rearrangement following initial protein precipitation.

The diffraction properties of these crystals were improved by overnight soaking in a precipitating solution containing 14% (w/v) PEG 6000 prior to flash-freezing in liquid nitrogen, extending the diffraction limits from 3.5 Å to 1.9 Å.

### Data collection and structure solution

X-ray diffraction data were collected at SLS (Vilingen, Switzerland) and processed using iMosflm [24]. The structure was determined by molecular replacement using the BALBES [25] server with the DGC structure from *Methylococcus capsulatus* (PDB code: 3ICL), which has 40% amino-acid sequence identity with tDGC, as a search probe. The structure was built by iterative cycles of model building at the computer graphics using COOT [26] and refinement using Autobuster [27]. The geometrical parameters for c-di-GMP were generated using the c-di-GMP coordinates 2V0N from the PDB (www.rcsb.org). Buried solvent accessible surface areas upon dimer formation were calculated using the PDBePISA server [28]. The quality of the structures was assessed using the MOLPROBITY [29] server and the figures were generated using the Pymol software [30]. Data collection and structure refinement parameters are summarized in **Tables 1 and 2**. The atomic coordinates and structure factors are deposited with the Protein Data Bank under accession codes 4URS (for the c-di-GMP bound inactive dimer), 4URG (for the c-di-GMP bound active-like dimer) and 4URQ (for the apo R158A mutant).

### Site-directed mutagenesis

For thermostability analysis, point mutations were introduced into the wild type tDGC GGDEF domain in the pET22b vector. The forward primers: 5' ACC ATT CTG CTG TAT GCC ATG AAA GAA GAA TAT 3' (for D177A; mutant codon underlined), 5' CTG AGC ACC TTT CGT GCC CCG GTG CGT GTG GAA 3' (for E196A) and 5' GAA AGT GGC GGA TAT GGC TAT GTA TAA AGC CAA AG 3' (for R233A) along with the
<table><caption>**Table 1.** Data collection statistics.</caption><thead><tr><th>tDGC</th><th>c-di-GMP bound inactive dimer</th><th>c-di-GMP bound active-like dimer</th><th>c-di-GMP free R158A</th></tr></thead><tbody><tr><td>Detector Type</td><td>Pilatus</td><td>Pilatus</td><td>MAR CCD</td></tr><tr><td>Synchrotron</td><td>SLS PXIII</td><td>SLS PXIII</td><td>NSRRC 13B1</td></tr><tr><td>Wavelength (Å)</td><td>1.00</td><td>1.00</td><td>1.00</td></tr><tr><td>Resolution (Å)</td><td>32.57–2.27 (2.45–2.27)</td><td>36.49–1.9 (36.49–1.9)</td><td>38.38–2.5 (2.56–2.5)</td></tr><tr><td>Space Group</td><td>P1</td><td>C2</td><td>C2</td></tr><tr><td>Molecules per asymmetric unit</td><td>2</td><td>2</td><td>6</td></tr><tr><td>Cell parameters (Å)</td><td>37.15, 37.14, 62.28</td><td>108.89, 52.12, 73.71</td><td>159.41 91.98 87.59</td></tr><tr><td>α/β/γ (°)</td><td>95.48, 100.23, 107.78</td><td>90, 124.05, 90</td><td>90, 90.04, 90</td></tr><tr><td>Measured reflections <sup>a</sup></td><td>62,597 (8,725)</td><td>184,283 (27,210)</td><td>323,030 (47,601)</td></tr><tr><td>Unique reflections <sup>a</sup></td><td>13,630 (1,971)</td><td>27,197 (3,951)</td><td>43,759 (6,374)</td></tr><tr><td>Redundancy</td><td>4.6 (4.4)</td><td>6.8 (6.9)</td><td>7.4 (7.5)</td></tr><tr><td>I/σ <sup>a</sup></td><td>6.5 (2.6)</td><td>11.4 (5.9)</td><td>15.5 (3.0)</td></tr><tr><td>Completeness (%)<sup>a</sup></td><td>99.6 (99.5)</td><td>100 (100)</td><td>99.7 (100)</td></tr><tr><td>R<sub>merge</sub> <sup>a</sup></td><td>0.143 (0.508)</td><td>0.118 (0.325)</td><td>0.088 (0.694)</td></tr></tbody></table>

respective complementary primers were used. The mutations were confirmed by sequencing and the constructs were overexpressed and purified following the same protocol that was used for the wild type tDGC protein.

## Melting temperature assay

Circular dichroism (CD) measurements were performed using a Jasco spectropolarimeter (J-810) with a 1 mm path length quartz cuvette. To determine the melting temperature, 200 µl protein

samples at a concentration of 0.8 mg/ml, in 20 mM Tris-HCl (pH 8.5), 250 mM NaCl, 5% (v/v) glycerol, 10 mM MgCl₂ were used. The temperature was gradually increased from 25 to 95°C at the rate of 1°C/min and ellipticity was recorded at a wavelength of 220 nm for every 0.2°C increment in temperature. The melting temperature was determined by plotting the CD data at 220 nm against temperature and fitting to a Boltzmann sigmoidal curve (between 60 and 95°C) using the curve fitting software XL Fit4.

<table><caption>**Table 2.** Refinement statistics.</caption><thead><tr><th>Tdgc</th><th>c-di-GMP bound inactive dimer</th><th>c-di-GMP bound active-like dimer</th><th>c-di-GMP free R158A</th></tr></thead><tbody><tr><td>Reflections used for refinement</td><td>13,625</td><td>27,172</td><td>43,743</td></tr><tr><td>Non-hydrogen atoms</td><td>2,980</td><td>2764</td><td>7531</td></tr><tr><td>Water</td><td>93</td><td>144</td><td>127</td></tr><tr><td>R<sub>work</sub> (%)</td><td>21.80 (19.50)</td><td>20.73 (20.79)</td><td>22.20 (23.21)</td></tr><tr><td>R<sub>free</sub> (%)</td><td>26.80 (27.13)</td><td>24.13 (21.09)</td><td>24.67 (25.43)</td></tr><tr><td>r.m.s.d. bond lengths (Å)</td><td>0.013</td><td>0.014</td><td>0.014</td></tr><tr><td>r.m.s.d. bond angles (°)</td><td>1.71</td><td>1.47</td><td>1.63</td></tr><tr><td>Mean B-factor (Å<sup>2</sup>)</td><td></td><td></td><td></td></tr><tr><td>Protein</td><td>38.8</td><td>25.2</td><td>53.3</td></tr><tr><td>Ligands:</td><td></td><td></td><td></td></tr><tr><td>C-di-GMP</td><td>45.6</td><td>10.2</td><td></td></tr><tr><td>Glycerol</td><td>47.5</td><td></td><td></td></tr><tr><td>MES</td><td>63.2</td><td></td><td>43.4</td></tr><tr><td>Water</td><td>38.9</td><td>31.1</td><td></td></tr><tr><td>Ramachandran Plot</td><td></td><td></td><td></td></tr><tr><td>Most favored (%)</td><td>98.42</td><td>98.69</td><td>95.72</td></tr><tr><td>Allowed(%)</td><td>1.58</td><td>1.31</td><td>4.17</td></tr><tr><td>Outliers (%)</td><td>0.99</td><td>0</td><td>0.11</td></tr><tr><td>PDB code</td><td>4URS</td><td>4URG</td><td>4URQ</td></tr></tbody></table>

<footnote><sup>a</sup>The values for the highest resolution shell are shown in parenthesis.</footnote>

<footnote>doi:10.1371/journal.pone.0110912.t001</footnote>

<footnote>doi:10.1371/journal.pone.0110912.t002</footnote>
# Results

## Monomer structure

We expressed and crystallized a domain of the TM1788 protein from *Thermotoga maritima* (248 residues, WP_010865403) spanning residues 82–248, hereafter named tDGC. The structures of tDGC and of its R158A mutant were determined in three crystal forms (i) tDGC in an active-like dimeric conformation (ii) tDGC in an inactive dimeric conformation (iii) R158A mutant in a monomeric conformation (**Tables 1 and 2**). The overall structure of the tDGC monomer features a central antiparallel β-sheet with the topology β2-β3-β1-β6, surrounded by five α-helices (**Fig. 1a**). The location of the conserved catalytic 167-GGDEF-171 motif (A-site) and the inhibitory I-site (158-RxxD-161) on loops connecting β2 and β3 and between α3 and β2 respectively, are depicted in **Fig. 1b**. tDGC shares a fold similar to adenylyl cyclases and the palm domain of polymerases, with a reaction mechanism involving two metal ions, which are observed in the structure of PleD, crystallized with GTPαS, but not in the present structures. According to this reaction scheme, Asp169 deprotonates the 3' hydroxyl group of the ribose of the GTP substrate to initiate an in-line nucleophilic attack on the α-phosphate of the target GTP, followed by elimination of the pyrophosphate moiety (**Fig. 1c**). The four residues responsible for coordinating the Mg²⁺ ion (Asp126 and Asp169) and for making contact with the guanine moiety (Asn134 and Asp143) are strictly conserved in tDGC, as shown in **Fig. 1c**.

## Inactive dimer mediated by c-di-GMP bound at the I site

The wild type tDGC protein co-purifies with c-di-GMP and co-crystallizes without any extraneous addition of ligand. The structure was determined at a resolution of 2.27 Å (**Tables 1 and 2**). In this crystal form, tDGC adopts a dimeric conformation, where the two monomers are related by a non-crystallographic dyad (**Fig. 2**). The buried accessible surface area between the two

monomers is 1,235 Å² and their interface is stabilized by one salt bridge and thirteen hydrogen bonds. I-site residues 158-RxxD-161 from the two monomers interact with two molecules of c-di-GMP that are mutually intercalated (**Fig. 2a**). This symmetrical arrangement is similar to the PleD structure [7], with the four guanyl bases of the two c-di-GMP molecules bound to the primary inhibitory sites of one molecule comprising the 158-RxxD-161 motif and a secondary inhibitory site R115 contributed by the second monomer. Residue R158 plays a key role in the interaction with both c-di-GMP molecules, via its guanidinium moiety that completes the set of stacking interactions (**Fig. 2b**). In this arrangement, the two monomers are locked in an inactive orientation with their A sites facing away from each other in a catalytically unproductive mode (**Fig. 2a**). The crystallographic dimer observed in this crystal form is consistent with gel filtration analysis showing that tDGC forms a dimer in solution (**Fig. 3a**). To demonstrate that dimerization is mediated by c-di-GMP, which co-purifies with tDGC, we incubated the protein with RocR, a PDE from *Pseudomonas aeruginosa*, to remove the bound c-di-GMP. Remarkably, following enzymatic treatment with RocR, the wild type tDGC protein elutes as a monomer (**Fig. 3b**).

## I-site mutant (R158A) structure

The R158A mutant structure was determined at a resolution of 2.5 Å (**Tables 1 and 2**). The protein does not establish large contacts in the crystal packing, which is consistent with its monomeric form in solution (**Fig. 3a**). When superimposed with the wild-type tDGC protein, the r.m.s.d. is 0.99 Å for 152 Cα atoms, suggesting the absence of large conformational changes between the two proteins. The R158A mutation abolishes binding of the c-di-GMP product and eliminates feedback inhibition [22]. From their gel filtration elution profile, the R158A mutant elutes as a monomer, whilst tDGC elutes as a dimer (**Fig. 3a**). The sole

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Domain architecture of tDGC. (a) Topology of the GGDEF domain of tDGC. The N- and C- termini of the polypeptide chain are indicated. The β-strands are depicted by pink arrows and α-helices by blue tubes. (b) Location of the A- and I-site on the structure of tDGC: The loops bearing the “GGDEF” motif at the A-site and “RxxD” motif at the I-site are shown in red and blue respectively. (c) Superposition of the A site of tDGC (cyan) and PleD (dark blue)⁷. GTPαS-Mg²⁺, bound to PleD is shown as sticks and green spheres respectively. Residues involved in metal ion binding and base recognition (represented as sticks and labeled according to tDGC and PleD numbering schemes) are strictly conserved between the two proteins.</figcaption></figure>

doi:10.1371/journal.pone.0110912.g001
<figure><img src="image_4.png" /><figcaption>**Figure 2.** The inactive tDGC dimer (a) The two monomers present in the asymmetric unit are colored in grey and yellow respectively with the two bridging c-di-GMP molecules shown as sticks. The views are along (left) and perpendicular (right) to the non-crystallographic dyad. (b) Stereoview of the dimerization mediated by two c-di-GMP molecules (labeled c-di-GMP1 and c-di-GMP2), bound at the I-site of tDGC, forcing the GGDEF domain in an inhibited conformation, with both A-sites facing away from each other. An omit map (blue mesh) with Fourier coefficients $2F_o - F_c$ where the c-di-GMP ligand was omitted from phase calculation is shown at 1σ contour level. Residues from the RxxD motif at the I-site forming hydrogen bonds with the bound c-di-GMP molecules are displayed as sticks and the distances between interacting atoms are displayed. The same color code as in panel **a** is used.
doi:10.1371/journal.pone.0110912.g002</figcaption></figure>

mutation R158A is thus able to eliminate the binding of the c-di-GMP ligand to the I-site, restoring the protein to its monomeric form, and preventing feedback inhibition (Fig. 3).

## tDGC in active-like dimer conformation

The structure of tDGC at pH 4, in an alternative dimeric form, was determined to a resolution of 1.9 Å (Tables 1 and 2). In this crystal form, two c-di-GMP molecules are bound to the I-site of a single tDGC monomer and the A-sites are brought close to each other (Fig. 4a). The absence of a second c-di-GMP molecule bound to the other monomer of the asymmetric unit, eliminates the locking of the monomers in a non-productive orientation. The buried interface area between the two monomers is 905 Å². In this novel dimeric arrangement, the two monomers are related by a rotation of 167°, a value close to the transformation thought to lead to an active 2-fold symmetric dimer, poised to react with two GTP substrate molecules to perform the cyclization reaction [7] (Fig. 4a). According to the postulated mechanism for the formation of c-di-GMP, each A-site binds to a GTP molecule to form a 2-fold symmetric active dimer with the A sites in proximity [7]. Despite several attempts using soaking or cocrystallization, we could not observe bound GTPαS in this crystal form. The crystals were obtained at pH of 4 at which the carboxylic groups of the

aspartate residues involved in metal binding and base recognition (Fig. 1c) are likely to be protonated. The rather acidic conditions of crystallization are therefore likely to hamper metal and substrate binding although the protein conformation looks active (Fig. 4). The structure of PldD⁷ bound to GTPαS was therefore used as a guide to position two GTP molecules in an optimized 2-fold symmetric dimer (Fig. 4b). Only minor structural movements are then required to positions the α-phosphate groups of the two GTP molecules near the 3' oxygen atoms of the ribose of the adjacent GTP, showing how nucleophilic attack can initiate the cyclization reaction (Fig. 4c).

## Comparison of thermostable tDGC with mesophilic homologues

The thermostable GGDEF domain of tDGC shares a conserved structure with homologous proteins from three mesophilic species: *Pseudomonas aeruginosa*, *Marinobacter aquaeolei* and *Geobacter sulfurreducens* (Table 1 in File S1) and one thermo-tolerant organism *Methylococcus capsulatus*¹⁶. Several molecular factors that potentially contribute to thermostability were examined (Table 1 in File S1). This comparison revealed that the increased thermostability of tDGC could be predominantly due to an excess of salt bridges (total of 5) compared to mesophilic homologues (1 or
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Size Exclusion Chromatograms (a) Gel filtration chromatogram using a S75 column, of the wild type tDGC protein (red) that elutes as a c-di-GMP-bound dimer and of the R158A mutant (green) that purifies as a monomer and is devoid of nucleotides.
(b) Following incubation of the wild type tDGC bound to c-di-GMP with RocR (in the gel filtration buffer containing 25 mM Tris-HCl, 300 mM NaCl, 5% (v/v) glycerol, 10 mM MgCl₂ for 2 hrs at room temperature), c-diGMP is converted to pGpG and tDGC elutes predominantly as a monomer. SDS PAGE analysis of the dimeric and monomeric fractions demonstrates that both peaks contain tDGC.
doi:10.1371/journal.pone.0110912.g003</figcaption></figure>

2) (Table 2 in File S1). Only salt bridges that connect residues far apart in the tDGC amino-acid sequence were considered, as protein unfolding would not disrupt salt bridges involving residues adjacent in the sequence. A sequence alignment (Fig. 5a) reveals that of the five salt bridges present, two are conserved across protein homologues (D91-R165 and D126-K238), whilst the three salt bridges: K118-D177, R152-E196 and R233-D219 (Fig. 5) are unique to the thermostable enzyme. This is consistent with the observation that several *Thermotoga maritima* proteins possess an excess of salt bridges, compared to their mesophilic counterparts [31,32].

## Contribution of salt bridges to thermostability

To analyze the contribution of each salt bridge to tDGC thermostability, single mutants D177A, E196A and R233A were made to disrupt the three unique solvent exposed salt-bridges (Fig. 5b). The three mutants eluted as dimers in gel filtration, indicating that c-di-GMP is bound to their I-site. To evaluate the stabilization effects brought about by nucleotide binding and dimerization in the thermal denaturation experiments, the nucleotide-free monomeric form of the wild type protein was generated by incubating with RocR [23], in order to digest the c-

di-GMP molecule bridging the two tDGC monomers. This digestion step was followed by gel filtration to remove RocR and pGpG. The melting temperatures of the wild type tDGC and the mutants were measured by circular dichroism at 220 nm (Fig. 6). The monomeric and dimeric forms of the wild type tDGC protein exhibit comparable melting temperatures with T<sub>m</sub> values of 85.5°C and 85.0°C respectively. Disruption of the R152-E196 salt bridge has the smallest effect on tDGC thermostability (decrease of 4.2°C). However, disruption of the other two salt bridges, K118-D177 and D219-R233, lead to large reductions in the melting temperature of tDGC by 12.3°C and 8.6°C, respectively (Fig. 6).

## Discussion

We report the structure of tDGC from *Thermotoga maritima*, an important enzyme for the large scale enzymatic synthesis of c-di-GMP and the first GGDEF domain from a hyperthermophile. In its inactive form, the tDGC dimer is closely superimposable to the GGDEF dimers from *Marinobacter aquaeolei* (3IGN)[18], PleD from *Caulobacter crescentus* (2V0N)⁷, WspR from *Pseudomonas syringae* (3I5A) [13] and WspR from *Pseudomonas aeruginosa* (3BRE) [6], suggesting a common mechanism of inhibition by
<figure><img src="image_6.png" /><figcaption>**Figure 4.** The active-like tDGC dimer and the cyclization reaction. (a) Structure of tDGC crystallized in a dimeric active-like conformation, with the two half active sites (loops containing the GGDEF motif colored in blue) facing each other. One monomer is colored yellow and the other grey. The single c-di-GMP molecule bound to the I-site of the monomer colored in yellow, is shown as sticks. (b) "Optimized" dimer where the residual transformation (Figure S2 in File S1) was applied to the monomer colored grey to generate a dimer with exact two-fold symmetry. The view is along the non-crystallographic dyad that runs across the two GTP molecules displayed as sticks. (c) Magnified view of the tDGC-GTP-Mg Michaelis complex, modeled on the basis of the "optimized" 2-fold symmetric dimer. The arrows indicate the nucleophilic attack of the 3' oxygen atom on the α-phosphate of the adjacent GTP. (d) Superposition of the "optimized" active dimer (this work, yellow and grey) and the c-diGMP cross-linked YdeH dimer (PDB code: 3TVK, Zaehringer and Schirmer) (colored in pink and teal).</figcaption></figure>

doi:10.1371/journal.pone.0110912.g004

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Structure-based sequence alignment of selected DGCs and salt-bridges unique to tDGC. (a) Sequence alignment between tDGC from *Thermotoga maritima* (this work), mesophilic DGC domains (3ICL, 4IOB, 3IGN and 3EZU). The sequence of the structurally characterized PleD (2V0N) from *Caulobacter crescentus* is also included for comparison. The secondary structure elements of tDGC are displayed above the sequences. Strictly conserved residues are in red boxes and chemically similar residues are colored in red. Residues D126 and D169 of tDGC are involved in divalent metal coordination. The A- and I-site residues are highlighted with green and black triangles respectively. The blue triangles indicate salt-bridge forming residues that are conserved across the five proteins. Red triangles indicate residues unique to tDGC involved in salt bridge formation. (b) Mapping of the salt bridges unique to the thermostable enzyme onto the 3D structure of tDGC: Lys118-Asp177, Arg233-Asp219 and Arg152-Glu196 were disrupted in the site-directed mutagenesis study.</figcaption></figure>

doi:10.1371/journal.pone.0110912.g005
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Temperature-dependent unfolding of tDGC and of tDGC single mutants recorded by CD spectropolarimetry. The mean residue molar ellipticity vs temperature for proteins tDGC (dimeric form), tDGC (monomeric form, following treatment with RocR), D177A, E196A and R233A is plotted. A significant difference in the melting profiles of the wild type and the single tDGC mutants (disrupting individual salt bridge mutants) is visible. The melting temperature T~m~, was determined by a Boltzmann sigmoid analysis (see text). doi:10.1371/journal.pone.0110912.g006</figcaption></figure>

domain immobilization. This dimeric form is mediated by two intercalated c-di-GMP molecules bridging the two monomers, such that its two half-active sites are placed in a non-productive orientation. The arginine and aspartate residues of the RxxD primary I-site motif are conserved in 57% of diguanylate cyclases [7]. In addition, 64% of proteins that contain this inhibitory motif have a secondary secondary binding site (R115) on an adjacent subunit that participates in locking the domains together. Thus, this mode of restricting interdomain mobility could constitute a widespread mechanism of regulation and feedback inhibition of GGDEF proteins devoid of other regulatory domains. In the case of tDGC, however, we did not examine possible contributions from its N-terminal 83 residues in regulating the formation of an active or inactive GGDEF dimer.

Based on simulations of bound and free forms of the GGDEF domains, it was previously proposed that product binding at the I-site could induce rearrangements of A-site residues, through a allosteric signal transmitted between these two distant sites [11]. To investigate if c-di-GMP binding at the I-site causes conformational changes at the A-site, the structures of the active-like and inactive dimer were compared (Figure S1 in File S1). The backbone and side-chain conformations are closely similar, suggesting the absence of signal propagation between the A- and I-sites, upon c-di-GMP binding. We therefore conclude that stand-alone GGDEF domains are inhibited predominantly through domain immobilization mediated by c-di-GMP bound at the I-site.

Using the tDGC active-like dimer as a guide, we were also able to propose a 2-fold symmetrical dimer conducive to GTP cyclization with a nucleophilic attack by the 3' oxygen groups of the ribose. Further structural and enzymatic studies are now needed to confirm this model.

Although various structural features could account for thermo-stability, a study that compared 94 proteins structures from the hyperthermophile *Thermotoga maritima* with protein homologues

from mesophilic organisms showed an increased number of salt bridges and compactness of the protein [33]. Another study reported that 73% of proteins from *Thermotoga maritima* have a higher “contact order” (average sequence distance between all pairs of contacting residues) compared to their mesophilic counterparts [34]. tDGC shares most structural parameters with its mesophilic counterparts (Table 1 in File S1), but it has three exclusive salt bridges. From the denaturation assay, the difference in the melting temperature of the wild type protein in its dimeric form and monomeric form are similar (T~m~ ca 85.0°C), indicating that the nucleotide-mediated dimerization has minimal influence on protein thermostability. The disruption of salt bridges K118-D177 and D219-R233 decreases T~m~ by 12.3 and 8.6°C, respectively. Interestingly, residues D177 and D219 belong to loop regions of tDGC, suggesting a greater contribution to thermostability of salt bridges located in flexible regions of the protein. In contrast, salt bridge R152-E196, which connects alpha helices α3 and α4 has a more modest influence on thermostability (T~m~ decrease of 4.2°C). In addition to salt-bridges, a contribution to the thermostability of tDGC could derive from a larger number of hydrophobic interactions [35] (Table 1 in File S1). This aspect was not studied here.

In conclusion, our study identified a dimeric active-like conformation for a thermostable GGDEF domain, where the relative orientation of both active sites is compatible with c-di-GMP formation. The R158A mutation in the I-site abolishes feedback inhibition by preventing formation of the inactive c-di-GMP-locked dimer. Our work also points to a few structural features accounting for enhanced thermostability of tDGC including the formation of three surface exposed unique salt bridges. Point mutations introducing salt bridges in solvent-exposed loops could constitute one possible approach to engineer more thermostable proteins.
# Supporting Information

**File S1 Table 1, Table 2, Figure S1 and Figure S2.**
(DOCX)

## References

1. Ross P, Mayer R, Weinhouse H, Amikam D, Huggirat Y, et al. (1990) The cyclic diguanylic acid regulatory system of cellulose synthesis in Acetobacter xylinum. Chemical synthesis and biological activity of cyclic nucleotide dimer, trimer, and phosphothioate derivatives. J Biol Chem 265: 18933-18943.

2. Romling U, Gomelsky M, Galperin MY (2005) C-di-GMP: the dawning of a novel bacterial signalling system. Mol Microbiol 57: 629-639.

3. Jenal U, Malone J (2006) Mechanisms of cyclic-di-GMP signaling in bacteria. Annu Rev Genet 40: 385-407.

4. Romling U, Galperin MY, Gomelsky M (2013) Cyclic di-GMP: the first 25 years of a universal bacterial second messenger. Microbiol Mol Biol Rev 77: 1-52.

5. Chan C, Paul R, Samoray D, Amiot NC, Giese B, et al. (2004) Structural basis of activity and allosteric control of diguanylate cyclase. Proc Natl Acad Sci U S A 101: 17084-17089.

6. De N, Pirruccello M, Krasteva PV, Bae N, Raghavan RV, et al. (2008) Phosphorylation-independent regulation of the diguanylate cyclase WspR. Plos Biology 6: 601-617.

7. Wassmann P, Chan C, Paul R, Beck A, Heerklotz H, et al. (2007) Structure of BeF3- modified response regulator PleD: implications for diguanylate cyclase activation, catalysis, and feedback inhibition. Structure 15: 915-927.

8. Christen M, Christen B, Folcher M, Schauerte A, Jenal U (2005) Identification and characterization of a cyclic di-GMP-specific phosphodiesterase and its allosteric control by GTP. J Biol Chem 280: 30829-30837.

9. Schmidt AJ, Ryjenkov DA, Gomelsky M (2005) The ubiquitous protein domain EAL is a cyclic diguanylate-specific phosphodiesterase: enzymatically active and inactive EAL domains. J Bacteriol 187: 4774-4781.

10. Galperin MY, Nikolskaya AN, Koonin EV (2001) Novel domains of the prokaryotic two-component signal transduction systems. FEMS Microbiol Lett 203: 11-21.

11. Christen B, Christen M, Paul R, Schmid F, Folcher M, et al. (2006) Allosteric control of cyclic di-GMP signaling. J Biol Chem 281: 32015-32024.

12. Zahringer F, Lacanna E, Jenal U, Schirmer T, Boehm A (2013) Structure and signaling mechanism of a zinc-sensory diguanylate cyclase. Structure 21: 1149-1157.

13. De N, Navarro MV, Raghavan RV, Sondermann H (2009) Determinants for the activation and autoinhibition of the diguanylate cyclase response regulator WspR. J Mol Biol 393: 619-633.

14. Giardina G, Paiardini A, Fernicola S, Franceschini S, Rinaldo S, et al. (2013) Investigating the Allosteric Regulation of YfiN from Pseudomonas aeruginosa: Clues from the Structure of the Catalytic Domain. PLoS One 8: e81324.

15. Yang CY, Chin KH, Chuah ML, Liang ZX, Wang AH, et al. (2011) The structure and inhibition of a GGDEF diguanylate cyclase complexed with (c-di-GMP)(2) at the active site. Acta Crystallogr D Biol Crystallogr 67: 997-1008.

16. Galperin MY (2006) Structural classification of bacterial response regulators: diversity of output domains and domain combinations. J Bacteriol 188: 4169-4182.

17. Schirmer T, Jenal U (2009) Structural and mechanistic determinants of c-di-GMP signalling. Nat Rev Microbiol 7: 724-735.

18. Vorobiev SM, Neely H, Yu B, Seetharaman J, Xiao R, et al. (2012) Crystal structure of a catalytically active GG(D/E)EF diguanylate cyclase domain from Marinobacter aquaeolei with bound c-di-GMP product. J Struct Funct Genomics 13: 177-183.

## Author Contributions

Conceived and designed the experiments: JL KS. Performed the experiments: AD CWL JL. Analyzed the data: AD JL. Contributed reagents/materials/analysis tools: ZL. Wrote the paper: AD JL.

19. Karaolis DK, Newstead MW, Zeng X, Hyodo M, Hayakawa Y, et al. (2007) Cyclic di-GMP stimulates protective innate immunity in bacterial pneumonia. Infect Immun 75: 4942-4950.

20. Karaolis DK, Means TK, Yang D, Takahashi M, Yoshimura T, et al. (2007) Bacterial c-di-GMP is an immunostimulatory molecule. J Immunol 178: 2171-2181.

21. Kiburu I, Shurer A, Yan L, Sintim HO (2008) A simple solid-phase synthesis of the ubiquitous bacterial signaling molecule, c-di-GMP and analogues. Mol Biosyst 4: 518-520.

22. Rao F, Pasunooti S, Ng Y, Zhuo W, Lim L, et al. (2009) Enzymatic synthesis of c-di-GMP using a thermophilic diguanylate cyclase. Anal Biochem 389: 138-142.

23. Chen MW, Kotaka M, Vonrhein C, Bricogne G, Rao F, et al. (2012) Structural insights into the regulatory mechanism of the response regulator RocR from Pseudomonas aeruginosa in cyclic Di-GMP signaling. J Bacteriol 194: 4837-4846.

24. Battye TG, Kontogiannis L, Johnson O, Powell HR, Leslie AG (2011) iMOSFLM: a new graphical interface for diffraction-image processing with MOSFLM. Acta Crystallogr D Biol Crystallogr 67: 271-281.

25. Long F, Vagin AA, Young P, Murshudov GN (2008) BALBES: a molecular-replacement pipeline. Acta Crystallogr D Biol Crystallogr 64: 125-132.

26. Emsley P, Lohkamp B, Scott WG, Cowtan K (2010) Features and development of Coot. Acta Crystallogr D Biol Crystallogr 66: 486-501.

27. Smart OS, Womack TO, Flensburg C, Keller P, Paciorek W, et al. (2012) Exploiting structure similarity in refinement: automated NCS and target-structure restraints in BUSTER. Acta Crystallogr D Biol Crystallogr 68: 368-380.

28. Krissinel E, Henrick K (2007) Inference of macromolecular assemblies from crystalline state. J Mol Biol 372: 774-797.

29. Davis IW, Leaver-Fay A, Chen VB, Block JN, Kapral GJ, et al. (2007) MolProbity: all-atom contacts and structure validation for proteins and nucleic acids. Nucleic Acids Res 35: W375-383.

30. Schrodinger LLC (2010) The PyMOL Molecular Graphics System, Version 1.3r1.

31. Merz A, Knochel T, Jansonius JN, Kirschner K (1999) The hyperthermostable indoleglycerol phosphate synthase from Thermotoga maritima is destabilized by mutational disruption of two solvent-exposed salt bridges. J Mol Biol 288: 753-763.

32. Pappenberger G, Schurig H, Jaenicke R (1997) Disruption of an ionic network leads to accelerated thermal denaturation of D-glyceraldehyde-3-phosphate dehydrogenase from the hyperthermophilic bacterium Thermotoga maritima. J Mol Biol 274: 676-683.

33. Robinson-Rechavi M, Alibes A, Godzik A (2006) Contribution of electrostatic interactions, compactness and quaternary structure to protein thermostability: lessons from structural genomics of Thermotoga maritima. J Mol Biol 356: 547-557.

34. Robinson-Rechavi M, Godzik A (2005) Structural genomics of thermotoga maritima proteins shows that contact order is a major determinant of protein thermostability. Structure 13: 857-860.

35. Pace CN (1992) Contribution of the hydrophobic effect to globular protein stability. J Mol Biol 226: 29-35.