Non-linear effects of temperature and urea on the thermodynamics and kinetics of folding and unfolding of hisactophilin

Extensive measurements and analysis of thermodynamic stability and kinetics of urea-induced unfolding and folding of hisactophilin are reported for 5-50 degrees C, at pH 6.7. Under these conditions hisactophilin has moderate thermodynamic stability, and equilibrium and kinetic data are well fit by a two-state transition between the native and the denatured states. Equilibrium and kinetic m values decrease with increasing temperature, and decrease with increasing denaturant concentration. The betaF values at different temperatures and urea concentrations are quite constant, however, at about 0.7. This suggests that the transition state for hisactophilin unfolding is native-like and changes little with changing solution conditions, consistent with a narrow free energy profile for the transition state. The activation enthalpy and entropy of unfolding are unusually low for hisactophilin, as is also the case for the corresponding equilibrium parameters. Conventional Arrhenius and Eyring plots for both folding and unfolding are markedly non-linear, but these plots become linear for constant DeltaG/T contours. The Gibbs free energy changes for structural changes in hisactophilin have a non-linear denaturant dependence that is comparable to non-linearities observed for many other proteins. These non-linearities can be fit for many proteins using a variation of the Tanford model, incorporating empirical quadratic denaturant dependencies for Gibbs free energies of transfer of amino acid constituents from water to urea, and changes in fractional solvent accessible surface area of protein constituents based on the known protein structures. Noteworthy exceptions that are not well fit include amyloidogenic proteins and large proteins, which may form intermediates. The model is easily implemented and should be widely applicable to analysis of urea-induced structural transitions in proteins.     

Putative interhelix ion pairs involved in the stability of myoglobin.

An earlier theoretical study predicted that specific ion pair interactions between neighboring helices should be important in stabilizing myoglobin. To measure these interactions in sperm whale myoglobin, single mutations were made to disrupt them. To obtain reliable DeltaG values, conditions were found in which the urea induced unfolding of holomyoglobin is reversible and two-state. The cyanomet form of myoglobin satisfies this condition at pH 5, 25 degrees C. The unfolding curves monitored by far-UV CD and Soret absorbance are superimposable and reversible. None of the putative ion pairs studied here makes a large contribution to the stability of native myoglobin. The protein stability does decrease somewhat between 0 and 0.1 M NaCl, however, indicating that electrostatic interactions contribute favorably to myoglobin stability at pH 5.0. A previous mutational study indicated that the net positive charge of the A[B]GH subdomain of myoglobin is an important factor affecting the stability of the pH 4 folding intermediate and potential ion pairs within the subdomain do not contribute significantly to its stability. One of the assumptions made in that study is tested here: replacement of either positively or negatively charged residues outside the A[B]GH subdomain has no significant effect on the stability of the pH 4 molten globule.     

Role of the charge-charge interactions in defining stability and halophilicity of the CspB proteins

Charge-charge interactions on the surface of native proteins are important for protein stability and can be computationally redesigned in a rational way to modulate protein stability. Such computational effort led to an engineered protein, CspB-TB that has the same core as the mesophilic cold shock protein CspB-Bs from Bacillus subtilis, but optimized distribution of charge-charge interactions on the surface. The CspB-TB protein shows an increase in the transition temperature by 20 degrees C relative to the unfolding temperature of CspB-Bs. The CspB-TB and CspB-Bs protein pair offers a unique opportunity to further explore the energetics of charge-charge interactions as the substitutions at the same sequence positions are done in largely similar structural but different electrostatic environments. In particular we addressed two questions. What is the contribution of charge-charge interactions in the unfolded state to the protein stability and how amino acid substitutions modulate the effect of increase in ionic strength on protein stability (i.e. protein halophilicity). To this end, we experimentally measured the stabilities of over 100 variants of CspB-TB and CspB-Bs proteins with substitutions at charged residues. We also performed computational modeling of these protein variants. Analysis of the experimental and computational data allowed us to conclude that the charge-charge interactions in the unfolded state of two model proteins CspB-Bs and CspB-TB are not very significant and computational models that are based only on the native state structure can adequately, i.e. qualitatively (stabilizing versus destabilizing) and semi-quantitatively (relative rank order), predict the effects of surface charge neutralization or reversal on protein stability. We also show that the effect of ionic strength on protein stability (protein halophilicity) appears to be mainly due to the screening of the long-range charge-charge interactions.     

Use of site-directed mutagenesis to destabilize native apomyoglobin relative to folding intermediates

