Structural genomic projects envision almost routine protein structure determinations, which are currently imaginable only for small proteins with molecular weights below 25,000 Da. For larger proteins, structural insight can be obtained by breaking them into small segments of amino acid sequences that can fold into native structures, even when isolated from the rest of the protein. Such segments are autonomously folding units (AFU) and have sizes suitable for fast structural analyses. Here, we propose to expand an intuitive procedure often employed for identifying biologically important domains to an automatic method for detecting putative folded protein fragments. The procedure is based on the recognition that large proteins can be regarded as a combination of independent domains conserved among diverse organisms. We thus have developed a program that reorganizes the output of BLAST searches and detects regions with a large number of similar sequences. To automate the detection process, it is reduced to a simple geometrical problem of recognizing rectangular shaped elevations in a graph that plots the number of similar sequences at each residue of a query sequence. We used our program to quantitatively corroborate the premise that segments with conserved sequences correspond to domains that fold into native structures. We applied our program to a test data set composed of 99 amino acid sequences containing 150 segments with structures listed in the Protein Data Bank, and thus known to fold into native structures. Overall, the fragments identified by our program have an almost 50% probability of forming a native structure, and comparable results are observed with sequences containing domain linkers classified in SCOP. Furthermore, we verified that our program identifies AFU in libraries from various organisms, and we found a significant number of AFU candidates for structural analysis, covering an estimated 5 to 20% of the genomic databases. Altogether, these results argue that methods based on sequence similarity can be useful for dissecting large proteins into small autonomously folding domains, and such methods may provide an efficient support to structural genomics projects.

Physicochemical properties are potentially useful in predicting functional differences between aligned protein subfamilies. We present a method that considers physicochemical properties from ancestral sequences predicted to have given rise to the subfamilies of interest by gene duplication. Comparison between two map kinases subfamilies, p38 and ERK, revealed a region that had an excess of change in properties after gene duplication followed by conservation within the two subfamilies. This region corresponded to that experimentally defined as important for substrate and pathway specificity. The derived scores for the region of interest were found to differ significantly in their distribution compared to the rest of the protein when the Kolmogorov–Smirnov test was applied (p = 0.005). Thus, the incorporation of ancestral physicochemical properties is useful in predicting functional differences between protein subfamilies. In addition, the method was applied to the MKK and MAPK components of the p38 and JNK pathways. These proteins showed a similar pattern in their evolution and regions predicted to confer functional differences are discussed.

In a general approach to the understanding of protein adaptation to high temperature, molecular models of the closely related mesophilic Streptomyces sp. S38 Xyl1 and thermophilic Thermomonospora fusca TfxA family 11 xylanases were built and compared with the three-dimensional (3D) structures of homologous enzymes. Some of the structural features identified as potential contributors to the higher thermostability of TfxA were introduced in Xyl1 by site-directed mutagenesis in an attempt to improve its thermostability and thermophilicity. A new Y11–Y16 aromatic interaction, similar to that present in TfxA and created in Xyl1 by the T11Y mutation, improved both the thermophilicity and thermostability. Indeed, the optimum activity temperature (70 vs. 60 °C) and the apparent Tm were increased by about 9 °C, and the mutant was sixfold more stable at 57 °C. The combined mutations A82R/F168H/N169D/Δ170 potentially creating a R82–D169 salt bridge homologous to that present in TfxA improved the thermostability but not the thermophilicity. Mutations R82/D170 and S33P seemed to be slightly destabilizing and devoid of influence on the optimal activity temperature of Xyl1. Structural analysis revealed that residues Y11 and Y16 were located on β-strands B1 and B2, respectively. This interaction should increase the stability of the N-terminal part of Xyl1. Moreover, Y11 and Y16 seem to form an aromatic continuum with five other residues forming putative subsites involved in the binding of xylan (+3, +2, +1, −1, −2). Y11 and Y16 might represent two additional binding subsites (−3, −4) and the T11Y mutation could thus improve substrate binding to the enzyme at higher temperature and thus the thermophilicity of Xyl1.

