A cross-strand Trp Trp pair stabilizes the hPin1 WW domain at the expense of function

Using the human Pin1 WW domain (hPin1 WW), we show that replacement of two nearest neighbor non-hydrogen-bonded residues on adjacent beta-strands with tryptophan (Trp) residues increases beta-sheet thermodynamic stability by 4.8 kJ mol(-1) at physiological temperature. One-dimensional NMR studies confirmed that introduction of the Trp-Trp pair does not globally perturb the structure of the triple-stranded beta-sheet, while circular dichroism studies suggest that the engineered cross-strand Trp-Trp pair adopts a side-chain conformation similar to that first reported for a designed "Trp-zipper" beta-hairpin peptide, wherein the indole side chains stack perpendicular to each other. Even though the mutated side chains in wild-type hPin1 WW are not conserved among WW domains and compose the beta-sheet surface opposite to that responsible for ligand binding, introduction of the cross-strand Trp-Trp pair effectively eliminates hPin1 WW function as assessed by the loss of binding affinity toward a natural peptide ligand. Maximizing both thermodynamic stability and the domain function of hPin1 WW by the above mentioned approach appears to be difficult, analogous to the situation with loop 1 optimization explored previously. That introduction of a non-hydrogen-bonded cross-strand Trp-Trp pair within the hPin1 WW domain eliminates function may provide a rationale for why this energetically favorable pairwise interaction has not yet been identified in WW domains or any other biologically evolved protein with known three-dimensional structure.     

Influence of hPin1 WW N-terminal domain boundaries on function, protein stability, and folding

An N-terminally truncated and cooperatively folded version (residues 6-39) of the human Pin1 WW domain (hPin1 WW hereafter) has served as an excellent model system for understanding triple-stranded beta-sheet folding energetics. Here we report that the negatively charged N-terminal sequence (Met1-Ala-Asp-Glu-Glu5) previously deleted, and which is not conserved in highly homologous WW domain family members from yeast or certain fungi, significantly increases the stability of hPin1 WW (approximately 4 kJ mol(-1) at 65 degrees C), in the context of the 1-39 sequence based on equilibrium measurements. N-terminal truncations and mutations in conjunction with a double mutant cycle analysis and a recently published high-resolution X-ray structure of the hPin1 cis/trans-isomerase suggest that the increase in stability is due to an energetically favorable ionic interaction between the negatively charged side chains in the N terminus of full-length hPin1 WW and the positively charged epsilon-ammonium group of residue Lys13 in beta-strand 1. Our data therefore suggest that the ionic interaction between Lys13 and the charged N terminus is the optimal solution for enhanced stability without compromising function, as ascertained by ligand binding studies. Kinetic laser temperature-jump relaxation studies reveal that this stabilizing interaction has not formed to a significant extent in the folding transition state at near physiological temperature, suggesting a differential contribution of the negatively charged N-terminal sequence to protein stability and folding rate. As neither the N-terminal sequence nor Lys13 are highly conserved among WW domains, our data further suggest that caution must be exercised when selecting domain boundaries for WW domains for structural, functional, or thermodynamic studies.     

The role of backbone conformational heat capacity in protein stability: temperature dependent dynamics of the B1 domain of Streptococcal protein G

The contributions of backbone NH group dynamics to the conformational heat capacity of the B1 domain of Streptococcal protein G have been estimated from the temperature dependence of 15N NMR-derived order parameters. Longitudinal (R1) and transverse (R2) relaxation rates, transverse cross-relaxation rates (eta(xy)), and steady state [1H]-15N nuclear Overhauser effects were measured at temperatures of 0, 10, 20, 30, 40, and 50 degrees C for 89-100% of the backbone secondary amide nitrogen nuclei in the B1 domain. The ratio R2/eta(xy) was used to identify nuclei for which conformational exchange makes a significant contribution to R2. Relaxation data were fit to the extended model-free dynamics formalism, incorporating an axially symmetric molecular rotational diffusion tensor. The temperature dependence of the order parameter (S2) was used to calculate the contribution of each NH group to conformational heat capacity (Cp) and a characteristic temperature (T*), representing the density of conformational energy states accessible to each NH group. The heat capacities of the secondary structure regions of the B1 domain are significantly higher than those of comparable regions of other proteins, whereas the heat capacities of less structured regions are similar to those in other proteins. The higher local heat capacities are estimated to contribute up to approximately 0.8 kJ/mol K to the total heat capacity of the B1 domain, without which the denaturation temperature would be approximately 9 degrees C lower (78 degrees C rather than 87 degrees C). Thus, variation of backbone conformational heat capacity of native proteins may be a novel mechanism that contributes to high temperature stabilization of proteins.     

