Absence of kinetic thermal stabilization in a hyperthermophile rubredoxin indicated by 40 microsecond folding in the presence of irreversible denaturation

Temperature-dependent iron motion in extremophile rubredoxins - no need for 'corresponding states'

Extremophile organisms are known that can metabolize at temperatures down to - 25 °C (psychrophiles) and up to 122 °C (hyperthermophiles). Understanding viability under extreme conditions is relevant for human health, biotechnological applications, and our search for life elsewhere in the universe. Information about the stability and dynamics of proteins under environmental extremes is an important factor in this regard. Here we compare the dynamics of small Fe-S proteins - rubredoxins - from psychrophilic and hyperthermophilic microorganisms, using three different nuclear techniques as well as molecular dynamics calculations to quantify motion at the Fe site. The theory of 'corresponding states' posits that homologous proteins from different extremophiles have comparable flexibilities at the optimum growth temperatures of their respective organisms. Although 'corresponding states' would predict greater flexibility for rubredoxins that operate at low temperatures, we find that from 4 to 300 K, the dynamics of the Fe sites in these homologous proteins are essentially equivalent.

Protein folding at extreme temperatures: Current issues

The range of temperatures compatible with life is currently estimated from -25°C, as exemplified by metabolically active bacteria between sea ice crystals, and up to 122°C in hydrothermal vents as exemplified by the archaeon Methanopyrus kandleri. In the context of protein folding, as soon as a polypeptide emerges from the ribosome, it is exposed to the effects of environmental temperatures. Recent investigations have shown that the rate of protein folding is not adapted to extreme temperatures and should be very fast at high temperature and low in cold environments. This lack of adaptation is driven by kinetic constraints on protein stability. To counteract the deleterious effects of fast protein folding in hyperthermophiles, chaperones such as the Trigger Factor hold and slow down the rate of folding intermediates. Prolyl isomerization, a rate-limiting step in the folding of many proteins, is strongly temperature-dependent and impairs folding of psychrophilic proteins in the cold. This is compensated by reduction of the proline content in cold-adapted proteins, by an increased number of prolyl isomerases encoded in the genome of psychrophilic microorganisms and by overexpression of prolyl isomerases under low temperature cultivation. After folding, the native state is reached and although extremophilic proteins share the same fold, dramatic differences in stability have been recorded by differential scanning calorimetry.

Solution-state NMR structure and biophysical characterization of zinc-substituted rubredoxin B (Rv3250c) from Mycobacterium tuberculosis

Owing to the evolution of multi-drug-resistant and extremely drug-resistant Mycobacterium tuberculosis strains, there is an urgent need to develop new antituberculosis strategies to prevent TB epidemics in the industrial world. Among the potential new drug targets are two small nonheme iron-binding proteins, rubredoxin A (Rv3251c) and rubredoxin B (Rv3250c), which are believed to play a role in electron-transfer processes. Here, the solution structure and biophysical properties of one of these two proteins, rubredoxin B (Mt-RubB), determined in the zinc-substituted form are reported. The zinc-substituted protein was prepared by expressing Mt-RubB in minimal medium containing excess zinc acetate. Size-exclusion chromatography and NMR spectroscopy indicated that Mt-RubB was a monomer in solution. The structure (PDB entry 2kn9) was generally similar to those of other rubredoxins, containing a three-stranded antiparallel β-sheet (β2-β1-β3) and a metal tetrahedrally coordinated to the S atoms of four cysteine residues (Cys9, Cys12, Cys42 and Cys45). The first pair of cysteine residues is at the C-terminal end of the first β-strand and the second pair of cysteine residues is towards the C-terminal end of the loop between β2 and β3. The structure shows the metal buried deeply within the protein, an observation that is supported by the inability to remove the metal with excess EDTA at room temperature. Circular dichroism spectroscopy shows that this stability extends to high temperature, with essentially no change being observed in the CD spectrum of Mt-RubB upon heating to 353 K.

Thermostability in rubredoxin and its relationship to mechanical rigidity

The source of increased stability in proteins from organisms that thrive in extreme thermal environments is not well understood. Previous experimental and theoretical studies have suggested many different features possibly responsible for such thermostability. Many of these thermostabilizing mechanisms can be accounted for in terms of structural rigidity. Thus a plausible hypothesis accounting for this remarkable stability in thermophilic enzymes states that these enzymes have enhanced conformational rigidity at temperatures below their native, functioning temperature. Experimental evidence exists to both support and contradict this supposition. We computationally investigate the relationship between thermostability and rigidity using rubredoxin as a case study. The mechanical rigidity is calculated using atomic models of homologous rubredoxin structures from the hyperthermophile Pyrococcus furiosus and mesophile Clostridium pasteurianum using the FIRST software. A global increase in structural rigidity (equivalently a decrease in flexibility) corresponds to an increase in thermostability. Locally, rigidity differences (between mesophilic and thermophilic structures) agree with differences in protection factors.

