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

Essentiality of core hydrophobicity to the structure and function of archaeal chromatin protein Cren7

Studies on the structure-function relationship of protein greatly help to understand not only the principles of protein folding but also the rationales of protein engineering. Crenarchaeal chromatin protein Cren7 provides an excellent research model for this issue. The small protein adopts a 'β-barrel' fold, formed by the double-stranded antiparallel β-sheet 1 tightly packing with the triple-stranded antiparallel β-sheet 2. The simple structure of Cren7 is stabilized by the hydrophobic core between the β-sheets, consisting of the side chains of V8, V10, L20, V25, F41 and F50. In the present work, mutation analyses by alanine substitution of each of the residues in the hydrophobic core were performed. Circular dichroism spectra and nuclear magnetic resonance analyses showed that mutation of F41 led to a significant misfolding of Cren7 through disruption of the β-sheets. Meanwhile, the mutant F41A showed a reduced thermostatility (Tm of 53.2 °C), as compared with the wild-type Cren7 (Tm > 80 °C). Biolayer interferometry and nick-closure assays showed the largely unchanged activities in DNA binding and supercoiling of F41A, indicating the DNA interface of Cren7 was generally retained in F41A. However, F41A was unable to mediate DNA bridging, probably due to the impairment in forming oligomers/polymers on DNA. Atomic force microscopic images of the F41A-DNA complexes also revealed that F41A nearly completely lost the ability to compact DNA into highly condensed structures. Our results not only reveal the critical role of F41 in protein folding of Cren7 but also provide new insights into the structure-function relationships of thermostable proteins.

Evaluation of Multi-Objective Optimization Algorithms for NMR Chemical Shift Assignment

An automated NMR chemical shift assignment algorithm was developed using multi-objective optimization techniques. The problem is modeled as a combinatorial optimization problem and its objective parameters are defined separately in different score functions. Some of the heuristic approaches of evolutionary optimization are employed in this problem model. Both, a conventional genetic algorithm and multi-objective methods, i.e., the non-dominated sorting genetic algorithms II and III (NSGA2 and NSGA3), are applied to the problem. The multi-objective approaches consider each objective parameter separately, whereas the genetic algorithm followed a conventional way, where all objectives are combined in one score function. Several improvement steps and repetitions on these algorithms are performed and their combinations are also created as a hyper-heuristic approach to the problem. Additionally, a hill-climbing algorithm is also applied after the evolutionary algorithm steps. The algorithms are tested on several different datasets with a set of 11 commonly used spectra. The test results showed that our algorithm could assign both sidechain and backbone atoms fully automatically without any manual interactions. Our approaches could provide around a 65% success rate and could assign some of the atoms that could not be assigned by other methods.

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.

Switching an anti-IgG binding site between archaeal extremophilic proteins results in Affitins with enhanced pH stability

As a useful reagent for biotechnological applications, a scaffold protein needs to be as stable as possible to ensure longer lifetimes. We have developed archaeal extremophilic proteins from the “7 kDa DNA-binding” family as scaffolds to derive affinity proteins (Affitins). In this study, we evaluated a rational structure/sequence-guided approach to stabilize an Affitin derived from Sac7d by transferring its human IgG binding site onto the framework of the more thermally stable Sso7d homolog. The chimera obtained was functional, well expressed in Escherichia coli, but less thermally stable than the original Affitin (T(m) = 74.2 °C vs. T(m) = 80.4 °C). Two single mutations described as thermally stabilizing wild type Sso7d were introduced into chimeras. Only the double mutation nearly restored thermal stability (T(m) = 76.9 °C). Interestingly, the chimera and its double mutant were stable from pH 0 up to at least pH 13. Our results show that it is possible to increase further the stability of Affitins toward alkaline conditions (+2 pH units) while conserving their advantageous properties. As Affitins are based on a growing family of homologs from archaeal extremophiles, we conclude that this approach offers new potential for their improvement, which will be useful in demanding biotechnological applications.

Insights into thermoadaptation and the evolution of mesophily from the bacterial phylum Thermotogae

