Protein design, a minimalist approach

Database algorithm for generating protein backbone and side-chain co-ordinates from a C alpha trace application to model building and detection of co-ordinate errors

The problem of constructing all-atom model co-ordinates of a protein from an outline of the polypeptide chain is encountered in protein structure determination by crystallography or nuclear magnetic resonance spectroscopy, in model building by homology and in protein design. Here, we present an automatic procedure for generating full protein co-ordinates (backbone and, optionally, side-chains) given the C alpha trace and amino acid sequence. To construct backbones, a protein structure database is first scanned for fragments that locally fit the chain trace according to distance criteria. A best path algorithm then sifts through these segments and selects an optimal path with minimal mismatch at fragment joints. In blind tests, using fully known protein structures, backbones (C alpha, C, N, O) can be reconstructed with a reliability of 0.4 to 0.6 A root-mean-square position deviation and not more than 0 to 5% peptide flips. This accuracy is sufficient to identify possible errors in protein co-ordinate sets. To construct full co-ordinates, side-chains are added from a library of frequently occurring rotamers using a simple and fast Monte Carlo procedure with simulated annealing. In tests on X-ray structures determined at better than 2.5 A resolution, the positions of side-chain atoms in the protein core (less than 20% relative accessibility) have an accuracy of 1.6 A (r.m.s. deviation) and 70% of chi 1 angles are within 30 degrees of the X-ray structure. The computer program MaxSprout is available on request.

Repacking protein interiors

Several goals of protein engineering may be achieved through redesign and repacking of protein interiors. The effects of interior apolar substitutions on protein stability depend strongly on the site of the substitution. One reason for this is that protein interiors have properties both of apolar liquids and of crystalline solids. Substitutions at interior sites affect the stability of a protein by changing the hydrophobicity, but each site in a protein has a characteristic energy associated with introducing packing changes, and the net stability depends on both of these factors.

A thermodynamic scale for the helix-forming tendencies of the commonly occurring amino acids

Amino acids have distinct conformational preferences that influence the stabilities of protein secondary and tertiary structures. The relative thermodynamic stabilities of each of the 20 commonly occurring amino acids in the alpha-helical versus random coil states have been determined through the design of a peptide that forms a noncovalent alpha-helical dimer, which is in equilibrium with a randomly coiled monomeric state. The alpha helices in the dimer contain a single solvent-exposed site that is surrounded by small, neutral amino acid side chains. Each of the commonly occurring amino acids was substituted into this guest site, and the resulting equilibrium constants for the monomer-dimer equilibrium were determined to provide a list of free energy difference (delta delta G degree) values.

Ion-water and ion-polypeptide correlations in a gramicidin-like channel. A molecular dynamics study

This work describes a molecular dynamics study of ion-water and ion-polypeptide correlation in a model gramicidin-like channel (the polyglycine analogue) based upon interaction between polarizable, multipolar groups. The model suggests that the vicinity of the dimer junction and of the ethanolamine tail are regions of unusual flexibility. Cs+ binds weakly in the mouth of the channel: there it coordinates five water molecules and the #11CO group with which it interacts strongly and is ideally aligned. In the channel interior it is generally pentacoordinate; at the dimer junction, because of increased channel flexibility, it again becomes essentially hexacoordinate. The ion is also strongly coupled to the #13 CO but not to either #9 or #15, consistent with 13C NMR data. Water in the channel interior is strikingly different from bulk water; it has a much lower mean dipole moment. This correlates with our observation (which differs from that of previous studies) that water-water angular correlations do not persist within the channel, a result independent of ion occupancy or ionic polarity. In agreement with streaming potential measurements, there are seven single file water molecules associated with Cs+ permeation; one of these is always in direct contact with bulk water. At the mouth of an ion-free channel, there is a pattern of dipole moment alteration among the polar groups. Due to differential interaction with water, exo-carbonyls have unusually large dipole moments whereas those of the endo-carbonyls are low. The computed potential of mean force for CS+ translocation is qualitatively reasonable. However, it only exhibits a weakly articulated binding site and it does not quantitatively account for channel energetics. Correction for membrane polarization reduces, but does not eliminate, these problems.

Studies of synthetic helical peptides using circular dichroism and nuclear magnetic resonance