NMR analysis of native-state protein conformational flexibility by hydrogen exchange

The rate of hydrogen exchange for the most protected amides of a protein is widely used to provide an estimate of global conformational stability by analyzing the exchange kinetics in the unfolded state in terms of model peptide exchange rates. The exchange behavior of the other amides of the protein which do not exchange via a global unfolding mechanism can provide insight into the smaller-scale conformational transitions that facilitate access to solvent as required for the exchange reaction. However, since the residual tertiary structure in the exchange-competent conformation can modulate the chemistry of the exchange reaction, equilibrium values estimated from normalization with model peptide rates are open to question. To overcome this limitation, the most robust approaches utilize differential analyses as a function of experimental variables such as denaturant concentration, temperature, pH, and mutational variation. Practical aspects of these various differential analysis techniques are considered with illustrations drawn from the literature.

A billion-fold range in acidity for the solvent-exposed amides of Pyrococcus furiosus rubredoxin

The exchange rates of the static solvent-accessible amide hydrogens of Pyrococcus furiosus rubredoxin range from near the diffusion-limited rate to a billion-fold slower for the non-hydrogen-bonded Val 38 (eubacterial numbering). Hydrogen exchange directly monitors the kinetic acidity of the peptide nitrogen. Electrostatic solvation free energies were calculated by Poisson-Boltzmann methods for the individual peptide anions that form during the hydroxide-catalyzed exchange reaction to examine how well the predicted thermodynamic acidities match the experimentally determined kinetic acidities. With the exception of the Ile 12 amide, the differential exchange rate constant for each solvent-exposed amide proton that is not hydrogen bonded to a backbone carbonyl can be predicted within a factor of 6 (10 (0.78)) root-mean-square deviation (rmsd) using the CHARMM22 electrostatic parameter set and an internal dielectric value of 3. Under equivalent conditions, the PARSE parameter set yields a larger rmsd value of 1.28 pH units, while the AMBER parm99 parameter set resulted in a considerably poorer correlation. Either increasing the internal dielectric value to 4 or reducing it to a value of 2 significantly degrades the quality of the prediction. Assigning the excess charge of the peptide anion equally between the peptide nitrogen and the carbonyl oxygen also reduces the correlation to the experimental data. These continuum electrostatic calculations were further analyzed to characterize the specific structural elements that appear to be responsible for the wide range of peptide acidities observed for these solvent-exposed amides. The striking heterogeneity in the potential at sites along the protein-solvent interface should prove germane to the ongoing challenge of quantifying the contribution that electrostatic interactions make to the catalytic acceleration achieved by enzymes.

NMR and X-ray analysis of structural additivity in metal binding site-swapped hybrids of rubredoxin

Background

Chimeric hybrids derived from the rubredoxins of Pyrococcus furiosus (Pf) and Clostridium pasteurianum (Cp) provide a robust system for the characterization of protein conformational stability and dynamics in a differential mode. Interchange of the seven nonconserved residues of the metal binding site between the Pf and Cp rubredoxins yields a complementary pair of hybrids, for which the sum of the thermodynamic stabilities is equal to the sum for the parental proteins. Furthermore, the increase in amide hydrogen exchange rates for the hyperthermophile-derived metal binding site hybrid is faithfully mirrored by a corresponding decrease for the complementary hybrid that is derived from the less thermostable rubredoxin, indicating a degree of additivity in the conformational fluctuations that underlie these exchange reactions.

Results

Initial NMR studies indicated that the structures of the two complementary hybrids closely resemble "cut-and-paste" models derived from the parental Pf and Cp rubredoxins. This protein system offers a robust opportunity to characterize differences in solution structure, permitting the quantitative NMR chemical shift and NOE peak intensity data to be analyzed without recourse to the conventional conversion of experimental NOE peak intensities into distance restraints. The intensities for 1573 of the 1652 well-resolved NOE crosspeaks from the hybrid rubredoxins were statistically indistinguishable from the intensities of the corresponding parental crosspeaks, to within the baseplane noise level of these high sensitivity data sets. The differences in intensity for the remaining 79 NOE crosspeaks were directly ascribable to localized dynamical processes. Subsequent X-ray analysis of the metal binding site-swapped hybrids, to resolution limits of 0.79 Å and 1.04 Å, demonstrated that the backbone and sidechain heavy atoms in the NMR-derived structures lie within the range of structural variability exhibited among the individual molecules in the crystallographic asymmetric unit (~0.3 Å), indicating consistency with the "cut-and-paste" structuring of the hybrid rubredoxins in both crystal and solution.

Conclusion

