Elongation of the BH8 beta-hairpin peptide: Electrostatic interactions in beta-hairpin formation and stability

An elongated version of the de novo designed beta-hairpin peptide, BH8, has allowed us to gain insight into the role of electrostatic interactions in beta-hairpin stability. A Lys-Glu electrostatic pair has been introduced by adding a residue at the beginning and at the end of the N-terminal and C-terminal strands, respectively, of the beta-hairpin structure, in both orientations. The two resulting peptides and controls having Ala residues at these positions and different combinations of Ala with Lys, or Glu residues, have been analyzed by nuclear magnetic resonance (NMR), under different pH and ionic strength conditions. All of the NMR parameters, in particular the conformational shift analysis of Calpha protons and the coupling constants, (3)J(HNalpha), correlate well and the population estimates are in reasonable agreement among the different methods used. In the most structured peptides, we find an extension of the beta-hairpin structure comprising the two extra residues. Analysis of the pH and salt dependence shows that ionic pairs contribute to beta-hairpin stability. The interaction is electrostatic in nature and can be screened by salt. There is also an important salt-independent contribution of negatively charged groups to the stability of this family of beta-hairpin peptides.

Computational estimation of specific side chain interaction energies in alpha helices

We have used a structure energy-based computer program developed for protein design, Perla, to provide theoretical estimates of all specific side chain-side chain interaction energies occurring in alpha helices. The computed side chain-side chain interaction energies were used as substitutes for the corresponding values used by the helix/coil transition algorithm, AGADIR. Predictions of peptide helical contents were nearly as successful as those obtained with the originally calibrated set of parameters; a correlation to experimentally observed alpha-helical populations of 0.91 proved that our theoretical estimates are reasonably correct for amino acid pairs that are frequent in our database of peptides. Furthermore, we have determined experimentally the previously uncharacterized interaction energies for Lys-Ile, Thr-Ile, and Phe-Ile amino acid pairs at i,i + 4 positions. The experimental values compare favorably with the computed theoretical estimates. Importantly, the computed values for Thr-Ile and Phe-Ile interactions are better than the energies based on chemical similarity, whereas for Lys-Ile they are similar. Thus, computational techniques can be used to provide precise energies for amino acid pairwise interactions, a fact that supports the development of structure energy-based computational tools for structure predictions and sequence design.

Conformational exchange on the microsecond time scale in alpha-helix and beta-hairpin peptides measured by 13C NMR transverse relaxation

13C-NMR relaxation experiments (T(1), T(2), T(1)(rho), and NOE) were performed on selectively enriched residues in two peptides, one hydrophobic staple alpha-helix-forming peptide GFSKAELAKARAAKRGGY and one beta-hairpin-forming peptide RGITVNGKTYGR, in water and in water/trifluoroethanol (TFE). Exchange contributions, R(ex), to spin-spin relaxation rates for (13)C(alpha) and (13)C(beta) groups were derived and were ascribed to be mainly due to peptide folding-unfolding. To evaluate the exchange time, tau(ex), from R(ex), the chemical shift difference between folded and unfolded states, Deltadelta, and the populations of these states, p(i), were determined from the temperature dependence of (13)C chemical shifts. For both peptides, values for tau(ex) fell in the 1 micros to 10 micros range. Under conditions where the peptides are most folded (water/TFE, 5 degrees C), tau(ex) values for all residues in each respective peptide were essentially the same, supporting the presence of a global folding-unfolding exchange process. Rounded-up average tau(ex) values were 4 micros for the helix peptide and 9 micros for the hairpin peptide. This 2-3-fold difference in exchange times between helix and hairpin peptides is consistent with that observed for folding-unfolding of other small peptides.

Atomic resolution structure of a mutant of the spectrin SH3 domain

The crystal structure of an alpha-spectrin Src-homology 3 (SH3) domain mutant has been refined at 1.12 A resolution. This X-ray structure is at the highest resolution achieved so far for an SH3 domain. The structure allows the identification of a complete set of specific interactions and is useful for the elucidation of relations between structure and pH-dependent thermodynamic stability in a series of SH3 domain mutants.

Protein design based on folding models

The basic rules governing the folding of small, single-domain proteins are being discovered. New algorithms that can predict the major features of the folding process give the opportunity to design and optimise protein folding in a rational way. Recent experimental works suggest that sequence-specific features should be integrated in folding models to improve their performance.

Improving the refolding yield of interleukin-4 through the optimization of local interactions

