Designed molecules that fold to mimic protein secondary structures

Molecules that fold to mimic protein secondary structures have emerged as important targets of bioorganic chemistry. Recently, a variety of compounds that mimic helices, turns, and sheets have been developed, with notable advances in the design of beta-peptides that mimic each of these structures. These compounds hold promise as a step toward synthetic molecules with protein-like properties and as drugs that block protein-protein interactions.     

Linkers for improved cleavage of fusion proteins with an engineered alpha-lytic protease.

Addition of an N-terminal fusion partner can greatly aid the expression and purification of a recombinant protein in Escherichia coli. We investigated two genetically engineered proteases designed to remove the fusion partner after the protein of interest has been expressed. Recombinant human insulin-like growth factor-II (hIGF-II) has been produced from E. coli-derived fusion proteins using a novel enzymatic cleavage system that uses a mutant of alpha-lytic protease. Initially, two potential fusion protein linkers were designed, Pro-Ala-Pro-His (PAPH) and Pro-Ala-Pro-Met (PAPM), and were tested as substrates in the form of synthetic dodecapeptides. Using mass spectrometry and reverse-phase HPLC, the position of cleavage was confirmed and the kinetics of synthetic peptide cleavage were examined. Use of the linkers in hIGF-II fusion proteins produced in E. coli was then evaluated. The fusion proteins constructed consist of the first 11 amino acids of porcine growth hormone linked N-terminally to hIGF-II by six amino acids that include the dipeptide Val-Asn followed by a variable tetrapeptide protease cleavage motif. Mass spectrometry and N-terminal sequencing confirmed that proteolytic cleavage of the fusion proteins had occurred at the predicted sites. Using the fusion proteins as substrates, the cleavage of the rationally designed motifs by the alpha-lytic protease mutant was compared. The fusion protein containing the motif PAPM had a k(cat)/K(M) ratio indicating a 1.6-fold preference over the PAPH fusion protein for cleavage by this enzyme. Furthermore, when hIGF-II fusion proteins containing the designed cleavable linkers were processed with the engineered alpha-lytic protease, they gave greatly improved yields of native hIGF-II compared to an analogous fusion protein cleaved by H64A subtilisin. Comparison of the peptide and protein cleavage studies shows that the efficient proteolysis of the cleavage motifs is an inherent property of the designed sequences and is not determined by secondary or tertiary structure in the fusion proteins.     

Analysis of interactive packing of secondary structural elements in alpha/beta units in proteins

An alpha-helix and a beta-strand are said to be interactively packed if at least one residue in each of the secondary structural elements loses 10% of its solvent accessible contact area on association with the other secondary structural element. An analysis of all such 5,975 nonidentical alpha/beta units in protein structures, defined at < or = 2.5 A resolution, shows that the interaxial distance between the alpha-helix and the beta-strand is linearly correlated with the residue-dependent function, log[(V/nda)/n-int], where V is the volume of amino acid residues in the packing interface, nda is the normalized difference in solvent accessible contact area of the residues in packed and unpacked secondary structural elements, and n-int is the number of residues in the packing interface. The beta-sheet unit (beta u), defined as a pair of adjacent parallel or antiparallel hydrogen-bonded beta-strands, packing with an alpha-helix shows a better correlation between the interaxial distance and log(V/nda) for the residues in the packing interface. This packing relationship is shown to be useful in the prediction of interaxial distances in alpha/beta units using the interacting residue information of equivalent alpha/beta units of homologous proteins. It is, therefore, of value in comparative modeling of protein structures.     

Thermal folding and mechanical unfolding pathways of protein secondary structures

Mechanical stretching of secondary structures is studied through molecular dynamics simulations of a Go-like model. Force versus displacement curves are studied as a function of the stiffness and velocity of the pulling device. The succession of stretching events, as measured by the order in which contacts are ruptured, is compared to the sequencing of events during thermal folding and unfolding. Opposite cross-correlations are found for an alpha-helix and a beta-hairpin structure. In a tandem of two alpha-helices, the two constituent helices unravel nearly simultaneously. A simple condition for simultaneous versus sequential unraveling of repeat units is presented.     