Thermophiles are extremophiles that grow optimally at temperatures >45 °C. To survive and maintain function of their biological molecules, they have a suite of characteristics not found in organisms that grow at moderate temperature (mesophiles). At the cellular level, thermophiles have mechanisms for maintaining their membranes, nucleic acids, and other cellular structures. At the protein level, each of their proteins remains stable and retains activity at temperatures that would denature their mesophilic homologs. Conversely, cellular structures and proteins from thermophiles may not function optimally at moderate temperatures. These differences between thermophiles and mesophiles presumably present a barrier for evolutionary transitioning between the 2 lifestyles. Therefore, studying closely related thermophiles and mesophiles can help us determine how such lifestyle transitions may happen. The bacterial phylum Thermotogae contains hyperthermophiles, thermophiles, mesophiles, and organisms with temperature ranges wide enough to span both thermophilic and mesophilic temperatures. Genomic, proteomic, and physiological differences noted between other bacterial thermophiles and mesophiles are evident within the Thermotogae. We argue that the Thermotogae is an ideal group of organisms for understanding of the response to fluctuating temperature and of long-term evolutionary adaptation to a different growth temperature range.

Affinity transfer to the archaeal extremophilic Sac7d protein by insertion of a CDR

Artificially transforming a scaffold protein into binders often consists of introducing diversity into its natural binding region by directed mutagenesis. We have previously developed the archaeal extremophilic Sac7d protein as a scaffold to derive affinity reagents (Affitins) by randomization of only a flat surface, or a flat surface and two short loops with natural lengths. Short loops are believed to contribute to stability of extremophilic proteins, and loop extension has been reported detrimental for the thermal and chemical stabilities of mesophilic proteins. In this work, we wanted to evaluate the possibility of designing target-binding proteins based on Sac7d by using a complementary determining region (CDR). To this aim, we inserted into three different loops a 10 residues CDR from the cAb-Lys3 anti-lysozyme camel antibody. The chimeras obtained were as stable as wild-type (WT) Sac7d at extreme pH and their structural integrity was supported. Chimeras were thermally stable, but with T(m)s from 60.9 to 66.3°C (cf. 91°C for Sac7d) which shows that loop extension is detrimental for thermal stability of Sac7d. The loop 3 enabled anti-lysozyme activity. These results pave the way for the use of CDR(s) from antibodies and/or extended randomized loop(s) to increase the potential of binding of Affitins.

Tolerance of the archaeal Sac7d scaffold protein to alternative library designs: characterization of anti-immunoglobulin G Affitins

Engineered protein scaffolds have received considerable attention as alternatives to antibodies in both basic and applied research, as they can offer superior biophysical properties often associated with a simpler molecular organization. Sac7d has been demonstrated as an effective scaffold for molecular recognition. Here, we used the initial L1 'flat surface' library constructed by randomization of 14 residues, to identify ligands specific for human immunoglobulin G. To challenge the plasticity of the Sac7d protein scaffold, we designed the alternative L2 'flat surface & loops' library whereof only 10 residues are randomized. Representative binders (Affitins) of the two libraries exhibited affinities in the low nanomolar range and were able to recognize different epitopes within human immunoglobulin G. These Affitins were stable up to pH 12 while largely conserving other favorable properties of Sac7d protein, such as high expression yields in Escherichia coli, solubility, thermal stability up to 80.7°C, and acidic stability (pH 0). In agreement with our library designs, mutagenesis study revealed two distinct binding areas, one including loops. Together, our results indicate that the Sac7d scaffold tolerates alternative library designs, which further expands the diversity of Affitins and may provide a general way to create tailored affinity tools for demanding applications.

Experimental assessment of the importance of amino acid positions identified by an entropy-based correlation analysis of multiple-sequence alignments

The analysis of a multiple-sequence alignment (MSA) with correlation methods identifies pairs of residue positions whose occupation with amino acids changes in a concerted manner. It is plausible to assume that positions that are part of many such correlation pairs are important for protein function or stability. We have used the algorithm H2r to identify positions k in the MSAs of the enzymes anthranilate phosphoribosyl transferase (AnPRT) and indole-3-glycerol phosphate synthase (IGPS) that show a high conn(k) value, i.e., a large number of significant correlations in which k is involved. The importance of the identified residues was experimentally validated by performing mutagenesis studies with sAnPRT and sIGPS from the archaeon Sulfolobus solfataricus. For sAnPRT, five H2r mutant proteins were generated by replacing nonconserved residues with alanine or the prevalent residue of the MSA. As a control, five residues with conn(k) values of zero were chosen randomly and replaced with alanine. The catalytic activities and conformational stabilities of the H2r and control mutant proteins were analyzed by steady-state enzyme kinetics and thermal unfolding studies. Compared to wild-type sAnPRT, the catalytic efficiencies (k(cat)/K(M)) were largely unaltered. In contrast, the apparent thermal unfolding temperature (T(M)(app)) was lowered in most proteins. Remarkably, the strongest observed destabilization (ΔT(M)(app) = 14 °C) was caused by the V284A exchange, which pertains to the position with the highest correlation signal [conn(k) = 11]. For sIGPS, six H2r mutant and four control proteins with alanine exchanges were generated and characterized. The k(cat)/K(M) values of four H2r mutant proteins were reduced between 13- and 120-fold, and their T(M)(app) values were decreased by up to 5 °C. For the sIGPS control proteins, the observed activity and stability decreases were much less severe. Our findings demonstrate that positions with high conn(k) values have an increased probability of being important for enzyme function or stability.

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.

