Thermal unfolding of the DNA-binding protein Sso7d from the hyperthermophile Sulfolobus solfataricus

Beyond antibody engineering: directed evolution of alternative binding scaffolds and enzymes using yeast surface display

Pioneered exactly 20 years ago, yeast surface display (YSD) continues to take a major role in protein engineering among the high-throughput display methodologies that have been developed to date. The classical yeast display technology relies on tethering an engineered protein to the cell wall by genetic fusion to one subunit of a dimeric yeast-mating agglutination receptor complex. This method enables an efficient genotype–phenotype linkage while exploiting the benefits of a eukaryotic expression machinery. Over the past two decades, a plethora of protein engineering efforts encompassing conventional antibody Fab and scFv fragments have been reported. In this review, we will focus on the versatility of YSD beyond conventional antibody engineering and, instead, place the focus on alternative scaffold proteins and enzymes which have successfully been tailored for purpose with regard to improving binding, activity or specificity.

DNA Modulates the Interaction of Genetically Engineered DNA-Binding Proteins and Gold Nanoparticles: Diagnosis of High-Risk HPV Infection

Gene detection has an important role in diagnosing several serious diseases and genetic defects in modern clinical medicine. Herein, we report a fast and convenient gene detection method based on the modulation of the interaction between a heat-resistant double-stranded DNA (dsDNA)-binding protein (Sso7d) and gold nanoparticles (Au NPs). We prepared a recombinant Cys-Sso7d, which is Sso7d with an extra cysteine (Cys) residue in the N-terminus, through protein engineering to control the interaction between Sso7d and Au NPs. Cys-Sso7d exhibited a stronger affinity for Au NPs and more easily induced the aggregation of Au NPs than Sso7d. In addition, Cys-Sso7d retained its ability to bind with dsDNA. The aggregation of Au NPs induced by Cys-Sso7d was diminished in the presence of dsDNA, which could be utilized as a transduction mechanism for the detection of the polymerase chain reaction (PCR) products of human papillomavirus (HPV) gene fragments (HPV types 16 and 18). The Cys-Sso7d/Au NP probe could detect as few as 1 copy of the HPV gene. The sensitivity and specificity of the Cys-Sso7d/Au NP probe for Pap smear clinical specimens (n = 52) for HPV 16 and HPV 18 detection were 85.7%/100.0% and 85.7%/91.7%, respectively. Our results demonstrate that the Cys-Sso7d/Au NP probe can be used to diagnose high-risk HPV types in Pap smear samples with high sensitivity, specificity, and accuracy.

Protein folding at extreme temperatures: Current issues

The range of temperatures compatible with life is currently estimated from -25°C, as exemplified by metabolically active bacteria between sea ice crystals, and up to 122°C in hydrothermal vents as exemplified by the archaeon Methanopyrus kandleri. In the context of protein folding, as soon as a polypeptide emerges from the ribosome, it is exposed to the effects of environmental temperatures. Recent investigations have shown that the rate of protein folding is not adapted to extreme temperatures and should be very fast at high temperature and low in cold environments. This lack of adaptation is driven by kinetic constraints on protein stability. To counteract the deleterious effects of fast protein folding in hyperthermophiles, chaperones such as the Trigger Factor hold and slow down the rate of folding intermediates. Prolyl isomerization, a rate-limiting step in the folding of many proteins, is strongly temperature-dependent and impairs folding of psychrophilic proteins in the cold. This is compensated by reduction of the proline content in cold-adapted proteins, by an increased number of prolyl isomerases encoded in the genome of psychrophilic microorganisms and by overexpression of prolyl isomerases under low temperature cultivation. After folding, the native state is reached and although extremophilic proteins share the same fold, dramatic differences in stability have been recorded by differential scanning calorimetry.

Thermodynamics of globular proteins

The analysis of temperature-induced unfolding of proteins in aqueous solutions was performed. Based on the data of thermodynamic parameters of protein unfolding and using the method of semi-empirical calculations of hydration parameters at reference temperature 298 K, we obtained numerical values of enthalpy, free energy, and entropy which characterize the unfolding of proteins in the 'gas phase'. It was shown that specific values of the energy of weak intramolecular bonds (∆H^int), conformational free energy (∆G^conf) and entropy (∆S^conf) are the same for proteins with molecular weight 7-25 kDa. Using the energy value (∆H^int) and the proposed approach for estimation of the conformational entropy of native protein (S^NC), numerical values of the absolute free energy (G^NC) were obtained.

The archaeal "7 kDa DNA-binding" proteins: extended characterization of an old gifted family

The “7 kDa DNA-binding” family, also known as the Sul7d family, is composed of chromatin proteins from the Sulfolobales archaeal order. Among them, Sac7d and Sso7d have been the focus of several studies with some characterization of their properties. Here, we studied eleven other proteins alongside Sac7d and Sso7d under the same conditions. The dissociation constants of the purified proteins for binding to double-stranded DNA (dsDNA) were determined in phosphate-buffered saline at 25 °C and were in the range from 11 μM to 22 μM with a preference for G/C rich sequences. In accordance with the extremophilic origin of their hosts, the proteins were found highly stable from pH 0 to pH 12 and at temperatures from 85.5 °C to 100 °C. Thus, these results validate eight putative “7 kDa DNA-binding” family proteins and show that they behave similarly regarding both their function and their stability among various genera and species. As Sac7d and Sso7d have found numerous uses as molecular biology reagents and artificial affinity proteins, this study also sheds light on even more attractive proteins that will facilitate engineering of novel highly robust reagents.

Strong Enrichment of Aromatic Residues in Binding Sites from a Charge-neutralized Hyperthermostable Sso7d Scaffold Library

