Amyloid-forming peptides selected proteolytically from phage display library

We demonstrated that amyloid-forming peptides could be selected from phage-displayed library via proteolysis-based selection protocol. The library of 28-residue peptides based on a sequence of the second zinc finger domain of Zif268, and computationally designed betabetaalpha peptide, FSD-1, was presented monovalently on the surface of M13 phage. The library coupled the infectivity of phage particles to proteolytic stability of a peptide introduced into the coat protein III linker. It was designed to include variants with a strong potential to fold into betabetaalpha motif of zinc finger domains, as expected from secondary structure propensities, but with no structure stabilization via zinc ion coordination. As our primary goal was to find novel monomeric betabetaalpha peptides, the library was selected for stable domains with the assumption that folded proteins are resistant to proteolysis. After less than four rounds of proteolytic selection with trypsin, chymotrypsin, or proteinase K, we obtained a number of proteolysis-resistant phage clones containing several potential sites for proteolytic attack with the proteinases. Eight peptides showing the highest proteolysis resistance were expressed and purified in a phage-free form. When characterized, the peptides possessed proteolytic resistance largely exceeding that of the second zinc finger domain of Zif268 and FSD-1. Six of the characterized peptides formed fibrils when solubilized at high concentrations. Three of them assembled into amyloids as determined through CD measurements, Congo red and thioflavin T binding, and transmission electron microscopy.     

Chimeric Protein Engineering

Protein stability can be enhanced by the incorporation of non-natural amino acids and semi-rigid peptidomimetics to lower the entropic penalty upon protein folding through preorganization. An example is the incorporation of aminoisobutyric acid (Aib, α-methylalanine) into proteins to restrict the Φ and Ψ backbone angles adjacent to Aib to those associated with helix formation. Reverse-turn analogs were introduced into the sequences of HIV protease and ribonuclease A that enhanced their stability and retained their native enzymatic activity. In this work, a chimeric protein, design_4, was engineered, in silico, by replacing the C-terminal helix of full sequence design protein (FSD-1) with a semi-rigid helix mimetic. Residues 1–16 of FSD-1 was ligated in silico with the N-terminus of a phenylbipyridyl-based helix mimetic to form design_4. The designed chimeric protein was stable and maintained the designed fold in a 100-nanosecond molecular dynamics simulation at 280 K. Its β-hairpin adopted conformations that formed three additional hydrogen bonds. Compared to FSD-1, design_4 contained fewer peptide bonds and internal degrees of freedom; it should, therefore, be more resistant to proteolytic degradation and denaturation.     

Modulating calmodulin binding specificity through computational protein design

We report the computational redesign of the protein-binding interface of calmodulin (CaM), a small, ubiquitous Ca(2+)-binding protein that is known to bind to and regulate a variety of functionally and structurally diverse proteins. The CaM binding interface was optimized to improve binding specificity towards one of its natural targets, smooth muscle myosin light chain kinase (smMLCK). The optimization was performed using optimization of rotamers by iterative techniques (ORBIT), a protein design program that utilizes a physically based force-field and the Dead-End Elimination theorem to compute sequences that are optimal for a given protein scaffold. Starting from the structure of the CaM-smMLCK complex, the program considered 10(22) amino acid residue sequences to obtain the lowest-energy CaM sequence. The resulting eightfold mutant, CaM_8, was constructed and tested for binding to a set of seven CaM target peptides. CaM_8 displayed high binding affinity to the smMLCK peptide (1.3nM), similar to that of the wild-type protein (1.8nM). The affinity of CaM_8 to six other target peptides was reduced, as intended, by 1.5-fold to 86-fold. Hence, CaM_8 exhibited increased binding specificity, preferring the smMLCK peptide to the other targets. Studies of this type may increase our understanding of the origins of binding specificity in protein-ligand complexes and may provide valuable information that can be used in the design of novel protein receptors and/or ligands.     

Understanding the folding and stability of a zinc finger-based full sequence design protein with replica exchange molecular dynamics simulations