Sac7d unfolds at low pH in the absence of salt, with the greatest extent of unfolding obtained at pH 2. We have previously shown that the acid unfolded protein is induced to refold by decreasing the pH to 0 or by addition of salt (McCrary BS, Bedell J, Edmondson SP, Shriver JW, 1998, J Mol Biol 276:203–224). Both near-ultraviolet circular dichroism spectra and ANS fluorescence enhancements indicate that the acid- and salt-induced folded states have a native fold and are not molten globular. 1H,15N heteronuclear single quantum coherence NMR spectra confirm that the native, acid-, and salt-induced folded states are essentially identical. The most significant differences in amide 1H and 15N chemical shifts are attributed to hydrogen bonding to titrating carboxyl side chains and through-bond inductive effects. The 1H NMR chemical shifts of protons affected by ring currents in the hydrophobic core of the acid- and salt-induced folded states are identical to those observed in the native. The radius of gyration of the acid-induced folded state at pH 0 is shown to be identical to that of the native state at pH 7 by small angle X-ray scattering. We conclude that acid-induced collapse of Sac7d does not lead to a molten globule but proceeds directly to the native state. The folding of Sac7d as a function of pH and anion concentration is summarized with a phase diagram that is similar to those observed for other proteins that undergo acid-induced folding except that the A-state is encompassed by the native state. These results demonstrate that formation of a molten globule is not a general property of proteins that are refolded by acid.

The potential for engineering stable proteins with multiple amino acid substitutions was explored. Eleven lysine, five methionine, two tryptophan, one glycine, and three threonine substitutions were simultaneously made in barley chymotrypsin inhibitor-2 (CI-2) to substantially improve the essential amino acid content of the protein. These substitutions were chosen based on the three-dimensional structure of CI-2 and an alignment of homologous sequences. The initial engineered protein folded into a wild-type-like structure, but had a free energy of unfolding of only 2.2 kcal/mol, considerably less than the wild-type value of 7.5 kcal/mol. Restoration of the lysine mutation at position 67 to the wild-type arginine increased the free energy of unfolding to 3.1 kcal/mol. Subsequent cysteine substitutions at positions 22 and 82 resulted in disulfide bond formation and a protein with nearly wild-type thermodynamic stability (7.0 kcal/mol). None of the engineered proteins retained inhibitory activity against chymotrypsin or elastase, and all had substantially reduced inhibitory activity against subtilisin. The proteolytic stabilities of the proteins correlated with their thermodynamic stabilities. Reduction of the disulfide bond resulted in substantial loss of both thermodynamic and proteolytic stabilities, confirming that the disulfide bond, and not merely the cysteine substitutions, was responsible for the increased stability. We conclude that it is possible to replace over a third of the residues in CI-2 with minimal disruption of stability and structural integrity.

The ecdysone receptor (ECR), a nuclear transcription factor controlling insect development, is a novel target for insecticides such as dibenzoylhydrazines with low environmental and toxicological impacts. To understand the high selectivity of such synthetic molecules toward ECR, two homology models of the Chironomus tentans ECR ligand-binding domain (LDB) have been constructed by taking as templates the known LBD crystal structures of the retinoic acid and vitamin D receptors. Docking of 20-hydroxyecdysone (20E) and dibenzoylhydrazines to the receptor suggests a novel superposition of the natural and synthetic molecules; the N-tert-butyl substituent of the dibenzoylhydrazines extends significantly beyond the 20E volume. Our ECR–LBD protein models rationalize how 20E and dibenzoylhydrazines interact with the ligand-binding pocket. The homology model complexes provide new insights that can be exploited in the rational design of new environmentally safe insecticides.