Interleukin-4 (IL-4) is a multifunctional cytokine that plays an important role in the regulation of various immune responses. However, the development of IL-4 or IL-4 variants into potential therapeutic drugs is hindered by the low efficiency of the in vitro refolding process of this protein. In this work, we have investigated the improvement of the refolding yield of IL-4 using two different rational design approaches. The first one is based on the so-called inverse hydrophobic effect and involved the replacement of a solvent exposed, non-conserved, hydrophobic residue (W91) by serine. This led to an increase in stability of 1.4 kcal mol(-1) and shifted the midpoint transition temperature (Tm) from 62 to 70 degrees C. The second approach is based on the stabilization of alpha-helices through the introduction of favorable local interactions. This strategy resulted in the following helix sequence for helix C of IL-4, 68ASAAEANRHKQLIRFLKRLDRNLWGLAG95. The mutant protein was stabilized by 0.5 kcal mol(-1), the Tm shifted to 68 degrees C, and a two-fold increase in the refolding yield was consistently observed. Our results make the large-scale production of IL-4 derivatives economically more viable, suggest that a similar approach can be applied to other related proteins, and may represent a general strategy to improve in vitro refolding yields through the selective optimization of the stability of alpha-helices.

The SH3-fold family: experimental evidence and prediction of variations in the folding pathways

To investigate the relationships between protein topology, amino acid sequence and folding mechanisms, the folding transition state of the Sso7d protein has been characterised both experimentally and theoretically. Although Sso7d protein has a similar topology to that of the SH3 domains, the structure of its transition state is different from that of alpha-spectrin and src SH3 domains previously studied. The folding algorithm, Fold-X, including an energy function with specific sequence features, accounts for these differences and reproduces with a good agreement the set of experimental phi(double dagger-U) values obtained for the three proteins. Our analysis shows that taking into account sequence features underlying protein topology is critical for an accurate prediction of the folding process.

Kinesin subfamily UNC104 contains a FHA domain: boundaries and physicochemical characterization

By sequence analysis we show that the U104 domain found in the UNC104 subfamily of kinesins is a forkhead homology-associated domain (FHA). A combination of limited proteolysis, mass spectroscopy, and physicochemical analysis define this domain as a genuine autonomously folding domain. Our data show that the FHA domain is shorter than previously reported since the C-terminal alpha-helix is not part of its minimum core. Key amino acids postulated to recognize phosphorylated residues are conserved. These data suggest that the kinesin FHA domains are functional domains involved in protein-protein interactions regulated by phosphorylation.

Thermodynamic and structural characterization of Asn and Ala residues in the disallowed II' region of the Ramachandran plot

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

Internal motional amplitudes and correlated bond rotations in an alpha-helical peptide derived from 13C and 15N NMR relaxation

Peptide GFSKAELAKARAAKRGGY folds in an alpha-helical conformation that is stabilized by formation of a hydrophobic staple motif and an N-terminal capping box (Munoz V. Blanco FJ, Serrano L, 1995, Struct Biol 2:380-385). To investigate backbone and side-chain internal motions within the helix and hydrophobic staple, residues F2, A5, L7, A8, and A10 were selectively 13C- and 15N-enriched and NMR relaxation experiments were performed in water and in water/trifluoroethanol (TFE) solution at four Larmor frequencies (62.5, 125, 150, and 200 MHz for 13C). Relaxation data were analyzed using the model free approach and an anisotropic diffusion model. In water, angular variances of motional vectors range from 10 to 20 degrees and backbone phi,psi bond rotations for helix residues A5, L7, A8, and A10 are correlated indicating the presence of Calpha-H, Calpha-Cbeta, and N-H rocking-type motions along the helix dipole axis. L7 side-chain CbetaH2 and CgammaH motions are also correlated and as motionally restricted as backbone CalphaH, suggesting considerable steric hindrance with neighboring groups. In TFE which stabilizes the fold, internal motional amplitudes are attenuated and rotational correlations are increased. For the side chain of hydrophobic staple residue F2, wobbling-in-a-cone type motions dominate in water, whereas in TFE, the Cbeta-Cgamma bond and phenyl ring fluctuate more simply about the Calpha-Cbeta bond. These data support the Daragan-Mayo model of correlated bond rotations (Daragan VA, Mayo KH, 1996, J Phys Chem 100:8378-8388) and contribute to a general understanding of internal motions in peptides and proteins.