Each of the significant energetic interactions in the metal binding site-swapped hybrids appears to exhibit a 1-to-1 correspondence with the interactions present in the corresponding parental rubredoxin structure, thus providing a structural basis for the observed additivity in conformational stability and dynamics. The congruence of these X-ray and NMR experimental data offers additional support for the interpretation that the conventional treatment of NOE distance restraints contributes substantially to the systematic differences that are commonly reported between NMR- and X-ray-derived protein structures.

Normal carbon acid referencing for protein amide hydrogen exchange

Measurement of the exchange kinetics for amide hydrogens along the protein backbone continues to offer valuable insights into structural stability and conformational dynamics. Since such studies routinely compare samples that differ in solution conditions or mechanical handling, normalization of the relative exchange rates can present a potentially significant source of experimental uncertainty. The carbon acids 1,3-dimethylimidazolium cation and thiomethylacetonitrile exhibit base catalyzed exchange rates similar to those of the slowly exchanging amides, under conditions typical for protein studies. With 13C enrichment at the acidic carbon position to facilitate selective observation, such carbon acids offer practical internal calibration of exchange.

Residue cluster additivity of thermodynamic stability in the hydrophobic core of mesophile vs. hyperthermophile rubredoxins

The branched sidechain residues 24 and 33 in the hydrophobic core of rubredoxin differ between the Clostridium pasteurianum (Cp) and Pyrococcus furiosus (Pf) sequences. Their X-ray structures indicate that these two sidechains are in van der Waals contact with each other, while neither appears to significantly interact with the other nonconserved residues. The simultaneous interchange of residues 24 and 33 between the Cp and Pf rubredoxin sequences yield a complementary pair of hybrid proteins for which the sum of their thermodynamic stabilities equals that of the parental rubredoxins. The 1.2 kcal/mol change arising from this two residues interchange accounts for 21% of the differential thermodynamic stability between the mesophile and hyperthermophile proteins. The additional interchange of the sole nonconserved aromatic residue in the hydrophobic core yields a 0.78 kcal/mol deviation from thermodynamic additivity.

Electrostatic potential energy within a protein monitored by metal charge-dependent hydrogen exchange

Hydrogen exchange measurements on Zn(II)-, Ga(III)-, and Ge(IV)-substituted Pyrococcus furiosus rubredoxin demonstrate that the log ratio of the base-catalyzed rate constants (Delta log k(ex)) varies inversely with the distance out to at least 12 A from the metal. This pattern is consistent with the variation of the amide nitrogen pK values with the metal charge-dependent changes in the electrostatic potential. Fifteen monitored amides lie within this range, providing an opportunity to assess the strength of electrostatic interactions simultaneously at numerous positions within the structure. Poisson-Boltzmann calculations predict an optimal effective internal dielectric constant of 6. The largest deviations between the experimentally estimated and the predicted DeltapK values appear to result from the conformationally mobile charged side chains of Lys-7 and Glu-48 and from differential shielding of the peptide units arising from their orientation relative to the metal site.

Tetrathiolate coordination of germanium(IV) in a protein active site

The tetracysteine metal coordination site of the rubredoxins from Clostridium pasteurianum (Cp) and Pyrococcus furiosus (Pf) are shown to stably bind the inorganic Ge(IV) ion. This is the first characterized coordination complex of tetravalent germanium with a biological macromolecule. Zn(II), Ga(III) and Ge(IV) substitution yields differential NMR chemical shifts for the 1H and 15N amide resonances throughout much of the protein structure. The differential shifts for the six backbone amides that hydrogen bond to the metal-coordinated sulfurs indicate that the pseudo 2-fold symmetry of the active site is more closely maintained in the hyperthermophile Pf rubredoxin than in its mesophile Cp homolog. These three metal substitutions form an isoelectronic series of small diamagnetic proteins for which reference structures are known to 1A resolution. These series provide a promising system to analyze theoretical predictions of the effects of differential charge distribution on chemical shifts from both proximal and long range interactions.

Enhanced thermal stability achieved without increased conformational rigidity at physiological temperatures: Spatial propagation of differential flexibility in rubredoxin hybrids

