Temperature range of thermodynamic stability for the native state of reversible two-state proteins

Loops around the Heme Pocket Have a Critical Role in the Function and Stability of BsDyP from Bacillus subtilis

Bacillus subtilis BsDyP belongs to class I of the dye-decolorizing peroxidase (DyP) family of enzymes and is an interesting biocatalyst due to its high redox potential, broad substrate spectrum and thermostability. This work reports the optimization of BsDyP using directed evolution for improved oxidation of 2,6-dimethoxyphenol, a model lignin-derived phenolic. After three rounds of evolution, one variant was identified displaying 7-fold higher catalytic rates and higher production yields as compared to the wild-type enzyme. The analysis of X-ray structures of the wild type and the evolved variant showed that the heme pocket is delimited by three long conserved loop regions and a small α helix where, incidentally, the mutations were inserted in the course of evolution. One loop in the proximal side of the heme pocket becomes more flexible in the evolved variant and the size of the active site cavity is increased, as well as the width of its mouth, resulting in an enhanced exposure of the heme to solvent. These conformational changes have a positive functional role in facilitating electron transfer from the substrate to the enzyme. However, they concomitantly resulted in decreasing the enzyme’s overall stability by 2 kcal mol⁻¹, indicating a trade-off between functionality and stability. Furthermore, the evolved variant exhibited slightly reduced thermal stability compared to the wild type. The obtained data indicate that understanding the role of loops close to the heme pocket in the catalysis and stability of DyPs is critical for the development of new and more powerful biocatalysts: loops can be modulated for tuning important DyP properties such as activity, specificity and stability.

Consensus protein engineering on the thermostable histone-like bacterial protein HUs significantly improves stability and DNA binding affinity

Consensus-based protein engineering strategy has been applied to various proteins and it can lead to the design of proteins with enhanced biological performance. Histone-like HUs comprise a protein family with sequence variety within a highly conserved 3D-fold. HU function includes compacting and regulating bacterial DNA in a wide range of biological conditions in bacteria. To explore the possible impact of consensus-based design in the thermodynamic stability of HU proteins, the approach was applied using a dataset of sequences derived from a group of 40 mesostable, thermostable, and hyperthermostable HUs. The consensus-derived HU protein was named HUBest, since it is expected to perform best. The synthetic HU gene was overexpressed in E. coli and the recombinant protein was purified. Subsequently, HUBest was characterized concerning its correct folding and thermodynamic stability, as well as its ability to interact with plasmid DNA. A substantial increase in HUBest stability at high temperatures is observed. HUBest has significantly improved biological performance at ambience temperature, presenting very low K_d values for binding plasmid DNA as indicated from the Gibbs energy profile of HUBest. This K_d may be associated to conformational changes leading to decreased thermodynamic stability and, therefore, higher flexibility at ambient temperature.

Advancing Therapeutic Protein Discovery and Development through Comprehensive Computational and Biophysical Characterization

Therapeutic protein candidates should exhibit favorable properties that render them suitable to become drugs. Nevertheless, there are no well-established guidelines for the efficient selection of proteinaceous molecules with desired features during early stage development. Such guidelines can emerge only from a large body of published research that employs orthogonal techniques to characterize therapeutic proteins in different formulations. In this work, we share a study on a diverse group of proteins, including their primary sequences, purity data, and computational and biophysical characterization at different pH and ionic strength. We report weak linear correlations between many of the biophysical parameters. We suggest that a stability comparison of diverse therapeutic protein candidates should be based on a computational and biophysical characterization in multiple formulation conditions, as the latter can largely determine whether a protein is above or below a certain stability threshold. We use the presented data set to calculate several stability risk scores obtained with an increasing level of analytical effort and show how they correlate with protein aggregation during storage. Our work highlights the importance of developing combined risk scores that can be used for early stage developability assessment. We suggest that such scores can have high prediction accuracy only when they are based on protein stability characterization in different solution conditions.