Synthesis, folding, and structure of the beta-turn mimic modified B1 domain of streptococcal protein G

The mechanism of beta-sheet formation remains a fundamental issue in our understanding of the protein folding process, but is hampered by the often encountered kinetic competition between folding and aggregation. The role of local versus nonlocal interactions has been probed traditionally by mutagenesis of both turn and strand residues. Recently, rigid organic molecules that impose a correct chain reversal have been introduced in several small peptides to isolate the importance of the long-range interactions. Here, we present the incorporation of a well-studied beta-turn mimic, designated as the dibenzofuran-based (DBF) amino acid, in the B1 domain of streptococcal protein G (B1G), and compare our results with those obtained upon insertion of the same mimic into the N-terminal beta-hairpin of B1G (O Melnyk et al., 1998, Lett Pept Sci 5:147-150). The DBF-B1G domain conserves the structure and the functional and thermodynamical properties of the native protein, whereas the modified peptide does not adopt a native-like conformation. The nature of the DBF flanking residues in the modified B1G domain prevents the beta-turn mimic from acting as a strong beta-sheet nucleator, which reinforces the idea that the native beta-hairpin formation is not driven by the beta-turn formation, but by tertiary interactions.     

Factors involved in the stability of isolated beta-sheets: Turn sequence, beta-sheet twisting, and hydrophobic surface burial

We have recently reported on the design of a 20-residue peptide able to form a significant population of a three-stranded up-and-down antiparallel beta-sheet in aqueous solution. To improve our beta-sheet model in terms of the folded population, we have modified the sequences of the two 2-residue turns by introducing the segment DPro-Gly, a sequence shown to lead to more rigid type II' beta-turns. The analysis of several NMR parameters, NOE data, as well as Deltadelta(CalphaH), DeltadeltaC(beta), and Deltadelta(Cbeta) values, demonstrates that the new peptide forms a beta-sheet structure in aqueous solution more stable than the original one, whereas the substitution of the DPro residues by LPro leads to a random coil peptide. This agrees with previous results on beta-hairpin-forming peptides showing the essential role of the turn sequence for beta-hairpin folding. The well-defined beta-sheet motif calculated for the new designed peptide (pair-wise RMSD for backbone atoms is 0.5 +/- 0.1 A) displays a high degree of twist. This twist likely contributes to stability, as a more hydrophobic surface is buried in the twisted beta-sheet than in a flatter one. The twist observed in the up-and-down antiparallel beta-sheet motifs of most proteins is less pronounced than in our designed peptide, except for the WW domains. The additional hydrophobic surface burial provided by beta-sheet twisting relative to a "flat" beta-sheet is probably more important for structure stability in peptides and small proteins like the WW domains than in larger proteins for which there exists a significant contribution to stability arising from their extensive hydrophobic cores.     

A molecular dynamics study of the 41-56 beta-hairpin from B1 domain of protein G.