Conformational transitions of human calcitonin (hCT) during fibril formation in the acidic and neutral conditions were investigated by high-resolution solid-state 13C NMR spectroscopy. In aqueous acetic acid solution (pH 3.3), a local α-helical form is present around Gly10, whereas a random coil form is dominant as viewed from Phe22, Ala26, and Ala31 in the monomer form on the basis of the 13C chemical shifts. On the other hand, a local β-sheet form as viewed from Gly10 and Phe22, and both β-sheet and random coil as viewed from Ala26 and Ala31 were detected in the fibril at pH 3.3. The results indicate that conformational transitions from α-helix to β-sheet, and from random coil to β-sheet forms occurred in the central and C-terminus regions, respectively, during the fibril formation. The increased 13C resonance intensities of fibrils after a certain delay time suggests that the fibrillation can be explained by a two-step reaction mechanism in which the first step is a homogeneous association to form a nucleus, and the second step is an autocatalytic heterogeneous fibrillation. In contrast to the fibril at pH 3.3, the fibril at pH 7.5 formed a local β-sheet conformation at the central region and exhibited a random coil at the C-terminus region. Not only a hydrophobic interaction among the amphiphilic α-helices, but also an electrostatic interaction between charged side chains can play an important role for the fibril formation at pH 7.5 and 3.3 acting as electrostatically favorable and unfavorable interactions, respectively. These results suggest that hCT fibrils are formed by stacking antiparallel β-sheets at pH 7.5 and a mixture of antiparallel and parallel β-sheets at pH 3.3.

Pigeon liver malic enzyme was inactivated and cleaved at Asp141, Asp194, and Asp464 by the Cu2+-ascorbate system in acidic environment. Site-specific mutagenesis was performed at these putative metal-binding sites. Three point mutants, D141N, D194N, and D464N; three double mutants, D(141,194)N, D(194,464)N, and D(141,464)N; and a triple mutant, D(141,194,464)N; as well as the wild-type malic enzyme (WT) were successfully cloned and expressed in Escherichia coli cells. All recombinant enzymes, except the triple mutant, were purified to apparent homogeneity by successive Q-Sepharose and adenosine-2′,5′-bisphosphate-agarose columns. The mutants showed similar apparent Km,NADP values to that of the WT. The Km,Mal value was increased in the D141N and D194N mutants. The Km,Mn value, on the other hand, was increased only in the D141N mutant by 14-fold, corresponding to ∼1.6 kcal/mol for the Asp141-Mn2+ binding energy. Substrate inhibition by l-malate was only observed in WT, D464N, and D(141,464)N. Initial velocity experiments were performed to derive the various kinetic parameters. The possible interactions between Asp141, Asp194, and Asp464 were analyzed by the double-mutation cycles and triple-mutation box. There are synergistic weakening interactions between Asp141 and Asp194 in the metal binding that impel the D(141,194)N double mutant to an overall specificity constant [kcat/(Kd,MnKm,MalKm,NADP)] at least four orders of magnitude smaller than the WT value. This difference corresponds to an increase of 6.38 kcal/mol energy barrier for the catalytic efficiency. Mutation at Asp464, on the other hand, has partial additivity on the mutations at Asp141 and Asp194. The overall specificity constants for the double mutants D(194,464)N and D(141,464)N or the triple mutant D(141,194,464)N were decreased by only 10- to 100-fold compared to the WT. These results strongly suggest the involvement of Asp141 in the Mn2+-l-malate binding for the pigeon liver malic enzyme. The Asp194 and Asp464, which may be oxidized by nonspecific binding of Cu2+, are involved in the Mn2+-l-malate binding or catalysis indirectly by modulating the binding affinity of Asp141 with the Mn2+.

