Oxyanion Binding Alters Conformation and Quaternary Structure of the C-terminal Domain of the Transcriptional Regulator ModE

IMPLICATIONS FOR MOLYBDATE-DEPENDENT REGULATION, SIGNALING, STORAGE, AND TRANSPORT*

David G. GourleyDagger §, Alexander W. SchüttelkopfDagger §, Lisa A. AndersonDagger , Nicholas C. Price||, David H. BoxerDagger , and William N. HunterDagger **

From Dagger  The Wellcome Trust Biocentre, University of Dundee, Dundee, DD1 5EH and the  Department of Biological Sciences, University of Stirling, Stirling, FK9 4LA, United Kingdom

Received for publication, January 31, 2001, and in revised form, March 19, 2001

    ABSTRACT
TOP
ABSTRACT
INTRODUCTION
MATERIALS AND METHODS
RESULTS AND DISCUSSION
REFERENCES

The molybdate-dependent transcriptional regulator ModE of Escherichia coli functions as a sensor of intracellular molybdate concentration and a regulator for the transcription of several operons that control the uptake and utilization of molybdenum. We present two high-resolution crystal structures of the C-terminal oxyanion-binding domain in complex with molybdate and tungstate. The ligands bind between subunits at the dimerization interface, and analysis reveals that oxyanion selectivity is determined primarily by size. The relevance of the structures is indicated by fluorescence measurements, which show that the oxyanion binding properties of the C-terminal domain of ModE are similar to those of the full-length protein. Comparisons with the apoprotein structure have identified structural rearrangements that occur on binding oxyanion. This molybdate-dependent conformational switch promotes a change in shape and alterations to the surface of the protein and may provide the signal for recruitment of other proteins to construct the machinery for transcription. Sequence and structure-based comparisons lead to a classification of molybdate-binding proteins.

    INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
MATERIALS AND METHODS
RESULTS AND DISCUSSION
REFERENCES

Molybdenum is an essential trace element required for the catalytic activity of several enzymes in animals, plants, and bacteria. In some cases the transition metal is complexed with a unique pterin forming the molybdopterin cofactor, and in others it forms part of an iron-molybdenum cluster cofactor (1, 2). Escherichia coli acquires molybdenum in the form of MoO<UP><SUB>4</SUB><SUP>2−</SUP></UP> by using a high-affinity ABC-type molybdate transporter system encoded by the modABCD operon (3). ModA, the structure of which has been determined (4, 5), is similar to the sulfate-binding protein of Salmonella typhimurium (SBP)1 (6) and a member of the periplasmic-binding protein family. ModA binds and transfers molybdate to ModB at the outer surface of the cytoplasmic membrane. ModB is an integral membrane protein and ModC a membrane associated protein; together they transport the oxyanion across the membrane using an ATP-dependent mechanism. The role of ModD is at present unknown. Exposure of E. coli to high levels of molybdate leads to repression of modABCD by a mechanism that involves the molybdate-dependent transcriptional regulator known as ModE (7-9). The modE gene, situated immediately upstream from the modABCD operon, codes for a transcriptional regulator able to control the uptake and utilization of this transition metal. ModE binds molybdate, and the complex can function as a repressor by binding at a site that overlaps the transcription start of the modABCD operon (10-12). In addition, the ModE-molybdate complex acts as a positive regulator of several genes including those coding for the molybdoenzymes dimethyl sulfoxide reductase (13) and nitrate reductase A (14), as well as for hyc, the hydrogenase 3 structural operon. It also mediates expression of the moaABCDE operon, which encodes the first enzymes required for the biosynthesis of molybdopterin (15).

E. coli ModE functions as a homodimer of subunits of 262 amino acids. The protein binds MoO<UP><SUB>4</SUB><SUP>2−</SUP></UP> and WO<UP><SUB>4</SUB><SUP>2−</SUP></UP> with a high affinity (Kd ~0.8 µM in each case) and stoichiometry of 2 oxyanions/dimer. Oxyanion binding produces a large quench (50%) of the intrinsic protein fluorescence accompanied by a blue-shift of the emission spectrum maximum, which indicates that a conformational change is induced by these ligands (11). The crystal structure of apo-ModE revealed that the protein forms distinct N- and C-terminal domains (Fig. 1, Ref. 16). The N-terminal domain, comprising 121 residues, is mostly alpha -helical (60% alpha -helix, 20% beta -strand) with a winged helix-turn-helix DNA binding motif. The C-terminal molybdate-binding domains, which display a pronounced asymmetry, comprise residues 122-262 and contain 60% beta -strands, which form two beta -barrels. The first barrel is constructed from a combination of strands beta 6-beta 9 with beta 15 and the second from strands beta 10-beta 14. Each barrel constitutes a subdomain (Fig. 1) of similar folds, although the order of secondary structure elements with respect to the sequence varies. The subdomain fold is related to the oligomer binding fold (OB-fold (17)), a five-stranded Greek-key beta -barrel capped by an alpha -helix. The subdomains share significant sequence and structural homology with each other and also with the 7-kDa molybdenum-containing Mop proteins from Clostridium pasteurianum (18, 19) and Sporomusa ovata (20). These proteins are implicated in the molybdenum metabolism of a variety of microorganisms (3, 11, 21). We have previously shown that each subdomain corresponds to the structure of a single Mop unit and termed the C-terminal domain of ModE the DiMop domain (16). The term molbindin has been introduced (19) to identify Mop-like proteins, and we made the further distinction that this refers to proteins composed of only Mop domains.


View larger version (41K):
[in this window]
[in a new window]
 