Exciton transport in the PE545 complex: insight from atomistic QM/MM-based quantum master equations and elastic network models

In this paper, we work out a parameterization of environmental noise within the Haken-Strobl-Reinenker (HSR) model for the PE545 light-harvesting complex, based on atomic-level quantum mechanics/molecular mechanics (QM/MM) simulations. We use this approach to investigate the role of various auto- and cross-correlations in the HSR noise tensor, confirming that site-energy autocorrelations (pure dephasing) terms dominate the noise-induced exciton mobility enhancement, followed by site energy-coupling cross-correlations for specific triplets of pigments. Interestingly, several cross-correlations of the latter kind, together with coupling-coupling cross-correlations, display clear low-frequency signatures in their spectral densities in the 30-70 [Formula: see text] region. These slow components lie at the limits of validity of the HSR approach, which requires that environmental fluctuations be faster than typical exciton transfer time scales. We show that a simple coarse-grained elastic-network-model (ENM) analysis of the PE545 protein naturally spotlights collective normal modes in this frequency range that represent specific concerted motions of the subnetwork of cysteines covalenty linked to the pigments. This analysis strongly suggests that protein scaffolds in light-harvesting complexes are able to express specific collective, low-frequency normal modes providing a fold-rooted blueprint of exciton transport pathways. We speculate that ENM-based mixed quantum classical methods, such as Ehrenfest dynamics, might be promising tools to disentangle the fundamental designing principles of these dynamical processes in natural and artificial light-harvesting structures.

In-silico prediction of concentration-dependent viscosity curves for monoclonal antibody solutions

Early stage developability assessments of monoclonal antibody (mAb) candidates can help reduce risks and costs associated with their product development. Forecasting viscosity of highly concentrated mAb solutions is an important aspect of such developability assessments. Reliable predictions of concentration-dependent viscosity behaviors for mAb solutions in platform formulations can help screen or optimize drug candidates for flexible manufacturing and drug delivery options. Here, we present a computational method to predict concentration-dependent viscosity curves for mAbs solely from their sequence-structural attributes. This method was developed using experimental data on 16 different mAbs whose concentration-dependent viscosity curves were experimentally obtained under standardized conditions. Each concentration-dependent viscosity curve was fitted with a straight line, via logarithmic manipulations, and the values for intercept and slope were obtained. Intercept, which relates to antibody diffusivity, was found to be nearly constant. In contrast, slope, the rate of increase in solution viscosity with solute concentration, varied significantly across different mAbs, demonstrating the importance of intermolecular interactions toward viscosity. Next, several molecular descriptors for electrostatic and hydrophobic properties of the 16 mAbs derived using their full-length homology models were examined for potential correlations with the slope. An equation consisting of hydrophobic surface area of full-length antibody and charges on V_H, V_L, and hinge regions was found to be capable of predicting the concentration-dependent viscosity curves of the antibody solutions. Availability of this computational tool may facilitate material-free high-throughput screening of antibody candidates during early stages of drug discovery and development.

On the role of flexibility in protein-ligand interactions: the example of p53 tetramerization domain

The recognition of protein surfaces by designed ligands has become an attractive approach in drug discovery. However, the variable nature and irregular behavior of protein surfaces defy this new area of research. The easy to understand "lock-and-key" model is far from being the ideal paradigm in biomolecular interactions and, hence, any new finding on how proteins and ligands behave in recognition events paves a step of the way. Herein, we illustrate a clear example on how an increase in flexibility of both protein and ligand can result in an increase in the stability of the macromolecular complex. The biophysical study of the interaction between a designed flexible tetraguanidinium-calix[4]arene and the tetramerization domain of protein p53 (p53TD) and its natural mutant p53TD-R337H shows how the floppy mutant domain interacts more tightly with the ligand than the well-packed wild-type protein. Moreover, the flexible calixarene ligand interacts with higher affinity to both wild-type and mutated protein domains than a conformationally rigid calixarene analog previously reported. These findings underscore the crucial role of flexibility in molecular recognition processes, for both small ligands and large biomolecular surfaces.