Towards understanding a molecular switch mechanism: thermodynamic and crystallographic studies of the signal transduction protein CheY

The signal transduction protein CheY displays an alpha/beta-parallel polypeptide folding, including a highly unstable helix alpha4 and a strongly charged active site. Helix alpha4 has been shown to adopt various positions and conformations in different crystal structures, suggesting that it is a mobile segment. Furthermore, the instability of this helix is believed to have functional significance because it is involved in protein-protein contacts with the transmitter protein kinase CheA, the target protein FliM and the phosphatase CheZ. The active site of CheY comprises a cluster of three aspartic acid residues and a lysine residue, all of which participate in the binding of the Mg(2+) needed for the protein activation. Two steps were followed to study the activation mechanism of CheY upon phosphorylation: first, we independently substituted the three aspartic acid residues in the active site with alanine; second, several mutations were designed in helix alpha 4, both to increase its level of stability and to improve its packing against the protein core. The structural and thermodynamic analysis of these mutant proteins provides further evidence of the connection between the active-site area and helix alpha 4, and helps to understand how small movements at the active site are transmitted and amplified to the protein surface.

Thermodynamic analysis of helix-engineered forms of the activation domain of human procarboxypeptidase A2

Thermodynamic characterization of the activation domain of human procarboxypeptidase A2, ADA2h, and its helix-engineered mutants was carried out by differential scanning calorimetry. The mutants were engineered by changing residues in the exposed face of the two alpha helices in order to increase their stability. At neutral and alkaline pH the three mutants, alpha-helix 1 (M1), alpha-helix 2 (M2) and alpha-helix 1 and alpha-helix 2 (DM), were more stable than the wild-type domain, in the order DM, M2, M1 and wild-type. Under these conditions the CD and NMR spectra of all the variants are very similar, indicating that this increase in stability is not the result of gross structural changes. Calorimetric analysis shows that the stabilizing effect of mutating the water-exposed surfaces of the helices seems to be mainly entropic, because the mutations do not change the enthalpy or the increase in heat capacity of denaturation. The unfolding behavior of all variants changes under acidic conditions: whereas wild-type and M1 have a strong tendency to aggregate, giving rise to a beta conformation upon unfolding, M2 and DM unfold reversibly, M2 being more stable than DM. CD and NMR experiments at pH 3.0 suggest that a region involving residues of the second and third beta strands as well as part of alpha-helix 1 changes its conformation. It seems that the enhanced stability of the altered conformation of M2 and DM reduces the aggregation tendency of ADA2h at acidic pH.

Protein engineering as a strategy to avoid formation of amyloid fibrils

The activation domain of human procarboxypeptidase A2 (ADA2h) aggregates following thermal or chemical denaturation at acidic pH. The aggregated material contains well-defined ordered structures with all the characteristics of the fibrils associated with amyloidotic diseases. Variants of ADA2h containing a series of mutations designed to increase the local stability of each of the two helical regions of the protein have been found to have a substantially reduced propensity to form fibrils. This arises from a reduced tendency of the denatured species to aggregate rather than from a change in the overall stability of the native state. The reduction in aggregation propensity may result from an increase in the stability of local relative to longer range interactions within the polypeptide chain. These findings show that the intrinsic ability of a protein to form amyloid can be altered substantially by protein engineering methods without perturbing significantly its overall stability or activity. This suggests new strategies for combating diseases associated with the formation of aggregated proteins and for the design of novel protein or peptide therapeutics.

Evaluation of a coagulation analyzer

The CA-6000 compared favorably to the well-established 1600C in terms of clinical utility. Its enhanced reagent capacity and open nature make it a wise instrument selection option.

Stabilisation of alpha-helices by site-directed mutagenesis reveals the importance of secondary structure in the transition state for acylphosphatase folding

