| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||



* Instituto de Biología Molecular y Celular, Universidad Miguel Hernández, Alicante, Spain;
Lehrstuhl für Mikrobiologies, Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen, Germany; and
Biocomputation and Complex Systems Physics Institute, Zaragoza, Spain
Correspondence: Address reprints requests to José L. Neira, Instituto de Biología Molecular y Celular, Edificio Torregaitán, Universidad Miguel Hernández, Avda. del Ferrocarril s/n, 03202, Elche (Alicante), Spain. Tel.: 34-966658459; Fax: 34-966658758; E-mail: jlneira{at}umh.es.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
EI proteins are composed of two domains (6
,7
): a proteolytically inert N-terminal (EIN), and a proteolytically unstable C-terminal domain (EIC). The isolated EIN is a monomer, which cannot be phosphorylated by PEP, but participates in a reversible phosphorylation reaction with phosphorylated HPr (2
,6
). The structure of EIN from Escherichia coli has been solved by x-ray (8
) and NMR (9
), and its complex with HPr has been described by NMR (10
). EIN is composed of an HPr-binding
-helical subdomain and an
/ß subdomain containing the active site histidine. On the other hand, the EIC domain contains the PEP-binding region and the dimerization interface, as shown by the crystal structure of the EIC from Thermoanaerobacter tengcongensis (11
). The folding features of both domains are different. EIN shows a reversible thermal two-state transition; whereas EIC has a broad-scan-rate, concentration-dependent thermal transition (2
,12
). Since the folding process is complex in multidomain proteins, where each domain may be able to refold independently and interdomain interactions may affect the overall folding process, EI can be used, in addition, as a model protein to improve our knowledge of protein folding and assembly of multidomain-oligomeric proteins.
Streptomyces are soil-dwelling actinomycetes which grow on a vast variety of carbon sources, such as mono- and disaccharides, as well as the most abundant biopolymers on Earth (cellulose, chitin, and xylan). These actinomycetes are the origin of approximately two-thirds of all natural antibiotics currently produced by the pharmaceutical industry. The complete genome of Streptomyces coelicolor has been sequenced, and it shows the largest number of genes found in any bacterium (13
). The presence of the different components of the PTS in S. coelicolor has been reported, and the corresponding proteins cloned and expressed (14
16
).
Although there is a wealth of information on the PTS in Gram-negative and low-G+C Gram-positive bacteria, information on PTS proteins present in high-G+C Gram-positive bacteria is very limited. We have undertaken an extensive description of the structures and conformational stabilities of the PTS proteins of S. coelicolor, first, to understand in depth the mechanism of the PTS; and, second, to provide new clues on the protein folding problem. For instance, we have shown that the thermodynamic folding parameters (
Cp and m-values) of the HPr of S. coelicolor, HPrsc, are different from those of the other members of the HPr family (17
19
). Interestingly enough, HPrsc, in contrast to the HPrs in Gram-negative bacteria and low-G+C Gram-positive bacteria (1
), is not involved in general carbon regulation (20
,21
). Then, one is tempted to ask, is there any relationship among the different functions and structures of the PTS proteins on one side, and the thermodynamic parameters and folding properties on the other? To answer this question, it is necessary to describe the conformational preferences of proteins, where studies of different members of the family have been described. There is now increasing evidence that looking at the structure alone is not sufficient to answer these scientific questions (22
). Rather, structures need to be associated to other bodies of data, putting them into specific biological context. Those other sets of data include folding classification, taxonomic identification, their role in metabolic routes, biophysical data, and conformational stability.
Preliminary Fourier transform infrared spectroscopy (FTIR) and circular dichroism (CD) thermal denaturations at physiological pH on the EI enzyme of S. coelicolor, EIsc, indicate that its thermal denaturation behavior is different from that of other EI proteins, since only one transition is observed (23
). In this work, we further investigate the stability and structure of EIsc at different pH in the absence and in the presence of its effectors (PEP and Mg2+), by using theoretical and experimental approaches. Our results show that dimeric EIsc acquired a folded conformation at physiological pH; however, 8-anilino-1-naphthalene-sulfonic acid (ANS)-binding experiments, FTIR, and CD measurements suggested that the protein populated a molten-globule-like species at low pH. Further, we modeled the structure of EIsc, by using the structures of the N-terminal domain of EI from E. coli, EINec (9
), and the C-terminal domain of T. tengcongensis, EICtt (11
), as templates. This modeled structure has allowed us to explain the experimental conformational and thermodynamic features of the protein. Finally, we studied the conformational and stability properties of EIsc in the presence of the effectors (i.e., Mg2+ and PEP). Our results show that the secondary structure of EIsc was not substantially altered, as suggested by FTIR experiments, but quenching experiments and proteolysis patterns suggest that there were small conformational changes upon binding of PEP; in addition, the thermal stability of EIsc was slightly modified when any effector was added. In conclusion, our approach suggests that thermodynamical analysis might reveal subtle conformational changes in otherwise similar folds.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Protein expression and purification
The EIsc clone comprises residues 1576, with an extra methionine and a His6-tag at the N-terminus (15
). The protein was expressed and purified as described (23
). Protein was >95% pure as judged by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE). EIsc was dephosphorylated by overnight dialysis at 278 K against phosphate potassium buffer (10 mM, pH 7.5), 10 mM sodium pyruvate, 1 mM MgCl2, followed by extensive dialysis against the same buffer in the absence of MgCl2 and pyruvate. Samples were frozen in liquid nitrogen and stored at 80°C. Protein concentration was calculated by using the extinction coefficients of Trp, Tyr, Cys, and Phe residues (24
). The yield of the protein was 58 mg per liter of culture.
CD measurements in the far-ultraviolet (UV) region
Spectra were collected on a Jasco J810 (Tokyo, Japan) spectropolarimeter fitted with a thermostated cell holder and interfaced with a Peltier cell, working at 298 K. The instrument was periodically calibrated with (+) 10-camphorsulphonic acid. The final concentrations were 5 mM MgCl2 and 2.5 mM PEP. The mean residue ellipticity, [
], was obtained from the raw ellipticity data,
, as described (23
), and all spectra were corrected with the proper baseline.
In the pH-denaturant experiments followed either by CD or fluorescence (see below), the pH was measured before and after completion of experiments with an ultrathin Aldrich electrode in a Radiometer pH-meter (Copenhagen, Denmark). Three-point calibration of the pH-meter was performed by using standards from Radiometer. The salts and acids used in buffer preparation were: phosphoric acid, pH 2.03.0; formic acid, pH 3.04.0; acetic acid, pH 4.05.5; monosodium dihydrogen phosphate, pH 6.07.0; Tris basic, pH 7.59.0; sodium carbonate, pH 9.511.0; and sodium phosphate, pH 11.513.0.
Steady-state measurements
Spectra were acquired in a 0.1-cm-pathlength cell with a speed of 50 nm/min and a response time of 4 s, and were averaged over four scans. Protein concentration was 4 µM.
Thermal denaturation measurements
Experiments were performed with a response time of 8 s and at constant heating rates of 15, 30, and 60 K/h to test for microreversibility (6
), in the presence of 1 mM ß-mercaptoethanol. No variations in the Tm were observed among the three scan rates, and thus, the criterion for microreversibility was satisfied. Measurements were acquired every 0.2 K. Thermal scans were followed by changes in the ellipticity at 222 nm from 298 to 353 K. The reversibility of thermal transitions was tested by recording a new scan after cooling the thermally denatured sample to 298 K, and comparing the thermal denaturation curve with that obtained in the first scan.
Fluorescence measurements
Fluorescence spectra were collected in a Cary Eclipse spectrofluorometer (Varian, San Carlos, CA) interfaced with a Peltier cell. A 1-cm-pathlength quartz cell (Hellma, Muellheim/Baden, Germany) was used. The slit widths were 5 nm for the excitation and emission wavelengths. Protein concentration was 4 µM, unless otherwise indicated. Steady-state and quenching experiments were acquired at 298 K.
Steady-state intrinsic fluorescence
Protein samples were excited at 280 and 295 nm between pH 212. Experiments were recorded between 300 and 400 nm. The signal was acquired every second, and the increment of wavelength was set to 1 nm. Blank corrections were made in all spectra.
Steady-state ANS binding
ANS binding was detected by collecting fluorescence spectra at different pH in the presence of 100 µM of dye. Excitation wavelength was 370 nm, and the emission intensity was collected from 430 to 700 nm. Stock solutions of ANS were prepared in water by using a molar extinction coefficient of 6.8 x 103 M1 cm1 at 370 nm (25
). Blank corrections were made in all spectra.
Quenching experiments
Quenching of intrinsic tryptophan and tyrosine fluorescence by either iodide or acrylamide was examined at different pH values, by excitation at 280 and 295 nm. In the experiments employing KI as a quencher, ionic strength was kept constant by addition of KCl; Na2S2O3 was added to a final concentration of 0.1 M to avoid formation of I3. The dynamic and static quenching constants were obtained by fitting the fluorescence intensity at 338 nm to the Stern-Volmer equation, which includes an exponential term to account for static quenching (26
):
![]() | (1) |
is the static quenching constant, F0 is the fluorescence in the absence of quencher, and [X] is the concentration of the quencher. The range of quencher concentrations was 00.7 M. Iodide quenching did not show a significant static component, and then the exponential term was not included. In experiments aimed to detect protein self-association, the concentration of EIsc was varied at pH 7 from 0.5 to 4 µM, and only KI was used as the quencher.
In experiments aimed at detecting the burial of the aromatic residues upon binding to PEP and/or Mg2+, the concentration of EIsc was 4 µM at pH 7.0, and only acrylamide was used as quencher. Concentrations of effectors were 2.5 mM of PEP and 5 mM of Mg2+. Experiments were only carried out at 280 nm, where the signal/noise ratio of the Stern-Volmer plots was higher.
Thermal denaturation measurements
Thermal unfolding curves were determined in the presence of 1 mM ß-mercaptoethanol by following the emission fluorescence at 315, 335, or 350 nm, after excitation at 280 and 295 nm, with a scan rate of 30 K/h. The three wavelengths yielded the same unfolding curves (data not shown).
Thermal denaturation experiments were also carried out by following the emission fluorescence of ANS either at 480 or 520 nm, after excitation at 360, 380, or 390 nm, at three different pH, 2.0, 3.5, and 8.0 The scan rates were 30 K/h.
Trypsin digestion experiments
EIsc was mixed with trypsin (at final concentrations of 1 µg/µl for EIsc and 0.002 µg/µl for trypsin, that is, the final concentration rate between both proteins was 500:1) in a total volume of 17 µl in 100 mM of Tris buffer at pH 7 (27
). Lower concentration rates of protein to trypsin were used (100:1, 200:1, and 300:1), but in all these experiments the amount of nondigested EIsc in the absence of effectors was insignificant after 1 min of reaction (data not shown). Samples were incubated at 25°C for 10 min. Digestion was stopped by the addition of 17 µl of SDS-PAGE loading buffer, and samples were immediately run on a gel. The intensities of the bands at different times were measured by densitometry. Experiments were carried out in the absence and in the presence of the effectors (at final concentrations of 2.5 mM PEP and 5 mM Mg2+), and every digestion was repeated three times.
Differential scanning calorimetry
Differential scanning calorimetry (DSC) experiments were performed with a VP-DSC differential scanning calorimeter (MicroCal, Northampton, MA) by using a protein concentration of 1.27 mg/ml (
20 µM). Protein was extensively dialyzed against 2 L of the working buffer (5 mM HEPES buffer containing 10 mM NaCl and 2 mM ß-mercaptoethanol, pH 7) at 278 K. In experiments with the effectors (5 mM Mg2+ and 2.5 mM PEP), the dialyzed protein solution was supplemented with such effectors taken from a stock solution freshly prepared in the same buffer. The sample with the protein and the effectors was incubated for a minimum of 2 h at room temperature. Samples were degassed under vacuum for 1015 min with gentle stirring before being loaded into the calorimetric cell. DSC experiments were performed under a constant external pressure of 30 psi to avoid bubble formation and samples were heated at a constant scan rate of 60 K/h. Once the first scan was completed, the samples were cooled down in situ to 283 K for 15 min and rescanned under the same experimental conditions to test for the reversibility of the denaturation reaction. Calorimetric denaturations were reversible in all cases. Experimental data were corrected from small mismatches between the two cells by subtracting the buffer baseline before data analysis. Once thermograms were normalized to protein concentration, a chemical baseline calculated from the progress of the unfolding transition was subtracted. For experiments involving PEP, control scans were acquired with PEP added to buffer versus buffer in the reference cell, and baselines were found to be the same as in buffer-buffer baselines, suggesting that under these conditions hydrolysis of PEP during heating was not significant. The excess heat capacity functions were then analyzed using the software package Origin 7.0 (Microcal). The best fit was obtained by using the two-state model.
FTIR
The protein was speed-vacuum dried (Speed-vac rotatory evaporator, Savant Instruments, Farmingdale, NY) and dissolved in the deuterated buffer containing 100 mM NaCl. Protein was still functional after the concentration step, as shown by binding experiments (28
). No pH corrections were done for the isotope effects. Samples of EIsc at a final concentration of 44.7 mg/ml (
1 mM) were placed between a pair of CaF2 windows separated by a 50-µm-thick spacer in a Harrick demountable cell. Spectra were acquired on a Bruker (Madison, WI) FTIR instrument equipped with a deuterated triglycine sulfate detector and thermostated with a Braun (Kronberg, Germany) water bath. The cell holder was continuously filled with dry air. When effectors were added, the final concentrations were 5 mM MgCl2 and 2.5 mM PEP at pH 7.
Steady-state measurements
Three-hundred scans per sample were taken, averaged, apodized with a Happ-Genzel function, and Fourier transformed to give a final resolution of 2 cm1. The signal/noise ratio of the spectra was better than 1000:1. Buffer contributions were subtracted, and the resulting spectra were used for analysis. To quantify the secondary structure components, the amide I' band was decomposed into its constituents by curve-fitting (based on a combination of Gaussian and Lorentzian functions). The number and position of bands obtained from the deconvoluted (with a Lorentzian bandwidth of 18 cm1 and a resolution enhancement factor of 2) and the Fourier derivative (with a power of 3 and a breakpoint of 0.3) of the spectra (29
31
) were used for analysis of the percentages of secondary structure. The experiments recorded at 298 K, and acquired at pH 3.0, 7.0, 8.5, and 10.0. Experiments between pH 4 and 5.5, where an acidic titration was observed (see Results), could not be carried out since the acetic moiety used in the buffer absorbs strongly in the amide I' region of the infrared spectrum.
Thermal denaturation measurements
Thermal denaturation experiments were carried out with a scanning rate of 50 K/h, and acquired every 3 K. Fifty scans per temperature were acquired and averaged.
Analysis of the pH- and thermal-denaturation curves
The pH-denaturation experiments were analyzed assuming that both species, protonated and deprotonated, contributed equally to the fluorescence intensity, Y:
![]() | (2) |
Thermal-denaturation parameters were obtained by following the change in the physical property, Y (the fluorescence intensity, the
and the half-height width of the amide I' band of the FTIR experiments), to:
![]() | (3) |
N + ßN[D] and YD =
D + ßD[D] are the baselines of the folded and unfolded states, respectively, for which a linear relationship with denaturant is assumed; R is the gas constant; and T is the temperature in K.
The change in free energy is given by (32
):
![]() | (4) |
Hm is the van't Hoff enthalpy change, Tm is the midpoint of the thermal denaturation,
Cp is the heat capacity change, and Ct is the total molar concentration of the protein expressed in dimer equivalents. The fittings to Eqs. 24 were carried out by using the general curve-fit option of Kaleidagraph (Abelbeck software) working on a PC computer.
Modeling of the EIsc structure
The structures of the N- and C-terminal domains of S. coelicolor were obtained by homology modeling. Sequence alignments between species were made with CLUSTALW at the EBI website (http://www.ebi.ac.uk) by using the Needleman-Wunsch global alignment algorithm (with the following parameters: matrix: EBLOSUM62, gap penalty value: 10.0, extended penalty value: 5.0) to find the optimum alignment (a 40% identity and a 52% similarity).
The models were constructed with the SWISS MODEL server (33
) by using the structures of EINec (PDB accession code 1ZYM) and EICtt (PDB accession code 2BG5) as templates. The orientation and optimization of side chains were carried out in two steps. Firstly, those residues making van der Waals clashes were selected and fitted with "Quick and Dirty" algorithms; and secondly, models were energy minimized (100 steps of steepest descent and 100 conjugate gradients, with a cut-off of 10 Å for nonbonded interactions) with Insight II (Biosym/MSI). Structure edition was made with Swiss PDB viewer v3.7 (34
) and Insight II. The model was tested in terms of energy with Fold-X (http://foldx.embl.de) (35
,36
) which evaluates the properties of the structure, such as the atomic contact map, the accessibility of atoms and residues, the backbone dihedral angles, and the hydrogen bond and electrostatic networks of the protein. Fold-X was also used to evaluate the interaction energy in the dimer formed by both C-terminal domains by calculating the difference in stability energy between the dimer and the isolated domains. Further, the model was evaluated with PROCHECK (37
), showing a Ramachandran plot with 89.1% of the residues in the most favorable regions and 10.3% in additional generously allowed regions. Molecular graphics were created with PyMOL (38
).
| RESULTS |
|---|
|
|
|---|
Then, the experiments monitoring protein stability under different conditions were carried out at protein concentrations
4 µM.
pH-induced structural changes
Steady-state intrinsic fluorescence measurements.
Fluorescence spectroscopy reports the changes in the tertiary structure of a protein (26
). The emission fluorescence spectrum of EIsc, by excitation at either 280 or 295 nm, showed a maximum at 335 nm, and thus, it was dominated by the emission of the sole tryptophan residue, Trp461. The changes in the fluorescence intensity at 335 nm, at both excitation wavelengths, showed a bell-shaped profile, with a maximum at neutral pH (Fig. 1 A, open squares). At acidic pH, a sigmoidal curve with a pKa of 4.3 ± 0.6 and n = 1.4 ± 0.3 was observed. The large fitting errors are probably due to the steepness of the experimental data between pH 5 and 6.
|
Fluorescence quenching.
The solvent-accessibility of tryptophan and tyrosine residues was examined by iodide and acrylamide quenching at several pHs. Three pHs were chosen based on the steady-state fluorescence results: pH 5, close to the midpoint of the titration observed in ANS-binding and fluorescence intrinsic experiments; neutral pH, where the intensity of steady-state fluorescence experiments was at its maximum; and pH 12, in the basic region. Control experiments were carried out at 7 M GdmHCl, pH 7, where the protein was fully unfolded. Experiments at pH <4 led to sample precipitation.
Acrylamide quenching experiments, measured by excitation at either 280 or 295 nm, showed exponential Stern-Volmer plots (Table 1), which yielded similar results: the Ksv were larger at acidic than at neutral pH, and the experiments in GdmHCl had the largest values of Ksv. At basic pHs, the Ksvs were slightly larger than at neutral pH, but smaller than those at acidic pHs. The use of KI yielded similar results at both excitation wavelengths (Table 1): the values of Ksv were larger at acidic pHs and in GdmHCl than at neutral pH. At basic pHs, the Ksvs were similar, within the error, to those observed at pH 7.
|
Far-UV CD experiments.
The far-UV CD spectrum of EIsc at neutral pH showed the features of an
-helical protein with minima at 208 and 222 nm (41
,42
) (Fig. 1 B). This shape did not change as the pH was modified, except in the region between 3.0 to 5.5. In this pH range, the minimum at 208 nm disappeared, whereas that at 222 nm was still present. In the whole pH range explored, the absolute values of the [
] at 222 nm showed a bell-shaped tendency (Fig. 1 C), with a maximum shifted to lower pH (at pH 6) than that found in fluorescent experiments (Fig. 1 A, open squares). This could be due to the fact that CD is reporting global changes in the secondary structure, and fluorescence is monitoring the changes in tertiary structure around fluorescent residues.
FTIR experiments.
At any of the explored pH, the percentages of ß-hairpin (at 1693 cm1) turns/loops/(0,
) ß-sheet (at 1679 cm1) and antiparallel ß-sheet (at 1633 cm1) remained essentially unaltered. Conversely, at acidic pH, the percentage of ß-turns appearing at 1668 cm1 was decreased, and the percentages of helical structure, loops, and disordered structures (bands at 1648 and 1657 cm1, which could not be independently assigned at pH 3) were higher (50.9%) than those at neutral pH (pH 7 (44.5%) and 8.5 (40.3%)), but similar, within the error, to those observed at pH 10 (55%). These data suggest that the percentages of helical structure in EIsc at the extremes of pH were increased.
Stability of EIsc as monitored by DSC and thermal denaturations
From the spectroscopic results above, it can be concluded that EIsc acquired a folded conformation around physiological pH, and that this structure was basically unaltered (except for changes in the ellipticity, Fig. 1 C) between pH 7 and 8. Since, EIsc shows a sigmoidal thermal transition at neutral pH, with a Tm = 330 K (23
), we can wonder whether there are sigmoidal thermal transitions at other pH, and if so, whether the Tms are the same as at pH 7. To address these questions, we carried out calorimetric and thermal measurements by using DSC and several spectroscopic techniques (Table 3).
|
|
Far-UV CD experiments.
Experiments were carried out at pH 2, 3, 4, 7, 9, and 11. In all cases, the absolute value of the ellipticity decreased linearly, as the temperature was raised. Sigmoidal transitions were only observed at pH 7 (23
) and 9, with Tm close to 330 K (Fig. 2 B).
DSC experiments.
Experiments were carried out at pH 4, 7, and 10. At pH 4, the sample precipitated, and at pH 10, no endotherm was observed. The endotherm at pH 7 (Fig. 3 A, lowest line) could be fitted to a two-state mechanism, with Tm = 332 K (Table 3).
|
Structure of EIsc by FTIR (steady-state measurements)
In Fig. 3 B, the amide I' band of EIsc in the presence of both effectors is shown. The deconvoluted bands in the FTIR spectra suggested that, within the experimental uncertainty, the percentages of secondary structure remained unaltered upon binding of both effectors (Table 4). Similar findings were obtained when any of the effectors was added separately.
|
= 0.5 ± 0.1 M1) similar to those measured in isolated EIsc (Table 1), and smaller when PEP (1.6 ± 0.3 M1 and
= 0.5 ± 0.1 M1) or both effectors (1.8 ± 0.2 M1 and
= 0.4 ± 0.2 M1) were added.
Trypsin digestion experiments
The digestion patterns in the presence of Mg2+ showed the same behavior as in the absence of any effector: most of the protein was rapidly digested within the time of the experiments (open and solid squares). Conversely, in the presence of PEP (either alone or with Mg2+), the enzyme was more resistant to proteolysis (open and solid circles) (Fig. 4).
|
The stability of EIsc (thermal and calorimetric measurements)
Thermal denaturation experiments in the presence of the effectors were carried out using far-UV CD and FTIR spectroscopies, and DSC measurements (Figs. 2 C and 3 A) (Table 3). The addition of PEP slightly destabilized EIsc, and the presence of Mg2+ stabilized it. In general, the presence of Mg2+, PEP, or both effectors together had a minor effect on the thermal stability of the protein as judged by the slight changes in Tm (the 
Tms from DSC measurements equal +1.5 K for Mg2+, 1.5 K for PEP, and 0.7 K when both effectors were present).
The EIsc modeled structure
We have modeled the EIsc structure to find a framework in which our experimental results could be adequately explained. Further, this modeled structure will make it possible to explain whether the structures alone can account for the differences in the thermodynamical parameters among the members of the family.
The N-terminal domain showed basically the same structure as the EINec (8
,9
), since the alignments for homology modeling showed only minor deletions in loops between
-helix 1 and
-helix 2,
-helix 2 and
-helix 3, and
-helix 3 and
-helix 4 (Fig. 5 A). In essence, the structure consisted of an
-helix subdomain, which is a four-helix bundle, comprising
-helices 14; an
/ß subdomain, comprising a ß-sandwich, formed by a four-stranded parallel ß-sheet (ß-strands 14) and a three-stranded antiparallel ß-sheet (ß-strands 1, 5, and 6), and three short helices (
-helices 57) that contain the phosphorylation site (His186) (Fig. 5 B). In the modeled structure, there was a long
-helix that is supposed to serve as a linker to the EIC domain.
|
8 barrel fold with an additional helix-loop-helix (formed by
-helix 9-
-helix 10) at its C-terminal. This helix-loop-helix motif folded against the barrel, covering a hydrophobic patch (which includes Trp461); the motif was linked to the barrel by two salt bridges, namely, Arg537-Asp467 and Arg546-Asp459 (Fig. 5 C). The modeled structure of the EICsc monomer differed from that of EICtt (11
-helix 4 and ß-strand 5, which is part of neither the active center nor the dimerization surface.
However, the Fold-X analysis of the modeled EICsc dimer showed differences in the stability and in the interaction energy when compared to those of EICtt (11
). Similarly to that observed for EICec, the electrostatic, as well as the van der Waals interactions, were decreased in EICsc, losing several salt bridges and hydrogen bonds when compared with the structure of EICtt. The helix-loop-helix motif formed by the
-helix 9-
-helix 10 shows intermonomer salt bridges: Asp540 in the motif with Arg422' in
-helix 5 from the other monomer (the prime is used to designate residues in the other monomer). The dimerization surface included several salt-bridge networks, such as Asp467-Arg424', Arg456-Asp358', Asp540-Arg422', Asn345-Asp448', as well as backbone-backbone hydrogen bonds (Fig. 5 C). When compared with the EICtt, the dimeric surface was losing hydrogen bonds (which yielded 
G = + 5.42 kcal mol1) and hydrophobic contacts (
G = + 4.72 kcal mol1), making the monomer-monomer interaction less stable. These differences could be due to the thermophilic nature of EICtt.
The structure of the modeled dimer nicely explained why the solvent accessibility of the Trp461 was reduced as the concentration of EIsc increased (see before). Upon dimer formation, the higher part of the groove, formed by the helix-loop-helix motif and the rest of the protein (and which traps Trp461), was partly fixed by the other monomer via several salt bridges, thus hampering solvent-accessibility to the indole moiety.
We attempted, by docking of the EINsc and EICsc domains, to model the whole structure of EIsc, trying to keep the active-site histidine (in the EINsc domain) close to the PEP binding site (in the EICsc domain). Our models were in all cases (500 runs) unsuccessful, since the distance between both sites was very large to allow for the transference of the phosphate group (>10 Å). This result was also reported in a study of docking of EICtt with EINec (11
), and suggests that the structure of some, if not both, domains should be different in the intact EIsc from in the isolated domains.
| DISCUSSION |
|---|
|
|
|---|
Above pH 8, the intensity of the fluorescence spectra decreased (Fig. 1 A) and the maximum wavelength of the emission spectra moved toward higher values (data not shown), indicating solvent exposure of the fluorescence residues. Furthermore, the Ksv values were higher than those measured at neutral pH but smaller than those at acidic pH and in 7 M GdmHCl (Table 1). However, the molar ellipticity at 222 nm at basic pH was larger (in absolute value) than that at pH 7.0 (Fig. 1 C), thus indicating that the helical content of EIsc increased. This increase in the helicity was further confirmed by the FTIR experiments carried out at pH 10 (Table 2). Taken together, these data suggest that there were structural changes, which partially solvent-exposed the indole moiety and/or tyrosine residues, and modified the ellipticity at 222 nm (41
,42
).
|
] at 222 nm (Fig. 1 C) and the FTIR data acquired at pH 3 (Table 2). Furthermore, the shape of the CD spectrum was altered between pH 3.0 and 5.5 (Fig. 1 B) at wavelengths
208 nm. Changes in the CD signal at 208 nm (where the parallel polarized
-
* transition occurs) are the result of a flexibility of the structure (46
The biological relevance of the partially folded species at acidic pH is unknown, and we do not know whether the residual secondary structure detected by CD and FTIR (Table 2) was native-like. We can speculate, however, based on the modeled three-dimensional structure of the EICsc, that the large amphipathic surface in the dimer interface could be exposed to solvent upon monomer formation; furthermore, there is a network of aspartic residues involved in salt bridges at the dimeric interface (Fig. 5 C). The partial solvent exposure of this hydrophobic surface (due to protonation of some aspartic residue) could then cause binding to ANS, and thus the molten-globule of EIsc could have some of its dimeric interface partially disrupted. However, our experimental observations at low pH also can be explained, for instance, by the solvent exposure of the helix-loop-helix formed by
-helix 9 and
-helix 10, whose movement from the rest of the protein would expose a large hydrophobic surface able to bind ANS, and would alter the fluorescence of Trp461 (Fig. 1 A, open squares). To test this hypothesis, we carried out thermal denaturation experiments followed by the fluorescence of ANS at three different pH (2.0, 3.5, and 8.0) and in all cases, the intensity of the fluorescence at 480 nm was very high at low temperatures (between 288 and 293 K), but decreased as the temperature was raised (data not shown). However, the absence of native baseline at low temperatures did not yield reliable results when fitting the data to Eq. 4.
Folding and stability of EIsc in the absence of effectors
Two unfolding mechanisms have been described in oligomeric proteins exposed to high temperatures or high concentrations of chemical denaturants at equilibrium (49
51
): 1), dissociation followed by unfolding of the native or partially unfolded monomeric species; and 2), denaturation without the accumulation of the native dissociated species. In this work, between pH 6 and 8, EIsc was a well-folded dimer that could be heat-denatured in a single sigmoidal transition, and thus it followed the second mechanism proposed. All the thermal and the calorimetric transitions were reversible and they could be fitted to a two-state model (Table 3).
Converse to observations in other EI proteins by DSC (6
) and CD (52
), the DSC experiments with EIsc did not show two transitions. In the EI protein from Salmonella typhimurium, EIst, the DSC scans show a single peak between pH 6.5 and 7 (6
), due to simultaneous unfolding between the C- and N-terminal domains. This simultaneous unfolding is strongly pH-dependent, and at high pH, two peaks with different thermodynamical properties are observed, one corresponding to EIN (with pH-independent Tm values) and a second corresponding to EIC (with pH-dependent Tm-values) (6
). Interestingly enough, two peaks also appear in the unfolding of EIec (4
). In EIsc, a single peak was observed in CD and DSC experiments, but in the pH range where the thermal transitions could be followed (pH 79), two different peaks were not observed; further, the Tm of the sole transition observed was pH-independent between pH 7 and 9 (Fig. 2 B). As the EINst and EINec have pH-independent Tm values close to 328 K (4
,5
), and thus similar to the Tm value measured in EIsc (Table 2), it is tempting to suggest that the sole peak observed in the thermal scans of EIsc was due to the EINsc. This different behavior in the thermograms among the members of a protein family is not new, and it has been reported, for instance, in the thermodynamic unfolding of dimeric histone-like proteins (53
).
An estimation of the Tm of EICsc, based on the fluorescence results (Fig. 2 A), was also obtained. This value (311 ± 5 K, at 4 µM protein concentration, 315 ± 9 K, at 6 µM protein concentration) was similar, within the experimental uncertainty, to that obtained in the intact EIec (4
), where two tryptophan residues are located at the C-terminal region of the protein. The absence of a long native baseline at this and other pH in EIsc precluded any reliable conclusion on the pH independence of this value of Tm.
We carried out measurements at different protein concentrations either by CD and fluorescence in the concentration range 48 µM at pH 7. The transitions measured by CD became broader as the concentration increased, and although a small increase was observed when the concentration was raised (3 µM yielded Tm = 328.8 ± 0.2 K; 4 µM, Tm = 329.9 ± 0.4 K (Table 1); and 8 µM, Tm = 330.5 ± 0.6 K), all the measurements are within experimental error, and the differences are probably due to the broadness of the transition. The fact that the CD thermal transitions of EIsc became broad as the concentration was increased suggests that the thermal transition of the EICsc is also being observed, as happens in EIec (52
). Either we could not unambiguously determine the protein-concentration independence of the thermal midpoint from fluorescence measurements (see above, for the values, and Fig. 2 A), since the native baselines showed large slopes, probably due to thermoquenching of the tryptophan (26
), leading to large fitting errors. In EIec, the thermal midpoint of the transition monitored by fluorescence is concentration-dependent (52
), but the slopes of the native baseline are smaller due to the fact that two tryptophan residues are observed in the C-terminus of EIec.
Structure and stability of EIsc in the presence of effectors
The addition of the effectors to EIsc did not alter substantially the secondary structure of the protein, as shown by FTIR (Table 4). Similar results have been observed in 1), the enzyme malate synthase G (54
), a close structural analog of pyruvate dikinase and EIC (11
), where binding to pyruvate does not cause large structural changes (55
); and 2), the pyruvate dikinase from maize, where only small changes are observed in selected regions upon binding of PEP (56
). These small changes upon binding of the effectors were also monitored by the changes in the accessibility of fluorescent residues and the trypsin digestion patterns. Upon binding, the accessibilities of the fluorescent residues, when compared to those in EIsc alone, were decreased in the presence of PEP (either alone or with Mg2+), and remained unaltered in the presence of Mg2+ alone. Furthermore, trypsin digestion was reduced in the presence of PEP (either alone or with Mg2+), but not in the presence of Mg2+ alone. Taken together, both results suggest the presence of small conformational changes, which otherwise are not detected by FTIR (Table 4), when PEP was bound to EIsc. These changes led to burial of the Trp and/or Tyr and other regions susceptible to trypsin digestion.
The presence of Mg2+ slightly increased the Tm of EIsc alone, probably due to formation of an octahedral-coordination-geometry complex to the native state of EIsc, as occurs in other homologous proteins (56
,57
). The addition of PEP alone slightly destabilized the protein (Table 3), suggesting that the conformational changes monitored by trypsin digestion and quenching experiments did not result in large thermodynamic variations. Interestingly enough, the addition of both effectors to EIsc partially compensated for the stabilizing effect of Mg2+. We hypothesize that this compensatory effect could be due to a change in the coordination sphere of the Mg2+ by directly involving the phosphoryl group of the PEP moiety, as it occurs in a member of the pyruvate dikinase family (57
).
Thermodynamic and folding properties versus structure among the members of a protein family
Early studies suggested that folding and equilibrium pathways of homologous proteins were similar, and that the folding of a certain structure was thus conserved through evolution (58
,59
). However, there are examples of homologous proteins that unfold following different patterns and, in addition, have different thermodynamical parameters (60
,61
). In EIsc, although its sequence homology and its structure are similar to those of other EI proteins, the calorimetric equilibrium unfolding, as shown by DSC, is different. These changes in the thermodynamical parameters point to subtle conformational variations in the packing of some residues. The differences among the modeled structure and that of EICtt in the hydrogen-bond pattern, and the van der Waals contacts in the dimerization interface could explain those changes. At the moment, however, we do not know whether these small changes in the structure are also responsible for the observed considerable specificity of the phosphotransfer reaction between EI and HPr among the different species and the differences in their enthalpic variations in their binding reactions (45
,62
).
A similar reasoning can be used to explain the outcome of the effectors. The secondary structure of EIsc upon binding was not largely affected (Table 4), but there were subtle changes in the stability of the protein (Table 3), which suggest that small rearrangements occurred upon binding, as further confirmed by quenching fluorescence and trypsin digestion experiments.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
This work was supported by Projects CTQ2004-07716 from the Spanish Ministerio de Educación y Ciencia (MEC) to J.G., and Projects GV04B-402 from Generalitat Valenciana, and CTQ2004-04474 and CTQ2005-00360/BQU from MEC, to J.L.N. F.T. and H.N. were supported by grant SFB473 of the Deutsche Forschungsgemeinschaft. E.H.G. is a recipient of a predoctoral fellowship from Generalitat Valenciana.
Submitted on October 31, 2005; accepted for publication March 7, 2006.
| REFERENCES |
|---|
|
|
|---|
2. Dimitrova, M. N., A. Peterkofsky, and A. Ginsburg. 2003. Opposing effects of phosphoenolpyruvate and pyruvate with Mg2+ on the conformational stability and dimerization of phosphotransferase enzyme I from Escherichia coli. Protein Sci. 12:20472056.
3. Chauvin, F., L. Brand, and S. Roseman. 1994. Sugar transport by the bacterial phosphotransferase system. Studies on the molecular weight and association of Enzyme I. Characterization of the Escherichia coli enzyme I monomer/dimer transition kinetics by fluorescence anisotropy. J. Biol. Chem. 269:2027020274.
4. Dimitrova, M. N., R. H. Szczepanowski, S. B. Ruvinov, A. Peterkofsky, and A. Ginsburg. 2002. Interdomain interaction and substrate coupling effects of dimerization and conformational stability of Enzyme I of the Escherichia coli phosphoenolpyruvate:sugar phosphotransferase system. Biochemistry. 41:906913.[CrossRef][Medline]
5. Lux, R., K. Jahreis, K. Bettenbrock, J. S. Parkinson, and J. W. Lengeler. 1995. Coupling the phosphotransferase system and the methyl-accepting chemotaxis protein-dependent chemotaxis signaling pathways of Escherichia coli. Proc. Natl. Acad. Sci. USA. 92:1158311587.
6. LiCalsi, C., T. S. Crocenzi, E. Freire, and S. Roseman. 1991. Sugar transport by the bacterial phosphotransferase system. Structural and thermodynamics domains of Enzyme I of Salmonella typhimurium. J. Biol. Chem. 266: 1951919527.
7. Lee, B. R., P. Lecchi, L. Panell, H. Jaffe, and A. Peterkofsky. 1994. Identification of the N-terminal domain of enzyme I of the Escherichia coli phosphoenolpyruvate:sugar phosphotransferase system produced by proteolytic digestion. Arch. Biochem. Biophys. 312:121124.[CrossRef][Medline]
8. Liao, D. I., E. Silverton, Y. J. Seok, B. R. Lee, A. Peterkofsky, and D. R. Davies. 1996. The first step in sugar transport: crystal structure of the amino terminal domain of enzyme I of the E. coli PEP:sugar phosphotransferase system and a model of the phosphotransfer complex with HPr. Structure. 4:861872.[Medline]
9. Garret, D. S., Y. J. Seok, D. I. Liao, A. Peterkofsky, A. M. Gronenborg, and G. M. Clore. 1997. Solution structure of the 30 kDa N-terminal domain of enzyme I of the Escherichia coli phosphoenolpyruvate:sugar phosphotransferase system by multidimensional NMR. Biochemistry. 36:25172530.[CrossRef][Medline]
10. Garret, D. S., Y. J. Seok, A. Peterkofsky, A. M. Gronenborg, and G. M. Clore. 1999. Solution structure of the 40,000 Mr phosphoryl transfer complex between the N-terminal domain of enzyme I and HPr. Nat. Struct. Biol. 6:166173.[CrossRef][Medline]
11. Oberholzer, A. E., M. Bumann, P. Schneider, C. Bächler, C. Siebold, U. Baumann, and B. Erni. 2005. Crystal structure of the phosphoenolpyruvate-binding enzyme I domain from the Thermoanaerobacter tengcongensis PEP:sugar phosphotransferase system. J. Mol. Biol. 346:521532.[CrossRef][Medline]
12. Seok, Y. J., B. R. Lee, P. P. Zhu, and A. Peterkofsky. 1996. Importance of the carboxyl-terminal domain of enzyme I of the Escherichia coli phosphoenolpyruvate:sugar phosphotransferase system for phosphoryl donor specificity. Proc. Natl. Acad. Sci. USA. 93:347351.
13. Bentley, S. D., K. F. Chater, A. M. Cerdeno-Tarraga, G. L. Challis, N. R. Thomson, K. D. James, D. E. Harris, M. A. Quail, H. Kieser, D. Harper, A. Bateman, S. Brown, et al. 2002. Complete genome sequence of the model actinomycete Streptomyces coelicolor A3(2). Nature. 417:141147.[CrossRef][Medline]
14. Titgemeyer, F., J. Walkenhorst, X. Cui, J. Reizer, and M. H. Jr. Saier. 1994. Proteins of the phosphoenolpyruvate:sugar phosphotransferase system in Streptomyces: possible involvement in the regulation of antibiotic production. Res. Microbiol. 145:8992.[Medline]
15. Nothaft, H., D. Dresel, A. Willimek, K. Mahr, M. Niederweis, and F. Titgemeyer. 2003. The phosphotransferase system of Streptomyces coelicolor is biased for N-acetylglucosamine metabolism. J. Bacteriol. 185:70197023.
16. Parche, S., R. Schmid, and F. Titgemeyer. 1999. The PTS system of Streptomyces coelicolor: identification and biochemical analysis of a histidine phosphocarrier protein HPr encoded by the gene ptsH. Eur. J. Biochem. 265:308317.[Medline]
17. Fernandez-Ballester, G., J. Maya, A. Martin, S. Parche, J. Gomez, F. Titgemeyer, and J. L. Neira. 2003. The histidine-phosphocarrier protein of Streptomyces coelicolor folds by a partially folded species at low pH. Eur. J. Biochem. 270:22542267.[Medline]
18. Neira, J. L., and J. Gomez. 2004. The conformational stability of the Streptomyces coelicolor histidine-phosphocarrier protein. Characterization of cold denaturation and urea-protein interactions. Eur. J. Biochem. 271:21652181.[Medline]
19. Azuaga, A. I., J. L. Neira, and N. A. van Nuland. 2005. HPr as a model protein in structure, interaction, folding and stability studies. Protein Pept. Lett. 12:123137.[CrossRef][Medline]
20. Butler, M. J., J. Deutscher, P. W. Postma, T. J. Wilson, A. Galinier, and M. J. Bibb. 1999. Analysis of a ptsH homologue from Streptomyces coelicolor A3(2). FEMS Microbiol. Lett. 177:279288.