A cavity-forming mutation in insulin induces segmental unfolding of a surrounding α-helix

To investigate the cooperativity of insulin's structure, a cavity-forming substitution was introduced within the hydrophobic core of an engineered monomer. The substitution, Ile(A2)-->Ala in the A1-A8 alpha-helix, does not impair disulfide pairing between chains. In accord with past studies of cavity-forming mutations in globular proteins, a decrement was observed in thermodynamic stability (DeltaDeltaG(u) 0.4-1.2 kcal/mole). Unexpectedly, CD studies indicate an attenuated alpha-helix content, which is assigned by NMR spectroscopy to selective destabilization of the A1-A8 segment. The analog's solution structure is otherwise similar to that of native insulin, including the B chain's supersecondary structure and a major portion of the hydrophobic core. Our results show that (1) a cavity-forming mutation in a globular protein can lead to segmental unfolding, (2) tertiary packing of Ile(A2), a residue of low helical propensity, stabilizes the A1-A8 alpha-helix, and (3) folding of this segment is not required for native disulfide pairing or overall structure. We discuss these results in relation to a hierarchical pathway of protein folding and misfolding. The Ala(A2) analog's low biological activity (0.5% relative to the parent monomer) highlights the importance of the A1-A8 alpha-helix in receptor recognition.     

Hierarchical protein "un-design": insulin's intrachain disulfide bridge tethers a recognition alpha-helix

A hierarchical pathway of protein folding can enable segmental unfolding by design. A monomeric insulin analogue containing pairwise substitution of internal A6-A11 cystine with serine [[Ser(A6),Ser(A11),Asp(B10),Lys(B28),Pro(B29)]insulin (DKP[A6-A11](Ser))] was previously investigated as a model of an oxidative protein-folding intermediate [Hua, Q. X., et al. (1996) J. Mol. Biol. 264, 390-403]. Its structure exhibits local unfolding of an adjoining amphipathic alpha-helix (residues A1-A8), leading to a 2000-fold reduction in activity. Such severe loss of function, unusual among mutant insulins, is proposed to reflect the cost of induced fit: receptor-directed restoration of the alpha-helix and its engagement in the hormone's hydrophobic core. To test this hypothesis, we have synthesized and characterized the corresponding alanine analogue [[Ala(A6),Ala(A11),Asp(B10),Lys(B28), Pro(B29)]insulin (DKP[A6-A11](Ala))]. Untethering the A6-A11 disulfide bridge by either amino acid causes similar perturbations in structure and dynamics as probed by circular dichroism and (1)H NMR spectroscopy. The analogues also exhibit similar decrements in thermodynamic stability relative to that of the parent monomer as probed by equilibrium denaturation studies (Delta Delta G(u) = 3.0 +/- 0.5 kcal/mol). Despite such similarities, the alanine analogue is 50 times more active than the serine analogue. Enhanced receptor binding (Delta Delta G = 2.2 kcal/mol) is in accord with alanine's greater helical propensity and more favorable hydrophobic-transfer free energy. The success of an induced-fit model highlights the applicability of general folding principles to a complex binding process. Comparison of DKP[A6-A11](Ser) and DKP[A6-A11](Ala) supports the hypothesis that the native A1-A8 alpha-helix functions as a preformed recognition element tethered by insulin's intrachain disulfide bridge. Segmental unfolding by design provides a novel approach to dissecting structure-activity relationships.     

Crystal structure of allo-Ile(A2)-insulin, an inactive chiral analogue: implications for the mechanism of receptor binding