The effects of stabilising mutations on the folding process of common-type acylphosphatase have been investigated. The mutations were designed to increase the helical propensity of the regions of the polypeptide chain corresponding to the two alpha-helices of the native protein. Various synthetic peptides incorporating the designed mutations were produced and their helical content estimated by circular dichroism. The most substantial increase in helical content is found for the peptide carrying five mutations in the second alpha-helix. Acylphosphatase variants containing the corresponding mutations display, to different extents, enhanced conformational stabilities as indicated by equilibrium urea denaturation experiments monitored by changes of intrinsic fluorescence. All the protein variants studied here refold with apparent two-state kinetics. Mutations in the first alpha-helix are responsible for a small increase in the refolding rate, accompanied by a marked decrease in the unfolding rate. On the other hand, multiple mutations in the second helix result in a considerable increase in the refolding rate without any significant effect on the unfolding rate. Addition of trifluoroethanol was found to accelerate the folding of the acylphosphatase variants, the extent of the acceleration being inversely proportional to the intrinsic rate of folding of the corresponding mutant. The trifluoroethanol-induced acceleration is far less marked for those variants whose alpha-helical structure is efficiently stabilised by amino acid replacements. This observation suggests that trifluoroethanol acts in a similar manner to the stabilising mutations in promoting native-like secondary structure. Analysis of the kinetic data indicates that the second helix is fully consolidated in the transition state for folding of acylphosphatase, whereas the first helix is only partially formed. These data suggest that the second helix is an important element in the folding process of the protein.

Engineering stability in gene networks by autoregulation

The genetic and biochemical networks which underlie such things as homeostasis in metabolism and the developmental programs of living cells, must withstand considerable variations and random perturbations of biochemical parameters. These occur as transient changes in, for example, transcription, translation, and RNA and protein degradation. The intensity and duration of these perturbations differ between cells in a population. The unique state of cells, and thus the diversity in a population, is owing to the different environmental stimuli the individual cells experience and the inherent stochastic nature of biochemical processes (for example, refs 5 and 6). It has been proposed, but not demonstrated, that autoregulatory, negative feedback loops in gene circuits provide stability, thereby limiting the range over which the concentrations of network components fluctuate. Here we have designed and constructed simple gene circuits consisting of a regulator and transcriptional repressor modules in Escherichia coli and we show the gain of stability produced by negative feedback.

Similarities between the spectrin SH3 domain denatured state and its folding transition state

We have expanded our description of the energy landscape for folding of the SH3 domain of chicken alpha-spectrin by a detailed structural characterization of its denatured state ensemble (DSE). This DSE is significantly populated under mildly acidic conditions in equilibrium with the folded state. Evidence from heteronuclear nuclear magnetic resonance (NMR) experiments on (2)H, (15)N-labeled protein suggests the presence of conformers whose residual structure bears some resemblence to the structure of the folding transition state of this protein. NMR analysis in a mutant with an engineered, non-native alpha-helical tendency shows a significant amount of local non-native structure in the mutant, while the overall characteristics of the DSE are unchanged. Comparison with recent theoretical predictions of SH3 domain folding reactions reveals an interesting correlation with the predicted early events. Based on these results and recent data from other systems, we propose that the DSE of a protein will resemble the intermediate or transition state of its nearest rate-limiting step, as a consequence of simple energetic and kinetic principles.

Model for stathmin/OP18 binding to tubulin

Stathmin/OP18 is a regulatory phosphoprotein that controls microtubule (MT) dynamics. The protein does not have a defined three-dimensional structure, although it contains three distinct regions (an unstructured N-terminus, N: 1-44; a region with high helix propensity, H 1: 44-89; and a region with low helix propensity, H 2: 90-142). The full protein and a combination of H 1 and H 2 inhibits tubulin polymerization, while the combination of H 1 and the N-terminus is less efficient. None of the individual three regions alone are functional in this respect. However, all of them cross-link to alpha-tubulin, but only full-length stathmin produces high-molecular-weight products. Mass spectrometry analysis of alpha-tubulin-stathmin/OP18 and its truncation products shows that full-length stathmin/OP18 binds to the region around helix 10 of alpha-tubulin, a region that is involved in longitudinal interactions in the MT, sequestering the dimer and possibly linking two tubulin heterodimers. In the absence of the N-terminus, stathmin/OP18 binds to only one molecule of alpha-tubulin, at the top of the free tubulin heterodimer, preventing polymerization.

The folding transition state between SH3 domains is conformationally restricted and evolutionarily conserved

The protein engineering analysis of the alpha-spectrin SH3 domain at three different stability conditions (pH 7.0, 3.5 and 2.5) reveals a folding transition state structured around the distal loop beta-hairpin and the 310-helix. This region is impervious to overall changes in protein stability, suggesting a transition state ensemble with little conformational variability. Comparison with the Src SH3 domain (36% sequence homology) indicates that the transition state in this protein family may be conserved. Discrepancies at some positions can be rationalized in terms of the different interactions made by the different side chains in both domains. Brønsted plot analysis confirms the straight phi(doubledagger-U) results and shows two folding subdomains for this small protein. These results, together with previous data on circular permutants of the alpha-spectrin SH3 domain, indicate that polypeptide topology and chain connectivity play a major role in the folding reaction of this protein family.