Stabilization of a metabolic enzyme by library selection in Thermus thermophilus

The anthranilate phosphoribosyl transferase from the hyperthermophilic archaeon Sulfolobus solfataricus (sAnPRT, encoded by strpD), which catalyzes the third step in tryptophan biosynthesis, is a thermostable homodimer with low enzymatic activity at room temperature. We have combined two mutations leading to the monomerization and two mutations leading to the activation of sAnPRT. The resulting "activated monomer" sAnPRT-I36E-M47D+D83G-F149S, which is much more labile than wild-type sAnPRT, was stabilized by a combination of random mutagenesis and metabolic library selection using the extremely thermophilic bacterium Thermus thermophilus as host. This approach led to the identification of five mutations that individually increased the thermal stability of sAnPRT-I36E-M47D+D83G-F149S by 1 to 8 °C, and by 13 °C when combined. The beneficial exchanges were located in different parts of the protein structure, but none of them led to the "re-dimerization" of the enzyme. We observed a negative correlation between thermal stability and catalytic activity of the mutants; this suggests that conformational flexibility is required for catalysis by sAnPRT.

Contribution of hydrophobic interactions to protein stability

Our goal was to gain a better understanding of the contribution of hydrophobic interactions to protein stability. We measured the change in conformational stability, Δ(ΔG), for hydrophobic mutants of four proteins: villin headpiece subdomain (VHP) with 36 residues, a surface protein from Borrelia burgdorferi (VlsE) with 341 residues, and two proteins previously studied in our laboratory, ribonucleases Sa and T1. We compared our results with those of previous studies and reached the following conclusions: (1) Hydrophobic interactions contribute less to the stability of a small protein, VHP (0.6±0.3 kcal/mol per -CH(2)- group), than to the stability of a large protein, VlsE (1.6±0.3 kcal/mol per -CH(2)- group). (2) Hydrophobic interactions make the major contribution to the stability of VHP (40 kcal/mol) and the major contributors are (in kilocalories per mole) Phe18 (3.9), Met13 (3.1), Phe7 (2.9), Phe11 (2.7), and Leu21 (2.7). (3) Based on the Δ(ΔG) values for 148 hydrophobic mutants in 13 proteins, burying a -CH(2)- group on folding contributes, on average, 1.1±0.5 kcal/mol to protein stability. (4) The experimental Δ(ΔG) values for aliphatic side chains (Ala, Val, Ile, and Leu) are in good agreement with their ΔG(tr) values from water to cyclohexane. (5) For 22 proteins with 36 to 534 residues, hydrophobic interactions contribute 60±4% and hydrogen bonds contribute 40±4% to protein stability. (6) Conformational entropy contributes about 2.4 kcal/mol per residue to protein instability. The globular conformation of proteins is stabilized predominantly by hydrophobic interactions.

Molecular simulations on the thermal stabilization of DNA by hyperthermophilic chromatin protein Sac7d, and associated conformational transitions

Sac7d belongs to a family of chromosomal proteins, which are crucial for thermal stabilization of DNA at higher growth temperatures. It is capable of binding DNA nonspecifically, and is responsible for the increase in the melting temperature of DNA in the bound form up to 85 °C. Molecular dynamics (MD) simulations were performed at different temperatures on two protein-DNA complexes of Sac7d. Various structural and energetic parameters were calculated to examine the DNA stability and to investigate the conformational changes in DNA and the protein-DNA interactions. Room temperature simulations indicated very good agreement with the experimental structures. The protein structure is nearly unchanged at both 300 and 360 K, and only up to five base pairs of the DNA are stabilized by Sac7d at 360 K. However, the MD simulations on DNA alone systems show that they lose their helical structures at 360 K further supporting the role of Sac7d in stabilizing the oligomers. At higher temperatures (420 and 480 K), DNA undergoes denaturation in the presence and the absence of the protein. The DNA molecules were found to undergo B- to A-form transitions consistent with experimental studies, and the extent of these transitions are examined in detail. The extent of sampling B- and A-form regions was found to show temperature and sequence dependence. Multiple MD simulations yielded similar results validating the proposed model. Interaction energy calculations corresponding to protein-DNA binding indicates major contribution due to DNA backbone, explaining the nonspecific interactions of Sac7d.