The crystal structure of an inactive chiral analogue of insulin containing nonstandard substitution allo-Ile(A2) is described at 2.0 A resolution. In native insulin, the invariant Ile(A2) side chain anchors the N-terminal alpha-helix of the A-chain to the hydrophobic core. The structure of the variant protein was determined by molecular replacement as a T(3)R(3) zinc hexamer. Whereas respective T- and R-state main-chain structures are similar to those of native insulin (main-chain root-mean-square deviations (RMSD) of 0.45 and 0.54 A, respectively), differences in core packing are observed near the variant side chain. The R-state core resembles that of the native R-state with a local inversion of A2 orientation (core side chain RMSD 0.75 A excluding A2); in the T-state, allo-Ile(A2) exhibits an altered conformation in association with the reorganization of the surrounding side chains (RMSD 0.98 A). Surprisingly, the core of the R-state is similar to that observed in solution nuclear magnetic resonance (NMR) studies of an engineered T-like monomer containing the same chiral substitution (allo-Ile(A2)-DKP-insulin; Xu, B., Hua, Q. X., Nakagawa, S. H., Jia, W., Chu, Y. C., Katsoyannis, P. G., and Weiss, M. A. (2002) J. Mol. Biol. 316, 435-441). Simulation of NOESY spectra based on crystallographic protomers enables the analysis of similarities and differences in solution. The different responses of the T- and R-state cores to chiral perturbation illustrates both their intrinsic plasticity and constraints imposed by hexamer assembly. Although variant T- and R-protomers retain nativelike protein surfaces, the receptor-binding activity of allo-Ile(A2)-insulin is low (2% relative to native insulin). This seeming paradox suggests that insulin undergoes a change in conformation to expose Ile(A2) at the hormone-receptor interface.     

Nonlocal structural perturbations in a mutant human insulin: sequential resonance assignment and 13C-isotope-aided 2D-NMR studies of [PheB24-->Gly]insulin with implications for receptor recognition.

Insulin's mechanism of receptor binding is not well understood despite extensive study by mutagenesis and X-ray crystallography. Of particular interest are "anomalous" analogues whose bioactivities are not readily rationalized by crystal structures. Here the structure and dynamics of one such analogue (GlyB24-insulin) are investigated by circular dichroism (CD) and isotope-aided 2D-NMR spectroscopy. The mutant insulin retains near-native receptor-binding affinity despite a nonconservative substitution (PheB24-->Gly) in the receptor-binding surface. Relative to native insulin, GlyB24-insulin exhibits reduced dimerization; the monomer (the active species) exhibits partial loss of ordered structure, as indicated by CD studies and motional narrowing of selected 1H-NMR resonance. 2D-NMR studies demonstrate that the B-chain beta-turn (residues B20-23) and beta-strand (residues B24-B28) are destabilized; essentially native alpha-helical secondary structure (residues A3-A8, A13-A18, and B9-B19) is otherwise maintained. 13C-Isotope-edited NOESY studies demonstrate that long-range contacts observed between the B-chain beta-strand and the alpha-helical core in native insulin are absent in the mutant. Implications for the mechanism of insulin's interaction with its receptor are discussed.     

Importance of the solvent-exposed residues of the insulin B chain alpha-helix for receptor binding

Conjointly, the solvent-exposed residues of the central alpha-helix of the B chain form a well-defined ridge, which is flanked and partly overlapped by the two described insulin receptor binding surfaces on either side of the insulin molecule. To evaluate the importance of this interface in insulin receptor binding, we developed a new powerful method that allows us to introduce all the naturally occurring amino acids into a given position and subsequently determine the receptor binding affinities of the resulting insulin analogues. The total amino acid scanning mutagenesis was performed at positions B9, B10, B12, B13, B16, and B17, and the vast majority of the insulin analogue precursors were expressed and secreted in amounts close to that of the wild-type (human insulin) precursor. The analogue binding data revealed that positions B12 and B16 were the two positions most affected by the amino acid substitutions. Interestingly, the receptor binding affinities of the B13 analogues were also markedly affected by the amino acid substitutions, suggesting that GluB13 indeed is a part of insulin's binding surface. The B10 library screen generated analogues covering a wide range of (20-340%) of relative binding affinities, and the results indicated that a structural stabilization of the central alpha-helix and thereby a more rigid presentation of the binding epitope at the insulin receptor is important for receptor recognition. In conclusion, systematic amino acid scanning mutagenesis allowed us to confirm the importance of the B chain alpha-helix as a central recognition element serving as a linker of a continual binding surface.     

Toward understanding the role of insulin alpha-helix II in the growth-promoting activity of insulin.