The Sso7d protein from the hyperthermophilic archaeon Sulfolobus solfataricus is an attractive binding scaffold because of its small size (7 kDa), high thermal stability (T_m of 98 °C), and absence of cysteines and glycosylation sites. However, as a DNA-binding protein, Sso7d is highly positively charged, introducing a strong specificity constraint for binding epitopes and leading to nonspecific interaction with mammalian cell membranes. In the present study, we report charge-neutralized variants of Sso7d that maintain high thermal stability. Yeast-displayed libraries that were based on this reduced charge Sso7d (rcSso7d) scaffold yielded binders with low nanomolar affinities against mouse serum albumin and several epitopes on human epidermal growth factor receptor. Importantly, starting from a charge-neutralized scaffold facilitated evolutionary adaptation of binders to differentially charged epitopes on mouse serum albumin and human epidermal growth factor receptor, respectively. Interestingly, the distribution of amino acids in the small and rigid binding surface of enriched rcSso7d-based binders is very different from that generally found in more flexible antibody complementarity-determining region loops but resembles the composition of antibody-binding energetic hot spots. Particularly striking was a strong enrichment of the aromatic residues Trp, Tyr, and Phe in rcSso7d-based binders. This suggests that the rigidity and small size of this scaffold determines the unusual amino acid composition of its binding sites, mimicking the energetic core of antibody paratopes. Despite the high frequency of aromatic residues, these rcSso7d-based binders are highly expressed, thermostable, and monomeric, suggesting that the hyperstability of the starting scaffold and the rigidness of the binding surface confer a high tolerance to mutation.

aKMT Catalyzes Extensive Protein Lysine Methylation in the Hyperthermophilic Archaeon Sulfolobus islandicus but is Dispensable for the Growth of the Organism

Protein methylation is believed to occur extensively in creanarchaea. Recently, aKMT, a highly conserved crenarchaeal protein lysine methyltransferase, was identified and shown to exhibit broad substrate specificity in vitro Here, we have constructed an aKMT deletion mutant of the hyperthermophilic crenarchaeon Sulfolobus islandicus The mutant was viable but showed a moderately slower growth rate than the parental strain under non-optimal growth conditions. Consistent with the moderate effect of the lack of aKMT on the growth of the cell, expression of a small number of genes, which encode putative functions in substrate transportation, energy metabolism, transcriptional regulation, stress response proteins, etc, was differentially regulated by more than twofold in the mutant strain, as compared with that in the parental strain. Analysis of the methylation of total cellular protein by mass spectrometry revealed that methylated proteins accounted for ∼2/3 (1,158/1,751) and ∼1/3 (591/1,757) of the identified proteins in the parental and the mutant strains, respectively, indicating that there is extensive protein methylation in S. islandicus and that aKMT is a major protein methyltransferase in this organism. No significant sequence preference was detected at the sites of methylation by aKMT. Methylated lysine residues, when visible in the structure, are all located on the surface of the proteins. The crystal structure of aKMT in complex with S-adenosyl-l-methionine (SAM) or S-adenosyl homocysteine (SAH) reveals that the protein consists of four α helices and seven β sheets, lacking a substrate recognition domain found in PrmA, a bacterial homolog of aKMT, in agreement with the broad substrate specificity of aKMT. Our results suggest that aKMT may serve a role in maintaining the methylation status of cellular proteins required for the efficient growth of the organism under certain non-optimal conditions.

Post-translational methylations of the archaeal Mre11:Rad50 complex throughout the DNA damage response

The Mre11:Rad50 complex is central to DNA double strand break repair in the Archaea and Eukarya, and acts through mechanical and nuclease activities regulated by conformational changes induced by ATP binding and hydrolysis. Despite the widespread use of Mre11 and Rad50 from hyperthermophilic archaea for structural studies, little is known in the regulation of these proteins in the Archaea. Using purification and mass spectrometry approaches allowing nearly full sequence coverage of both proteins from the species Sulfolobus acidocaldarius, we show for the first time post-translational methylation of the archaeal Mre11:Rad50 complex. Under basal growth conditions, extensive lysine methylations were identified in Mre11 and Rad50 dynamic domains, as well as methylation of a few aspartates and glutamates, including a key Mre11 aspartate involved in nuclease activity. Upon γ-irradiation induced DNA damage, additional methylated residues were identified in Rad50, notably methylation of Walker B aspartate and glutamate residues involved in ATP hydrolysis. These findings strongly suggest a key role for post-translational methylation in the regulation of the archaeal Mre11:Rad50 complex and in the DNA damage response.

The functional diversity of protein lysine methylation

Large‐scale characterization of post‐translational modifications (PTMs), such as phosphorylation, acetylation and ubiquitination, has highlighted their importance in the regulation of a myriad of signaling events. While high‐throughput technologies have tremendously helped cataloguing the proteins modified by these PTMs, the identification of lysine‐methylated proteins, a PTM involving the transfer of one, two or three methyl groups to the ε‐amine of a lysine side chain, has lagged behind. While the initial findings were focused on the methylation of histone proteins, several studies have recently identified novel non‐histone lysine‐methylated proteins. This review provides a compilation of all lysine methylation sites reported to date. We also present key examples showing the impact of lysine methylation and discuss the circuitries wired by this important PTM.

The entropic nature of protein thermal stabilization

We performed thermodynamic analysis of temperature-induced unfolding of mesophilic and thermophilic proteins. It was shown that the variability in protein thermostability associated with pH-dependent unfolding or linked to the substitution of amino acid residues on the protein surface is evidence of the governing role of the entropy factor. Numerical values of conformational components in enthalpy, entropy and free energy which characterize protein unfolding in the "gas phase" were obtained. Based on the calculated absolute values of entropy and free energy, a model of protein unfolding is proposed in which the driving force is the conformational entropy of native protein, as an energy of the heat motion (T·S(NC)) increasing with temperature and acting as an factor devaluating the energy of intramolecular weak bonds in the transition state.

Avidity-mediated virus separation using a hyperthermophilic affinity ligand

