Modeling the structure of Pyrococcus furiosus rubredoxin by homology to other X-ray structures

Dynamics and unfolding pathways of a hyperthermophilic and a mesophilic rubredoxin

Molecular dynamics simulations in solution are performed for a rubredoxin from the hyperthermophilic archaeon Pyrococcus furiosus (RdPf) and one from the mesophilic organism Desulfovibrio vulgaris (RdDv). The two proteins are simulated at four temperatures: 300 K, 373 K, 473 K (two sets), and 500 K; the various simulations extended from 200 ps to 1,020 ps. At room temperature, the two proteins are stable, remain close to the crystal structure, and exhibit similar dynamic behavior; the RMS residue fluctuations are slightly smaller in the hyperthermophilic protein. An analysis of the average energy contributions in the two proteins is made; the results suggest that the intraprotein energy stabilizes RdPf relative to RdDv. At 373 K, the mesophilic protein unfolds rapidly (it begins to unfold at 300 ps), whereas the hyperthermophilic does not unfold over the simulation of 600 ps. This is in accord with the expected stability of the two proteins. At 473 K, where both proteins are expected to be unstable, unfolding behavior is observed within 200 ps and the mesophilic protein unfolds faster than the hyperthermophilic one. At 500 K, both proteins unfold; the hyperthermophilic protein does so faster than the mesophilic protein. The unfolding behavior for the two proteins is found to be very similar. Although the exact order of events differs from one trajectory to another, both proteins unfold first by opening of the loop region to expose the hydrophobic core. This is followed by unzipping of the beta-sheet. The results obtained in the simulation are discussed in terms of the factors involved in flexibility and thermostability.

A model for the unusual kinetics of thermal denaturation of rubredoxin

The thermal denaturation of the simple, redox-active iron protein rubredoxin is characterized by a slow, irreversible decay of the characteristic red color of the iron center at elevated temperatures in the presence of oxygen at pH 7.8. The denaturation rate is essentially constant and the time period for complete bleaching is nearly independent of protein concentration. These two characteristics of the kinetics can be fit by a simple self-catalyzed kinetics model consisting of the combination of a first-order decay and catalysis by some product of that decay, i.e., dP/dt = k1[A] + (k2[P][A])/(K(m) + [A]), where A is native rubredoxin, P, is unspecified product, k1 is a first-order rate constant, and k2 and K(m) are the catalytic constants. In order for the second term to be of this simple form over the full course of a decay, the model must include the condition that the reaction is effectively irreversible. This model has properties which suggest other biological roles in regulation (changes in k1 or k2 can dramatically modulate the kinetics), in timing (titer-independent fixed reaction time), and in self-activation reactions. At one extreme (k1 >> k2) the kinetics becomes exponential, but at the other extreme (k2 >> k1) they show a dramatic and rapid terminal increase after a lag period. Some obvious possible roles in the kinetics of programmed cell death, prion disease, and protease autoactivation are discussed.

Biochemical diversity among sulfur-dependent, hyperthermophilic microorganisms

