Stability and dynamics in a hyperthermophilic protein with melting temperature close to 200 degrees C

Identification and characterization of thermophilic amylase producing bacterial isolates from the brick kiln soil

The present experiment was designed to isolate bacterial strains from the brick kiln soil and to check the activity and enzyme kinetics of amylase from these isolates. The bacterial colonies were isolated from soil samples through the serial dilution method. The bacterial isolates were identified through morphological, electron microscopic and molecular analysis. The 16S ribosomal RNA sequences of the isolates IR-1, IR-2, IR-3, IR-8, and IR-9 showed high similarities with Bacillus tequilensis, Bacillus paramycoides, Proteus alimentorum, Bacillus wiedmannii, and Pseudomonas aeruginosa, respectively. All of the bacterial isolates showed a positive catalase activity except IR-9. Furthermore, the isolates showed variable antagonistic effects against different bacterial pathogens. All of the strains produced indole acetic acid (IAA), and the concentrations increased in the presence of tryptophan application. The isolates showed the amylase enzyme activity and maximum activity of isolates was achieved in 4% starch concentration. The IR-9 isolate showed the highest amylase activity of 5.9 U/ml. The V max values of the extracellular amylase from different bacterial isolates ranged between 12.90 and 50.00 IU ml-1. The lowest K m value of 6.33 mg starch was recorded for IR-8 and the maximum K cat value of 2.50 min-1 was observed for IR-3. The amylase activity of the isolates was significantly affected by a range of different incubation time, temperature, and pH values. Further tests are required before the potential utilization of these isolates for amylase production, and in the biopesticide and biofertilizer applications.

Ligand interactions and the protein order-disorder energetic continuum

Intrinsically disordered proteins as computationally predicted account for ∼1/3 of eukaryotic proteomes, are involved in a plethora of biological functions, and have been linked to several human diseases as a result of their dysfunctions. Here, we present a picture wherein an energetic continuum describes protein structural and conformational propensities, ranging from the hyperstable folded proteins on one end to the hyperdestabilized and sometimes functionally disordered proteins on the other. We distinguish between proteins that are folding-competent but disordered because of marginal stability and those that are disordered due mainly to the absence of folding code-completing structure-determining interactions, and postulate that disordered proteins that are unstructured by way of partial population of protein denatured states represent a sizable proportion of the proteome.

Differential Effects of Hydrophobic Core Packing Residues for Thermodynamic and Mechanical Stability of a Hyperthermophilic Protein

Proteins from organisms that have adapted to environmental extremes provide attractive systems to explore and determine the origins of protein stability. Improved hydrophobic core packing and decreased loop-length flexibility can increase the thermodynamic stability of proteins from hyperthermophilic organisms. However, their impact on protein mechanical stability is not known. Here, we use protein engineering, biophysical characterization, single-molecule force spectroscopy (SMFS), and molecular dynamics (MD) simulations to measure the effect of altering hydrophobic core packing on the stability of the cold shock protein TmCSP from the hyperthermophilic bacterium Thermotoga maritima. We make two variants of TmCSP in which a mutation is made to reduce the size of aliphatic groups from buried hydrophobic side chains. In the first, a mutation is introduced in a long loop (TmCSP L40A); in the other, the mutation is introduced on the C-terminal β-strand (TmCSP V62A). We use MD simulations to confirm that the mutant TmCSP L40A shows the most significant increase in loop flexibility, and mutant TmCSP V62A shows greater disruption to the core packing. We measure the thermodynamic stability (ΔGD-N) of the mutated proteins and show that there is a more significant reduction for TmCSP L40A (ΔΔG = 63%) than TmCSP V62A (ΔΔG = 47%), as might be expected on the basis of the relative reduction in the size of the side chain. By contrast, SMFS measures the mechanical stability (ΔG*) and shows a greater reduction for TmCSP V62A (ΔΔG* = 8.4%) than TmCSP L40A (ΔΔG* = 2.5%). While the impact on the mechanical stability is subtle, the results demonstrate the power of tuning noncovalent interactions to modulate both the thermodynamic and mechanical stability of a protein. Such understanding and control provide the opportunity to design proteins with optimized thermodynamic and mechanical properties.

The key to the extraordinary thermal stability of P. furiosus holo-rubredoxin: iron binding-guided packing of a core aromatic cluster responsible for high kinetic stability of the native structure