Immunoaffinity separation of large multivalent species such as viruses is limited by the stringent elution conditions necessary to overcome their strong and highly avid interaction with immobilized affinity ligands on the capture surface. Here we present an alternate strategy that harnesses the avidity effect to overcome this limitation. Red clover necrotic mosaic virus (RCNMV), a plant virus relevant to drug delivery applications, was chosen as a model target for this study. An RCNMV binding protein (RBP) with modest binding affinity (K(D) ~100 nM) was generated through mutagenesis of the Sso7d protein from Sulfolobus solfataricus and used as the affinity ligand. In our separation scheme, RCNMV is captured by a highly avid interaction with RBP immobilized on a nickel surface through a hexahistidine (6xHis) tag. Subsequently, disruption of the multivalent interaction and release of RCNMV is achieved by elution of RBP from the nickel surface. Finally, RCNMV is separated from RBP by exploiting the large difference in their molecular weights (~8 MDa vs. ~10 kDa). Our strategy not only eliminates the need for harsh elution conditions, but also bypasses chemical conjugation of the affinity ligand to the capture surface. Stable non-antibody affinity ligands to a wide spectrum of targets can be generated through mutagenesis of Sso7d and other hyperthermophilic proteins. Therefore, our approach may be broadly relevant to cases where capture of large multivalent species from complex mixtures and subsequent release without the use of harsh elution conditions is necessary.

Insights into the thermal stabilization and conformational transitions of DNA by hyperthermophile protein Sso7d: molecular dynamics simulations and MM-PBSA analysis

In the assembly of DNA-protein complex, the DNA kinking plays an important role in nucleoprotein structures and gene regulation. Molecular dynamics (MD) simulations were performed on specific protein-DNA complexes in this study to investigate the stability and structural transitions of DNA depending on temperature. Furthermore, we introduced the molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) approach to analyze the interactions between DNA and protein in hyperthermophile. Focused on two specific Sso7d-DNA complexes (PDB codes: 1BNZ and 1BF4), we performed MD simulations at four temperatures (300, 360, 420, and 480 K) and MM-PBSA at 300 and 360 K to illustrate detailed information on the changes of DNA. Our results show that Sso7d stabilizes DNA duplex over a certain temperature range and DNA molecules undergo B-like to A-like form transitions in the binary complex with the temperature increasing, which are consistent with the experimental data. Our work will contribute to a better understanding of protein-DNA interaction.

Thermal stability and unfolding pathways of Sso7d and its mutant F31A: insight from molecular dynamics simulation

The thermo-stability and unfolding behaviors of a small hyperthermophilic protein Sso7d as well as its single-point mutation F31A are studied by molecular dynamics simulation at temperatures of 300 K, 371 K and 500 K. Simulations at 300 K show that the F31A mutant displays a much larger flexibility than the wild type, which implies that the mutation obviously decreases the protein's stability. In the simulations at 371 K, although larger fluctuations were observed, both of these two maintain their stable conformations. High temperature simulations at 500 K suggest that the unfolding of these two proteins evolves along different pathways. For the wild-type protein, the C-terminal alpha-helix is melted at the early unfolding stage, whereas it is destroyed much later in the unfolding process of the F31A mutant. The results also show that the mutant unfolds much faster than its parent protein. The deeply buried aromatic cluster in the F31A mutant dissociates quickly relative to the wild-type protein at high temperature. Besides, it is found that the triple-stranded antiparallel β-sheet in the wild-type protein plays an important role in maintaining the stability of the entire structure.

Highly stable binding proteins derived from the hyperthermophilic Sso7d scaffold

We have shown that highly stable binding proteins for a wide spectrum of targets can be generated through mutagenesis of the Sso7d protein from the hyperthermophilic archaeon Sulfolobus solfataricus. Sso7d is a small (~7 kDa, 63 amino acids) DNA-binding protein that lacks cysteine residues and has a melting temperature of nearly 100 °C. We generated a library of 10(8) Sso7d mutants by randomizing 10 amino acid residues on the DNA-binding surface of Sso7d, using yeast surface display. Binding proteins for a diverse set of model targets could be isolated from this library; our chosen targets included a small organic molecule (fluorescein), a 12 amino acid peptide fragment from the C-terminus of β-catenin, the model proteins hen egg lysozyme and streptavidin, and immunoglobulins from chicken and mouse. Without the application of any affinity maturation strategy, the binding proteins isolated had equilibrium dissociation constants in the nanomolar to micromolar range. Further, Sso7d-derived binding proteins could discriminate between closely related immunoglobulins. Mutant proteins based on Sso7d were expressed at high yields in the Escherichia coli cytoplasm. Despite extensive mutagenesis, Sso7d mutants have high thermal stability; five of six mutants analyzed have melting temperatures >89 °C. They are also resistant to chemical denaturation by guanidine hydrochloride and retain their secondary structure after extended incubation at extreme pH values. Because of their favorable properties, such as ease of recombinant expression, and high thermal, chemical and pH stability, Sso7d-derived binding proteins will have wide applicability in several areas of biotechnology and medicine.

Thermodynamic basis for the stabilities of three CutA1s from Pyrococcus horikoshii,Thermus thermophilus, and Oryza sativa, with unusually high denaturation temperatures

In order to elucidate the stabilization mechanism of CutA1 from Pyrococcus horikoshii (PhCutA1) with a denaturation temperature of nearly 150 degrees C, GuHCl denaturation and heat denaturation were examined at neutral and acidic pHs. As a comparison, CutA1 proteins from Thermus thermophilus (TtCutA1) and Oryza sativa (OsCutA1) were also examined, which have lower optimum growth temperatures of 75 and 28 degrees C, respectively, than that (98 degrees C) of P. horikoshii. GuHCl-induced unfolding and refolding curves of the three proteins showed hysteresis effects due to an unusually slow unfolding rate. The midpoints of refolding for PhCutA1, TtCutA1 and OsCutA1 were 5.7 M, 3.3 M, and 2.3 M GuHCl, respectively, at pH 8.0 and 37 degrees C. DSC experiments with TtCutA1 and OsCutA1 showed that the denaturation temperatures were remarkably high, 112.8 and 97.3 degrees C, respectively, at pH 7.0 and that the good heat reversibility was amenable to thermodynamic analyses. At acidic pH, TtCutA1 showed higher stability to both heat and denaturant than PhCutA1. Combined with the data for DSC and denaturant denaturation, the unfolding Gibbs energy of PhCutA1 could be depicted as a function of temperature. It was experimentally revealed that (1) the unusually high stability of PhCutA1 basically originates from a common trimer structure of the three proteins, (2) the stability of PhCutA1 is superior to those of the other two CutA1s over all temperatures above 0 degrees C at neutral pH, due to the decrease in both enthalpy and entropy, and (3) ion pairs of PhCutA1 contribute to the unusually high stability at neutral pH.