Full sequence design protein FSD-1 is a designed protein based on the motif of zinc finger protein. In this work, its folding mechanism and thermal stability are investigated using the replica exchange molecular dynamics model with the water molecules being treated explicitly. The results show that the folding of the FSD-1 is initiated by the hydrophobic collapse, which is accompanied with the formation of the C-terminal alpha-helix. Then the folding proceeds with the formation of the beta-hairpin and the further package of the hydrophobic core. Compared with the beta-hairpin, the alpha-helix has much higher stability. It is also found that the N-capping motif adopted by the FSD-1 contributes to the stability of the alpha-helix dramatically. The hydrophobic contacts made by the side chain of Tyr3 in the native state are essential for the stabilization of the beta-hairpin. It is also found that the folding of the N-terminal beta-hairpin and the C-terminal alpha-helix exhibits weak cooperativity, which is consistent with the experimental data. Meanwhile, the folding pathway is compared between the FSD-1 and the target zinc finger peptide, and the possible role of the zinc ion on the folding pathway of zinc finger is proposed. Proteins 2007. (c) 2007 Wiley-Liss, Inc.     

Computational design of calmodulin mutants with up to 900-fold increase in binding specificity

Calmodulin (CaM) is a ubiquitous second messenger protein that regulates a variety of structurally and functionally diverse targets in response to changes in Ca²⁺ concentration. CaM-dependent protein kinase II (CaMKII) and calcineurin (CaN) are the prominent CaM targets that play an opposing role in many cellular functions including synaptic regulation. Since CaMKII and CaN compete for the available Ca²⁺/CaM, the differential affinity of these enzymes for CaM is crucial for achieving a balance in Ca²⁺ signaling. We used the computational protein design approach to modify CaM binding specificity for these two targets. Starting from the X-ray structure of CaM in complex with the CaM-binding domain of CaMKII, we optimized CaM interactions with CaMKII by introducing mutations into the CaM sequence. CaM optimization was performed with a protein design program, ORBIT, using a modified energy function that emphasized intermolecular interactions in the sequence selection procedure. Several CaM variants were experimentally constructed and tested for binding to the CaMKII and CaN peptides using the surface plasmon resonance technique. Most of our CaM mutants demonstrated small increase in affinity for the CaMKII peptide and substantial decrease in affinity for the CaN peptide compared to that of wild-type CaM. Our best CaM design exhibited an about 900-fold increase in binding specificity towards the CaMKII peptide, becoming the highest specificity switch achieved in any protein-protein interface through the computational protein design approach. Our results show that computational redesign of protein-protein interfaces becomes a reliable method for altering protein binding affinity and specificity.     

A Second Look at Mini-Protein Stability: Analysis of FSD-1 Using Circular Dichroism, Differential Scanning Calorimetry, and Simulations

Mini-proteins that contain <50 amino acids often serve as model systems for studying protein folding because their small size makes long timescale simulations possible. However, not all mini-proteins are created equal. The stability and structure of FSD-1, a 28-residue mini-protein that adopted the ββα zinc-finger motif independent of zinc binding, was investigated using circular dichroism, differential scanning calorimetry, and replica-exchange molecular dynamics. The broad melting transition of FSD-1, similar to that of a helix-to-coil transition, was observed by using circular dichroism, differential scanning calorimetry, and replica-exchange molecular dynamics. The N-terminal β-hairpin was found to be flexible. The FSD-1 apparent melting temperature of 41°C may be a reflection of the melting of its α-helical segment instead of the entire protein. Thus, despite its attractiveness due to small size and purposefully designed helix, sheet, and turn structures, the status of FSD-1 as a model system for studying protein folding should be reconsidered.     

NMR study and molecular dynamics simulations of optimized beta-hairpin fragments of protein G

The stability and structure of several beta-hairpin peptide variants derived from the C-terminus of the B1 domain of protein G were investigated by a number of experimental and computational techniques. Our analysis shows that the structure and stability of this hairpin can be greatly affected by one or a few simple mutations. For example, removing an unfavorable charge near the N-terminus of the peptide (Glu42 to Gln or Thr) or optimization of the N-terminal charge-charge interactions (Gly41 to Lys) both stabilize the peptide, even in water. Furthermore, a simple replacement of a charged residue in the turn (Asp47 to Ala) changes the beta-turn conformation. Finally, we show that the effects of combining these single mutations are additive, suggesting that independent stabilizing interactions can be isolated and evaluated in a simple model system. Our results indicate that the structure and stability of this beta-hairpin peptide can be modulated in numerous ways and thus contributes toward a more complete understanding of this important model beta-hairpin as well as to the folding and stability of larger peptides and proteins.     

In silico identification and design of potent peptide inhibitors against PDZ-3 domain of Postsynaptic Density Protein (PSD-95)