Pyrococcus furiosus rubredoxin (PfRd), a small, monomeric, 53 residues-long, iron-containing, electron-transfer protein of known structure is sometimes referred to as being the most structurally-stable protein known to man. Here, using a combination of mutational and spectroscopic (CD, fluorescence, and NMR) studies of differently made holo- and apo-forms of PfRd, we demonstrate that it is not the presence of iron, or even the folding of the PfRd chain into a compact well-folded structure that causes holo-PfRd to display its extraordinary thermal stability, but rather the correct iron binding-guided packing of certain residues (specifically, Trp3, Phe29, Trp36, and also Tyr10) within a tight aromatic cluster of six residues in PfRd's hydrophobic core. Binding of the iron atom appears to play a remarkable role in determining subtle details of residue packing, forcing the chain to form a hyper-thermally stable native structure which is kinetically stable enough to survive (subsequent) removal of iron. On the other hand, failure to bind iron causes the same chain to adopt an equally well-folded native-like structure which, however, has a differently-packed aromatic cluster in its core, causing it to be only as stable as any other ordinary mesophile-derived rubredoxin. Our studies demonstrate, perhaps for the very first time ever that hyperthermal stability in proteins can owe to subtle differences in residue packing vis a vis mesostable proteins, without there being any underlying differences in either amino acid sequence, or bound ligand status.

Conformational dynamics and stability of HP35 studied with 2D IR vibrational echoes

Two-dimensional infrared (2D IR) vibrational echo spectroscopy was used to measure the fast dynamics of two variants of chicken villin headpiece 35 (HP35). The CN of cyanophenylalanine residues inserted in the hydrophobic core were used as a vibrational probe. Experiments were performed on both singly (HP35-P) and doubly CN-labeled peptide (HP35-P(2)) within the wild-type sequence, as well as on HP-35 containing a singly labeled cyanophenylalanine and two norleucine mutations (HP35-P NleNle). There is a remarkable similarity between the dynamics measured in singly and doubly CN-labeled HP35, demonstrating that the presence of an additional CN vibrational probe does not significantly alter the dynamics of the small peptide. The substitution of two lysine residues by norleucines markedly improves the stability of HP35 by replacing charged with nonpolar residues, stabilizing the hydrophobic core. The results of the 2D IR experiments reveal that the dynamics of HP35-P are significantly faster than those of HP35-P NleNle. These observations suggest that the slower structural fluctuations in the hydrophobic core, indicating a more tightly structured core, may be an important contributing factor to HP35-P NleNle's increased stability.

Electrostatics of hydrogen exchange for analyzing protein flexibility

Electrostatic interactions at the protein-aqueous interface modulate the reactivity of solvent-exposed backbone amides by a factor of at least a billion fold. The brief (∼10 ps) lifetime of the peptide anion formed during the hydroxide-catalyzed exchange reaction helps enable the experimental rates to be robustly predictable by continuum dielectric methods. Since this ability to predict the structural dependence of exchange reactivity also applies to the protein amide hydrogens that are only rarely exposed to the bulk solvent phase, electrostatic analysis of the experimental exchange rates provides an effective assessment of whether a given model ensemble is consistent with the properly weighted Boltzmann conformational distribution of the protein native state.

Protein stability and enzyme activity at extreme biological temperatures

Psychrophilic microorganisms thrive in permanently cold environments, even at subzero temperatures. To maintain metabolic rates compatible with sustained life, they have improved the dynamics of their protein structures, thereby enabling appropriate molecular motions required for biological activity at low temperatures. As a consequence of this structural flexibility, psychrophilic proteins are unstable and heat-labile. In the upper range of biological temperatures, thermophiles and hyperthermophiles grow at temperatures > 100 °C and synthesize ultra-stable proteins. However, thermophilic enzymes are nearly inactive at room temperature as a result of their compactness and rigidity. At the molecular level, both types of extremophilic proteins have adapted the same structural factors, but in opposite directions, to address either activity at low temperatures or stability in hot environments. A model based on folding funnels is proposed accounting for the stability-activity relationships in extremophilic proteins.

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.

Dynamics in Thermotoga neapolitana adenylate kinase: 15N relaxation and hydrogen-deuterium exchange studies of a hyperthermophilic enzyme highly active at 30 degrees C

Backbone conformational dynamics of Thermotoga neapolitana adenylate kinase in the free form (TNAK) and inhibitor-bound form (TNAK*Ap5A) were investigated at 30 degrees C using (15)N NMR relaxation measurements and NMR monitored hydrogen-deuterium exchange. With kinetic parameters identical to those of Escherichia coli AK (ECAK) at 30 degrees C, TNAK is a unique hyperthermophilic enzyme. These catalytic properties make TNAK an interesting and novel model to study the interplay between protein rigidity, stability, and activity. Comparison of fast time scale dynamics (picosecond to nanosecond) in the open and closed states of TNAK and ECAK at 30 degrees C reveals a uniformly higher rigidity across all domains of TNAK. Within this framework of a rigid TNAK structure, several residues located in the AMP-binding domain and in the core-lid hinge regions display high picosecond to nanosecond time scale flexibility. Together with the recent comparison of ECAK dynamics with those of hyperthermophilic Aquifex aeolicus AK (AAAK), our results provide strong evidence for the role of picosecond to nanosecond time scale fluctuations in both stability and activity. In the slow time scales, TNAK's increased rigidity is not uniform but localized in the AMP-binding and lid domains. The core domain amides of ECAK and TNAK in the open and closed states show comparable protection against exchange. Significantly, the hinges framing the lid domain show similar exchange data in ECAK and TNAK open and closed forms. Our NMR relaxation and hydrogen-deuterium exchange studies therefore suggest that TNAK maintains high activity at 30 degrees C by localizing flexibility to the hinge regions that are key to facilitating conformational changes.

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.

