Folding Mechanism of Proteins Im7 and Im9: Insight from All-Atom Simulations in Implicit and Explicit Solvent.

Im7 and Im9 are evolutionary related proteins with almost identical native structures. In spite of their structural similarity, experiments show that Im7 folds through a long-lived on-pathway intermediate, while Im9 folds according to two-state kinetics. In this work, we use a recently developed enhanced path sampling method to generate many folding trajectories for these proteins, using realistic atomistic force fields, in both implicit and explicit solvent. Overall, our results are in good agreement with the experimental ϕ values and with the result of ϕ-value-restrained molecular dynamics (MD) simulations. However, our implicit solvent simulations fail to predict a qualitative difference in the folding pathways of Im7 and Im9. In contrast, our simulations in explicit solvent correctly reproduce the fact that only protein Im7 folds through a on-pathway intermediate. By analyzing our atomistic trajectories, we provide a physical picture which explains the observed difference in the folding kinetics of these chains.

Patterns of adaptation in a laboratory evolved thermophilic enzyme.

The heat sensitive psychrophilic protease subtilisin S41 was previously subjected to three rounds of mutagenesis/recombination and screening, resulting in variant 3-2G7, whose half-life at 60 degrees C is approx. 500 times that of wild-type. Here we report the results of five additional generations of laboratory evolution starting from 3-2G7. The half-life of 8th generation enzyme 8-4A9 at 60 degrees C is 1200 times that of wild-type, and slightly more than twice that of 3-2G7. This half-life is >20-fold greater than those of homologous mesophilic subtilisins SSII and BPN'. Circular dichroism melting curves indicate that subtilisin 8-4A9 unfolds at temperatures approx. 25 degrees C higher than wild-type. It is also substantially more resistant to proteolysis at 30 degrees C. Nearly half of the 13 amino acid substitutions accumulated in 8-4A9 involve the mutation of serine residues. This mirrors a pattern observed in natural proteins, where serines are statistically less prevalent in thermophilic enzymes compared to mesophilic ones.

How enzymes adapt: lessons from directed evolution.

Enzymes that are adapted to widely different temperature niches are being used to investigate the molecular basis of protein stability and enzyme function. However, natural evolution is complex: random noise, historical accidents and ignorance of the selection pressures at work during adaptation all cloud comparative studies. Here, we review how adaptation in the laboratory by directed evolution can complement studies of natural enzymes in the effort to understand stability and function. Laboratory evolution experiments can attempt to mimic natural evolution and identify different adaptive mechanisms. However, laboratory evolution might make its biggest contribution in explorations of nonnatural functions, by allowing us to distinguish the properties nutured by evolution from those dictated by the laws of physical chemistry.

Cold adaptation of a mesophilic subtilisin-like protease by laboratory evolution.

Enzymes isolated from organisms native to cold environments generally exhibit higher catalytic efficiency at low temperatures and greater thermosensitivity than their mesophilic counterparts. In an effort to understand the evolutionary process and the molecular basis of cold adaptation, we have used directed evolution to convert a mesophilic subtilisin-like protease from Bacillus sphaericus, SSII, into its psychrophilic counterpart. A single round of random mutagenesis followed by recombination of improved variants yielded a mutant, P3C9, with a catalytic rate constant (k(cat)) at 10 degrees C 6.6 times and a catalytic efficiency (k(cat)/K(M)) 9.6 times that of wild type. Its half-life at 70 degrees C is 3.3 times less than wild type. Although there is a trend toward decreasing stability during the progression from mesophile to psychrophile, there is not a strict correlation between decreasing stability and increasing low temperature activity. A first generation mutant with a >2-fold increase in k(cat) is actually more stable than wild type. This suggests that the ultimate decrease in stability may be due to random drift rather than a physical incompatibility between low temperature activity and high temperature stability. SSII shares 77. 4% identity with the naturally psychrophilic protease subtilisin S41. Although SSII and S41 differ at 85 positions, four amino acid substitutions were sufficient to generate an SSII whose low temperature activity is greater than that of S41. That none of the four are found in S41 indicates that there are multiple routes to cold adaptation.

Directed evolution study of temperature adaptation in a psychrophilic enzyme.

We have used laboratory evolution methods to enhance the thermostability and activity of the psychrophilic protease subtilisin S41, with the goal of investigating the mechanisms by which this enzyme can adapt to different selection pressures. A combined strategy of random mutagenesis, saturation mutagenesis and in vitro recombination (DNA shuffling) was used to generate mutant libraries, which were screened to identify enzymes that acquired greater thermostability without sacrificing low-temperature activity. The half-life of seven-amino acid substitution variant 3-2G7 at 60 degrees C is approximately 500 times that of wild-type and far surpasses those of homologous mesophilic subtilisins. The dependence of half-life on calcium concentration indicates that enhanced calcium binding is largely responsible for the increased stability. The temperature optimum of the activity of 3-2G7 is shifted upward by approximately 10 degrees C. Unlike natural thermophilic enzymes, however, the activity of 3-2G7 at low temperatures was not compromised. The catalytic efficiency, k(cat)/K(M), was enhanced approximately threefold over a wide temperature range (10 to 60 degrees C). The activation energy for catalysis, determined by the temperature dependence of k(cat)/K(M) in the range 15 to 35 degrees C, is nearly identical to wild-type and close to half that of its highly similar mesophilic homolog, subtilisin SSII, indicating that the evolved S41 enzyme retained its psychrophilic character in spite of its dramatically increased thermostability. These results demonstrate that it is possible to increase activity at low temperatures and stability at high temperatures simultaneously. The fact that enzymes displaying both properties are not found in nature most likely reflects the effects of evolution, rather than any intrinsic physical-chemical limitations on proteins.

