A binding free energy hot spot in the ankyrin repeat protein GABPbeta mediated protein-protein interaction

The frequently observed ankyrin repeat motif represents a structural scaffold evolved for mediating protein-protein interactions. As such, these repeats modulate a diverse range of cellular functions. We thermodynamically characterized the heterodimeric GA-binding protein (GABP) alphabeta complex and focused specifically on the interaction mediated by the ankyrin repeat domain of the GABPbeta. Our isothermal titration calorimetric analysis of the interaction between the GABP subunits determined an association constant (K(A)) of 6.0 x 10(8) M(-1) and that the association is favorably driven by a significant change in enthalpy (DeltaH) and a minor change in entropy (-TDeltaS). A total of 16 GABPbeta interface residues were chosen for alanine scanning mutagenesis. The calorimetrically measured differences in the free energy of binding were compared to computationally calculated values resulting in a correlation coefficient r = 0.71. We identified three spatially contiguous hydrophobic and aromatic residues that form a binding free energy hot spot (DeltaDeltaG > 2.0 kcal/mol). One residue provides structural support to the hot spot residues. Three non-hot spot residues are intermediate contributors (DeltaDeltaG approximately 1.0 kcal/mol) and create a canopy-like structure over the hot spot residues to possibly occlude solvent and orientate the subunits. The remaining interface residues are located peripherally and have weak contributions. Finally, our mutational analysis revealed a significant entropy-enthalpy compensation for this interaction.     

SUPREX (Stability of Unpurified Proteins from Rates of H/D Exchange) analysis of the thermodynamics of synergistic anion binding by ferric-binding protein (FbpA), a bacterial transferrin

SUPREX (stability of unpurified proteins from rates of H/D exchange) is a H/D exchange- and matrix-assisted laser desorption/ionization (MALDI)-based technique for characterizing the equilibrium unfolding/refolding properties of proteins and protein-ligand complexes. Here, we describe the application of SUPREX to the thermodynamic analysis of synergistic anion binding to iron-loaded ferric-binding protein (Fe(3+)FbpA-X, X = synergistic anion). The in vivo function of FbpA is to transport unchelated Fe(3+) across the periplasmic space of certain Gram-negative bacteria, a process that requires simultaneous binding of a synergistic anion. Our results indicate that Fe(3+)FbpA-X is not a so-called "ideal" protein system for SUPREX analyses because it does not exhibit two-state folding properties and it does not exhibit EX2 H/D exchange behavior. However, despite these nonideal properties of the Fe(3+)FbpA-X protein-folding/unfolding reaction, we demonstrate that the SUPREX technique is still amenable to the quantitative thermodynamic analysis of synergistic anion binding to Fe(3+)FbpA. As part of this work, the SUPREX technique was used to evaluate the DeltaDeltaG(f) values of four synergistic anion-containing complexes of Fe(3+)FbpA (i.e., Fe(3+)FbpA-PO(4), Fe(3+)FbpA-citrate, Fe(3+)FbpA-AsO(4), and Fe(3+)FbpA-SO(4)). The DeltaDeltaG(f) value obtained for Fe(3+)FbpA-citrate relative to Fe(3+)FbpA-PO(4) (1.45 +/- 0.44 kcal/mol), is in good agreement with that reported previously (1.98 kcal/mol). The value obtained for Fe(3+)FbpA-AsO(4) (0.58 +/- 0.45 kcal/mol) was also consistent with that reported previously (0.68 kcal/mol), but the measurement error is very close to the magnitude of the value. This work (i) demonstrates the utility of the SUPREX method for studying anion binding by FbpA, (ii) provides the first evaluation of a DeltaDeltaG(f) value for Fe(3+)FbpA-SO(4), -1.43 +/- 0.17 kcal/mol, and (iii) helps substantiate our hypothesis that the synergistic anion plays a role in controlling the lability of iron bound to FbpA in the transport process.     

Alanine scanning mutagenesis of Abeta(1-40) amyloid fibril stability