Site-directed mutagenesis has been used to study the effect on the stability of human apomyoglobin (apoMb) of modifying the size, hydrophobicity, and charge of a central residue in the G.B helix-helix packing interface. Some stability measurements have also been made on the corresponding holomyoglobins (heme present). Cys-110, a central helix pairing residue in the G helix, has been changed to Ala, Ser, Asp, and Leu. Stability to low-pH-induced unfolding has been measured for both native apoMb and the compact folding intermediate discovered by Griko et al. [Griko, Y. V., Privalov, P. L., Venyaminov, S. Y., & Kutyshenko, V. P. (1988) J. Mol. Biol. 202, 127-138]. As judged by its circular dichroism spectrum, this intermediate has a substantial helix content (about 35%). Whether or not this inferred helical structure is closely related to the myoglobin structure is not yet known. The mutational evidence shows that integrity of G.B helix pairing is important for the stability of apoMb as well as of myoglobin and that this helix pairing site is very sensitive to both steric and electrostatic disruption. Our results also suggest that G.B helix pairing does not stabilize the compact intermediate; hence, disrupting this site destabilizes the native protein relative to the compact intermediate. Such selective destabilization of the native state relative to equilibrium folding intermediates is not restricted to acid denaturation: urea denaturation of the Leu mutant appears to display at least one stable intermediate, while wild-type and the remaining mutant apoMbs undergo two-state urea unfolding transitions.     

Rapid unfolding of a domain populates an aggregation-prone intermediate that can be recognized by GroEL

Some amino acid substitutions in phage P22 coat protein cause a temperature-sensitive folding (tsf) phenotype. In vivo, these tsf amino acid substitutions cause coat protein to aggregate and form intracellular inclusion bodies when folded at high temperatures, but at low temperatures the proteins fold properly. Here the effects of tsf amino acid substitutions on folding and unfolding kinetics and the stability of coat protein in vitro have been investigated to determine how the substitutions change the ability of coat protein to fold properly. The equilibrium unfolding transitions of the tsf variants were best fit to a three-state model, N if I if U, where all species concerned were monomeric, a result confirmed by velocity sedimentation analytical ultracentrifugation. The primary effect of the tsf amino acid substitutions on the equilibrium unfolding pathway was to decrease the stability (DeltaG) and the solvent accessibility (m-value) of the N if I transition. The kinetics of folding and unfolding of the tsf coat proteins were investigated using tryptophan fluorescence and circular dichroism (CD) at 222 nm. The tsf amino acid substitutions increased the rate of unfolding by 8-14-fold, with little effect on the rate of folding, when monitored by tryptophan fluorescence. In contrast, when folding or unfolding reactions were monitored by CD, the reactions were too fast to be observed. The tsf coat proteins are natural substrates for the molecular chaperones, GroEL/S. When native tsf coat protein monomers were incubated with GroEL, they bound efficiently, indicating that a folding intermediate was significantly populated even without denaturant. Thus, the tsf coat proteins aggregate in vivo because of an increased propensity to populate this unfolding intermediate.     

Correlation between mutational destabilization of phage T4 lysozyme and increased unfolding rates

The thermodynamics and kinetics of unfolding of 28 bacteriophage T4 lysozyme variants were compared by using urea gradient gel electrophoresis. The mutations studied cause a variety of sequence changes at different residues throughout the polypeptide chain and result in a wide range of thermodynamic stabilities. A striking relationship was observed between the thermodynamic and kinetic effects of the amino acid replacements: All the substitutions that destabilized the native protein by 2 kcal/mol or more also increased the rate of unfolding. The observed increases in unfolding rate corresponded to a decrease in the activation energy of unfolding (delta Gu) at least 35% as large as the decrease in thermodynamic stability (delta Gu). Thus, the destabilizing lesions bring the free energy of the native state closer to that of both the unfolded state and the transition state for folding and unfolding. Since a large fraction of the mutational destabilization is expressed between the transition state and the native conformation, the changes in folding energetics cannot be accounted for by effects on the unfolded state alone. The results also suggest that interactions throughout much of the folded structure are altered in the formation of the transition state during unfolding.     

Changing Folding and Binding Stability in a Viral Coat Protein: A Comparison between Substitutions Accessible through Mutation and Those Fixed by Natural Selection

