Core-packing constraints, hydrophobicity and protein design

Structural studies of intrinsically disordered MLL-fusion protein AF9 in complex with peptidomimetic inhibitors

AF9 (MLLT3) and its paralog ENL(MLLT1) are members of the YEATS family of proteins with important role in transcriptional and epigenetic regulatory complexes. These proteins are two common MLL fusion partners in MLL-rearranged leukemias. The oncofusion proteins MLL-AF9/ENL recruit multiple binding partners, including the histone methyltransferase DOT1L, leading to aberrant transcriptional activation and enhancing the expression of a characteristic set of genes that drive leukemogenesis. The interaction between AF9 and DOT1L is mediated by an intrinsically disordered C-terminal ANC1 homology domain (AHD) in AF9, which undergoes folding upon binding of DOT1L and other partner proteins. We have recently reported peptidomimetics that disrupt the recruitment of DOT1L by AF9 and ENL, providing a proof-of-concept for targeting AHD and assessing its druggability. Intrinsically disordered proteins, such as AF9 AHD, are difficult to study and characterize experimentally on a structural level. In this study, we present a successful protein engineering strategy to facilitate structural investigation of the intrinsically disordered AF9 AHD domain in complex with peptidomimetic inhibitors by using maltose binding protein (MBP) as a crystallization chaperone connected with linkers of varying flexibility and length. The strategic incorporation of disulfide bonds provided diffraction-quality crystals of the two disulfide-bridged MBP-AF9 AHD fusion proteins in complex with the peptidomimetics. These successfully determined first series of 2.1-2.6 Å crystal complex structures provide high-resolution insights into the interactions between AHD and its inhibitors, shedding light on the role of AHD in recruiting various binding partner proteins. We show that the overall complex structures closely resemble the reported NMR structure of AF9 AHD/DOT1L with notable difference in the conformation of the β-hairpin region, stabilized through conserved hydrogen bonds network. These first series of AF9 AHD/peptidomimetics complex structures are providing insights of the protein-inhibitor interactions and will facilitate further development of novel inhibitors targeting the AF9/ENL AHD domain.

Determination of the optimal connector length to enhance stability of backbone-circularized granulocyte colony-stimulating factor

Improving protein stability is important for industrial applications, and one promising method for achieving this is backbone circularization. As connector length affects stability, predicting and elucidating a more stable connector length is necessary for development of the backbone circularization method. However, the relationship between connector length and protein stability has not been completely elucidated. Here, we determined the most stable connector length for granulocyte colony-stimulating factor by changing one residue at a time to produce connector length variants and then measuring their thermal stability. Analysis of the local structures obtained from the predicted structures of the circularized variants revealed that an approach using helix length, dihedral backbone angle, and number of unbonded hydrogen bond donors and acceptors is suitable for identifying connector lengths with higher stability.

Modularity of the hydrophobic core and evolution of functional diversity in fold A glycosyltransferases

Hydrophobic cores are fundamental structural properties of proteins typically associated with protein folding and stability; however, how the hydrophobic core shapes protein evolution and function is poorly understood. Here, we investigated the role of conserved hydrophobic cores in fold-A glycosyltransferases (GT-As), a large superfamily of enzymes that catalyze formation of glycosidic linkages between diverse donor and acceptor substrates through distinct catalytic mechanisms (inverting versus retaining). Using hidden Markov models and protein structural alignments, we identify similarities in the phosphate-binding cassette (PBC) of GT-As and unrelated nucleotide-binding proteins, such as UDP-sugar pyrophosphorylases. We demonstrate that GT-As have diverged from other nucleotide-binding proteins through structural elaboration of the PBC and its unique hydrophobic tethering to the F-helix, which harbors the catalytic base (xED-Asp). While the hydrophobic tethering is conserved across diverse GT-A fold enzymes, some families, such as B3GNT2, display variations in tethering interactions and core packing. We evaluated the structural and functional impact of these core variations through experimental mutational analysis and molecular dynamics simulations and find that some of the core mutations (T336I in B3GNT2) increase catalytic efficiency by modulating the conformational occupancy of the catalytic base between “D-in” and acceptor-accessible “D-out” conformation. Taken together, our studies support a model of evolution in which the GT-A core evolved progressively through elaboration upon an ancient PBC found in diverse nucleotide-binding proteins, and malleability of this core provided the structural framework for evolving new catalytic and substrate-binding functions in extant GT-A fold enzymes.

Free energy simulations to understand the effect of Met → Ala mutations at positions 205, 206 and 213 on stability of human prion protein

Prion diseases are a family of infectious amyloid diseases affecting human and animals. Prion propagation in transmissible spongiform encephalopathies is associated with the unfolding and conversion of normal cellular prion protein into its pathogenic scrapie form. Understanding the fundamentals of prion protein aggregation caused by mutations is crucial to unravel the pathology of prion diseases. To help understand the contributions of individual residues to the stability of the human prion protein, we have carried out free energy simulations based on atomistic molecular dynamics trajectories. We focus on Met → Ala mutations at positions 205, 206 and 213, which are mostly buried residues located on helix 3 of the protein. The simulations predicted that all three mutations destabilize the prion protein. Changes in unfolding free energy upon mutation, ∆∆G, are 3.10 ± 0.79, 2.00 ± 0.26 and 3.06 ± 0.66 kcal/mol for M205A, M206A and M213A, respectively, in excellent agreement with the corresponding experimental values of 3.09 ± 0.28, 1.50 ± 0.34 and 3.12 ± 0.27 kcal/mol [T. Hart et al. (2009) PNAS 106, 5651-5656]. Component analysis indicates that the major contributions to the loss of protein stability arise from van der Waals interactions for the M205A and M206A mutations, and from van der Waals and covalent energy terms for M213A. Interestingly, while free energy contributions from a majority of residues neighboring the mutation sites tend to stabilize the wild type, there are a few residues stabilizing the mutant side chains. Our results show that this approach to free energy calculation can be very useful for understanding the detailed mechanism of human prion protein stability.

Funneled energy landscape unifies principles of protein binding and evolution

Most proteins have evolved to spontaneously fold into native structure and specifically bind with their partners for the purpose of fulfilling biological functions. According to Darwin, protein sequences evolve through random mutations, and only the fittest survives. The understanding of how the evolutionary selection sculpts the interaction patterns for both biomolecular folding and binding is still challenging. In this study, we incorporated the constraint of functional binding into the selection fitness based on the principle of minimal frustration for the underlying biomolecular interactions. Thermodynamic stability and kinetic accessibility were derived and quantified from a global funneled energy landscape that satisfies the requirements of both the folding into the stable structure and binding with the specific partner. The evolution proceeds via a bowl-like evolution energy landscape in the sequence space with a closed-ring attractor at the bottom. The sequence space is increasingly reduced until this ring attractor is reached. The molecular-interaction patterns responsible for folding and binding are identified from the evolved sequences, respectively. The residual positions participating in the interactions responsible for folding are highly conserved and maintain the hydrophobic core under additional evolutionary constraints of functional binding. The positions responsible for binding constitute a distributed network via coupling conservations that determine the specificity of binding with the partner. This work unifies the principles of protein binding and evolution under minimal frustration and sheds light on the evolutionary design of proteins for functions.