A helix-turn-strand structural motif common in alpha-beta proteins.

By exhaustive structural comparisons, we have found that about one-third of the alpha-helix-turn-beta-strand polypeptides in alpha-beta barrel domains share a common structural motif. The chief characteristics of this motif are that first, the geometry of the turn between the alpha-helix and the beta-strand is somewhat constrained, and second, the beta-strand contains a hydrophobic patch that fits into a hydrophobic pocket on the alpha-helix. The geometry of the turn does not seem to be a major determinant of the alpha-beta unit, because the turns vary in length from four to six residues. However, the motif does not occur when there are few constraints on the geometry of the turn-for instance, when the turns between the alpha-helix and the beta-strands are very long. It also occurs much less frequently in flat-sheet alpha-beta proteins, where the topology is much less regular and the amount of twist on the sheet varies considerably more than in the barrel proteins. The motif may be one of the basic building blocks from which alpha-beta barrels are constructed.     

Energetics of the structure and chain tilting of antiparallel beta-barrels in proteins

The preferred structural pattern of antiparallel beta-barrels in proteins, described as the right-handed tilting of the peptide strands with respect to the axis of the barrel, is accounted for in terms of intra- and interchain interaction energies. It is related to the preference of beta-sheets for right-handed twisting. Conformational energy computations have been carried out on three eight-stranded antiparallel beta-barrels composed of six-residue strands, in which L-Val and Gly alternate, and having a right-handed, a left-handed, or no tilt. After energy minimization, the relative energies of these structures were 0.0, 8.6, and 46.1 kcal/mol, respectively; i.e., the right-tilted beta-barrel is favored energetically, in agreement with anti-parallel beta-barrels observed in proteins. Tilting of the barrel is favored, relative to the nontilted structure, by both intra- and interstrand interactions, because tilting allows better packing of the bulky side chains. On the other hand, the energy difference between the left- and right-tilted barrels arises essentially from intrachain interactions. This is a consequence of the preference of beta-sheets for a right-handed twist. Space limitations inside the barrel are satisfied if there is an alternation of bulky residues and residues with small or no side chain (preferably Gly) in neighboring positions on adjacent strands. Such a pattern is seen frequently in antiparallel beta-barrels of globular proteins. The computations indicate that a structure with Val...Gly pairs can be accommodated in a beta-barrel with no distortion.     

Influence of local and residual structures on the scaling behavior and dimensions of unfolded proteins

Although recent spectroscopic studies of chemically denatured proteins hint at significant nonrandom residual structure, the results of extensive small angle X-ray scattering studies suggest random coil behavior, calling for a coherent understanding of these seemingly contradicting observations. Here, we report the results of a Monte Carlo study of the effects of two types of local structures, alpha helix and Polyproline II (PPII) helix, on the dimensions of random coil polyalanine chains viewed as a model of highly denatured proteins. We find that although Flory's power law scaling, long regarded as a signature of random coil behavior, holds for chains containing up to 90% alpha or PPII helix, the absolute magnitude of the chain dimensions is sensitive to helix content. As residual alpha helix content increases, the chain contracts until it reaches a minimum radius at approximately 70% helix, after which the chain dimensions expand rapidly. With an alpha helix content of approximately 20%, corresponding to the Ramachandran probability of being in the helical basin, experimentally observed radii of gyration are recovered. Experimental radii are similarly recovered at an alpha helix content of approximately 87%, providing an explanation for the previously puzzling experimental finding that the dimensions of the highly helical methanol-induced unfolded state are experimentally indistinguishable from those of the helix-poor urea-unfolded state. In contrast, the radius of gyration increases monotonically with increasing PPII content, and is always more expanded than the dimensions observed experimentally. These results suggest that PPII is unlikely the sole, dominant preferred conformation for unfolded proteins.     

