共查询到20条相似文献，搜索用时 15 毫秒
1.
2.
Stefan Knapp Pekka T. Mattson Petya Christova Kurt D. Berndt Andrej Karshikoff Mauno Vihinen C.I. Edvard Smith Rudolf Ladenstein 《Proteins》1998,31(3):309-319
The thermal unfolding of three SH3 domains of the Tec family of tyrosine kinases was studied by differential scanning calorimetry and CD spectroscopy. The unfolding transition of the three protein domains in the acidic pH region can be described as a reversible two-state process. For all three SH3 domains maximum stability was observed in the pH region 4.5 < pH < 7.0 where these domains unfold at temperatures of 353K (Btk), 342K (Itk), and 344K (Tec). At these temperatures an enthalpy change of 196 kJ/mol, 178 kJ/mol, and 169 kJ/mol was measured for Btk-, Itk-, and Tec-SH3 domains, respectively. The determined changes in heat capacity between the native and the denatured state are in an usual range expected for small proteins. Our analysis revealed that all SH3 domains studied are only weakly stabilized and have free energies of unfolding which do not exceed 12–16 kJ/mol but show quite high melting temperatures. Comparing unfolding free energies measured for eukaryotic SH3 domains with those of the topologically identical Sso7d protein from the hyperthermophile Sulfolobus solfataricus, the increased melting temperature of the thermostable protein is due to a broadening as well as a significant lifting of its stability curve. However, at their physiological temperatures, 310K for mesophilic SH3 domains and 350K for Sso7d, eukaryotic SH3 domains and Sso7d show very similar stabilities. Proteins 31:309–319, 1998. © 1998 Wiley-Liss, Inc. 相似文献
3.
Kalyanashis Jana Rukmankesh Mehra Budheswar Dehury Tom L. Blundell Kasper P. Kepp 《Proteins》2020,88(9):1233-1250
Protein thermostability is important to evolution, diseases, and industrial applications. Proteins use diverse molecular strategies to achieve stability at high temperature, yet reducing the entropy of unfolding seems required. We investigated five small α-proteins and five β-proteins with known, distinct structures and thermostability (Tm) using multi-seed molecular dynamics simulations at 300, 350, and 400 K. The proteins displayed diverse changes in hydrogen bonding, solvent exposure, and secondary structure with no simple relationship to Tm. Our dynamics were in good agreement with experimental B-factors at 300 K and insensitive to force-field choice. Despite the very distinct structures, the native-state (300 + 350 K) free-energy landscapes (FELs) were significantly broader for the two most thermostable proteins and smallest for the three least stable proteins in both the α- and β-group and with both force fields studied independently (tailed t-test, 95% confidence level). Our results suggest that entropic ensembles stabilize proteins at high temperature due to reduced entropy of unfolding, viz., ΔG = ΔH − TΔS. Supporting this mechanism, the most thermostable proteins were also the least kinetically stable, consistent with broader FELs, typified by villin headpiece and confirmed by specific comparison to a mesophilic ortholog of Thermus thermophilus apo-pyrophosphate phosphohydrolase. We propose that molecular strategies of protein thermostabilization, although diverse, tend to converge toward highest possible entropy in the native state consistent with the functional requirements. We speculate that this tendency may explain why many proteins are not optimally structured and why molten-globule states resemble native proteins so much. 相似文献
4.
L. Wang Y. Duan R. Shortle B. Imperiali P. A. Kollman 《Protein science : a publication of the Protein Society》1999,8(6):1292-1304
BBA1 is a designed protein that has only 23 residues. It is the smallest protein without disulfide bridges that has a well-defined tertiary structure in solution. We have performed unfolding molecular dynamics simulations on BBA1 and some of its mutants at 300, 330, 360, and 400 K to study their kinetic stability as well as the unfolding mechanism of BBA1. It was shown that the unfolding simulations can provide insights into the forces that stabilize the protein. Packing, hydrophobic interactions, and a salt bridge between Asp12 and Lys16 were found to be important to the protein's stability. The unfolding of BBA1 goes through two major steps: (1) disruption of the hydrophobic core and (2) unfolding of the helix. The beta-hairpin remains stable in the unfolding because of the high stability of the type II' turn connecting the two beta-strands. 相似文献
5.
Biswajit Gorai Arasu Prabhavadhni 《Journal of biomolecular structure & dynamics》2013,31(9):2037-2047
Unfolding stabilities of two homologous proteins, cardiotoxin III and short-neurotoxin (SNTX) belonging to three-finger toxin (TFT) superfamily, have been probed by means of molecular dynamics (MD) simulations. Combined analysis of data obtained from steered MD and all-atom MD simulations at various temperatures in near physiological conditions on the proteins suggested that overall structural stabilities of the two proteins were different from each other and the MD results are consistent with experimental data of the proteins reported in the literature. Rationalization for the differential structural stabilities of the structurally similar proteins has been chiefly attributed to the differences in the structural contacts between C- and N-termini regions in their three-dimensional structures, and the findings endorse the ‘CN network’ hypothesis proposed to qualitatively analyse the thermodynamic stabilities of proteins belonging to TFT superfamily of snake venoms. Moreover, the ‘CN network’ hypothesis has been revisited and the present study suggested that ‘CN network’ should be accounted in terms of ‘structural contacts’ and ‘structural strengths’ in order to precisely describe order of structural stabilities of TFTs. 相似文献
6.
7.
The beta-hairpin trpzip2 can be tuned continuously from a two-state folder to folding on a rough energy landscape without a dominant refolding barrier. At high denaturant concentration, this extremely stable peptide exhibits a single apparent "two-state" transition temperature when monitored by different spectroscopic techniques. However, under optimal folding conditions the hairpin undergoes an unusual folding process with three clusters of melting transitions ranging from 15 degrees C to 160 degrees C, as monitored by 12 different experimental and computational observables. We explain this behavior in terms of a rough free energy landscape of the unfolded peptide caused by multiple tryptophan interactions and alternative backbone conformations. The landscape is mapped out by potentials of mean force derived from replica-exchange molecular dynamics simulations. Implications for deducing cooperativity from denaturant titrations, for the origin of folding cooperativity, and for the folding of thermophilic proteins are pointed out. trpzip is an excellent small tunable model system for the glass-like folding transitions predicted by landscape theory. 相似文献
8.
In this work, the conformational behavior of ribonuclease Sso7d is studied as a function of chirality of its constituting amino acids. Both optimized structures (using molecular mechanics with the CHARMM force field) and dynamic behavior (obtained by molecular dynamic simulations) are compared. 相似文献
9.
Sso7d from the extreme thermophilic crenarchaeon Sulfolobus solfataricus is a multifunctional protein in in vitro assays, whose in vivo role is still puzzling. Crystals of Sso7d in complex with DNA elucidated the protein surface involved in the binding to the nucleic acid, whereas the locations of the Sso7d regions responsible for a chaperone activity in renaturing protein aggregates (i.e., the protein-binding surface and the site of ATPase activity) are still unknown. We identified the regions of Sso7d involved in protein-binding by limited proteolysis experiments associated to advanced mass spectrometric procedures performed on isolated Sso7d and Sso7d in complex with the peptide melittin. By affinity labeling of Sso7d with the ATP analogue 5'-p-fluorosulfonylbenzoyl adenosine and characterization of the labeled tryptic peptides by tandem mass spectrometry, we found that Y7 and K39 are residues involved in ATP binding/hydrolysis. Insights into the positions of the ligands melittin and ATP were achieved by a molecular modeling study; the models obtained were in agreement with most experimental data. A comparison among the complexes of Sso7d with DNA, with melittin, and with ATP showed that the DNA-binding surface and the protein-binding surface overlap, whereas the ATPase site is mostly independent of the binding sites for the nucleic acid and melittin. 相似文献
10.
B.L. de Groot D.M.F. van Aalten R.M. Scheek A. Amadei G. Vriend H.J.C. Berendsen 《Proteins》1997,29(2):240-251
A method is presented that generates random protein structures that fulfil a set of upper and lower interatomic distance limits. These limits depend on distances measured in experimental structures and the strength of the interatomic interaction. Structural differences between generated structures are similar to those obtained from experiment and from MD simulation. Although detailed aspects of dynamical mechanisms are not covered and the extent of variations are only estimated in a relative sense, applications to an IgG-binding domain, an SH3 binding domain, HPr, calmodulin, and lysozyme are presented which illustrate the use of the method as a fast and simple way to predict structural variability in proteins. The method may be used to support the design of mutants, when structural fluctuations for a large number of mutants are to be screened. The results suggest that motional freedom in proteins is ruled largely by a set of simple geometric constraints. Proteins 29:240–251, 1997. © 1997 Wiley-Liss, Inc. 相似文献
11.
Transient thermal unfolding of ubiquitin is investigated using nonlinear infrared spectroscopy after a nanosecond laser temperature jump (T-jump). The abrupt change in the unfolding free energy surface and the ns time resolution allow us to observe a fast response on ns to micros time-scales, which we attribute to downhill unfolding, before a cross-over to ms kinetics. The downhill unfolding by a sub-population of folded proteins is induced through a shift of the barrier toward the native state. By adjusting the T-jump width, the effect of the initial (T(i)) and final (T(f)) temperature on the unfolding dynamics can be separated. From the amplitude of the fast downhill unfolding, the fractional population prepared at the unfolding transition state is obtained. This population increases with both T(i) and with T(f). A two-state kinetic analysis of the ms refolding provides thermodynamic information about the barrier height. By a combination of the fast and slow unfolding and folding parameters, a quasi-two-state kinetic analysis is performed to calculate the time-dependent population changes of the folded state. This calculation coincides with the experimentally obtained population changes at low temperature but deviations are found in the T-jump from 67 to 78 degrees C. Using temperature-dependent barrier height changes, a temperature Phi value analysis is performed. The result shows a decreasing trend of Phi(T) with temperature, which indicates an increase of the heterogeneity of the transition state. We conclude that ubiquitin unfolds along a well-defined pathway at low temperature which expands with increasing temperature to include multiple routes. 相似文献
12.
Thermal stability of mutant proteins has been investigated using temperature dependent molecular dynamics (MD) simulations in vacuo. The numerical modeling was aimed at mimicking protein expansion upon heating. After the conditions for an expanding protein accessible surface area were established for T4 lysozyme and barnase wild-type proteins, MD simulations were carried out under the same conditions using the crystal structures of several mutant proteins. The computed thermal expansion of the accessible surface area of mutant proteins was found to be strongly correlated with their experimentally measured stabilities. A similar, albeit weaker, correlation was observed for model mutant proteins. This opens the possibility of obtaining stability information directly from protein structure. 相似文献
13.
The impacts of three charged‐residue‐involved mutations, E46A, R3E, and R3E/L66E, on the thermostability and folding behavior of the cold shock protein from the themophile Bacillus caldolyticus (Bc‐Csp) were investigated by using a modified Gō‐like model, in which the nonspecific electrostatic interactions of charged residues were taken into account. Our simulation results show that the wild‐type Bc‐Csp and its three mutants are all two‐sate folders, which is consistent with the experimental observations. It is found that these three mutations all lead to a decrease of protein thermodynamical stability, and the effect of R3E mutation is the strongest. The lower stability of these three mutants is due to the increase of the enthalpy of the folded state and the entropy of the unfolded state. Using this model, we also studied the folding kinetics and the folding/unfolding pathway of the wild‐type Bc‐Csp as well as its three mutants and then discussed the effects of electrostatic interactions on the folding kinetics. The results indicate that the substitutions at positions 3 and 46 largely decrease the folding kinetics, whereas the mutation of residue 66 only slightly decreases the folding rate. This result agrees well with the experimental observations. It is also found that these mutations have little effects on the folding transition state and the folding pathway, in which the N‐terminal β sheet folds earlier than the C‐terminal region. We also investigated the detailed unfolding pathway and found that it is really the reverse of the folding pathway, providing the validity of our simulation results. Proteins 2010. © 2010 Wiley‐Liss, Inc. 相似文献
14.
MacDonald JT Purkiss AG Smith MA Evans P Goodfellow JM Slingsby C 《Protein science : a publication of the Protein Society》2005,14(5):1282-1292
The thermodynamic and kinetic stabilities of the eye lens family of betagamma-crystallins are important factors in the etiology of senile cataract. They control the chance of proteins unfolding, which can lead to aggregation and loss of transparency. betaB2-Crystallin orthologs are of low stability and comprise two typical betagamma-crystallin domains, although, uniquely, the N-terminal domain has a cysteine in one of the conserved folded beta-hairpins. Using high-temperature (500 K) molecular dynamics simulations with explicit solvent on the N-terminal domain of rodent betaB2-crystallin, we have identified in silico local flexibility in this folded beta-hairpin. We have shown in vitro using two-domain human betaB2-crystallin that replacement of this cysteine with a more usual aromatic residue (phenylalanine) results in a gain in conformational stability and a reduction in the rate of unfolding. We have used principal components analysis to visualize and cluster the coordinates from eight separate simulated unfolding trajectories of both the wild-type and the C50F mutant N-terminal domains. These data, representing fluctuations around the native well, show that although the mutant and wild-type appear to behave similarly over the early time period, the wild type appears to explore a different region of conformational space. It is proposed that the advantage of having this low-stability cysteine may be correlated with a subunit-exchange mechanism that allows betaB2-crystallin to interact with a range of other beta-crystallin subunits. 相似文献
15.
The molecular basis of conformational diseases frequently resides in mutant proteins constituting a subset of the vast mutational space. While the subtleties of protein structure point to molecular dynamics (MD) techniques as promising tools for an efficient exploration of such a space, the average size of proteins and the time scale of unfolding events make this goal difficult with present computational capabilities. We show here, nevertheless, that an efficient approach is already feasible for modular proteins. Familial hypercholesterolemia (FH) is a conformational disease linked to mutations in the gene encoding the low density lipoprotein receptor. A high percentage of these mutations has been found in the seven small modular binding repeats of the receptor. Taking advantage of its small size, we have performed an in depth MD study of the fifth binding repeat. Fast unfolding dynamics have been observed in the absence of a structural bound calcium ion, which agrees with its reported essential role in the stability of the module. In addition, several mutations detected in FH patients have been analyzed, starting from the native conformation. Our results indicate that in contrast with the wild type protein and an innocuous control mutant, disease-related mutants experience, in short simulation times (2-8 ns), gross departures from the native state that lead to unfolded conformations and, in some cases, to binding site desorganization deriving in calcium release. Computational diagnosis of mutations leading to conformational diseases seems thus feasible, at least for small or modular pathogenic proteins. 相似文献
16.
The protein moiety of the Braun's E. coli outer membrane lipoprotein (Lpp-56) is an attractive object of biophysical investigation in several aspects. It is a homotrimeric, parallel coiled coil, a class of coiled coils whose stability and folding have been studied only occasionally. Lpp-56 possesses unique structural properties and exhibits extremely low rates of folding and unfolding. It is natural to ask how the specificity of the structure determines the extraordinary physical chemical properties of this protein. Recently, a seemingly controversial data on the stability and unfolding rate of Lpp-56 have been published (Dragan et al., Biochemistry 2004;43: 14891-14900; Bjelic et al., Biochemistry 2006;45:8931-8939). The unfolding rate constant measured using GdmCl as the denaturing agent, though extremely low, was substantially higher than that obtained on the basis of thermal unfolding. If this large difference arises from the effect of screening of electrostatic interactions induced by GdmCl, electrostatic interactions would appear to be an important factor determining the unusual properties of Lpp-56. We present here a computational analysis of the electrostatic properties of Lpp-56 combining molecular dynamics simulations and continuum pK calculations. The pH-dependence of the unfolding free energy is predicted in good agreement with the experimental data: the change in DeltaG between pH 3 and pH 7 is approximately 60 kJ mol(-1). The results suggest that the difference in the stability of the protein observed using different experimental methods is mainly because of the effect of the reduction of electrostatic interactions when the salt (GdmCl) concentration increases. We also find that the occupancy of the interhelical salt bridges is unusually high. We hypothesize that electrostatic interactions, and the interhelical salt bridges in particular, are an important factor determining the low unfolding rate of Lpp-56. 相似文献
17.
To provide evidence for the interpretation of temperature‐dependent unfolding kinetics and the downhill unfolding scenario presented in the accompanying experimental article (Part I), the free energy surface of ubiquitin unfolding is calculated using statistical mechanical models of the Muñoz‐Eaton (ME) form. The models allow only two states for each amino acid residue, folded or unfolded, and permutations of these states generate an ensemble of microstates. One‐dimensional free energy curves are calculated using the number of folded residues as a reaction coordinate. The proposed sequential unfolding of ubiquitin's β‐sheet is tested by mapping the free energy onto two reaction coordinates inspired by the experiment as follows: the number of folded residues in ubiquitin's stable β‐strands I and II and those of the less stable strands III–V. Although the original ME model successfully captures folding features of zipper‐like one‐dimensional folders, it misses important tertiary interactions between residues that are far from each other in primary sequence. To take tertiary contacts into account, partially folded microstates based on a spherical growth model are included in the calculation and compared with the original model. By calculating the folding probability of each residue for a given point on the free energy surface, the unfolding pathway of ubiquitin is visualized. At low temperature, thermal unfolding occurs along a sequential unfolding pathway as follows: disruption of the β‐strands III–V followed by unfolding of the strands I and II. At high temperature, multiple unfolding routes are formed. The heterogeneity of the transition state explains the global nonexponential unfolding observed in the T‐jump experiment at high temperature. The calculation also reports a high stability for the α‐helix of ubiquitin. Proteins 2008. © 2008 Wiley‐Liss, Inc. 相似文献
18.
We have traditionally relied on extremely elevated temperatures (498K, 225 degrees C) to investigate the unfolding process of proteins within the timescale available to molecular dynamics simulations with explicit solvent. However, recent advances in computer hardware have allowed us to extend our thermal denaturation studies to much lower temperatures. Here we describe the results of simulations of chymotrypsin inhibitor 2 at seven temperatures, ranging from 298K to 498K. The simulation lengths vary from 94ns to 20ns, for a total simulation time of 344ns, or 0.34 micros. At 298K, the protein is very stable over the full 50ns simulation. At 348K, corresponding to the experimentally observed melting temperature of CI2, the protein unfolds over the first 25ns, explores partially unfolded conformations for 20ns, and then refolds over the last 35ns. Above its melting temperature, complete thermal denaturation occurs in an activated process. Early unfolding is characterized by sliding or breathing motions in the protein core, leading to an unfolding transition state with a weakened core and some loss of secondary structure. After the unfolding transition, the core contacts are rapidly lost as the protein passes on to the fully denatured ensemble. While the overall character and order of events in the unfolding process are well conserved across temperatures, there are substantial differences in the timescales over which these events take place. We conclude that 498K simulations are suitable for elucidating the details of protein unfolding at a minimum of computational expense. 相似文献
19.
玉米赤霉烯酮是世界上污染最为广泛的一种镰刀菌(Fusarium)毒素,严重危害牲畜以及人类的健康。来源于粉红螺旋聚孢霉(Clonostachys rosea)的玉米赤霉烯酮水解酶(zearalenone hydrolase,ZHD)能有效降解饲料中的玉米赤霉烯酮,然而饲料加工中的高温环境限制了该酶的应用。基于结构特征的理性设计可为酶的热稳定性改造提供指导。本研究首先基于蛋白质结构比对(multiple structure alignment, MSTA)筛选ZHD的结构灵活区,随后基于序列保守性打分以及构象自由能计算设计突变文库,得到基于136号和220号残基的9个单点突变设计。结果表明,9个突变体的热熔融温度(Tm)提高了0.4–5.6℃,其中S220R和S220W热稳定性表现最好,Tm分别提高了5.6℃和4.0℃,45℃下的热半失活时间分别延长了15.4倍和3.1倍,相对酶活分别为野生型的70.6%和57.3%。分子动力学模拟分析表明突变位点及附近区域的作用力得到了增强,突变体S220R和S220W的220-K130氢键成键概率分别增加了37.1%和19.3%、K130-D223盐桥成键概率分别增加了30.1%和12.5%,为ZHD热稳定性的提高作出了贡献。这项工作表明结合天然酶的结构比对、序列分析及自由能计算的热稳定性改造策略的可行性,并获得了热稳定性增强的ZHD变体,为ZHD在工业上的应用打下基础。 相似文献
20.
The analysis of the dynamic behavior of enzymes is fundamental to structural biology. A direct relationship between protein flexibility and biological function has been shown for bovine pancreatic ribonuclease (RNase A) (Rasmussen et al., Nature 1992;357:423-424). More recently, crystallographic studies have shown that functional motions in RNase A involve the enzyme beta-sheet regions that move concertedly on substrate binding and release (Vitagliano et al., Proteins 2002;46:97-104). These motions have been shown to correspond to intrinsic dynamic properties of the native enzyme by molecular dynamics (MD) simulations. To unveil the occurrence of these collective motions in other members of pancreatic-like superfamily, we carried out MD simulations on human angiogenin (Ang). Essential dynamics (ED) analyses performed on the trajectories reveal that Ang exhibits collective motions similar to RNase A, despite the limited sequence identity (33%) of the two proteins. Furthermore, we show that these collective motions are also present in ensembles of experimentally determined structures of both Ang and RNase A. Finally, these subtle concerted beta-sheet motions were also observed for other two members of the pancreatic-like superfamily by comparing the ligand-bound and ligand-free structures of these enzymes. Taken together, these findings suggest that pancreatic-like ribonucleases share an evolutionary conserved dynamic behavior consisting of subtle beta-sheet motions, which are essential for substrate binding and release. 相似文献