Stabilization of backbone-circularized protein is attained by synergistic gains in enthalpy of folded structure and entropy of unfolded structure

Backbone circularization is an effective technique for protein stabilization. Here, we investigated the effect of a connector, an engineered segment that connects two protein termini, on the conformational stability of previously designed circularized variants of granulocyte colony-stimulating factor (G-CSF). Heat tolerance and chemical denaturation analyses revealed that aggregation resistance and thermodynamic stability of the circularized variants were superior to those of linear G-CSF. Crystal structure and molecular dynamics (MD) simulation of the most thermodynamically stable variant (C166) revealed a high number of intramolecular hydrogen bonds in both the connector region and Helix D adjacent to the connector region in the folded structure. MD simulations and theoretical calculations involving different force fields indicated a reduction in the main chain entropy of C166 in the unfolded state and increase in the intramolecular hydrogen bond energy of C166 in the folded structure. Although backbone circularization is usually considered to alter chain entropy of the unfolded state, the data indicated that it could also improve the conformational enthalpy of the folded state. Further structural examination of the connector region confirmed that protein design based on a statistical analysis of local structures is an effective approach for predicting an optimum connector length to improve the conformational stability of backbone-circularized proteins. Protein design using backbone circularization with an optimum connector length will be useful for the development of effective and safe protein therapeutics. DATABASE: Structural data are available in Protein Data Bank under the accession number 5ZO6.

Structural Perturbations due to Mutation (H1047R) in Phosphoinositide-3-kinase (PI3Kα) and Its Involvement in Oncogenesis: An in Silico Insight

PI3Kα is a heterodimer protein consisting of two subunits (p110α and p85α) which promotes various signaling pathways. Oncogenic mutation in the catalytic subunit p110α of PI3Kα at the 1047 position in the kinase domain substitutes the histidine with arginine. This mutation brings about conformational transitions in the protein complex. These transitions expose the membrane binding region of PI3Kα, and then it independently binds to the cell membrane through its kinase domain without the involvement of the membrane-bound protein RAS. We observed notable changes between the protein complexes (p110α-p85α) of native and mutant structures at the atomic level using molecular dynamics simulations. Simulation results revealed formation of a less number of hydrogen bonds between the two subunits in the mutant protein complex which led the two subunits to move away from each other. This increase in distance between the subunits led to an expanded structure, thereby increasing the flexibility of the protein complex. Furthermore, a study of secondary structure elements and the electrostatic potential of the protein also gave a molecular insight into the change in interaction patterns of the protein with the plasma membrane. Our finding clearly indicates the role of mutation in oncogenesis and provides an insight into considering the structural aspects to handle this mutation.

Methionine in a protein hydrophobic core drives tight interactions required for assembly of spider silk

Web spiders connect silk proteins, so-called spidroins, into fibers of extraordinary toughness. The spidroin N-terminal domain (NTD) plays a pivotal role in this process: it polymerizes spidroins through a complex mechanism of dimerization. Here we analyze sequences of spidroin NTDs and find an unusually high content of the amino acid methionine. We simultaneously mutate all methionines present in the hydrophobic core of a spidroin NTD from a nursery web spider’s dragline silk to leucine. The mutated NTD is strongly stabilized and folds at the theoretical speed limit. The structure of the mutant is preserved, yet its ability to dimerize is substantially impaired. We find that side chains of core methionines serve to mobilize the fold, which can thereby access various conformations and adapt the association interface for tight binding. Methionine in a hydrophobic core equips a protein with the capacity to dynamically change shape and thus to optimize its function.

Spider silk is of interest in material science research. Here the authors show that the tight binding of a spider silk protein domain relies on the amino acid methionine, which is abundant in the domain core where it facilitates dynamic shape adaption of the binding interface.

Superfunneled Energy Landscape of Protein Evolution Unifies the Principles of Protein Evolution, Folding, and Design

Evolution is essential for shaping the biological functions. Darwin proposed the selection as the driving force for evolution upon mutations. While mutations are clear, the quantification of the selection force is still challenging. In this study, we identified and quantified both thermodynamic stability and kinetic accessibility as the selection forces for protein evolution. The protein evolution can be viewed and quantified as a trajectory moving along a superfunneled energy landscape with a line attractor at the bottom. The resulting evolved sequences and structures show strong protein characteristics including the hydrophobic core, high designability, and fast folding. The evolution principle uncovered here is validated on real proteins and sheds light on the protein design.

Allosteric Modulation of Binding Specificity by Alternative Packing of Protein Cores

Hydrophobic cores are often viewed as tightly packed and rigid, but they do show some plasticity and could thus be attractive targets for protein design. Here we explored the role of different functional pressures on the core packing and ligand recognition of the SH3 domain from human Fyn tyrosine kinase. We randomized the hydrophobic core and used phage display to select variants that bound to each of three distinct ligands. The three evolved groups showed remarkable differences in core composition, illustrating the effect of different selective pressures on the core. Changes in the core did not significantly alter protein stability, but were linked closely to changes in binding affinity and specificity. Structural analysis and molecular dynamics simulations revealed the structural basis for altered specificity. The evolved domains had significantly reduced core volumes, which in turn induced increased backbone flexibility. These motions were propagated from the core to the binding surface and induced significant conformational changes. These results show that alternative core packing and consequent allosteric modulation of binding interfaces could be used to engineer proteins with novel functions.

pH Dependence of Chitosan Enzymolysis

As a means of making chitosan more useful in biotechnological applications, it was hydrolyzed using pepsin, chitosanase and α-amylase. The enzymolysis behavior of these enzymes was further systematically studied for its effectiveness in the production of low-molecular-weight chitosans (LMWCs) and other derivatives. The study showed that these enzymes depend on ion hydronium (H3O+), thus on pH with a pH dependence fitting R2 value of 0.99. In y = 1.484[H^+] + 0.114, the equation of pH dependence, when [H^+] increases by one, y (k_0/k_m) increases by 1.484. From the temperature dependence study, the activation energy (Ea) and pre-exponential factor (A) were almost identical for two of the enzymes, but a considerable difference was observed in comparison with the third enzyme. Chitosanase and pepsin had nearly identical Ea, but α-amylase was significantly lower. This serves as evidence that the hydrolysis reaction of α-amylase relies on low-barrier hydrogen bonds (LBHBs), which explains its low Ea in actual conditions. The confirmation of this phenomenon was further derived from a similarly considerable difference in the order magnitudes of A between α-amylase and the other two enzymes, which was more than five. Variation of the rate constants of the enzymatic hydrolysis of chitosan with temperature follows the Arrhenius equation.