Fig. 1.   Ribbon diagram depicting the architecture of intact apo-ModE (16). Trp18, shown as a ball-and-stick representation, is labeled A186 and B186 to identify the subunit. The N-terminal DNA-binding domains are red and orange for subunits A and B, respectively, and the winged helix-turn-helix motifs are labeled (wHTHA, wHTHB). The molybdate-binding domains are divided into Mop subdomains labeled Mop1A, Mop2A, Mop1B, and Mop2B according to subunit.

Although the structure of ModE has been determined, there are a paucity of data concerning how the oxyanion ligands are bound, since we have been unable to determine the structure of the complete ModE-molybdate complex. Recent analysis of the Mop protein from S. ovata provided the first structural details of oxyanion binding to a molbindin, in this case tungstate (20), the results of which will be described later. To study the protein-oxyanion interactions relevant to the function of ModE we cloned and characterized the C-terminal (DiMop) domain. We show that this domain retains the ability to bind oxyanions selectively with properties similar to the intact protein, having determined high-resolution crystal structures of both the molybdate and tungstate complexes. These structures provide high-resolution information on the protein-ligand interactions and allow us to compare the apo-DiMop component of ModE with the oxyanion-loaded domain, a comparison that identifies a conformational change induced by oxyanion binding. Sequence and structural comparisons with other Mop proteins led us to propose a classification of the proteins that utilize the Mop domain.

    MATERIALS AND METHODS
TOP
ABSTRACT
INTRODUCTION
MATERIALS AND METHODS
RESULTS AND DISCUSSION
REFERENCES

Expression of the C-terminal DiMop Domain of ModE-- A fragment of the modE gene coding for the C-terminal domain was amplified by polymerase chain reaction using pHW121 as the template plasmid, which carries a fragment covering the modE region in the vector pUC18 (8). The oligonucleotide primer used for the 5'-end of the C-terminal domain-encoding DNA was 5'-G CGC CAT ATG CAA ACC AGC-3', which incorporated an NdeI recognition site (underlined) and covered the first 3 codons of the C-terminal domain. The primer used for the 3'-end of the C-terminal domain was 5'-GCG CGG ATC CCA CGC TTA GCG CAG-3', which covered the last 6 codons of modE and contained a BamHI recognition site (underlined). The polymerase chain reaction-amplified fragment was cloned as an 0.45-kilobase NdeI-BamHI insert in the N-terminal hexa-His tag expression vector pET15b (Novagen) to create plasmid Plaa7. The expression of protein is under the control of the T7-Lac promoter, and the recombinant DiMop domain carries a thrombin cleavage site between the protein and the His tag.

Plasmid Plaa7 was heat shock-transformed into E. coli BL21(DE3), and transformants were selected on Luria Bertani (LB) agar plates containing 100 µg ml-1 ampicillin. A transformant was grown at 37 °C in LB broth with ampicillin to mid-log phase at which point 0.4 mM (final concentration) isopropyl-beta -D-thiogalactopyranoside was added and growth continued with vigorous aeration for a further 4 h. Cells were harvested by centrifugation (2500 × g) at 4 °C and were then resuspended in 50 mM Tris-HCl, pH 7.6, 250 mM NaCl, 5 mM benzamidine and stored at -20 °C. Bacterial cells were broken by passage through a French press, the insoluble cell debris pelleted by centrifugation at 4 °C (18000 × g) for 15 min, and the cell extract passed through a 0.2-µm filter and applied to a 5-ml metal chelate affinity column (Hi-Trap, Amersham Pharmacia Biotech) charged using nickel chloride. The column was washed with four column volumes of 50 mM Tris-HCl, pH 7.6, 250 mM NaCl, and 5 mM benzamidine. The His-tagged protein was eluted with a 0-500 mM imidazole gradient in the same buffer and then incubated with thrombin (Amersham Pharmacia Biotech) for 12 h at 20 °C to remove the histidine tag. The DiMop domain was separated from the thrombin, uncleaved fusion protein, and N-terminal peptide by strong anion exchange chromatography using a RESOURCETM Q (Amersham Pharmacia Biotech) column on a BioCAD 700E work station. Pooled fractions were concentrated to 20 mg ml-1 in 50 mM Tris-HCl, pH 7.6, and sample purity was assessed by SDS-polyacrylamide gel electrophoresis and MALDI-TOF (matrix-assisted laser desorption time-of-flight) mass spectrometry (Voyager DE STR; PerSeptive Biosystems). Protein concentration was estimated spectrophotometrically at 280 nm using the calculated molar extinction coefficient of 12660 M-1 cm-1 (22). The yield of purified protein was ~40 mg liter-1 of bacterial culture. Following thrombin cleavage the DiMop domain consisted of residues 124-262 of ModE with a Gly-Ser-His extension at the N terminus. For reasons of consistency we maintained the amino acid numbering of the intact ModE.

Fluorescence Measurements-- Intrinsic protein fluorescence spectra were recorded using a PerkinElmer LS50B spectrofluorimeter maintained at 25 °C. Excitation was at 295 nm (5 nm bandpass for both excitation and emission) at a scan rate of 100 nm min-1 with a protein concentration of 29 µM of monomer in 50 mM Tris-HCl, pH 7.6. Binding parameters were determined by changes in fluorescence at 347 nm on successive additions of sodium molybdate to a solution of protein (38 µM monomer). Corrections were made for dilution effects but not for the inner filter effect because molybdate does not absorb significantly at the wavelengths used.

