|
Advertisement | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 280, Issue 17, 17397-17407, April 29, 2005
Energetics of Structural Transitions of the Addiction Antitoxin MazE
IS A PROGRAMMED BACTERIAL CELL DEATH DEPENDENT ON THE INTRINSICALLY FLEXIBLE NATURE OF THE ANTITOXINS?*
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|---|
|
|
|---|
4.5 and T < 50 °C, the unstructured polypeptide chains of the MazE dimer fold into (pre)molten globule-like conformations that thermally stabilize the dimeric form of the protein. The simulation based on the thermodynamic and structural information on various addiction modules suggests that both the conformational adaptability of the dimeric antitoxin form (binding to the toxins and DNA) and the reversible transformation to the more flexible monomeric form are essential for the regulation of bacterial cell life and death. | INTRODUCTION |
|---|
|
|
|---|
Toxin-antitoxin systems, some of which are homologous to these extrachromosomal addiction modules, have been found on the chromosomes of several bacteria. For example, in E. coli, there are several such pairs of genes, including mazEF (2325), chpBIK (24), relBE (2628), yefM-yoeB (2931), and dinJ-yafQ (32). The first toxin-antitoxin system contained in a bacterial chromosome that was described as regulatable and responsible for programmed cell death was the E. coli mazEF module (25), located in the relA operon (23). The product of mazF (MazF) is a stable toxin, and the product of mazE (MazE) is a labile protein, degraded by the protease ClpAP (25). Therefore, preventing MazF-mediated death requires a continuous cellular production of MazE. In contrast to the addiction modules that are triggered by the loss of the element, death mediated by chromosomal mazEF is achieved by several stressful conditions that prevent mazEF expression. Initially, the mazEF module was found to be under the control of ppGpp (25), the amino acid starvation signal molecule produced by the RelA protein (33). Overproduction of ppGpp leads to inhibition of the expression of mazEF and thereby to cell death (25, 34). However, as reported recently, inhibiting the expression of mazEF and thus inducing cell death can be also achieved by several additional stressful conditions. These include the inhibition of transcription and/or translation (3537), DNA damage (38, 39), and oxidative stress (39).
In the last few years, the chromosomal mazEF toxin-antitoxin systems have attracted additional attention, in particular in the following four directions. (a) The first is studies on the mechanism of MazF action. MazF has been shown to inhibit translation by cleaving mRNAs at specific sites (40, 41). However, two laboratories have reported contradictory results for the mechanism of cleavage. According to one group, MazF inhibits translation by a ribosome-dependent mechanism (40), and according to a second group, MazF is a sequence-specific endoribonuclease that cleaves mRNAs at ACA sequences in a codon- and ribosome-independent manner (41). This second version was further confirmed by a more recent report from the same laboratory showing that, at least in vitro, MazF functions enzymatically in a similar manner to RNase A, although with different sequence specificity (42). (b) The second direction is studies on the mazEF-mediated cell death network. There is a "point of no return" in MazF-mediated cell death; overexpression of MazE can only reverse MazF lethality over a short window of time (43), suggesting that MazF is a mediator rather than an executioner of cell lethality (44). (c) The third direction is studies on the function of mazEF showing that it can prevent spreading of phage infection (45). (d) The fourth is molecular studies characterizing the forces that drive the unfolding (folding) of proteins (MazE, MazF) and the interactions that drive MazE-MazF, MazE-DNA, and MazE-MazF-DNA complex formation. The first step toward this goal was achieved by crystallization of the antitoxin MazE and by determination of part of its three-dimensional structure (46). The resulting structure is the first of any addiction antitoxin, the crystallization of which has been hampered because of its susceptibility to proteases and its high content of unstructured polypeptide. We circumvented this problem by crystallizing MazE in complex with the specific dromedary antibody fragment VHH. Furthermore, the energetics of MazE interactions with VHH was investigated in detail. It has been suggested that, under physiological conditions in solution, a significant part of the MazE polypeptide chain may be unstructured, as in the crystal of the MazE-VHH complex (47). Almost simultaneously, Kamada et al. (48) determined the crystal structure of the MazE-MazF complex, in which MazE chains appear to be less unstructured than in the complex with VHH.
It is known that proteins that lack intrinsic globular structure (intrinsically flexible, intrinsically unstructured, natively unfolded, pliable) are frequently involved in some of the most important regulatory functions in the cell (30, 4955). Numerous examples of domains that are unstructured in solution but become structured upon binding to the target have been identified. Thus, the intrinsic lack of structure can offer functional advantages to a protein, including its ability to bind several different targets (30, 4957).
In this study, we investigated to what extent MazE belongs to this category of proteins by studying its conformational changes induced by variation in temperature, pH, protein concentration, and urea concentration. The energetics of MazE order-disorder transitions was studied by various techniques (differential scanning calorimetry, spectropolarimetry, and fluorometry) and is discussed in terms of its three-dimensional structures (46, 48). An attempt was made to predict some functional features of the addiction modules from the equilibrium thermodynamic model, which involves antitoxin and toxin molecules and DNA. To the best of our knowledge, the research described here represents for the first time a complete thermodynamic study of structural transitions of a protein with an intrinsically flexible nature.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
Differential Scanning Calorimetry (DSC)
The thermally induced transitions of MazE were measured between pH 1.5 and 9.4 using a Nano-II DSC differential scanning calorimeter (Calorimetry Sciences Corp.). The heating rate was 1 °C/min1, and the concentration of proteins in the measuring cell was
1.2 mg ml1. To obtain the presented thermograms ((Cp Cp,N) versus temperature curves) (see Fig. 2), the heat capacity of the protein in the initial (native) state was subtracted from the raw signal corrected for buffer contribution. The transition enthalpies (
Hcal) were obtained by integration of the (Cp Cp,N) versus temperature curves.
CD Spectropolarimetry
CD spectra were recorded using a Model 62A DS spectropolarimeter (Aviv, Lakewood, NJ) equipped with a programmable thermoelectrically controlled cell holder. The presented mean residue ellipticities ([
], degrees cm2 dmol1) were obtained from raw data (ellipticities) by subtracting the corresponding spectra of the buffer solution and taking into account the mass concentration (c), molecular mass (Mw), number of MazE amino acid residues (N), and optical path length (l) through the relation [
] =
·Mw/(N·c·l).
Thermal DenaturationChanges in MazE secondary (far-UV, 205255 nm; 1-mm cuvette) and tertiary (near-UV, 255310 nm; 1-cm cuvette) structure were followed by recording CD spectra at a MazE concentration of 0.5 mg ml1 (see Fig. 3, a and b) at different temperatures with a temperature step of 3 °C. At a lower concentration (0.05 mg ml1), temperature scans in the far-UV region were performed also at
= 230 nm in a 1-cm cuvette.
Urea DenaturationThe changes in MazE secondary structure upon increasing urea concentrations (08 M) were followed by measuring the ellipticity at 225 nm in a 1-cm cuvette at a MazE concentration of 0.05 mg ml1 (see Fig. 5, b and c). The measurements were performed at a number of different temperatures.
Fluorescence Spectroscopy
Fluorescence emission spectra were recorded using a PerkinElmer Life Sciences LS 50 luminescence spectrometer equipped with a thermally controlled cell holder and a cuvette with a 1-cm path length. Thermal denaturation of MazE was followed by the emission fluorescence of the external fluorescent probe 1-anilino-8-naphthalenesulfonic acid ammonium salt (ANS; Fluka, Buchs, Switzerland). The emission spectra (415615 nm) were recorded at different temperatures the range 890 °C using
ex = 400 nm (see Fig. 4). The concentrations of MazE and ANS in the buffer solutions were 0.05 and 0.126 mg ml1, respectively. Urea-induced denaturation of MazE was followed by measuring the changes in MazE intrinsic (Trp) emission fluorescence between 310 and 410 nm (
ex = 280 nm) at urea concentrations of 08 M. The concentration of MazE was 0.05 mg ml1 (see Fig. 5b), and the measurements were performed at a number of different temperatures.
Model Analysis of Thermally Induced and Urea-induced Transitions
Our results reveal that the majority of the monitored structural transitions can be adequately described in terms of a reversible two-state process in which the protein molecule can exist either in macro state N or in macro state D. Such process can be presented as shown in Equation 1,
![]() | (Eq. 1) |
G0 is the corresponding standard Gibbs free energy change per mol of monomer; n is the number of subunits to which each protein dissociates upon denaturation;
is the fraction of protein in the denatured state; and CT is the total protein concentration given in moles of monomer liter1. It should be noted that, under some specific conditions, two successive transitions that do not overlap significantly were observed (see "Results"). In such cases, the two transitions can be treated separately: the first one that corresponds to the transition of the native dimer N2 into the intermediate dimer I2 (1/2N2
1/2I2) is a special case of Equation 1 with n = 1, and the second one that describes the transition of the intermediate dimer I2 into the denatured monomer D (1/2I2
D) is a common example of Equation 2 with n = 2. The thermodynamic parameters of denaturation obtained for the two successive transitions are practically the same as those obtained from the model that takes into account all three states simultaneously (1/2N2
1/2I2
D) (see the supplemental material).
G0 can always be expressed by the integrated Gibbs-Helmholtz equation (Equation 2),
![]() | (Eq. 2) |
H0(T
) is the standard enthalpy of denaturation at the reference temperature T
(transition temperature at
= 0.5) and
Cp0 is the corresponding standard heat capacity change assumed to be independent of temperature. According to the model (Equation 1), one can express an average of a physical property, F (in our case, F is the partial molar enthalpy of the protein, the mean residue ellipticity ([
]), or the fluorescence intensity (FL)), in terms of the corresponding contributions FN and FD, which characterize states N and D (47, 59), respectively (Equation 3).
![]() | (Eq. 3) |
In the case of transitions followed by spectroscopic methods (CD and fluorescence), FN and FD are assumed to be linear functions of temperature (thermal denaturation) or urea concentration (urea denaturation). The model function (Cp Cp,N) for the DSC signal is the partial molar heat capacity of the protein relative to state N. It can be derived from the first partial derivative of Equation 3 on temperature at constant pressure (6062) (Equation 4).
![]() | (Eq. 4) |
Taking into account Equations 1, 2, 3, 4, the observed temperature profiles (melting curves) can be described in terms of the parameters
H0(T
),
Cp0, and T
. Their values were obtained from fitting the model function (CD = Equation 3 and DSC = Equation 4) to the experimental temperature profiles using the Levenberg-Marquardt nonlinear
2 regression procedure (63). In the case of urea-induced denaturation at a given temperature, the denaturation profiles were described by Equations 1 and 3 combined with the well known empirical relation
G0 =
Go0 m[urea], where
Go0 is the standard Gibbs free energy in the absence of urea and m is the proportionality coefficient. At each temperature, the parameters
Go0 and m were adjusted according to the
2 regression procedure mentioned above. These values were then fitted by Equation 2 (for
G0 =
Go0) to obtain the corresponding parameters
H0(T
),
Cp0, and T
in the absence of urea.
Structure-based Thermodynamic Calculations
The non-polar (AN) and polar (AP) solvent-accessible surface areas of proteins were calculated with NACCESS Version 2.1 (64). AN and AP of native MazE were obtained from the crystal structure of the MazE-VHH complex (probe size of 1.4 Å), and the unstructured residues were assumed to be exposed in the same way as they are in the extended Ala-X-Ala tripeptide. AN and AP contained by the denatured MazE were estimated as the sum of the corresponding accessibilities of the protein residues in an extended Ala-X-Ala tripeptide.
Cp0 accompanying MazE denaturation and dissociation (1/2N2
D) was calculated from the corresponding changes in non-polar and polar accessible areas using the expression introduced by Murphy and Freire (65) (Equation 5).
![]() | (Eq. 5) |
Cp0 values for MazE denaturation and dissociation (1/2N2
D). The enthalpy change accompanying the dissociation and denaturation of MazE (1/2N2
D) was calculated using the parameterized expression for
H0 introduced by Xie and Freire (69) (Equation 6),
![]() | (Eq. 6) |
H0 value observed with most globular proteins at their median transition temperature of 60 °C. The corresponding entropy change was calculated as a sum of four contributions (6971) (Equation 7).
![]() | (Eq. 7) |
Ssol0, which describes the exposure of polar and non-polar groups to the solvent upon protein denaturation and dissociation, was obtained as
Ssol0 =
Cp0 ln(T/385.15) (71, 72). The second term that reflects the change in the side chain conformational entropy upon dissociation (1/2N2
N) was calculated as
Ssc0 =
(
AscSsc0/AAla-X-Ala). The sum was taken over each amino acid in the protein-protein interface (excluding alanine, proline, and disulfide-bonded cysteine), and its side chain conformational entropy was scaled by its change in the accessible surface area (
Asc) normalized to its A in the corresponding Ala-X-Ala tripeptide (70). The Ssc0 values were those reported by Lee et al. (73). The third contribution (
Srt0) is due to mixing of the solvent and solute molecules when the number of solute molecules in the initial state differs from that in the final state. For the reaction 1/2N2
N carried out in a hypothetical 1 M standard state, this contribution amounts to
Srt0 = (1/2)R ln(1/55.5) (74). There is a considerable debate in the literature whether this contribution accounts for the entropy change due to the increase in translational/rotational degrees of freedom upon dissociation (67, 70, 71, 7585). Nevertheless, its value has been found as the most appropriate for description of protein-protein rigid body dissociations (70, 71, 75). The fourth contribution (
Sconf0) refers to the configurational entropy change that accompanies the unfolding of the protein (N
D). It is estimated as
Sconf0 =
N
·4.3 cal K1 (mol residue)1, where
N
is the number of amino acid residues participating in the unfolding process and 4.3 cal K1 (mol residue)1 is the average configurational entropy change associated with the passage of an amino acid from a buried to a solvent-exposed state (65). To correlate the energetics of denaturation of MazE with its structural features, the
AN and
AP values were estimated also from Equations 5 and 6 using the experimentally obtained
H0 and
Cp0 values (see Fig. 8c). Moreover,
N
was estimated from Equation 7 as
N
= (
S0 (
Ssol0 +
Ssc0 +
Srt0))/4.3 cal K (mol residue)1, where
S0 is the measured standard entropy change (see Fig. 9d).
Mean Hydrophobicity and Mean Net Charge
The hydrophobicity of each amino acid sequence was calculated by the Kyte and Doolittle approximation (86) using a window size of five amino acids. The hydrophobicity of individual residues was normalized to a scale of 01 in these calculations. The mean hydrophobicity
H
b is defined as the sum of the normalized hydrophobicities of all residues divided by the number of residues in the polypeptide. The mean net charge
R
is defined as the net charge at pH 7 divided by the total number of residues. The calculations for various addiction antitoxin and toxin sequences were performed using the Swiss Institute of Bioinformatics server ExPASy (available at www.expasy.ch) (87).
Distribution of Molecular Species Involved in Addiction Modules
The model for calculation of fractions of antitoxin and toxin species is defined by the postulates that were extracted from the thermodynamic and structural information on addiction modules (10, 11, 4648, 8892) and from the physical properties of an E. coli cell (available at redpoll.pharmacy.ualberta.ca/CCDB/cgi-bin/STAT_NEW.cgi). The postulates and the details on the numerical procedure are given in the Supplemental Material.
| RESULTS |
|---|
|
|
|---|
The extent of reversibility of the observed thermally induced transitions was checked by performing two consecutive temperature scans (DSC and CD) and/or by measuring spectra (CD, intrinsic, and ANS fluorescence) after cooling down the MazE sample to the pre-transitional temperature. Reversibility of urea denaturation was examined by dialyzing concentrated urea solutions against the appropriate buffer and measuring an optical property (CD and fluorescence) before and after the dialysis. The extent of reversibility of all monitored transitions (except for the high temperature one at pH 1.5) is >80% which makes the laws of reversible thermodynamics applicable for their description. By contrast, the model-dependent parameters (Equation 1) that describe the high temperature transition at pH 1.5 (reversibility of
30%) can serve only for comparative purposes (Table I).
|
25 °C (Fig. 2). The high temperature transitions can be described in terms of the equilibrium two-state model (Equations 1, 2, and 4) with n = 2, whereas for the description of low temperature transitions, a two-state model with n = 1 should be used. The low temperature transitions are characterized at pH 1.5 by the model-independent (
Hcal = 9 kcal mol1) and model-dependent (
H0(T
) = 9 kcal mol1 and T
= 24 °C) quantities, whereas at pH 3.1, the corresponding values are
Hcal = 12 kcal mol1,
H0(T
) = 13 kcal mol1, and T
= 28 °C.
Cp0 values were too small to be determined accurately and were therefore neglected in our model analysis. The successive low and high temperature transitions of MazE observed at pH 1.5 and 3.1 were described also by a three-state model (1/2N2
1/2I2
D; see the Supplemental Material). The resulting
H0(T
),
Cp0, and T
values for both transitions were practically the same as those obtained when each transition was considered as an independent two-state transition. This agreement clearly indicates that the overlap of the two successive transitions is insignificant, and therefore, MazE denaturation can be discussed in terms of one or two independent two-state processes. Considering these facts, the thermodynamic stability of MazE at different conditions was described by the corresponding parameters (
H0(T
),
Cp0, and T
) derived from the model analysis of the measured DSC data (see Figs. 6 and 7 and Table I).
CD SpectropolarimetryThe far-UV CD spectra measured at pH 7.1 and T = 37 °C (Fig. 1c) show that the fractions of
-helix,
-sheet, and the remainder of the structure estimated according to Provencher and Glöckner (93) are consistent with those observed in the structure of MazE within the crystal of the MazE-VHH complex (Fig. 1a). Furthermore, the far-UV CD spectra of MazE (Figs. 1c and 8a) show an increase in the CD signal upon lowering the pH from 7.1 to 1.5, which is typical for intrinsically flexible proteins (51). One may speculate that the increase in the CD signal at pH 1.5 observed at
< 205 nm can be ascribed to the increased amount of the coil-like conformation, whereas at
> 210, the increased CD intensity may be due to some other type of structuring. According to the spectra measured after dialysis of MazE solution at pH 1.5 against the appropriate buffer solution at pH 7.1, this low pH-induced transition is highly reversible (
90%). The far-UV CD intensity at
> 210 nm increased with temperature at all pH values studied (Figs. 1c, 3a, and 8a). We believe that such behavior that resembles structuring observed within intrinsically flexible proteins (51) may be ascribed to the structuring of the disordered part of MazE molecules. A relatively high intensity CD signal was observed at low temperature and pH >1.5 also in the near-UV range (asymmetry of aromatic residues in MazE), which indicates the presence of a substantial amount of tertiary structure in MazE dimers (Fig. 3b). With increasing temperature, this amount decreased at all monitored pH values; however, a substantial amount of tertiary structure remains present even at high temperature (see Fig. 8b). Moreover, Fig. 3, a and b, shows that, also upon chemical denaturation (8 M urea at 25 °C), neither the secondary nor tertiary structure of MazE was fully destroyed.
|
|
= 230 and 280 nm, respectively, show at pH 1.59.4 one transition at T > 60 °C, whereas at pH
4.5, another transition at
25 °C was detected (Fig. 3, c and d). The melting curves describing the high temperature transitions show concentration dependence (Fig. 3, a, c, and d), indicating that these processes are not monomolecular (n > 1; see Equations 1, 2, 3). By contrast, low temperature transitions observed at pH
4.5 seem to be concentration-independent and can be characterized as monomolecular (n = 1). Some melting curves measured at higher concentrations do not level off at T < 100 °C (limitation of the method), which made the determination of the end (denatured) state as a linear function of temperature questionable (Equation 3). Consequently, the end states determined at lower concentrations were used in the model analysis of these curves. Fig. 3c shows a typical temperature profile at pH
4.5 measured at two different concentrations. The analysis of the three-state model (1/2N2
1/2I2
D; see the Supplemental Material) confirms the observations from DSC measurements that the two transitions are rather independent (Fig. 3d). Wherever possible, the melting curves were analyzed in terms of an appropriate model (Equations 1, 2, 3: low temperature transitions, n = 1; and high temperature transitions, n = 2). In these cases, the model function (Equation 3) displays good agreement with experimentally measured mean residue ellipticities (Fig. 3), and the resulting
H0(T
),
Cp0, and T
values agree well with those derived from DSC experiments (Table I).
ANS FluorescenceChanges in ANS fluorescence are frequently used to detect non-native intermediate conformations of globular proteins (9496). This is because such intermediates are characterized by the presence of solvent-exposed non-polar clusters that bind ANS molecules. As a result, an increase in ANS fluorescence intensity accompanied by a pronounced blue shift of the emission maximum is observed (95, 96). At pH
4.5, the observed changes in the intensity at T < 50 °C (Fig. 4) show that, under these conditions, at least a part of the MazE molecule exists in a (pre)molten globule-like conformation (52, 97). Furthermore, at pH 1.5, temperature-induced transitions at
25 °C (red shift of
3 nm) and 65 °C (red shift of
5 nm) were observed that correspond well to the low and high temperature transitions observed with DSC and CD (Fig. 4b). At pH >4.5, no significant changes in the intensity or position of the emission maximum with increasing temperature were observed, which indicates that, at pH >4.5, neither the native dimeric state nor the denatured monomeric state of MazE contains exposed clusters of non-polar residues.
Urea-induced Transitions
Urea-induced denaturation of Maze at pH 7.1 was monitored by intrinsic fluorescence and far-UV CD. In Fig. 5a, an insert of the MazE structure is presented that shows that Trp is the only aromatic amino acid residue in the structured part of the MazE molecule (there are two more Trp residues in the unstructured part) that is exposed to a large extent to the solvent. This is reflected in the MazE fluorescence emission spectra, which display, in the measured solution, (
em)max values that are close to those observed with pure Trp in water. Our measurements also show that urea-induced denaturation of MazE is characterized, contrary to many globular proteins, by an increase in the intrinsic fluorescence intensity accompanied by a small blue shift of
3 nm (Fig. 5b). This result is consistent with the measured CD spectra of MazE in urea solutions (Fig. 3, a and b), which show that, even in 8 M urea, neither the secondary nor tertiary structure of MazE disappears (98). Moreover, the observed blue shift and increase in fluorescence intensity may be indicative of a local structuring of residues around one, two, or all three Trp residues. Interestingly, a similar blue shift was observed also with the thermal denaturation of MazE in urea-free aqueous solutions. Fig. 5b also shows that the far-UV CD intensity decreased with increasing urea concentration and that urea-induced transitions, followed by either CD or fluorescence, can be described in terms of an appropriate two-state model (Equations 1 and 3 with n = 2).
|
|
| DISCUSSION |
|---|
|
|
|---|
G0 versus temperature) of MazE at pH 7.1 obtained from thermal and urea-induced denaturation are presented in Fig. 6 together with the corresponding enthalpy and entropy contributions. At T = 585 °C, the
G0 values obtained from urea-induced denaturation (
G0 =
Go0; see "Experimental Procedures") are slightly higher than those obtained from thermal denaturation. Because, in the urea-denatured state, the fraction of secondary and tertiary structure is substantially lower than in the thermally denatured state (Fig. 3), the corresponding difference in
G0 may be ascribed mainly to the enhanced hydrophobic effect accompanying the urea-induced denaturation, i.e. to the increased number of unfolded residues accompanied by an increase in the solvent-accessible surface area. This explanation is supported by the experimental observation that
Cp0 obtained from the urea denaturation (0.46 kcal mol1 K1) is significantly higher than
Cp0 determined from the thermal denaturation (0.23 kcal mol1 K1). In their excellent publications, Makhatadze and Privalov (66), Makhatadze and co-workers (99), Murphy and Freire (65), Xie and Freire (69), Spolar and Record (67), and Myers et al. (68) correlated protein structural features to the corresponding thermodynamics of unfolding. In their structure-based calculations of thermodynamic quantities of unfolding, they used empirical parameterization based largely upon changes in polar (
AP) and non-polar (
AN) surface areas exposed to the solvent.
AN and
AP are usually calculated from the structural data. Thus, their values depend on the structural definition of the initial (native) and final (unfolded) states of the protein. The crystal structure of MazE (46) may be a good approximation of the free MazE dimer native form in solution because VHH does not interrupt MazE-DNA association (46, 47) and does not interfere with MazF binding (see the figure in the Supplemental Material). This approximation is supported by the analysis of the far-UV CD spectra, which, under physiological conditions, predict a MazE structure (46) that is close to the one observed in the crystal of the MazE-VHH complex (Fig. 1). The final (unfolded) state of MazE may be defined as a state in which each MazE residue (X) is exposed to the solvent in the same way as in the corresponding Ala-X-Ala tripeptide. Using the described initial and final states, we obtained the thermodynamic profile (Equations 5, 6, 7) for the complete unfolding and dissociation of MazE (1/2N2
D) (Fig. 6), in which all 44 structured residues are involved in the changes of the solvent-exposed surfaces (
AN,U = 2907 Å2,
AP,U = 1498 Å2, and
N
U = 44). This calculation leads to much too high
Cp0 (0.92 kcal K1 mol1) and
Sconf0 values. As a result, the calculated
G0 value is considerably lower than the experimental one.
|
|
AN,
AP, and
N
that accompany MazE denaturation were calculated from experimental thermodynamic data (Equations 5, 6, 7). The calculations show that, upon denaturation and monomerization of MazE induced by urea (
AN = 1793 Å2,
AP = 1320 Å2, and
N
= 31), more residues unfolded, and more non-polar surface was exposed to solvent than in the case of thermal denaturation (
AN = 1113 Å2,
AP = 1035 Å2, and
N
= 21). The estimates are in complete agreement with our explanation of the discrepancy between
G0(urea) and
G0(thermal) and show that, at physiological pH, both thermally denatured and urea-denatured states are far from being completely unfolded (
N
< the total number of structured residues,
N
U = 44). We would like to point out that these
AN,
AP, and
N
values comprise errors of the empirical parameterization and those of the measured quantities and can therefore be considered only as reasonably good approximations (see Fig. 8, c and d). Finally, we estimated the average number of unfolded residues (
N
) simply as
N
= (
AN +
AP)·
N
U/(
AN,U +
AP,U), where
AN and
AP refer to the thermal or urea-induced denaturation of MazE, and
AN,U,
AP,U, and
N
U are defined as described above. Fig. 8d shows that
N
values determined in this way are surprisingly close to the corresponding
N
values calculated from Equation 7, in which the configurational entropy change upon actual unfolding is taken as
Sconf0 =
N
·4.3 cal K1 (mol residue)1 (65). We believe that the observed agreement between the
N
values obtained from independent parameterization of
Cp0 and
H0 on the one side and
S0 on the other strongly supports the parameterization of the experimental quantities
Cp0,
H0, and
S0 suggested by Equations 5, 6, 7. In other words, the suggested parameterization seems to be successful provided that the "true"
AN and
AP values are used.
As shown in Fig. 7a, the stability of the MazE dimer conformation strongly depends on pH. In the studied temperature range, the measured
G0 versus pH curve exhibits a maximum at pH
4.5. Such behavior may be explained in terms of different structural characteristics of the final (high temperature) monomeric state and the initial (low temperature) dimeric state. At physiological pH, the low temperature dimeric state contains the unstructured C-terminal half characterized by low mean hydrophobicity and low negative mean net charge (see Fig. 10). Lowering the pH results in reduction of the net charge, and consequently, the unstructured MazE chains start organizing themselves into a sort of a compact molten globule-like conformation. This structuring effect is opposed by a simultaneous change in the conformation of the structured N-terminal half (low mean net charge and high mean hydrophobicity) (see Fig. 10), which is responsible for dimer formation. Because of the strong favorable contribution of the hydrophobic effect (mainly hydrophobic residues are involved in dimer formation) (Fig. 1b), a surplus of a positive charge on the N-terminal half at low pH does not break the dimer apart; however, it reduces the fraction of the secondary and tertiary structure. Inspection of data describing the high temperature monomeric state of MazE shows that, relative to the low temperature state, the secondary structure complexity of the high temperature state is about the same at all measured pH values (Fig. 8a). By contrast, the amount of the high temperature tertiary structure to the low temperature state is reduced to about the same extent only at pH
4.5. At lower pH values, this reduction appears to be much smaller (Fig. 8b). It appears that the thermodynamic parameters of denaturation result from the compensating effect of the partial folding of the unstructured C-terminal half and monomerization coupled with the partial unfolding of the structured N-terminal half. Inspection of Table I shows that lowering the pH thermally stabilizes the dimeric form of MazE (T
at a given MazE concentration increases with decreasing pH). Because, at pH 1.57.1,
Cp0 and
H0 values decrease upon lowering the pH, the maximum stability reached at pH
4.5 is most likely due to the simultaneous lowering of
Sconf0. The suggested explanation is in full agreement with the calculations of
AN,
AP, and
N
(Fig. 8, c and d) based on the experimental
Cp0,
H0, and
S0 values. They show that the net effect of pH and temperature is such that, at high temperature, fewer residues unfold, and fewer non-polar and polar surfaces become exposed to the solvent at pH
4.5 than at physiological pH. Moreover,
AN,
AP, and
N
values are much lower than those calculated for the complete unfolding of MazE, thus indicating that, irrespective of pH, the high temperature state of MazE is not a fully unfolded state.
|
50% of the structured polypeptide chains that contain a substantial amount of tertiary structure at physiological temperature and pH, we cannot characterize the MazE dimer strictly as an intrinsically flexible (natively unfolded) protein (4955, 97). However, to describe the changes in its spectral and thermodynamic properties, we have to take into account the intrinsically flexible nature of its C-terminal half (4955, 97, 100). Furthermore, the flexible nature of MazE is hidden in the dimer-monomer equilibrium because the fractions of dimer and monomer depend on the total (dimer + monomer) MazE concentration, CT (Fig. 7b).
Functional Stability of the Addiction AntitoxinsAddiction antitoxins and toxins can exist in the cell as free monomers and free dimers and can also form various antitoxin-toxin, antitoxin-DNA, and antitoxin-toxin-DNA complexes. If we want to know which species are important for regulation of cell death (or life), we have to understand its molecular mechanism. This can be achieved only by knowledge of the corresponding thermodynamic data on protein unfolding (folding) and specific protein-protein and protein-DNA interactions. Information of this type is rather scarce, but nevertheless, some thermodynamic characteristics common to addiction modules are available (10, 11, 47, 8991). For example, the measured
G0 for the MazE dimer-monomer transition is very similar to the corresponding
G0 values observed for the antitoxins CcdA (89) and ParD (91). Toxins have much higher thermodynamic stability (89). Moreover, it is known that the dimeric antitoxin (A2) has two binding sites for the dimeric toxin (T2) and vice versa, so they can form complexes of the type... A2T2A2T2... = A2iT2j (46, 48, 90, 92). The expression of the antitoxin and toxin can be inhibited either by low affinity binding of the antitoxin to the promoter/operator DNA (10, 11, 47, 89, 90) or by binding of the antitoxin-toxin complexes to the same part of the DNA with higher affinity (10, 48, 90, 92). According to these observations, one may expect that the strength of the specific interactions (binding affinities) within various addiction modules would be of similar magnitude. This led us to the construction of the equilibrium model that is able to predict some features essential for understanding the function of addiction modules (for details, see the Supplemental Material). The first aim of its application was to estimate the total concentrations of antitoxin ([A]T) and toxin ([T]T) under ordinary cellular conditions (in the absence of stress). It was assumed that, under steady-state conditions (the rate of antitoxin and toxin synthesis is the same as the rate of their degradation), the total concentrations of the antitoxin and toxin are about the same and that no toxin is bound to its cellular target. Because the total concentration of DNA (encoding for the antitoxin and toxin) in the cell can be estimated as 2.5 nM, the antitoxin and toxin concentrations were varied simultaneously from 108 to 105 M. Fig. 9a shows that the modeled addiction module showed the highest sensitivity to the changes in [A]T and [T]T in the 107 to 106 M range (maximum sensitivity (d(fraction)/d log[A]T) is at log-[A]T = log[T]T
6.6
[A]T = [T]T
2.5 x 107M). Consequently, it is not reasonable to expect that the cell would function at [A]T
[T]T > 106M, where practically all antitoxin- and toxin-antitoxin-binding sites on DNA are occupied (inhibition of expression). On the other hand, in the absence of stress, the cell would probably not function at [A]T
[T]T < 107 M, where almost all antitoxin and toxin molecules are in the unbound state and their total number (N) is very small (N = 4 at [A]T
[T]T = 108M). Fig. 9b shows that, in this range, the most populated antitoxin-toxin complexes contain more T2 than A2 molecules (A2T4 and A4T6). The A2T4 form was also observed experimentally (48, 90, 92). By contrast, DNA is in the 107 to 106 M range occupied mainly by complexes in which A2 and T2 molecules are at a 1:1 ratio (Fig. 9c) and favors binding of longer A2iT2j associates, which is in accordance with recently proposed structural models (46, 48). Fig. 9d indicates how the fractions of various species change with a [T]T/[A]T ratio at fixed [T]T = 2.5 x 107M. This simulation approximately shows the response of the addiction module under stress conditions when antitoxin and toxin expression is inhibited. With increasing [T]T/[A]T ratios, the concentration of free dimeric toxin rises, and thus, it becomes available for interaction with its cellular target. If the rate of degradation of the antitoxin were known (for CcdA, T
60 min in the presence of CcdB) (101), it would not be difficult to convert the [T]T/[A]T ratio to the time scale and to monitor the dying of the cell.
|
[T]T
2.5 x 107 M) dealing with
100 antitoxin and
100 toxin molecules. If this is the case, practically every single molecule is important. To keep the cell alive, the regulation of an addiction module has to be very precise and rapid because it must compensate for even very small fluctuations in the number of antitoxin and toxin molecules. Our model shows that, in a long enough period of time (cell division cycle of E. coli of
20 min), a single antitoxin molecule spends
40% of the time in complexes with toxin and DNA,
50% of time as a free dimer, and
10% of time in the monomeric form (Fig. 9a).
|
|
To the best of our knowledge, the numerical simulation presented here is the first attempt to combine thermodynamic and structural information on addiction modules in terms of a mathematical model. We are aware that some properties defined in the model differ from module to module. From the numerical simulation point of view (convergence), it is not a problem to increase the size (number) of the A2iT2j complexes and the number of DNA-binding sites or/and to include more specific binding modes (more binding constants) in the model. It is also possible to include the kinetics (degradation rates for antitoxin and toxin) and the binding of the toxin to the specific target. The model has been developed to incorporate specific properties of a particular addiction module and to simulate some of its actions. At present, there are too little thermodynamic and structural data available to propagate the model. However, we believe that even this limited size of the model comprehends the basic physical meaning of the interplay between the antitoxin, toxin, and DNA.
| FOOTNOTES |
|---|
The on-line version of this article (available at http://www.jbc.org) contains supplemental "Experimental Procedures," supplemental Equations 115, and a supplemental figure. ![]()
To whom correspondence should be addressed. Tel.: 386-1-241-9414; Fax: 386-1-241-9437; E-mail: jurij.lah{at}fkkt.uni-lj.si.
1 The abbreviations used are: HEPPS, 4-(2-hydroxyethyl)-1-piperazinepropanesulfonic acid; DSC, differential scanning calorimetry; ANS, 1-anilino-8-naphthalenesulfonic acid ammonium salt. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. Simic, N. De Jonge, R. Loris, G. Vesnaver, and J. Lah Driving Forces of Gyrase Recognition by the Addiction Toxin CcdB J. Biol. Chem., July 24, 2009; 284(30): 20002 - 20010. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. M. Simon, F. J. R. Sousa, R. Mohana-Borges, and G. C. Walker Regulation of Escherichia coli SOS mutagenesis by dimeric intrinsically disordered umuD gene products PNAS, January 29, 2008; 105(4): 1152 - 1157. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. M. Moritz and P. J. Hergenrother Toxin-antitoxin systems are ubiquitous and plasmid-encoded in vancomycin-resistant enterococci PNAS, January 2, 2007; 104(1): 311 - 316. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Engelberg-Kulka, R. Hazan, and S. Amitai mazEF: a chromosomal toxin-antitoxin module that triggers programmed cell death in bacteria J. Cell Sci., October 1, 2005; 118(19): 4327 - 4332. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Zhang, L. Zhu, J. Zhang, and M. Inouye Characterization of ChpBK, an mRNA Interferase from Escherichia coli J. Biol. Chem., July 15, 2005; 280(28): 26080 - 26088. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| All ASBMB Journals | Molecular and Cellular Proteomics |
| Journal of Lipid Research | ASBMB Today |