![]()
|
|
||||||||
J. Biol. Chem., Vol. 280, Issue 6, 4684-4695, February 11, 2005
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||


¶


¶

¶¶
From the
Laboratory of Bioinorganic Chemistry, Department of Agro-Environmental Science and Technology, University of Bologna, I-40127 Bologna, Italy, ||Laboratory of Protein Biochemistry and Protein Engineering, Department of Biochemistry, Physiology, and Microbiology, University of Gent, B-9000 Gent, Belgium, 
Magnetic Resonance Center and Department of Chemistry, University of Firenze, I-50019 Firenze, Italy, and the 
Department of Biochemistry and Molecular Biology, the Pennsylvania State University, University Park, Pennsylvania 16802
Received for publication, July 27, 2004 , and in revised form, November 9, 2004.
| ABSTRACT |
|---|
|
|
|---|
-barrel surrounded by
-helices and a P-loop at the N terminus. The protein dynamic behavior observed in solution is critically discussed relative to the structural model, using algorithms for disorder predictions. The results suggest that UreG proteins belong to the class of intrinsically unstructured proteins that need the interaction with cofactors or other protein partners to perform their function. It is also proposed that metal ions such as Zn2+ could have important structural roles in the urease activation process. | INTRODUCTION |
|---|
|
|
|---|
|
3
3
3 quaternary structure with the three different subunits encoded by the ureC, ureB, and ureA genes, respectively, is synthesized in the apo-form that is devoid of nickel. The incorporation of Ni2+ into the active site, leading to the activation of the enzyme, is still poorly understood and is thought to occur in vivo as a stepwise assembly process (6).
The assembly of the active site in vitro can be achieved by using high, nonphysiological concentrations of Ni2+ ions and bicarbonate as the source of CO2, which is needed for the carbamoylation of the active site lysine residue (7). At physiological concentrations of Ni2+ and CO2, the participation of four accessory proteins is required (812). These proteins (UreD, UreF, UreG, and UreE) are encoded by four genes, which are also present in the urease operon together with ureA, ureB, and ureC (12).
Many functional studies have been carried out on the urease accessory proteins from K. aerogenes. KaUreD (
30 kDa) binds to apourease and appears to induce a conformational change required for the next steps of the activation process (13, 14). KaUreF (
25 kDa) binds the KaUreD-apourease complex and seems to facilitate carbamoylation of the nickel-bridging lysine residue and to prevent Ni2+ binding to the noncarbamoylated apourease (15). Cross-linking experiments showed an interaction between KaUreD and the
and
subunits of apourease, whereas KaUreF was shown to interact with the
subunit and to induce a conformational change capable of increasing the accessibility of the nickel ions and CO2 to residues in the active site (16). The interaction of UreD with UreF and with the
-subunit of apourease was suggested by immunoprecipitation and two-hybrid studies carried out on these proteins from Proteus mirabilis (17). In H. pylori a similar experiment indicated the interaction between UreF and UreH, the latter corresponding to UreD in other bacteria (18). KaUreG (
22 kDa) can form a quaternary complex with KaUreDF-apourease, suggesting that such large aggregate is the minimum competent species for the in vitro urease activation (7, 19) and could be required for the process occurring in vivo (20). Finally, KaUreE (a homodimer of
35 kDa) is thought to bind the KaUreDFG-apourease complex, acting as a nickel transporter that delivers Ni2+ to the active site of the enzyme (2125).
Among the four accessory proteins, UreG plays an essential role in coupling cellular metabolism and bioenergetics to the assembly of urease. This protein contains a fully conserved P-loop motif, which is also present in many nucleotide-binding proteins and which is probably related to the in vivo GTP requirement for assembly of the urease active site (26). A direct relationship between UreG and GTP requirement is proven by the evidence that GTP is needed for activation of the Ka-UreDFG-apourease complex, although it has an inhibitory effect on the nickel reconstitution of the apourease and of the KaUreD-apourease and KaUreDF-apourease complexes (26). UreG is also involved in delivering CO2 necessary for the carbamoylation of the nickel-bridging lysine; the curves that correlate the urease activation to different bicarbonate concentrations indicate a higher rate and level of enzymatic activation in the presence of UreDFG-apourease complex (26). In vitro, optimal levels of apourease activation require 0.5 mM GTP. Larger GTP concentrations lead to a decrease of urease activity, probably caused by the chelation of Ni2+ by the nucleotide; this is consistent with the fact that elevated levels of Mg2+ ions partially restore the activation (26). In the presence of UreE, this inhibitory effect is lost, and the urease activation occurs at significantly lower GTP concentrations (27). These observations indicate a correlation between GTP hydrolysis by UreG and Ni2+ transfer to apourease by UreE. Additional evidence of the UreE-UreG interaction in vivo has been obtained by using two-hybrid systems and immunoprecipitation experiments in H. pylori (18). It has been proposed that UreG may induce GTP-dependent structural changes of the apourease, increasing the accessibility of both Ni2+ and CO2 to the developing active site. Alternatively, UreG may use GTP and CO2 to synthesize carboxyphosphate, which could serve as an excellent CO2 donor to the Lys residue (26).
The thorough functional studies described above have paved the way to an understanding of the mechanism of the urease accessory proteins and, in particular, of UreG in Ni2+ trafficking and metabolic regulation. The next echelon in the comprehension of the urease chemistry implies the study of the interaction mechanisms between the accessory proteins and the enzyme at the molecular level, in order to understand how the urease active site is assembled. This goal cannot be achieved without detailed structural information on each accessory protein and its biochemical properties. The only crystal structures available for the urease chaperones are those of UreE from K. aerogenes (28) and B. pasteurii (29). The high degree of similarity between these two structures and between the structures of Ka- and Bp-urease are indications of a conserved molecular mechanism of urease activity, activation, and metal-center building in different species (30). No structural detail is available for the other three chaperones. A study, published in 1997 (20), reported that KaUreG is monomeric in solution and that it does not, by itself, hydrolyze GTP or ATP. Another report in 2003 (31) and concerning HpUreG essentially confirmed these results.
This paper describes a thorough study performed on the recombinant UreG from B. pasteurii. In particular, the cloning, expression, and purification of the protein in its native and His-tagged forms are described, together with evidence confirming the identity of the isolated protein. The oligomerization state and the hydrodynamic properties of the protein in solution were examined by using size exclusion chromatography coupled with light scattering experiments, and the protein folding was checked using circular dichroism, mass spectrometry, and NMR. The metal-binding capability and enzymatic GTPase activity of BpUreG were established for the first time and discussed. Finally, threading (fold recognition) algorithms were applied to calculate a model for the protein structure that is consistent with its solution properties and enzymatic activity. The results represent a significant contribution to the understanding of the role of this metallochaperone in the urease active site assembly.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
Taq polymerase for the PCR and dNTPs were from Display System Biotech. The oligonucleotide primers were synthesized in the Nucleic Acid Facility at the Pennsylvania State University (University Park, PA). The PCR product obtained by using these primers was digested by a combination of NdeI and BamHI restriction enzymes (New England Biolabs), purified by electrophoresis on a 1% (w/v) agarose gel, extracted, and precipitated. By using T4 DNA ligase (Promega), this DNA fragment was ligated at a 2-fold excess of insert to vector, into plasmid pET3a (Novagen), which had been digested with NdeI and BamHI, treated with alkaline phosphatase, and purified by electrophoresis. For the His-tagged BpUreG, the plasmid pET15b was used. Plasmid DNA was isolated from transformants of Escherichia coli strain DH5
(Invitrogen) by the rapid alkaline extraction method, as described (33), digested with appropriate restriction enzymes, and analyzed by agarose gel electrophoresis. The resulting pET3a::ureG and pET15b::ureG plasmids were purified using the StrataPrepTM Plasmid MiniPrep kit (Stratagene). The sequence of the cloned BpUreG gene was confirmed by DNA sequencing. The constructs for both plasmids were inserted by electroporation (Bio-Rad GenePulser II) into the E. coli BL21(DE3) expression host (Novagen) grown in shaking flasks at 37 °C in a medium with the Luria-Bertani (LB) composition (Amersham Biosciences) or on agar (1.5%) plates with the same composition.
BpUreG ExpressionBased on the T7 system (34), large scale expression of BpUreG and His-tagged BpUreG was achieved in 2.5-liter batches of minimum M9 liquid media (1 liter contained6gofNa2HPO4, 3 g of KH2PO4, 0.5 g of NaCl, 1.25 g of (NH4)2SO4, 0.246 g of MgSO4) supplemented with 4 g of glucose per liter of culture. The 15N-enriched proteins were obtained using a medium containing (15NH4)2SO4. Transformed E. coli BL21(DE3) cells were grown at 37 °C (28 °C for the His-tagged protein) with vigorous stirring, until the A600 reached 0.50.8. Expression was induced by addition of isopropyl
-thiogalactopyranoside (IPTG) to a final concentration of 0.5 mM. The cells were harvested 4 h after induction by centrifugation at 8,000 x g for 10 min at 4 °C. The cells were resuspended in 25 ml of 50 mM Tris-HCl buffer, pH 8, containing 5 mM EDTA and lysozyme (200 µg/ml). After incubation at 30 °C for 20 min, followed by the addition of DNase I (20 µg ml-1) and additional incubation at 37 °C for 20 min, the cells were disrupted by two passages through a French pressure cell (SLM, Aminco) at 20,000 pounds/square inch. The cell pellet was separated from the supernatant by centrifugation at 15,000 x g for 15 min at 4 °C.
BpUreG PurificationIn the case of native BpUreG, the pellet was washed (resuspended in 25 ml of buffer using a mixer homogenizer and centrifuged at 15,000 x g for 15 min at 4 °C) three times with 50 mM Tris-HCl buffer, pH 8, containing 5 mM EDTA, 1 mM DTT, and 2% (w/v) Triton X-100 and three times with the same buffer without Triton X-100. The pellet was subsequently resuspended and incubated overnight at 4 °C in 50 mM Tris-HCl buffer, pH 8, containing 1 mM DTT and 2 M urea. The soluble fraction, obtained after removal of the precipitated material by centrifugation (15,000 x g, 15 min), was loaded onto a Q-Sepharose XK 26/10 column (Amersham Biosciences) that had been pre-equilibrated with 2 volumes of 50 mM Tris-HCl buffer, pH 8, containing 1 mM DTT and 2 M urea. The column was washed using a flow rate of 3 ml min-1 with the starting buffer until the base line was stable. The protein was eluted from the column with a 400-ml linear gradient of NaCl (01 M). Fractions containing BpUreG were combined, diluted with the elution buffer to a protein concentration of 0.3 mg/ml, and dialyzed (5-kDa cut-off membrane) overnight at 4 °C against 50 mM Tris-HCl buffer, pH 8. The resulting solution of BpUreG was concentrated by using 5-kDa cut-off membrane Amicon and Centricon ultrafiltration units (Millipore), to a final volume of 5 ml, and centrifuged (15 min at 14,000 x g) to remove the precipitated material. The resulting solution was loaded onto a Superdex 75 XK 26/60 column conditioned with 50 mM Tris-HCl buffer, pH 8, containing 0.15 M NaCl and 1 mM DTT. BpUreG was eluted at a flow rate of 2 ml min-1, and the purified protein, amounting to
50 mg per liter of culture, was concentrated to 2.5 mg ml-1 and stored at -80 °C.
In the case of His-tagged BpUreG, the supernatant after pellet separation was loaded onto a column containing 8 ml of the nickel-nitrilotriacetic acid Superflow affinity resin (Qiagen) pre-equilibrated with 40 ml of 50 mM Tris-HCl buffer, pH 8.0, containing 5 mM imidazole, washed with 30 ml of the same buffer containing 20 mM imidazole, and eluted using the same buffer containing 100 mM imidazole.
Protein purity, as well as the molecular mass of BpUreG in denaturing conditions, was estimated by SDS-PAGE according to the method of Laemmli (35), by using a Bio-Rad Mini-Protean II apparatus. Proteins were separated on 15% (w/v) acrylamide-bisacrylamide separating gels that were stained using either Coomassie Brilliant Blue R-250 or silver staining.
Protein concentration was measured by using a Jasco 7800 spectro-photometer and a value for the extinction coefficient (
280 = 10,810 M-1 cm-1) calculated from the amino acid sequence using the ProtParam web site (au.expasy.org/tools/protparam.html). This value is in good agreement with that obtained by using the Bio-Rad assay that is based on the Bradford colorimetric method (36).
Protein Sequence DeterminationAll reagents, solvents and instruments were obtained from Applied Biosystems. N-terminal sequence analysis was performed in the gas-pulsed liquid phase using a model 476A protein sequencer with a micro-reaction chamber and an on-line high pressure liquid chromatography system for phenylthiohydantoin analysis. Absorbance was monitored at 269 nm. C-terminal sequence analysis was performed on a Procise 494C protein sequencer using C-terminal sequencing chemistry (37, 38). Prior to this analysis, the sample was adsorbed on a Prosorb sample preparation cartridge and, after subsequent washes with MilliQ-filtered water, treated with 200 mM phenylisocyanate in acetonitrile under basic conditions (124 mM diisopropylethylamine/acetonitrile) in order to modify the
-amino group of the lysine residues into stable phenylureas. The alkylated thiohydantoin (ATH)-amino acids were analyzed on-line using a thermostated (38 °C) C18 reverse-phase column (2.1 x 220 mm, 5 µM). A linear gradient with a flow rate of 300 µl min-1 was formed using a 140C microgradient system with the following solvents: solvent A, 35 mM sodium acetate buffer, 3.5% tetrahydrofuran/MQ water, and solvent B, 100% acetonitrile. The ATH-amino acid derivatives were monitored using a 785A absorbance detector set at 254 nm, and quantified relative to a 100 pmol of ATH-amino acid standard. The methyl naphthylthiohydantoin amino acid standards were obtained from the supplier.
Mass SpectrometryAll mass spectrometric analyses were performed on a quadrupole TOF mass spectrometer (Micromass), interfaced to a chip-based nano-ESI source (NanoMate100, Advion Biosciences). The mass spectra were processed using MassLynx version 3.1 software of Micromass. Before analysis, the buffer was changed to 50 mM NH4 acetate, pH 6.5. The denatured protein (1 µM) was measured in 50% acetonitrile, 0.1% formic acid. For native measurements, the protein (5 µM) was kept in 50 mM NH4 acetate, pH 6.5. Acquisition conditions were a spraying voltage of 1.5 (denatured protein) or 1.7 kV (native protein), gas pressure of 0.3 pounds/square inch, and an acquisition time of 3 (denatured protein) or 10 min (native protein) across an m/z range of 5003000.
Hydrodynamic PropertiesThe molecular mass and hydrodynamic radius of the native protein were estimated by standard size exclusion chromatography. A small amount (100 µl, 2.5 mg/ml) of the purified protein solution was applied to a Superdex-75 HR 10/30 FPLC column that had been equilibrated with 50 mM Tris-HCl, pH 8.0, containing 0.15 M NaCl, at a flow rate of 0.5 ml min-1 in order to estimate the apparent hydrodynamic volume of the purified protein. The column was calibrated using an Amersham Biosciences low molecular weight gel filtration calibration kit.
Absolute estimates of molecular mass and hydrodynamic radius of BpUreG were determined using a combination of size exclusion chromatography, multiple angle light scattering (MALS), and quasi-elastic light scattering (QELS). BpUreG (100 µl, 2.5 mg/ml) in 20 mM Tris-HCl, pH 7.5, 150 mM NaCl was loaded onto an S-200 16/60 column (Amersham Biosciences), pre-equilibrated with the same buffer, and eluted at room temperature at a flow rate of 1 ml/min. The column was connected downstream to a multiangle laser light (690.0 nm) scattering DAWN EOS photometer (Wyatt Technology). Quasi-elastic (dynamic) light scattering data were collected at a 90° angle by using a Wyatt-QELS device. The concentration of the eluted protein was determined using a refractive index detector (Optilab DSP, Wyatt). Values of 0.185 for the refractive index increment (dn/dc) and 1.330 for the solvent refractive index were used. Molecular weights were determined from a Zimm plot. Data were analyzed using the Astra 4.90.07 software (Wyatt Technology), following the manufacturer's indications.
Circular Dichroism SpectroscopyThe CD spectra of BpUreG and its His-tagged analogue were measured at 20 °C, using a Jasco 710 spectropolarimeter flushed with N2, and a cuvette with 0.01-cm path length. The buffer was 20 mM phosphate, pH 7.5, containing 0.15 M NaCl. The spectra were registered from 190 to 300 nm at 0.2-nm intervals. Ten spectra were accumulated at room temperature and averaged to achieve an appropriate signal-to-noise ratio. The spectrum of the buffer was subtracted. The secondary structure composition of BpUreG was evaluated with the tool available on the Dichroweb server of the Centre for Protein and Membrane Structure and Dynamics, www.cryst.bbk.ac.uk/cdweb/html/home.html (39) using the reference sets 4 and 7.
NMR SpectroscopyNMR spectra of 15N-enriched KaUreG and His-tagged BpUreG were recorded at pH 8.0 and 298 K on a Bruker Avance 800 spectrometer operating at 800.13 MHz. The KaUreG spectrum was recorded using a 5-mm reverse detection probe on a 1-mM sample, whereas the spectrum of His-tagged BpUreG was obtained using a TXI cryoprobe on a 0.45-mM sample. The spectrum of BpUreG isolated from inclusion bodies was recorded on a 0.45-mm 15N-enriched sample at pH 8.0 and 298 K using a Bruker DRX Avance 500 spectrometer operating at 500.13 MHz and equipped with a TXO cryoprobe. 1H,15N-HSQC spectra were acquired using sensitivity improvement (4042) and consisted of 848 scans, spectral windows of 1116 ppm in the proton dimension, and 3040 ppm in the nitrogen dimension, with the carrier set at the water frequency and 118 ppm, respectively. Relaxation delays (including acquisition time) in the range 0.91.2 s were employed. Matrices of 1024 x 256 points or 2048 x 128 points were acquired and transformed into 1024 x 512 or 2048 x 512 points.
Measurement of GTPase ActivityGTP hydrolyzing activity was measured using a colorimetric method. The reaction mixture, containing 20 mM Tris-HCl, pH 8.0, 0.075 M NaCl, 5 mM MgCl2,2mM GTP, and 5 µM BpUreG, in the absence or presence of 25 µM ZnSO4, was incubated at 37 °C. Aliquots (90 µl) were removed at different incubation times and added to 30 µl of a 35% trichloroacetic acid/water solution. Phosphate concentration was determined by the malachite green assay (43).
Metal Binding ExperimentsIn all operations, care was taken to avoid exogenous metal contamination. Ni2+ and Zn2+ nitrate salts solutions were prepared starting from ICP 1000 ppm standard solutions (CPI International) diluted to 1 mM with buffer A, containing NaCl 0.15 M. Equal volumes of BpUreG and metal solutions were mixed in 1:1 ratio to yield a constant concentration of protein (40 and 20 µM BpUreG for nickel and zinc, respectively) and an increasing concentration of metal ion. No precipitation was observed during the titration. The resulting mixtures were incubated for1hat37 °Cand overnight at 4 °C and then filtered by centrifugation using 0.5-ml Centricon filter (cut-off membrane 5 kDa). 400 µl of the filtered solution was diluted to 8 ml with MilliQ water. Metal analysis was performed using a Spectro Ciros CCD ICP optical emission spectrometer (Spectro Analytical Instruments) in combination with a Lichte nebulizer and a peristaltic pump for sample introduction. The ICP emission spectrometer system was calibrated by serial dilutions of appropriate single and multielement standards (CPI International). The standardization curve was made by using standard solutions in the range 0500 µM of nickel and zinc in 2.5 mM Tris-HCl, pH 8, and 7.5 mM NaCl with a linear fitting. An Rf power of 1400 watts, a nebulizer gas flow of 0.8 liter min-1, and a plasma gas flow of 14 liters min-1 were used. The sample uptake was set at 2 ml min-1 for 24 s and a wash time of 15 s at 4 ml min-1 plus 45 s at 2 ml min-1 for each sample. Quality control was established by evaluation of buffer containing standards. In order to estimate the total metal added to every protein sample, 200 µl of every metal solution were mixed in a 1:1 ratio with the blank buffer, diluted to 8 ml with MilliQ water, and measured as the filtered samples. The 221.648 and 231.604 nm lines for nickel and the 202.548, 206.191, and 213.604 nm lines for zinc were used for analysis. The measured metal ion concentrations were corrected with the value obtained for the filtered solution of the protein incubated with the buffer blank, without metal ion. The experimental points were fitted using the MacCurveFit software, and the fit was optimized using a Quasi-Newton algorithm.
Calculation of BpUreG Molecular StructureA similarity search of the protein sequences related to BpUreG was carried out using the program FASTA3 (44, 45) applied to the SwissProt data base. Multiple alignment of all sequences was performed using ClustalW (www.ebi.ac.uk/clustalw) (46). The alignment was optimized using information deriving from secondary structure predictions provided by the program JPRED (www.compbio.dundee.ac.uk/~www-jpred) (47). The PONDR VL-XT algorithm (48, 49) for the prediction of disordered regions of the UreG sequences was accessed through the web site www.pondr.com/(Molecular Kinetics, Inc.).
The sequence of BpUreG was used to search for templates using the 3D-Jury predictor meta-server (50, 51) available at the address bioinfo.pl/Meta/. Only the templates with a 3D-Jury score higher than 80% of the best score were selected. Multiple sequence alignment of the sequence of BpUreG with the selected templates was performed using ClustalW (46). Alignment optimization was carried out comparing information deriving from the secondary structure calculated for the templates using DSSP (52) with BpUreG secondary structure prediction provided by the program JPRED (47).
Model structures were calculated using the program MODELER 6.2 (53) with the model-default options. The program PROSA II (version 3.0, 1994) (54) was used for selecting the best models provided by MODELER and for protein structure analysis to test the coherency and validity of the model structures. The Z score reported in this work is derived through the standard "hide and seek" procedure of the program, by which the score is correlated to the difference in potential energy, calculated using mean field potentials, between the input structure and other randomly assigned folds for its amino acid sequence. A lower Z score corresponds to a more favorable potential energy associated with the structure under examination.
Structure validation was performed using PROCHECK (55) and WHATIF (56). The calculated final structure was deposited in the www.postgenomicnmr.net site. The molecular surface and the electrostatic color coding was generated by the program GRASP (57) using a probe radius of 1.4 Å. The electrostatic potential was calculated using a simple version of a Poisson-Boltzmann solver with the GRASP full charge set. All the histidine residues were considered neutral, and the N- and C-terminal residues were charged. Dielectric constants of 80 and 2 were used for the solvent and protein interior, respectively. The topological diagram was drawn using the program Topology of Protein Structure (TOPS) (58) available at www.tops.leeds.ac.uk.
| RESULTS |
|---|
|
|
|---|
|
Another approach to solve the solubility problem of BpUreG was attempted by cloning and expressing the His-tagged protein using the pET15b plasmid. This method, similar to the pET3a system, produced large amounts of protein in the inclusion bodies but allowed the purification of the protein from the supernatant in a single chromatographic step using a Ni2+-affinity column. The yield of His-tagged BpUreG was significantly lower than in the case of the protein purified from the insoluble pellet (10 mg per liter), but this amount allowed us to compare the fold properties of the protein isolated from the soluble and insoluble fraction, using CD and NMR spectroscopy, as well as the activity assay.
Electrospray ionization, quadrupole TOF mass spectra of BpUreG in denaturing conditions (Fig. 1B) confirmed the high purity of the isolated protein and indicated a mass of 23,084.5 Da, in agreement with the theoretical mass (23,084.1 Da).
As expected from the gene sequence, Edman degradation of the N terminus of BpUreG yielded the sequence MKTIHL. The protein was also subjected to C-terminal amino acid sequence analysis, which provided the expected sequence ESK. In the case of the latter experiment, the sensitivity dropped rather rapidly, yielding only the last three residues, most likely because of the presence of consecutive "problematic amino acids" (59). However, this information, together with the mass of the recombinant protein, is sufficient to demonstrate that the recombinant-purified BpUreG protein is intact and unmodified.
Hydrodynamic Properties and Oligomerization of BpUreG Gel filtration experiments performed on BpUreG revealed the presence of two peaks corresponding, according to their retention volumes, to hydrodynamic radii (Rh) of 2.8 (minor peak) and 3.5 nm (predominant form), and apparent molecular mass of 39 and 59 kDa. The gel filtration profile for the soluble His-tagged BpUreG analogue was very similar to that of BpUreG isolated from inclusion bodies, with a smaller amount of the minor peak. These masses are not easily related to multimeric forms of the BpUreG monomer (23,084 Da). These data, based on the calibration of the column with globular standard proteins, assuming a globular structure for BpUreG, are strongly dependent on the properties of the protein and the column, as well as on the possible interactions between the protein and the solid phase. Therefore, a better estimate of the molecular masses and hydrodynamic radii of the protein eluted in the two different peaks was obtained by using a combination of SEC, MALS, and QELS (60, 61). The advantage of such a system is that the molecular mass determination is based on the fundamental light scattering properties of the macromolecule. The results obtained are shown in Fig. 2. The elution profile and the MALS data are consistent with the presence of a predominant dimeric form of the protein in solution (
95%, mass = 40.0 ± 0.5 kDa), whereas the monomeric BpUreG represents only a minor fraction (
5%, mass = 28.6 ± 4.9 kDa). The MALS data confirm the anomalous chromatographic behavior of the protein and suggest either a nonspherical (prolate or oblate) shape of the protein or the presence of nonspecific interactions between the protein and the column phase. Because the scattering of the applied light (wavelength 690 nm) by such small molecules is isotropic, no information on the protein size can be obtained from the static MALS. However, such information can be derived from dynamic light scattering (QELS) experiments performed simultaneously with the MALS measurements. The hydrodynamic radius (Rh) of the dimeric BpUreG, as determined by QELS, is 2.00 ± 0.02 nm, whereas the value of Rh for the monomer is larger (2.70 ± 0.20 nm). The apparent inconsistency between the mass and the volume of the monomer and dimer of BpUreG can be resolved by considering that the monomer is largely unfolded. Considering the larger value of Rh for the monomer, the retention volume for the monomeric form in the size exclusion column is expected to be smaller than that for the dimeric protein; the aberrant behavior of the monomer in the gel filtration experiment could be explained with unspecific interactions occurring between the unfolded monomer and the solid phase of the column. The mass spectrum of BpUreG under nondenaturing conditions revealed the presence of the monomer only, an indication that the dimerization does not involve the formation of covalent bonds (not shown). In the ESI-MS experiments, the injected sample must be maintained at low ionic strength, a condition that could destabilize the dimer, especially if hydrophobic forces are the main responsibility for the dimerization of the protein. Moreover, hydrophobic interactions are weak in gas phase, further causing the dimer to dissociate in the TOF stage. The mass spectrum also confirms that the monomer is largely unfolded, with the presence of populations featuring charge states (from 9+ to 29+) reflecting high solvent accessibility of protonation sites.
|
208 nm and a pronounced shoulder at
220 nm, as well as a maximum positive deflection at
193 nm, typical for the presence of both
-helix and
-strand regions in the protein. The spectrum was quantitatively analyzed in the range 190240 nm, on a per amide basis calculated from protein content and sequence, using all the different fitting programs at the Dichroweb server (39) and all the possible reference sets. The best fit was selected on the basis of the normalized root mean square deviation (0.038) between the experimental and calculated data, obtained by using the variable selection method program CDSSTR (62) and the reference set n = 4. From this analysis, a secondary structure composition of 15%
-helix, 29%
-strand, 26% turns, and 30% random coil was estimated for BpUreG. The CD spectrum of His-tagged BpUreG confirmed these data, yielding 18%
-helix, 25%
-strand, 29% turns, and 30% random coil.
|
|
|
![]() | (Eq. 1) |
|
-strands. Two additional helices (H4 and H5) and one strand (S4) are found in some sequences but are not fully conserved. These structural elements are separated by large portions of turns or coils of variable length. The regions predicted to be in a strand conformation are mainly characterized by the presence of hydrophobic residues, suggesting that they constitute the hydrophobic core of the protein. On the other hand, the helices are mostly amphipathic, suggesting their involvement in intermediate regions between the hydrophobic core and the solvent-exposed surface of the protein. The similarity of the secondary structure elements is paralleled by a high degree of sequence identity among all UreG proteins (between 49 and 62%). Large portions of the sequences are fully conserved, and these regions mostly occur in the loops and coils rather than in the helices and strands. This suggests that these less structured regions are functionally most important, whereas the helices and strands are only needed to confer the necessary overall structure to the protein. The P-loop, needed for GTP binding, is found between strand S1 and helix H1 in the N-terminal region of the protein and consists of a conserved GPVGXGKT motif, where X is usually Ser or, rarely, Ala.
|
In order to calculate a model structure for the fully folded BpUreG, the first choice involved the use of protein templates having a known structure and high homology with the target protein (homology building). However, a data base search for homologous proteins of known structure using FASTA3 (44, 45) and PSI-BLAST (64, 65) did not result in any hit with a sequence identity higher than 30%, preventing the use of standard homology modeling protocols for structure prediction (66). Therefore, algorithms based on fold recognition were attempted. The 3D-Jury approach was singled out for its performance in structure prediction in the most recent CASP5 experiment (67), and the use of this protocol resulted in five template structures with a score higher than 81.3 (i.e. higher than 80% of the best), as reported in Table I. The PDB 1FFH [PDB] structure has been refined (68) to a higher resolution (1.10 Å) and the PDB code 1LS1 [PDB] was used in its place. All the found templates fall into the c.37.1 SCOP class (69), indicating that BpUreG is similar to proteins that contain a P-loop and are involved in GTP-based metabolic processes.
|
|
The fold of the BpUreG monomer is characterized by a central open
-barrel, formed by seven parallel and two anti-parallel strands, surrounded by six
-helices connected with loops (Fig. 8). The 36% of the modeled residues of BpUreG are involved in
-helices, and the 22% are part of the central
-strands, whereas the remaining 42% comprises turns and coils. The P-loop is located between strand S1 and helix H1 and is found on one side of a deep pocket defined, on the other side, by the loop located between strand S5 and helix H4. Fig. 8 also report the positions of the so-called Switch I (between strand S2 and helix H2) and Switch II (between strand S5 and helix H4) regions, found to be important for the binding of the GTP-Mg2+ adduct in G-proteins (70). The binding of Mg2+ is generally needed for GTPase function, as it is involved in the GTP phosphate chain binding to the protein at typical consensus sequence ([DE]XXG) (70). The topological diagram for the model of BpUreG is also shown in Fig. 8. The C-terminal portion of BpUreG, comprising residues 189211, was not modeled due to the absence of a structural template for this region. This peptide sequence is predicted to exist in a
-helical conformation (H7 in Fig. 7), whose topology and orientation with respect to the rest of the protein could not be determined. The presence in the model of strand S9 at the C terminus, in contrast to the predicted presence of an helix (H7), may be explained by an artifact because of a termination effect of the modeling process.
In order to gain more information on BpUreG model structure characteristics, the solid surface representation of the electrostatic potential was calculated (Fig. 9). The molecular surface of the model structure is characterized by large patches of negative charge, with the exception of a large neutral zone situated in the region of helixes H4, H5, and H6. The large pocket observed in the C-terminal region is probably due to the absence, in the model, of the C-terminal helix, which comprises the last 24 residues. This suggests that such helix could fill this pocket in a folded form of the protein.
|
| DISCUSSION |
|---|
|
|
|---|
By using a combination of SEC and MALS (60, 61), BpUreG was shown to exist in solution as a dimer. In the past, hydrodynamic studies on KaUreG reported the presence of only the monomer in solution (20), whereas for HpUreG the molecular mass under native conditions was not determined (31). In the case of KaUreG only chromatographic experiments were used to establish its molecular aggregation form in solution. This approach is partially uncorrected and subject to errors, because it is strictly dependent on the shape of the protein and on its interaction properties with the solid chromatographic phase. This possibility was proven for BpUreG, shown to produce an aberrant chromatographic profile that could have lead to wrong conclusions were it not for the use of light scattering techniques, which yielded an incontrovertible value for the molecular mass of the protein in solution. The aggregation state of BpUreG is similar to that found for H. pylori HypB, a GTPase involved in the activation of Ni2+-dependent hydrogenase and reported to exist as a mixture of the monomeric and dimeric forms (76). The hydrodynamic radii measured by dynamic light scattering (QELS) for the dimer and the monomer BpUreG (2.00 ± 0.02 and 2.70 ± 0.20 nm, respectively) suggest that the latter is largely unfolded, whereas the dimer is present with a compact behavior, therefore representing the actual functional form of the protein.
In solution, BpUreG showed a well defined secondary structure, with both
-helices (15%) and
-strand regions (29%), as determined using CD spectroscopy (Fig. 3). These data do not contain any information regarding the tertiary structure of BpUreG, which has therefore been evaluated using NMR spectroscopy. 1H,15N-HSQC NMR spectra revealed that BpUreG does not possess a rigid tertiary structure but exists in solution in fast equilibrium among different conformations (Fig. 4) and therefore contains large portions of unfolded backbone. The intrinsically unfolded state of BpUreG in solution is not an artifact due to the purification method, which involved the use of a small amount of urea to solubilize the protein from inclusion bodies. This is proven by the essentially identical spectroscopic and hydrodynamic properties of the His-tagged form of BpUreG, purified in a single step from the soluble cellular extract.
The large similarity of the NMR spectral properties of Ka-UreG and BpUreG suggests that the presence of large portions of the protein backbone undergoing conformational changes, and therefore causing a substantial intrinsic protein unfolding, is a general feature for all UreG chaperones. This evidence indicates that UreG belongs to the ever-growing class of intrinsically unstructured proteins. These are natively unfolded polypeptides that undergo disorder-order transitions among the random coil, pre-molten globule, molten globule, or fully folded states during or prior to their biological function (63, 77, 78). The observed conformational plasticity of UreG proteins is consistent with their role as chaperones with GTPase activity, assuming a fully active conformation only in the presence of a preformed complex with other urease accessory proteins, as experimentally proven (20). This characteristic could be related to the functional role of UreG in vivo, because it could permit us to minimize the unwanted hydrolysis of GTP unless the protein is ready to perform its role together with its partner chaperones in the preformed UreDFG-apourease complex. The high homology found in the predicted unstructured parts of the UreG sequences could well be related to this need, indicating the evolutionary functional importance of these loop regions, whereas the
and
elements perform a structural role. Large unstructured portions were predicted to be present in the central part of all UreG sequences found in the data base search (Fig. 7). This characteristic could explain the fluxional behavior detected by NMR spectroscopy for both BpUreG and KaUreG, confirming that this characteristic is a general feature of this class of proteins.
The structure of BpUreG in the fully folded state has been modeled by using a fold recognition procedure. Considering that BpUreG is an intrinsically unstructured protein, we can presume that this model predicts the structure of the protein when other cofactors or protein partners force it to assume its functional conformation. The structural prediction shows a fold typical for a GTPase protein (Figs. 8 and 9) with an identifiable negatively charged P-loop region, likely capable of containing a GTP molecule, and the Switch I and Switch II regions, putatively involved in the Mg2+ binding and generally required for the GTPase function. The overall fold comprises an internal open
-barrel surrounded by
-helices. The hydrophobic and amphipathic composition of the sequences, predicted respectively as
-strands and
-helices in the multiple sequence alignment (Fig. 7), confirms this architecture. The secondary structure composition derived from the model is calculated as 36% for the helix and 22% for
-strands, whereas the remaining 42% is in coils or turns. The partial contrast between this prediction and the calculation of the secondary structure from the CD spectrum (15% of
-helices and 29% of
-strands) can be explained considering that the contradiction involves mostly the more solvent-exposed
-helices, likely more subject to conformational fluctuations, rather than the
-strands, situated in the internal and protected hydrophobic protein core.
The observation of the presence, in the UreG sequences, of a fully conserved P-loop motif, has led us to consider the possible GTPase activity for this protein. A study, published in 1997 (20), reported that KaUreG does not, by itself, hydrolyze GTP or ATP. Indeed, no native or added GTP or ATP was found to be associated with isolated KaUreG, and no interaction of Ka-UreG with ATP- or GTP-linked resins was observed. Similar results were obtained for HpUreG, which showed negligible GTPase activity for the isolated protein (31). The present study demonstrates that BpUreG (and its His-tagged analogue) features a clear GTPase activity, even if low (kcat = 0.04 min-1) compared with other GTPases (79) (Fig. 5). This activity is lower than (but comparable with) the one showed by HypB from B. japonicum (kcat = 0.18 min-1) (80) and from E. coli (kcat = 0.17 min-1) (81). These results indicate that BpUreG, although present in solution as a fluxional, partially unfolded molecule, displays a level of enzymatic activity indicating that a significant fraction of UreG molecules is in the correct fold for catalysis. Alternatively, the fold around the catalytic site could well be correct, and the registered activity intrinsic for this protein, with the unfolded conformation possibly involving a different protein region, as the one indicated by PONDR prediction.
BpUreG has been demonstrated to bind two Zn2+ ions for dimer, whereas the affinity is 10-fold lower for Ni2+ ions (Fig. 6). The analysis of the structural model of BpUreG reveals the presence of a putative metal-binding site, rich in residues commonly found in zinc-binding protein sites and could reasonably be proposed as implicated in the binding of one Zn2+ per monomer. The fully conserved residues likely involved in the metal binding (Glu-64, Cys-68, and His-70) would fall well within the predicted disordered region of UreG (Fig. 7). In this regard, it is interesting to show that this putative metal-binding site is in a region (residues 40-80 in BpUreG numeration) that has been very difficult to model, as described above, probably because of its natively unfolded trait. Therefore, we could speculate that the binding of the Zn2+ ion could induce a conformational change and contribute to the stabilization of the protein backbone in this region. In this case, the metal ion would assume a structural role, stabilizing a certain protein conformation upon binding, rather than a role in the catalytic activity of the protein. Indeed, no change in GTPase activity was observed in the presence of a 5-fold excess Zn2+ ion in the assay mixture.
The specific metal-binding capability of UreG has never been observed before and could be related to its direct role in the assembly of the urease active site (10, 20, 26). HypB, the counterpart of UreG in the [Ni,Fe]-hydrogenase system, has often shown a nickel-sequestering ability, probably due to the presence of an N-terminal His tag, as in the case of Rhizobium leguminosarum (82) and Bradyrhizobium japonicum (80), whereas HpHypB and the HypB from E. coli, lacking the N-terminal His tag, do not bind nickel (76, 81). Our study revealed that BpUreG is able to bind zinc even in the absence of a His tag. In this regard, it is important to consider that UreE, the best characterized accessory protein in this system, is able to bind both Ni2+ and Zn2+ with comparable affinity (