The structural and dynamical behavior of the 41-56 beta-hairpin from the protein G B1 domain (GB1) has been studied at different temperatures using molecular dynamics (MD) simulations in an aqueous environment. The purpose of these simulations is to establish the stability of this hairpin in view of its possible role as a nucleation site for protein folding. The conformation of the peptide in the crystallographic structure of the protein GB1 (native conformation) was lost in all simulations. The new equilibrium conformations are stable for several nanoseconds at 300K (>10 ns), 350 K (>6.5 ns), and even at 450 K (up to 2.5 ns). The new structures have very similar hairpin-like conformations with properties in agreement with available experimental nuclear Overhauser effect (NOE) data. The stability of the structure in the hydrophobic core region during the simulations is consistent with the experimental data and provides further evidence for the role played by hydrophobic interactions in hairpin structures. Essential dynamics analysis shows that the dynamics of the peptide at different temperatures spans basically the same essential subspace. The main equilibrium motions in this subspace involve large fluctuations of the residues in the turn and ends regions. Of the six interchain hydrogen bonds, the inner four remain stable during the simulations. The space spanned by the first two eigenvectors, as sampled at 450 K, includes almost all of the 47 different hairpin structures found in the database. Finally, analysis of the hydration of the 300 K average conformations shows that the hydration sites observed in the native conformation are still well hydrated in the equilibrium MD ensemble.     

Conformational studies of the C-terminal 16-amino-acid-residue fragment of the B3 domain of the immunoglobulin binding protein G from Streptococcus

The structure and stability of the 16-amino-acid-residue fragment [IG(46-61)] corresponding to the C-terminal beta-hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus was investigated by means of CD and NMR spectroscopy and by differential scanning calorimetry. The CD and 2D NMR experiments were carried out (i) in water at different temperatures and (ii) at one temperature (305 K), with only CD, at different TFE concentrations. Our results show that the IG(46-61) peptide possesses organized three-dimensional structure at all investigated temperatures. The three-dimensional structure of the IG(46-61) peptide resembles the general shape of a beta-hairpin that is also observed for this peptide in the experimental structure of the B3 domain in the whole G protein; the structure is stabilized by hydrophobic interactions between nonpolar side chains. Our study shows that the melting temperature of the IG(46-61) peptide is about 320 K which supports the hypothesis that the investigated peptide can serve as a folding initiation site of the B3 domain of the immunoglobulin binding protein G.     

The role of side chain conformational flexibility in surface recognition by Tenebrio molitor antifreeze protein

Two-dimensional nuclear magnetic resonance spectroscopy was used to investigate the flexibility of the threonine side chains in the beta-helical Tenebrio molitor antifreeze protein (TmAFP) at low temperatures. From measurement of the (3)J(alphabeta) (1)H-(1)H scalar coupling constants, the chi(1) angles and preferred rotamer populations can be calculated. It was determined that the threonines on the ice-binding face of the protein adopt a preferred rotameric conformation at near freezing temperatures, whereas the threonines not on the ice-binding face sample many rotameric states. This suggests that TmAFP maintains a preformed ice-binding conformation in solution, wherein the rigid array of threonines that form the AFP-ice interface matches the ice crystal lattice. A key factor in binding to the ice surface and inhibition of ice crystal growth appears to be the close surface-to-surface complementarity between the AFP and crystalline ice, and the lack of an entropic penalty associated with freezing out motions in a flexible ligand.     

Molecular dynamics as a tool to detect protein foldability. A mutant of domain B1 of protein G with non-native secondary structure propensities

The usefulness of molecular dynamics to assess the structural integrity of mutants containing several mutations has been investigated. Our goal was to determine whether molecular dynamics would be able to discriminate mutants of a protein having a close-to-wild-type fold, from those that are not folded under the same conditions. We used as a model the B1 domain of protein G in which we replaced the unique central alpha-helix by the sequence of the second beta-hairpin, which has a strong intrinsic propensity to form this secondary structure in solution. In the resulting protein, one-third of the secondary structure has been replaced by a non-native one. Models of the mutants were built based on the three-dimensional structure of the wild-type GB1 domain. During 2 ns of molecular dynamics simulations on these models, mutants containing up to 10 mutations in the helix retained the native fold, while another mutant with an additional mutation unfolded. This result is in agreement with our circular dichroism and NMR experiments, which indicated that the former mutants fold into a structure similar to the wild-type, as opposed to the latter mutant which is partly unfolded. Additionally, a mutant containing six mutations scattered through the surface of the domain, and which is unfolded, was also detected by the simulation. This study suggests that molecular dynamics calculations could be performed on molecular models of mutants of a protein to evaluate their foldability, prior to a mutagenesis experiment.     