We describe here an alanine scanning mutational analysis of the Abeta(1-40) amyloid fibril monitored by fibril elongation thermodynamics derived from critical concentration values for fibril growth. Alanine replacement of most residues in the amyloid core region, residues 15-36, leads to destabilization of the elongation step, compared to wild-type, by about 1kcal/mol, consistent with a major role for hydrophobic packing in Abeta(1-40) fibril assembly. Where comparisons are possible, the destabilizing effects of Ala replacements are generally in very good agreement with the effects of Ala replacements of the same amino acid residues in an element of parallel beta-sheet in the small, globular protein Gbeta1. We utilize these Ala-WT DeltaDeltaG values to filter previously described Pro-WT DeltaDeltaG values, creating Pro-Ala DeltaDeltaG values that specifically assess the sensitivity of a sequence position, in the structural context of the Abeta fibril, to replacement by proline. The results provide a conservative view of the energetics of Abeta(1-40) fibril structure, indicating that positions 18-21, 25-26, and 32-33 within amyloid structure are particularly sensitive to the main-chain disrupting effects of Pro replacements. In contrast, residues 14-17, 22, 24, 27-31, and 34-39 are relatively insensitive to Pro replacements; most N-terminal residues were not tested. The results are discussed in terms of amyloid fibril structure and folding energetics, in particular focusing on how the data compare to those from other structural studies of Abeta(1-40) amyloid fibrils grown in phosphate-buffered saline at 37 degrees C under unstirred ("quiescent") conditions.     

DNA helicase RepA: cooperative ATPase activity and binding of nucleotides

The steady-state kinetic parameters of the ATPase activity of the homohexameric DNA helicase RepA and the binding of the fluorescent analogue epsilonADP to RepA have been studied. ssDNA stimulates RepA ATPase activity optimally at acidic pH 5.3-6.0. The sigmoidal kinetic curves in both the absence and presence of ssDNA show strong positive cooperativity for ATP hydrolysis, with oligonucleotides longer than 10mer optimal for ssDNA-stimulated ATPase activity. Fluorescence titrations show that, at 25 degrees C and in the absence of DNA, the binding of epsilonADP to RepA is biphasic with three high (K(1) = 1.54 x 10(6) M(-1)) and three low (K(2) = 4.71 x 10(4) M(-)(1)) affinity binding sites differing by 30-40-fold in binding constants. In the absence of cofactors, RepA melts cooperatively at T(m) = 65.8 +/- 0.1 degrees C and is more stable in the presence of ATPgammaS, T(m) = 68.1 +/- 0.2 degrees C (DeltaDeltaG 0.95 kcal/mol), than in the presence of ADP, T(m) = 66. 5 +/- 0.1 degrees C (DeltaDeltaG 0.29 kcal/mol), indicating that the additional phosphate group in ATPgammaS has a significant influence on RepA structure. A model is proposed in which individual subunits of RepA sequentially and cooperatively perform a multistep ATP hydrolytic cycle.     

Analysis of the interactions of human ribonuclease inhibitor with angiogenin and ribonuclease A by mutagenesis: importance of inhibitor residues inside versus outside the C-terminal "hot spot"

Ribonuclease inhibitor (RI) binds diverse mammalian RNases with extraordinary avidity. Here, we have investigated the structural basis for this tight binding and broad specificity by mutational analysis of the complexes of RI with angiogenin (Ang) and RNase A (K(D)=0.5 fM and 43 fM, respectively). Both crystal structures are known; the interfaces are large, and the ligands dock similarly, although few of the specific interactions formed are analogous. Our previous mutagenesis studies focused primarily on one contact region, containing RI 434-438 and the enzymatic active site. Many single-residue replacements produced extensive losses of binding energy (2.3-5.9 kcal/mol), suggesting that this region constitutes a "hot spot" in both cases. We have now explored the roles of most of the remaining RI residues that interact with Ang and/or RNase A. One major cluster in each complex lies in a Trp-rich area of RI, containing Trp261, Trp263, Trp318, and Trp375. Although the energy losses from individual replacements in this portion of the Ang complex were small-to-moderate (0-1.5 kcal/mol), the changes from multiple substitutions were much greater than additive, and the binding energy provided by this region is estimated to be approximately 6 kcal/mol (30 % of total). Effects of replacing combinations of hot spot components had also been found to be superadditive, and this negative cooperativity is now shown to extend to the neighboring contact residue RI Ser460. The overall contribution of the hot spot, taking superadditivity into account, is then approximately 14-15 kcal/mol. The hot spot and Trp-rich regions, although spatially well separated, are themselves functionally linked. No other parts of the RI-Ang interface appear to be energetically important. Binding of RNase A is more sensitive to substitutions throughout the interface, with free energy losses>/=1 kcal/mol produced by nearly all replacements examined, so that the sum of losses greatly exceeds the binding energy of the complex. This discrepancy can be explained, in part, by positive cooperativity, as evident from the subadditive effects observed when combinations of residues in either the hot spot or Trp-rich region are replaced. These findings suggest that the binding energy may be more widely distributed in the RNase A complex than in the Ang complex.     