For further understanding the contribution of the alpha-helix II (alphaII) in the growth-promoting activity of insulin, the residues A2Ile, A5Gln, and A8Thr located in alphaII were mutated to Leu, Glu, and Tyr, respectively. Three mutant insulins, [A2Leu]human insulin, [A5Glu]human insulin, and [A8Tyr]human insulin, were prepared by means of site-directed mutagenesis. The in vitro growth-promoting activities of the three mutant insulins, measured using GR2H6 cells, were 7.5%, 291%, and 250% of that of native insulin, respectively. Their receptor-binding activities to the insulin receptor were 2.3%, 46.7%, and 138.7%, respectively, compared with native insulin. Both the growth-promoting and receptor-binding activities of [A2Leu]human insulin and [A3Leu]insulin (Shi et al., 1997) were parallel and greatly decreased compared with native insulin. The results demonstrate that the residues A2Ile and A3Val in the alphaII are essential for the growth-promoting activity of insulin, and the growth-promoting function of insulin might be performed through, or mainly through, binding to the insulin receptor. The growth-promoting activities of [A5Glu]human insulin and [A8Tyr]human insulin were increased 6-fold and 2-fold, respectively, compared with native insulin, indicating that their growth-promoting activities might be expressed by, or mainly by, binding to the IGF-1 receptor.     

Mechanism of insulin chain combination. Asymmetric roles of A-chain alpha-helices in disulfide pairing

The A and B chains of insulin combine to form native disulfide bridges without detectable isomers. The fidelity of chain combination thus recapitulates the folding of proinsulin, a precursor protein in which the two chains are tethered by a disordered connecting peptide. We have recently shown that chain combination is blocked by seemingly conservative substitutions in the C-terminal alpha-helix of the A chain. Such analogs, once formed, nevertheless retain high biological activity. By contrast, we demonstrate here that chain combination is robust to non-conservative substitutions in the N-terminal alpha-helix. Introduction of multiple glycine substitutions into the N-terminal segment of the A chain (residues A1-A5) yields analogs that are less stable than native insulin and essentially without biological activity. (1)H NMR studies of a representative analog lacking invariant side chains Ile(A2) and Val(A3) (A chain sequence GGGEQCCTSICSLYQLENYCN; substitutions are italicized and cysteines are underlined) demonstrate local unfolding of the A1-A5 segment in an otherwise native-like structure. That this and related partial folds retain efficient disulfide pairing suggests that the native N-terminal alpha-helix does not participate in the transition state of the reaction. Implications for the hierarchical folding mechanisms of proinsulin and insulin-like growth factors are discussed.     

Inactive conformation of an insulin despite its wild-type sequence.

The peptide group between residues B24 and B25 of insulin was replaced by an ester bond. This modification only in the backbone was meant to eliminate a structurally important H-bond between the amide proton of B25 and the carbonyl oxygen of A19, and consequently to enhance detachment of the C-terminal B-chain from the body of the molecule, exposing the underlying A-chain. According to a model derived from the effects of side-chain substitutions, main-chain shortening, and crosslinking, this conformational change is prerequisite for receptor binding. Contrary to the expectation that increased flexibility would increase receptor binding and activity, depsi-insulin ([B24-B25 CO-O]insulin) has turned out be only 3-4% potent. In search of an explanation for this observation, the solution structure of depsi-insulin was determined by two-dimensional 1H-NMR spectroscopy. It was found that the loss of the B25-A19 H-bond does not entail detachment of the C-terminal B-chain. On the contrary, it is overcompensated by a gain in hydrophobic interaction achieved by insertion of the Phe B25 side chain into the molecule's core. This is possible because of increased rotational freedom in the backbone owing to the ester bond. Distortion of the B20-B23 turn and an altered direction of the distal B-chain are consequences that also affect self-association. The exceptional position of the B25 side chain is thus the key feature of the depsi-insulin structure. Being buried in the interior, it is not available for guiding the interaction with the receptor, a crucial role attributed to it by the model. This seems to be the main reason why the structure of depsi-insulin represents an inactive conformation.     

Extracellular loops 1 and 3 and their associated transmembrane regions of the calcitonin receptor-like receptor are needed for CGRP receptor function