Slow unfolding of monomeric proteins from hyperthermophiles with reversible unfolding

Based on the differences in their optimal growth temperatures microorganisms can be classified into psychrophiles, mesophiles, thermophiles, and hyperthermophiles. Proteins from hyperthermophiles generally exhibit greater stability than those from other organisms. In this review, we collect data about the stability and folding of monomeric proteins from hyperthermophilies with reversible unfolding, from the equilibrium and kinetic aspects. The results indicate that slow unfolding is a general strategy by which proteins from hyperthermophiles adapt to higher temperatures. Hydrophobic interaction is one of the factors in the molecular mechanism of the slow unfolding of proteins from hyperthermophiles.

Ligand-binding interactions and stability

The reversible interaction or binding of ligands to biological macromolecules is fundamental to nearly every aspect of biochemistry and cell biology. Binding events typically do not occur in isolation in biochemistry, and are almost always coupled or linked to other reactions such as protonation changes, other ligand-binding interactions, structural transitions, and folding. It is rarely sufficient to simply state that something binds. An understanding of binding requires a measure of affinity, stoichiometry, and the contributions of linked reactions. Emphasis is placed here on defining binding and the influence of linkage on binding and stability using both spectroscopic and calorimetric data.

Carboxyl pK(a) values, ion pairs, hydrogen bonding, and the pH-dependence of folding the hyperthermophile proteins Sac7d and Sso7d

Sac7d and Sso7d are homologous, hyperthermophile proteins with a high density of charged surface residues and potential ion pairs. To determine the relative importance of specific amino acid side-chains in defining the stability and function of these Archaeal chromatin proteins, pK(a) values were measured for the acidic residues in both proteins using (13)C NMR chemical shifts. The stability of Sso7d enabled titrations to pH 1 under low-salt conditions. Two aspartate residues in Sso7d (D16 and D35) and a single glutamate residue (G54) showed significantly perturbed pK(a) values in low salt, indicating that the observed pH-dependence of stability was primarily due to these three residues. The pH-dependence of backbone amide NMR resonances demonstrated that perturbation of all three pK(a) values was primarily the result of side-chain to backbone amide hydrogen bonds. Few of the significantly perturbed acidic pK(a) values in Sac7d and Sso7d could be attributed to primarily ion pair or electrostatic interactions. A smaller perturbation of E48 (E47 in Sac7d) was ascribed to an ion pair interaction that may be important in defining the DNA binding surface. The small number (three) of significantly altered pK(a) values was in good agreement with a linkage analysis of the temperature, pH, and salt-dependence of folding. The linkage of the ionization of two or more side-chains to protein folding led to apparent cooperativity in the pH-dependence of folding, although each group titrated independently with a Hill coefficient near unity. These results demonstrate that the acid pH-dependence of protein stability in these hyperthermophile proteins is due to independent titration of acidic residues with pK(a) values perturbed primarily by hydrogen bonding of the side-chain to the backbone. This work demonstrates the need for caution in using structural data alone to argue the importance of ion pairs in stabilizing hyperthermophile proteins.

Residue cluster additivity of thermodynamic stability in the hydrophobic core of mesophile vs. hyperthermophile rubredoxins

The branched sidechain residues 24 and 33 in the hydrophobic core of rubredoxin differ between the Clostridium pasteurianum (Cp) and Pyrococcus furiosus (Pf) sequences. Their X-ray structures indicate that these two sidechains are in van der Waals contact with each other, while neither appears to significantly interact with the other nonconserved residues. The simultaneous interchange of residues 24 and 33 between the Cp and Pf rubredoxin sequences yield a complementary pair of hybrid proteins for which the sum of their thermodynamic stabilities equals that of the parental rubredoxins. The 1.2 kcal/mol change arising from this two residues interchange accounts for 21% of the differential thermodynamic stability between the mesophile and hyperthermophile proteins. The additional interchange of the sole nonconserved aromatic residue in the hydrophobic core yields a 0.78 kcal/mol deviation from thermodynamic additivity.

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.