Identification,classification, and analysis of beta-bulges in proteins.

A beta-bulge is a region of irregularity in a beta-sheet involving two beta-strands. It usually involves two or more residues in the bulged strand opposite to a single residue on the adjacent strand. These irregularities in beta-sheets were identified and classified automatically, extending the definition of beta-bulges given by Richardson et al. (Richardson, J.S., Getzoff, E.D., & Richardson, D.C., 1978, Proc. Natl. Acad. Sci. USA 75, 2574-2578). A set of 182 protein chains (170 proteins) was used, and a total of 362 bulges were extracted. Five types of beta-bulges were found: classic, G1, wide, bent, and special. Their characteristic amino acid preferences were found for most classes of bulges. Basically, bulges occur frequently in proteins; on average there are more than two bulges per protein. In general, beta-bulges produce two main changes in the structure of a beta-sheet: (1) disrupt the normal alternation of side-chain direction; (2) accentuate the twist of the sheet, altering the direction of the surrounding strands.     

Fishing new proteins in the twilight zone of genomes: the test case of outer membrane proteins in Escherichia coli K12, Escherichia coli O157:H7, and other Gram-negative bacteria

We address the problem of clustering the whole protein content of genomes into three different categories-globular, all-alpha, and all-beta membrane proteins-with the aim of fishing new membrane proteins in the pool of nonannotated proteins (twilight zone). The focus is then mainly on outer membrane proteins. This is performed by using an integrated suite of programs (Hunter) specifically developed for predicting the occurrence of signal peptides in proteins of Gram-negative bacteria and the topography of all-alpha and all-beta membrane proteins. Hunter is tested on the well and partially annotated proteins (2160 and 760, respectively) of Escherichia coli K 12 scoring as high as 95.6% in the correct assignment of each chain to the category. Of the remaining 1253 nonannotated sequences, 1099 are predicted globular, 136 are all-alpha, and 18 are all-beta membrane proteins. In Escherichia coli 0157:H7 we filtered 1901 nonannotated proteins. Our analysis classifies 1564 globular chains, 327 inner membrane proteins, and 10 outer membrane proteins. With Hunter, new membrane proteins are added to the list of putative membrane proteins of Gram-negative bacteria. The content of outer membrane proteins per genome (nine are analyzed) ranges from 1.5% to 2.4%, and it is one order of magnitude lower than that of inner membrane proteins. The finding is particularly relevant when it is considered that this is the first large-scale analysis based on validated tools that can predict the content of outer membrane proteins in a genome and can allow cross-comparison of the same protein type between different species.     

Application of a directed conformational search for generating 3-D coordinates for protein structures from alpha-carbon coordinates.

A directed conformational search algorithm using the program CONGEN (ref. 3), which samples backbone conformers, is described. The search technique uses information from the partially built structures to direct the search process and is tested on the problem of generating a full set of backbone Cartesian coordinates given only alpha-carbon coordinates. The method has been tested on six proteins of known structure, varying in size and classification, and was able to generate the original backbone coordinates with RMSs ranging from 0.30-0.87A for the alpha-carbons and 0.5-0.99A RMSs for the backbone atoms. Cis peptide linkages were also correctly identified. The procedure was also applied to two proteins available with only alpha-carbon coordinates in the Brookhaven Protein Data Bank; thioredoxin (SRX) and triacylglycerol acylhydrolase (TGL). All-atom models are proposed for the backbone of both these proteins. In addition, the technique was applied to randomized coordinates of flavodoxin to assess the effects of irregularities in the data on the final RMS. This study represents the first time a deterministic conformational search was used on such a large scale.     

Effects of detergents on the secondary structures of prion protein peptides as studied by CD spectroscopy.