Preliminary neutron crystallographic analysis of selectively CH3-protonated deuterated rubredoxin from Pyrococcus furiosus

Neutron crystallography is used to locate H atoms in biological materials and can distinguish between negatively scattering hydrogen-substituted and positively scattering deuterium-substituted positions in isomorphous neutron structures. Recently, Hauptman & Langs (2003; Acta Cryst. A59, 250-254) have shown that neutron diffraction data can be used to solve macromolecular structures by direct methods and that solution is aided by the presence of negatively scattering H atoms in the structure. Selective-labeling protocols allow the design and production of H/D-labeled macromolecular structures in which the ratio of H to D atoms can be precisely controlled. Methyl selective-labeling protocols were applied to introduce (1H-delta methyl)-leucine and (1H-gamma methyl)-valine into deuterated rubredoxin from Pyrococcus furiosus (PfRd). Here, the production, crystallization and preliminary neutron analysis of a selectively CH3-protonated deuterated PfRd sample, which provided a high-quality neutron data set that extended to 1.75 A resolution using the new LADI-III instrument at the Institut Laue-Langevin, are reported. Preliminary analysis of neutron density maps allows unambiguous assignment of the positions of H atoms at the methyl groups of the valine and leucine residues in the otherwise deuterated rubredoxin structure.

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.

Thermodynamic stability and folding of proteins from hyperthermophilic organisms

Life grows almost everywhere on earth, including in extreme environments and under harsh conditions. Organisms adapted to high temperatures are called thermophiles (growth temperature 45-75 degrees C) and hyperthermophiles (growth temperature >or= 80 degrees C). Proteins from such organisms usually show extreme thermal stability, despite having folded structures very similar to their mesostable counterparts. Here, we summarize the current data on thermodynamic and kinetic folding/unfolding behaviors of proteins from hyperthermophilic microorganisms. In contrast to thermostable proteins, rather few (i.e. less than 20) hyperthermostable proteins have been thoroughly characterized in terms of their in vitro folding processes and their thermodynamic stability profiles. Examples that will be discussed include co-chaperonin proteins, iron-sulfur-cluster proteins, and DNA-binding proteins from hyperthermophilic bacteria (i.e. Aquifex and Theromotoga) and archea (e.g. Pyrococcus, Thermococcus, Methanothermus and Sulfolobus). Despite the small set of studied systems, it is clear that super-slow protein unfolding is a dominant strategy to allow these proteins to function at extreme temperatures.

Spatial distribution of dielectric shielding in the interior of Pyrococcus furiosus rubredoxin as sampled in the subnanosecond timeframe by hydrogen exchange

Experimental pK values of ionizable sidechains provide the most direct test for models representing dielectric shielding within the interior of a protein. However, only the strongly shifted pK values are particularly useful for discriminating among models. NMR titration studies have usually found only one or two such shifted pK values in each protein, so that the fitting of the experimental data to a uniform internal dielectric (epsilon(int)) model is not well constrained. The observed variation among proteins for such epsilon(int) estimates may reflect nonuniformity of dielectric shielding within each protein interior or qualitative differences between individual proteins. The differential amide kinetic acidities for a series of metal-substituted rubredoxins are shown to be consistent with Poisson-Boltzmann predictions of dielectric shielding that is relatively uniform for all of the amides that are sensitive to the metal charge, a region which corresponds to roughly 1/3 of the internal volume. The effective epsilon(int) values near 6 that are found in this study are significantly lower than many such estimates derived from sidechain pK measurements. The differing timeframes in which dielectric relaxation can respond to the highly transient peptide anion as compared to the longer lived states of the charged sidechains offers an explanation for the lower apparent dielectric constant deduced from these measurements.

Lessons in stability from thermophilic proteins

Studies that compare proteins from thermophilic and mesophilic organisms can provide insights into ability of thermophiles to function at their high habitat temperatures and may provide clues that enable us to better define the forces that stabilize all proteins. Most of the comparative studies have focused on thermal stability and show, as expected, that thermophilic proteins have higher Tm values than their mesophilic counterparts. Although these comparisons are useful, more detailed thermodynamic analyses are required to reach a more complete understanding of the mechanisms thermophilic protein employ to remain folded over a wider range of temperatures. This complete thermodynamic description allows one to generate a stability curve for a protein that defines how the conformational stability (DeltaG) varies with temperature. Here we compare stability curves for many pairs of homologous proteins from thermophilic and mesophilc organisms. Of the basic methods that can be employed to achieve enhanced thermostability, we find that most thermophilic proteins use the simple method that raises the DeltaG at all temperatures as the principal way to increase their Tm. We discuss and compare this thermodynamic method with the possible alternatives. In addition we propose ways that structural alterations and changes to the amino acid sequences might give rise to varied methods used to obtain thermostability.