Atomic-resolution crystal structure of Borrelia burgdorferi outer surface protein A via surface engineering

Outer surface protein A (OspA) from Borrelia burgdorferi has an unusual dumbbell-shaped structure in which two globular domains are connected with a "single-layer" beta-sheet (SLB). The protein is highly soluble, and it has been recalcitrant to crystallization. Only OspA complexes with Fab fragments have been successfully crystallized. OspA contains a large number of Lys and Glu residues, and these "high entropy" residues may disfavor crystal packing because some of them would need to be immobilized in forming a crystal lattice. We rationally designed a total of 13 surface mutations in which Lys and Glu residues were replaced with Ala or Ser. We successfully crystallized the mutant OspA without a bound Fab fragment and extended structure analysis to a 1.15 Angstroms resolution. The new high-resolution structure revealed a unique backbone hydration pattern of the SLB segment in which water molecules fill the "weak spots" on both faces of the antiparallel beta-sheet. These well-defined water molecules provide additional structural links between adjacent beta-strands, and thus they may be important for maintaining the rigidity of the SLB that inherently lacks tight packing afforded by a hydrophobic core. The structure also revealed new information on the side-chain dynamics and on a solvent-accessible cavity in the core of the C-terminal globular domain. This work demonstrates the utility of extensive surface mutation in crystallizing recalcitrant proteins and dramatically improving the resolution of crystal structures, and provides new insights into the stabilization mechanism of OspA.     

Analysis of the factors that stabilize a designed two-stranded antiparallel beta-sheet

Autonomously folding beta-hairpins (two-strand antiparallel beta-sheets) have become increasingly valuable tools for probing the forces that control peptide and protein conformational preferences. We examine the effects of variations in sequence and solvent on the stability of a previously designed 12-residue peptide (1). This peptide adopts a beta-hairpin conformation containing a two-residue loop (D-Pro-Gly) and a four-residue interstrand sidechain cluster that is observed in the natural protein GB1. We show that the conformational propensity of the loop segment plays an important role in beta-hairpin stability by comparing 1 with (D)P--> N mutant 2. In addition, we show that the sidechain cluster contributes both to conformational stability and to folding cooperativity by comparing 1 with mutant 3, in which two of the four cluster residues have been changed to serine. Thermodynamic analysis suggests that the high loop-forming propensity of the (D)PG segment decreases the entropic cost of beta-hairpin formation relative to the more flexible NG segment, but that the conformational rigidity of (D)PG may prevent optimal contacts between the sidechains of the GB1-derived cluster. The enthalpic favorability of folding in these designed beta-hairpins suggests that they are excellent scaffolds for studying the fundamental mechanisms by which amino acid sidechains interact with one another in folded proteins.     

Crystal structures and increased stabilization of the protein G variants with switched folding pathways NuG1 and NuG2

We recently described two protein G variants (NuG1 and NuG2) with redesigned first hairpins that were almost twice as stable, folded 100-fold faster, and had a switched folding mechanism relative to the wild-type protein. To test the structural accuracy of our design algorithm and to provide insights to the dramatic changes in the kinetics and thermodynamics of folding, we have now determined the crystal structures of NuG1 and NuG2 to 1.8 A and 1.85 A, respectively. We find that they adopt hairpin structures that are closer to the computational models than to wild-type protein G; the RMSD of the NuG1 hairpin to the design model and the wild-type structure are 1.7 A and 5.1 A, respectively. The crystallographic B factor in the redesigned first hairpin of NuG1 is systematically higher than the second hairpin, suggesting that the redesigned region is somewhat less rigid. A second round of structure-based design yielded new variants of NuG1 and NuG2, which are further stabilized by 0.5 kcal/mole and 0.9 kcal/mole.     

Structural and dynamic characterization of an unfolded state of poplar apo-plastocyanin formed under nondenaturing conditions