Simple approach for ranking structure determining residues

Mutating residues has been a common task in order to study structural properties of the protein of interest. Here, we propose and validate a simple method that allows the identification of structural determinants; i.e., residues essential for preservation of the stability of global structure, regardless of the protein topology. This method evaluates all of the residues in a 3D structure of a given globular protein by ranking them according to their connectivity and movement restrictions without topology constraints. Our results matched up with sequence-based predictors that look up for intrinsically disordered segments, suggesting that protein disorder can also be described with the proposed methodology.

Automated selection of stabilizing mutations in designed and natural proteins

The ability to engineer novel protein folds, conformations, and enzymatic activities offers enormous potential for the development of new protein therapeutics and biocatalysts. However, many de novo and redesigned proteins exhibit poor hydrophobic packing in their predicted structures, leading to instability or insolubility. The general utility of rational, structure-based design would greatly benefit from an improved ability to generate well-packed conformations. Here we present an automated protocol within the RosettaDesign framework that can identify and improve poorly packed protein cores by selecting a series of stabilizing point mutations. We apply our method to previously characterized designed proteins that exhibited a decrease in stability after a full computational redesign. We further demonstrate the ability of our method to improve the thermostability of a well-behaved native protein. In each instance, biophysical characterization reveals that we were able to stabilize the original proteins against chemical and thermal denaturation. We believe our method will be a valuable tool for both improving upon designed proteins and conferring increased stability upon native proteins.

Conservation and divergence on plant seed 11S globulins based on crystal structures

The crystal structures of two pro-11S globulins namely: rapeseed procruciferin and pea prolegumin are presented here. We have extensively compared them with the other known structures of plant seed 11S and 7S globulins. In general, the disordered regions in the crystal structures among the 11S globulins correspond to their five variable regions. Variable region III of procruciferin is relatively short and is in a loop conformation. This region is highly disordered in other pro-11S globulin crystals. Local helical and strand variations also occur across the group despite general structure conservation. We showed how these variations may alter specific physicochemical, functional and physiological properties. Aliphatic hydrophobic residues on the molecular surface correlate well with Tm values of the globulins. We also considered other structural features that were reported to influence thermal stability but no definite conclusion was drawn since each factor has additive or subtractive effect. Comparison between proA3B4 and mature A3B4 revealed an increase in r.m.s.d. values near variable regions II and IV. Both regions are on the IE face. Secondary structure based alignment of 11S and 7S globulins revealed 16 identical residues. Based on proA3B4 sequence, Pro60, Gly128, Phe163, Phe208, Leu213, Leu227, Ile237, Pro382, Val404, Pro425 and Val 466 are involved in trimer formation and stabilization. Gly28, Gly74, Asp135, Gly349 and Gly397 are involved in correct globular folding.

Orientation-dependent backbone-only residue pair scoring functions for fixed backbone protein design

Background

Empirical scoring functions have proven useful in protein structure modeling. Most such scoring functions depend on protein side chain conformations. However, backbone-only scoring functions do not require computationally intensive structure optimization and so are well suited to protein design, which requires fast score evaluation. Furthermore, scoring functions that account for the distinctive relative position and orientation preferences of residue pairs are expected to be more accurate than those that depend only on the separation distance.

Results

Residue pair scoring functions for fixed backbone protein design were derived using only backbone geometry. Unlike previous studies that used spherical harmonics to fit 2D angular distributions, Gaussian Mixture Models were used to fit the full 3D (position only) and 6D (position and orientation) distributions of residue pairs. The performance of the 1D (residue separation only), 3D, and 6D scoring functions were compared by their ability to identify correct threading solutions for a non-redundant benchmark set of protein backbone structures. The threading accuracy was found to steadily increase with increasing dimension, with the 6D scoring function achieving the highest accuracy. Furthermore, the 3D and 6D scoring functions were shown to outperform side chain-dependent empirical potentials from three other studies. Next, two computational methods that take advantage of the speed and pairwise form of these new backbone-only scoring functions were investigated. The first is a procedure that exploits available sequence data by averaging scores over threading solutions for homologs. This was evaluated by applying it to the challenging problem of identifying interacting transmembrane alpha-helices and found to further improve prediction accuracy. The second is a protein design method for determining the optimal sequence for a backbone structure by applying Belief Propagation optimization using the 6D scoring functions. The sensitivity of this method to backbone structure perturbations was compared with that of fixed-backbone all-atom modeling by determining the similarities between optimal sequences for two different backbone structures within the same protein family. The results showed that the design method using 6D scoring functions was more robust to small variations in backbone structure than the all-atom design method.

Conclusions

Backbone-only residue pair scoring functions that account for all six relative degrees of freedom are the most accurate and including the scores of homologs further improves the accuracy in threading applications. The 6D scoring function outperformed several side chain-dependent potentials while avoiding time-consuming and error prone side chain structure prediction. These scoring functions are particularly useful as an initial filter in protein design problems before applying all-atom modeling.

Protein stabilization by the rational design of surface charge-charge interactions

The design of proteins with increased stability has many important applications in biotechnology. In recent years, strategies involving directed evolution, sequence-based design, or computational design have proven successful for generating stabilized proteins. A brief overview of the various methods that have been used to increase protein stability is presented, followed by a detailed example of how the rational design of surface charge-charge interactions has provided a robust method for protein stabilization.

A computational approach for the rational design of stable proteins and enzymes: optimization of surface charge-charge interactions

The design of stable proteins and enzymes is not only of particular biotechnological importance, but also addresses some important fundamental questions. While there are a number of different options available for designing or engineering stable proteins, the field of computational design provides fast and universal methods for stabilizing proteins of interest. One of the successful computational design strategies focuses on stabilizing proteins through the optimization of charge-charge interactions on the protein surface. By optimizing surface interactions, it is possible to alleviate some of the challenges that accompany efforts to redesign the protein core. The rational design of surface charge-charge interactions also allows one to optimize only the interactions that are distant from binding sites or active sites, making it possible to increase stability without adversely affecting activity. The optimization of surface charge-charge interactions is discussed in detail along with the experimental evidence to demonstrate that this is a robust and universal approach to designing proteins with enhanced stability.

