Conformational stability and activity of ribonuclease T1 and mutants

Structure and dynamics of the fish eye lens protein, γM7-crystallin

The vertebrate eye lens contains high concentrations of crystallins. The dense lenses of fish are particularly abundant in a class called γM-crystallin whose members are characterized by an unusually high methionine content and partial loss of the four tryptophan residues conserved in all γ-crystallins from mammals which are proposed to contribute to protection from UV-damage. Here, we present the structure and dynamics of γM7-crystallin from zebrafish (Danio rerio). The solution structure shares the typical two-domain, four-Greek-key motif arrangement of other γ-crystallins, with the major difference noted in the final loop of the N-terminal domain, spanning residues 65-72. This is likely due to the absence of the conserved tryptophans. Many of the methionine residues are exposed on the surface but are mostly well-ordered and frequently have contacts with aromatic side chains. This may contribute to the specialized surface properties of these proteins that exist under high molecular crowding in the fish lens. NMR relaxation data show increased backbone conformational motions in the loop regions of γM7 compared to those of mouse γS-crystallin and show that fast internal motion of the interdomain linker in γ-crystallins correlates with linker length. Unfolding studies monitored by tryptophan fluorescence confirm results from mutant mouse γS-crystallin and show that unfolding of a βγ-crystallin domain likely starts from unfolding of the variable loop containing the more fluorescently quenched tryptophan residue, resulting in a native-like unfolding intermediate.

Monomeric banana lectin at acidic pH overrules conformational stability of its native dimeric form

Banana lectin (BL) is a homodimeric protein categorized among jacalin-related family of lectins. The effect of acidic pH was examined on conformational stability of BL by using circular dichroism, intrinsic fluorescence, 1-anilino-8-napthalene sulfonate (ANS) binding, size exclusion chromatography (SEC) and dynamic light scattering (DLS). During acid denaturation of BL, the monomerization of native dimeric protein was found at pH 2.0. The elution profile from SEC showed two different peaks (59.65 ml & 87.98 ml) at pH 2.0 while single peak (61.45 ml) at pH 7.4. The hydrodynamic radii (R h) of native BL was 2.9 nm while at pH 2.0 two species were found with R h of 1.7 and 3.7 nm. Furthermore at, pH 2.0 the secondary structures of BL remained unaltered while tertiary structure was significantly disrupted with the exposure of hydrophobic clusters confirming the existence of molten globule like state. The unfolding of BL with different subunit status was further evaluated by urea and temperature mediated denaturation to check their stability. As inferred from high Cm and ΔG values, the monomeric form of BL offers more resistance towards chemical denaturation than the native dimeric form. Besides, dimeric BL exhibited a Tm of 77°C while no loss in secondary structures was observed in monomers even up to 95°C. To the best of our knowledge, this is the first report on monomeric subunit of lectins showing more stability against denaturants than its native dimeric state.

Tryptophan conformations associated with partial unfolding in ribonuclease T1

The origin of the biexponential fluorescence decay of Trp in ribonuclease T1 under mildly destabilizing conditions, such as increased pH or temperature, or the presence of detergent, is still not understood. We have performed two extended replica-exchange molecular dynamics simulations to obtain a detailed representation of the native state at two protonation states corresponding to a high and low pH. At high pH, the appearance of partially unfolded states is evident. We found that this pH-induced destabilization originates from increased global repulsion as well as reduced local favorable electrostatic interactions and reduced H-bonding strength of His(27), His(40), and His(92). At high pH, alternative tryptophan rotamers appear and are linked to a distorted environment of the tryptophan, which also acts as a separate source of ground-state heterogeneity. The total population of these alternative conformations agrees well with the amplitude of the experimentally observed secondary fluorescence lifetime.

SIMPLE estimate of the free energy change due to aliphatic mutations: Superior predictions based on first principles

The bioinformatics revolution of the last decade has been instrumental in the development of empirical potentials to quantitatively estimate protein interactions for modeling and design. Although computationally efficient, these potentials hide most of the relevant thermodynamics in 5-to-40 parameters that are fitted against a large experimental database. Here, we revisit this longstanding problem and show that a careful consideration of the change in hydrophobicity, electrostatics, and configurational entropy between the folded and unfolded state of aliphatic point mutations predicts 20-30% less false positives and yields more accurate predictions than any published empirical energy function. This significant improvement is achieved with essentially no free parameters, validating past theoretical and experimental efforts to understand the thermodynamics of protein folding. Our first principle analysis strongly suggests that both the solute-solute van der Waals interactions in the folded state and the electrostatics free energy change of exposed aliphatic mutations are almost completely compensated by similar interactions operating in the unfolded ensemble. Not surprisingly, the problem of properly accounting for the solvent contribution to the free energy of polar and charged group mutations, as well as of mutations that disrupt the protein backbone remains open.

Activating an Enzyme by an Engineered Coiled Coil Switch

We designed a de novo protein based on a circular permutant of RNaseT1, in which the enzymatic activity can be manipulated by engineered peptide binding. The circular permutant of RNaseT1 was obtained by tethering the original C- and N-termini with a GPAG linker and cleaving the molecule between Glu82 and Asn83. This mutant lacked enzymatic activity, due to the destabilization of entire protein structure. We previously reported the construction of ABC-type heterotrimeric coiled coil peptides, in which the A- and B-type peptides cannot form the folded trimeric structure without the C-type peptide. The introduction of the A- and B-type coiled coil peptides to the C- and N-termini of the circular permutant of RNaseT1, respectively, and the subsequent addition of the C-type coiled coil peptide enabled the RNaseT1 domain to refold properly, thus, restoring the enzymatic activity. The formation of the trimeric coiled coil structure should bring the cleaved sites of RNaseT1 close enough to refold the RNaseT1 domain spontaneously.

Asp79 makes a large, unfavorable contribution to the stability of RNase Sa