The first and third extracellular loops (ECL) of G protein-coupled receptors (GPCRs) have been implicated in ligand binding and receptor function. This study describes the results of an alanine/leucine scan of ECLs 1 and 3 and loop-associated transmembrane (TM) domains of the secretin-like GPCR calcitonin receptor-like receptor which associates with receptor activity modifying protein 1 to form the CGRP receptor. Leu195Ala, Val198Ala and Ala199Leu at the top of TM2 all reduced αCGRP-mediated cAMP production and internalization; Leu195Ala and Ala199Leu also reduced αCGRP binding. These residues form a hydrophobic cluster within an area defined as the “minor groove” of rhodopsin-like GPCRs. Within ECL1, Ala203Leu and Ala206Leu influenced the ability of αCGRP to stimulate adenylate cyclase. In TM3, His219Ala, Leu220Ala and Leu222Ala have influences on αCGRP binding and cAMP production; they are likely to indirectly influence the binding site for αCGRP as well as having an involvement in signal transduction. On the exofacial surfaces of TMs 6 and 7, a number of residues were identified that reduced cell surface receptor expression, most noticeably Leu351Ala and Glu357Ala in TM6. The residues may contribute to the RAMP1 binding interface. Ile360Ala impaired αCGRP-mediated cAMP production. Ile360 is predicted to be located close to ECL2 and may facilitate receptor activation. Identification of several crucial functional loci gives further insight into the activation mechanism of this complex receptor system and may aid rational drug design.     

Double-sided ubiquitin binding of Hrs-UIM in endosomal protein sorting

Hrs has an essential role in sorting of monoubiquitinated receptors to multivesicular bodies for lysosomal degradation, through recognition of ubiquitinated receptors by its ubiquitin-interacting motif (UIM). Here, we present the structure of a complex of Hrs-UIM and ubiquitin at 1.7-A resolution. Hrs-UIM forms a single alpha-helix, which binds two ubiquitin molecules, one on either side. These two ubiquitin molecules are related by pseudo two-fold screw symmetry along the helical axis of the UIM, corresponding to a shift by two residues on the UIM helix. Both ubiquitin molecules interact with the UIM in the same manner, using the Ile44 surface, with equal binding affinities. Mutational experiments show that both binding sites of Hrs-UIM are required for efficient degradative protein sorting. Hrs-UIM belongs to a new subclass of double-sided UIMs, in contrast to its yeast homolog Vps27p, which has two tandem single-sided UIMs.     

Receptor recognition by a hepatitis B virus reveals a novel mode of high affinity virus-receptor interaction

The duck hepatitis B virus model system was used to elucidate the characteristics of receptor (carboxypeptidase D, gp180) interaction with polypeptides representing the receptor binding site in the preS part of the large viral surface protein. We demonstrate the pivotal role of carboxypeptidase D for virus entry and show its C-domain represents the virus attachment site, which binds preS with extraordinary affinity. Combining results from surface plasmon resonance spectroscopy and two-dimensional NMR analysis we resolved the contribution of preS sequence elements to complex stability and show that receptor binding potentially occurs in two steps. Initially, a short alpha-helix in the C-terminus of the receptor binding domain facilitates formation of a primary complex. This complex is stabilized sequentially, involving approximately 60 most randomly structured amino acids preceding the helix. Thus, hepadnaviruses exhibit a novel mechanism of high affinity receptor interaction by conserving the potential to adapt structure during binding rather than to preserve it per se. We propose that this process represents an alternative strategy to escape immune surveillance and the evolutionary pressure inherent in the compact hepadnaviral genome organization.     

Mechanism of the reaction catalyzed by mandelate racemase. 2. Crystal structure of mandelate racemase at 2.5-A resolution: identification of the active site and possible catalytic residues