Unique intrinsic properties of peptides like low toxicity, high biological activity, and specificity make them attractive therapeutic agents. PDZ-binding peptide inhibitors have been demonstrated for curing of Alzheimer, Parkinson, Dementia, and other central nervous system ailments. In this article, we report the successful use of an integrated computational protocol to analyze the structural basis of how peptides bind to the shallow groove of the third PDZ domain (PDZ-3) from the postsynaptic density (PSD-95) protein. This protocol employs careful and precise computational techniques for design of new strategy for predicting novel and potent peptides against PDZ protein. We attempted to generate a pharmacophore model using crystal structure of peptide inhibitor bound to the PDZ-3. A highly specific and sensitive generated pharmacophore model was used for screening virtual database generated using different combination of amino acid substitutions as well as decoy peptide database for its sensitivity and specificity. Identified hit peptides were further analyzed by docking studies, and their stability analyzed using solvated molecular dynamics. Quantum Mechanics/Molecular Mechanics (QM/MM) interaction energy and GMX-PBSA scoring schemes were used for ranking of stable peptides. Computational approach applied here generated encouraging results for identifying peptides against PDZ interaction model. The workflow can be further exercised as a virtual screening technique for reducing the search space for candidate target peptides against PDZ domains.     

Cation-pi interactions studied in a model coiled-coil peptide

Cation-pi interactions between aromatic amino acids and the positively charged residues lysine and arginine have been proposed to play an important role in stabilizing protein structure. We have used a peptide that adopts a coiled coil structure as a model system to evaluate the energetic contribution of cation-pi interactions to protein folding. Peptides were designed in which phenylalanine, tyrosine, and tryptophan were placed at a solvent-exposed position of the helix, one turn removed from an arginine residue that could provide a favorable cation-pi interaction. Only the arginine-phenylalanine pairing provided significant stabilization of the peptide structure and it appears that hydrophobic packing, rather than the cation-pi effect, is more likely to be responsible for the stability of this peptide. We conclude that any stabilizing effect of cation-pi interactions in these peptides is much smaller than that predicted from computational studies.     

Molecular design of Stat3‐derived peptide selectivity between BET proteins Brd2 and Brd4 in ovarian cancer

Stat3 signaling has been recognized as a potential therapeutic target of human ovarian cancer. The signaling is transducted through the peptide‐medicated interaction of Stat3 with BET family members Brd2 and Brd4 –– 2 highly homologous proteins involved in differential downstream pathways. Here, we reported a successful design of peptide selectivity between the Brd2 and Brd4. The design resulted in 3 linear peptides SMSLQCX YLGVA, QSKVLTX SYWGA, and RQCNLGX LYMNY with high or moderate selectivity for Brd2 over Brd4 (S = 3.3‐fold, 6.8‐fold, and 4.2‐fold, respectively) as compared with the native Stat3 peptide ²⁸¹HNLLRIX QFLQS²⁹² (S = 2.5‐fold). Structural analysis revealed that peptide N‐terminus and hydrogen bonds play important roles in the peptide interaction stability and specificity with Brd2 and Brd4. This study would help to establish an integrated in silico‐in vitro method for rational molecular design of peptide ligand selectivity between homologous protein receptors.     

Designed protein G core variants fold to native-like structures: sequence selection by ORBIT tolerates variation in backbone specification

The solution structures of two computationally designed core variants of the beta 1 domain of streptococcal protein G (G beta 1) were solved by (1)H NMR methods to assess the robustness of amino acid sequence selection by the ORBIT protein design package under changes in protein backbone specification. One variant has mutations at three of 10 core positions and corresponds to minimal perturbations of the native G beta 1 backbone. The other, with mutations at six of 10 positions, was calculated for a backbone in which the separation between G beta 1's alpha-helix and beta-sheet was increased by 15% relative to native G beta 1. Exchange broadening of some resonances and the complete absence of others in spectra of the sixfold mutant bespeak conformational heterogeneity in this protein. The NMR data were sufficiently abundant, however, to generate structures of similar, moderately high quality for both variants. Both proteins adopt backbone structures similar to their target folds. Moreover, the sequence selection algorithm successfully predicted all core chi(1) angles in both variants, five of six chi(2) angles in the threefold mutant and four of seven chi(2) angles in the sixfold mutant. We conclude that ORBIT calculates sequences that fold specifically to a geometry close to the template, even when the template is moderately perturbed relative to a naturally occurring structure. There are apparently limits to the size of acceptable perturbations: In this study, the larger perturbation led to undesired dynamic behavior.     

Designing a 20-residue protein