The two most buried carboxyl groups in ribonuclease Sa (RNase Sa) are Asp33 (99% buried; pK 2.4) and Asp79 (85% buried; pK 7.4). Above these pK values, the stability of the D33A variant is 6kcal/mol less than wild-type RNase Sa, and the stability of the D79A variant is 3.3kcal/mol greater than wild-type RNase Sa. The key structural difference between the carboxyl groups is that Asp33 forms three intramolecular hydrogen bonds, and Asp79 forms no intramolecular hydrogen bond. Here, we focus on Asp79 and describe studies of 11 Asp79 variants. Most of the variants were at least 2kcal/mol more stable than wild-type RNase Sa, and the most interesting was D79F. At pH 3, below the pK of Asp79, RNase Sa is 0.3kcal/mol more stable than the D79F variant. At pH 8.5, above the pK of Asp79, RNase Sa is 3.7kcal/mol less stable than the D79F variant. The unfavorable contribution of Asp79 to the stability appears to result from the Born self-energy of burying the charge and, more importantly, from unfavorable charge-charge interactions. To counteract the effect of the negative charge on Asp79, we prepared the Q94K variant and the crystal structure showed that the amino group of the Lys formed a hydrogen-bonded ion pair (distance, 2.71A; angle, 100 degrees ) with the carboxyl group of Asp79. The stability of the Q94K variant was about the same as the wild-type at pH 3, where Asp79 is uncharged, but 1kcal/mol greater than that of wild-type RNase Sa at pH 8.5, where Asp79 is charged. Differences in hydrophobicity, steric strain, Born self-energy, and electrostatic interactions all appear to contribute to the range of stabilities observed in the variants. When it is possible, replacing buried, non-hydrogen bonded, ionizable side-chains with non-polar side-chains is an excellent means of increasing protein stability.

Effect of Salt on the Urea-Unfolded Form of Barstar Probed bymValue Measurements†

To probe for residual structure present in the urea-unfolded form of the small protein barstar, to determine how salt might modulate such structure, and to determine how such structure might affect the stability of the protein, mutant variants that display m values different from that of the wild-type protein have been studied. The mutant proteins were obtained by site-directed mutagenesis at residue positions located on the surface of the folded protein. The m value, which represents the preferential free energy of interaction of urea with the unfolded form in comparison to that with the folded state, was determined from equilibrium urea-induced unfolding curves. Mutant proteins for which the m values were significantly greater than (m(+) mutant forms), significantly smaller than (m(-) mutant forms), or similar to (m(0) mutant forms) the m value determined for the wild-type protein were studied. The unfolded forms of the m(0), m(+) and m(-) mutant proteins represent different components within the unfolded form ensemble, which differ from each other in their solvent-exposed surface areas. Hence, the m value has been used as a measure of residual structure in the unfolded form. To further understand the nature of structures present in the unfolded form ensemble, the effects of the salt KCl on the stabilities of the wild-type and the mutant proteins, as well as on the structures present in the unfolded form ensemble, were also studied. It was found that the m values of the m(0), m(+) and m(-) mutant proteins all converge to the wild-type m value in the presence of KCl. This result indicates that the salt modulates residual structure in the unfolded form by screening electrostatic interactions that maintain compact and expanded components in the unfolded protein ensemble. The use of free energy cycles has allowed the effect of salt on the structure and free energy of the unfolded protein to be related to the stability of the protein.

Determination of the NMR structure of Gln25-ribonuclease T1

Ribonuclease (RNase) T1 is a guanyloribonuclease, having two isozymes in nature, Gln25- and Lys25-RNase T1. Between these two isozymes, there is no difference in catalytic activity and three-dimensional structure; however, Lys25-RNase T1 is slightly more stable than Gln25-RNase T1. Recently, it has been suggested that the existence of a salt bridge between Lys25 and Asp29/Glu31 in Lys25-RNase T1 contributes to the stability. To elucidate the effects of the replacement of Lys25 with a Gln on the conformation and microenvironments of RNase T1 in detail, the three-dimensional solution structure of Gln25-RNase T1 was determined by simulated-annealing calculations. As a result, the topology of the overall folding was shown to be very similar to that of the Lys25-isozyme except for some differences. In particular, there were two differences in the property of torsion angles of the two disulfide bonds and the conformations of the residues 11-13, 63-66, and 92-93. With regard to the residues 11-13, the lack of the above-mentioned salt bridge in Gln25-RNase T1 was thought to induce the conformational difference of this segment as compared with the Lys25-isozyme. Furthermore, it was proposed that the perturbation of this segment might transfer to the residues 92-93 via the two disulfide bonds.

Structurally homologous all beta-barrel proteins adopt different mechanisms of folding

Acidic fibroblast growth factors from human (hFGF-1) and newt (nFGF-1) (Notopthalamus viridescens) are 16-kDa, all beta-sheet proteins with nearly identical three-dimensional structures. Guanidine hydrochloride (GdnHCl)-induced unfolding of hFGF-1 and nFGF-1 monitored by fluorescence and far-UV circular dichroism (CD) shows that the FGF-1 isoforms differ significantly in their thermodynamic stabilities. GdnHCl-induced unfolding of nFGF-1 follows a two-state (Native state to Denatured state(s)) mechanism without detectable intermediate(s). By contrast, unfolding of hFGF-1 monitored by fluorescence, far-UV circular dichroism, size-exclusion chromatography, and NMR spectroscopy shows that the unfolding process is noncooperative and proceeds with the accumulation of stable intermediate(s) at 0.96 M GdnHCl. The intermediate (in hFGF-1) populated maximally at 0.96 M GdnHCl has molten globule-like properties and shows strong binding affinity to the hydrophobic dye, 1-Anilino-8-naphthalene sulfonate (ANS). Refolding kinetics of hFGF-1 and nFGF-1 monitored by stopped-flow fluorescence reveal that hFGF-1 and nFGF-1 adopts different folding mechanisms. The observed differences in the folding/unfolding mechanisms of nFGF-1 and hFGF-1 are proposed to be either due to differential stabilizing effects of the charged denaturant (Gdn(+) Cl(-)) on the intermediate state(s) and/or due to differences in the structural interactions stabilizing the native conformation(s) of the FGF-1 isoforms.