The crystal structure of mandelate racemase (MR) has been solved at 3.0-A resolution by multiple isomorphous replacement and subsequently refined against X-ray diffraction data to 2.5-A resolution by use of both molecular dynamics refinement (XPLOR) and restrained least-squares refinement (PROLSQ). The current crystallographic R-factor for this structure is 18.3%. MR is composed of two major structural domains and a third, smaller, C-terminal domain. The N-terminal domain has an alpha + beta topology consisting of a three-stranded antiparallel beta-sheet followed by an antiparallel four alpha-helix bundle. The central domain is a singly wound parallel alpha/beta-barrel composed of eight central strands of beta-sheet and seven alpha-helices. The C-terminal domain consists of an irregular L-shaped loop with several short sections of antiparallel beta-sheet and two short alpha-helices. This C-terminal domain partially covers the junction between the major domains and occupies a region of the central domain that is filled by an eight alpha-helix in all other known parallel alpha/beta-barrels except for the barrel domain in muconate lactonizing enzyme (MLE) [Goldman, A., Ollis, D. L., & Steitz, T. A. (1987) J. Mol. Biol. 194, 143] whose overall polypeptide fold and amino acid sequence are strikingly similar to those of MR [Neidhart, D. J., Kenyon, G. L., Gerlt, J. A., & Petsko, G. A. (1990) Nature 347, 692]. In addition, the crystal structure reveals that, like MLE, MR is tightly packed as an octamer of identical subunits. The active site of MR is located between the two major domains, at the C-terminal ends of the beta-strands in the alpha/beta-barrel domain. The catalytically essential divalent metal ion is ligated by three side-chain carboxyl groups contributed by residues of the central beta-sheet. A model of a productive substrate complex of MR has been constructed on the basis of difference Fourier analysis at 3.5-A resolution of a complex between MR and (R,S)-p-iodomandelate, permitting identification of residues that may participate in substrate binding and catalysis. The ionizable groups of both Lys 166 and His 297 are positioned to interact with the chiral center of substrate, suggesting that both of these residues may function as acid/base catalysts.(ABSTRACT TRUNCATED AT 400 WORDS)     

The crystal structure of the human (S100A8/S100A9)2 heterotetramer, calprotectin, illustrates how conformational changes of interacting alpha-helices can determine specific association of two EF-hand proteins

The EF-hand proteins S100A8 and S100A9 are important calcium signalling proteins that are involved in wound healing and provide clinically relevant markers of inflammatory processes, such as rheumatoid arthritis and inflammatory bowel disease. Both can form homodimers via distinct modes of association, probably of lesser stability in the case of S100A9, whereas in the presence of calcium S100A8 and S100A9 associate to calprotectin, the physiologically active heterooligomer. Here we describe the crystal structure of the (S100A8/S100A9)(2) heterotetramer at 1.8 A resolution. Its quaternary structure illustrates how specific heteroassociation is energetically driven by a more extensive burial of solvent accessible surface areas in both proteins, most pronounced for S100A9, thus leading to a dimer of heterodimers. A major contribution to tetramer association is made by the canonical calcium binding loops in the C-terminal halves of the two proteins. The mode of heterodimerisation in calprotectin more closely resembles the subunit association previously observed in the S100A8 homodimer and provides trans stabilisation for S100A9, which manifests itself in a significantly elongated C-terminal alpha-helix in the latter. As a consequence, two different putative zinc binding sites emerge at the S100A8/S100A9 subunit interface. One of these corresponds to a high affinity arrangement of three His residues and one Asp side-chain, which is unique to the heterotetramer. This structural feature explains the well known Zn(2+) binding activity of calprotectin, whose overexpression can cause strong dysregulation of zinc homeostasis with severe clinical symptoms.     

NMR solution structure of the receptor binding domain of human alpha(2)-macroglobulin

Human alpha(2)-macroglobulin-proteinase complexes bind to their receptor, the low density lipoprotein receptor-related protein (LRP), through a discrete 138-residue C-terminal receptor binding domain (RBD), which also binds to the beta-amyloid peptide. We have used NMR spectroscopy on recombinantly expressed uniformly (13)C/(15)N-labeled human RBD to determine its three-dimensional structure in solution. Human RBD is a sandwich of two antiparallel beta-sheets, one four-strand and one five-strand, and also contains one alpha-helix of 2.5 turns and an additional 1-turn helical region. The principal alpha-helix contains two lysine residues on the outer face that are known to be essential for receptor binding. A calcium binding site (K(d) approximately 11 mM) is present in the loop region at one end of the beta-sandwich. Calcium binding principally affects this loop region and does not significantly perturb the stable core structure of the domain. The structure and NMR assignments will enable us to examine in solution specific binding of RBD to domains of the receptor and to beta-amyloid peptide.     