Conformational dynamics of the tetrameric hemoglobin molecule as revealed by hydrogen exchange: I. Effects of pH, temperature, and ligand binding

Transcriptional regulation depends on the appropriate set of positive and negative signals that provide correct gene expression. Studies in eukaryotic gene expression have shown NC2 to be a general repressor of transcription which blocks the interaction between TFIIB and TBP. However, during the last few years NC2 has been found to bind the transcriptionally active promoters and interact with several positive transcription factors. These data suggested a controversial role of NC2 in transcription. Using in vitro transcription on minimal LTR promoter of HIV-1, we show that removal of NC2 from HeLa nuclear extract increases the yield of transcripts as well as unspecific transcription initiation in a template amount dependent manner. Fractions of HeLa nuclear extract containing NC2 can restore basal transcription repression and precise selection of transcription initiation point. This points to a new role for NC2 as a repressor of inaccurate transcription initiation that allows specific transcription to take place.

Slow irreversible unfolding of Pyrococcus furiosus triosephosphate isomerase: Separation and quantitation of conformers through a novel electrophoretic approach

The thermostability of hyperthermophile proteins is not easily studied because such proteins tend to be extremely recalcitrant to unfolding. Weeks of exposure to structurally destabilizing conditions are generally required to elicit any evidence of conformational change(s). The main reason for this extreme kinetic stability would appear to be the dominance of local unfolding transitions that occur within different parts of the structures of these molecules; put differently, local sub structural unfolding transitions that occur autonomously and reversibly are thought to fail to cooperate to bring about global unfolding in a facile manner, leading to a low overall observed rate of unfolding. For reasons that are not yet fully understood, unfolding is also reported to occur irreversibly in hyperthermophile proteins. Therefore, conventional experimental approaches are often unsuited to the study of their unfolding. Here, we describe a novel electrophoretic approach that facilitates separation, direct visualization, and quantitation of the folded, partially folded, and unfolded forms of the hyperthermophile protein triosephosphate isomerase from Pyrococcus furiosus, produced in the course of its irreversible structural destabilization by the combined action of heat and chemical agents. Our approach exploits (i) the irreversibility of global unfolding effected by heat and denaturants such as urea or guanidine hydrochloride, (ii) the stability of the native form of the protein to unfolding by the anionic detergent sodium dodecyl sulfate, (iii) the differential susceptibilities of various protein conformations to being bound by SDS, and (iv) the differential electrophoretic migration behavior displayed as a consequence of differential SDS binding.

Calorimetric and spectroscopic investigations of the thermal denaturation of wild type nitrite reductase

Nitrite reductase (NiR) is a multicopper protein, with a trimeric structure containing two types of copper site: type 1 is present in each subunit whereas type 2 is localized at the subunits interface. The paper reports on the thermal behaviour of wild type NiR from Alcaligenes faecalis S-6. The temperature-induced changes of the copper centres are characterized by optical spectroscopy and electron paramagnetic resonance spectroscopy, and by establishing the thermal stability by differential scanning calorimetry. The calorimetric profile of the enzyme shows a single endothermic peak with maximum heat absorption at T(m) approximately 100 degrees C, revealing an exceptional thermal stability. The thermal transition is irreversible and the scan rate dependence of the calorimetric trace indicates that the denaturation of NiR is kinetically controlled. The divergence of the activation energy values determined by different methods is used as a criterion for the inapplicability of the one-step irreversible model. The best fit of the DSC profiles is obtained when the classical Lumry-Eyring model, N<-->U-->F, is considered. The simulation results indicate that the irreversible step prevails on the reversible one. Moreover, it is found that the conformational changes within the type-1 copper environments precede the denaturation of the whole protein. No evidence of protein dissociation within the temperature range investigated was observed.

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.

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.

Increased peptide deformylase activity for N-formylmethionine processing of proteins overexpressed in Escherichia coli: application to homogeneous rubredoxin production

Deformylation of the initiator N-formylmethionine does not always proceed to completion for proteins overexpressed in Escherichia coli. To overcome this limitation, the def gene encoding the Escherichia coli peptide deformylase was cloned into the plysS plasmid under the tetracycline (Tc) promoter control. The efficiency of this constitutive level of peptide deformylase expression was demonstrated for the case of the rubredoxins from both mesophilic and hyperthermophilic organisms which normally retain a majority of their N-formyl terminal form. Indicating the potential structural/functional significance of residual formylation, the presence of a highly solvent exposed N-formyl group in rubredoxin is discernable in the amide NMR chemical shifts for the active site metal-coordinating cysteines more than 21A away.