Testing hypotheses about determinants of protein structure with high-precision, high-throughput stability measurements and statistical modeling

Statistical modeling provides the mathematics to use data from large numbers of mutant proteins to generate information about hypotheses concerning protein structure not easily obtained from anecdotal studies on small numbers of mutants. Here we use the unfolding free energies of 303 unique eglin c mutant proteins obtained from high-precision, high-throughput chemical denaturation measurements to assess models concerning helix stability. A model with helix propensity as the sole determinant of stability accounts for 83% of the mutant-to-mutant variation in stability for 99% of the mutant proteins (three outliers). When position effects and side chain-side chain interactions are added to the model, the fraction of variation explained increases to 92%. The propensity parameters in this model are identical to helix propensity values derived from other approaches. Measurement error accounts for another 1% of the mutant-to-mutant variation in stability. While the data support terms for several of the expected stabilizing/destabilizing effects, it does not support terms for several others, including i, i + 3 effects in the center of the helix and helix-dipole effects. In addition, the model does better with terms for several stabilizing/destabilizing effects for which we cannot identify the physical basis. The precision of our unfolding stability measurements (+/-0.087 kcal/mol) allows us to conclude that the 7% of variation in stabilities of the mutant proteins not accounted for by the model or by measurement variation is both real and large with respect to the nonpropensity terms in the model. The analysis also shows that the common practice of using C(m)m(av) instead of C(m)m(mut) to calculate DeltaG(HOH,N-D) values for each mutant protein results in a loss of information. We see no correlation between the residuals derived from the full model and m(mut) - m(wt), and hence it is unlikely our m(mut) values reflect mutant-to-mutant differences in the denatured state.

Thermodynamic analysis of the contributions of the copper ion and the disulfide bridge to azurin stability: synergism among multiple depletions

The stabilizing potential of the copper ion and the disulfide bridge in azurin has been explored with the aim of inspecting the ways in which these two factors influence one another. Specifically, whether copper and disulfide contributions to protein stability are additive has been examined. To this aim, the thermal unfolding of a copper-depleted mutant lacking the disulfide bridge between Cys3 and Cys26 (apo C3A/C26A azurin) was studied by differential scanning calorimetry. A comparison of the unfolding parameters of holo and apo C3A/C26A azurin with the apo C3A/C26A protein has shown that the effects of simultaneous copper and disulfide depletion are additive only at two temperatures: T=15 degrees C and T=67 degrees C. Within this range the presence of the copper ion and the disulfide bridge has a positive synergistic effect on azurin stability. These findings might have implications for the rational use of the stabilizing potential of copper and disulfides in copper protein engineering.

Instability, stabilization, and formulation of liquid protein pharmaceuticals

One of the most challenging tasks in the development of protein pharmaceuticals is to deal with physical and chemical instabilities of proteins. Protein instability is one of the major reasons why protein pharmaceuticals are administered traditionally through injection rather than taken orally like most small chemical drugs. Protein pharmaceuticals usually have to be stored under cold conditions or freeze-dried to achieve an acceptable shelf life. To understand and maximize the stability of protein pharmaceuticals or any other usable proteins such as catalytic enzymes, many studies have been conducted, especially in the past two decades. These studies have covered many areas such as protein folding and unfolding/denaturation, mechanisms of chemical and physical instabilities of proteins, and various means of stabilizing proteins in aqueous or solid state and under various processing conditions such as freeze-thawing and drying. This article reviews these investigations and achievements in recent years and discusses the basic behavior of proteins, their instabilities, and stabilization in aqueous state in relation to the development of liquid protein pharmaceuticals.

A self-consistent, microenvironment modulated screened coulomb potential approximation to calculate pH-dependent electrostatic effects in proteins

An improved approach is presented for calculating pH-dependent electrostatic effects in proteins using sigmoidally screened Coulomb potentials (SCP). It is hypothesized that a key determinant of seemingly aberrant behavior in pKa shifts is due to the properties of the unique microenvironment around each residue. To help demonstrate this proposal, an approach is developed to characterize the microenvironments using the local hydrophobicity/hydrophilicity around each residue of the protein. The quantitative characterization of the microenvironments shows that the protein is a complex mosaic of differing dielectric regions that provides a physical basis for modifying the dielectric screening functions: in more hydrophobic microenvironments the screening decreases whereas the converse applies to more hydrophilic regions. The approach was applied to seven proteins providing more than 100 measured pKa values and yielded a root mean square deviation of 0.5 between calculated and experimental values. The incorporation of the local hydrophobicity characteristics into the algorithm allowed the resolution of some of the more intractable problems in the calculation of pKa. Thus, the divergent shifts of the pKa of Glu-35 and Asp-66 in hen egg white lysozyme, which are both about 90% buried, was correctly predicted. Mechanistically, the divergence occurs because Glu-35 is in a hydrophobic microenvironment, while Asp-66 is in a hydrophilic microenvironment. Furthermore, because the calculation of the microenvironmental effects takes very little CPU time, the computational speed of the SCP formulation is conserved. Finally, results from different crystal structures of a given protein were compared, and it is shown that the reliability of the calculated pKa values is sufficient to allow identification of conformations that may be more relevant for the solution structure.

Monitoring the sizes of denatured ensembles of staphylococcal nuclease proteins: implications regarding m values, intermediates, and thermodynamics