Calorimetric examination of high-affinity Src SH2 domain-tyrosyl phosphopeptide binding: dissection of the phosphopeptide sequence specificity and coupling energetics

SH2 domains are protein modules which interact with specific tyrosine phosphorylated sequences in target proteins. The SH2 domain of the Src kinase binds with high affinity to a tyrosine phosphorylated peptide containing the amino acids Glu, Glu, and Ile (EEI) at the positions +1, +2, and +3 C-terminal to the phosphotyrosine, respectively. To investigate the degree of selectivity of the Src SH2 domain for each amino acid of the EEI motif, the binding thermodynamics of a panel of substitutions at the +1 (Gln, Asp, Ala, Gly), +2 (Gln, Asp, Ala, Gly), and +3 (Leu, Val, Ala, Gly) positions were examined using titration microcalorimetry. It was revealed that the Src SH2 domain is insensitive (DeltaDeltaG degrees 相似文献   

Impact of proline residues on parvalbumin stability

Despite its higher net charge and reduced opportunities for favorable tertiary interactions, Ca(2+)-free rat beta-parvalbumin is more stable than rat alpha-parvalbumin. Under conditions wherein alpha denatures at 45.8 degrees C, beta denatures at 53.6 degrees. The homologous chicken beta isoform known as CPV3 also exhibits heightened stability-prompting an inquiry into the stabilizing influence of Pro-21 and Pro-26. Individual P21A and P26A mutations lower the T(m) of rat beta by 3.2 degrees, decreasing conformational stability by 0.74 kcal/mol. Simultaneous replacement of Pro-21 and Pro-26 essentially abolishes the excess stability (DeltaT(m) = -7.6 degrees; DeltaDeltaG(conf) = -1.77 kcal/mol). Significantly, the P21A/P26A variant displays Ca(2+) affinity virtually indistinguishable from wild-type beta, implying that structural alterations in the AB domain do not necessarily influence the divalent ion affinity of the CD-EF domain. The consequences of introducing proline at positions 21 and 26 in rat alpha were also examined. Whereas the H26P mutation raises the T(m) by 5.6 degrees (DeltaDeltaG(conf) = 1.25 kcal/mol), A21P lowers the T(m) by 8.5 degrees (DeltaDeltaG(conf) = -1.9 kcal/mol). Replacement of Ala-21 by proline in an alpha AB/beta CD-EF chimera increases the T(m) by 5.8 degrees (DeltaDeltaG(conf) = 0.95 kcal/mol), implying that the destabilization of alpha by Pro-21 results from steric conflict with a residue in the CD-EF domain. Consistent with that hypothesis, the K80S mutation markedly stabilizes alpha A21P, yielding a protein with a T(m) 2.0 degrees higher than wild-type alpha. The observed differences in stability resulting from proline addition/removal are largely consistent with alterations in main-chain and side-chain conformational entropy.     

How protein stability and new functions trade off

Numerous studies have noted that the evolution of new enzymatic specificities is accompanied by loss of the protein's thermodynamic stability (DeltaDeltaG), thus suggesting a tradeoff between the acquisition of new enzymatic functions and stability. However, since most mutations are destabilizing (DeltaDeltaG>0), one should ask how destabilizing mutations that confer new or altered enzymatic functions relative to all other mutations are. We applied DeltaDeltaG computations by FoldX to analyze the effects of 548 mutations that arose from the directed evolution of 22 different enzymes. The stability effects, location, and type of function-altering mutations were compared to DeltaDeltaG changes arising from all possible point mutations in the same enzymes. We found that mutations that modulate enzymatic functions are mostly destabilizing (average DeltaDeltaG = +0.9 kcal/mol), and are almost as destabilizing as the "average" mutation in these enzymes (+1.3 kcal/mol). Although their stability effects are not as dramatic as in key catalytic residues, mutations that modify the substrate binding pockets, and thus mediate new enzymatic specificities, place a larger stability burden than surface mutations that underline neutral, non-adaptive evolutionary changes. How are the destabilizing effects of functional mutations balanced to enable adaptation? Our analysis also indicated that many mutations that appear in directed evolution variants with no obvious role in the new function exert stabilizing effects that may compensate for the destabilizing effects of the crucial function-altering mutations. Thus, the evolution of new enzymatic activities, both in nature and in the laboratory, is dependent on the compensatory, stabilizing effect of apparently "silent" mutations in regions of the protein that are irrelevant to its function.     

Modeling implicit reorganization in continuum descriptions of protein-protein interactions

The determination of free energies that govern protein-protein recognition is essential for a detailed molecular understanding of biological specificity. Continuum models of macromolecular interactions, in which the solvent is treated by an implicit representation and the proteins are treated semi-microscopically, are computationally tractable for estimating free energies, yet many questions remain concerning their accuracy. This article reports a continuum analysis of the free-energy changes underlying the binding of 31 interfacial alanine substitutions of two complexes of the antihen egg white lysozyme (HEL) antibody D1.3 bound with HEL or the antibody E5.2. Two implicit schemes for modeling the effects of protein and solvent relaxation were examined, in which the protein environment was treated as either homogeneous with a "protein dielectric constant" of epsilon(p) = 4 or inhomogeneous, with epsilon(p) = 4 for neutral residues and epsilon(p) = 25 for ionized residues. The results showed that the nonuniform dielectric model reproduced the experimental differences better, with an average absolute error of +/-1.1 kcal/mol, compared with +/-1.4 kcal/mol for the uniform model. More importantly, the error for charged residues in the nonuniform model is +/-0.8 kcal/mol and is nearly half of that corresponding to the uniform model. Several substitutions were clearly problematic in determining qualitative trends and probably required explicit structural reorganization at the protein-protein interface.     

The stability of transmembrane helix interactions measured in a biological membrane

Despite some promising progress in the understanding of membrane protein folding and assembly, there is little experimental information regarding the thermodynamic stability of transmembrane helix interactions and even less on the stability of transmembrane helix-helix interactions in a biological membrane. Here we describe an approach that allows quantitative measurement of transmembrane helix interactions in a biological membrane, and calculation of changes in the interaction free energy resulting from substitution of single amino acids. Dimerization of several variants of the glycophorin A transmembrane domain are characterized and compared to the wild-type (wt) glycophorin A transmembrane helix dimerization. The calculated DeltaDeltaG(app) values are further compared with values found in the literature. In addition, we compare interactions between the wt glycophorin A transmembrane domain and helices in which critical glycine residues are replaced by alanine or serine, respectively. The data demonstrate that replacement of the glycine residues by serine is less destabilizing than replacement by alanine with a DeltaDeltaG(app) value of about 0.4 kcal/mol. Our study comprises the first measurement of a transmembrane helix interaction in a biological membrane, and we are optimistic that it can be further developed and applied.     

A general method for the quantitative analysis of functional chimeras: applications from site-directed mutagenesis and macromolecular association

Two new parameters, I: and C:, are introduced for the quantitative evaluation of functional chimeras: I: (impact) and C: (context dependence) are the free energy difference and sum, respectively, of the effects on a given property measured in forward and retro chimeras. The forward chimera is made by substitution of a part "a" from ensemble A into the analogous position of homologous ensemble B (S:(B --> A)). The C: value is a measure of the interaction of the interrogated position with its surroundings, whereas I: is an expression of the quantitative importance of the probed position. Both I: and C: vary with the evaluated property, for example, kinetics, binding, thermostability, and so forth. The retro chimera is the reverse substitution of the analogous part "b" from B into A, S:(A --> B). The I: and C: values derived from original data for forward and retro mutations in aspartate and tyrosine aminotransferase, from literature data for quasi domain exchange in oncomodulin and for the interaction of Tat with bovine and human TAR are evaluated. The most salient derived conclusions are, first, that Thr 109 (AATase) or Ser 109 (TATase) is an important discriminator for dicarboxylic acid selectivity by these two enzymes (I: < -2.9 kcal/mol). The T109S mutation in AATase produces a nearly equal and opposite effect to S109T in TATase (C: < 0.4 kcal/mol). Second, an I: value of 5.5 kcal/mol describes the effects of mirror mutations D94S (site 1) and S55D (site 2) in the Ca(2+) binding sites of oncomodulin on Ca(2+) affinity. The second mirror set, G98D (site 1) and D59G (site 2), yields a smaller impact (I: = -3.4 kcal/mol) on Ca(2+) binding; however, the effect is significantly more nearly context independent (C: = -0.6 versus C: = -2.7 kcal/mol). Third, the stem and loop regions of HIV and BIV TAR are predominantly responsible for the species specific interaction with BIV Tat(65-81) (I: = -1.5 to -1.6 kcal/mol), whereas I: = 0.1 kcal/mol for bulge TAR chimeras. The C: values are from -0.3 to -1.2 kcal/mol. The analysis described should have important applications to protein design.     

Quantitative evaluation of the chicken lysozyme epitope in the HyHEL-10 Fab complex: free energies and kinetics.

The hen (chicken) egg-white lysozyme (HEWL) epitope for the monoclonal antibody HyHEL-10 Fab (Fab-10) was investigated by alanine scan mutagenesis. The association rate constants (k(on)) for the HEWL Fab-10 complexes were obtained from the homogenous solution method described in the preceding paper (Taylor et al., 1998). A new method for determining the dissociation rate constant (k(off)) for the complex, by trapping nascent free antibody with an inactive HEWL mutant is described. The values of k(on) fall within a factor of 2 of the wild-type (WT) HEWL value (1.43+/-0.13 X 10(6)M(-1)s(-1)), while the increases in k(off)more nearly reflect the total change in free energies of the complex (deltadeltaG(D)). The dissociation constants (K(D)) were measured directly in those cases where satisfactory kinetic data could not be obtained. The Y20A, K96A, and K97A HEWL.Fab-10 complexes are destabilized by more than 4 kcal/mol compared to the WT complex. The R21A, L75A, and D101A antibody complexes are moderately destabilized (0.7 < deltadeltaG(D)< or = 1.0 kcal/mol). Additional mutations of the "hotspot" residues (Tyr20, Lys96, Lys97) were constructed to probe, more precisely, the nature of their contributions to complex formation. The results show that the entire hydrocarbon side chains of Tyr20 and Lys97, and only the epsilon-amino group of Lys96, contribute to the stability of the complex. The value of deltadeltaG(D) for the R21A mutant complex is a distinct outlier in the Arg21 replacement series demonstrating the importance of supplementing alanine scan mutagenesis with additional mutations.     

Evaluation of direct and cooperative contributions towards the strength of buried hydrogen bonds and salt bridges

An experimental approach to evaluate the net binding free energy of buried hydrogen bonds and salt bridges is presented. The approach, which involves a modified multiple-mutant cycle protocol, was applied to selected interactions between TEM-1-beta-lactamase and its protein inhibitor, BLIP. The selected interactions (two salt bridges and two hydrogen bonds) all involving BLIP-D49, define a distinct binding unit. The penta mutant, where all side-chains constructing the binding unit were mutated to Ala, was used as a reference state to which combinations of side-chains were introduced. At first, pairs of interacting residues were added allowing the determination of interaction energies in the absence of neighbors, using double mutant cycles. Addition of neighboring residues allowed the evaluation of their cooperative effects on the interaction. The two isolated salt bridges were either neutral or repulsive whereas the two hydrogen bonds contribute 0.3 kcal mol(-1 )each. Conversely, a double mutant cycle analysis of these interactions in their native environment showed that they all stabilize the complex by 1-1.5 kcal mol(-1). Examination of the effects of neighboring residues on each of the interactions revealed that the formation of a salt bridge triad, which involves two connected salt bridges, had a strong cooperative effect on stabilizing the complex independent of the presence or absence of additional neighbors. These results demonstrate the importance of forming net-works of buried salt bridges. We present theoretical electrostatic calculations which predict the observed mode of cooperativity, and suggest that the cooperative networking effect results from the favorable contribution of the protein to the interaction. Furthermore, a good correlation between calculated and experimentally determined interaction energies for the two salt bridges, and to a lesser extent for the two hydrogen bonds, is shown. The data analysis was performed on values of DeltaDeltaG(double dagger)K(d) which reflect the strength of short range interactions, while DeltaDeltaG(o)K(D) values which include the effects of long range electrostatic forces that alter specifically DeltaDeltaG(double dagger)k(a) were treated separately.     

In-silico epitope identification and design of Uricase mutein with reduced immunogenicity

The clinical utilization of Uricase against gout is limited due to the immunogenicity. In the present article, we identified the antigenic determinants of Uricase and reduced their immunogenicity via in-silico mutagenesis. Multiple sequence alignment and motif analysis were carried out to identify the conserved residues in evolutionary process. Emini surface accessibility, Parker hydrophilicity, and Karplus & Schulz flexibility methods were employed to predict the linear B-cell epitopes of both Ag-Uricase and Bf-Uricase. Deimmunization approach identified T-cell epitopes and the hot spot residues. Reduced antigenic probability was obtained in case of T159W, D169C, N264W and Y203D mutations for Ag-Uricase, while S139 V, K215W, G216 F and I172 P mutations for Bf-Uricase. The binding affinity values of uric acid towards the catalytic pocket of Ag-Uricase and Bf-Uricase models were found to be -48.71 kcal/mol and -40.93 kcal/mol, respectively. This energy is further stabilized in the mutant model by -6.36 kcal/mol and -1.45 kcal/mol for Ag-Uricase and Bf-Uricase, respectively. About 100 ns molecular dynamics simulation was performed to evaluate the conformational stability of both native and mutated Uricase. Insights obtained from this study provide guidelines for experimental design of Uricase muteins with reduced antigenicity.     

The vnd/NK-2 homeodomain: thermodynamics of reversible unfolding and DNA binding for wild-type and with residue replacements H52R and H52R/T56W in helix III

The conformational stabilities of the vnd (ventral nervous system defective)/NK-2 homeodomain [HD(wt); residues 1-80 that encompass the 60-residue homeodomain] and those harboring mutations in helix III of the DNA recognition site [HD(H52R) and HD(H52R/T56W)] have been investigated by differential scanning calorimetry (DSC) and ellipticity changes at 222 nm. Thermal unfolding reactions at pH 7.4 are reversible and repeatable in the presence of 50-500 mM NaCl with DeltaC(p) = 0.52 +/- 0.04 kcal K(-1) mol(-1). A substantial stabilization of HD(wt) is produced by 50 mM phosphate or by the addition of 100-500 mM NaCl to 50 mM Hepes, pH 7.4, buffer (from T(m) = 35.5 degrees C to T(m) 43-51 degrees C; DeltaH(vH) congruent with 47 +/- 5 kcal mol(-1)). The order of stability is HD(H52R/T56W) > HD(H52R) > HD(wt), irrespective of the anions present. Progress curves for ellipticity changes at 222 nm as a function of increasing temperature are fitted well by a two-state unfolding model, and the cooperativity of secondary structure changes is greater for mutant homeodomains than for HD(wt) and also is increased by adding 100 mM NaCl to Hepes buffer. A 33% quench of the intrinsic tryptophanyl residue fluorescence of HD(wt) by phosphate binding (K(D)' = 2.6 +/- 0.3 mM phosphate) is reversed approximately 60% by DNA binding. Thermodynamic parameters for vnd/NK-2 homeodomain proteins binding sequence-specific 18 bp DNA have been determined by isothermal titration calorimetry (10-30 degrees C). Values of DeltaC(p) are +0.25, -0.17, and -0.10 +/- 0.04 kcal K(-1) mol(-1) for HD(wt), HD(H52R), and HD(H52R/T56W) binding duplex DNA, respectively. Interactions of homeodomains with DNA are enthalpically controlled at 298 K and pH 7.4 with corresponding DeltaH values of -6.6 +/- 0.5, -10.8 +/- 0.1, and -9.0 +/- 0.6 kcal mol(-1) and DeltaG' values of -11.0 +/- 0.1, -11.0 +/- 0.1, and -11.3 +/- 0.3 kcal mol(-1) with a binding stoichiometry of 1.0 +/- 0.1. Thermodynamic parameters for DNA binding are not predicted from homeodomain structural changes that occur upon complexing to DNA and must reflect also solvent and possibly DNA rearrangements.     

Thermodynamic stability of the bacteriorhodopsin lattice as measured by lipid dilution

To determine the strength of noncovalent interactions that stabilize a membrane protein complex, we have developed an in vitro method for quantifying the dissociation of the bacteriorhodopsin (BR) lattice, a naturally occurring two-dimensional crystal. A lattice suspension was titrated with a short- and long-chain phosphatidylcholine mixture to dilute BR within the lipid bilayer. The fraction of BR in the lattice form as a function of added lipid was determined by visible circular dichroism spectroscopy and fit with a cooperative self-assembly model to obtain a critical concentration for lattice assembly. Critical concentration values of wild-type and mutant proteins were used to calculate the change in lattice stability upon mutation (DeltaDeltaG). By using this method, a series of mutant proteins was examined in which residues at the BR-BR interface were replaced with smaller amino acids, either Ala or Gly. Most of the mutant lattices were destabilized, with DeltaDeltaG values of 0.2-1.1 kcal/mol at 30 degrees C, consistent with favorable packing of apolar residues in the membrane. One mutant, I45A, was stabilized by approximately 1.0 kcal/mol, possibly due to increased lipid entropy. The DeltaDeltaG values agreed well with previous in vivo measurements, except in the case of I45A. The ability to measure the change in stability of mutant protein complexes in a lipid bilayer may provide a means of determining the contributions of specific protein-protein and protein-lipid interactions to membrane protein structure.     

Contribution of structural peculiarities of onconase to its high stability and folding kinetics

Onconase (ONC) from Rana pipiens is the smallest member of the ribonuclease A (RNase A) superfamily. Despite a tertiary structure similar to RNase A, ONC is distinguished by an extremely high thermodynamic stability. In the present paper we have probed the significance of three structural regions, which exhibit structural peculiarities in comparison to RNase A, for the stability of ONC to temperature and guanidine hydrochloride induced denaturation: (i) the N-terminal pyroglutamate residue, (ii) the hydrophobic cluster between helix I and the first beta-sheet, and (iii) the C-terminal disulfide bond. For this purpose, the enzyme variants 相似文献   

DeltaG-based prediction and experimental confirmation of SYCRP1-binding sites on the Synechocystis genome

DNA-binding sites for SYCRP1, which is a regulatory protein of the cyanobacterium Synechocystissp. PCC6803, were predicted for the whole genome sequence by estimating changes in the binding free energy () for SYCRP1 for those sites. The values were calculated by summing DeltaDeltaG values derived from systematic single base-pair substitution experiments (symmetrical and cooperative binding model). Of the calculated binding sites, 23 sites with a value <3.9kcal.mol(-1) located upstream or between the ORFs were selected as putative binding sites for SYCRP1. In order to confirm whether SYCRP1 actually binds to these binding sites or not, 11 sites with the lowest values were tested experimentally, and we confirmed that SYCRP1 binds to ten of the 11 sites with a DeltaDeltaG(total) value <3.9kcal.mol(-1). The best correlation coefficient between and the observed DeltaDeltaG(total) for binding of SYCRP1 to those sites was 0.78. These results suggest that the DeltaDeltaG values derived from systematic single base-pair experiments may be used to screen for potential binding sites of a regulatory protein in the genome sequence.     

Electrostatic contributions to T4 lysozyme stability: solvent-exposed charges versus semi-buried salt bridges

We carried our Poisson-Boltzmann (PB) calculations for the effects of charge reversal at five exposed sites (K16E, R119E, K135E, K147E, and R154E) and charge neutralization and proton titration of the H31-D70 semi-buried salt bridge on the stability of T4 lysozyme. Instead of the widely used solvent-exclusion (SE) surface, we used the van der Waals (vdW) surface as the boundary between the protein and solvent dielectrics (a protocol established in our earlier study on charge mutations in barnase). By including residual charge-charge interactions in the unfolded state, the five charge reversal mutations were found to have DeltaDeltaG(unfold) from -1.6 to 1.3 kcal/mol. This indicates that the variable effects of charge reversal observed by Matthews and co-workers are not unexpected. The H31N, D70N, and H31N/D70N mutations were found to destabilize the protein by 2.9, 1.3, and 1.6 kcal/mol, and the pK(a) values of H31 and D70 were shifted to 9.4 and 0.6, respectively. These results are in good accord with experimental data of Dahlquist and co-workers. In contrast, if the SE surface were used, the H31N/D70N mutant would be more stable than the wild-type protein by 1.3 kcal/mol. From these and additional results for 27 charge mutations on five other proteins, we conclude that 1) the popular view that electrostatic interactions are generally destabilizing may have been based on overestimated desolvation cost as a result of using the SE surface as the dielectric boundary; and 2) while solvent-exposed charges may not reliably contribute to protein stability, semi-buried salt bridges can provide significant stabilization.