Previous studies have shown that most random amino acid substitutions destabilize protein folding (i.e. increase the folding free energy). No analogous studies have been carried out for protein-protein binding. Here we use a structure-based model of the major coat protein in a simple virus, bacteriophage φX174, to estimate the free energy of folding of a single coat protein and binding of five coat proteins within a pentameric unit. We confirm and extend previous work in finding that most accessible substitutions destabilize both protein folding and protein-protein binding. We compare the pool of accessible substitutions with those observed among the φX174-like wild phage and in experimental evolution with φX174. We find that observed substitutions have smaller effects on stability than expected by chance. An analysis of adaptations at high temperatures suggests that selection favors either substitutions with no effect on stability or those that simultaneously stabilize protein folding and slightly destabilize protein binding. We speculate that these mutations might involve adjusting the rate of capsid assembly. At normal laboratory temperature there is little evidence of directional selection. Finally, we show that cumulative changes in stability are highly variable; sometimes they are well beyond the bounds of single substitution changes and sometimes they are not. The variation leads us to conclude that phenotype selection acts on more than just stability. Instances of larger cumulative stability change (never via a single substitution despite their availability) lead us to conclude that selection views stability at a local, not a global, level.     

Kinetic stability of designed glycosylation mutants of Coprinus cinereus peroxidase

The effect of glycans and surface mutations on protein unfolding induced by heat or urea has been studied. Removal of the only native high mannose type glycan in the N142P, N142T, and N142D CIP mutants reduced the lifetime to half of that of wtCIP at irreversible conditions of unfolding. The effect was moderate at reversible conditions. Five glycomutants designed to have 0, 1, 2, 4 and 6N glycans showed a correlation between increased carbohydrate mass and increased stability toward irreversible unfolding. The results are in agreement with a dampening effect of glycans on backbone fluctuation in both the native and the unfolded states. However, experiments in reversible conditions were less clear because of additional effects of an increasing number of amino acid substitutions and aggregation. Examples of strong effects from minor surface changes were also observed.     

An improved experimental system for determining small folding entropy changes resulting from proline to alanine substitutions

Changes in protein stability can be achieved by making substitutions that increase or decrease the available conformations of the unfolded protein without altering the conformational freedom of the folded protein. Matthews and coworkers (1987) proposed that proline to alanine (P --> A) substitution would achieve this type of entropic destabilization. By comparing the Ramachandran area associated with alanine and proline residues, Matthews et al. estimated the unfolding entropy change resulting from P --> A substitution to be 4.8 cal mol(-1) K(-1). Although such an entropy difference would produce a substantial free energy change, accurately resolving such free energy changes into entropic and enthalpic components has been difficult. Here, we attempt to quantify the unfolding entropy change produced by P --> A substitution by amplifying the effect through multiple substitutions, and by decreasing the uncertainty in determining the unfolding entropy. Variants of a repeat protein, the Drosophila Notch ankyrin domain, were constructed with a varying number of P --> A substitutions at structurally conserved positions. Unfolding entropy values of the variants were determined from free energy measurements taken over a common temperature range using chemical denaturation. Our findings confirm the prediction that increasing the number of proline residues present in similar local environments increases the unfolding entropy. The average value of this increase in unfolding entropy is 7.7 +/- 4.2 cal mol(-1) K(-1), which is within error of the value estimated by Matthews et al. (1987).     

Filling small, empty protein cavities: structural and energetic consequences

Most proteins contain small cavities that can be filled by replacing cavity-lining residues by larger ones. Since shortening mutations in hydrophobic cores tend to destabilize proteins, it is expected that cavity-filling mutations may conversely increase protein stability. We have filled three small cavities in apoflavodoxin and determined by NMR and equilibrium unfolding analysis their impact in protein structure and stability. The smallest cavity (14 A3) has been filled, at two different positions, with a variety of residues and, in all cases, the mutant proteins are locally unfolded, their structure and energetics resembling those of an equilibrium intermediate of the thermal unfolding of the wild-type protein. In contrast, two slightly larger cavities of 20 A3 and 21 A3 have been filled with Val to Ile or Val to Leu mutations and the mutants preserve both the native fold and the equilibrium unfolding mechanism. From the known relationship, observed in shortening mutations, between stability changes and the differential hydrophobicity of the exchanged residues and the volume of the cavities, the filling of these apoflavodoxin cavities is expected to stabilize the protein by approximately 1.5 kcal mol(-1). However, both urea and thermal denaturation analysis reveal much more modest stabilizations, ranging from 0.0 kcal mol(-1) to 0.6 kcal mol(-1), which reflects that the accommodation of single extra methyl groups in small cavities requires some rearrangement, necessarily destabilizing, that lowers the expected theoretical stabilization. As the size of these cavities is representative of that of the typical small, empty cavities found in most proteins, it seems unlikely that filling this type of cavities will give rise to large stabilizations.     