Fluorescence and size-exclusion chromatography (SEC) are used to monitor urea denaturation of wild-type staphylococcal nuclease (SN) as well as the m+ and m- mutants A69T and V66W, respectively. It is found that the SEC partition coefficient, 1/Kd, is directly proportional to the Stokes radii of proteins. From the Stokes radii, the denatured ensembles of the three proteins are found to be highly compact in the limit of low urea concentration and expand significantly with increasing urea concentration. The m values from fluorescence-detected denaturation of the SN proteins are generally considered to reflect the relative sizes of denatured ensembles. However, the rank order of m values of the SN proteins studied do not correspond to the rank order of denatured ensemble sizes detected by 1/Kd, suggesting that m values reflect more than just surface area increases on denaturation. SEC provides two complementary ways to demonstrate the existence of intermediates in urea denaturation and illustrates that V66W undergoes a three-state transition. Fluorescence-detected urea denaturations of A69T and wt SN do not correspond with 1/Kd-detected denaturation profiles, a result that would ordinarily mean that the transitions are non-two-state. However, this interpretation fails to recognize the rapidly changing size and thermodynamic character of the denatured ensembles of these proteins both within and outside of the transition zone. The implications of the changing sizes and thermodynamic character of the denatured ensembles for SN proteins are manifold, requiring a reconsideration of the thermodynamics of proteins whose denatured ensembles behave as those of SN proteins.

A core folding model for catalysis by the hammerhead ribozyme accounts for its extraordinary sensitivity to abasic mutations

Introducing abasic nucleotides at each of 13 positions in the conserved core of the hammerhead ribozyme causes a large decrease in the extent of catalysis [Peracchi, A., et al. (1996) Proc. Natl. Acad. Sci. U.S.A. 93, 11522]. This extreme sensitivity to structural defects is in contrast to the behavior of protein enzymes and larger ribozymes. Several additional differences in the behavior of the hammerhead relative to that of protein enzymes and larger ribozymes are described herein. The deleterious effects of the abasic mutations are not relieved by lowering the temperature, by increasing the concentration of monovalent or divalent metal ions, or by adding polyamines, in contrast to effects observed with protein enzymes and large RNA enzymes. In addition, the abasic mutations do not significantly weaken substrate binding. These results and previous observations are all accounted for by a "core folding" model in which the stable ground state structure of the hammerhead ribozyme complexed with the substrate is a partially folded state that must undergo an additional folding event to achieve its catalytic conformation. We propose that the peculiar behavior of the hammerhead arises because the limited structural interconnections in a small RNA enzyme do not allow the ground state to stably adopt the catalytic conformation; within the globally folded catalytic conformation, limited structural interconnections may further impair catalysis by hampering the precise alignment of active site functional groups. This behavior represents a basic manifestation of the well-recognized interconnection between folding and catalysis.

Thermodynamic and kinetic characterization of calf lens gammaF-crystallin

Gamma-crystallin is reported to be conformationally stable because of its internal structural symmetry, and gammaF (gammaIVa) is the most stable among the various gamma-crystallin gene products. However, there is no detailed report on its thermodynamic and kinetic stability. In the present study, detailed unfolding of gammaF-crystallin was investigated by equilibrium and kinetics methods with fluorescence and far-UV CD spectroscopic measurements. The GdnHCl-induced unfolding curves probed by Trp emission maximum and intensity showed a sharp single-step transition. Upon widening the unfolding transition with the use of urea in 1.5 M GdnHCl, a more proper fit for thermodynamic analysis was obtained. GammaF-Crystallin underwent a straightforward two-state process (N <==> U) without showing any measurable amount of intermediate. The conformational stability, as measured by deltaG(D)H2O (approximately 9 kcal/mol), indicates that gammaF-crystallin is a very stable protein. The high activation energy deltaG++H2O (approximately 24 kcal/mol), calculated from unfolding kinetics monitored by far-UV CD at 218 nm, also indicates that the native and unfolded states are separated by a high activation energy barrier.

Conformational stability and thermodynamics of folding of ribonucleases Sa, Sa2 and Sa3

Ribonucleases Sa, Sa2, and Sa3 are three small, extracellular enzymes produced by different strains of Streptomyces aureofaciens with amino acid sequences that are 50% identical. We have studied the unfolding of these enzymes by heat and urea to determine the conformational stability and its dependence on temperature, pH, NaCl, and the disulfide bond. All three of the Sa ribonucleases unfold reversibly by a two-state mechanism with melting temperatures, Tm, at pH 7 of 48.4 degrees C (Sa), 41.1 degrees C (Sa2), and 47.2 degrees C (Sa3). The Tm values are increased in the presence of 0.5 M NaCl by 4.0 deg. C (Sa), 0.1 deg. C (Sa2), and 7.2 deg. C (Sa3). The Tm values are decreased by 20.0 deg. C (Sa), 31.5 deg. C (Sa2), and 27.0 deg. C (Sa3) when the single disulfide bond in the molecules is reduced. We compare these results with similar studies on two other members of the microbial ribonuclease family, RNase T1 and RNase Ba (barnase), and with a member of the mammalian ribonuclease family, RNase A. At pH 7 and 25 degrees C, the conformational stabilities of the ribonucleases are (kcal/mol): 2.9 (Sa2), 5.6 (Sa3), 6.1 (Sa), 6.6 (T1), 8.7 (Ba), and 9.2 (A). Our analysis of the stabilizing forces suggests that the hydrophobic effect contributes from 90 to 110 kcal/mol and that hydrogen bonding contributes from 70 to 105 kcal/mol to the stability of these ribonucleases. Thus, we think that the hydrophobic effect and hydrogen bonding make large but comparable contributions to the conformational stability of these proteins.

A decade of protein engineering on ribonuclease T1--atomic dissection of the enzyme-substrate interactions