The dihydrolipoamide succinyltransferase (E2o) component of the α-ketoglutarate dehydrogenase complex catalyzes the transfer of a succinyl group from the S-succinyldihydrolipoyl moiety to coenzyme A. E2o is normally a 24-mer, but is found as a trimer when E2o is expressed with a C-terminal [His]6 tag. The crystal structure of the trimeric form of the catalytic domain (CD) of the Escherichia coli E2o has been solved to 3.0 Å resolution using the Molecular Replacement method. The refined model contains an intact trimer in the asymmetric unit and has an R-factor of 0.257 (Rfree = 0.286) for 18,699 reflections between 10.0 and 3.0 Å resolution. The core of tE2oCD (residues 187–396) superimposes onto that of the cubic E2oCD with an RMS difference of 0.4 Å for all main-chain atoms. The C-terminal end of tE2oCD (residues 397–404) rotates by an average of 37° compared to cubic E2oCD, disrupting the normal twofold interface. Despite the alteration of quaternary structure, the active site of tE2oCD shows no significant differences from that of the cubic E2oCD, although several side chains in the active site are more ordered in the trimeric form of E2oCD. Analysis of the available sequence data suggests that the majority of E2 components have active sites that resemble that of E. coli E2oCD. The remaining E2 components can be divided into three groups based on active-site sequence similarity. Analysis of the surface properties of both crystal forms of E. coli E2oCD suggests key residues that may be involved in the protein–protein contacts that occur between the catalytic and lipoyl domains of E2o.

Calbindin D28k is a highly conserved Ca2+-binding protein abundant in brain and sensory neurons. The 261-residue protein contains six EF-hands packed into one globular domain. In this study, we have reconstituted calbindin D28k from two fragments containing three EF-hands each (residues 1–132 and 133–261, respectively), and from other combinations of small and large fragments. Complex formation is studied by ion-exchange and size-exclusion chromatography, electrophoresis, surface plasmon resonance, as well as circular dichroism (CD), fluorescence, and NMR spectroscopy. Similar chromatographic behavior to the native protein is observed for reconstituted complexes formed by mixing different sets of complementary fragments, produced by introducing a cut between EF-hands 1, 2, 3, or 4. The C-terminal half (residues 133–261) appears to have a lower intrinsic stability compared to the N-terminal half (residues 1–132). In the presence of Ca2+, NMR spectroscopy reveals a high degree of structural similarity between the intact protein and the protein reconstituted from the 1–132 and 133–261 fragments. The affinity between these two fragments is 2 × 107 M−1, with association and dissociation rate constants of 2.7 × 104 M−1 s−1 and 1.4 × 10−3 s−1, respectively. The complex formed in the presence of Ca2+ is remarkably stable towards unfolding by urea and heat. Both the complex and intact protein display cold and heat denaturation, although residual α-helical structure is seen in the urea denatured state at high temperature. In the absence of Ca2+, the fragments do not recombine to yield a complex resembling the intact apo protein. Thus, calbindin D28k is an example of a protein that can only be reconstituted in the presence of bound ligand. The α-helical CD signal is increased by 26% after addition of Ca2+ to each half of the protein. This suggests that Ca2+-induced folding of the fragments is important for successful reconstitution of calbindin D28k.

α-Macroglobulin inhibits a broad spectrum of proteinases by forming macromolecular cages inside which proteinases are cross-linked and trapped. Upon formation of a complex with proteinase, α-macroglobulin undergoes a large conformational change that results in the exposure of its receptor-binding domain (RBD). Engagement of this domain by α-macroglobulin receptor permits clearance of the α-macroglobulin: proteinase complex from circulation. The crystal structure of rat α1-macroglobulin RBD has been determined at 2.3 Å resolution. The RBD is composed of a nine-stranded β-sandwich and a single α-helix that has been implicated as part of the receptor binding site and that lies on the surface of the β-sandwich. The crystallographic asymmetric unit contains a dimer of RBDs related by approximate twofold symmetry such that the putative receptor recognition sites of the two monomers are contiguous. By gel filtration and ultracentrifugation, it is shown that RBD dimers form in solution with a dissociation constant of ∼50 μM. The structure of the RBD dimer might mimic a conformation of transformed α-macroglobulin in which the proposed receptor binding residues are exposed on one face of the dimer. A pair of phenylalanine residues replaces a cystine that is conserved in other members of the macroglobulin family. These residues participate in a network of aromatic side-chain interactions that appears to stabilize the dimer interface.