How do thermophilic proteins and proteomes withstand high temperature?

We attempt to understand the origin of enhanced stability in thermophilic proteins by analyzing thermodynamic data for 116 proteins, the largest data set achieved to date. We compute changes in entropy and enthalpy at the convergence temperature where different driving forces are maximally decoupled, in contrast to the majority of previous studies that were performed at the melting temperature. We find, on average, that the gain in enthalpy upon folding is lower in thermophiles than in mesophiles, whereas the loss in entropy upon folding is higher in mesophiles than in thermophiles. This implies that entropic stabilization may be responsible for the high melting temperature, and hints at residual structure or compactness of the denatured state in thermophiles. We find a similar trend by analyzing a homologous set of proteins classified based only on the optimum growth temperature of the organisms from which they were extracted. We find that the folding free energy at the temperature of maximal stability is significantly more favorable in thermophiles than in mesophiles, whereas the maximal stability temperature itself is similar between these two classes. Furthermore, we extend the thermodynamic analysis to model the entire proteome. The results explain the high optimal growth temperature in thermophilic organisms and are in excellent quantitative agreement with full thermal growth rate data obtained in a dozen thermophilic and mesophilic organisms.

Unfolding analysis of the mature and unprocessed forms of Bacillus licheniformis γ-glutamyltranspeptidase

Bacillus licheniformis γ-glutamyltranspeptidase (BlGGT) undergoes an autocatalytic process to generate 44.9 and 21.7 kDa subunits; however, a mutant protein (T399A) loses completely the processing ability and mainly exists as a precursor. For a comprehensive understanding of their structural features, the biophysical properties of these two proteins were investigated by circular dichroism and fluorescence spectroscopy. Tryptophan fluorescence and circular dichroism spectra were nearly identical for BlGGT and T399A, but unfolding analyses revealed that these two proteins had a different sensitivity towards temperature- and guanidine hydrochloride (GdnHCl)-induced denaturation. BlGGT and the unprocessed T399A displayed T(m) values of 61.4°C and 68.1°C, respectively, and thermal unfolding of both proteins was found to be highly irreversible. Fluorescence quenching analysis showed that T399A had a dynamic quenching constant similar to that of the wild-type enzyme. BlGGT started to unfold beyond ∼2.14 M GdnHCl and reached an unfolded intermediate, [GdnHCl](0.5, N - U), at 2.85 M, corresponding to free energy change [Formula: see text] of 12.34 kcal mol( - 1), whereas the midpoint of the denaturation curve for T399A was approximately 3.94 M, corresponding to a [Formula: see text] of 4.45 kcal mol( - 1). Taken together, it can be concluded that the structural stability of BlGGT is superior to that of T399A.

Analysis of the urinary peptidome associated with Helicobacter pylori infection

AIM: To investigate the relationship between urinary peptide changes and Helicobacter pylori (H. pylori) infection using urinary peptidome profiling.

METHODS: The study was performed in volunteers (n = 137) who gave informed consent. Urinary peptides were enriched by magnetic beads based weak cation exchange chromatography and spectrums acquired by matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) mass spectrometry (MS). ClinProTools bioinformatics software was used for statistical analysis and the recognition of peptide patterns. The marker peptides were identified by LTQ Obitrap XL tandem MS.

RESULTS: Approximately 50 proteins or peptides which loaded onto the magnetic beads were detected by MALDI-TOF MS. By optimizing the parameters of the model, the Genetic Algorithm model had good recognition capability (97%) and positive predictive value (94%). Based on the model, 2 markers with molecular masses of 6788 and 1912 Da were found that differentiated between H. pylori positive and negative volunteers. The m/z 1912 sequence was parsed as SKQFTSSTSYNRGDSTF. The peptide was identified as isoform 1 of the fibrinogen α chain precursor, whose concentration in urine was markedly higher in H. pylori infected volunteers than in H. pylori non-infected ones.