During the last decade, protein engineering has been used to identify the residues that contribute to the ribonuclease-T1-catalyzed transesterification. His40, Glu58 and His92 accelerate the associative nucleophilic displacement at the phosphate atom by the entering 2'-oxygen downstream guanosines in a highly cooperative manner. Glu58, assisted by the protonated His40 imidazole, abstracts a proton from the 2'-oxygen, while His92 protonates the leaving group. Tyr38, Arg77 and Phe100 further stabilize the transition state of the reaction. A functionally independent subsite, including Asn36 and Asn98, contributes to chemical turnover by aligning the substrate relative to the catalytic side chains upon binding of the leaving group. An invariant structural motive, involving residues 42-46, renders ribonuclease T1 guanine specific through a series of intermolar hydrogen bonds. Tyr42 contributes significantly to guanine binding through a parallel face-to-face stacking interaction. Tyr45, often referred to as the lid of the guanine-binding site, does not contribute to the binding of the base.

Self-consistent field approach to protein structure and stability. I: pH dependence of electrostatic contribution

Starting from the simple case of an external field acting on noninteracting particles, a formulation of the self-consistent field theory for treating proteins and unfolded protein chains with multiple interacting titratable groups is given. Electrostatic interactions between the titratable groups are approximated by a Debye-Huckel expression. Amino acid residues are treated as polarizable bodies with a single dielectric constant. Dielectric properties of protein molecules are described in terms of local dielectric constants determined by the space distribution of residue volume density around each ionized residue. Calculations are based on average charges of titratable groups, distance of separation between them, on their pKa's, residue volumes and on the local dielectric constant. A set of different residue volumes is used to analyze the influence of the permanent dipole of polar parts of the residue on calculated titration curves, electrostatic contribution to the free energy of protein stability, and pK shifts. Calculations with zero volumes--which means that charged portions of protein molecules are viewed as part of the high dielectric medium--give good agreement with experimental data. The theory was tested against most accurate approaches currently available for the calculation of the pKa's of ionizable groups based upon finite difference solutions of the Poisson-Boltzmann equation (FDPB). For 70 theoretically calculated pKa's in a total of six proteins the accuracy of the approach presented here is assessed by comparison of computed pKa's with that measured. The overall root-mean-square error is 0.79, compared to the value 0.89 obtained by FDPB approach given in the paper of Antosiewicz et al. (J. Mol. Biol. 238:415-436, 1994). The test of Debye-Huckel approximation for the electrostatic pair interactions shows that it is in excellent agreement with experimental data as well as the calculations of the FDPB and Tanford-Kirkwood methods on the pK shifts of His64 in the active site of subtilisin over the whole range of ionic strengths. (Gilson and Honig, Proteins 3:32-52, 1988; Russell et al., J.Mol.Biol. 193:803-813, 1987). The theory was also analytically and numerically tested on a simple models where the exact statistical mechanical treatment is still simple (Yang et al., Proteins 15:252-265, 1993; Bashford and karplus, J. Phys. Chem. 95:9556-9561, 1991).

Thermodynamic features of the thermal unfolding of human serum albumin

The unfolding process of human serum albumin (HSA) was studied by thermal effect on the native fluorescence of the protein, thermal inactivation of the hydrolase activity of albumin and differential scanning calorimetry using the high sensitive calorimeter developed by Privalov. The denaturation process can be described by an approximation of the model of Eyring and Lumry: native [symbol: see text] unfolded reversible [symbol: see text] unfolded irreversible. It was found that the rate of irreversible step was very slow (at temperatures below 74 degrees C), allowing the resolution of the denaturation process as a reversible one on the basis of two states approximation. However, the presence of intramolecular cooperation in the thermal denaturation process at temperatures above 74 degrees C cannot be discarded, which might be favoring the aggregation of albumin molecules. The midpoint temperature of unfolding obtained by differential scanning calorimetry was of 63.1 degrees C +/- 0.4 at pH 7.4. This value was independent of the rate of scanning and it is in agreement with those obtained by techniques such as thermal effect on the protein fluorescence and on the hydrolase activity of albumin. The enthalpy of unfolding at pH 7.4 was 88.9 +/- 4 Kcal/mol. This value was low compared with those obtained for other proteins, suggesting the presence of a molten globule in the unfolding pathway of albumin. The neutral-basic conformational change (pH 7.4) of albumin did not modify the thermal stability and the enthalpy of denaturation of the protein. A pH below 4.3 (transition acid-neutral) the presence of a second peak in the thermogram of albumin with a TM of 46.2 degrees C +/- 0.9 would be suggesting a lost of cooperativity between the various domains of albumin in the unfolding.

The determinants of pKas in proteins

Although validation studies show that theoretical models for predicting the pKas of ionizable groups in proteins are increasingly accurate, a number of important questions remain: (1) What factors limit the accuracy of current models? (2) How can conformational flexibility of proteins best be accounted for? (3) Will use of solution structures in the calculations, rather than crystal structures, improve the accuracy of the computed pKas? and (4) Why does accurate prediction of protein pKas seem to require that a high dielectric constant be assigned to the protein interior? This paper addresses these and related issues. Among the conclusions are the following: (1) computed pKas averaged over NMR structure sets are more accurate than those based upon single crystal structures; (2) use of atomic parameters optimized to reproduce hydration energies of small molecules improves agreement with experiment when a low protein dielectric constant is assumed; (3) despite use of NMR structures and optimized atomic parameters, pKas computed with a protein dielectric constant of 20 are more accurate than those computed with a low protein dielectric constant; (4) the pKa shifts in ribonuclease A that result from phosphate binding are reproduced reasonably well by calculations; (5) the substantial pKa shifts observed in turkey ovomucoid third domain result largely from interactions among ionized groups; and (6) both experimental data and calculations indicate that proteins tend to lower the pKas of Asp side chains but have little overall effect upon the pKas of other ionizable groups.

DsbA-mediated disulfide bond formation and catalyzed prolyl isomerization in oxidative protein folding