The contribution of the Ser45 hydrogen bond to biotin binding activation and equilibrium thermodynamics was investigated by biophysical and X-ray crystallographic studies. The S45A mutant exhibits a 1,700-fold greater dissociation rate and 907-fold lower equilibrium affinity for biotin relative to wild-type streptavidin at 37 °C, indicating a crucial role in binding energetics. The crystal structure of the biotin-bound mutant reveals only small changes from the wild-type bound structure, and the remaining hydrogen bonds to biotin retain approximately the same lengths. No additional water molecules are observed to replace the missing hydroxyl, in contrast to the previously studied D128A mutant. The equilibrium ΔG°, ΔH°, ΔS°, ΔC°P, and activation ΔG[Dagger] of S45A at 37 °C are −13.7 ± 0.1 kcal/mol, −21.1 ± 0.5 kcal/mol, −23.7 ± 1.8 cal/mol K, −223 ± 12 cal/mol K, and 20.0 ± 2.5 kcal/mol, respectively. Eyring analysis of the large temperature dependence of the S45A off-rate resolves the ΔH[Dagger] and ΔS[Dagger] of dissociation, 25.8 ± 1.2 kcal/mol and 18.7 ± 4.3 cal/mol K. The large increases of ΔH[Dagger] and ΔS[Dagger] in the mutant, relative to wild-type, indicate that Ser45 could form a hydrogen bond with biotin in the wild-type dissociation transition state, enthalpically stabilizing it, and constraining the transition state entropically. The postulated existence of a Ser45-mediated hydrogen bond in the wild-type streptavidin transition state is consistent with potential of mean force simulations of the dissociation pathway and with molecular dynamics simulations of biotin pullout, where Ser45 is seen to form a hydrogen bond with the ureido oxygen as biotin slips past this residue after breaking the native hydrogen bonds.

The chaperonin HSP60 (GroEL) proteins are essential in eubacterial genomes and in eukaryotic organelles. Functional regions inferred from mutation studies and the Escherichia coli GroEL 3D crystal complexes are evaluated in a multiple alignment across 43 diverse HSP60 sequences, centering on ATP/ADP and Mg2+ binding sites, on residues interacting with substrate, on GroES contact positions, on interface regions between monomers and domains, and on residues important in allosteric conformational changes. The most evolutionary conserved residues relate to the ATP/ADP and Mg2+ binding sites. Hydrophobic residues that contribute in substrate binding are also significantly conserved. A large number of charged residues line the central cavity of the GroEL–GroES complex in the substrate-releasing conformation. These span statistically significant intra- and inter-monomer three-dimensional (3D) charge clusters that are highly conserved among sequences and presumably play an important role interacting with the substrate. Unaligned short segments between blocks of alignment are generally exposed at the outside wall of the Anfinsen cage complex. The multiple alignment reveals regions of divergence common to specific evolutionary groups. For example, rickettsial sequences diverge in the ATP/ADP binding domain and Gram-positive sequences diverge in the allosteric transition domain. The evolutionary information of the multiple alignment proffers attractive sites for mutational studies.

The backbone dynamics of the four-helical bundle cytokine leukemia inhibitory factor (LIF) have been investigated using 15N NMR relaxation and amide proton exchange measurements on a murine–human chimera, MH35-LIF. For rapid backbone motions (on a time scale of 10 ps to 100 ns), as probed by 15N relaxation measurements, the dynamics parameters were calculated using the model-free formalism incorporating the model selection approach. The principal components of the inertia tensor of MH35-LIF, as calculated from its NMR structure, were 1:0.98:0.38. The global rotational motion of the molecule was, therefore, assumed to be axially symmetric in the analysis of its relaxation data. This yielded a diffusion anisotropy D∥/D⊥ of 1.31 and an effective correlation time (4D⊥ + 2D∥)−1 of 8.9 ns. The average values of the order parameters (S2) for the four helices, the long interhelical loops, and the N-terminus were 0.91, 0.84, and 0.65, respectively, indicating that LIF is fairly rigid in solution, except at the N-terminus. The S2 values for the long interhelical loops of MH35-LIF were higher than those of their counterparts in short-chain members of the four-helical bundle cytokine family. Residues involved in LIF receptor binding showed no consistent pattern of backbone mobilities, with S2 values ranging from 0.71 to 0.95, but residues contributing to receptor binding site III had relatively lower S2 values, implying higher amplitude motions than for the backbone of sites I and II. In the relatively slow motion regime, backbone amide exchange measurements showed that a number of amides from the helical bundle exchanged extremely slowly, persisting for several months in 2H2O at 37 °C. Evidence for local unfolding was considered, and correlations among various structure-related parameters and the backbone amide exchange rates were examined. Both sets of data concur in showing that LIF is one of the most rigid four-helical bundle cytokines.