Neutron crystallographic study on rubredoxin from Pyrococcus furiosus by BIX-3, a single-crystal diffractometer for biomacromolecules

The structure of a partially deuterated rubredoxin from the hyperthermophilic archaeon Pyrococcus furiosus, an organism that grows optimally at 100 degrees C, was determined by using the neutron single-crystal diffractometer dedicated for biological macromolecules (BIX-3) at the JRR-3M reactor of the Japan Atomic Energy Research Institute. Data were collected at room temperature up to a resolution of 1.5 A, and the completeness factor of the data set was 81.9%. The model contains 306 H and 50 D atoms. A total of 37 hydration water molecules were identified, with 15 having all three atoms fully located and the remaining D2O molecules partially defined. The model has been refined to final agreement factors of R = 18.6% and Rfree = 21.7%. Several orientations of the O-D bonds of side chains, whose assignments from x-ray data were previously ambiguous, were clearly visible in the neutron structure. Although most backbone N-H bonds had undergone some degree of H/D exchange throughout the rubredoxin molecule, 5 H atom positions still had distinctly negative (H) peaks. The neutron Fourier maps clearly showed the details of an extensive set of H bonds involving the ND3+ terminus that may contribute to the unusual thermostability of this molecule.

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

The striking kinetic stability of many proteins derived from hyperthermophilic organisms has led to the proposal that such stability may result from a heightened activation barrier for unfolding independent of a corresponding increase in the thermodynamic stability. This in turn implies a corresponding retardation of the folding reaction. A commonly cited model for kinetic thermal stabilization is the rubredoxin from Pyrococcus furiosus (Pf), which exhibits an irreversible denaturation lifetime at 100 degrees C of nearly a week. Utilizing protein resonances shifted well outside of the random coil chemical shift envelope, nuclear magnetic resonance (NMR) chemical exchange measurements on Pf rubredoxin as well as on the mesophile Clostridium pasteurianum (Cp) rubredoxin demonstrate reversible thermal transition temperatures of 144 degrees C (137 degrees C for the N-terminal modified A2K variant) and 104 degrees C, respectively, with similar (un)folding rates of approximately 25,000 s(-1), only modestly slower than the diffusion controlled rate. The absence of a substantial activation barrier to rubredoxin folding as well as the similar folding kinetics of the mesophile protein indicate that kinetic stabilization has not been utilized by the hyperthermophile rubredoxin in achieving its extreme thermal stability. The two-state folding kinetics observed for Pf rubredoxin contradict a previous assertion of multiphasic folding based on hydrogen exchange data extrapolated to an estimated midpoint of transition temperature (T(m)) of nearly 200 degrees C. This discrepancy is resolved by the observation that the base-catalyzed hydrogen exchange of the model dipeptide (N-acetyl-L-cysteine-N-methylamide)4-Cd2+ is 23-fold slower than that of the free cysteine model dipeptide used to normalize the Pf rubredoxin hydrogen exchange data.

Contribution of the [FeII(SCys)4] site to the thermostability of rubredoxins

The thermostabilities of Fe(2+) ligation in rubredoxins (Rds) from the hyperthermophile Pyrococcus furiosus (Pf) and the mesophiles Clostridium pasteurianum (Cp) and Desulfovibrio vulgaris (Dv) were compared. Residue 44 forms an NH.S(Cys) hydrogen bond to one of the cysteine ligands to the [Fe(SCys)(4)] site, and substitutions at this location affect the redox properties of the [Fe(SCys)(4)] site. Both Pf Rd and Dv Rd have an alanine residue at position 44, whereas Cp Fd has a valine residue. Wild-type proteins were examined along with V44A and A44V "exchange" mutants of Cp and Pf Rds, respectively, in order to assess the effects of the residue at position 44 on the stability of the [Fe(SCys)(4)] site. Stability of iron ligation was measured by temperature-ramp and fixed-temperature time course experiments, monitoring iron release in both the absence and presence of more thiophilic metals (Zn(2+), Cd(2+)) and over a range of pH values. The thermostability of the polypeptide fold was concomitantly measured by fluorescence, circular dichroism, and (1)H NMR spectroscopies. The A44V mutation strongly lowered the stability of the [Fe(II)(SCys)(4)] site in Pf Rd, whereas the converse V44A mutation of Cp Rd significantly raised the stability of the [Fe(II)(SCys)(4)] site, but not to the levels measured for wild-type Dv Rd. The region around residue 44 is thus a significant contributor to stability of iron coordination in reduced Rds. This region, however, made only a minor contribution to the thermostability of the protein folding, which was found to be higher for hyperthermophilic versus mesophilic Rds, and largely independent of the residue at position 44. These results, together with our previous studies, show that localized charge density, solvent accessibility, and iron site/backbone interactions control the thermostability of the [Fe(SCys)(4)] site. The iron site thermostability does make a minor contribution to the overall Rd thermostability. From a mechanistic standpoint, we also found that attack of displacing ions (H(+), Cd(2+)) on the Cys42 sulfur ligand at the [Fe(SCys)(4)] site occurs through the V8 side and not the V44 side of the iron site.