Exploring the Denatured State Ensemble by Single-Molecule Chemo-Mechanical Unfolding: The Effect of Force,Temperature, and Urea

While it is widely appreciated that the denatured state of a protein is a heterogeneous conformational ensemble, there is still debate over how this ensemble changes with environmental conditions. Here, we use single-molecule chemo-mechanical unfolding, which combines force and urea using the optical tweezers, together with traditional protein unfolding studies to explore how perturbants commonly used to unfold proteins (urea, force, and temperature) affect the denatured-state ensemble. We compare the urea m-values, which report on the change in solvent accessible surface area for unfolding, to probe the denatured state as a function of force, temperature, and urea. We find that while the urea- and force-induced denatured states expose similar amounts of surface area, the denatured state at high temperature and low urea concentration is more compact. To disentangle these two effects, we use destabilizing mutations that shift the T_m and C_m. We find that the compaction of the denatured state is related to changing temperature as the different variants of acyl-coenzyme A binding protein have similar m-values when they are at the same temperature but different urea concentration. These results have important implications for protein folding and stability under different environmental conditions.     

Denaturation induced aggregation in α‐crystallin: differential action of chaotropes

α‐Crystallin is a member of small heat shock proteins and is believed to play an exceptional role in the stability of eye lens proteins. The disruption or denaturation of the protein arrangement or solubility of the crystallin proteins can lead to vision problems including cataract. In the present study, we have examined the effect of chemical denaturants urea and guanidine hydrochloride (GdnHCl) on α‐crystallin aggregation, with special emphasis on protein conformational changes, unfolding, and amyloid fibril formation. GdnHCl (4 M) induced a 16 nm red shift in the intrinsic fluorescence of α‐crystallin, compared with 4 nm shift by 8 M urea suggesting a major change in α‐crystallin structure. Circular dichroism analysis showed marked increase in the ellipticity of α‐crystallin at 216 nm, suggesting gain in β‐sheet structure in the presence of GdnHCl (0.5–1 M) followed by unfolding at higher concentration (2–6 M). However, only minor changes in the secondary structure of α‐crystallin were observed in the presence of urea. Moreover, 8‐anilinonaphthalene‐1‐sulfonic acid fluorescence measurement in the presence of GdnHCl and urea showed changes in the hydrophobicity of α‐crystallin. Amyloid studies using thioflavin T fluorescence and congo red absorbance showed that GdnHCl induced amyloid formation in α‐crystallin, whereas urea induced aggregation in this protein. Electron microscopy studies further confirmed amyloid formation of α‐crystallin in the presence of GdnHCl, whereas only aggregate‐like structures were observed in α‐crystallin treated with urea. Our results suggest that α‐crystallin is susceptible to unfolding in the presence of chaotropic agents like urea and GdnHCl. The destabilized protein has increased likelihood to fibrillate. Copyright © 2016 John Wiley & Sons, Ltd.     

Energetics of aliphatic deletions in protein cores

Although core residues can sometimes be replaced by shorter ones without introducing significant changes in protein structure, the energetic consequences are typically large and destabilizing. Many efforts have been devoted to understand and predict changes in stability from analysis of the environment of mutated residues, but the relationships proposed for individual proteins have often failed to describe additional data. We report here 17 apoflavodoxin large-to-small mutations that cause overall protein destabilizations of 0.6-3.9 kcal.mol(-1). By comparing two-state urea and three-state thermal unfolding data, the overall destabilizations observed are partitioned into effects on the N-to-I and on the I-to-U equilibria. In all cases, the equilibrium intermediate exerts a "buffering" effect that reduces the impact of the overall destabilization on the N-to-I equilibrium. The performance of several structure-energetics relationships, proposed to explain the energetics of hydrophobic shortening mutations, has been evaluated by using an apoflavodoxin data set consisting of 14 mutations involving branching-conservative aliphatic side-chain shortenings and a larger data set, including similar mutations implemented in seven model proteins. Our analysis shows that the stability changes observed for any of the different types of mutations (LA, IA, IV, and VA) in either data set are best explained by a combination of differential hydrophobicity and of the calculated volume of the modeled cavity (as previously observed for LA and IA mutations in lysozyme T4). In contrast, sequence conservation within the flavodoxin family, which is a good predictor for charge-reversal stabilizing mutations, does not perform so well for aliphatic shortening ones.     