Previously, we determined the DNA and amino acid sequences as well as biochemical and biophysical properties of a series of fungal phytases. The amino acid sequences displayed 49–68% identity between species, and the catalytic properties differed widely in terms of specific activity, substrate specificity, and pH optima. With the ultimate goal to combine the most favorable properties of all phytases in a single protein, we attempted, in the present investigation, to increase the specific activity of Aspergillus fumigatus phytase. The crystal structure of Aspergillus niger NRRL 3135 phytase known at 2.5 Å resolution served to specify all active site residues. A multiple amino acid sequence alignment was then used to identify nonconserved active site residues that might correlate with a given favorable property of interest. Using this approach, Gln27 of A. fumigatus phytase (amino acid numbering according to A. niger phytase) was identified as likely to be involved in substrate binding and/or release and, possibly, to be responsible for the considerably lower specific activity (26.5 vs. 196 U·[mg protein]−1 at pH 5.0) of A. fumigatus phytase when compared to Aspergillus terreus phytase, which has a Leu at the equivalent position. Site-directed mutagenesis of Gln27 of A. fumigatus phytase to Leu in fact increased the specific activity to 92.1 U·(mg protein)−1, and this and other mutations at position 27 yielded an interesting array of pH activity profiles and substrate specificities. Analysis of computer models of enzyme–substrate complexes suggested that Gln27 of wild-type A. fumigatus phytase forms a hydrogen bond with the 6-phosphate group of myo-inositol hexakisphosphate, which is weakened or lost with the amino acid substitutions tested. If this hydrogen bond were indeed responsible for the differences in specific activity, this would suggest product release as the rate-limiting step of the A. fumigatus wild-type phytase reaction.

The folding kinetics of human common-type acylphosphatase (cAcP) from its urea- and TFE-denatured states have been determined by stopped-flow fluorescence techniques. The refolding reaction from the highly unfolded state formed in urea is characterized by double exponential behavior that includes a slow phase associated with isomerism of the Gly53–Pro54 peptide bond. However, this slow phase is absent when refolding is initiated by dilution of the highly α-helical denatured state formed in the presence of 40% trifluoroethanol (TFE). NMR studies of a peptide fragment corresponding to residues Gly53–Gly69 of cAcP indicate that only the native-like trans isomer of the Gly–Pro peptide bond is significantly populated in the presence of TFE, whereas both the cis and trans isomers are found in an ∼1:9 ratio for the peptide bond in aqueous solution. Molecular modeling studies in conjunction with NMR experiments suggest that the trans isomer of the Gly53–Pro54 peptide bond is stabilized in TFE by the formation of a nonnative-like hydrogen bond between the CO group of Gly53 and the NH group of Lys57. These results therefore reveal that a specific nonnative interaction in the denatured state can increase significantly the overall efficiency of refolding.