Protein stabilization by compatible solutes

Heteronuclear NMR relaxation measurements and hydrogen exchange data have been used to characterize protein dynamics in the presence or absence of stabilizing solutes from hyperthermophiles. Rubredoxin from Desulfovibrio gigas was selected as a model protein and the effect of diglycerol phosphate on its dynamic behaviour was studied. The presence of 100 mM diglycerol phosphate induces a fourfold increase in the half-life for thermal denaturation of D. gigas rubredoxin. A model-free analysis of the protein backbone relaxation parameters shows an average increase of generalized order parameters of 0.015 reflecting a small overall reduction in mobility of fast-scale motions. Hydrogen exchange data acquired over a temperature span of 20 degrees C yielded thermodynamic parameters for the structural opening reactions that allow for the exchange. This shows that the closed form of the protein is stabilized by an additional 1.6 kJ x mol(-1) in the presence of the solute. The results seem to indicate that the stabilizing effect is due mainly to a reduction in mobility of the slower, larger-scale motions within the protein structure with an associated increase in the enthalpy of interactions.

Protein hydrogen exchange mechanism: local fluctuations

Experiments were done to study the dynamic structural motions that determine protein hydrogen exchange (HX) behavior. The replacement of a solvent-exposed lysine residue with glycine (Lys8Gly) in a helix of recombinant cytochrome c does not perturb the native structure, but it entropically potentiates main-chain flexibility and thus can promote local distortional motions and large-scale unfolding. The mutation accelerates amide hydrogen exchange of the mutated residue by about 50-fold, neighboring residues in the same helix by less, and residues elsewhere in the protein not at all, except for Leu98, which registers the change in global stability. The pattern of HX changes shows that the coupled structural distortions that dominate exchange can be several residues in extent, but they expose to exchange only one amide NH at a time. This "local fluctuation" mode of hydrogen exchange may be generally recognized by disparate near-neighbor rates and a low dependence on destabilants (denaturant, temperature, pressure). In contrast, concerted unfolding reactions expose multiple neighboring amide NHs with very similar computed protection factors, and they show marked destabilant sensitivity. In both modes, ionic hydrogen exchange catalysts attack from the bulk solvent without diffusing through the protein matrix.

Submolecular cooperativity produces multi-state protein unfolding and refolding

Hydrogen exchange experiments show that cytochrome c and other proteins under native conditions reversibly unfold in a multi-step manner. The step from one intermediate to the next is determined by the intrinsically cooperative nature of secondary structural elements, which is retained in the native protein. Folding uses the same pathway in the reverse direction, moving from the unfolded to the native state through relatively discrete intermediate forms by the sequential addition of native-like secondary structural units.

Thermolysin-linearized microcin J25 retains the structured core of the native macrocyclic peptide and displays antimicrobial activity

Microcin J25 (MccJ25) is the single macrocyclic antimicrobial peptide belonging to the ribosomally synthesized class of microcins that are secreted by Enterobacteriaceae. It showed potent antibacterial activity against several Salmonella and Escherichia strains and exhibited a compact three-dimensional structure [Blond et al. (2001) Eur. J. Biochem., 268, 2124-2133]. The molecular mechanisms involved in the biosynthesis, folding and mode of action of MccJ25 are still unknown. We have investigated the structure and the antimicrobial activity of thermolysin-linearized MccJ25 (MccJ25-L1-21: VGIGTPISFY10GGGAGHVPEY20F), as well as two synthetic analogs, sMccJ25-L1-21 (sequence of the thermolysin-cleaved MccJ25) and sMccJ25-L12-11 (C-terminal sequence of the MccJ25 precursor: G12GAGHVPEYF21V1GIGTPISFYG11). The three-dimensional solution structure of MccJ25-L1-21, determined by two-dimensional NMR, consists of a boot-shaped hairpin-like well-defined 8-19 region flanked by disordered N and C termini. This structure is remarkably similar to that of cyclic MccJ25, and includes a short double-stranded antiparallel beta-sheet (8-10/17-19) perpendicular to a loop (Gly11-His16). The thermolysin-linearized MccJ25-L1-21 had antibacterial activity against E. coli and S. enteritidis strains, while both synthetic analogues lacked activity and organized structure. We show that the 8-10/17-19 beta-sheet, as well as the Gly11-His16 loop are required for moderate antibacterial activity and that the Phe21-Pro6 loop and the MccJ25 macrocyclic backbone are necessary for complete antibacterial activity. We also reveal a highly stable 8-19 structured core present in both the native MccJ25 and the thermolysin-linearized peptide, which is maintained under thermolysin treatment and resists highly denaturing conditions.