Truncation and mutation of a poorly folded 39-residue peptide has produced 20-residue constructs that are >95% folded in water at physiological pH. These constructs optimize a novel fold, designated as the 'Trp-cage' motif, and are significantly more stable than any other miniprotein reported to date. Folding is cooperative and hydrophobically driven by the encapsulation of a Trp side chain in a sheath of Pro rings. As the smallest protein-like construct, Trp-cage miniproteins should provide a testing ground for both experimental studies and computational simulations of protein folding and unfolding pathways. Pro Trp interactions may be a particularly effective strategy for the a priori design of self-folding peptides.     

Adaptive recognition in RNA complexes with peptides and protein modules

Recently, progress has been made towards the structural characterization of the novel folds of RNA-bound arginine-rich peptides and the architecture of their peptide-binding RNA pockets in viral and phage systems. These studies are based on an approach whereby the peptide and RNA components are minimalist modular domains that undergo adaptive structural transitions upon complex formation. Such complexes are characterized by recognition alignments in which the tertiary fold of the RNA generates binding pockets with the potential to envelop minimal elements of protein secondary structure. Strikingly, the peptides fold as isolated alpha-helical or beta-hairpin folds within their RNA major-groove targets, without the necessity of additional appendages for anchorage within the binding pocket. The RNA peptide-binding pocket architectures are sculptured through precisely positioned mismatches, triples and looped-out bases, which accommodate amino acid sidechains through hydrophobic, hydrogen bonding and ionic intermolecular contacts. By contrast, protein modules associated with the HIV-1 nucleocapsid and MS2 phage coat target their RNA binding sites through the insertion of specificity-determining RNA base residues within conserved hydrophobic pockets and crevices on the protein surface, with the bases anchored through hydrogen bonding interactions. These alternative strategies of RNA recognition at the peptide and protein module level provide novel insights into the principles, patterns and diversity of the adaptive transitions associated with the recognition process.     

Importance of alpha-helix N-capping motif in stabilization of betabetaalpha fold

FSD-1 (full sequence design 1) is a protein folded in a betabetaalpha motif, designed on the basis of the second zinc finger domain of Zif268 by a substitution of its metal coordination site with a hydrophobic core. In this work, we analyzed the possibility of introducing the DNA recognition motif of the template zinc finger (S(13)RSDH(17)) into FSD-1 sequence in order to obtain a small DNA-binding module devoid of cross-link(s) or metal cofactors. The hybrid protein was unfolded, as judged by CD and NMR criteria. To reveal the role of each of the five amino acids, which form the N-capping motif of the alpha-helix, we analyzed conformational and stability properties of eight FSD-1 mutants. We used a shielded methyl group of Leu 18 and a CD signal at 215 nm as a convenient measure of the folded state. Glu 17-->His substitution at the N(3) in S(13)NEKE(17) variant decreased the folded structure content from 90% to 25% (equivalent to 1.8 kcal * mole(-1) destabilization) by disruption of N-capping interactions, and had the most significant effect among single mutants studied here. The N(cap) Asn 14 substitution with Arg considerably decreased stability, reducing structure content from 90% to 40% (1.4 kcal * mole(-1) destabilization) by disruption of a helix-capping hydrogen bond and destabilization of a helix macrodipole. The N(1) Glu 15-->Ser mutation also produced a considerable effect (1.0 kcal * mole(-1) destabilization), again emphasizing the significance of electrostatic interactions in alpha-helix stabilization.     

Rational design and improvement of the dimerization‐disrupting peptide selectivity between ROCK‐I and ROCK‐II kinase isoforms in cerebrovascular diseases

Human rho‐associated coiled‐coil forming kinases (ROCKs) ROCK‐I and ROCK‐II have been documented as attractive therapeutic targets for cerebrovascular diseases. Although ROCK‐I and ROCK‐II share a high degree of structural conservation and are both present in classic rho/ROCK signaling pathway, their downstream substrates and pathological functions may be quite different. Selective targeting of the two kinase isoforms with traditional small‐molecule inhibitors is a great challenge due to their surprisingly high homology in kinase domain (~90%) and the full identity in kinase active site (100%). Here, instead of developing small‐molecule drugs to selectively target the adenosine triphosphate (ATP) site of two isoforms, we attempt to design peptide agents to selectively disrupt the homo‐dimerization event of ROCK kinases through their dimerization domains which have a relatively low conservation (~60%). Three helical peptides H1, H2, and H3 are split from the kinase dimerization domain, from which the isolated H2 peptide is found to have the best capability to rebind at the dimerization interface. A simulated annealing (SA) iteration method is used to improve the H2 peptide selectivity between ROCK‐I and ROCK‐II. The method accepts moderate degradation in peptide affinity in order to maximize the affinity difference between peptide binding to the two isoforms. Consequently, hundreds of parallel SA runs yielded six promising peptide candidates with ROCK‐I over ROCK‐II (I over II [IoII]) calculated selectivity and four promising peptide candidates with ROCK‐II over ROCK‐I (II over I [IIoI]) calculated selectivity. Subsequent anisotropy assays confirm that the selectivity values range between 13.2‐fold and 83.9‐fold for IoII peptides, and between 5.8‐fold and 21.2‐fold for IIoI peptides, which are considerably increased relative to wild‐type H2 peptide (2.6‐fold for IoII and 2.0‐fold for IIoI). The molecular origin of the designed peptide selectivity is also analyzed at structural level; it is revealed that the peptide residues can be classified into conserved, non‐conserved, and others, in which the non‐conserved residues play a crucial role in defining peptide selectivity, while conserved residues confer stability to kinase‐peptide binding.     