The relationship between the structure of a free ligand in solution and the structure of its bound form in a complex is of great importance to the understanding of the energetics and mechanism of molecular recognition and complex formation. In this study, we use a structure-based thermodynamic approach to study the dissociation of the complex between the toxin microcystin-LR (MLR) and the catalytic domain of protein phosphatase-1 (PP-1c) for which the crystal structure of the complex is known. We have calculated the thermodynamic parameters (enthalpy, entropy, heat capacity, and free energy) for the dissociation of the complex from its X-ray structure and found the calculated dissociation constant (4.0 × 10−11) to be in excellent agreement with the reported inhibitory constant (3.9 × 10−11). We have also calculated the thermodynamic parameters for the dissociation of 47 PP-1c:MLR complexes generated by docking an ensemble of NMR solution structures of MLR onto the crystal structure of PP-1c. In general, we observe that the lower the root-mean-square deviation (RMSD) of the docked complex (compared to the X-ray complex) the closer its free energy of dissociation (ΔG°d) is to that calculated from the X-ray complex. On the other hand, we note a significant scatter between the ΔG°d and the RMSD of the docked complexes. We have identified a group of seven docked complexes with ΔG°d values very close to the one calculated from the X-ray complex but with significantly dissimilar structures. The analysis of the corresponding enthalpy and entropy of dissociation shows a compensation effect suggesting that MLR molecules with significant structural variability can bind PP-1c and that substantial conformational flexibility in the PP-1c:MLR complex may exist in solution.

An important goal of protein design is to understand the forces that stabilize a particular fold in preference to alternative folds. Here, we describe an extension of earlier studies in which we successfully designed a stable, native-like helical protein that is 50% identical in sequence to a predominantly β-sheet protein, the B1 domain of Streptococcal IgG-binding protein G. We report the characteristics of a series of variants of our original design that have even higher sequence identity to the B1 domain. Their properties illustrate the extent to which protein stability and conformation can be modulated through careful manipulation of key amino acid residues. Our results have implications for understanding conformational change phenomena of central biological importance and in probing the malleability of the sequence/structure relationship.

The amino terminal domain of enzyme I (residues 1–258 + Arg; EIN) and full length enzyme I (575 residues; EI) harboring active-site mutations (H189E, expected to have properties of phosphorylated forms, and H189A) have been produced by protein bioengineering. Differential scanning calorimetry (DSC) and temperature-induced changes in ellipticity at 222 nm for monomeric wild-type and mutant EIN proteins indicate two-state unfolding. For EIN proteins in 10 mM K-phosphate (and 100 mM KCl) at pH 7.5, ΔH ≅ 140 ± 10 (160) kcal mol−1 and ΔCp ≅ 2.7 (3.3) kcal K−1 mol−1. Transition temperatures (Tm) are 57 (59), 55 (58), and 53 (56) °C for wild-type, H189A, and H189E forms of EIN, respectively. The order of conformational stability for dephospho-His189, phospho-His189, and H189 substitutions of EIN at pH 7.5 is: His > Ala > Glu > His-PO32− due to differences in conformational entropy. Although H189E mutants have decreased Tm values for overall unfolding the amino terminal domain, a small segment of structure (3 to 12%) is stabilized (Tm ∼ 66–68 °C). This possibly arises from an ion pair interaction between the γ-carboxyl of Glu189 and the ε-amino group of Lys69 in the docking region for the histidine-containing phosphocarrier protein HPr. However, the binding of HPr to wild-type and active-site mutants of EIN and EI is temperature-independent (entropically controlled) with about the same affinity constant at pH 7.5: K′A = 3 ± 1 × 105 M−1 for EIN and ∼1.2 × 105 M−1 for EI.

Residue Asn47 at position L1 of a type II′ β-turn of the α-spectrin SH3 domain is located in a disallowed region of the Ramachandran plot (φ = 56 ± 12, ψ = −118 ± 17). Therefore, it is expected that replacement of Asn47 by Gly should result in a considerable stabilization of the protein. Thermodynamic analysis of the N47G and N47A mutants shows that the change in free energy is small (∼0.7 kcal/mol; ∼3 kJ/mol) and comparable to that found when mutating a Gly to Ala in a α-helix or β-sheet. X-ray structural analysis of these mutants shows that the conformation of the β-turn does not change upon mutation and, therefore, that there is no relaxation of the structure, nor is there any gain or loss of interactions that could explain the small energy change. Our results indicate that the energetic definition of II′ region of the Ramachandran plot (φ = 60 ± 30, ψ = −115 ± 15) should be revised for at least Ala and Asn in structure validation and protein design.