Exceptional stability of a [2Fe-2S] ferredoxin from hyperthermophilic bacterium Aquifex aeolicus

Aquifex aeolicus is the only hyperthermophile that is known to contain a plant- and mammalian-type [2Fe-2S] ferredoxin (Aae Fd1). This unique protein contains two cysteines, in addition to the four that act as ligands of the [2Fe-2S] cluster, which form a disulfide bridge. We have investigated the stability of Aae Fd1 with (wild-type) and without (C87A variant) the disulfide bond, with respect to pH, thermal and chemical perturbation, and compared the results to those for the mesophilic [2Fe-2S] ferredoxin from spinach. Unfolding reactions of all three proteins are irreversible due to cluster decomposition in the unfolded state. Wild-type and C87A Aae Fd1 proteins are extremely stable: unfolding at 20 degrees C requires high concentrations of the chemical denaturant and long incubation times. Moreover, their thermal-unfolding midpoints are 40-50 degrees higher than that for spinach ferredoxin (pH 7). The stability of the Aae Fd1 protein is significantly lower at pH 2.5 than pH 7 and 10, suggesting that ionic interactions play a role in structural integrity. Interestingly, the iron-sulfur cluster in C87A Aae Fd1 rearranges into a transient species with absorption bands at 520 and 610 nm, presumably a linear three-iron cluster, in the high-pH unfolded state.

Cytochrome c folding pathway: kinetic native-state hydrogen exchange

Native-state hydrogen exchange experiments under EX1 conditions can distinguish partially unfolded intermediates by their formation rates and identify the amide hydrogens exposed and protected in each. Results obtained define a cytochrome c intermediate seen only poorly before and place it early on the major unfolding pathway. Four distinct unfolding steps are found to be kinetically ordered in the same pathway sequence inferred before.

The temperature dependence of the hydrogen exchange in the SH3 domain of alpha-spectrin

The amide hydrogen-deuterium exchange (HX) in the Src homology region 3 (SH3) domain of alpha-spectrin has been measured by nuclear magnetic resonance as a function of temperature between 8 and 46 degrees C. The analysis of the temperature dependence of HX from a statistical thermodynamic point of view has allowed us to estimate the enthalpies and entropies of the conformational processes leading to HX. The results indicate that under native conditions the domain undergoes a wide variety of conformational fluctuations, ranging from local motions, mainly located in loops, turns and chain ends and involving only low enthalpy and entropy, to extensive structural disruptions affecting its core and involving enthalpies and entropies that come fairly close to those observed during global unfolding.

O2 penetration and proton burial depth in proteins: applicability to fold family recognition

Paramagnetically induced relaxation effects of O2 and the nitroxide 4-hydroxy TEMPO were measured for the amide protons of perdeuterated rubredoxin from the hyperthermophilic archaeon Pyrococcus furiosus and the mesophilic bacterium Clostridium pasteurianum. For both O2 and the impermeant nitroxide, the induced relaxation at the static solvent inaccessible amide sites is dominated by long-range interactions with the paramagnetic species in the bulk aqueous phase. The upper bound of O2 solubility in the internal matrix of the rubredoxins is one-tenth that of the bulk aqueous phase. Furthermore, the difference between the oxygen solubilities inside the two rubredoxins is at most 1% that of bulk water O2 solubility, suggesting that there are only modest differences in this measure of fluidity for the mesophile vs hyperthermophile protein interiors. Calculations based on the assumption of a paramagnet uniformly distributed on the protein exterior yield accurate predictions at nearly all amide sites for the minimum relaxation value observed from either the O2 or nitroxide data. Model calculations indicate that the readily obtained paramagnetically induced relaxation effects should prove effective in recognition of structural homology for proteins that are too widely diverged for sequence-based recognition.

pH dependence of the hydrogen exchange in the SH3 domain of alpha-spectrin

Using nuclear magnetic resonance we have measured the hydrogen exchange (HX) in the Src homology region 3 (SH3) domain of alpha-spectrin as a function of pH*. At very acidic pH* values the exchange of most residues appears to occur via global unfolding, although several residues show abnormally large Gibbs energies of exchange, suggesting the presence of some residual structure in the unfolded state. At higher pH* HX occurs mainly via local or partial unfoldings. We have been able to characterize the coupling between the electrostatic interactions in this domain and the conformational fluctuations occurring under native conditions by analyzing the dependence upon pH* of the Gibbs energy of exchange. The SH3 domain seems to be composed of a central core, which requires large structural disruptions to become exposed to the solvent, surrounded by smaller subdomains, which fluctuate independently.