Hyperthermophiles are a recently discovered group of microorganisms that grow at and above 90 degrees C. They currently comprise over 20 different genera, and except for two novel bacteria, all are classified as Archaea. The majority of these organisms are obligately anaerobic heterotrophs that reduce elemental sulfur (S degree) to H2S. The best studied from a biochemical perspective are the archaeon, Pyrococcus furiosus, and the bacterium, Thermotoga maritima, both of which are saccharolytic. P. furiosus is thought to contain a new type of Entner-Doudoroff pathway for the conversion of carbohydrates ultimately to acetate, H2 and CO2. The pathway is independent of nicotinamide nucleotides and involves novel types of ferredoxin-linked oxidoreductases, one of which has tungsten, a rarely used element, as a prosthetic group. The only site of energy conservation is at the level of acetyl CoA, which is the presence of ADP and phosphate is converted to acetate and ATP in a single step. In contrast, T. maritima utilizes a conventional Embden-Meyerhof pathway for sugar oxidation. P. furiosus also utilizes peptides as a sole carbon and energy source. Amino acid oxidation is thought to involve glutamate dehydrogenase together with at least three types of novel ferredoxin-linked oxidoreductases which catalyze the oxidation of 2-ketoglutarate, aryl pyruvates and formaldehyde. One of these enzymes also utilizes tungsten. In P. furiosus, virtually all of the reductant that is generated during the catabolism of both carbohydrates and peptides is channeled to a cytoplasmic hydrogenase. This enzyme is now termed sulhydrogenase, as it reduces both protons to H2 and S degrees (or polysulfide) to H2S. S degrees reduction appears to lead to the conservation of energy in P. furiosus but not in T. maritima, although the mechanism by which this occurs is not known.

Thermostability and thermoactivity of enzymes from hyperthermophilic Archaea

Enzymes from hyperthermophilic microorganisms are characteristically thermostable and thermoactive at extremely high temperatures. Information about the basis for the structure and function of these novel proteins is beginning to emerge. However, there are very few generalizations that can be drawn at this point that can be derived from the limited number of studies that have focused on biocatalysis and thermostability at extremely high temperatures.

Computational chemistry and molecular modeling of electron-transfer proteins

The methods of computational chemistry and molecular modeling are becoming more and more accessible to biochemists with the advent of fast, inexpensive graphics workstations and well-tested computer programs. The state of the art in small molecules allows chemists to use these programs as "black boxes" and be confident of the results at an amazingly high level of precision. This is not the case, however, for biological macro-molecules at this time. Therefore, it is necessary before using the programs listed in Section I that we familiarize ourselves with their theoretical basis and limitations. It is also important that they be used in the context of the accumulated literature on their use. The survey given in this chapter is intended as an introduction to these tools and as a source for initiating the discovery process with the literature cited.

Knowledge-based protein modeling

Knowledge, both from the three-dimensional structures of homologous proteins and from the general analysis of protein structure, is of value in modeling a protein of known sequence but unknown structure. While many models are still constructed at least in part by manual methods on graphics devices, automated procedures have come into greater use. These procedures include those that assemble fragments of structure from other known structures and those that derive coordinates for the model from the satisfaction of restraints placed on atomic positions.

Investigations of the thermostability of rubredoxin models using molecular dynamics simulations

The affects of differences in amino acid sequence on the temperature stability of the three-dimensional structure of the small beta-sheet protein, rubredoxin (Rd), was revealed when a set of homology models was subjected to molecular dynamics simulations at relatively high temperatures. Models of Rd from the hyperthermophile, Pyrococcus furiosus (Pf), an organism that grows optimally at 100 degrees C, were compared to three mesophilic Rds of known X-ray crystal structure. Simulations covering the limits of known Rd thermostabilities were carried out at temperatures of 300 K, 343 K, 373 K, and 413 K. They suggest that Rd stability is correlated with structural dynamics. Because the dynamic behavior of three Pf Rd models was consistently different from the dynamic behavior of the three mesophilic Rd structures, detailed analysis of the temperature-dependent dynamic behavior was carried out. The major differences between the models of the protein from the hyperthermophile and the others were: (1) an obvious temperature-dependent transition in the mesophilic structures not seen with the Pf Rd models, (2) consistent AMBER energy for the Pf Rd due to differences in nonbonded interaction terms, (3) less variation in the average conformations for the Pf Rd models with temperature, and (4) the presence of more extensive secondary structure for the Pf Rd models. These unsolvated dynamics simulations support a simple, general hypothesis to explain the hyperthermostability of Pf Rd. Its structure simplifies the conformational space to give a single minimum accessible over an extreme range of temperatures, whereas the mesophilic proteins sample a more complex conformational space with two or more minima over the same temperature range.