Recognition of misfolded proteins by Lon, a AAA(+) protease

Proteins unfold constantly in cells, especially under stress conditions. Degradation of denatured polypeptides by Lon and related ATP-dependent AAA(+) proteases helps prevent toxic aggregates formation and other deleterious consequences, but how these destructive enzymatic machines distinguish between damaged and properly folded proteins is poorly understood. Here, we show that Escherichia coli Lon recognizes specific sequences -- rich in aromatic residues -- that are accessible in unfolded polypeptides but hidden in most native structures. Denatured polypeptides lacking such sequences are poor substrates. Lon also unfolds and degrades stably folded proteins with accessible recognition tags. Thus, protein architecture and the positioning of appropriate targeting sequences allow Lon degradation to be dependent or independent of the folding status of a protein. Our results suggest that Lon can recognize multiple signals in unfolded polypeptides synergistically, resulting in nanomolar binding and a mechanism for discriminating irreversibly damaged proteins from transiently unfolded elements of structure.

Spatially directed assembly of a heterotetrameric Cre-Lox synapse restricts recombination specificity

The pseudo-fourfold homotetrameric synapse formed by Cre protein and target DNA restricts site-specific recombination to sequences containing dyad-symmetric Cre-binding repeats. Mixtures of engineered altered-specificity Cre monomers can form heterotetramers that recombine nonidentical asymmetric sequences, allowing greater flexibility for target site selection in the genome of interest. However, the variety of tetramers allowed by random subunit association increases the chances of unintended reactivity at nontarget sites. This problem can be circumvented by specifying a unique spatial arrangement of heterotetramer subunits. By reconfiguring intersubunit protein-protein contacts, we directed the assembly of two different Cre monomers, each having a distinct DNA sequence specificity, in an alternating (ABAB) configuration. This designed heterotetramer preferentially recombined a particular pair of asymmetric Lox sites over other pairs, whereas a mixture of freely associating subunits showed little bias. Alone, the engineered monomers had reduced reactivity towards both dyad-symmetric and asymmetric sites. Specificity arose because the organization of Cre-binding repeats of the preferred substrate matched the programmed arrangement of the subunits in the heterotetrameric synapse. When this "spatial matching" principle is applied, Cre-mediated recombination can be directed to asymmetric DNA sequences with greater fidelity.

Recent advances in our understanding of protein conformational stability from a pharmaceutical perspective

The marginal conformational stability of proteins has made them in some cases less than ideal candidates for pharmaceutical agents. Recent progress in our understanding of protein structure and stability has provided the opportunity to design the desired degree of stability into protein drug candidates. Modifications such as the optimisation of interior side-chain packing, the introduction of new ion-pairs, as well as the design of stabilising disulfide bridges and ligand binding sites, all offer the opportunity to produce proteins with enhanced stability properties.

Towards a three-alpha-helix bundle protein that binds volatile general anesthetics

The general anesthetics halothane and chloroform are capable of binding to synthetic water-soluble four-alpha-helix bundles, which model the putative in vivo receptors. In this study, we investigate the binding of these anesthetics to synthetic water-soluble three-alpha-helix bundles. A series of variants containing up to four X-to-Ala and up to four X-to-Met substitutions was made; and the effect of these substitutions on structure, stability and anesthetic binding affinity was examined. Generally, the amount of alpha-helix and the stability of the three-alpha-helix bundles decreased as the number of X-to-Ala substitutions increased. A concomitant red-shift in tryptophan fluorescence lambdamax was seen, suggesting an increased flexibility of the native structure. Up to four X-to-Met substitutions had little effect on the amount of alpha-helix, but an increase in tryptophan lambdamax was seen for the variants with three and four methionine substitutions. The exceptions were a) a variant with a clustering of alanine and methionine residues at one end of the three-alpha-helix bundle, suggesting a gate structure that can admit ligand molecules; and b) a variant with a single Leu35Ala substitution, suggesting that at select positions, the size of the side chain is important for defining anesthetic binding affinity.

Replacement of Staphylococcal Nuclease Hydrophobic Core Residues with Those from Thermophilic Homologues Indicates Packing is Improved in Some Thermostable Proteins

The importance of tight hydrophobic core packing in stabilizing proteins found in thermophilic organisms has been vigorously disputed. Here, portions of the cores found in three thermophilic homologues were transplanted into the core of staphylococcal nuclease, a protein of modest stability. Packing of the core was evaluated by comparing interaction energy of the three mutants to the comprehensive mutant library built up previously at these same sites in staphylococcal nuclease. It was found that the interaction energy of one thermophilic sequence is extraordinarily favorable and the interaction energies of other two transplanted thermophilic sequences are good, comparable to the interaction energies of mutant cores based on cores found in mesophilic homologues. As expected when transferring just a portion of the core sequence, the mutant proteins were destabilized overall relative to wild-type staphylococcal nuclease. The overall conclusion is that improvement of packing interactions is a mechanism to confer stability employed in some proteins from thermophiles, but not all.

Selection of stably folded proteins by phage-display with proteolysis

To facilitate the process of protein design and learn the basic rules that control the structure and stability of proteins, combinatorial methods have been developed to select or screen proteins with desired properties from libraries of mutants. One such method uses phage-display and proteolysis to select stably folded proteins. This method does not rely on specific properties of proteins for selection. Therefore, in principle it can be applied to any protein. Since its first demonstration in 1998, the method has been used to create hyperthermophilic proteins, to evolve novel folded domains from a library generated by combinatorial shuffling of polypeptide segments and to convert a partially unfolded structure to a fully folded protein.

Combinatorial approaches to protein stability and structure

Why do proteins adopt the conformations that they do, and what determines their stabilities? While we have come to some understanding of the forces that underlie protein architecture, a precise, predictive, physicochemical explanation is still elusive. Two obstacles to addressing these questions are the unfathomable vastness of protein sequence space, and the difficulty in making direct physical measurements on large numbers of protein variants. Here, we review combinatorial methods that have been applied to problems in protein biophysics over the last 15 years. The effects of hydrophobic core composition, the most important determinant of structure and stability, are still poorly understood. Particular attention is given to core composition as addressed by library methods. Increasingly useful screens and selections, in combination with modern high-throughput approaches borrowed from genomics and proteomics efforts, are making the empirical, statistical correlation between sequence and structure a tractable problem for the coming years.

Crystallographic analysis of an "anticalin" with tailored specificity for fluorescein reveals high structural plasticity of the lipocalin loop region