The side chain interaction index as a tool for predicting fast-folding elements and the structure and stability of engineered peptides

The side chain interaction index (SCII) is a method of calculating the propensity for short-range interactions among side chains within a peptide sequence. Here, it is shown that the SCII values of secondary structure elements that have been shown to fold early and independently cluster separately from those of structures that fold later and/or are dependent on long-range interactions. In addition, the SCII values of engineered peptides that spontaneously adopt a particular desired fold in solution are significantly different from those of engineered peptides that fail to exhibit a stable conformation. Thus, the SCII, as a measure of local structural stability, constitutes a useful tool in folding prediction and in protein/peptide engineering. A program that allows rapid calculation of SCII values is presented.     

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.     

Reduced tendency to form a beta turn in peptides from the P22 tailspike protein correlates with a temperature-sensitive folding defect

A family of mutants of the P22 bacteriophage tailspike protein has been characterized as temperature sensitive for folding (tsf) by King and co-workers [King, J. (1986) Bio/Technology 4, 297-303]. There is substantial evidence that the tsf mutations alter the folding pathway but not the stability of the final folded protein. Several point mutations are known to cause the tsf phenotype; most of these occur in regions of the tailspike sequence likely to take up reverse turns. Hence, it has been hypothesized that the correct folding of the P22 tailspike protein requires formation of turns and that the mutations causing tsf phenotypes interfere at this critical stage. We have tested this hypothesis by study of isolated peptides corresponding to a region of the P22 tailspike harboring a tsf mutation. Comparison of the tendencies of wild-type and tsf sequences to adopt turn conformations was achieved by the synthesis of peptides with flanking cysteine residues and the use of a thiol-disulfide exchange assay. We find that the wild-type sequence, either as a decapeptide (Ac-CVKFPGIETC-CONH2) or as a dodecapeptide (Ac-CYVKFPGIETLC-CONH2), has a 3-5-fold greater tendency for its termini to approach closely enough to form the intramolecular disulfide than do the peptide sequences corresponding to the tsf mutant sequences, which have a Gly----Arg substitution (Ac-CVKFPRIETC-CONH2 or Ac-CYVKFPRIETLC-CONH2). A peptide with a D-Arg substituted for the Gly has a slightly higher turn propensity than does the wild type. Together with data from nuclear magnetic resonance analysis of the oxidized peptides, this suggests that a type II beta turn is favored by the wild-type sequence. Our results on isolated peptides from the P22 tailspike protein support the model for its folding that includes reverse turn formation as a critical step.     

Improving coiled-coil stability by optimizing ionic interactions

Alpha-helical coiled coils are a common protein oligomerization motif stabilized mainly by hydrophobic interactions occurring along the coiled-coil interface. We have recently designed and solved the structure of a two-heptad repeat coiled-coil peptide that is stabilized further by a complex network of inter- and intrahelical salt-bridges in addition to the hydrophobic interactions. Here, we extend and improve the de novo design of this two heptad-repeat peptide by four newly designed peptides characterized by different types of ionic interactions. The contribution of these different types of ionic interactions to coiled-coil stability are analyzed by CD spectroscopy and analytical ultracentrifugation. We show that all peptides are highly alpha-helical and two of them are 100% dimeric under physiological conditions. Furthermore, we have solved the X-ray structure of the most stable of these peptides and the rational design principles are verified by comparing this structure to the structure of the parent peptide. We show that by combining the most favorable inter- and intrahelical salt-bridge arrangements it is possible to design coiled-coil oligomerization domains with improved stability properties.