We have designed a set of 17-residue synthetic peptides to be monomeric helices in aqueous solution. Circular dichrosim experiments indicate the presence of helical structure in aqueous solution at low temperature and low pH. The two-dimensional nuclear magnetic resonance results for one of the peptides show a segment of ten residues which clearly meets all of the criteria for the existence of helical structure at both 5 degrees C and 15 degrees C. The first four residues of the peptide are in a largely extended conformation. Calculations suggest that residues 5 through 14 are significantly helical at 5 degrees C. When the temperature is increased, circular dichroism spectra indicate that the helical content decreases. At 15 degrees C, the 3JN alpha coupling constants increase in the helical region, indicating an increase in motion or conformational averaging in the helical segment. None of the peptides has pH titration behavior consistent with salt bridge stabilization of helical conformation. Our data lend themselves to interpretation with the helix dipole model and specific side-chain interactions. When the N and C termini charges are removed the helical content of the peptides increases. The amount of helicity increases as the pH is lowered, due to the ionization of His16. Much of the helical stabilization appears to be due to a specific side-chain interaction between His16 and Tyr12.

Solubility as a function of protein structure and solvent components

This review deals with ways of stabilizing proteins against aggregation and with methods to determine, predict, and increase solubility. Solvent additives (osmolytes) that stabilize proteins are listed with a description of their effects on proteins and on the solvation properties of water. Special attention is given to areas where solubility limitations pose major problems, as in the preparation of highly concentrated solutions of recombinant proteins for structural determination with NMR and X-ray crystallography, refolding of inclusion body proteins, studies of membrane protein dynamics, and in the formulation of proteins for pharmaceutical use. Structural factors relating to solubility and possibilities for protein engineering are analyzed.

5'-3' and 3'-5' translation of the same RNA results in hydropathically similar peptides that are antigenically related

When a single RNA sequence is read in either the 5'-3' or 3'-5 direction, the translated peptides often are hydropathically similar even though their sequences may be different. To investigate whether hydropathically similar peptides might also be antigenically related, two peptides were synthesized from the substance P anti-sense RNA transcript: CAU CAA UCC AAA GAA CUG CUG AGG CUU GGG UCG. Translation of this RNA in the 5'-3' direction and in the 3'-5' direction resulted in two different peptides. HQSKELLRLGS and AGFGVVKKPNY, respectively. As anticipated, both peptides shared similar hydropathic profiles but were quite different with respect to their sequences. To examine their antigenic relatedness, mice were immunized with either peptide, and monoclonal antibodies were produced. Using an enzyme-linked immunosorbent assay, it was possible to demonstrate that the majority of monoclonal antibodies, selected for reactivity against the original immunogen, also reacted with the other peptide. The observed binding was determined to be specific since reactivity could be blocked with either soluble peptide. Thus, we demonstrate that hydropathically similar peptides obtained from the same RNA but translated in opposite directions are antigenically related despite difference in amino acid sequences.

Haemopoietic receptors and helical cytokines

A bewitching interplay of proteins, variously clothed as chemical messengers and cellular receptors, control the pace of growth and the course of progressive differentiation in blood cell types. The messengers are lymphokines, interleukins, colony-stimulating factors, growth hormones and interferons: generically, the cytokines. The second components of the regulatory pairs are membrane-spanning receptor proteins: these molecules transduce the specific binding of cognate cytokines into a mitogenic cellular response. In this article, Fernando Bazan discusses a provocative structural model for cytokine-receptor interactions which, if correct, will alter perceptions of the evolutionary design of the haemopoietic system.

Design of DNA-binding peptides based on the leucine zipper motif

A class of transcriptional regulator proteins bind to DNA at dyad-symmetric sites through a motif consisting of (i) a "leucine zipper" sequence that associates into noncovalent, parallel, alpha-helical dimers and (ii) a covalently connected basic region necessary for binding DNA. The basic regions are predicted to be disordered in the absence of DNA and to form alpha helices when bound to DNA. These helices bind in the major groove forming multiple hydrogen-bonded and van der Waals contacts with the nucleotide bases. To test this model, two peptides were designed that were identical to natural leucine zipper proteins only at positions hypothesized to be critical for dimerization and DNA recognition. The peptides form dimers that bind specifically to DNA with their basic regions in alpha-helical conformations.

Crystal structure of alpha 1: implications for protein design