Effect of finite size on cooperativity and rates of protein folding

We analyze the dependence of cooperativity of the thermal denaturation transition and folding rates of globular proteins on the number of amino acid residues, N, using lattice models with side chains, off-lattice Go models, and the available experimental data. A dimensionless measure of cooperativity, Omega(c) (0 < Omega(c) < infinity), scales as Omega(c) approximately N(zeta). The results of simulations and the analysis of experimental data further confirm the earlier prediction that zeta is universal with zeta = 1 + gamma, where exponent gamma characterizes the susceptibility of a self-avoiding walk. This finding suggests that the structural characteristics in the denaturated state are manifested in the folding cooperativity at the transition temperature. The folding rates k(F) for the Go models and a dataset of 69 proteins can be fit using k(F) = k(F)0 exp(-cN(beta)). Both beta = 1/2 and 2/3 provide a good fit of the data. We find that k(F) = k(F)0 exp(-cN(1/2)), with the average (over the dataset of proteins) k(F)0 approximately (0.2 micros)(-1) and c approximately 1.1, can be used to estimate folding rates to within an order of magnitude in most cases. The minimal models give identical N dependence with c approximately 1. The prefactor for off-lattice Go models is nearly 4 orders of magnitude larger than the experimental value.

Thermochemistry of protein-DNA interaction studied with temperature-controlled nonequilibrium capillary electrophoresis of equilibrium mixtures

We introduce temperature-controlled nonequilibrium capillary electrophoresis of equilibrium mixtures (NECEEM) and demonstrate its use to study thermochemistry of protein-DNA interactions. Being a homogeneous kinetic method, temperature-controlled NECEEM uniquely allows finding temperature dependencies of equilibrium and kinetic parameters of complex formation without the immobilization of the interacting molecules on the surface of a solid substrate. In this work, we applied temperature-controlled NECEEM to study the thermochemistry of two protein-DNA pairs: (i) Taq DNA polymerase with its DNA aptamer and (ii) E. coli single-stranded DNA binding protein with a 20-base-long single-stranded DNA. We determined temperature dependencies of three parameters: the equilibrium binding constant (Kb), the rate constant of complex dissociation (k(off)), and the rate constant of complex formation (k(on)). The Kb(T) functions for both protein-DNA pairs had phase-transition-like points suggesting temperature-dependent conformational changes in structures of the interacting macromolecules. Temperature dependencies of k(on) and k(off) provided insights into how the conformational changes affected two opposite processes: binding and dissociation. Finally, thermodynamic parameters, DeltaH and DeltaS, for complex formation were found for different conformations. With its unique features and potential applicability to other macromolecular interactions, temperature-controlled NECEEM establishes a valuable addition to the arsenal of analytical methods used to study dynamic molecular complexes.

Lessons in stability from thermophilic proteins

Studies that compare proteins from thermophilic and mesophilic organisms can provide insights into ability of thermophiles to function at their high habitat temperatures and may provide clues that enable us to better define the forces that stabilize all proteins. Most of the comparative studies have focused on thermal stability and show, as expected, that thermophilic proteins have higher Tm values than their mesophilic counterparts. Although these comparisons are useful, more detailed thermodynamic analyses are required to reach a more complete understanding of the mechanisms thermophilic protein employ to remain folded over a wider range of temperatures. This complete thermodynamic description allows one to generate a stability curve for a protein that defines how the conformational stability (DeltaG) varies with temperature. Here we compare stability curves for many pairs of homologous proteins from thermophilic and mesophilc organisms. Of the basic methods that can be employed to achieve enhanced thermostability, we find that most thermophilic proteins use the simple method that raises the DeltaG at all temperatures as the principal way to increase their Tm. We discuss and compare this thermodynamic method with the possible alternatives. In addition we propose ways that structural alterations and changes to the amino acid sequences might give rise to varied methods used to obtain thermostability.

Slow irreversible unfolding of Pyrococcus furiosus triosephosphate isomerase: Separation and quantitation of conformers through a novel electrophoretic approach

The thermostability of hyperthermophile proteins is not easily studied because such proteins tend to be extremely recalcitrant to unfolding. Weeks of exposure to structurally destabilizing conditions are generally required to elicit any evidence of conformational change(s). The main reason for this extreme kinetic stability would appear to be the dominance of local unfolding transitions that occur within different parts of the structures of these molecules; put differently, local sub structural unfolding transitions that occur autonomously and reversibly are thought to fail to cooperate to bring about global unfolding in a facile manner, leading to a low overall observed rate of unfolding. For reasons that are not yet fully understood, unfolding is also reported to occur irreversibly in hyperthermophile proteins. Therefore, conventional experimental approaches are often unsuited to the study of their unfolding. Here, we describe a novel electrophoretic approach that facilitates separation, direct visualization, and quantitation of the folded, partially folded, and unfolded forms of the hyperthermophile protein triosephosphate isomerase from Pyrococcus furiosus, produced in the course of its irreversible structural destabilization by the combined action of heat and chemical agents. Our approach exploits (i) the irreversibility of global unfolding effected by heat and denaturants such as urea or guanidine hydrochloride, (ii) the stability of the native form of the protein to unfolding by the anionic detergent sodium dodecyl sulfate, (iii) the differential susceptibilities of various protein conformations to being bound by SDS, and (iv) the differential electrophoretic migration behavior displayed as a consequence of differential SDS binding.

Contribution of protein-surface ion pairs of a hyperthermophilic protein on thermal and thermodynamic stability

Hyperthermophilic proteins possess many ion pairs on their surface. To reveal the role of the ion pairs, O6-methylguanine-DNA methyltransferase from Thermococcus kodakaraensis KOD1 (Tk-MGMT) was studied as a model protein. The maximum free-energy changes of the protein in 0.1 and 0.5 M NaCl at pH 7.0 were 61.7 kJ mol(-1) at 31.5 degrees C and 77.4 kJ mol(-1) at 39.7 degrees C, respectively. On the other hand, mid points of the thermal unfolding temperatures in 0.1 and 0.5 M NaCl at pH 7.0 were 94.8 degrees C and 90.1 degrees C, respectively. The results suggest that the protein-surface ion pairs contribute to thermal stability (Tm), rather than thermodynamic stability (DeltaG).