The artificial lipocalin FluA with novel specificity toward fluorescein was derived via combinatorial engineering from the bilin-binding protein, BBP by exchange of 16 amino acids in the ligand pocket. Here, we describe the crystal structure of FluA at 2.0 A resolution in the space group P2(1) with two protein-ligand complexes in the asymmetric unit. In both molecules, the characteristic beta-barrel architecture with the attached alpha-helix is well preserved. In contrast, the four loops at one end of the beta-barrel that form the entrance to the binding site exhibit large conformational deviations from the wild-type protein, which can be attributed to the sidechain replacements. Specificity for the new ligand is furnished by hydrophobic packing, charged sidechain environment, and hydrogen bonds with its hydroxyl groups. Unexpectedly, fluorescein is bound in a much deeper cavity than biliverdin IX(gamma) in the natural lipocalin. Triggered by the substituted residues, unmutated sidechains at the bottom of the binding site adopt conformations that are quite different from those observed in the BBP, illustrating that not only the loop region but also the hydrophobic interior of the beta-barrel can be reshaped for molecular recognition. Particularly, Trp 129 participates in a tight stacking interaction with the xanthenolone moiety, which may explain the ultrafast electron transfer that occurs on light excitation of the bound fluorescein. These structural findings support our concept of using lipocalins as a scaffold for the engineering of so-called "anticalins" directed against prescribed targets as an alternative to recombinant antibody fragments.

How do helix-helix interactions help determine the folds of membrane proteins? Perspectives from the study of homo-oligomeric helical bundles

The final, structure-determining step in the folding of membrane proteins involves the coalescence of preformed transmembrane helices to form the native tertiary structure. Here, we review recent studies on small peptide and protein systems that are providing quantitative data on the interactions that drive this process. Gel electrophoresis, analytical ultracentrifugation, and fluorescence resonance energy transfer (FRET) are useful methods for examining the assembly of homo-oligomeric transmembrane helical proteins. These methods have been used to study the assembly of the M2 proton channel from influenza A virus, glycophorin, phospholamban, and several designed membrane proteins-all of which have a single transmembrane helix that is sufficient for association into a transmembrane helical bundle. These systems are being studied to determine the relative thermodynamic contributions of van der Waals interactions, conformational entropy, and polar interactions in the stabilization of membrane proteins. Although the database of thermodynamic information is not yet large, a few generalities are beginning to emerge concerning the energetic differences between membrane and water-soluble proteins: the packing of apolar side chains in the interior of helical membrane proteins plays a smaller, but nevertheless significant, role in stabilizing their structure. Polar, hydrogen-bonded interactions occur less frequently, but, nevertheless, they often provide a strong driving force for folding helix-helix pairs in membrane proteins. These studies are laying the groundwork for the design of sequence motifs that dictate the association of membrane helices.

Preorganization of molecular binding sites in designed diiron proteins

De novo protein design provides an attractive approach to critically test the features that are required for metalloprotein structure and function. Previously we designed and crystallographically characterized an idealized dimeric model for the four-helix bundle class of diiron and dimanganese proteins [Dueferri 1 (DF1)]. Although the protein bound metal ions in the expected manner, access to its active site was blocked by large bulky hydrophobic residues. Subsequently, a substrate-access channel was introduced proximal to the metal-binding center, resulting in a protein with properties more closely resembling those of natural enzymes. Here we delineate the energetic and structural consequences associated with the introduction of these binding sites. To determine the extent to which the binding site was preorganized in the absence of metal ions, the apo structure of DF1 in solution was solved by NMR and compared with the crystal structure of the di-Zn(II) derivative. The overall fold of the apo protein was highly similar to that of the di-Zn(II) derivative, although there was a rotation of one of the helices. We also examined the thermodynamic consequences associated with building a small molecule-binding site within the protein. The protein exists in an equilibrium between folded dimers and unfolded monomers. DF1 is a highly stable protein (K(diss) = 0.001 fM), but the dissociation constant increases to 0.6 nM (deltadeltaG = 5.4 kcalmol monomer) as the active-site cavity is increased to accommodate small molecules.

Unusual fluorescence of W168 in Plasmodium falciparum triosephosphate isomerase, probed by single-tryptophan mutants

Plasmodium falciparum triosephosphate isomerase (PfTIM) contains two tryptophan residues, W11 and W168. One is positioned in the interior of the protein, and the other is located on the active-site loop 6. Two single-tryptophan mutants, W11F and W168F, were constructed to evaluate the contributions of each chromophore to the fluorescence of the wild-type (wt) protein and to probe the utility of the residues as spectroscopic reporters. A comparative analysis of the fluorescence spectra of PfTIMwt and the two mutant proteins revealed that W168 possesses an unusual, blue-shifted emission (321 nm) and exhibits significant red-edge excitation shift of fluorescence. In contrast, W11 emits at 332 nm, displays no excitation dependence of fluorescence, and behaves like a normal buried chromophore. W168 has a much shorter mean lifetime (2.7 ns) than W11 (4.6 ns). The anomalous fluorescence properties of W168 are abolished on unfolding of the protein in guanidinium chloride (GdmCl) or at low pH. Analysis of the tryptophan environment using a 1.1-A crystal structure established that W168 is rigidly held by a complex network of polar interactions including a strong hydrogen bond from Y164 to the indole NH group. The environment is almost completely polar, suggesting that electrostatic effects determine the unusually low emission wavelength of W168. To our knowledge this is a unique observation of a blue-shifted emission from a tryptophan in a polar environment in the protein. The wild-type and mutant proteins show similar levels of enzymatic activity and secondary and tertiary structure. However, the W11F mutation appreciably destabilizes the protein to unfolding by urea and GdmCl. The fluorescence of W168 is shown to be extremely sensitive to binding of the inhibitor, 2-phosphoglycolic acid.

Attenuating functions of the C terminus of lambda integrase

The tyrosine family site-specific recombinases, in contrast to the related type I topoisomerases, which act as monomers on a single DNA molecule, rely on multi-protein complexes to synapse partner DNAs and coordinate two sequential strand exchanges involving four nicking-closing reactions. Here, we analyze three mutants of the catalytic domain of lambda integrase (Int), A241V, I353M and W350ter that are defective for normal recombination, but possess increased topoisomerase activity. The mutant enzymes can carry out individual DNA strand exchanges using truncated substrates or Holliday junctions, and they show more DNA-cleavage activity than wild-type Int on isolated att sites. Structural modeling predicts that the substituted residues may destabilize interactions between the C-terminal beta-strand (beta7) of Int and the core of the protein. The cleavage-competent state of Int requires the repositioning of the nucleophile (Y342) located on beta6 and the catalyst K235 located on the flexible beta2-beta3 loop, relative to their positions in a crystal structure of the inactive conformation. We propose that the anchoring of beta7 against the protein core restrains the movement of Tyr342 and/or Lys235, causing an attenuation of cleavage activity in most contexts. Within a bona fide recombination complex, the release of strand beta7 would allow Tyr342 and Lys235 to assume catalytically active conformations in coordination with other Int protomers in the complex. The loss of beta7 packing by misalignment or truncation in the mutant proteins described here causes a loss of regulated activity, thereby favoring DNA cleavage activity in monomeric complexes and forfeiting the coordination of strand-exchange necessary for efficient recombination.