CONCLUSION: The appearance of urinary fibrinogen degradation products is caused by an active H. pylori-induced process.

Protein stability, flexibility and function

Proteins rely on flexibility to respond to environmental changes, ligand binding and chemical modifications. Potentially, a perturbation that changes the flexibility of a protein may interfere with its function. Millions of mutations have been performed on thousands of proteins in quests for a delineation of the molecular details of their function. Several of these mutations interfered with the binding of a specific ligand with a concomitant effect on the stability of the protein scaffold. It has been ambiguous and not straightforward to recognize if any relationships exist between the stability of a protein and the affinity for its ligand. In this review, we present examples of proteins where changes in stability results in changes in affinity and of proteins where stability and affinity are uncorrelated. We discuss the possibility for a relationship between stability and binding. From the data presented is it clear that there are specific sites (flexibility hotspots) in proteins that are important for both binding and stability. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

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.

Urinary proteomics--a tool for biomarker discovery

The strong need for the discovery of novel disease markers together with the development of high-throughput techniques that provide highly sensitive analysis of protein content in tissues and bodily fluids, using proteomics, has opened the completely new chapter in biomarker discovery. The detection of biomarkers based on urinary proteome analysis is rapidly advancing and may provide new tools to improve non-invasive diagnostics, prognostics, and therapy enhancement. As a tool for biomarker discovery, urinary proteomics is especially fruitful in the area of early diagnostics and differentiation of renal damage, and it possesses enormous potential for improving and expanding non-invasive cancer diagnostics. An abundance of urinary proteins could provide a wide variety of biomarkers for the diagnosis and follow-up of many systemic diseases as well. This article reviews the utility of urinary proteomics for biomarker discovery from the perspective of clinical application. Despite huge potential and prompt development of urinary proteomics, many challenges are still in front of us. Research effort and financial investment have to be oriented on providing strategies for exceeding current methodological and technical obstacles in a way to ensure the successful validation and implementation of newly discovered urinary biomarkers. The result is expected to be the development of new non-invasive tests and procedures able to guarantee higher efficiency of patient care and provide needed personalized medical approach.

Mass spectrometry based proteomics profiling as diagnostic tool in oncology: current status and future perspective

Proteomics analysis has been heralded as a novel tool for identifying new and specific biomarkers that may improve diagnosis and monitoring of various disease states. Recent years have brought a number of proteomics profiling technologies. Although proteomics profiling has resulted in the detection of disease-associated differences and modification of proteins, current proteomics technologies display certain limitations that are hampering the introduction of these new technologies into clinical laboratory diagnostics and routine applications. In this review, we summarize current advances in mass spectrometry based biomarker discovery. The promises and challenges of this new technology are discussed with particular emphasis on diagnostic perspectives of mass-spectrometry based proteomics profiling for malignant diseases.

A water-explicit lattice model of heat-, cold-, and pressure-induced protein unfolding

We investigate the effect of temperature and pressure on polypeptide conformational stability using a two-dimensional square lattice model in which water is represented explicitly. The model captures many aspects of water thermodynamics, including the existence of density anomalies, and we consider here the simplest representation of a protein: a hydrophobic homopolymer. We show that an explicit treatment of hydrophobic hydration is sufficient to produce cold, pressure, and thermal denaturation. We investigate the effects of the enthalpic and entropic components of the water-protein interactions on the overall folding phase diagram, and show that even a schematic model such as the one we consider yields reasonable values for the temperature and pressure ranges within which highly compact homopolymer configurations are thermodynamically stable.

Urinary mesothelin provides greater sensitivity for early stage ovarian cancer than serum mesothelin, urinary hCG free beta subunit and urinary hCG beta core fragment

OBJECTIVES

Early detection of ovarian cancer should improve overall survival. Multiple serum markers have been evaluated as possible tests to detect early stage disease, but few urine markers have been studied. Mesothelin has been detected in serum from patients with ovarian cancer, but has not been previously reported in urine.

METHODS