The peculiar nature of unfolding of the human prion protein

Spontaneous conformational transition of the prion protein from an alpha-helical isoform to a beta-sheet-rich isoform underlies the pathogenesis of sporadic prion diseases. To study the rate-limiting steps of spontaneous conversion, the formation of amyloid fibrils by the recombinant human PrP C-terminal fragment spanning residues 90-231 (recPrP) was monitored in the presence of urea. The kinetics of spontaneous fibril formation displayed sigmoidal behavior involving a lag phase. The shortest lag phase was observed at partially denaturing conditions, close to the concentration of urea corresponding to the middle point of unfolding. This result indicates that unfolding intermediates may be important for the conversion. To test whether unfolding intermediates are formed, we employed size-exclusion chromatography and circular dichroism spectroscopy to monitor urea denaturation of recPrP. Both techniques showed a single sigmoidal transition with very similar thermodynamic parameters of denaturation and that the transition can be described by a simple equilibrium between folded and denatured states. Detailed analyses of data, however, revealed that the dimensions of both the native and denatured species gradually increases with urea. Expansion of the native species is also accompanied by an increase in efficiency of the energy transfer from a single Trp residue to 1-anilinonaphthalene-8-sulfonate dye as measured by fluorescence. These data illustrate that thermodynamic character of the native ensemble changes gradually with environmental conditions. Such behavior is consistent with the thermodynamically variable model, and may be linked to the ability of PrP to adopt distinct abnormal conformations under pathologic conditions.     

Heat capacity changes upon burial of polar and nonpolar groups in proteins

In this paper we address the question of whether the burial of polar and nonpolar groups in the protein locale is indeed accompanied by the heat capacity changes, DeltaC(p), that have an opposite sign, negative for nonpolar groups and positive for polar groups. To accomplish this, we introduced amino acid substitutions at four fully buried positions of the ubiquitin molecule (Val5, Val17, Leu67, and Gln41). We substituted Val at positions 5 and 17 and Leu at position 67 with a polar residue, Asn. As a control, Ala was introduced at the same three positions. We also replaced the buried polar Gln41 with Val and Leu, nonpolar residues that have similar size and shape as Gln. As a control, Asn was introduced at Gln41 as well. The effects of these amino acid substitutions on the stability, and in particular, on the heat capacity change upon unfolding were measured using differential scanning calorimetry. The effect of the amino acid substitutions on the structure was also evaluated by comparing the (1)H-(15)N HSQC spectra of the ubiquitin variants. It was found that the Ala substitutions did not have a considerable effect on the heat capacity change upon unfolding. However, the substitutions of aliphatic side chains (Val or Leu) with a polar residue (Asn) lead to a significant (> 30%) decrease in the heat capacity change upon unfolding. The decrease in heat capacity changes does not appear to be the result of significant structural perturbations as seen from the HSQC spectra of the variants. The substitution of a buried polar residue (Gln41) to a nonpolar residue (Leu or Val) leads to a significant (> 25%) increase in heat capacity change upon unfolding. These results indicate that indeed the heat capacity change of burial of polar and nonpolar groups has an opposite sign. However, the observed changes in DeltaC(p) are several times larger than those predicted, based on the changes in water accessible surface area upon substitution.     

Role of conserved residues in structure and stability: tryptophans of human serum retinol-binding protein, a model for the lipocalin superfamily

Serum retinol binding protein (RBP) is a member of the lipocalin family, proteins with up-and-down beta-barrel folds, low levels of sequence identity, and diverse functions. Although tryptophan 24 of RBP is highly conserved among lipocalins, it does not play a direct role in activity. To determine if Trp24 and other conserved residues have roles in stability and/or folding, we investigated the effects of conservative substitutions for the four tryptophans and some adjacent residues on the structure, stability, and spectroscopic properties of apo-RBP. Crystal structures of recombinant human apo-RBP and of a mutant with substitutions for tryptophans 67 and 91 at 1.7 A and 2.0 A resolution, respectively, as well as stability measurements, indicate that these relatively exposed tryptophans have little influence on structure or stability. Although Trp105 is largely buried in the wall of the beta-barrel, it can be replaced with minor effects on stability to thermal and chemical unfolding. In contrast, substitutions of three different amino acids for Trp24 or replacement of Arg139, a conserved residue that interacts with Trp24, lead to similar large losses in stability and lower yields of native protein generated by in vitro folding. The results and the coordinated nature of natural substitutions at these sites support the idea that conserved residues in functionally divergent homologs have roles in stabilizing the native relative to misfolded structures. They also establish conditions for studies of the kinetics of folding and unfolding by identifying spectroscopic signals for monitoring the formation of different substructures.     