X-ray diffraction shows the structure of a synthetic protein model, formed from noncovalent self-association of a 12-residue peptide and of sulfate ions at low pH. This peptide is a fragment of a 16-residue polypeptide that was designed to form an amphiphilic alpha helix with a ridge of Leu residues along one helical face. By interdigitation of the leucines of four such helices, the design called for self-association into a four-alpha-helical bundle. The crystal structure (2.7 angstrom resolution; R factor = 0.215) reveals a structure more complex than the design, with both a tetramer and a hexamer. The alpha-helical tetramer with leucine interior has more oblique crossing angles than most four-alpha-helical bundles; the hexamer has a globular hydrophobic core of 12 leucine residues and three associated sulfate ions. Computational analysis suggests that the hexameric association is tighter than the tetrameric one. The consistency of the structure with the design is discussed, as well as the divergence.

Calcium-induced peptide association to form an intact protein domain: 1H NMR structural evidence

The 70-residue carboxyl-terminal domain of the muscle contractile protein troponin-C contains two helix-loop-helix calcium (Ca)-binding sites that are related to each other by approximate twofold rotational symmetry. Hydrophobic residues from the helices and a short three residue beta sheet at the interface of the two sites act to stabilize the protein domain in the presence of Ca. A synthetic 34-residue peptide representing one of these sites (site III) has been synthesized and studied by H-1 nuclear magnetic resonance (NMR) spectroscopy. In solution this peptide undergoes a Ca-induced conformational change to form the helix-loop-helix Ca-binding motif. Two-dimensional nuclear Overhauser effect spectra have provided evidence for the formation of a beta sheet and interactions between several hydrophobic residues from opposing helices as found in troponin-C. It is proposed that a symmetric two-site dimer similar in tertiary structure to the carboxyl-terminal domain of troponin-C forms from the assembly of two site III peptides in the Ca-bound form.

Stabilization of helical structure in two 17-residue amphipathic analogues of the C-terminal peptide of cytochrome C

The conformations of two 17-residue peptide analogues derived from the C-terminal sequence of pigeon cytochrome c (native sequence = KAERADLIAYLKQATAK) were examined in aqueous and lipid environments by CD spectroscopy. The two analogues, KKLLKKLIAYLKQATAK (K peptide) and EELLEELIAYLKQATAK (E peptide), were made amphipathic with respect to helical segregation by substituting a 6-residue sequence at the N-terminus of the native peptide. Their structures were compared to the native peptide under aqueous conditions of varying pH and temperature, and in the presence of liposomes composed of phosphatidylcholine and phosphatidylserine in the ratio of 9:1. The results indicated that the native peptide remains unstructured under all the conditions examined even though this region of the native molecule is surface exposed and helical. The E peptide, however, was helical under aqueous conditions at 25 degrees C from pH 2-10 with a maximum helicity at pH 4 (54% helix from analysis of CD data). The ellipticity of the E peptide at pH 4 and 8 was concentration dependent, indicating an aggregation phenomenon. In studies in which the CD spectrum was measured at different temperatures, the E peptide became more helical at lower temperatures at pH 4 but not at pH 8. Upon interaction with a lipid membrane in the form of liposomes, there appeared to be a slight destabilization in the structure of the E peptide. The K peptide in an aqueous environment behaved like the native peptide in that it was structureless at all pHs and temperatures examined. In the presence of liposomes, however, this peptide had a high helical content (75% helix from analysis of CD data). These findings suggest that while stabilization of the helix dipole with negative charges at the N-terminus are important in inducing helical conformation in the E peptide, hydrophobic interactions created during aggregation appear to provide the principal stabilizing force. The results with the K peptide demonstrate that the positive N-terminal sequence of this peptide is able to interact with the negatively charged head groups in the phospholipid membrane in such a fashion as to stabilize a helical structure that is not apparent in an aqueous environment alone.

Ion-channel properties of mastoparan, a 14-residue peptide from wasp venom, and of MP3, a 12-residue analogue