Mesothelin was assayed in the serum and in the urine from 28 patients with early stage (I/II) invasive epithelial ovarian cancers, 111 with advanced stage (III/IV) invasive disease and 19 with tumors of low malignant potential. Marker values have been compared to those in healthy controls and 115 patients with benign pelvic masses. Thresholds were set to include 95% of mesothelin values for 127 sera and 89 urines from healthy women. Urine values were considered: (1) as assayed; (2) normalized using the ratio of serum to urine creatinine; and (3) normalized using the Cockroft-Gault formula for glomerular filtration rate (GFR). Urines were also assayed for human chorionic gonadotropin (hCG) free beta subunit and beta subunit core fragment and similarly normalized.

RESULTS

Optimal sensitivity for early stage disease was obtained when urine mesothelin was normalized using GFR. A greater fraction of patients with early stage disease was detected with the mesothelin urine assay (42%) than with the serum assay (12%). Similarly, 75% of patients with advanced ovarian cancer had elevated mesothelin in urine compared to 48% in serum. Serum and urine levels of mesothelin correlated for early (p=0.02) and late (p<0.001) disease. Urine mesothelin exhibited greater sensitivity for early stage ovarian cancer than did hCG free beta subunit or beta subunit core fragment and complementarity was not observed.

CONCLUSION

Urine mesothelin deserves further evaluation as a biomarker for detection of early stage ovarian cancer in combination with other urinary markers.

Life at low temperatures: is disorder the driving force?

The thermodynamic characterization of various biological systems from psychrophiles points to a larger entropic contribution when compared to the corresponding mesophilic or (hyper) thermophilic counterparts, either at the level of the macromolecules (thermodynamic and kinetic stabilities) or of their function (ligand binding, catalytic activity). It is suggested here that in an environment characterized by a low heat content (enthalpy) and at temperatures that strongly slowdown molecular motions, the cold-adapted biological systems rely on a larger disorder to maintain macromolecular dynamics and function. Such pre-eminent involvement of entropy is observed in the experimental results and, from a macroscopic point of view, is also reflected for instance by the steric hindrances introduced by cis-unsaturated and branched lipids to maintain membrane fluidity, by the loose conformation of psychrophilic proteins or by the local destabilization of tRNA by dihydrouridine in psychrophilic bacteria.

Purification and characterization of a novel pectinase from Acrophialophora nainiana with emphasis on its physicochemical properties

An extracellular pectinase (PECI) was purified to apparent homogeneity from liquid state cultures of the thermophilic fungus Acrophialophora nainiana by ultrafiltration and a combination of gel filtration and ion-exchange chromatographic procedures. The molecular masses of PECI were 35,500 and 30,749 Da, as determined by SDS-PAGE and mass spectrometry, respectively. It was more active at 60 degrees C and pH 8.0 and showed high stability at 50 degrees C with half-life of 7 days. However at 60 and 70 degrees C, PECI was much less stable with half lives of approximately 20 and 3 min, respectively. The thermostability of purified PECI was also investigated by fluorescence and circular dichroism spectroscopy. Fluorescence revealed that the unfolding transition region was observed between 45 and 70 degrees C. A major decrease in the stability was found at 70 degrees C. Circular dichroism measurements at pH between 5.0 and 9.0 showed a transition temperature (T(m)) range of 50-55 degrees . The thermodynamic analysis of these results showed that EPGI is thermal stable protein exhibiting maximum stability (DeltaG(25)) of 22.65 and 19.19 kcal/mol at pH 8.0 and 9.0, respectively. The apparent K(m) value on pectin from citrus fruits was 4.22 mgml(-1). PECI exhibited no detectable activity of pectin methylesterase, endo-polygalacturonase, mannanase, xylanase and cellulase. However, it showed exo-polygalacturonase and pectin lyase activities. The presence of carbohydrate was detected in the pure PECI. It was activated by l-tryptophan, DEPC, DTT, DTNB, DTP, l-cystein and beta-mercaptoethanol and inhibited by NBS, Fe(2+), Cu(2+), Zn(2+), Mn(2+), Al(3+) and Ca(2+). The enzyme showed homology with a pectin lyases from Xanthomonas campestris and Bacillus licheniformis.