Thermal and conformational stability of Ssh10b protein from archaeon Sulfolobus shibattae

The secondary structure of the DNA binding protein Ssh10b is largely unaffected by change in temperature between 25 degrees C and 85 degrees C, indicating that the protein is highly thermostable. Here, we report the temperature-dependent equilibrium denaturation of Ssh10b in the presence of guanidine hydrochloride (GdnHCl). It was found that the transition midpoint values of the temperature (T(m)), and changes of enthalpy (DeltaH(m)) and entropy (DeltaS(m)) of Ssh10b unfolding were linearly decreasing with increasing GdnHCl concentration. The true values of the thermodynamic parameters, T(m)=402 K, DeltaH(m)=590+/-40 kJ x mol(-1) and DeltaS(m)=1.4+/-0.15 kJ x T(-1) x mol(-1), were obtained by linear extrapolation to 0 M GdnHCl. The value of the heat capacity change of Ssh10b unfolding, DeltaC(p)=3.8+/-0.2 kJ x T(-1) x mol(-1) (approx. 19 J T(-1) x mol residue(-1)), was obtained from the measured thermodynamic parameters. This is significantly smaller than that of the average value for mesophilic proteins (50 J.K(-1) x mol residue(-1)) or the value calculated from the Ssh10b structural data (64 J T(-1) x mol residue(-1)). A consequence of the small DeltaC(p) is that the DeltaG of Ssh10b is larger than that of mesophilic proteins, while the values of DeltaH and T*DeltaS are smaller. The small DeltaC(p) of Ssh10b appears to result mainly from the presence of compactness in the denatured state.

Contribution of a Single-Turn α-Helix to the Conformational Stability and Activity of the Alkaline Proteinase Inhibitor of Pseudomonas aeruginosa

The alkaline proteinase inhibitor of Pseudomonas aeruginosa (APRin), a high-affinity inhibitor of the serralysin family of bacterial metalloproteinases, is folded into an eight-stranded beta-barrel with an N-terminal trunk linked to the barrel by a single-turn alpha-helix (helix A, residues 8-11). We show here that deletion or modification of helix A decreases the conformational stability of APRin as assessed by thermal and chemical denaturation with guanidinium chloride (GdmCl). The apparent melting temperature T(m) of the wild-type protein was 81.5 degrees C at pH 7.1 as assessed by circular dichroism and 87.5 degrees C by differential scanning calorimetry. Reduction of the single disulfide bond of APRin decreased T(m) by approximately 18 degrees C, while deletion of residues 6-10 or 1-10 lowered T(m) by approximately 8 and approximately 14 degrees C, respectively. DeltaG(u) as assessed by chemical denaturation was 7.2 kcal mol(-)(1) at 25 degrees C for wild-type APRin and was decreased by 3.4, 2.4, and 2.6 kcal mol(-)(1) by disulfide reduction, deletion of residues 6-10, and deletion of residues 1-10, respectively. In contrast, deletion of residues 1-5 had no significant effect on either T(m) or DeltaG(u). Substitution of five helix-breaking Gly or Pro residues in positions 6-10 as well as disruption of hydrogen bonds involving residues within helix A (mutants Asp10Pro and Trp15Phe) also decreased T(m) and DeltaG(u). The data suggest that a hydrogen-bonding network involving Leu11 in helix A and Trp15 located at the top of the barrel may prevent access of solvent to the interior of the barrel. Disruption of the helix could facilitate solvation of the nonpolar interior of the barrel, thereby destabilizing its folded structure. Kinetic studies with single amino acid mutants in helix A indicate that it modulates the affinity of APRin for APR primarily by influencing the dissociation rate of the inhibitor from the complex.

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.

Using protein folding rates to test protein folding theories

The fastest simple, kinetically two-state protein folds a million times more rapidly than the slowest. Here we review many recent theories of protein folding kinetics in terms of their ability to qualitatively rationalize, if not quantitatively predict, this fundamental experimental observation.

Insertion of the cytochrome b5 heme-binding loop into an SH3 domain. Effects on structure and stability, and clues about the cytochrome's architecture

Under native conditions, apocytochrome b(5) exhibits a stable core and a disordered heme-binding region that refolds upon association with the cofactor. The termini of this flexible region are in close proximity, suggesting that loop closure may contribute to the thermodynamic properties of the apocytochrome. A chimeric protein containing 43 residues encompassing the cytochrome loop was constructed using the cyanobacterial photosystem I accessory protein E (PsaE) from Synechococcus sp. PCC 7002 as a structured scaffold. PsaE has the topology of an SH3 domain, and the insertion was engineered to replace its 14-residue CD loop. NMR and optical spectroscopies showed that the hybrid protein (named EbE1) was folded under native conditions and that it retained the characteristics of an SH3 domain. NMR spectroscopy revealed that structural and dynamic differences were confined near the site of loop insertion. Variable-temperature 1D NMR spectra of EbE1 confirmed the presence of a kinetic unfolding barrier. Thermal and chemical denaturations of PsaE and EbE1 demonstrated cooperative, two-state transitions; the stability of the PsaE scaffold was found only moderately compromised by the insertion, with a DeltaT(m) of 8.3 degrees C, a DeltaC(m) of 1.5 M urea, and a DeltaDeltaG degrees of 4.2 kJ/mole. The data implied that the penalty for constraining the ends of the inserted region was lower than the approximately 6.4 kJ/mole calculated for a self-avoiding chain. Extrapolation of these results to cytochrome b(5) suggested that the intrinsic stability of the folded portion of the apoprotein reflected only a small detrimental contribution from the large heme-binding domain.

Unusually Slow Denaturation and Refolding Processes of Pyrrolidone Carboxyl Peptidase from a Hyperthermophile Are Highly Cooperative: Real-Time NMR Studies