Stability scale and atomic solvation parameters extracted from 1023 mutation experiments

The stability scale of 20 amino acid residues is derived from a database of 1023 mutation experiments on 35 proteins. The resulting scale of hydrophobic residues has an excellent correlation with the octanol-to-water transfer free energy corrected with an additional Flory-Huggins molar-volume term (correlation coefficient r = 0.95, slope = 1.05, and a near zero intercept). Thus, hydrophobic contribution to folding stability is characterized remarkably well by transfer experiments. However, no corresponding correlation is found for hydrophilic residues. Both the hydrophilic portion and the entire scale, however, correlate strongly with average burial accessible surface (r = 0.76 and 0.97, respectively). Such a strong correlation leads to a near uniform value of the atomic solvation parameters for atoms C, S, O/N, O(-0.5), and N(+0.5,1). All are in the range of 12-28 cal x mol(-1) A(-2), close to the original estimate of hydrophobic contribution of 25-30 cal x mol(-1) A(-2) to folding stability. Without any adjustable parameters, the new stability scale and new atomic solvation parameters yielded an accurate prediction of protein-protein binding free energy for a separate database of 21 protein-protein complexes (r = 0.80 and slope = 1.06, and r = 0.83 and slope = 0.93, respectively).

Thermodynamic consequences of burial of polar and non-polar amino acid residues in the protein interior

Effects of amino acid substitutions at four fully buried sites of the ubiquitin molecule on the thermodynamic parameters (enthalpy, Gibbs energy) of unfolding were evaluated experimentally using differential scanning calorimetry. The same set of substitutions has been incorporated at each of four sites. These substitutions have been designed to perturb packing (van der Waals) interactions, hydration, and/or hydrogen bonding. From the analysis of the thermodynamic parameters for these ubiquitin variants we conclude that: (i) packing of non-polar groups in the protein interior is favorable and is largely defined by a favorable enthalpy of van der Waals interactions. The removal of one methylene group from the protein interior will destabilize a protein by approximately 5 kJ/mol, and will decrease the enthalpy of a protein by 12 kJ/mol. (ii) Burial of polar groups in the non-polar interior of a protein is highly destabilizing, and the degree of destabilization depends on the relative polarity of this group. For example, burial of Thr side-chain in the non-polar interior will be less destabilizing than burial of Asn side-chain. This decrease in stability is defined by a large enthalpy of dehydration of polar groups upon burial. (iii) The destabilizing effect of dehydration of polar groups upon burial can be compensated if these buried polar groups form hydrogen bonding. The enthalpy of this hydrogen bonding will compensate for the unfavorable dehydration energy and as a result the effect will be energetically neutral or even slightly stabilizing.

Computer-assisted re-design of spectrin SH3 residue clusters

We have developed a protein design computer program, called Perla, which performs searches in sequence space to uncover optimal amino acid sequences for desired protein three-dimensional structures. Optimal sequences are localised at the minima of a sequence-structure energy landscape defined using a complex scoring function (an all-atom molecular mechanics force field plus statistical terms including entropy and solvation) measured with respect to a reference state simulating a denatured protein. Sequence choices eventually optimise side chain packing, secondary structure propensities, and hydrogen bonding and electrostatics interactions. Perla was used to re-design clusters of residues of the SH3 domain of alpha-spectrin. Several mutant proteins were produced and characterised. Some of our designed proteins have significantly higher stabilities (stability enhancements about 0.25, 0.70 and 1.0 kcal mol(-1)) than the wild-type protein. These successful protein re-designs, and similar examples found in the literature, establish the quality of the structure-based computational approach to protein design.

Four-body potentials reveal protein-specific correlations to stability changes caused by hydrophobic core mutations

Mutational experiments show how changes in the hydrophobic cores of proteins affect their stabilities. Here, we estimate these effects computationally, using four-body likelihood potentials obtained by simplicial neighborhood analysis of protein packing (SNAPP). In this procedure, the volume of a known protein structure is tiled with tetrahedra having the center of mass of one amino acid side-chain at each vertex. Log-likelihoods are computed for the 8855 possible tetrahedra with equivalent compositions from structural databases and amino acid frequencies. The sum of these four-body potentials for tetrahedra present in a given protein yields the SNAPP score. Mutations change this sum by changing the compositions of tetrahedra containing the mutated residue and their related potentials. Linear correlation coefficients between experimental mutational stability changes, Delta(DeltaG(unfold)), and those based on SNAPP scoring range from 0.70 to 0.94 for hydrophobic core mutations in five different proteins. Accurate predictions for the effects of hydrophobic core mutations can therefore be obtained by virtual mutagenesis, based on changes to the total SNAPP likelihood potential. Significantly, slopes of the relation between Delta(DeltaG(unfold)) and DeltaSNAPP for different proteins are statistically distinct, and we show that these protein-specific effects can be estimated using the average SNAPP score per residue, which is readily derived from the analysis itself. This result enhances the predictive value of statistical potentials and supports previous suggestions that "comparable" mutations in different proteins may lead to different Delta(DeltaG(unfold)) values because of differences in their flexibility and/or conformational entropy.

De novo design and structural characterization of proteins and metalloproteins

De novo protein design has recently emerged as an attractive approach for studying the structure and function of proteins. This approach critically tests our understanding of the principles of protein folding; only in de novo design must one truly confront the issue of how to specify a protein's fold and function. If we truly understand proteins, it should be possible to design receptors, enzymes, and ion channels from scratch. Further, as this understanding evolves and is further refined, it should be possible to design proteins and biomimetic polymers with properties unprecedented in nature.

Contribution of proton linkage to the thermodynamic stability of the major cold-shock protein of Escherichia coli CspA