Amino acid substitutions affecting protein dynamics in eglin C do not affect heat capacity change upon unfolding

The heat capacity change upon unfolding (deltaC(p)) is a thermodynamic parameter that defines the temperature dependence of the thermodynamic stability of proteins; however, physical basis of the heat capacity change is not completely understood. Although empirical surface area-based calculations can predict heat capacity changes reasonably well, accumulating evidence suggests that changes in hydration of those surfaces is not the only parameter contributing to the observed heat capacity changes upon unfolding. Because packing density in the protein interior is similar to that observed in organic crystals, we hypothesized that changes in protein dynamics resulting in increased rigidity of the protein structure might contribute to the observed heat capacity change upon unfolding. Using differential scanning calorimetry we characterized the thermodynamic behavior of a serine protease inhibitor eglin C and two eglin C variants with altered native state dynamics, as determined by NMR. We found no evidence of changes in deltaC(p) in either of the variants, suggesting that changes in rigidity do not contribute to the heat capacity change upon unfolding in this model system.     

Comparison of the folding mechanism of highly homologous proteins in the lipid‐binding protein family

The folding mechanism of two closely related proteins in the intracellular lipid‐binding protein family, human bile acid‐binding protein (hBABP), and rat bile acid‐binding protein (rBABP) were examined. These proteins are 77% identical (93% similar) in sequence. Both of these single domain proteins fit well to a two‐state model for unfolding by fluorescence and circular dichroism at equilibrium. Three phases were observed during the unfolding of rBABP by fluorescence but only one phase was observed during the unfolding of hBABP, suggesting that at least two kinetic intermediates accumulate during the unfolding of rBABP that are not observed during the unfolding of hBABP. Fluorine NMR was used to examine the equilibrium unfolding behavior of the W49 side chain in 6‐fluorotryptophan‐labeled rBABP and hBABP. The structure of rBABP appears to be more dynamic than that of hBABP in the vicinity of W49 in the absence of denaturant, and urea has a greater effect on this dynamic behavior for rBABP than for hBABP. As such, the folding behavior of highly sequence related proteins in this family can be quite different. These differences imply that moderately sized proteins with high sequence and structural similarity can still populate quite different structures during folding. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Effect of Substitutions in Surface Amino Acid on Energy Profile of Apomyoglobin

Studies on the process of spontaneous protein folding into a unique native state are an important issue of molecular biology. Apomyoglobin from the sperm whale is a convenient model for these studies in vitro. Here, we present the results of equilibrium and kinetic experiments carried out in a study on the folding and unfolding of eight mutant apomyoglobin forms of with hydrophobic amino acid substitutions on the protein surface. Calculated values of apparent constants of folding/unfolding rates, as well as the data on equilibrium conformational transitions in the urea concentration range of 0–6 M at 11°C are given. Based on the obtained information on the kinetic properties of the studied proteins, a Φ-value analysis of the transition state has been performed and values of urea concentrations corresponding to the midpoint of the transition from the native to intermediate state have been determined for the given forms of mutant apomyoglobin. It has been found that a significant increase in the stability of the native state can be achieved by a small number of amino acid substitutions on the protein surface. It has been shown that the substitution of only one amino acid residue exclusively affects the height of the energy barrier that separates different states of apomyoglobin.     

分子动力学模拟尿素对水溶液中蛋白质构象转变的影响

采用分子动力学方法和全原子模型研究尿素和水分子对模型蛋白S-肽链结构转化的影响。模拟结果显示S-肽链的变性速率常数k值随着尿素浓度的增加而先降低后升高,在尿素浓度为2.9 mol/L时达到最低值。模拟了不同尿素浓度下尿素-肽链、水-肽链以及肽链分子氢键的形成状况。结果表明:尿素浓度较低时,尿素分子与S-肽链的极性氨基酸侧链形成氢键,但不破坏其分子内的骨架氢键,尿素在S-肽链水化层外形成限制性空间,增强了S-肽链的稳定性。随着尿素的升高,尿素分子进入S-肽链内部并与其内部氨基酸残基形成氢键,导致S-肽链的骨架氢键丧失,S-肽链发生去折叠。上述模拟结果与文献报道的实验结果一致,从分子水平上揭示了尿素对蛋白质分子结构变化的影响机制,对于研究和发展蛋白质折叠及稳定化技术具有指导意义。