Plastocyanin is a predominantly beta-sheet protein containing a type I copper center. The conformational ensemble of a denatured state of apo-plastocyanin formed in solution under conditions of low salt and neutral pH has been investigated by multidimensional heteronuclear NMR spectroscopy. Chemical shift assignments were obtained by using three-dimensional triple-resonance NMR experiments to trace through-bond heteronuclear connectivities along the backbone and side chains. The (3)J(HN,Halpha) coupling constants, (15)N-edited proton-proton nuclear Overhauser effects (NOEs), and (15)N relaxation parameters were also measured for the purpose of structural and dynamic characterization. Most of the residues corresponding to beta-strands in the folded protein exhibit small upfield shifts of the (13)C(alpha) and (13)CO resonances relative to random coil values, suggesting a slight preference for backbone dihedral angles in the beta region of (phi,psi) space. This is further supported by the presence of strong sequential d(alphaN)(i, i + 1) NOEs throughout the sequence. The few d(NN)(i, i + 1) proton NOEs that are observed are mostly in regions that form loops in the native plastocyanin structure. No medium or long-range NOEs were observed. A short sequence, between residues 59 and 63, was found to populate a nonnative helical conformation in the unfolded state, as indicated by the shift of the (13)C(alpha), (13)CO, and (1)H(alpha) resonances relative to random coil values and by the decreased values of the (3)J(HN,Halpha) coupling constants. The (15)N relaxation parameters indicate restriction of motions on a nanosecond timescale in this region. Intriguingly, this helical conformation is present in a sequence that is close to but not in the same location as the single short helix in the native folded protein. The results are consistent with earlier NMR studies of peptide fragments of plastocyanin and confirm that the regions of the sequence that form beta-strands in the native protein spontaneously populate the beta-region of (phi,psi) space under folding conditions, even in the absence of stabilizing tertiary interactions. We conclude that the state of apo-plastocyanin present under nondenaturing conditions is a noncompact unfolded state with some evidence of nativelike and nonnative local structuring that may be initiation sites for folding of the protein.     

Simulation of pH-dependent edge strand rearrangement in human beta-2 microglobulin

Amyloid fibrils formed from unrelated proteins often share morphological similarities, suggesting common biophysical mechanisms for amyloidogenesis. Biochemical studies of human beta-2 microglobulin (beta2M) have shown that its transition from a water-soluble protein to insoluble aggregates can be triggered by low pH. Additionally, biophysical measurements of beta2M using NMR have identified residues of the protein that participate in the formation of amyloid fibrils. The crystal structure of monomeric human beta2M determined at pH 5.7 shows that one of its edge beta-strands (strand D) adopts a conformation that differs from other structures of the same protein obtained at higher pH. This alternate beta-strand arrangement lacks a beta-bulge, which may facilitate protein aggregation through intermolecular beta-sheet association. To explore whether the pH change may yield the observed conformational difference, molecular dynamics simulations of beta2M were performed. The effects of pH were modeled by specifying the protonation states of Asp, Glu, and His, as well as the C terminus of the main chain. The bulged conformation of strand D is preferred at medium pH (pH 5-7), whereas at low pH (pH < 4) the straight conformation is observed. Therefore, low pH may stabilize the straight conformation of edge strand D and thus increase the amyloidogenicity of beta2M.     

Conformational stability and thermodynamic characterization of the lipoic acid bearing domain of human mitochondrial branched chain alpha-ketoacid dehydrogenase

The lipoic acid bearing domain (hbLBD) of human mitochondrial branched chain alpha-ketoacid dehydrogenase (BCKD) plays important role of substrate channeling in oxidative decarboxylation of the branched chain alpha-ketoacids. Recently hbLBD has been found to follow two-step folding mechanism without detectable presence of stable or kinetic intermediates. The present study describes the conformational stability underlying the folding of this small beta-barrel domain. Thermal denaturation in presence of urea and isothermal urea denaturation titrations are used to evaluate various thermodynamic parameters defining the equilibrium unfolding. The linear extrapolation model successfully describes the two-step; native state <-->denatured state unfolding transition of hbLBD. The average temperature of maximum stability of hbLBD is estimated as 295.6 +/- 0.9 K. Cold denaturation of hbLBD is also predicted and discussed.     