The stability of protein is defined not only by the hydrogen bonding, hydrophobic effect, van der Waals interactions, and salt bridges. Additional, much more subtle contributions to protein stability can arise from surface residues that change their properties upon unfolding. The recombinant major cold shock protein of Escherichia coli CspA an all-beta protein unfolds reversible in a two-state manner, and behaves in all other respects as typical globular protein. However, the enthalpy of CspA unfolding strongly depends on the pH and buffer composition. Detailed analysis of the unfolding enthalpies as a function of pH and buffers with different heats of ionization shows that CspA unfolding in the pH range 5.5-9.0 is linked to protonation of an amino group. This amino group appears to be the N-terminal alpha-amino group of the CspA molecule. It undergoes a 1.6 U shift in pKa values between native and unfolded states. Although this shift in pKa is expected to contribute approximately 5 kJ/mol to CspA stabilization energy the experimentally observed stabilization is only approximately 1 kJ/mol. This discrepancy is related to a strong enthalpy-entropy compensation due, most likely, to the differences in hydration of the protonated and deprotonated forms of the alpha-amino group.

De novo protein design. II. Plasticity in sequence space

It is generally accepted that many different protein sequences have similar folded structures, and that there is a relatively high probability that a new sequence possesses a previously observed fold. An indirect consequence of this is that protein design should define the sequence space accessible to a given structure, rather than providing a single optimized sequence. We have recently developed a new approach for protein sequence design, which optimizes the complete sequence of a protein based on the knowledge of its backbone structure, its amino acid composition and a physical energy function including van der Waals interactions, electrostatics, and environment free energy. The specificity of the designed sequence for its template backbone is imposed by keeping the amino acid composition fixed. Here, we show that our procedure converges in sequence space, albeit not to the native sequence of the protein. We observe that while polar residues are well conserved in our designed sequences, non-polar amino acids at the surface of a protein are often replaced by polar residues. The designed sequences provide a multiple alignment of sequences that all adopt the same three-dimensional fold. This alignment is used to derive a profile matrix for chicken triose phosphate isomerase, TIM. The matrix is found to recognize significantly the native sequence for TIM, as well as closely related sequences. Possible application of this approach to protein fold recognition is discussed.

Effects of core-packing on the structure, function, and mechanics of a four-helix-bundle protein ROP

The effects of core-packing on the structure, function and mechanics of the RNA-binding 4-helix-bundle Rop have been studied by molecular dynamics simulations. The structural, dynamical and geometrical properties of the Rop homodimer, (formed by the antiparallel juxtaposition of two helix-turn-helix motifs), have been compared with those of three protein variants described by Munson et al. (Protein Sci, 5:1584-1593, 1996), where the core of the native protein has been systematically repacked using a two-amino acid alphabet: Ala(2)Leu(2)-8, Ala(2)Leu(2)-8-rev, and Leu(2)Ala(2)-8. The results showed that it was possible to readily distinguish the inactive protein Leu(2)Ala(2)-8 from the other functionally active systems based on tertiary and quaternary structure criteria. Structural properties such as native secondary structure content did not correlate with biological activity. Biological activity was related in part to the relative arrangement of the residues within the binding site. But, more global aspects, related to the overall topology of the helical bundle, accounted for the small functional differences between Ala(2)Leu(2)-8 and Ala(2)Leu(2)-8-rev. Mechanically, the 4-helix-bundle absorbed core mutations by altering the local structure at the sequence termini and in the turns that join the two helices of each monomer, and by changing the overall orientation and separation of the extremely rigid helices. Proteins 1999;36:436-446.

Experimental evolution of a dense cluster of residues in tyrosyl-tRNA synthetase: quantitative effects on activity, stability and dimerization

A dense cluster of eight residues was identified at the crossing of two alpha-helices in tyrosyl-tRNA synthetase (TyrRS) from the thermophile Bacillus stearothermophilus. Its mechanism of evolution was characterized. Four residues of this cluster are not conserved in TyrRS from the mesophile Escherichia coli. The corresponding mutations were constructed in TyrRS(Delta1), a derivative of TyrRS from B. stearothermophilus in which the anticodon binding domain is deleted. Mutations I52L (i.e. Ile52 into Leu), M55L and L105V did not affect the activity of TyrRS(Delta1) in the pyrophosphate exchange reaction whereas T51P increased it. The kinetic stabilities of TyrRS(Delta1) and its mutant derivatives at 68.5 degreesC were determined from experiments of irreversible thermal precipitation. They were in the order L105V<I52L<T51P<Wild Type</=M55L; mutation I52L partially compensated L105V in these experiments whereas M55L was coupled neither to I52L nor to L105V. Mutations I52L and L105V affected the stability of the dimeric TyrRS(Delta1) at different steps of its unfolding by urea, monitored under equilibrium conditions by spectrofluorometry or size exclusion chromatography. I52L destabilized the association between the subunits even though residue Ile52 is more than 20 A away from the subunit interface. L105V destabilized the monomeric intermediate of unfolding. The two mutational pathways, going from the wild-type TyrRS(Delta1) to the I52L-L105V double mutant through each of the single mutants were not equivalent for the stability of the monomeric intermediate and for the total stability of the dimer. One pathway contained two neutral steps whereas the other pathway contained a destabilizing step followed by a stabilizing step. Mutation I52L allowed L105V along the first pathway and compensated it along the second pathway. Thus, the effects of I52L and L105V on stability depended on the structural context. The gain in activity due to T51P was at the expense of a slight destabilization.

Separation used for purification of recombinant proteins

The purification of molecules from recombinant cells may be strongly influenced by the molecular biology of gene isolation and expression. At the beginning of the process there may be a demand for information on the minute amounts of proteins and thus for ever increasingly sensitive techniques. Purification of recombinant proteins can differ from conventional purifications in several ways, depending on the solubility of the protein, occurrence in inclusion bodies, creation of fusion proteins with tags that enable simpler purification. Sometimes a (re)naturation step is required to get a bioactive protein. On the other hand, the techniques used in separation are essentially the same as for purification from the natural source and environment.

Ideal architecture of residue packing and its observation in protein structures

A simple model of sphere packing has been investigated as an ideal model for long-range interactions for the packing of non-bonded residues in protein structures. By superposing all residues, the geometry of packing around a central residue is investigated. It is found that all residues conform almost perfectly to this lattice model for sphere packing when a radius of 6.5 A is used to define non-bonded (virtual) interacting residues. Side-chain positions with respect to sequential backbone segments are relatively regular as well. This lattice can readily be used in conformation simulations to reduce the conformational space.

Protein folding: the endgame