Crystal Growth and Data Collection-- Crystals were grown at 20 °C using the hanging drop method by mixing 3 µl of protein solution (16 mg ml-1 protein, 50 mM Tris-HCl, pH 7.6, 10 mM Na2MoO4 or Na2WO4) with 3 µl of a reservoir solution. Similar crystals were obtained with reservoirs of 400 µl of 1.0 M trisodium citrate, 100 mM HEPES, pH 7.6, and 10% mM ethylene glycol or 1.4 M ammonium sulfate in 100 mM Tris-HCl, pH 8.4. The crystals, which attained a size of 0.3 mm in all dimensions over a period of 2 days, were cryoprotected with ethylene glycol and cooled to -173 °C, and the data were then measured using an RAXIS IV image plate with a RU-200 rotating anode (Cu Kalpha lambda  '3d 1.54Å). Two data sets were measured, one from a crystal grown in the presence of WO<UP><SUB>4</SUB><SUP>2−</SUP></UP> using tri-sodium citrate as the precipitant and the other from a crystal obtained in the presence of MoO<UP><SUB>4</SUB><SUP>2−</SUP></UP> using ammonium sulfate as precipitant. The crystals are tetragonal with space group P41 with two DiMop domains in the asymmetric unit and about 42% solvent content. All diffraction data were processed, reduced, and scaled using the HKL suite (23) and were then analyzed using CCP4 software (24) with 5% of the reflections used for calculation of the R-free. Further details are given in Table I.

The structure of the DiMop-WO<UP><SUB>4</SUB><SUP>2−</SUP></UP> complex was solved by molecular replacement with the program AMORE (25) using a polyalanine search model constructed from a monomer of the DiMop domain from apo-ModE (16). An optimal position for one molecule was found by using an integration radius of 20 Å and all data between 15 and 3 Å resolution (correlation coefficient 0.57). Positioning the second molecule in the asymmetric unit gave a correlation coefficient of 0.56. The resulting electron density maps indicated significant structural changes between the search model and the oxyanion complex. The tungstate positions were identified and incorporated into the model, which was then used to generate phases for use in the warpNtrace mode of ARP/wARP (26). Data to 1.9 Å were used, and a new model containing 268 of 280 residues with an R-factor of 25.1% and an R-free of 27.7% was obtained.

Interactive computer graphics model building was carried out using the program O (27). A Ni2+ ion was located in the same position of the tungstate complex as seen in one of the subunits of the apo-ModE structure and is considered an artifact of the purification procedure (16). Refinement was completed using REFMAC (28) combined with ARPP (29) to place the water molecules. The molybdate complex was solved using the structure of the tungstate-bound form with the water molecules removed and the tungstate replaced with molybdate as a starting model for calculating phases for use in the auto model building procedure in ARP/wARP. The model was refined in a procedure identical to that used for the tungstate complex; relevant statistics are presented in Table I. The molybdate and tungstate complexes are essentially identical, with an r.m.s. deviation of 0.3 Å for superposition of all Calpha atoms. We concentrated on the molybdate-bound structure, although for completeness some details of the tungstate complex are provided. The figures were produced using ALSCRIPT (30), MOLSCRIPT (31), and RASTER-3D (32).

                              
View this table:
[in this window]
[in a new window]
 
Table I
Data collection, processing, refinement, and model geometry statistics for both complexes

Coordinates-- Atomic coordinates and structure factors have been deposited with the Research Collaboratory for Structural Bioinformatics Protein Data Bank (accession codes 1H9R/1H9S for the tungstate and molybdate complexes, respectively).

    RESULTS AND DISCUSSION
TOP
ABSTRACT
INTRODUCTION
MATERIALS AND METHODS
RESULTS AND DISCUSSION
REFERENCES

Properties of the DiMop ModE Fragment-- The C-terminal DiMop fragment of ModE was expressed and purified in high yield, and analytical gel filtration analysis (data not shown) revealed it to be a dimer in solution. The DiMop fragment retains two of the three tryptophan residues of ModE (Trp131, Trp186), and the intrinsic protein fluorescence emission spectrum is similar to that reported for intact ModE with lambda max at 347 nm. Upon addition of 0.5 mM sodium molybdate to the DiMop fragment, the fluorescence emission spectrum was reduced by more than 60% of its value in the absence of molybdate (data not shown). This molybdate-induced quench is slightly greater than that reported for intact ModE (50% (11)), clearly demonstrating that the DiMop fragment retains the ability to bind molybdate. Moreover, it suggests that the third tryptophan (Trp49) of ModE, located in the DNA-binding domain, makes a less important contribution to the molybdate-induced fluorescence quench of ModE.

Titration of the DiMop fluorescence quench provided an estimate of the stoichiometry of molybdate bound per monomer fragment of 1 with a Kd of ~0.5 µM (Fig. 2). These values are in excellent agreement with those reported for molybdate binding to intact ModE; stoichiometry of 1 molybdate per monomer with a Kd of 0.8 µM. The DiMop fragment effectively retains the molybdate binding capacity and properties of intact ModE; this observation validated our approach to the exploration of ModE molybdate interactions by crystallographic analysis of the ligand-bound DiMop fragment.


View larger version (8K):
[in this window]
[in a new window]
 
Fig. 2.   Titration of the DiMop fragment with molybdate. The continuous curved line is the calculated curve for a stoichiometry of 1 molybdate bound/monomer with a Kd of 0.5 µM.

The Molybdate Binding Site-- The DiMop dimer binds two oxyanions, in sites designated I and II, at the dimerization interface (Fig. 3) using residues from both polypeptide chains. Each ligand binding site is created by the turn between beta 5 and beta 6 from one DiMop and helix alpha 6 and the C-terminal region of strand beta 15 from the other (Fig. 3B). The oxyanion is held in position by accepting nine hydrogens bonds with each oxygen participating in at least two such interactions. Five of the hydrogen bonds are donated from one polypeptide, three from the other, and one from a water molecule (Table II, Fig. 4). The MoO<UP><SUB>4</SUB><SUP>2−</SUP></UP> interacts with a single basic residue, Lys183. Other interactions from the same subunit involve hydrogen bonds donated from the main chain amides of Arg128, Ala184, and the side chain OH of Ser126. The MoO<UP><SUB>4</SUB><SUP>2−</SUP></UP> is bound to the other subunit by hydrogen bonds donated from Ser166 OH and both the OH and the main chain amide of Thr163. Sequence comparisons showed that Ser126, Arg128, Thr163, and Lys183 are conserved in all Mop-like sequences (Fig. 5), and it was proposed that the residues were involved in oxyanion binding (19). Arg128 binds the oxyanion using the main chain amide. The side chain forms an intramolecular salt bridge to Glu218 (not shown) and contributes to the stabilization of the tertiary structure and formation of the ligand-binding pocket.


View larger version (41K):
[in this window]
[in a new window]
 
Fig. 3.   The dimeric assembly and architecture of the DiMop domain of ModE complexed with molybdate (shown as van der Waals spheres with molybdenum in cyan and oxygen in red) occupying binding sites I and II. The same color scheme is used as in Fig. 1. A, view with Mop1B and Mop2B in a similar orientation to the corresponding subunit shown in Fig. 1. B, view orthogonal to A with selected elements of secondary structure labeled.

                              
View this table:
[in this window]
[in a new window]
 
Table II
Hydrogen bond lengths (Å) formed between the oxyanions and the protein


View larger version (30K):
[in this window]
[in a new window]
 
Fig. 4.   Stereoview of oxyanion binding site I. Nitrogen and oxygen atoms of the protein are shown as blue and red spheres, respectively, and a water molecule is shown as a light blue sphere. Residues are colored blue if part of Mop1A and green where they are part of Mop1B. Dotted black lines represent hydrogen bonds, and the close contact between the carbonyl oxygen of Ser126 and molybdenum is depicted as a yellow dotted line. The single-letter code is used for amino acids and prime to further differentiate those residues from Mop1B.


View larger version (30K):
[in this window]
[in a new window]
 
Fig. 5.   Sequence alignment of 13 Mop domains. The residues of ModE directly involved in oxyanion binding are marked with triangles, shaded to indicate the use of a side chain functional group and open where a main chain amide is used for binding. Residues enclosed in black are conserved in at least 11 of the sequences, those enclosed in gray correspond to at least 8 of 13 being conserved. The order of the sequences shown is as follows: The ModE Mop1 from E. coli, H. influenzae, and the two Mop domains from A. vinelandii ModA (a ModE homologue), the single Mops from C. pasteurianum, H. influenzae, and S. ovata, A. vinelandii ModG Mop1 and 2, E. coli ModC, and A. vinelandii ModD. The final two sequences are the Mop2 domains of E. coli and H. influenzae ModE, respectively.

The charge compensation for the oxyanion is achieved by interactions with a single basic side chain (Lys183) and partial charges contributed from the three main chain NH groups in the binding site. A strong peak in the electron density near the molybdate was assigned as a water molecule, rather than a monovalent cation such as sodium, on the basis of a water-like hydrogen bonding pattern. This water forms hydrogen bonds with the side chain OD1 of Gln144 from one subunit and the carbonyl oxygen of Thr232 of the partner subunit. Of note is a particularly short contact, 3.0 and 3.1 Å, between the metals in each site and the carbonyl oxygens of Ser126 (Fig. 4). The carbonyl is approximately equidistant from three of the molybdate oxygens and may help to align the oxyanion for binding. In S. ovata Mop (20) we note that the same structural feature is observed in the corresponding tungstate binding sites.

Size Determines Oxyanion Selectivity in ModE-- ModE functions by binding molybdate or tungstate but not other tetrahedral oxyanions such as phosphate, sulfate, or vanadate (9). In this work, DiMop retained bound MoO<UP><SUB>4</SUB><SUP>2−</SUP></UP> in the presence of 1.4 M ammonium sulfate used in the crystallization solutions. For selective binding of MoO<UP><SUB>4</SUB><SUP>2−</SUP></UP>, ModE can exploit differences in charge (PO<UP><SUB>4</SUB><SUP>3−</SUP></UP> and VO<UP><SUB>4</SUB><SUP>3−</SUP></UP>). The highly charged anions like VO<UP><SUB>4</SUB><SUP>3−</SUP></UP> are disadvantaged because ModE has only a single positively charged residue in the binding site. Size will be an important consideration in the case of PO<UP><SUB>4</SUB><SUP>3−</SUP></UP> and SO<UP><SUB>4</SUB><SUP>2−</SUP></UP>. Previous analysis of oxyanion geometry in small molecule structures have shown the average metal-oxygen distances of molybdate and tungstate are 1.75 ± 0.04 Å and 1.76 ± 0.002 Å, respectively, whereas the average sulfur-oxygen bond length in sulfate is 1.47 ± 0.02Å (5). Consequently MoO<UP><SUB>4</SUB><SUP>2−</SUP></UP> and WO<UP><SUB>4</SUB><SUP>2−</SUP></UP> are similar in size but much larger than SO<UP><SUB>4</SUB><SUP>2−</SUP></UP>. The program VOIDOO (33) was used to calculate the volume of the binding pocket accessible to a rolling probe of radius 1.4 Å. The value obtained for the DiMop domain is 78 Å3, which is comparable with the 72 Å3 obtained for the molybdate binding site in E. coli ModA (4) and Azotobacter vinelandii ModA (5). The volume of these molybdate-binding pockets is larger than that calculated for the sulfate-binding pocket in SBP, 59 Å3. The size of the oxyanion must therefore be important for selectivity.

ModA also exhibits a similar selectivity of tungstate and molybdate over other oxyanions even though it has a fold and binding site similar to SBP (4). ModA binds the oxyanion in a deep cleft formed at the junction of two globular domains. There are seven direct hydrogen bonds, all donated from main chain amides and side chain hydroxyls of neutral protein ligands (4, 5), and a contribution to binding from N-terminal alpha -helix dipoles. Although the DiMop domain is a different fold compared with ModA, both proteins use three OH groups from the Ser or Thr residues in ligand binding and have residues that bind ligands via a chelate-type interaction using both main chain nitrogen and side chain OH groups.

Molybdate-specific proteins, irrespective of fold, have larger binding sites and use size in conjunction with charge compensation to determine oxyanion selectivity. The mode of dimerization imposes size restrictions because the binding of a smaller oxyanion would require readjustment of the protein backbone and dimerization interface to optimize the length and angles of the hydrogen bonds involved in ligand binding.

The Overall Structure of the Ligand-bound DiMop Domain and Comparisons with the Apo Form-- Like intact ModE, the ligand-bound DiMop forms a homodimer with a secondary structure similar to the apoprotein (Figs. 1 and 3). An obvious difference between the two structures is a decrease in the asymmetry between the partner subunits that was observed in the apoprotein structure (16). A least-squares superposition of the Calpha atoms of chain A onto chain B gives an r.m.s. deviation of 1.0 Å for the apoprotein, whereas for the DiMop-molybdate complex the r.m.s. deviation is lowered to 0.6 Å. Residues 138-144 were excluded from this superposition because they are disordered in chain B of apo-ModE.

The binding of molybdate induces a change in the quaternary structure that, as a first approximation, can be described as the rigid body movement of the four Mop domains against each other; when aligned on Mop1A (as described in the legend for Fig. 6), the center of mass of Mop1B moves about 3.0 Å between the apo and MoO<UP><SUB>4</SUB><SUP>2−</SUP></UP>-bound state. As this movement is almost perpendicular to the line connecting the centers of mass of Mop1A and -1B, the distance between the Mop1s remains practically constant. In the same alignment the center of mass of Mop2A moves about 1.9 Å, mostly toward Mop1B on the other chain. Overall these movements lead to the molybdate-bound protein being more compact than the corresponding part of the apo structure, with the buried surface area increasing from 2 × 750.5 Å2 for the apo-dimer to 2 × 1442.5 Å2 for the ligand-bound form. For consistency residues B138-B144 were removed from the complex structure prior to the calculation. The reduction of exposed surface is caused only in part by the rigid body movements; a significant contribution is also made by a number of loops that move relative to their surroundings. The largest change (with Calpha movements up to 12.6 Å) is observed for residues 138-144 (loop I, Fig. 6). As mentioned above, this loop is partially disordered in the apoprotein with no density observed for one chain. However, it becomes ordered upon molybdate binding, moving to cover another flexible loop formed by residues 163-166 (loop II). This smaller loop is important because its movement is crucial for the formation of the molybdate binding site, one-half of which is constituted by residues 163 and 166 (see above) after the loop has moved approx 2 Å toward the anion. This movement closes the molybdate binding site and leads to the formation of new hydrogen bonds across the intersubunit interface, e.g. between the backbone oxygen of Ser126 and the hydroxyl of Ser166', as well as the NH2 of Arg169' or between the guanidino moiety of Arg128 and the Thr163' hydroxyl (prime is used to denote a residue on the partner subunit). The latter two residues are also involved in molybdate binding and are conserved in virtually all known Mop domains (Fig. 5). Residues 152-155 form a third mobile loop (loop III). In the apoprotein it contacts the symmetry-related loop of the Mop2 of the same chain (residues 224-226), and its movement upon molybdate binding follows, to a limited extent, that of the Mop2. We noted that all mobile loops are located on the Mop1 domain and that Mop2 actually behaves like a rigid body.


View larger version (72K):
[in this window]
[in a new window]
 
Fig. 6.   Alignment of the ModE molybdate-binding domain in its apo form (blue) and oxyanion-bound form (green). Two "intermediate" states that were generated by Cartesian morphing (using the program LSQMAN (35)) are shown to illustrate the transition. The two proteins were aligned on the Calpha atom of the rigid core of Mop1A and -B comprising residues 124-137, 145-151, 156-162, and 167-184. The color gradient from blue-green to red reflects increasing Calpha distances based on the separate alignments of the rigid core for every Mop domain. Even though Mop2s are essentially rigid bodies, for comparability they were aligned on the Calpha atoms of residues 196-209, 217-223, 227-233, and 238-253. Unless stated otherwise all alignments and center-of-mass calculations use these sets of residues. The molybdate binding sites are labeled as in Fig. 3.

We can now describe the probable order in which molybdate binding and the subsequent structural changes take place. We propose that molybdate initially binds to the half-site formed by residues 126, 128, 183, and 184 of one protein chain. The bound molybdate then effects the movement of loop II of the partner chain and with that consequently the rigid body movement of the partner subunit. This seems a reasonable assumption, as the binding site components of residues 126, 128, 183, and 184 are practically unchanged in the two structures and can be considered a preformed binding site, whereas loop II moves and is clearly influenced by molybdate binding. Further support for this half-site being responsible for initial ligand binding comes from the unequal distribution of protein-ligand interactions; as described above this is where three of the four oxygens of the molybdate are bound, whereas the loop II component provided by the other subunit binds only one.

The closing of the binding site by the movement of loop II then could allow the hitherto solvent-exposed and flexible loop I to rearrange itself and cover loop II, which leads to the formation of a number of new hydrogen bonds, e.g. Gln145' backbone oxygen and nitrogen with Ile162' backbone nitrogen and oxygen, respectively, or Val143' backbone oxygen with Ala164' backbone N, and also improves the hydrophobic packing. The now ordered loop I might then attract the loop formed by residues 210-216 of the other chain, thereby effecting the observed movement of that Mop2 domain and leading to the formation of further hydrophobic contacts. Finally the loop III of the other chain moves to "follow" the Mop2 domain movement.

The fluorescence perturbation experiments indicate that the environments of one or both tryptophans in the DiMop domain change on binding MoO<UP><SUB>4</SUB><SUP>2−</SUP></UP>, and the crystal structures provide an explanation. When the apo and ligand-bound forms of the DiMop domains are compared, there is no significant change in the environment of Trp131 (not shown), which is located on strand beta 6 and positioned in a cavity formed by the residues linking beta 7 and beta 8. Trp186, however, is at the dimerization interface on a 310 helix. In apo-ModE the Trp186 of subunit A is aligned parallel to Trp186 from the B subunit and separated by about 6 Å (Fig. 1). The adjustment of the DiMop domains on binding MoO<UP><SUB>4</SUB><SUP>2−</SUP></UP> alters the side chain conformation, renders the tryptophans more buried, and reduces their distance to 3.3 Å (Fig. 3B). In apo-ModE, Trp186 contributes 1% of the DiMop dimerization interface, but in the ligand-bound form this now increases significantly to 5% of a greater surface area. We therefore attribute the spectral changes observed on binding oxyanion to alterations associated with the environment of Trp186. It follows then that changes in the environment of Trp186 on the addition of molybdate lead to the changes observed in the near UV CD spectrum noted for ModE (11).

Implications for Other Mop Proteins-- Mop domains occur in a variety of proteins other than ModE; the molbindins consist of either one (e.g. Mop from Hemeophilus influenzae or Mop I-III from C. pasteurianum) or two (e.g. ModG from A. vinelandii) Mop domains. The function of these small proteins is unclear, but it has been suggested that they store molybdate/tungstate and are thus involved in maintaining homeostasis of these anions. As molybdate binding occurs at the interface between two Mop units, it is obvious that functional molbindins must be at least dimeric. The structure of a single Mop protein from S. ovata (20) revealed a hexameric arrangement (Fig. 7A) that provides three binding interfaces, similar to those described for DiMop, at which six WO<UP><SUB>4</SUB><SUP>2−</SUP></UP> are bound. The highly symmetric assembly creates two additional binding sites located on a 3-fold axis and formed by residues from three neighboring Mop domains (A/C/E and B/D/F) and therefore do not occur in ModE (Fig. 7).


View larger version (37K):
[in this window]
[in a new window]
 
Fig. 7.   Schematic of the modular arrangements for four classes of Mop proteins. A, the hexameric organization of the single Mops as exemplified by S. ovata Mop. Each Mop domain is indicated by a large sphere and the molybdate ions by smaller cyan-colored spheres. The oxyanions are positioned at the Mop-Mop interfacial regions. Only those sites at the dimer interface that bind 6 oxyanions are shown; the other sites are buried at the trimer interface behind subunit A in one case and subunit D in the other. B, the likely trimeric arrangement of polypeptides in E. coli ModG. The only difference with respect to S. ovata Mop is that the Mop domains are linked covalently as indicated by the cylinder joining the Mop domains. C, the dimeric DiMop arrangement in ModE; and D, a postulated organization for the Mop domain of ModC.

A Mop-like sequence is also observed as a C-terminal extension in the ATP-binding subunit of the mod-encoded molybdate transporter (e.g. ModC from E. coli or ModD from A. vinelandii). Again the function of the Mop domain is unclear, but by analogy with similar extensions, e.g. of the maltose transporter ATP-binding protein MalK (34), we suggest that it could allow regulation of transporter activity depending on intracellular molybdate levels. Such rapid feedback regulation would complement the transcription level control exercised by ModE.

Fig. 5 shows a sequence alignment of Mop domains from a representative selection of Mop-carrying proteins. With the exception of some of the Mop2s from ModE proteins (see below), all sequences show conservation of the residues involved in molybdate binding. In a number of cases Ser166 (numbering for the Mop1 domain of E. coli ModE) is replaced with an alanine, but from the structural data it seems perfectly conceivable that in those cases a water could take the place of the serine hydroxyl group.

Based on the known structures, we propose that Mop-containing proteins can be divided into at least four groups depending on the quaternary arrangement of their Mop moieties. The first group comprises the single-Mop molbindins. These Mops dimerize, and the dimers then trimerize to form a hexameric assembly as shown in Fig. 7A. Each pair of dimers provides a molybdate binding interface like that observed in ModE, altogether supplying six binding sites. Two additional MoO<UP><SUB>4</SUB><SUP>2−</SUP></UP> binding sites are located on the 3-fold axis. Thus, the hexameric Mop proteins bind eight molecules of molybdate, all of them at subunit interfaces. DiMop molbindins with two Mop domains on a single polypeptide constitute the second structural group. Even though no structural data are currently available, it can be expected that DiMop proteins will trimerize to form an assembly similar to the Mop hexamer, the main difference being that in this case the participating dimers are covalently linked (Fig. 7B). Likewise the DiMop molbindins should also posses eight molybdate binding sites.

E. coli ModE is a member of the third group. Here only two DiMop chains dimerize to form an assembly that can be superimposed almost perfectly onto two "`dimers" from the S. ovata Mop structure (Fig. 7C). As the "missing" dimer participated in the formation of not only the two "3-fold" sites but also two of the binding interfaces, it is clear that ModE-like proteins can bind only 2 mol of molybdate/mol of (dimeric) protein. It is furthermore clear that two of the four Mop domains (one from each dimer, dark green and purple in Fig. 7C) are no longer involved in oxyanion binding. Accordingly it is not surprising that in some ModE-like proteins (e.g. the ones from E. coli or H. influenzae, Fig. 5) the Mop2 domain has undergone amino acid changes that render molybdate binding impossible. Consider the Mop2 domain of E. coli ModE; the residue corresponding to Mop1 Ser126 is Asn198, which, in the observed conformation, would for steric reasons be unable to bind molybdate. Arg128 is substituted by Asp200, its side chain facing toward the hypothetical binding site, preventing the approach of an oxyanion. Asn252 from the Mop2 domain is the equivalent residue to Lys183 of Mop1, and again a residue important for molybdate binding is lost. The other two binding site residues, Thr163 and Ser166, are replaced by Pro234 and Glu237, respectively, which abolishes all hydrogen bonding and introduces another negative charge incompatible with oxyanion binding.

We note that Mop proteins with a proposed storage function maximize the number of molybdates bound per Mop unit (=1.33), whereas the oligomerization state of other Mop proteins with a function in signaling/feedback provide only two binding sites, thus ensuring that no molybdate is "wasted" (e.g. E. coli ModE binds 0.5 molybdates/Mop).

The fourth group of Mop proteins are the ATP-binding components of the molybdate-transporting ABC transporter system, for example E. coli ModC and A. vinelandii ModD. On the basis of sequence alignments (Fig. 5), these proteins carry a Mop domain that should be able to bind molybdate. To do so it must at least dimerize with a Mop domain-containing protein, and it is possible that this involves another subunit of ModC itself.

Conclusions-- Our analysis shows that the molybdate binding properties of the C-terminal DiMop fragment of E. coli ModE are the same as for the intact protein; it has provided details of oxyanion binding together with the conformational and quaternary changes that result. Structural data have revealed that E. coli ModE discriminates between oxyanions based on size and charge. This molecular recognition is vital for the biological role of ModE in sensing and then controlling the internal molybdate concentration in E. coli by regulation of the transcription of the modABCD operon. The DiMop domain undergoes a ligand-induced conformational change that produces an alteration to the surface of the DiMop dimer and by implication to the surface of ModE. We previously suggested a role for the C-terminal domain of ModE in the recruitment of partner proteins to form the complexes necessary for the regulation of transcription (16). The conformational alteration that the C-terminal domain undergoes when oxyanion binds may represent a molecular switch that regulates the recruitment of such partner proteins necessary for the positive regulation of transcription. The possibility exists that this conformational switch in the C-terminal domain may extend to the N-terminal domain and influence the interactions with DNA, although further work will be required to investigate this hypothesis.

    ACKNOWLEDGEMENTS

We thank C. Bond, D. Hall, and J. Harrison for their help and useful discussions.

    FOOTNOTES

* This work was funded by the Biotechnology and Biological Sciences Research Council and The Wellcome Trust.The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

The atomic coordinates and the structure factors (code 1H9R/1H9S) have been deposited in the Protein Data Bank, Research Collaboratory for Structural Bioinformatics, Rutgers University, New Brunswick, NJ (http://www.rcsb.org/).

§ These authors made equal contributions to this paper.

|| Present address: IBLS, Div. of Biochemistry and Molecular Biology, Joseph Black Bldg., University of Glasgow, Glasgow G12 8QQ, UK.

** To whom correspondence should be addressed. Tel.: 1382-345745; Fax: 1382-345764; E-mail: w.n.hunter@dundee.ac.uk.

Published, JBC Papers in Press, March 20, 2001, DOI 10.1074/jbc.M100919200

    ABBREVIATIONS

The abbreviations used are: SPB, sulfate-binding protein of Salmonella typhimurium; r.m.s., root-mean-square; Mop, molybdate-binding protein.

    REFERENCES
TOP
ABSTRACT
INTRODUCTION
MATERIALS AND METHODS
RESULTS AND DISCUSSION
REFERENCES

1. Kisker, C., Schindelin, H., and Rees, D. C. (1997) Annu. Rev. Biochem. 66, 233-267[CrossRef][Medline] [Order article via Infotrieve]
2. Rajagopalan, K. V. (1996) in Escherichia coli and Salmonella: Cellular and Molecular Biology (Neidhardt, F. C., ed), 2nd Ed. , pp. 674-679, American Society of Microbiology Press, Washington, D. C.
3. Pau, R. N., Klipp, W., and Leimkuhler, S. (1998) in Iron and Related Transition Metals in Microbial Metabolism (Winkelmann, G. , and Carrano, C. J., eds) , pp. 217-234, Harwood Academic Publishers, London, UK
4. Hu, Y., Rech, S., Gunsalus, R. P., and Rees, D. C. (1997) Nat. Struct. Biol. 4, 703-707[Medline] [Order article via Infotrieve]
5. Lawson, D. M., Williams, C. E., Mitchenall, L. A., and Pau, R. N. (1999) Structure 6, 1529-3939
6. Pflugrath, J. W., and Quiocho, F. A. (1988) J. Mol. Biol. 200, 163-180[Medline] [Order article via Infotrieve]
7. Rech, S., Deppenmeier, V., and Gunsalus, R. P. (1995) J. Bacteriol. 77, 1023-1029
8. Walkenhorst, H. M., Hemschemeier, S. K., and Eichenlaub, R. (1995) Microbiol. Res. 150, 347-361[Medline] [Order article via Infotrieve]
9. Grunden, A. M., Self, W. T, Villain, M., Blalock, J. E., and Shanmugam, K. T. (1999) J. Biol. Chem. 274, 24308-24315[Abstract/Free Full Text]
10. McNicholas, P. M., Rech, S. A., and Gunsalus, R. P. (1997) Mol. Microbiol. 23, 515-524[Medline] [Order article via Infotrieve]
11. Anderson, L. A., Palmer, T., Price, N. C., Bornemann, S., Boxer, D. H., and Pau, R. N. (1997) Eur. J. Biochem. 246, 119-126[Abstract]
12. McNicholas, P. M., Mazzotta, M. M., Rech, S. A., and Gunsalus, R. P. (1998) J. Bacteriol. 180, 4638-4643[Abstract/Free Full Text]
13. McNicholas, P. M., Chiang, R. C., and Gunsalus, R. P. (1998) Mol. Microbiol. 27, 197-208[CrossRef][Medline] [Order article via Infotrieve]
14. Self, W. T., Grunden, A. M., Hasona, A., and Shanmugam, K. T. (1999) Microbiology 145, 41-55[Abstract]
15. Anderson, L. A., McNairn, E., Leubke, T., Pau, R. N., and Boxer, D. H. (2000) J. Bacteriol. 182, 7035-7043[Abstract/Free Full Text]
16. Hall, D. R., Gourley, D. G., Leonard, G. A., Duke, E. M. H., Anderson, L. A., Boxer, D. H., and Hunter, W. N. (1999) EMBO J. 18, 1435-1446[Abstract/Free Full Text]
17. Murzin, A. G. (1993) EMBO J. 12, 861-867[Abstract]
18. Hinton, S. M., and Freyer, G. (1986) Nucleic Acids Res. 14, 9371-9380[Abstract]
19. Lawson, D. M., Williams, C. E., White, D. J., Choay, A. P., Mitchenall, L. A., and Pau, R. N. (1997) J. Chem. Soc. Dalton Trans. 21, 3981-3984[CrossRef]
20. Wagner, U. G., Stupperich, E., and Kratky, C. (2000) Structure 8, 1127-1136[CrossRef][Medline] [Order article via Infotrieve]
21. Grunden, A. M., and Shanmugam, K. T. (1997) Arch. Microbiol. 168, 345-354[CrossRef][Medline] [Order article via Infotrieve]
22. Gill, S. C., and von Hippel, P. H. (1989) Anal. Bioch. 182, 319-326[Medline] [Order article via Infotrieve]
23. Otwinowski, Z., and Minor, W. (1997) Methods Enzymol. 276, 307-326
24. Collaborative Computational Project Number 4. (1994) Acta Crystallogr. Sect. D Biol. Crystallogr. 50, 760-763[CrossRef][Medline] [Order article via Infotrieve]
25. Navaza, J., and Saludjian, P. (1997) Methods Enzymol. 276, 581-594
26. Perrakis, A., Morris, R. J., and Lamzin, V. S. (1999) Nat. Struct. Biol. 6, 458-463[CrossRef][Medline] [Order article via Infotrieve]
27. Jones, T. A., Zou, J. Y., Cowan, S. W., and Kjeldgaard, M. (1991) Acta Crystallogr. Sect. A 47, 110-119[CrossRef][Medline] [Order article via Infotrieve]
28. Murshudov, G. N., Vagin, A. A., and Dodson, E. J. (1997) Acta Crystallogr. Sect. D Biol. Crystallogr. 53, 240-255[CrossRef][Medline] [Order article via Infotrieve]
29. Lamzin, V. S., and Wilson, K. S. (1993) Acta Crystallogr. Sect. D Biol. Crystallogr. 49, 129-147[CrossRef][Medline] [Order article via Infotrieve]
30. Barton, G. J. (1993) Protein Eng. 6, 37-40[Medline] [Order article via Infotrieve]
31. Kraulis, P. J. (1991) J. Appl. Crystallogr. 24, 946-950[CrossRef]
32. Merritt, E. A., and Murphy, M. E. P (1994) Acta Crystallogr. Sect. D Biol. Crystallogr. 50, 869-873[CrossRef][Medline] [Order article via Infotrieve]
33. Kleywegt, G. J., and Jones, T. A. (1994) Acta Crystallogr. Sect. D Biol. Crystallogr. 50, 178-185[CrossRef]
34. Dean, D. A., Reizer, J., Nikaido, H., and Saier, M. H., Jr. (1990) J. Biol. Chem. 265, 21005-21010[Abstract/Free Full Text]
35. Kleywegt, G. J., and Jones, T. A. (1995) Structure 3, 535-540[Medline] [Order article via Infotrieve]


Copyright © 2001 by The American Society for Biochemistry and Molecular Biology, Inc.