The refolding rate of heat-denatured cysteine-free pyrrolidone carboxyl peptidase (PCP-0SH) from Pyrococcus furiosus has been reported to be unusually slow under some conditions. To elucidate the structural basis of the unusually slow kinetics of the protein, the denaturation and refolding processes of the PCP-0SH were investigated using a real-time 2D (1)H-(15)N HSQC and CD experiments. At 2 M urea denaturation of the PCP-0SH in the acidic region, all of the native peaks in the 2D HSQC spectrum completely disappeared. The conformation of the PCP-0SH just after removal of 6 M GuHCl could be observed as a stable intermediate (D(1) state) in 2D HSQC and CD experiments, which is similar to a molten globule structure. The D(1) state of the PCP-0SH, which is the initial state of refolding, corresponded to the state at 2 M urea and seemed to be the denatured state in equilibrium with the native state under the physiological conditions. The refolding of PCP-0SH from the D(1) state to the native state could be observed to be highly cooperative without any intermediates between them, even if the refolding rate was quite slow. In the higher concentration of denaturants, PCP-0SH showed HSQC and CD spectra characteristic of completely unfolded proteins called the D(2) state. The unusually slow refolding rate was discussed as originating in the conformations in the transition state and/or the retardation of reorganization in an ensemble of nonrandom denatured structures in the D(1) state.

Thermodynamics of core hydrophobicity and packing in the hyperthermophile proteins Sac7d and Sso7d

The influence of core hydrophobicity and packing on the structure and stability of the hyperthermophile proteins Sac7d and Sso7d have been studied by calorimetry, circular dichroism, and NMR. Valine 30 is positioned in Sac7d to allow a cavity-filling Val --> Ile substitution which occurs naturally in the homologous more thermostable Sso7d. The cavity-filling mutation in Sac7d has been characterized and compared to the reciprocal Ile --> Val mutation in Sso7d. A detailed analysis of the stability of the proteins was obtained by globally fitting the variation of DSC parameters and circular dichroism intensities as a function of temperature (0-100 degrees C), salt (0-0.3 M), and pH (0-8). A global analysis over such a range of conditions permitted an unusually precise measure of the thermodynamic parameters, as well as the separation of the thermodynamics of the intrinsic unfolding reaction from the linked effects of protonation and chloride binding associated with acid-induced folding. The results indicate differences in the energetics of unfolding Sac7d and Sso7d that would not be apparent from an analysis of DSC data alone using conventional methods. The sign and magnitude of the changes in DeltaG, DeltaH, TDeltaS, and DeltaC(P) of unfolding resulting from core Ile/Val substitutions in the two proteins were consistent with differences in hydrophobicity of Val and Ile and negligible changes in packing (van der Waals) interactions. The benefit of increased hydrophobicity of the core increased with temperature, with maximal effect around 116 degrees C. Increased hydrophobicity of the core achieved not only an increase in the free energy of unfolding, but also a lateral shift of the temperature of maximal stability to higher temperature.

Guanidine-induced unfolding of the Sso7d protein from the hyperthermophilic archaeon Sulfolobus solfataricus

The unfolding induced by guanidine hydrochloride of the small protein Sso7d from the hyperthermophilic archaeon Sulfolobus solfataricus has been investigated by means of circular dichroism and fluorescence measurements. At neutral pH and room temperature the midpoint of the transition occurred at 4M guanidine hydrochloride. Thermodynamic information was obtained by means of both the linear extrapolation model and the denaturant binding model, in the assumption of a two-state N<==>D transition. A comparison with thermodynamic data determined from the thermal unfolding of Sso7d indicated that the denaturant binding model has to be preferred. Finally, it is shown that Sso7d is the most stable against both temperature and guanidine hydrochloride among a set of globular proteins possessing a very similar 3D structure.

HIV-1 gp41 and gp160 are hyperthermostable proteins in a mesophilic environment. Characterization of gp41 mutants

HIV gp41(24-157) unfolds cooperatively over the pH range of 1.0-4.0 with T(m) values of > 100 degrees C. At pH 2.8, protein unfolding was 80% reversible and the DeltaH(vH)/DeltaH(cal) ratio of 3.7 is indicative of gp41 being trimeric. No evidence for a monomer-trimer equilibrium in the concentration range of 0.3-36 micro m was obtained by DSC and tryptophan fluorescence. Glycosylation of gp41 was found to have only a marginal impact on the thermal stability. Reduction of the disulfide bond or mutation of both cysteine residues had only a marginal impact on protein stability. There was no cooperative unfolding event in the DSC thermogram of gp160 in NaCl/P(i), pH 7.4, over a temperature range of 8-129 degrees C. When the pH was lowered to 5.5-3.4, a single unfolding event at around 120 degrees C was noted, and three unfolding events at 93.3, 106.4 and 111.8 degrees C were observed at pH 2.8. Differences between gp41 and gp160, and hyperthermostable proteins from thermophile organisms are discussed. A series of gp41 mutants containing single, double, triple or quadruple point mutations were analysed by DSC and CD. The impact of mutations on the protein structure, in the context of generating a gp41 based vaccine antigen that resembles a fusion intermediate state, is discussed. A gp41 mutant, in which three hydrophobic amino acids in the gp41 loop were replaced with charged residues, showed an increased solubility at neutral pH.

The role of the salt concentration, proton, and phosphate binding on the thermal stability of wild and cloned DNA-binding protein Sso7d from Sulfolobus solfataricus