Electrostatic interactions in leucine zippers: thermodynamic analysis of the contributions of Glu and His residues and the effect of mutating salt bridges

Electrostatic interactions play a complex role in stabilizing proteins. Here, we present a rigorous thermodynamic analysis of the contribution of individual Glu and His residues to the relative pH-dependent stability of the designed disulfide-linked leucine zipper AB(SS). The contribution of an ionized side-chain to the pH-dependent stability is related to the shift of the pK(a) induced by folding of the coiled coil structure. pK(a)(F) values of ten Glu and two His side-chains in folded AB(SS) and the corresponding pK(a)(U) values in unfolded peptides with partial sequences of AB(SS) were determined by 1H NMR spectroscopy: of four Glu residues not involved in ion pairing, two are destabilizing (-5.6 kJ mol(-1)) and two are interacting with the positive alpha-helix dipoles and are thus stabilizing (+3.8 kJ mol(-1)) in charged form. The two His residues positioned in the C-terminal moiety of AB(SS) interact with the negative alpha-helix dipoles resulting in net stabilization of the coiled coil conformation carrying charged His (-2.6 kJ mol(-1)). Of the six Glu residues involved in inter-helical salt bridges, three are destabilizing and three are stabilizing in charged form, the net contribution of salt-bridged Glu side-chains being destabilizing (-1.1 kJ mol(-1)). The sum of the individual contributions of protonated Glu and His to the higher stability of AB(SS) at acidic pH (-5.4 kJ mol(-1)) agrees with the difference in stability determined by thermal unfolding at pH 8 and pH 2 (-5.3 kJ mol(-1)). To confirm salt bridge formation, the positive charge of the basic partner residue of one stabilizing and one destabilizing Glu was removed by isosteric mutations (Lys-->norleucine, Arg-->norvaline). Both mutations destabilize the coiled coil conformation at neutral pH and increase the pK(a) of the formerly ion-paired Glu side-chain, verifying the formation of a salt bridge even in the case where a charged side-chain is destabilizing. Because removing charges by a double mutation cycle mainly discloses the immediate charge-charge effect, mutational analysis tends to overestimate the overall energetic contribution of salt bridges to protein stability.     

Crystal structure of the third extracellular domain of CD5 reveals the fold of a group B scavenger cysteine-rich receptor domain

Scavenger receptor cysteine-rich (SRCR) domains are ancient protein modules widely found among cell surface and secreted proteins of the innate and adaptive immune system, where they mediate ligand binding. We have solved the crystal structure at 2.2 A of resolution of the SRCR CD5 domain III, a human lymphocyte receptor involved in the modulation of antigen specific receptor-mediated T cell activation and differentiation signals. The first structure of a member of a group B SRCR domain reveals the fold of this ancient protein module into a central core formed by two antiparallel beta-sheets and one alpha-helix, illustrating the conserved core at the protein level of genes coding for group A and B members of the SRCR superfamily. The novel SRCR group B structure permits the interpretation of site-directed mutagenesis data on the binding of activated leukocyte cell adhesion molecule (ALCAM/CD166) binding to CD6, a closely related lymphocyte receptor homologue to CD5.     

Effects of side-chain characteristics on stability and oligomerization state of a de novo-designed model coiled-coil: 20 amino acid substitutions in position "d"