Femtosecond dynamics of rubredoxin: tryptophan solvation and resonance energy transfer in the protein

We report here studies of tryptophan (Trp) solvation dynamics in water and in the Pyrococcus furiosus rubredoxin protein, including the native and its apo and denatured forms. We also report results on energy transfer from Trp to the iron-sulfur [Fe-S] cluster. Trp fluorescence decay with the onset of solvation dynamics of the chromophore in water was observed with femtosecond resolution ( approximately 160 fs; 65% component), but the emission extended to the picosecond range (1.1 ps; 35% component). In contrast, the decay is much slower in the native rubredoxin; the Trp fluorescence decay extends to 10 ps and longer, reflecting the local rigidity imposed by residues and by the surface water layer. The dynamics of resonance energy transfer from the two Trps to the [Fe-S] cluster in the protein was observed to follow a temporal behavior characterized by a single exponential (15-20 ps) decay. This unusual observation in a protein indicates that the resonance transfer is to an acceptor of a well-defined orientation and separation. From studies of the mutant protein, we show that the two Trp residues have similar energy-transfer rates. The critical distance for transfer (R(0)) was determined, by using the known x-ray data, to be 19.5 A for Trp-36 and 25.2 A for Trp-3, respectively. The orientation factor (kappa(2)) was deduced to be 0.13 for Trp-36, clearly indicating that molecular orientation of chromophores in the protein cannot be isotropic with kappa(2) value of 2/3. These studies of solvation and energy-transfer dynamics, and of the rotational anisotropy, of the wild-type protein, the (W3Y, I23V, L32I) mutant, and the fmetPfRd variant at various pH values reveal a dynamically rigid protein structure, which is probably related to the hyperthermophilicity of the protein.

High thermal stability of 3-isopropylmalate dehydrogenase from Thermus thermophilus resulting from low DeltaC(p) of unfolding

To characterize the thermal stability of 3-isopropylmalate dehydrogenase (IPMDH) from an extreme thermophile, Thermus thermophilus, urea-induced unfolding of the enzyme and of its mesophilic counterpart from Escherichia coli was investigated at various temperatures. The unfolding curves were analyzed with a three-state model for E.coli IPMDH and with a two-state model for T.thermophilus IPMDH, to obtain the free energy change DeltaG degrees of each unfolding process. Other thermodynamic parameters, enthalpy change DeltaH, entropy change DeltaS and heat capacity change DeltaC(p), were derived from the temperature dependence of DeltaG degrees. The main feature of the thermophilic enzyme was its lower dependence of DeltaG degrees on temperature resulting from a low DeltaC(p). The thermophilic IPMDH had a significantly lower DeltaC(p), 1.73 kcal/mol.K, than that of E.coli IPMDH (20.7 kcal/mol.K). The low DeltaC(p) of T.thermophilus IPMDH could not be predicted from its change in solvent-accessible surface area DeltaASA. The results suggested that there is a large structural difference between the unfolded state of T.thermophilus and that of E.coli IPMDH. Another responsible factor for the higher thermal stability of T.thermophilus IPMDH was the increase in the most stable temperature T(s). The DeltaG degrees maximum of T.thermophilus IPMDH was much smaller than that of E.coli IPMDH. The present results clearly demonstrated that a higher melting temperature T(m) is not necessarily accompanied by a higher DeltaG degrees maximum.

On the nature of conformational openings: native and unfolded-state hydrogen and thiol-disulfide exchange studies of ferric aquomyoglobin

Native-state amide hydrogen exchange (HX) of proteins in the presence of denaturant has provided valuable details on the structures of equilibrium folding intermediates. Here, we extend HX theory to model thiol group exchange (SX) in single cysteine-containing variants of sperm whale ferric aquomyoglobin. SX is complementary to HX in that it monitors conformational opening events that expose side-chains, rather than the main chain, to solvent. A simple two-process model, consisting of EX2-limited local structural fluctuations and EX1-limited global unfolding, adequately accounts for all HX data. SX is described by the same model except at very low denaturant concentrations and when the bulky labeling reagent 5,5'-dithiobis (2-nitrobenzoic acid) (DTNB) is used. Under these conditions SX can occur by a novel denaturant-dependent process. This anomalous behavior is not observed when the smaller labeling reagent methyl methanethiosulfonate is employed, suggesting that it reflects a denaturant-induced increase in the amplitudes of local structural fluctuations. It also is not seen in heme-free apomyoglobin, which may indicate that local openings are sufficiently large in the absence of denaturant to allow DTNB unhindered access. Differences in SX kinetics obtained using the two labeling reagents provide estimates of the sizes of local opening reactions at different sites in the protein. At all sequence positions examined except for position 73, the same opening event appears to facilitate exchange of both backbone amide and side-chain thiol groups. The C73 thiol group is exposed by a low-energy fluctuation that does not expose its amide group to exchange.