The acidic pH (1.5-7.0) and ionic strength (0.005-0.2M) dependence of thermodynamic functions of protein Sso7d from Sulfolobus solfataricus, cloned (c-Sso7d) and N-heptapeptide deleted [c-des(1-7)Sso7d] in glycine, and phosphate buffers was studied by means of adiabatic scanning calorimetry. The difference of proton binding was estimated from deltaHcal(pH), Td(pH), and (deltaTd/deltapH). It was found that a single group non-co-operative ionization with apparent pKa = 3.25 for both cloned and deleted proteins govern the thermal unfolding of two different (protonated and unprotonated) forms. deltaH degrees is found to be pH-independent and the changes in stability (deltaG degrees ) originate from changes in entropy terms. The apparent pKa measured at high salt concentrations decreases with 0.5 pH units from glycine to phosphate and the free energy of transfer at high ionic strength is 0.7 kcal/mol. The ionic strength dependence for the pH-dependent D-states is very different at pH 6.0 and 1.5. This is consistent with the property of denatured state to be more compacted or "closed" (Dc) at neutral or weak acidic pH and more random or "open" (Do) at acidic pH. From the Bjerrum's relation was found the number of screened charges important for the unfolding process. The main conclusions are: (1) the thermal stability of Sso7d has prominently entropic nature; (2) a single non-co-operative ionization controls the conformations in the D-state; and (3) pH-dependent conformational equilibrium could be functionally important in Sso7d-DNA recognition.

A novel strategy to engineer DNA polymerases for enhanced processivity and improved performance in vitro

Mechanisms that allow replicative DNA polymerases to attain high processivity are often specific to a given polymerase and cannot be generalized to others. Here we report a protein engineering-based approach to significantly improve the processivity of DNA polymerases by covalently linking the polymerase domain to a sequence non-specific dsDNA binding protein. Using Sso7d from Sulfolobus solfataricus as the DNA binding protein, we demonstrate that the processivity of both family A and family B polymerases can be significantly enhanced. By introducing point mutations in Sso7d, we show that the dsDNA binding property of Sso7d is essential for the enhancement. We present evidence supporting two novel conclusions. First, the fusion of a heterologous dsDNA binding protein to a polymerase can increase processivity without compromising catalytic activity and enzyme stability. Second, polymerase processivity is limiting for the efficiency of PCR, such that the fusion enzymes exhibit profound advantages over unmodified enzymes in PCR applications. This technology has the potential to broadly improve the performance of nucleic acid modifying enzymes.

Structural and dynamic effects of α-Helix deletion in Sso7d: Implications for protein thermal stability

Sso7d is a 62-residue protein from the hyperthemophilic archaeon Sulfolobus solfataricus with a denaturation temperature close to 100 degrees C around neutral pH. An engineered form of Sso7d truncated at leucine 54 (L54Delta) is significantly less stable, with a denaturation temperature of 53 degrees C. Molecular dynamics (MD) studies of Sso7d and its truncated form at two different temperatures have been performed. The results of the MD simulations at 300 K indicate that: (1) the flexibility of Sso7d chain at 300 K agrees with that detected from X-ray and NMR structural studies; (2) L54Delta remains stable in the native folded conformation and possesses an overall dynamic behavior similar to that of the parent protein. MD simulations performed at 500 K, 10 ns long, indicate that, while Sso7d is in-silico resistant to high temperature, the truncated variant partially unfolds, revealing the early phases of the thermal unfolding pathway of the protein. Analysis of the trajectories of L54Delta suggests that the unzipping of the N-terminal and C-terminal beta-strands should be the first event of the unfolding pathway, and points out the regions more resistant to thermal unfolding. These findings allow one to understand the role played by specific interactions connecting the two ends of the chain for the high thermal stability of Sso7d, and support recent hypotheses on its folding mechanism emerged from site-directed mutagenesis studies.

Thermal stability and DNA binding activity of a variant form of the Sso7d protein from the archeon Sulfolobus solfataricus truncated at leucine 54

Sso7d is a 62-residue, basic protein from the hyperthermophilic archaeon Sulfolobus solfataricus. Around neutral pH, it exhibits a denaturation temperature close to 100 degrees C and a non-sequence-specific DNA binding activity. Here, we report the characterization by circular dichroism and fluorescence measurements of a variant form of Sso7d truncated at leucine 54 (L54Delta). It is shown that L54Delta has a folded conformation at neutral pH and that its thermal unfolding is a reversible process, represented well by the two-state N <=> D transition model, with a denaturation temperature of 53 degrees C. Fluorescence titration experiments indicate that L54Delta binds tightly to calf thymus DNA, even though the binding parameters are smaller than those of the wild-type protein. Therefore, the truncation of eight residues at the C-terminus of Sso7d markedly affects the thermal stability of the protein, which nevertheless retains a folded structure and DNA binding activity.

Crystal structure of the hyperthermophilic archaeal DNA-binding protein Sso10b2 at a resolution of 1.85 Angstroms

The crystal structure of a small, basic DNA binding protein, Sso10b2, from the thermoacidophilic archaeon Sulfolobus solfataricus was determined by the Zn multiwavelength anomalous diffraction method and refined to 1.85 A resolution. The 89-amino-acid protein adopts a betaalphabetaalphabetabeta topology. The structure is similar to that of Sso10b1 (also called Alba) from the same organism. However, Sso10b2 contains an arginine-rich loop RDRRR motif, which may play an important role in nucleic acid binding. There are two independent Sso10b2 proteins in the asymmetric unit, and a plausible stable dimer could be deduced from the crystal structure. Topology comparison revealed that Sso10b2 is similar to several RNA-binding proteins, including IF3-C, YhhP, and DNase I. Models of the Sso10b2 dimer bound to either B-DNA or A-DNA have been constructed.

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