The last stage of protein folding, the "endgame," involves the ordering of amino acid side-chains into a well defined and closely packed configuration. We review a number of topics related to this process. We first describe how the observed packing in protein crystal structures is measured. Such measurements show that the protein interior is packed exceptionally tightly, more so than the protein surface or surrounding solvent and even more efficiently than crystals of simple organic molecules. In vitro protein folding experiments also show that the protein is close-packed in solution and that the tight packing and intercalation of side-chains is a final and essential step in the folding pathway. These experimental observations, in turn, suggest that a folded protein structure can be described as a kind of three-dimensional jigsaw puzzle and that predicting side-chain packing is possible in the sense of solving this puzzle. The major difficulty that must be overcome in predicting side-chain packing is a combinatorial "explosion" in the number of possible configurations. There has been much recent progress towards overcoming this problem, and we survey a variety of the approaches. These approaches differ principally in whether they use ab initio (physical) or more knowledge-based methods, how they divide up and search conformational space, and how they evaluate candidate configurations (using scoring functions). The accuracy of side-chain prediction depends crucially on the (assumed) positioning of the main-chain. Methods for predicting main-chain conformation are, in a sense, not as developed as that for side-chains. We conclude by surveying these methods. As with side-chain prediction, there are a great variety of approaches, which differ in how they divide up and search space and in how they score candidate conformations.

The denaturation and degradation of stable enzymes at high temperatures

Now that enzymes are available that are stable above 100 degrees C it is possible to investigate conformational stability at this temperature, and also the effect of high-temperature degradative reactions in functioning enzymes and the inter-relationship between degradation and denaturation. The conformational stability of proteins depends upon stabilizing forces arising from a large number of weak interactions, which are opposed by an almost equally large destabilizing force due mostly to conformational entropy. The difference between these, the net free energy of stabilization, is relatively small, equivalent to a few interactions. The enhanced stability of very stable proteins can be achieved by an additional stabilizing force which is again equivalent to only a few stabilizing interactions. There is currently no strong evidence that any particular interaction (e.g. hydrogen bonds, hydrophobic interactions) plays a more important role in proteins that are stable at 100 degrees C than in those stable at 50 degrees C, or that the structures of very stable proteins are systematically different from those of less stable proteins. The major degradative mechanisms are deamidation of asparagine and glutamine, and succinamide formation at aspartate and glutamate leading to peptide bond hydrolysis. In addition to being temperature-dependent, these reactions are strongly dependent upon the conformational freedom of the susceptible amino acid residues. Evidence is accumulating which suggests that even at 100 degrees C deamidation and succinamide formation proceed slowly or not at all in conformationally intact (native) enzymes. Whether this is the case at higher temperatures is not yet clear, so it is not known whether denaturation of degradation will set the upper limit of stability for enzymes.

Access of ligands to cavities within the core of a protein is rapid

We have investigated the magnitude and timescale of fluctuations within the core of a protein using the exchange kinetics of indole and benzene binding to engineered hydrophobic cavities in T4 lysozyme. The crystal structures of variant-benzene complexes suggest that relatively large scale fluctuations (1-2 angstrom) of backbone atoms are required for entry of these ligands into the core. Nonetheless, these ligands enter the cavities rapidly, with bimolecular rate constants of approximately 10(6)-10(7) M(-1) s(-1) and a low activation barrier, 2-5 kcal mol(-1). These results suggest that protein cores undergo substantial fluctuations on the millisecond to microsecond timescale and that entry of small molecules into protein interiors is not strongly limited by steric occlusion.

Modeling side-chain conformation

Over the past few years, a number of methods for the calculation of side-chain conformations in proteins have been described. More recent studies have considered the effect of combinatorial packing, derivations from idealized rotameric structures and, to a limited extent, backbone flexibility on the quality and efficiency of calculations of protein side-chain conformation. Although further work is needed to address the issue of backbone displacements, the recent progress solves the packing problem to a significant degree. This opens the way for fruitful incorporation of these methods into general procedures for homology modeling and studies of ligand-protein interactions.

Exploring the allowed sequence space of a membrane protein

We present a comprehensive view of the tolerance of a membrane protein to sequence substitution. We find that the protein, diacylglycerol kinase from Escherichia coli, is extremely tolerant to sequence changes with three-quarters of the residues tolerating non-conservative changes. The conserved residues are distributed with approximately the same frequency in the soluble and transmembrane portions of the protein, but the most critical active-site residues appear to residue in the second cytoplasmic domain. It is remarkable that a unique structure of the membrane embedded portion of the protein can be encoded by a sequence that is so tolerant to substitution.

Conservative substitutions in the hydrophobic core of Rhodobacter sphaeroides thioredoxin produce distinct functional effects

The internal residue Phe 25 in Rhodobacter sphaeroides thioredoxin was changed to five amino acids (Ala, Val, Leu, Ile, Tyr) by site-directed mutagenesis, and the mutant proteins were characterized in vitro and in vivo using the mutant trxA genes in an Escherichia coli TrxA- background. The substitution F25A severely impaired the functional properties of the enzyme. Strains expressing all other mutations can grow on methionine sulfoxide with growth efficiencies of 45-60% that of the wild type at 37 degrees, and essentially identical at 42 degrees. At both temperatures, however, strains harboring the substitutions F25V and F25Y had lower growth rates and formed smaller colonies. In another in vivo assay, only the wild type and the F25I substitution allowed growth of phage T3/7 at 37 degrees, demonstrating that subtle modifications of the protein interior at position 25 Ile/Leu or Phe/Tyr) can produce significant biological effects. All F25 mutants were good substrates for E. coli thioredoxin reductase. Although turnover rates and apparent Km values were significantly lower for all mutants compared to the wild type, catalytic efficiency of thioredoxin reductase was similar for all substrates. Determination of the free energy of unfolding showed that the aliphatic substitutions (Val, Leu, Ile) significantly destabilized the protein, whereas the F25Y substitution did not affect protein stability. Thus, thermodynamic stability of R. sphaeroides thioredoxin variants is not correlated with the distinct functional effects observed both in vivo and in vitro.

New strategies in protein design

Initially, it was hoped that very simple rules could be sued to design proteins that embody all the characteristics of natural proteins. Indeed, with single-domain proteins as targets, it has been possible to design proteins that adopt the desired global fold. Yet, designed proteins with well defined structures and properties that mimic those of natural proteins remain elusive. Recent efforts in protein design have been directed toward addressing the basis for non-native characteristics in most protein designs. Although it is clear that specific tertiary interactions between all residues in a protein contribute to the final folded state, much attention has been placed on optimizing the packing of side chains in the hydrophobic core, with substantial success.

Evidence on close packing and cavities in proteins

The packing of a protein's constituent atoms and the attendant constraints placed upon them form the basis of many attempts to understand and predict protein structure, stability, folding and even function. Although the significance of packing is yet to be fully comprehended, recent experimental and theoretical investigations have increased our understanding through the description of mutational effects on structure and stability, determination of the limits of packing constraints for both protein folding and structure prediction, and delineation of packing guidelines on the basis of observed cavities in the native protein folds. These advances and allowing protein modellers, engineers and designers to tackle their problems from a more rational perspective.