The interrelationship between the acquisition of ordered structure, prolyl isomerization, and the formation of the disulfide bonds in assisted protein folding was investigated by using a variant of ribonuclease T1 (C2S/C10N-RNase T1) with a single disulfide bond and two cis-prolyl bonds as a model protein. The thiol-disulfide oxidoreductase DsbA served as the oxidant for forming the disulfide bond and prolyl isomerase A as the catalyst of prolyl isomerization. Both enzymes are from the periplasm of Escherichia coli. Reduced C2S/C10N-RNase T1 is unfolded in 0 M NaCl, but native-like folded in > or = 2 M NaCl. Oxidation of 5 microM C2S/C10N-RNase T1 by 8 microM DsbA (at pH 7.0, 25 degrees C) is very rapid with a t1/2 of about 10 s (the second-order rate constant is 7 x 10(3) s-1 M-1), irrespective of whether the reduced molecules are unfolded or folded. When they are folded, the product of oxidation is the native protein. When they are denatured, first the disulfide bond is formed in the unfolded protein chains and then the native structure is acquired. This slow reaction is limited in rate by prolyl isomerization and catalyzed by prolyl isomerase. The efficiency of this catalysis is strongly decreased by the presence of the disulfide bond. Apparently, the rank order of chain folding, prolyl isomerization, and disulfide bond formation can vary in the oxidative folding of ribonuclease T1. Such a degeneracy could generally be an advantage for protein folding both in vitro and in vivo.

Thermal stabilization of ribonuclease T1 by carboxymethylation at Glu-58 as revealed by 1H nuclear magnetic resonance spectroscopy

Ribonuclease T1 (RNase T1) carboxymethylated at the gamma-carboxyl group of Glu-58 with iodoacetic acid is known to be completely inactive while it retains an almost full substrate-binding ability. In order to further clarify the effects of the carboxymethylation, the thermal stabilities of intact and Glu-58-carboxymethylated (CM-) RNase T1 were compared by measuring 1H NMR spectra at various temperatures. The transition curves of unfolding were obtained by plotting, as a function of temperature, the peak areas for the alpha and delta protons of Asn-81 and Ile-90, respectively, which are well apart from each other in the three-dimensional structure of the enzyme. For each of intact and CM-RNase T1, the transition curve of the Asn-81 alpha proton was identical with that of the Ile-90 delta methyl protons, suggesting that the thermal unfolding occurred simultaneously in every part of the molecule of CM-RNase T1 as well as of intact RNase T1. The midpoint of unfolding was 52 degrees C for intact RNase T1, and was increased by 9 degrees C upon carboxymethylation at Glu-58. This marked stabilization by carboxymethylation is thought to be due to formation of a salt bridge between the introduced carboxymethyl group and the neighboring guanidium group of Arg-77.

Site-directed mutagenesis of cysteine residues of hepatitis B surface antigen. Analysis of two single mutants and the double mutant

The structure of hepatitis B surface antigen (HBsAg) is mainly maintained by an intricate disulfide network responsible for most of its structural and antigenic properties. Characterization of three cysteine-replacement mutants of HBsAg has been performed by both structural and immunological methods. Replacement of Cys121 or Cys124 with serine results in mutant proteins that show diminished binding titres to both monoclonal antibodies and to a polyclonal serum, indicating that a structural change has taken place. Circular dichroism analysis shows that the substitution of either of these two residues also diminishes the helical content of the protein. However, the double mutant, in which both cysteine residues have been simultaneously changed, reverts the properties of the single mutations, and shows similar behaviour to the wild-type protein. Both the single and double cysteine mutants are efficiently glycosylated and secreted from Chinese hamster ovary cells and, in all cases, the mutant proteins assemble into spherical particles of similar buoyant density to both the wild-type and serum derived HBsAg.

pH-dependent thermostabilization of Escherichia coli ribonuclease HI by histidine to alanine substitutions

Thermal stabilities of mutant ribonuclease HI proteins from Escherichia coli, in which each of five histidine residues was replaced with alanine, were examined at various pHs. Increases in the Tm values were observed at pH 3.0 for four of the mutant proteins, in which each of the four histidine residues exposed to the solvent was mutated, as compared to the Tm of the wild-type protein. The thermostabilization of three of the mutant proteins was dependent on pH, and only observed at low pH. The thermostabilizing effects of the His-->Ala substitutions were cumulative. The temperature of the midpoint of the transition in the thermal unfolding curves, Tm, of the most stable mutant enzyme, in which His 62, His 83, His 124, and His 127 were replaced by Ala, was 5.5 degrees C higher than that of the wild-type enzyme at pH 3.0. The stability of the wild-type protein decreased as the pH was lowered below pH 4, a condition favoring the protonation of carboxyl groups, probably due to unfavorable electrostatic interactions introduced by the increase in positive charges on the protein. Since imidazole groups are positively charged at pH 3.0, it seems likely that thermal stabilization at pH 3.0 by a His-->Ala substitution would be the result of a reduction in such unfavorable electrostatic interactions. These results suggest that amino acid substitutions that cause a decrease in the number of positive charges on the surface of a protein can be used as a general strategy to enhance protein stability at pH values below pH 4.

Contributions of the polar, uncharged amino acids to the stability of staphylococcal nuclease: evidence for mutational effects on the free energy of the denatured state