Proteomic-based discovery and characterization of glycosylated eosinophil-derived neurotoxin and COOH-terminal osteopontin fragments for ovarian cancer in urine

PURPOSE

The objective was to identify and characterize low molecular weight proteins/peptides in urine and their posttranslational modifications that might be used as a screening tool for ovarian cancer.

EXPERIMENTAL DESIGN

Urine samples collected preoperatively from postmenopausal women with ovarian cancer and benign conditions and from nonsurgical controls were analyzed by surface-enhanced laser desorption/ionization mass spectrometry and two-dimensional gel electrophoresis. Selected proteins from mass profiles were purified by chromatography and followed by liquid chromatography-tandem mass spectrometry sequence analysis. Specific antibodies were generated for further characterization, including immunoprecipitation and glycosylation. Quantitative and semiquantitative ELISAs were developed for preliminary validation in patients of 128 ovarian cancer, 52 benign conditions, 44 other cancers, and 188 healthy controls.

RESULTS

A protein (m/z approximately 17,400) with higher peak intensities in cancer patients than in benign conditions and controls was identified and subsequently defined as eosinophil-derived neurotoxin (EDN). A glycosylated form of EDN was specifically elevated in ovarian cancer patients. A cluster of COOH-terminal osteopontin was identified from two-dimensional gels of urine from cancer patients. Modified forms EDN and osteopontin fragments were elevated in early-stage ovarian cancers and a combination of both resulted to 93% specificity and 72% sensitivity.

CONCLUSIONS

Specific elevated posttranslationally modified urinary EDN and osteopontin COOH-terminal fragments in ovarian cancer might lead to potential noninvasive screening tests for early diagnosis. Urine with less complexity than serum and relatively high thermodynamic stability of peptides or metabolites is a promising study medium for discovery of the novel biomarkers which may present in many non-urinary tract neoplastic diseases.

Is it always possible to distinguish two- and three-state systems by evaluating the van't Hoff enthalpy?

Many small globular proteins are traditionally classified as thermodynamical two-state systems, i.e., the protein is either in the native, active state (folded) or in the denatured state (unfolded). We challenge this view and show that there may exist (protein) systems for which a van't Hoff analysis of experimental data cannot determine whether the system corresponds to two or three thermodynamical states when only temperatures in a narrow temperature region around the transition are considered. We generalize a widely employed two-state protein folding model to include a third, transition state. For this three-state system we systematically study the deviation of the calorimetric enthalpy (heat of transition) from the van't Hoff enthalpy, a measure of the two-stateness of a transition. We show that under certain conditions the heat capacity of the three-state system can be almost indistinguishable from the heat capacity for the two-state system over a broad temperature interval. The consequence may be that some three-state (or even more than three-states) systems have been misinterpreted as two-state systems when the conclusion is drawn solely upon the van't Hoff enthalpy. These findings are important not only for proteins, but also for the interpretation of thermodynamical systems in general.

Thermodynamic and dynamic factors involved in the stability of native protein structure in amorphous solids in relation to levels of hydration

The internal, dynamical fluctuations of protein molecules exhibit many of the features typical of polymeric and bulk small molecule glass forming systems. The response of a protein's internal molecular mobility to temperature changes is similar to that of other amorphous systems, in that different types of motions freeze out at different temperatures, suggesting they exhibit the alpha-beta-modes of motion typical of polymeric glass formers. These modes of motion are attributed to the dynamic regimes that afford proteins the flexibility for function but that also develop into the large-scale collective motions that lead to unfolding. The protein dynamical transition, T(d), which has the same meaning as the T(g) value of other amorphous systems, is attributed to the temperature where protein activity is lost and the unfolding process is inhibited. This review describes how modulation of T(d) by hydration and lyoprotectants can determine the stability of protein molecules that have been processed as bulk, amorphous materials. It also examines the thermodynamic, dynamic, and molecular factors involved in stabilizing folded proteins, and the effects typical pharmaceutical processes can have on native protein structure in going from the solution state to the solid state.