The extreme thermal stability of proteins from hyperthermophilic organisms is widely believed to arise from an increased conformational rigidity in the native state. In apparent contrast to this paradigm, both Pyrococcus furiosus (Pf) rubredoxin, the most thermostable protein characterized to date, and its Clostridium pasteurianum (Cp) mesophile homolog undergo a transient conformational opening of their multi-turn segments, which is more favorable in hyperthermophile proteins below room temperature. Substitution of the hyperthermophile multi-turn sequence into the mesophile protein sequence yields a hybrid, (14-33(Pf)) Cp, that exhibits a 12 degrees increase in its reversible thermal unfolding transition midpoint. Nuclear magnetic resonance (NMR) magnetization transfer-based hydrogen exchange was used to monitor backbone conformational dynamics in the subsecond time regime. Despite the substantially increased thermostability, flexibility throughout the entire main chain of the more thermostable hybrid is equal to or greater than that of the wild type mesophile rubredoxin near its normal growth temperature. In comparison to the identical core residues of the (14-33(Pf)) Cp rubredoxin hybrid, six spatially clustered residues in the parental mesophile protein exhibit a substantially larger temperature dependence of exchange. The exchange behavior of these six residues closely matches that observed in the multi-turn segment, consistent with a more extensive conformational process. These six core residues exhibit a much weaker temperature dependence of exchange in the (14-33(Pf)) Cp hybrid, similar to that observed for the multi-turn segment in its parental Pf rubredoxin. These results suggest that differential temperature dependence of flexibility can underlie variations in thermostability observed for mesophile versus hyperthermophile homologs.

Additivity in Both Thermodynamic Stability and Thermal Transition Temperature for Rubredoxin Chimeras via Hybrid Native Partitioning

Given any operational definition of pairwise interaction, the set of residues that differ between two structurally homologous proteins can be uniquely partitioned into subsets of clusters for which no such interactions occur between clusters. Although hybrid protein sequences that preserve such clustering are consistent with tertiary structures composed of only parental native-like interactions, the stability of such predicted structures will depend upon the physical robustness of the assumed interaction potential. A simple distance cutoff criterion was applied to the most thermostable protein known to predict such a seven-residue cluster in the metal binding site region of Pyrococcus furiosus rubredoxin and a mesophile homolog. Both conformational stability and thermal transition temperature measurements demonstrate that 39% of the differential stability arises from these seven residues.

Contribution of the multi-turn segment in the reversible thermal stability of hyperthermophile rubredoxin: NMR thermal chemical exchange analysis of sequence hybrids

Pyrococcus furiosus (Pf) rubredoxin is the most thermostable protein characterized to date. Reflecting the complications arising from irreversible denaturation of this protein, predictions of which structural regions confer differential thermal stability have utilized kinetic stability measurements, hydrogen exchange protection factors, long range hydrogen bond NMR spin couplings, and molecular dynamics simulations, and have primarily implicated the three-stranded beta-sheet and the adjacent metal binding site. Herein, NMR chemical exchange experiments demonstrate reversible two-state unfolding at the thermal transition temperature (T(m)) for hybrids of Pf and the mesophile Clostridium pasteurianum (Cp) rubredoxins which interchange residues 14-33, the so-called multi-turn segment. This complementary pair of hybrid rubredoxins exhibits largely additive incremental thermal stabilizations vs. the parental proteins. Both stabilization free energy measurements as well as incremental T(m) values indicate that a minimum of 37% of the total differential thermal stability resides in this multi-turn segment. Such a proportionality between DeltaDeltaG and incremental T(m) values is predicted for hybrid pairs exhibiting thermodynamic additivity in which the differential stability is predominantly enthalpic.

Structural determinants of protein stabilization by solutes

Despite their high sequence homology, rubredoxins from Desulfovibrio gigas and D. desulfuricans are stabilized to very different extents by compatible solutes such as diglycerol phosphate, the major osmolyte in the hyperthermophilic archaeon Archaeoglobus fulgidus[Lamosa P, Burke A, Peist R, Huber R, Liu M Y, Silva G, Rodrigues-Pousada C, LeGall J, Maycock C and Santos H (2000) Appl Environ Microbiol66, 1974-1979]. The principal structural difference between these two proteins is the absence of the hairpin loop in the rubredoxin from D. desulfuricans. Therefore, mutants of D. gigas rubredoxin bearing deletions in the loop region were constructed to investigate the importance of this structural feature on protein intrinsic stability, as well as on its capacity to undergo stabilization by compatible solutes. The three-dimensional structure of the mutant bearing the largest deletion, Delta17/29, was determined by 1H-NMR, demonstrating that, despite the drastic deletion, the main structural features were preserved. The dependence of the NH chemical shifts on temperature and solute concentration (diglycerol phosphate or mannosylglycerate) provide evidence of subtle conformational changes induced by the solute. The kinetic stability (as assessed from the absorption decay at 494 nm) of six mutant rubredoxins was determined at 90 degrees C and the stabilizing effect exerted by both solutes was assessed. The extent of protection conferred by each solute was highly dependent on the specific mutant examined: while the half-life for iron release in the wild-type D. gigas rubredoxin increased threefold in the presence of 0.1 M diglycerol phosphate, mutant Delta23/29 was destabilized. This study provides evidence for solute-induced compaction of the protein structure and occurrence of weak, specific interactions with the protein surface. The relevance of these findings to our understanding of the molecular basis for protein stabilization is discussed.