Pathogenic prion proteins (PrP(Sc)) are thought to be produced by alpha-helical to beta-sheet conformational changes in the normal cellular prion proteins (PrP(C)) located solely in the caveolar compartments. In order to inquire into the possible conformational changes due to the influences of hydrophobic environments within caveolae, the secondary structures of prion protein peptides were studied in various kinds of detergents by CD spectra. The peptides studied were PrP(129-154) and PrP(192-213); the former is supposed to assume beta-sheets and the latter alpha-helices, in PrP(Sc). The secondary structure analyses for the CD spectra revealed that in buffer solutions, both PrP(129-154) and PrP(192-213) mainly adopted random-coils (approximately 60%), followed by beta-sheets (30%-40%). PrP(129-154) showed no changes in the secondary structures even in various kinds of detergents such as octyl-beta-D-glucopyranoside (OG), octy-beta-D-maltopyranoside (OM). sodium dodecyl sulfate (SDS), Zwittergent 3-14 (ZW) and dodecylphosphocholine (DPC). In contrast, PrP(192-213) changed its secondary structure depending on the concentration of the detergents. SDS, ZW, OG and OM increased the alpha-helical content, and decreased the beta-sheet and random-coil contents. DPC also increased the alpha-helical content, but to a lesser extent than did SDS, ZW, OG or OM. These results indicate that PrP(129-154) has a propensity to adopt predominantly beta-sheets. On the other hand, PrP(192-213) has a rather fickle propensity and varies its secondary structure depending on the environmental conditions. It is considered that the hydrophobic environments provided by these detergents may mimic those provided by gangliosides in caveolae, the head groups of which consist of oligosaccharide chains containing sialic acids. It is concluded that PrP(C) could be converted into a nascent PrP(Sc) having a transient PrP(Sc) like structureunder the hydrophobic environments produced by gangliosides.     

Positional preference of proline in alpha-helices.

Conformational free energy calculations have been carried out for proline-containing alanine-based pentadecapeptides with the sequence Ac-(Ala)n-Pro-(Ala)m-NHMe, where n + m = 14, to figure out the positional preference of proline in alpha-helices. The relative free energy of each peptide was calculated by subtracting the free energy of the extended conformation from that of the alpha-helical one, which is used here as a measure of preference. The highest propensity is found for the peptide with proline at the N-terminus (i.e., Ncap + 1 position), and the next propensities are found at Ncap, N' (Ncap - 1), and C' (Ccap + 1) positions. These computed results are reasonably consistent with the positional propensities estimated from X-ray structures of proteins. The breaking in hydrogen bonds around proline is found to play a role in destabilizing alpha-helical conformations, which, however, provides the favored hydration of the corresponding N-H and C=O groups. The highest preference of proline at the beginning of alpha-helix appears to be due to the favored electrostatic and nonbonded energies between two residues preceding proline and the intrinsic stability of alpha-helical conformation of the proline residue itself as well as no disturbance in hydrogen bonds of alpha-helix by proline. The average free energy change for the substitution of Ala by Pro in a alpha-helix is computed to be 4.6 kcal/mol, which is in good agreement with the experimental value of approximately 4 kcal/mol estimated for an oligopeptide dimer and proteins of barnase and T4 lysozyme.     

Recurrent alpha beta loop structures in TIM barrel motifs show a distinct pattern of conserved structural features.