In order to quantitate the contributions of the polar, uncharged amino acids to the stability of the native state of staphylococcal nuclease, each of the 13 alanines, 9 glycines, 9 threonines, 6 prolines, 6 glutamines, 6 asparagines, and 3 serines was substituted, either with both alanine and glycine or with 1 of these 2 amino acids plus valine. For each mutant, the stability to reversible denaturation (delta GH2O) was quantitated by determining the Kapp for this reaction as a function of guanidine hydrochloride concentration. In addition, the parameter mGuHCl (= d(delta G)/d[GuHCl]) was calculated from the data. To identify the local structural features responsible for the relatively large and variable changes in delta GH2O and mGuHCl observed for the same type of substitution at different locations in nuclease, statistical correlations were sought between delta GH2O, mGuHCl, and a number of descriptors of the local structure. As with substitutions of the large hydrophobic amino acids [Shortle, D., Stites, W. E., & Meeker, A. K. (1990) Biochemistry 29, 8033-8041], mutation of polar, uncharged residues to Gly leads to a change in stability that, on average, correlates well with the degree to which the wild-type residue is buried. This correlation is especially significant for threonine, an amino acid with both polar and hydrophobic character, but is not demonstrated for the more typically hydrophobic residue alanine. As reported in the previous study of alanine/glycine substitutions of hydrophobic residues, a significant correlation between changes in stability and changes in the value of mGuHCl is again observed, strengthening the conclusion that the putative structural changes in the denatured state which lead to increases or decreases in mGuHCl are responsible for a significant fraction of the stability loss for an average mutant. The existence of this correlation is consistent with the denatured state of wild-type staphylococcal nuclease having evolved to a relatively high free energy via optimization of a balance between a maximal exposure of hydrophobic surface and a minimal gain in chain entropy. On average, mutations are less stable in proportion to the extent of which they perturb this balance. A new and puzzling correlation is reported between the extent of buriedness of a residue in the wild-type native state versus the difference in mGuHCl between the Ala mutation and the Gly mutation at that position.

Effect of cavity-modulating mutations on the stability of Escherichia coli ribonuclease HI

The size of the cavity around Ser68 of Escherichia coli ribonuclease HI was modulated by amino acid substitutions to examine the effects on the stability of the enzyme. Five mutant proteins, Ser68----Gly, Ser68----Ala, Ser68----Thr, Ser68----Val and Ser68----Leu, were constructed. Each of the mutant proteins exhibited at least 40% of the enzyme activity of the wild-type protein. The stabilities of the mutant proteins were determined from urea-denaturation and thermal-denaturation curves. Among the five mutations, only the Ser----Val mutation resulted in an increase in the stability of the enzyme. The melting temperature, tm, at pH 3.0 of the mutant protein Ser68----Val was increased by 1.9 degrees C. Its free-energy change of unfolding in the absence of urea, delta G(H2O), and the midpoint of the denaturation curve, [D]1/2, were also increased by 5.4 kJ/mol and 0.18 M, respectively. The increase in the stability of the enzyme is probably due to the filling of the cavity space around Ser68 by valine. However, the mutation of Ser68 to glycine or leucine residues resulted in a considerable decrease in stability. In these cases, some conformational changes occur, as suggested by the CD and 1H-NMR spectra of these mutant proteins.

Mutational studies of protein structures and their stabilities

The fundamental relationship between structure and function has served to guide investigations into the workings of living systems at all levels - from the whole organism to individual cells on down to individual molecules. When X-ray crystallography began to reveal the three-dimensional structures of proteins like myoglobin, lysozyme and RNase A, protein chemists were well prepared to draw inferences about functional mechanisms from the precise positioning of amino acid residues they could see. The close proximity between an amino acid side chain and a chemical group on a bound ligand strongly suggests a functional role for that side chain in binding affinity and specificity. Likewise, the nearly universal finding of large clusters of hydrophobic side chains buried in the core of proteins strongly supports a major functional role of hydrophobic interactions in protein folding and stability. Even though eminently plausible hypotheses like these, grounded in the most fundamental principles of chemistry and the logic of structure–function relationships, become widely accepted and make their way into textbooks, protein chemists have felt compelled to search for ways to test them and put them on a more quantitative basis.

Protein stability and molecular adaptation to extreme conditions

Proteins, due to the delicate balance of stabilizing and destabilizing interactions, are only marginally stable. Adaptation to extreme environments tends to shift the 'mesophilic' characteristics of proteins to the respective extremes of temperature, hydrostatic pressure, pH and salinity, such that, under the mutual physiological conditions, the molecular properties are similar regarding overall topology, flexibility and solvation. Enhanced intrinsic stability requires only minute local structural changes so that general strategies of stabilization cannot be established. Apart from mutative changes of amino-acid sequences, extrinsic factors (or cellular components) may be involved in 'extremophilic adaptation'. The molecular basis of acidophilic, alkalophilic and barophilic adaptation is still obscure. Mechanisms of enhanced thermal stability involve improved packing density, as well as specific local interactions. In halophiles, water and salt binding of the intrinsically stable protein inventory is accomplished by favoring acidic over basic amino acid residues and decreased hydrophobicity. General limits of viability are: (a) the susceptibility of the covalent structure of the polypeptide chain toward hydrolysis or hydrothermal degradation; (b) the competition of extreme solvent parameters with the weak electrostatic and hydrophobic interactions involved in protein stabilization; (c) perturbations of the folding and assembly of proteins; and (d) 'dislocation' of biochemical pathways due to effects of extreme conditions on the intricate network of metabolic reactions.

Computer modelling studies on the mechanism of action of ribonuclease T1

The mechanism of action of ribonuclease (RNase) T1 is still a matter of considerable debate as the results of x-ray, 2-D nmr and site-directed mutagenesis studies disagree regarding the role of the catalytically important residues. Hence computer modelling studies were carried out by energy minimisation of the complexes of RNase T1 and some of its mutants (His40Ala, His40Lys, and Glu58Ala) with the substrate guanyl cytosine (GpC), and of native RNase T1 with the reaction intermediate guanosine 2',3'-cyclic phosphate (G greater than p). The puckering of the guanosine ribose moiety in the minimum energy conformer of the RNase T1-GpC (substrate) complex was found to be O4'-endo and not C3'-endo as in the RNase T1-3'-guanylic acid (inhibitor/product) complex. A possible scheme for the mechanism of action of RNase T1 has been proposed on the basis of the arrangement of the catalytically important amino acid residues His40, Glu58, Arg77, and His92 around the guanosine ribose and the phosphate moiety in the RNase T1-GpC and RNase T1-G greater than p complexes. In this scheme, Glu58 serves as the general base group and His92 as the general acid group in the transphosphorylation step. His40 may be essential for stabilising the negatively charged phosphate moiety in the enzyme-transition state complex.