pH dependence thermal stability of a chymotrypsin inhibitor from Schizolobium parahyba seeds

The thermal stability of a Schizolobium parahyba chymotrypsin inhibitor (SPCI) as a function of pH has been investigated using fluorescence, circular dichroism, and differential scanning calorimetry (DSC). The thermodynamic parameters derived from all methods are remarkably similar and strongly suggest the high stability of SPCI under a wide range of pH. The transition temperature (T(m)) values ranging from 57 to 85.3 degrees C at acidic, neutral, and alkaline pH are in good agreement with proteins from mesophilic and thermophilic organisms and corroborate previous data regarding the thermal stability of SPCI. All methods gave transitions curves adequately fitted to a two-state model of the unfolding process as judged by the cooperative ratio between the van't Hoff and the calorimetric enthalpy energies close to unity in all of the pH conditions analyzed, except at pH 3.0. Thermodynamic analysis using all these methods reveals that SPCI is thermally a highly stable protein, over the wide range of pH from 3.0 to 8.8, exhibiting high stability in the pH region of 5.0-7.0. The corresponding maximum stabilities, DeltaG(25), were obtained at pH 7.0 with values of 15.4 kcal mol(-1) (combined fluorescence and circular dichroism data), and 15.1 kcal mol(-1) (DSC), considering a DeltaC(p) of 1.72 +/- 0.24 kcal mol(-1) K(-1). The low histidine content ( approximately 1.7%) and the high acidic residue content ( approximately 22.5%) suggests a flat pH dependence of thermal stability in the region 2.0-8.8 and that the decrease in thermal stability at low pH can be due to the differences in pK values of the acidic groups.

Computational analysis of sequence selection mechanisms

Mechanisms leading to gene variations are responsible for the diversity of species and are important components of the theory of evolution. One constraint on gene evolution is that of protein foldability; the three-dimensional shapes of proteins must be thermodynamically stable. We explore the impact of this constraint and calculate properties of foldable sequences using 3660 structures from the Protein Data Bank. We seek a selection function that receives sequences as input, and outputs survival probability based on sequence fitness to structure. We compute the number of sequences that match a particular protein structure with energy lower than the native sequence, the density of the number of sequences, the entropy, and the "selection" temperature. The mechanism of structure selection for sequences longer than 200 amino acids is approximately universal. For shorter sequences, it is not. We speculate on concrete evolutionary mechanisms that show this behavior.

Experiment-guided thermodynamic simulations on reversible two-state proteins: implications for protein thermostability

Here, we perform protein thermodynamic simulations within a set of boundary conditions, effectively blanketing the experimental data. The thermodynamic parameters, melting temperature (TG), enthalpy change at the melting temperature (DeltaHG) and heat capacity change (DeltaCp) were systematically varied over the experimentally observed ranges for small single domain reversible two-state proteins. Parameter sets that satisfy the Gibbs-Helmholtz equation and yield a temperature of maximal stability (TS) around room temperature were selected. The results were divided into three categories by arbitrarily chosen TG ranges. The TG ranges in these categories correspond to typical values of the melting temperatures observed for the majority of the proteins from mesophilic, thermophilic and hyperthermophilic organisms. As expected, DeltaCp values tend to be high in mesophiles and low in hyperthermophiles. An increase in TG is accompanied by an up-shift and broadening of the protein stability curves, however, with a large scatter. Furthermore, the simulations reveal that the average DeltaHG increases with TG up to approximately 360 K and becomes constant thereafter. DeltaCp decreases with TG with different rates before and after approximately 360 K. This provides further justification for the separate grouping of proteins into thermophiles and hyperthermophiles to assess their thermodynamic differences. This analysis of the Gibbs-Helmholtz equation has allowed us to study the interdependence of the thermodynamic parameters TG, DeltaHG and DeltaCp and their derivatives in a more rigorous way than possible by the limited experimental protein thermodynamics data available in the literature. The results provide new insights into protein thermostability and suggest potential strategies for its manipulation.