Mastoparan, a 14-residue peptide, has been investigated with respect to its ability to form ion channels in planar lipid bilayers. In the presence of 0.3-3.0 microM mastoparan, two types of activity are seen. Type I activity is characterized by discrete channel openings, exhibiting multiple conductance levels in the range 15-700 pS. Type II activity is characterized by transient increases in bilayer conductance, up to a maximum of about 650 pS. Both type I and type II activities are voltage dependent. Channel activation occurs if the compartment containing mastoparan is held at a positive potential; channel inactivation if the same compartment is held at a negative potential. Channel formation is dependent on ionic strength; channel openings are only observed at KCl concentrations of 0.3 M or above. Furthermore, raising the concentration of KCl to 3.0 M stabilizes the open form of the channel. Mastoparan channels are weakly cation selective, PK/Cl approximately 2. A 12-residue analogue, des-Ile1,Asn2-mastoparan, preferentially forms type I channels. The ion channels formed by these short peptides may be modelled in terms of bundles of transmembrane alpha-helices.

Design of crystalline helices of short oligopeptides as a possible model for nucleation of alpha-helix: role of water molecules in stabilizing helices

We have designed, synthesized, crystallized, and performed x-ray analysis of several hydrophobic tripeptides that show an extended near alpha-helical structure in the crystalline state. All of the tripeptides that show this remarkably stable helix crystallize with two or three water molecules; they all have glycine at the N terminus and have increasing hydrophobicity as one moves from the N to C terminus. Even though three residues in the oligomer are not sufficient to complete a turn, one of the water molecules acts as an added residue and links up adjacent tripeptide segments along the helix axis so that in the crystal, the helix appears effectively as one long continuous helix. Two of these tripeptides are stabilized by two water molecules that enable the peptides to complete a turn of the helix and extend the helical structure throughout the crystal by linking translationally related peptides by hydrogen bonds. In two other peptides, these roles are played by three rather than two water molecules. Though these tripeptides have different crystal symmetry, they all show the basic pattern of hydrated helix and packing, indicating the strong conformational preference for a stable structure even for these tripeptides. Such conformationally stable hydrated structures for short specific related sequences illustrate their possible importance in nucleating protein folding and in the role water molecules play in such events.

Synthesis and biophysical characterization of engineered topographic immunogenic determinants with alpha alpha topology

Model peptides with predetermined secondary, tertiary, and quaternary conformation have been successfully designed, synthesized, and characterized in an attempt to mimic the three-dimensional structure of an antigenic determinant. This work is a continuing effort to map the antigenic structure of the protein antigen lactate dehydrogenase C4 (LDH-C4) to develop a contraceptive vaccine. A putative topographic determinant with alpha alpha topology which associates into four-helix bundles was designed on the basis of the framework model of protein folding. An idealized amphiphilic 18-residue sequence (alpha 1) and a 40-residue alpha alpha fold (alpha 3) have been shown to form stable 4-helix structures in solution with a free energy of association on the order of -20.8 kcal/mol (tetramerization of alpha 1) and -7.8 kcal/mol (dimerization of alpha 3). Both alpha 1 and alpha 3 form stable monolayers at the air-water interface. The CD spectra of Langmuir-Blodgett monolayers are characteristically alpha-helical. Both CD and FTIR spectroscopic studies reval a high degree of secondary structure. The SAXS data strongly suggest that the helices are arranged in a four-helix bundle since the radius of gyration of 17.2 A and the vector distribution function are indicative of a prolate ellipsoid of axial dimensions and molecular weight appropriate for the four-helix bundle. The major contribution to the formation and stabilization of alpha 1 and alpha 3 is believed to be hydrophobic interaction between the amphiphilic alpha-helices. The displayed heptad repeat, helix dipole, ion pairs, and the loop sequence may have also contributed to the overall stability and antiparallel packing of the helices. A detailed structural analysis of a relevant topographic immunogenic determinant will elucidate the nature of antigen-antibody interactions as well as provide insight into protein folding intermediates.

Synthesis and structural stability of helichrome as an artificial hemeproteins

A detailed procedure is described for the synthesis of helichrome, which is the first successful example of polypeptide-based artificial hemeprotein. The segment synthesis-condensation approach used for the assembly of small proteins has proven to be extremely useful for protein mimetics as well. The final deprotection was performed using the TMSOTf-thioanisole method instead of the less-convenient hydrogen fluoride method. The unfolding transition of the alpha-helical conformation of helichrome induced by guanidine hydrochloride was studied to understand the stability and dynamics of the folded structure. The resulting parameters (C0.5 = 5.2 M and delta GH2O = -4.4 kcal mol-1) characterizing helichrome denaturation were comparable to that of native globular proteins.