A systematic survey of seven parallel alpha/beta barrel protein domains, based on exhaustive structural comparisons, reveals that a sizable proportion of the alpha beta loops in these proteins--20 out of a total of 49--belong to either one of two loop types previously described by Thornton and co-workers. Six loops are of the alpha beta 1 type, with one residue between the alpha-helix and beta-strand, and 13 are of the alpha beta 3 type, with three residues between the helix and the strand. Protein fragments embedding the identified loops, and termed alpha beta connections since they contain parts of the flanking helix and strand, have been analyzed in detail revealing that each type of connection has a distinct set of conserved structural features. The orientation of the beta-strand relative to the helix and loop portions is different owing to a very localized difference in backbone conformation. In alpha beta 1 connections, the chain enters the beta-strand via a residue adopting an extended conformation, while in alpha beta 3 it does so via a residue in a near alpha-helical conformation. Other conserved structural features include distinct patterns of side chain orientation relative to the beta-sheet surface and of main chain H-bonds in the loop and the beta-strand moieties. Significant differences also occur in packing interactions of conserved hydrophobic residues situated in the last turn of the helix. Yet the alpha-helix surface of both types of connections adopts similar orientations relative to the barrel sheet surface. Our results suggest furthermore that conserved hydrophobic residues along the sequence of the connections, may be correlated more with specific patterns of interactions made with neighboring helices and sheet strands than with helix/strand packing within the connection itself. A number of intriguing observations are also made on the distribution of the identified alpha beta 1 and alpha beta 3 loops within the alpha/beta-barrel motifs. They often occur adjacent to each other; alpha beta 3 loops invariably involve even numbered beta-strands, while alpha beta 1 loops involve preferentially odd beta-strands; all the analyzed proteins contain at least one alpha beta 3 loop in the first half of the eightfold alpha/beta barrel. Possible origins of all these observations, and their relevance to the stability and folding of parallel alpha/beta barrel motifs are discussed.     

Stabilization of discordant helices in amyloid fibril-forming proteins

Several proteins and peptides that can convert from alpha-helical to beta-sheet conformation and form amyloid fibrils, including the amyloid beta-peptide (Abeta) and the prion protein, contain a discordant alpha-helix that is composed of residues that strongly favor beta-strand formation. In their native states, 37 of 38 discordant helices are now found to interact with other protein segments or with lipid membranes, but Abeta apparently lacks such interactions. The helical propensity of the Abeta discordant region (K16LVFFAED23) is increased by introducing V18A/F19A/F20A replacements, and this is associated with reduced fibril formation. Addition of the tripeptide KAD or phospho-L-serine likewise increases the alpha-helical content of Abeta(12-28) and reduces aggregation and fibril formation of Abeta(1-40), Abeta(12-28), Abeta(12-24), and Abeta(14-23). In contrast, tripeptides with all-neutral, all-acidic or all-basic side chains, as well as phosphoethanolamine, phosphocholine, and phosphoglycerol have no significant effects on Abeta secondary structure or fibril formation. These data suggest that in free Abeta, the discordant alpha-helix lacks stabilizing interactions (likely as a consequence of proteolytic removal from a membrane-associated precursor protein) and that stabilization of this helix can reduce fibril formation.     

The noncatalytic triad of alpha-amylases: a novel structural motif involved in conformational stability

Chloride-activated alpha-amylases contain a noncatalytic triad, independent of the glycosidic active site, perfectly mimicking the catalytic triad of serine-proteases and of other active serine hydrolytic enzymes. Mutagenesis of Glu, His, and Ser residues in various alpha-amylases shows that this pattern is a structural determinant of the enzyme conformation that cannot be altered without losing the intrinsic stability of the protein. (1)H-(15)N NMR spectra of a bacterial alpha-amylase reveal proton signals that are identical with the NMR signature of catalytic triads and especially a deshielded proton involving a protonated histidine and displaying properties similar to that of a low barrier hydrogen bond. It is proposed that the H-bond between His and Glu of the noncatalytic triad is an unusually strong interaction, responsible for the observed NMR signal and for the weak stability of the triad mutants. Furthermore, a stringent template-based search of the Protein Data Bank demonstrated that this motif is not restricted to alpha-amylases, but is also found in 80 structures from 33 different proteins, amongst which SH2 domain-containing proteins are the best representatives.     

Mutational analysis of the active site of indoleglycerol phosphate synthase from Escherichia coli.