Position dependence of amino acid intrinsic helical propensities II: non-charged polar residues: Ser, Thr, Asn, and Gln

The assumption that the intrinsic alpha-helical propensities of the amino acids are position independent was critical in several helix/coil transition theories. In the first paper of these series, we reported that this is not the case for Gly and nonpolar aliphatic amino acids (Val, Leu, Met, and Ile). Here we have analyzed the helical intrinsic propensities of noncharged polar residues (Ser, Thr, Asn, and Gln) at different positions of a model polyalanine-based peptide. We found that Thr is more favorable (by approximately 0.3 kcal/mol) at positions N1 and N2 than in the helix center, although for Ser, Asn, and Gln the differences are smaller (+/-0.2 kcal/mol), and in many cases within the experimental error. There is a reasonable agreement (+/-0.2 kcal/mol) between the calculated free energies, using the ECEPP/2 force field equipped with a hydration potential, and the experimental data, except at position N1.

Local water bridges and protein conformational stability

Recent studies have pointed out the important role of local water structures in protein conformational stability. Here, we present an accurate and computationally effective way to estimate the free energy contribution of the simplest water structure motif--the water bridge. Based on the combination of empirical parameters for accessible protein surface area and the explicit consideration of all possible water bridges with the protein, we introduce an improved protein solvation model. We find that accounting for water bridge formation in our model is essential to understand the conformational behavior of polypeptides in water. The model formulation, in fact, does not depend on the polypeptide nature of the solute and is therefore applicable to other flexible biomolecules (i.e., DNAs, RNAs, polysaccharides, etc.).

A tale of two secondary structure elements: when a beta-hairpin becomes an alpha-helix

In this work, we have analyzed the relative importance of secondary versus tertiary interactions in stabilizing and guiding protein folding. For this purpose, we have designed four different mutants to replace the alpha-helix of the GB1 domain by a sequence with strong beta-hairpin propensity in isolation. In particular, we have chosen the sequence of the second beta-hairpin of the GB1 domain, which populates the native conformation in aqueous solution to a significant extent. The resulting protein has roughly 30 % of its sequence duplicated and maintains the 3D-structure of the wild-type protein, but with lower stability (up to -5 kcal/mol). The loss of intrinsic helix stability accounts for about 80 % of the decrease in free energy, illustrating the importance of local interactions in protein stability. Interestingly enough, all the mutant proteins, included the one with the duplicated beta-hairpin sequence, fold with similar rates as the GB1 domain. Essentially, it is the nature of the rate-limiting step in the folding reaction that determines whether a particular interaction will speed up, or not, the folding rates. While local contacts are important in determining protein stability, residues involved in tertiary contacts in combination with the topology of the native fold, seem to be responsible for the specificity of protein structures. Proteins with non-native secondary structure tendencies can adopt stable folds and be as efficient in folding as those proteins with native-like propensities.

Stable proline box motif at the N-terminal end of alpha-helices

We describe a novel N-terminal alpha-helix local motif that involves three hydrophobic residues and a Pro residue (Pro-box motif). Database analysis shows that when Pro is the N-cap of an alpha-helix the distribution of amino acids in adjacent positions changes dramatically with respect to the average distribution in an alpha-helix, but not when Pro is at position N1. N-cap Pro residues are usually associated to Ile and Leu, at position N', Val at position N3 and a hydrophobic residue (h) at position N4. The side chain of the N-cap Pro packs against Val, while the hydrophobic residues at positions N' and N4 make favorable interactions. To analyze the role of this putative motif (sequence fingerprint hPXXhh), we have synthesized a series of peptides and analyzed them by circular dichroism (CD) and NMR. We find that this motif is formed in peptides, and that the accompanying hydrophobic interactions contribute up to 1.2 kcal/mol to helix stability. The fact that some of the residues in this fingerprint are not good N-cap and helix formers results in a small overall stabilization of the alpha-helix with respect to other peptides having Gly as the N-cap and Ala at N3 and N4. This suggests that the Pro-box motif will not specially contribute to protein stability but to the specificity of its fold. In fact, 80% of the sequences that contain the fingerprint sequence in the protein database are adopting the described structural motif, and in none of them is the helix extended to place Pro at the more favorable N1 position.