Conformationally constrained alpha-helical peptide models for protein ion channels

Recently we described the design, synthesis, and characterization of some simple amphiphilic alpha-helical models for protein ion channels. These peptides, composed of only Leu and Ser residues, are hypothesized to form helical bundles capable of passing ions across phospholipid bilayers. In an effort to demonstrate that the peptides are, in fact, helical in their active ion-conducting state, the conformationally constrained amino acid, C alpha, C alpha-dimethylglycine (alpha-aminoisobutyric acid, Aib), was introduced simultaneously at three positions into one of the model peptides, H2N-(Leu-Ser-Leu-Leu-Leu-Ser-Leu)3-CONH2, giving H2N-(Leu-Ser-Leu-Aib-Leu-Ser-Leu)3-CONH2. Examination of a tetrameric model for the channel suggested that this substitution should have a minimal effect on conductance. CD spectroscopy of the Aib-modified and original peptide in phospholipid vesicles indicated that both were highly alpha-helical. Furthermore, the Aib-containing peptide formed proton channels nearly identical in conductance to the original peptide.

The design of idealized alpha/beta-barrels: analysis of beta-sheet closure requirements

The 8-fold parallel alpha/beta-barrel topology is encountered in proteins that display an impressive variety of functions, suggesting that this topology may be a rather nonspecific and stable folding motif. Consequently, this motif can be considered as an interesting framework to design novel proteins. It has been shown that the shape of the beta-sheet portion of the barrel can be approximated by a hyperboloid. This geometric object may therefore be used as a scaffold to construct an idealized eight-stranded beta-barrel. To facilitate the de novo design of such structures, a collection of modeling tools has been developed allowing secondary structure elements to be mapped onto the scaffold surface and rotation and translation operations to be performed about user defined axes while evaluating their contribution to the conformational energy of the system. These tools have been applied in a systematic study assessing the phi, psi requirements to design symmetric eight stranded beta barrels with optimal hydrogen bonding between adjacent beta-strands. It is observed that: (a) the beta-sheet structure can be closed without introducing irregular stagger between beta-strands and (b) the region of phi, psi dihedral angle space compatible with the formation of regular symmetric eight stranded beta-barrels coincides with the phi, psi region corresponding to average beta-strands in known protein structures, suggesting that barrel closure does not impose gross constraints on beta-strand geometry.

A molecular editor program

A program called MOLEDITOR has been developed for manipulating images of molecules and various geometric objects in real time and in three dimensions using Evans and Sutherland PS300 graphics displays and Digital VAX computers. In addition to the normal viewing operations of rotation, translation, clipping and scaling, the program allows molecules and objects to be placed in local frames for motion relative to other molecules or objects and allows several modes of motion about individual bonds. MOLEDITOR also allows the display and manipulation of molecules and objects within crystal lattices and of the lattices themselves via manipulation of the unit cell parameters. The program also supports display and manipulation of animated frames from dynamics simulations. It is also possible to dynamically monitor coordinates of and distances between items. An easy-to-use scripting capability allows the user to generate sequences of motion that can be recorded directly from the display.

Interleukin 2: prototype for a new generation of immunoactive pharmaceuticals

Molecular biological techniques have revealed the interleukin-2 receptor to be a dimer composed of one alpha-subunit and one beta-subunit which interact noncovalently in a cooperative manner. Site-directed mutagenesis, in conjunction with structural analysis, is beginning to clarify the relationship between structural components of the receptor and their function, and Thomas Ciardelli and Kendall Smith explain why this is bringing drug developers closer to the design of IL-2 agonists and antagonists.

Complementary peptides as interactive sites for protein binding

The observation that pairs of peptides encoded from complementary strands of nucleic acids can bind one another suggested to us that proteins might use similar "complementary peptide" sequences to facilitate their interactions with other proteins or peptides. Here we review evidence to support this hypothesis. Investigations by us and by others showed that receptors and their ligands may use complementary peptide sequences or their analogs to facilitate binding. In addition, antibodies made against pairs of complementary peptides bind one another in an idiotypic/anti-idiotypic fashion. Taken together, these findings suggest that complementary peptide sequences may be used as recognition units by pairs of proteins which interact.