Thermodynamic analysis of the unfolding and stability of the dimeric DNA-binding protein HU from the hyperthermophilic eubacterium Thermotoga maritima and its E34D mutant

We have studied the stability of the histone-like, DNA-binding protein HU from the hyperthermophilic eubacterium Thermotoga maritima and its E34D mutant by differential scanning microcalorimetry and CD under acidic conditions at various concentrations within the range of 2-225 micro m of monomer. The thermal unfolding of both proteins is highly reversible and clearly follows a two-state dissociation/unfolding model from the folded, dimeric state to the unfolded, monomeric one. The unfolding enthalpy is very low even when taking into account that the two disordered DNA-binding arms probably do not contribute to the cooperative unfolding, whereas the quite small value for the unfolding heat capacity change (3.7 kJ.K(-1).mol(-1)) stabilizes the protein within a broad temperature range, as shown by the stability curves (Gibbs energy functions vs. temperature), even though the Gibbs energy of unfolding is not very high either. The protein is stable at pH 4.00 and 3.75, but becomes considerably less so at pH 3.50 and below, to the point that a simple decrease in concentration will lead to unfolding of both the wild-type and the mutant protein at pH 3.50 and low temperatures. This indicates that various acid residues lose their charges leaving uncompensated positively charged clusters. The wild-type protein is more stable than its E34D mutant, particularly at pH 4.00 and 3.75 although less so at 3.50 (1.8, 1.6 and 0.6 kJ.mol(-1) at 25 degrees C for DeltaDeltaG at pH 4.00, 3.75 and 3.50, respectively), which seems to be related to the effect of a salt bridge between E34 and K13.

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.

Temperature dependence of the hydrophobic hydration and interaction of simple solutes: An examination of five popular water models

We examine five different popular rigid water models (SPC, SPCE, TIP3P, TIP4P, and TIP5P) using molecular dynamics simulations in order to investigate the hydrophobic hydration and interaction of apolar Lennard-Jones solutes as a function of temperature in the range between 275 and 375 K along the 0.1 MPa isobar. For all investigated models and state points we calculate the excess chemical potential for the noble gases and methane employing the Widom particle insertion technique. All water models exhibit too small hydration entropies, but show a clear hierarchy. TIP3P shows poorest agreement with experiment, whereas TIP5P is closest to the experimental data at lower temperatures and SPCE is closest at higher temperatures. As a first approximation, this behavior can be rationalized as a temperature shift with respect to the solvation behavior found in real water. A rescaling procedure inspired by the information theory model of Hummer et al. [Chem. Phys. 258, 349 (2000)] suggests that the different solubility curves for the different models and real water can be largely explained on the basis of the different density curves at constant pressure. In addition, the models that give a good representation of the water structure at ambient conditions (TIP5P, SPCE, and TIP4P) show considerably better agreement with the experimental data than the ones which exhibit less structured O-O correlation functions (SPC and TIP3P). In the second part of the paper we calculate the hydrophobic interaction between xenon particles directly from a series of 60 ns simulation runs. We find that the temperature dependence of the association is to a large extent related to the strength of the solvation entropy. Nevertheless, differences between the models seem to require a more detailed molecular picture. The TIP5P model shows by far the strongest temperature dependence. The suggested density rescaling is also applied to the chemical potential in the xenon-xenon contact-pair configuration, indicating the presence of a temperature where the hydrophobic interaction turns into purely repulsive. The predicted association for xenon in real water suggests the presence of a strong variation with temperature, comparable to the behavior found for TIP5P water. Comparing different water models and experimental data we conclude that a proper description of density effects is an important requirement for a water model to account correctly for the correct description of the hydrophobic effects. A water model exhibiting a density maximum at the correct temperature is desirable.