Contribution of the multi-turn segment in the reversible thermal stability of hyperthermophile rubredoxin: NMR thermal chemical exchange analysis of sequence hybrids

Pyrococcus furiosus (Pf) rubredoxin is the most thermostable protein characterized to date. Reflecting the complications arising from irreversible denaturation of this protein, predictions of which structural regions confer differential thermal stability have utilized kinetic stability measurements, hydrogen exchange protection factors, long range hydrogen bond NMR spin couplings, and molecular dynamics simulations, and have primarily implicated the three-stranded beta-sheet and the adjacent metal binding site. Herein, NMR chemical exchange experiments demonstrate reversible two-state unfolding at the thermal transition temperature (T(m)) for hybrids of Pf and the mesophile Clostridium pasteurianum (Cp) rubredoxins which interchange residues 14-33, the so-called multi-turn segment. This complementary pair of hybrid rubredoxins exhibits largely additive incremental thermal stabilizations vs. the parental proteins. Both stabilization free energy measurements as well as incremental T(m) values indicate that a minimum of 37% of the total differential thermal stability resides in this multi-turn segment. Such a proportionality between DeltaDeltaG and incremental T(m) values is predicted for hybrid pairs exhibiting thermodynamic additivity in which the differential stability is predominantly enthalpic.

The HPr proteins from the thermophile Bacillus stearothermophilus can form domain-swapped dimers

The study of proteins from extremophilic organisms continues to generate interest in the field of protein folding because paradigms explaining the enhanced stability of these proteins still elude us and such studies have the potential to further our knowledge of the forces stabilizing proteins. We have undertaken such a study with our model protein HPr from a mesophile, Bacillus subtilis, and a thermophile, Bacillus stearothermophilus. We report here the high-resolution structures of the wild-type HPr protein from the thermophile and a variant, F29W. The variant proved to crystallize in two forms: a monomeric form with a structure very similar to the wild-type protein as well as a domain-swapped dimer. Interestingly, the structure of the domain-swapped dimer for HPr is very different from that observed for a homologous protein, Crh, from B.subtilis. The existence of a domain-swapped dimer has implications for amyloid formation and is consistent with recent results showing that the HPr proteins can form amyloid fibrils. We also characterized the conformational stability of the thermophilic HPr proteins using thermal and solvent denaturation methods and have used the high-resolution structures in an attempt to explain the differences in stability between the different HPr proteins. Finally, we present a detailed analysis of the solution properties of the HPr proteins using a variety of biochemical and biophysical methods.

Thermodynamics of DNA binding and distortion by the hyperthermophile chromatin protein Sac7d

Sac7d is a hyperthermophile chromatin protein which binds non-specifically to the minor groove of duplex DNA and induces a sharp kink of 66 degrees with intercalation of valine and methionine side-chains. We have utilized the thermal stability of Sac7d and the lack of sequence specificity to define the thermodynamics of DNA binding over a wide temperature range. The binding affinity for poly(dGdC) was moderate at 25 degrees C (Ka = 3.5(+/-1.6) x 10(6) M(-1)) and increased by nearly an order of magnitude from 10 degrees C to 80 degrees C. The enthalpy of binding was unfavorable at 25 degrees C, and decreased linearly from 5 degrees C to 60 degrees C. A positive binding heat at 25 degrees C is attributed in part to the energy of distorting DNA, and ensures that the temperature of maximal binding affinity (75.1+/-5.6 degrees C) is near the growth temperature of Sulfolobus acidocaldarius. Truncation of the two intercalating residues to alanine led to a decreased ability to bend and unwind DNA at 25 degrees C with a small decrease in binding affinity. The energy gained from intercalation is slightly greater than the free energy penalty of bending duplex DNA. Surprisingly, reduced distortion from the double alanine substitution did not lead to a significant decrease in the heat of binding at 25 degrees C. In addition, an anomalous positive DeltaCp of binding was observed for the double alanine mutant protein which could not be explained by the change in polar and apolar accessible surface areas. Both the larger than expected binding enthalpy and the positive heat capacity can be explained by a temperature dependent structural transition in the protein-DNA complex with a Tm of 15-20 degrees C and a DeltaH of 15 kcal/mol. Data are discussed which indicate that the endothermic transition in the complex is consistent with DNA distortion.