The difference between the heat (T(G)) and the cold (T(G)') denaturation temperatures defines the temperature range (T(Range)) over which the native state of a reversible two-state protein is thermodynamically stable. We have performed a correlation analysis for thermodynamic parameters in a selected data set of structurally nonhomologous single-domain reversible two-state proteins. We find that the temperature range is negatively correlated with the protein size and with the heat capacity change (DeltaC(p)) but is positively correlated with the maximal protein stability [DeltaG(T(S))]. The correlation between the temperature range and maximal protein stability becomes highly significant upon normalization of the maximal protein stability with protein size. The melting temperature (T(G)) also shows a negative correlation with protein size. Consistently, T(G) and T(G)' show opposite correlations with DeltaC(p), indicating a dependence of the T(Range) on the curvature of the protein stability curve. Substitution of proteins in our data set with their homologues and arbitrary addition or removal of a protein in the data set do not affect the outcome of our analysis. Simulations of the thermodynamic data further indicate that T(Range) is more sensitive to variations in curvature than to the slope of the protein stability curve. The hydrophobic effect in single domains is the principal reason for these observations. Our results imply that larger proteins may be stable over narrower temperature ranges and that smaller proteins may have higher melting temperatures, suggesting why protein structures often differentiate into multiple substructures with different hydrophobic cores. Our results have interesting implications for protein thermostability.

Toward the physical basis of thermophilic proteins: linking of enriched polar interactions and reduced heat capacity of unfolding

The enrichment of salt bridges and hydrogen bonding in thermophilic proteins has long been recognized. Another tendency, featuring lower heat capacity of unfolding (DeltaC(p)) than found in mesophilic proteins, is emerging from the recent literature. Here we present a simple electrostatic model to illustrate that formation of a salt-bridge or hydrogen-bonding network around an ionized group in the folded state leads to increased folding stability and decreased DeltaC(p). We thus suggest that the reduced DeltaC(p) of thermophilic proteins could partly be attributed to enriched polar interactions. A reduced DeltaC(p) might serve as an indicator for the contribution of polar interactions to folding stability.

Maximal stabilities of reversible two-state proteins

The hydrophobic effect is the major force driving protein folding. Around room temperature, small organic solutes and hydrophobic amino acids have low solubilities in water and the hydrophobic effect is the strongest. These facts suggest that globular proteins should be maximally stable around room temperature. While this fundamental paradigm has been expected, it has not actually been shown to hold. Toward this goal, we have collected and analyzed experimental thermodynamic data for 31 proteins that show reversible two-state folding <--> unfolding transitions at or near neutral pH. Twenty-six of these are unique, and 20 of the 26 are maximally stable around room temperature irrespective of their structural properties, the melting temperature, or the living temperatures of their source organisms. Their average temperature of maximal stability is 293 +/- 8 K (20 +/- 8 degrees C). These proteins differ in size, fold, and number of domains, hydrophobic folding units, and oligomeric states. They derive from the cold-loving psychrophiles, from mesophiles, and from thermophiles. Analysis of the single-domain proteins present in this set shows that the variations in their thermodynamic parameters are correlated in a way which may explain the adaptation of the proteins to the living temperatures of the organisms from which they derive. The average energetic contribution of the individual amino acids toward protein stability decreases with an increase in protein size, suggesting that there may be an upper limit for protein maximal thermodynamic stability. For the remaining proteins, deviation of the maximal stability temperatures from room temperature may be due to greater uncertainties in their heat capacity change (DeltaC(p)) values, a weaker hydrophobic effect, and/or a stronger electrostatic contribution.

The unusually slow relaxation kinetics of the folding-unfolding of pyrrolidone carboxyl peptidase from a hyperthermophile, Pyrococcus furiosus

In order to understand the thermodynamic and kinetic basis of the intrinsic stability of proteins from hyperthermophiles, the folding-unfolding reactions of cysteine-free pyrrolidone carboxyl peptidase (Cys142/188Ser) (PCP-0SH) from Pyrococcus furiosus were examined using circular dichroism (CD) and differential scanning calorimetry (DSC) at pH 2.3, where PCP-0SH exists in monomeric form. DSC showed a strong dependence of the shape and position of the unfolding profiles on the scan rate, suggesting the stability of PCP-0SH under kinetic control. On DSC timescales, even at a scan rate of 1 deg. C/hour, heat denaturation of PCP-0SH was non-equilibrium. However, over a long period of incubation of the heat-denatured PCP-0SH at pre-transition temperatures, it refolded completely, indicating reversibility with very slow relaxation kinetics. The rates of refolding of the heat-denatured PCP-0SH determined from the time-resolved DSC and CD spectroscopic progress curves were found to be similar within experimental error, confirming the mechanism of refolding to be a two-state process. The equilibrium established with a relaxation time of 5080 seconds (at t(m)=46.5 degrees C), which is unusually higher than the relaxation times observed for mesophilic and hyperthermophilic proteins. The long relaxation time may lead to the apparent irreversibility of an unfolding process occurring on the DSC experiment timescale. The refolding rate (9.8 x 10(-5) s(-1)) peaked near the t(m) (=46.5 degrees C), whereas the stability profile reached maxima (11.8 kJ mol(-1)) at 17 degrees C. The results clearly indicate the unusual mode of protein destabilization via a drastic decrease in the rate of folding at low pH and still maintaining a high activation energy barrier (284 kJ mol(-1)) for unfolding, which provides an effective kinetic advantage to unusually stable proteins from hyperthermophiles.

High thermal stability of 3-isopropylmalate dehydrogenase from Thermus thermophilus resulting from low DeltaC(p) of unfolding

To characterize the thermal stability of 3-isopropylmalate dehydrogenase (IPMDH) from an extreme thermophile, Thermus thermophilus, urea-induced unfolding of the enzyme and of its mesophilic counterpart from Escherichia coli was investigated at various temperatures. The unfolding curves were analyzed with a three-state model for E.coli IPMDH and with a two-state model for T.thermophilus IPMDH, to obtain the free energy change DeltaG degrees of each unfolding process. Other thermodynamic parameters, enthalpy change DeltaH, entropy change DeltaS and heat capacity change DeltaC(p), were derived from the temperature dependence of DeltaG degrees. The main feature of the thermophilic enzyme was its lower dependence of DeltaG degrees on temperature resulting from a low DeltaC(p). The thermophilic IPMDH had a significantly lower DeltaC(p), 1.73 kcal/mol.K, than that of E.coli IPMDH (20.7 kcal/mol.K). The low DeltaC(p) of T.thermophilus IPMDH could not be predicted from its change in solvent-accessible surface area DeltaASA. The results suggested that there is a large structural difference between the unfolded state of T.thermophilus and that of E.coli IPMDH. Another responsible factor for the higher thermal stability of T.thermophilus IPMDH was the increase in the most stable temperature T(s). The DeltaG degrees maximum of T.thermophilus IPMDH was much smaller than that of E.coli IPMDH. The present results clearly demonstrated that a higher melting temperature T(m) is not necessarily accompanied by a higher DeltaG degrees maximum.