Indoleglycerol phosphate synthase catalyzes the ring closure of 1-(2-carboxyphenylamino)-1-deoxyribulose 5''-phosphate to indoleglycerol phosphate, the fifth step in the pathway of tryptophan biosynthesis from chorismate. Because chemical synthesis of indole derivatives from arylamino ketones requires drastic solvent conditions, it is interesting by what mechanism the enzyme catalyzes the same condensation reaction. Seven invariant polar residues in the active site of the enzyme from Escherichia coli have been mutated directly or randomly, to identify the catalytically essential ones. A strain of E. coli suitable for selecting and classifying active mutants by functional complementation was constructed by precise deletion of the trpC gene from the genome. Judged by growth rates of transformants on selective media, mutants with either S58 or S60 replaced by alanine were indistinguishable from the wild-type, but R186 replaced by alanine was still partially active. Saturation random mutagenesis of individual codons showed that E53 was partially replaceable by aspartate and cysteine, whereas K114, E163, and N184 could not be replaced by any other residue. Partially active mutant proteins were purified and their steady-state kinetic and inhibitor binding constants determined. Their relative catalytic efficiencies paralleled their relative complementation efficiencies. These results are compatible with the location of the essential residues in the active site of the enzyme and support a chemically plausible catalytic mechanism. It involves two enzyme-bound intermediates and general acid-base catalysis by K114 and E163 with the support of E53 and N184.     

Contribution of secondary structure to the heat capacity and enthalpy distribution of the unfolded state in proteins.

We have recently shown that one can construct the enthalpy distribution for protein molecules from experimental knowledge of the temperature dependence of the heat capacity. For many proteins the enthalpy distribution evaluated at the midpoint of the denaturation transition (corresponding to the maximum in the heat capacity vs temperature curve) is broad and biphasic, indicating two different populations of molecules (native and unfolded) with distinctly different enthalpies. At temperatures above the denaturation point, the heat capacity for the unfolded state in many proteins is quite large and using the analysis just mentioned, we obtain a gaussian-like enthalpy distribution that is very broad. A large value of the heat capacity indicates that there are structural changes going on in the unfolded state above the transition temperature. In the present paper we investigate the origin of this large heat capacity by considering the presence of changing amounts of secondary structure (specifically, alpha-helix) in the unfolded state. For this purpose we use the empirical estimates of the Zimm-Bragg sigma and s factors for all of the native amino acids in water as determined by Scheraga and co-workers. Using myoglobin as an example, we calculate probability profiles and distribution functions for the total number of helix states in the specific-sequence molecule. Given the partition function for the specific-sequence molecule, we can then calculate a set of enthalpy moments for the molecule from which we obtain a good estimate of the enthalpy distribution in the unfolded state. This distribution turns out to be quite narrow when compared with the distribution obtained from the raw heat capacity data. We conclude that there must be other major structural changes (backbone and solvent) that are not accounted for by the inclusion of alpha-helix in the unfolded state.     

Prediction of secondary structure by evolutionary comparison: application to the alpha subunit of tryptophan synthase

The amino acid sequences of the a subunits of tryptophan synthase from ten different microorganisms were aligned by standard procedures. The alpha helices, beta strands and turns of each sequence were predicted separately by two standard prediction algorithms and averaged at homologous sequence positions. Additional evidence for conserved secondary structure was derived from profiles of average hydropathy and chain flexibility values, leading to a joint prediction. There is good agreement between (1) predicted beta strands, maximal hydropathy and minimal flexibility, and (2) predicted loops, great chain flexibility, and protein segments that accept insertions of various lengths in individual sequences. The a subunit is predicted to have eight repeated beta-loop-alpha-loop motifs with an extra N-terminal alpha helix and an intercalated segment of highly conserved residues. This pattern suggests that the territory structure of the a subunit is an eightfold alpha/beta barrel. The distribution of conserved amino acid residues and published data on limited proteolysis, chemical modification, and mutagenesis are consistent with the alpha/beta barrel structure. Both the active site of the a subunit and the combining site for the beta 2 subunit are at the end of the barrel formed by the carboxyl-termini of the beta strands.