We describe the de novo design and biophysical characterization of a model coiled-coil protein in which we have systematically substituted 20 different amino acid residues in the central "d" position. The model protein consists of two identical 38 residue polypeptide chains covalently linked at their N termini via a disulfide bridge. The hydrophobic core contained Val and Ile residues at positions "a" and Leu residues at positions "d". This core allowed for the formation of both two-stranded and three-stranded coiled-coils in benign buffer, depending on the substitution at position "d". The structure of each analog was analyzed by CD spectroscopy and their relative stability determined by chemical denaturation using GdnHCI (all analogs denatured from the two-stranded state). The oligomeric state(s) was determined by high-performance size-exclusion chromatography and sedimentation equilibrium analysis in benign medium. Our results showed a thermodynamic stability order (in order of decreasing stability) of: Leu, Met, Ile, Tyr, Phe, Val, Gln, Ala, Trp, Asn, His, Thr, Lys, Ser, Asp, Glu, Arg, Orn, and Gly. The Pro analog prevented coiled-coil formation. The overall stability range was 7.4 kcal/mol from the lowest to the highest analog, indicating the importance of the hydrophobic core and the dramatic effect a single substitution in the core can have upon the stability of the protein fold. In general, the side-chain contribution to the level of stability correlated with side-chain hydrophobicity. Molecular modelling studies, however, showed that packing effects could explain deviations from a direct correlation. In regards to oligomerization state, eight analogs demonstrated the ability to populate exclusively one oligomerization state in benign buffer (0.1 M KCl, 0.05 M K(2)PO(4)(pH 7)). Ile and Val (the beta-branched residues) induced the three-stranded oligomerization state, whereas Tyr, Lys, Arg, Orn, Glu and Asp induced the two-stranded state. Asn, Gln, Ser, Ala, Gly, Phe, Leu, Met and Trp analogs were indiscriminate and populated two-stranded and three-stranded states. Comparison of these results with similar substitutions in position "a" highlights the positional effects of individual residues in defining the stability and numbers of polypeptide chains occurring in a coiled-coil structure. Overall, these results in conjunction with other work now generate a relative thermodynamic stability scale for 19 naturally occurring amino acid residues in either an "a" or "d" position of a two-stranded coiled-coil. Thus, these results will aid in the de novo design of new coiled-coil structures, a better understanding of their structure/function relationships and the design of algorithms to predict the presence of coiled-coils within native protein sequences.     

Contribution of Thr29 to the thermodynamic stability of goat alpha-lactalbumin as determined by experimental and theoretical approaches

The Thr29 residue in the hydrophobic core of goat alpha-lactalbumin (alpha-LA) was substituted with Val (Thr29Val) and Ile (Thr29Ile) to investigate the contribution of Thr29 to the thermodynamic stability of the protein. We carried out protein stability measurements, X-ray crystallographic analyses, and free energy calculations based on molecular dynamics simulation. The equilibrium unfolding transitions induced by guanidine hydrochloride demonstrated that the Thr29Val and Thr29Ile mutants were, respectively, 1.9 and 3.2 kcal/mol more stable than the wild-type protein (WT). The overall structures of the mutants were almost identical to that of WT, in spite of the disruption of the hydrogen bonding between the side-chain O-H group of Thr29 and the main-chain C=O group of Glu25. To analyze the stabilization mechanism of the mutants, we performed free energy calculations. The calculated free energy differences were in good agreement with the experimental values. The stabilization of the mutants was mainly caused by solvation loss in the denatured state. Furthermore, the O-H group of Thr29 favorably interacts with the C=O group of Glu25 to form hydrogen bonds and, simultaneously, unfavorably interacts electrostatically with the main-chain C=O group of Thr29. The difference in the free energy profile of the unfolding path between WT and the Thr29Ile mutant is discussed in light of our experimental and theoretical results.     

Protein recognition motifs: design of peptidomimetics of helix surfaces

Helices represent one of the most common recognition motifs in proteins. The design of nonpeptidic scaffolds, such as the 3,2',2'-tris-substituted terphenyl, that can imitate the side-chain orientation along one face of an alpha-helix potentially provides an effective means to modulate helix-recognition functions. Here, based on theoretical arguments, we described novel alpha-helix mimetics which are more effective than the terphenyl at constraining the aryl-aryl torsion angles to those associated with structures suitable for mimicking the alpha-helical twist for side-chain orientation and for superimposing the side chains of residues i, i + 3 or i + 4, i + 7 when compared with the alpha-beta side-chain vectors of the regular alpha-helix with an improved root mean square deviation (RMSD) of approximately 0.5 A. In addition, this study suggests that rotamer distributions around the C(alpha)--C(beta) bonds of these helix mimetics are similar to those of alpha-helices, except that these rotamer distributions show an approximately 60 degrees shift compared to those of alpha-helices when the mimetic axis is superimposed upon the helix axis. This change in rotamer orientation complicates mimicry of the helix surface.