Structural Biology and Crystallization Communications the Structure of Jann_2411 (duf1470) from Jannaschia Sp. at 1.45 a ˚ Resolution Reveals a New Fold (the Abate Domain) and Suggests Its Possible Role as a Transcription Regulator

The crystal structure of Jann_2411 from Jannaschia sp. strain CCS1, a member of the Pfam PF07336 family classified as a domain of unknown function (DUF1470), was solved to a resolution of 1.45 A ˚ by multiple-wavelength anomalous dispersion (MAD). This protein is the first structural representative of the DUF1470 Pfam family. Structural analysis revealed a two-domain organization, with the N-terminal domain presenting a new fold called the ABATE domain that may bind an as yet unknown ligand. The C-terminal domain forms a treble-clef zinc finger that is likely to be involved in DNA binding. Analysis of the Jann_2411 protein and the broader ABATE-domain family suggests a role as stress-induced transcriptional regulators.


Introduction
The complete genome sequences of hundreds of organisms are now known, each of which contains thousands of genes that have evolved to create the bewildering diversity of life. To understand this complexity at the molecular level requires the investigation of the function and structure of a vast number of proteins. A major goal of the Protein Structure Initiative (PSI; http://www.nigms.nih.gov/Initiatives/ PSI/) is to expand our knowledge of the protein universe by solving the structures of representative members of large as yet uncharacterized protein families. The Pfam database (Finn et al., 2008) contains over 2000 such families, termed domains of unknown function (DUFs), and understanding their structure will help to elucidate their function. Thus, to extend the structural coverage of proteins with uncharacterized biological function, we targeted the Pfam protein family DUF1470 and have determined the structure of the Jann_2411 gene product from Jannaschia sp. strain CCS1, an ecologically relevant marine proteobacterium found in coastal and open surface waters. Jann_2411 has a molecular weight of 20.7 kDa (residues 1-187) and a calculated isoelectric point of 6.6 and its crystal structure was determined using the semiautomated high-throughput pipeline of the Joint Center for Structural Genomics (JCSG; http://www.jcsg.org; Lesley et al., 2002) as part of the NIGMS Protein Structure Initiative. Structural analysis revealed a two-domain organization, with the Nterminal domain consisting of a new fold that we call the ABATE (for Alpha-Beta-hairpin-Alpha TandEm) domain and the C-terminal domain forming a treble-clef zinc finger that we have termed the CGNR zinc-finger domain after a characteristic sequence motif that is conserved in this family. Jan_2411 forms a dimer, with both monomers implicated in the formation of a putative DNA binding site, and analysis of its genomic context suggests a role for the ABATE-domain family in stress-induced transcriptional regulation.
The cloning junctions were confirmed by DNA sequencing. Protein expression was performed in a selenomethionine-containing medium, with suppression of normal methionine synthesis, using the Escherichia coli strain GeneHogs (Invitrogen). At the end of fermentation, lysozyme was added to the culture to a final concentration of 250 mg ml À1 and the cells were harvested. After one freeze-thaw cycle, the cells were homogenized in lysis buffer [50 mM HEPES pH 8.0, 50 mM NaCl, 10 mM imidazole, 1 mM Tris(2-carboxyethyl)phosphine hydrochloride (TCEP)] and passed through a Microfluidizer (Microfluidics). The lysate was clarified by centrifugation at 32 500g for 30 min and loaded onto nickel-chelating resin (GE Healthcare) pre-equilibrated with lysis buffer. The resin was washed with wash buffer [50 mM HEPES pH 8.0, 300 mM NaCl, 40 mM imidazole, 10%(v/v) glycerol, 1 mM TCEP] and the protein was eluted with elution buffer [20 mM HEPES pH 8.0, 300 mM imidazole, 10%(v/v) glycerol, 1 mM TCEP]. The eluate was bufferexchanged with HEPES crystallization buffer (20 mM HEPES pH 8.0, 200 mM NaCl, 40 mM imidazole, 1 mM TCEP) using a PD-10 column (GE Healthcare) and treated with 1 mg TEV protease per 15 mg eluted protein. The digested eluate was passed over nickelchelating resin (GE Healthcare) pre-equilibrated with HEPES crystallization buffer and the resin was washed with the same buffer. The flowthrough and wash fractions were combined and concentrated to 12.5 mg ml À1 by centrifugal ultrafiltration (Millipore) for crystallization assays. Jann_2411 was crystallized using the nanodroplet vapor-diffusion method (Santarsiero et al., 2002) with standard JCSG crystallization protocols (Lesley et al., 2002). Sitting drops composed of 200 nl protein solution mixed with 200 nl crystallization solution were equilibrated against a 50 ml reservoir at 277 K for 15 d prior to harvest. Initial screening for diffraction was carried out using the Stanford Automated Mounting system (SAM; http:// smb.slac.stanford.edu/facilities/hardware/SAM/UserInfo;  at the Stanford Synchrotron Radiation Lightsource (SSRL, Menlo Park, California, USA). The crystallization reagent consisted of 1.4 M sodium acetate and 0.1 M sodium cacodylate pH 6.5. A rodshaped crystal of approximate dimensions 150 Â 50 Â 50 mm was harvested for data collection. Glycerol was added to the crystal as a cryoprotectant to a final concentration of 20%(v/v). The diffraction data were indexed in monoclinic space group C2 (Table 1). The oligomeric state of Jann_2411 in solution was determined using a 1 Â 30 cm Superdex 200 size-exclusion column (GE Healthcare) coupled with miniDAWN static light scattering and Optilab differential refractive index detectors (SEC/SLS; Wyatt Technology). The mobile phase consisted of 20 mM Tris pH 8.0, 150 mM NaCl, and 0.02%(w/v) sodium azide. The molecular weight was calculated using ASTRA 5.1.5 software (Wyatt Technology).

Data collection, structure solution and refinement
Multiple-wavelength anomalous diffraction (MAD) data were collected at SSRL on beamline BL11-1 at wavelengths corresponding to the remote ( 1 ), inflection ( 2 ) and peak ( 3 ) of a selenium MAD experiment. The data sets were collected at 100 K with a MAR Mosaic 325 mm CCD detector (Rayonix, Evanston, Illinois, USA) using the Blu-Ice data-collection environment (McPhillips et al., 2002). The MAD data were integrated and reduced using XDS and scaled with XSCALE (Kabsch, 1993). Initial substructure solution was performed with SHELX (Sheldrick, 2008) and the phases were refined with SOLVE (Terwilliger & Berendzen, 1999), with a mean figure of merit of 0.38 (0.59-2.0 Å ) with two selenium sites. Density modification with RESOLVE (Terwilliger, 2003) was followed by automated model building using ARP/wARP (Cohen et al., 2004). Model completion and refinement were performed with Coot (Emsley & Cowtan, 2004) and REFMAC 5.5 (Winn et al., 2003) using the remote ( 1 ) data. The refinement included phase restraints from SOLVE and TLS refinement with four TLS groups per chain. Datacollection and refinement statistics are summarized in Table 1.

Identification of metal-binding sites
X-ray fluorescence emission peaks for selenium, arsenic, zinc and nickel were observed when the crystal was excited with X-rays 500 eV above the Se edge on SSRL beamline 11-1. In order to determine the identity of the metals at the individual sites in the structure, four additional data sets were collected on SSRL beamline 1-5. These data sets were collected to 2.9 Å resolution at wavelengths of 1.278, 1.292, 1.480 and 1.497 Å , which are above and below the zinc and nickel absorption edges. Data statistics are described in Table 2. Anomalous difference Fourier maps were calculated for each wavelength using the density-modified experimental MAD phases. The large changes in peak heights across the zinc or nickel absorption edge clearly showed that one site contained zinc and the other contained nickel. The integrated peak heights at the metal sites are shown with the peak heights for the selenium and sulfur sites as a reference in Table 3.  Table 1 Summary of crystal parameters, data-collection and refinement statistics for Jann_2411 (PDB code 3h0n).
Values in parentheses are for the highest resolution shell.
where F calc and F obs are the calculated and observed structure-factor amplitudes, respectively. § R free is the same as R cryst but for 5.0% of the total reflections that were chosen at random and omitted from refinement. } This value represents the total B that includes TLS and residual B components. † † Estimated overall coordinate error (Collaborative Computational Project, Number 4, 1994;Tickle et al., 1998).
The theoretical f 00 values at each wavelength are listed for comparison in Table 4.

Validation and deposition
The quality of the crystal structure was analyzed using the JCSG Quality Control server (http://smb.slac.stanford.edu/jcsg/QC). This server processes the coordinates and data using a variety of validation tools including AutoDepInputTool (Yang et al., 2004), MolProbity (Davis et al., 2007), WHATIF 5.0 (Vriend, 1990), RESOLVE (Terwilliger, 2003) and MOLEMAN2 (Kleywegt, 2000) as well as several in-house scripts and summarizes the output. Fig. 1(b) was adapted from an analysis using PDBsum (Laskowski et al., 2005) and all other figures were prepared with PyMOL (DeLano Scientific). Atomic coordinates and experimental structure factors for Jann_2411 have been deposited in the PDB (http://www.pdb.org) under accession code 3h0n.

Overall structure
The crystal structure of Jann_2411 (Fig. 1a) was determined to 1.45 Å resolution using the MAD method. Data collection and refinement statistics are summarized in Table 1. The final model included one protomer (residues 1-184), three acetate molecules, two glycerol molecules, one zinc ion, one nickel ion, two sodium ions and 240 water molecules in the asymmetric unit. Arg185, Ala186 and Thr187 at the C-terminus and Gly0 remaining after the cleavage of the expression/purification tag at the N-terminus were disordered and not modeled. Poor electron density was observed for the side chains of Lys45, Arg161, Asn162, Lys177 and Arg182. The side chains of the following residues were modeled in two conformations: Ile37, Asp42, Arg67, Asp74, His90, Gln103, Leu116, Glu119, Leu120 and Met123. The Matthews coefficient (V M ; Matthews, 1968) was 2.5 Å 3 Da À1 and the estimated solvent content was 50.9%. The Ramachandran plot produced by MolProbity (Davis et al., 2007) showed that 98.9% of the residues were in favored regions, with no outliers. Jann_2411 belongs to the Pfam family known as DUF1470, which accounts for the entire length of the protein sequence. However, the structure shows that Jann_2411 is actually comprised of two domains (Figs. 1a and 2a). The first domain (residues 1-142) can be visualized as two subdomains (H2-H4, 1-2 and H5-H7, 3-4) that share similar topology and secondary-structure elements, namely a helix--hairpin-helix motif (H2-1-2-H3 in the first subdomain; H6-3-4-H7 in the second subdomain), with an additional helix (H4 from the first subdomain and H5 from the second subdomain) linking the two motifs. We have therefore named this region the ABATE domain, representing the Alpha-Beta-hairpin-Alpha TandEm motif.
The helices of the first subdomain (H2 and H3) are stacked perpendicular to helices H5 and H7 which, together with H6, form a helical bundle capped on one end by the 3 10 -helix H8 and on the other by helix H4. In both subdomains, the -hairpin is oriented orthogonally with respect to the long axes of the helices in the ABATE motif. Superposition of the two ABATE motifs (residues 9-46 for ABATE1 and 80-139 for ABATE2) results in an r.m.s.d. of 3.2 Å over 27 residues (3% identity), which is non-significant. Additionally, sequence alignment shows little residue conservation in this domain, with no strictly conserved residues observed between representative ABATE-family sequences (Fig. 3), suggesting that this domain evolved as a single unit as opposed to the gene-duplication event that might be suggested by the presence of two ABATE motifs.
The second domain (residues 143-187; H8, 5-6, H9) forms a treble-clef zinc finger (Fig. 2). The zinc ion is coordinated by two cysteines (Cys147 and Cys152) from a loop termed the zinc knuckle, located between the strands of the third -hairpin (5-6), and two cysteines from the N-terminus of helix H9 (Cys168 and Cys172) (Fig. 2a). This arrangement of zinc-coordinating residues is typical of treble-clef zinc fingers (Grishin, 2001;Krishna et al., 2003). Other strictly conserved residues (Fig. 3) in this domain include Asp158 and Arg175. A high degree of conservation is observed for a number of positively charged residues (Arg143, Arg161, Arg165, Lys177, Arg182 and Arg184), suggesting that this region could present a nucleic acid binding site. Furthermore, residues 146 (a hydrophobic residue) and 167 (an aromatic residue) are highly conserved (Fig. 3) and could intercalate between the DNA bases. Based on the most conserved motif found in the C-terminal -helix in this family of proteins, we have named this domain the CGNR zinc finger. The actual amino-acid sequence in Jann_2411 is CQNR.  Table 3 Anomalous difference Fourier integrated peak heights.
The maps were calculated with data from 20 to 2.9 Å . The signal listed is the value reported by MAPMAN (Kleywegt & Jones, 1996) after integration of a sphere of radius 2 Å around the atom center from the final refined model. Sulfur and selenium sites are listed to provide a reference for differences in scale between different maps.   The above findings led to the re-evaluation of the Pfam DUF1470 family which, as a result of our study, will now be split into two entries in the next Pfam release (the current release is Pfam 23.0, July 2008). The original DUF1470 entry has been truncated and renamed to represent the ABATE domain, while a new Pfam family has been created for the C-terminal CGNR zinc-finger domain called zf-CGNR (Pfam accession PF11706).

Similarity to other proteins
SCOP classifies Jann_2411 as an + protein with an unusual fold (http://scop. mrc-lmb.cam.ac.uk/scop/data/scop.b.e.dda.b.b.b.html). A search with FATCAT (Ye & Godzik, 2004) using the N-terminal domain of Jann_2411 gave some hits involving helices H5-H7, but these structures [a three-helix bundle from a bacterial ATPase (PDB code 2v6y) and a histidine phosphotransferase domain (PDB code 1sr2)] involved less than half of the domain and displayed diverse functionalities. No single hit was found for the N-terminal domain in its entirety, leading us to propose that this domain represents a new fold. The C-terminal domain was structurally similar (main-chain r.m.s.d. of 2.5 Å over 40 residues with a sequence identity of 11%) to a plant homeodomain (PHD) finger from yeast (PDB code 2jmi), confirming the identity of this domain as a treble-clef zinc finger. Superposition of H9 onto the corresponding helix of the yeast structure revealed that the arrangement of the zinc ion and coordinating cysteines is conserved between the two structures (Fig. 2b).
Analysis of the crystallographic packing of Jann_2411 using the PISA server (Krissinel & Henrick, 2007) and analytical size-exclusion chromatography in combination with static light scattering indicate that a dimer is the likely quaternary form. The crystallographic dimer interface mainly involves hydrophobic contacts from the second -hairpin (strands 3-4), helices H6 and H7 and the intervening loops, and buries 990 Å 2 of surface area per monomer (Fig. 2a) Crystal structure of Jann_2411 from Jannaschia sp. strain CCS1. (a) Stereo ribbon diagram of the Jann_2411 monomer color-coded from the N-terminus (blue) to the according to the CastP server; Binkowski et al., 2003) along the dimer interface, delimited by the long loop connecting the last hairpin (strands 5-6) and helix H9. However, sequence and structure conservation in this domain is very weak, making quaternary states difficult to infer for the rest of the family.
Treble-clef zinc fingers are usually incorporated into larger structures and are found in proteins with a wide range of functions, many of which involve transcriptional regulation (Grishin, 2001;Krishna et al., 2003). Jannaschia sp. CCS1 is a member of the Roseobacter lineage, a taxon of marine bacteria. CCS1 is a phototroph that uses bacteriochlorophyll to harvest energy from light without the formation of oxygen. Genes predicted to have functional associations with Jann_2411 in the STRING database (http://string.embl.de) include a transmembrane protein of unknown function (Jann_2410) and the transcriptional regulator Jann_2412, a member of the Asr gene family. The Asr gene family is widespread in higher plants and most members of this family are up-regulated under a range of environmental stress conditions; their products are thought to function as transcriptional regulators (Frankel et al., 2006).
Other members of this newly defined ABATE protein family are found in plant symbionts (Rhizobium, Bradyrhizobium) and plant pathogens (Streptomyces, Ralstonia, Agrobacterium); they are around 180 residues in length and also contain the newly designated zf-CGNR (Pfam accession PF11706). Genome-location analysis of representative ABATE sequences shows co-occurrence with putative DNA-binding proteins, transcriptional regulators and membrane proteins implicated in ABC transport. In Streptomyces, several of the proteins co-occurring with ABATE homologs, such as sporulationspecific cell-division proteins and RNA polymerase sigma factors, are implicated in the control of aerial mycelium development (Dalton et al., 2007;Gordon et al., 2008) and are activated under conditions of cell-envelope stress, such as hyperosmolarity (Kormanec & Sevcikova, 2002).
A number of molecular mechanisms are shared between symbionts and pathogens, especially those involving host colonization and adaptation to a particular ecological niche (Hentschel et al., 2000). In both instances, changes in ecological and host environments necessitate fast adaptation strategies on the part of the microorganism. Zinc fingers possess the functional versatility necessary for this adaptation and have been exploited may times by nature and by the pharmaceutical industry (Papworth et al., 2006). The molecular function of the ABATE domain remains elusive. However, given the The C-terminal domain of Jann_2411 forms a zinc finger. (a) Ribbon representation of the Jann_2411 dimer. The zinc-finger domains are depicted in blue, the N-terminal domains are depicted in magenta and green and the zinc ions are shown as gray spheres. The zinc-coordinating cysteines are shown in ball-and-stick representation and labeled. (b) Stereoview of the structural superposition of the C-terminal domain of Jann_2411 (PDB code 3h0n, residues 144-187, gray) and a PHD finger fragment from yeast Yng1 protein (PDB code 2jmi, residues 38-83, blue). Zinc ions are shown as spheres and side chains of coordinating residues are indicated.
prediction that the C-terminal domain binds DNA, then the N-terminal domain may allow the protein to act as a signal-dependent transcriptional regulator with the ABATE domain conferring sensitivity to some as yet undefined ligand. This combination of a DNAbinding domain with a ligand-sensing domain is a prevalent form of regulation of operons in bacteria, such as the lactose or arabinose operons (Anantharaman et al., 2001). The likely dimeric nature of the protein and the long loop delimiting the suggested DNA-binding region hint at the possibility of an allosteric mechanism that might abolish DNA binding upon binding to an as yet unknown ligand. The LacI protein represents a canonical example for such an allosteric transition, with ligand-binding inducing a hinge-like motion that alters the relative subdomain orientations within the dimer, thereby changing DNA affinity (Lewis et al., 1996). Many other similar examples exist, such as TetR (Henssler et al., 2005;Premkumar et al., 2007;Koclega et al., 2007), all of which form dimers or higher order oligomers.
The availability of additional ABATE sequences and structures should shed light on the evolutionary history of this protein family. The information presented here, in combination with further biochemical and biophysical studies, should yield valuable insights into the functional role of Jann_2411. Models of Jann_2411 homologs can be accessed at http://www1.jcsg.org/cgi-bin/models/get_mor.pl? key=3h0nA.

Conclusions
The first structural representative of the DUF1470 family revealed a two-domain organization, with the N-terminal domain presenting a new fold and the C-terminal domain consisting of a treble-clef zinc finger. The structure additionally allowed a re-evaluation of the Pfam signature and the Pfam assignment and suggests a role for this family in stress-induced transcriptional regulation. Multiple sequence alignment of Jann_2411 and representative ABATE-family sequences from related species. Sequences were chosen from the DUF1470 Pfam seed alignment. The alignment was derived from the Pfam full alignment. UniProt abbreviations are as follows: Q28PN4_JANSC, gene locus Jann_2411 from Jannaschia sp. (strain CCS1); Q92QX7_RHIME, gene locus R01166 from Rhizobium meliloti; Q98CF3_RHILO, gene locus mlr5173 from R. loti; Q9L267_STRCO, gene locus SCO1542 from Streptomyces coelicolor; Q8UGB4_AGRT5, gene locus Atu1124 from Agrobacterium tumefaciens; Q8XT11_RALSO, gene locus RSp0306 from Ralstonia solonacearum; Q9RJT3_STRO, gene locus SCO0403 from S. coelicolor; Q89FX4_BRAJA, gene locus bII6575 from Bradyrhizobium japonicum; Q93J49_STRCO, gene locus SCO3054 from S. coelicolor. Residues are colored by conservation using the CHROMA software with default settings (Goodstadt & Ponting, 2001).