Histidine-40 of ribonuclease T1 acts as base catalyst when the true catalytic base, glutamic acid-58, is replaced by alanine

Mechanisms for the ribonuclease T1 (RNase T1; EC 3.1.27.3) catalyzed transesterification reaction generally include the proposal that Glu58 and His92 provide general base and general acid assistance, respectively [Heinemann, U., & Saenger, W. (1982) Nature (London) 299, 27-31]. This view was recently challenged by the observation that mutants substituted at position 58 retain high residual activity; a revised mechanism was proposed in which His40, and not Glu58, is engaged in catalysis as general base [Nishikawa, S., Morioka, H., Kim, H., Fuchimura, K., Tanaka, T., Uesugi, S., Hakoshima, T., Tomita, K., Ohtsuka, E., & Ikehara, M. (1987) Biochemistry 26, 8620-8624]. To clarify the functional roles of His40, Glu58, and His92, we analyzed the consequences of several amino acid substitutions (His40Ala, His40Lys, His40Asp, Glu58Ala, Glu58Gln, and His92Gln) on the kinetics of GpC transesterification. The dominant effect of all mutations is on Kcat, implicating His40, Glu58, and His92 in catalysis rather than in substrate binding. Plots of log (Kcat/Km) vs pH for wild-type, His40Lys, and Glu58Ala RNase T1, together with the NMR-determined pKa values of the histidines of these enzymes, strongly support the view that Glu58-His92 acts as the base-acid couple. The curves also show that His40 is required in its protonated form for optimal activity of wild-type enzyme. We propose that the charged His40 participates in electrostatic stabilization of the transition state; the magnitude of the catalytic defect (a factor of 2000) from the His40 to Ala replacement suggests that electrostatic catalysis contributes considerably to the overall rate acceleration. For Glu58Ala RNase T1, the pH dependence of the catalytic parameters suggests an altered mechanism in which His40 and His92 act as base and acid catalyst, respectively. The ability of His40 to adopt the function of general base must account for the significant activity remaining in Glu58-mutated enzymes.

Stability of recombinant Lys25-ribonuclease T1

The conformational stability of recombinant Lys25-ribonuclease T1 has been determined by differential scanning microcalorimetry (DSC), UV-monitored thermal denaturation measurements, and isothermal Gdn.HCl unfolding studies. Although rather different extrapolation procedures are involved in calculating the Gibbs free energy of stabilization, there is fair agreement between the delta G degrees values derived from the three different experimental techniques at pH 5, theta = 25 degrees C: DSC, 46.6 +/- 2.1 kJ/mol; UV melting curves, 48.7 +/- 5 kJ/mol; Gdn.HCl transition curves, 40.8 +/- 1.5 kJ/mol. Thermal unfolding of the enzyme is a reversible process, and the ratio of the van't Hoff and calorimetric enthalpy, delta HvH/delta Hcal, is 0.97 +/- 0.06. This result strongly suggests that the unfolding equilibrium of Lys25-ribonuclease T1 is adequately described by a simple two-state model. Upon unfolding the heat capacity increases by delta Cp degrees = 5.1 +/- 0.5 kJ/(mol.K). Similar values have been found for the unfolding of other small proteins. Surprisingly, this denaturational heat capacity change practically vanishes in the presence of moderate NaCl concentrations. The molecular origin of this effect is not clear; it is not observed to the same extent in the unfolding of bovine pancreatic ribonuclease A, which was employed in control experiments. NaCl stabilizes Lys25-ribonuclease T1. The transition temperature varies with NaCl activity in a manner that suggests two limiting binding equilibria to be operative. Below approximately 0.2 M NaCl activity unfolding is associated with dissociation of about one ion, whereas above that concentration about four ions are released in the unfolding reaction.(ABSTRACT TRUNCATED AT 250 WORDS)

Contribution of histidine residues to the conformational stability of ribonuclease T1 and mutant Glu-58----Ala

The pK values of the histidine residues in ribonuclease T1 (RNase T1) are unusually high: 7.8 (His-92), 7.9 (His-40), and 7.3 (His-27) [Inagaki et al. (1981) J. Biochem. 89, 1185-1195]. In the RNase T1 mutant Glu-58----Ala, the first two pK values are reduced to 7.4 (His-92) and 7.1 (His-40). These lower pKs were expected since His-92 (5.5 A) and His-40 (3.7 A) are in close proximity to Glu-58 at the active site. The conformational stability of RNase T1 increases by over 4 kcal/mol between pH 9 and 5, and this can be entirely accounted for by the greater affinity for protons by the His residues in the folded protein (average pK = 7.6) than in the unfolded protein (pk approximately 6.6). Thus, almost half of the net conformational stability of RNase T1 results from a difference between the pK values of the histidine residues in the folded and unfolded conformations. In the Glu-58----Ala mutant, the increase in stability between pH 9 and 5 is halved (approximately 2 kcal/mol), as expected on the basis of the lower pK values for the His residues in the folded protein (average pK = 7.1). As a consequence, RNase T1 is more stable than the mutant below pH 7.5, and less stable above pH 7.5. These results emphasize the importance of measuring the conformational stability as a function of pH when comparing proteins differing in structure.

Purification of recombinant ribonuclease T1 expressed in Escherichia coli

A protocol for the rapid purification of ribonuclease T1 expressed from a chemically synthesized gene cloned into Escherichia coli is described. QAE ion-exchange and Sephadex G-50 chromatography are used to give over 300 mg (88% yield) of pure ribonuclease T1 from 61 of liquid culture in 3 days. We also report a new absorption coefficient for RNase T1: E1%278 nm = 15.4.

Conformational stability of globular proteins

The conformational stability of ribonuclease T1 has been measured as a function of the variables of most interest to biochemists: temperature, pH, salt concentration, disulfide-bond content and amino acid sequence. The results provide insight into the forces that stabilize globular proteins.