Energetics of target peptide recognition by calmodulin: a calorimetric study.

Calmodulin is a small protein involved in the regulation of a wide variety of intracellular processes. The cooperative binding of Ca2+ to calmodulin's two Ca2+ binding domains induces conformational changes which allow calmodulin to activate specific target enzymes. The association of calmodulin with a peptide corresponding to the calmodulin binding site of rabbit smooth muscle myosin light chain kinase (smMLCKp) was studied using isothermal titration microcalorimetry. The dependence of the binding energetics on temperature, pH, Ca2+ concentration, and NaCl concentration were determined. It is found that the binding of calmodulin to smMLCKp proceeds with negative changes in enthalpy (deltaH), heat capacity (deltaCp), and entropy (deltaS) near room temperature, indicating that it is an enthalpically driven process that is entropically unfavorable. From these results it is concluded that the hydrophobic effect, an entropic effect which favors the removal of non-polar protein groups from water, is not a major driving force in calmodulin-smMLCKp recognition. Although a large number of non-polar side-chains are buried upon binding, these stabilize the complex primarily by forming tightly packed van der Waals interactions with one another. Binding at acidic pH was studied in order to assess the contribution of electrostatic interactions to binding. It is found that moving to acidic pH results in a large decrease in the Gibbs free energy of binding but no change in the enthalpy, indicating that electrostatic interactions contribute only entropically to the binding energetics. The accessible surface area and atomic packing density of the calmodulin-smMLCKp crystal structure are analyzed, and the results discussed in relation to the experimental data.

Structural energetics of barstar studied by differential scanning microcalorimetry.

The energetics of barstar denaturation have been studied by CD and scanning microcalorimetry in an extended range of pH and salt concentration. It was shown that, upon increasing temperature, barstar undergoes a transition to the denatured state that is well approximated by a two-state transition in solutions of high ionic strength. This transition is accompanied by significant heat absorption and an increase in heat capacity. The denaturational heat capacity increment at approximately 75 degrees C was found to be 5.6 +/- 0.3 kJ K-1 mol-1. In all cases, the value of the measured enthalpy of denaturation was notably lower than those observed for other small globular proteins. In order to explain this observation, the relative contributions of hydration and the disruption of internal interactions to the total enthalpy and entropy of unfolding were calculated. The enthalpy and entropy of hydration were found to be in good agreement with those calculated for other proteins, but the enthalpy and entropy of breaking internal interactions were found to be among the lowest for all globular proteins that have been studied. Additionally, the partial specific heat capacity of barstar in the native state was found to be 0.37 +/- 0.03 cal K-1 g-1, which is higher than what is observed for most globular proteins and suggests significant flexibility in the native state. It is known from structural data that barstar undergoes a conformational change upon binding to its natural substrate barnase.(ABSTRACT TRUNCATED AT 250 WORDS)

The encephalomyocarditis virus 3C protease is a substrate for the ubiquitin-mediated proteolytic system.

The encephalomyocarditis virus 3C protease has been shown to be rapidly degraded in infected cells and in vitro in rabbit reticulocyte lysate. The in vitro degradation, at least, is accomplished by a virus-independent, ATP-dependent proteolytic system. Here we identify this proteolytic system as the ubiquitin-mediated system. Incubation of the 3C protease in rabbit reticulocyte or cultured mouse cell lysate preparations, alone or in the presence of added ubiquitin or methylated ubiquitin, resulted in the generation of new higher molecular weight species. These new products were shown to be 3C protease-ubiquitin conjugates by their ability to bind antibodies against both the 3C protease and ubiquitin. Supplemental ubiquitin also stimulated the degradation of the 3C protease in these preparations. Large 3C protease-polyubiquitin conjugates were observed to accumulate in reticulocyte lysate in the presence of adenosine 5'-O-(3-thiotriphosphate), an inhibitor of the 26 S multicatalytic protease. This, combined with the fact that the proteolytic activity could be removed from the lysate by sedimentation, implicates the multicatalytic protease in the degradation of the 3C protease-ubiquitin conjugates. It was also found that the slow rate of degradation of a model polyprotein, which resembles the stable viral 3CD diprotein produced in vivo, is likely due to the fact that the polyprotein is a poor substrate for the ubiquitin-conjugating system.

Thermodynamics of ubiquitin unfolding.

The energetics of ubiquitin unfolding have been studied using differential scanning microcalorimetry. For the first time it has been shown directly that the enthalpy of protein unfolding is a nonlinear function of temperature. Thermodynamic parameters of ubiquitin unfolding were correlated with the structure of the protein. The enthalpy of hydrogen bonding in ubiquitin was calculated and compared to that obtained for other proteins. It appears that the energy of hydrogen bonding correlates with the average length of the hydrogen bond in a given protein structure.