Molecular dynamics simulations of alanine rich beta-sheet oligomers: Insight into amyloid formation

The aggregation observed in protein conformational diseases is the outcome of significant new beta-sheet structure not present in the native state. Peptide model systems have been useful in studies of fibril aggregate formation. Experimentally, it was found that a short peptide AGAAAAGA is one of the most highly amyloidogenic peptides. This peptide corresponds to the Syrian hamster prion protein (ShPrP) residues 113-120. The peptide was observed to be conserved in all species for which the PrP sequence has been determined. We have simulated the stabilities of oligomeric AGAAAAGA and AAAAAAAA (A8) by molecular dynamic simulations. Oligomers of both AGAAAAGA and AAAAAAAA were found to be stable when the size is 6 to 8 (hexamer to octamer). Subsequent simulation of an additional alpha-helical AAAAAAAA placed on the A8-octamer surface has revealed molecular events related to conformational change and oligomer growth. Our study addresses both the minimal oligomeric size of an aggregate seed and the mechanism of seed growth. Our simulations of the prion-derived 8-residue amyloidogenic peptide and its variant have indicated that an octamer is stable enough to be a seed and that the driving force for stabilization is the hydrophobic effect.     

The Bacillus licheniformis BlaP beta-lactamase as a model protein scaffold to study the insertion of protein fragments

Using genetic engineering technologies, the chitin-binding domain (ChBD) of the human macrophage chitotriosidase has been inserted into the host protein BlaP, a class A beta-lactamase produced by Bacillus licheniformis. The product of this construction behaved as a soluble chimeric protein that conserves both the capacity to bind chitin and to hydrolyze beta-lactam moiety. Here we describe the biochemical and biophysical properties of this protein (BlaPChBD). This work contributes to a better understanding of the reciprocal structural and functional effects of the insertion on the host protein scaffold and the heterologous structured protein fragments. The use of BlaP as a protein carrier represents an efficient approach to the functional study of heterologous protein fragments.     

Factors contributing to decreased protein stability when aspartic acid residues are in beta-sheet regions

Asp residues are significantly under represented in beta-sheet regions of proteins, especially in the middle of beta-strands, as found by a number of studies using statistical, modeling, or experimental methods. To further understand the reasons for this under representation of Asp, we prepared and analyzed mutants of a beta-domain. Two Gln residues of the immunoglobulin light-chain variable domain (V(L)) of protein Len were replaced with Asp, and then the effects of these changes on protein stability and protein structure were studied. The replacement of Q38D, located at the end of a beta-strand, and that of Q89D, located in the middle of a beta-strand, reduced the stability of the parent immunoglobulin V(L) domain by 2.0 kcal/mol and 5.3 kcal/mol, respectively. Because the Q89D mutant of the wild-type V(L)-Len domain was too unstable to be expressed as a soluble protein, we prepared the Q89D mutant in a triple mutant background, V(L)-Len M4L/Y27dD/T94H, which was 4.2 kcal/mol more stable than the wild-type V(L)-Len domain. The structures of mutants V(L)-Len Q38D and V(L)-Len Q89D/M4L/Y27dD/T94H were determined by X-ray diffraction at 1.6 A resolution. We found no major perturbances in the structures of these Q-->D mutant proteins relative to structures of the parent proteins. The observed stability changes have to be accounted for by cumulative effects of the following several factors: (1) by changes in main-chain dihedral angles and in side-chain rotomers, (2) by close contacts between some atoms, and, most significantly, (3) by the unfavorable electrostatic interactions between the Asp side chain and the carbonyls of the main chain. We show that the Asn side chain, which is of similar size but neutral, is less destabilizing. The detrimental effect of Asp within a beta-sheet of an immunoglobulin-type domain can have very serious consequences. A somatic mutation of a beta-strand residue to Asp could prevent the expression of the domain both in vitro and in vivo, or it could contribute to the pathogenic potential of the protein in vivo.