The role of a beta-bulge in the folding of the beta-hairpin structure in ubiquitin

Generative β-hairpin design using a residue-based physicochemical property landscape

De novo peptide design is a new frontier that has broad application potential in the biological and biomedical fields. Most existing models for de novo peptide design are largely based on sequence homology that can be restricted based on evolutionarily derived protein sequences and lack the physicochemical context essential in protein folding. Generative machine learning for de novo peptide design is a promising way to synthesize theoretical data that are based on, but unique from, the observable universe. In this study, we created and tested a custom peptide generative adversarial network intended to design peptide sequences that can fold into the β-hairpin secondary structure. This deep neural network model is designed to establish a preliminary foundation of the generative approach based on physicochemical and conformational properties of 20 canonical amino acids, for example, hydrophobicity and residue volume, using extant structure-specific sequence data from the PDB. The beta generative adversarial network model robustly distinguishes secondary structures of β hairpin from α helix and intrinsically disordered peptides with an accuracy of up to 96% and generates artificial β-hairpin peptide sequences with minimum sequence identities around 31% and 50% when compared against the current NCBI PDB and nonredundant databases, respectively. These results highlight the potential of generative models specifically anchored by physicochemical and conformational property features of amino acids to expand the sequence-to-structure landscape of proteins beyond evolutionary limits.

PyProtif: a PyMol plugin to retrieve and visualize protein motifs for structural studies

Proteins often possess several motifs and the ones with similar motifs were found to have similar biochemical properties and thus related biological functions. Thereby, multiple databases were developed to store information on such motifs in proteins. For instance, PDBsum stores the results of Promotif's generated structural motifs and Pfam stores pre-computed patterns of functional domains. In addition to the fact that all this stored information is extremely useful, we can further augment its importance if we ought to integrate these motifs into visualization software. In this work, we have developed PyProtif, a plugin for the PyMOL molecular visualization program, which automatically retrieves protein structural and functional motifs from different databases and integrates them in PyMOL for visualization and analyses. Through an expendable menu and a user-friendly interface, the plugin grants the users the ability to study simultaneously multiple proteins and to select and manipulate each motif separately. Thus, this plugin will be of great interest for structural, evolutionary and classification studies of proteins.

Potential of Raman spectroscopic techniques to study proteins

Proteins are large, complex molecules responsible for various biological processes. However, protein misfolding may lead to various life-threatening diseases. Therefore, it is vital to understand the shape and structure of proteins. Despite numerous techniques, a mechanistic understanding of the protein folding process is still unclear. Therefore, new techniques are continually being explored. In the present article, we have discussed the importance of Raman spectroscopy, Raman Optical Activity (ROA) and various other advancements in Raman spectroscopy to understand protein structure and conformational changes based on the review of our earlier work and recent literature. A Raman spectrum of a protein provides unique signatures for various secondary structures like helices, beta-sheets, turns, random structures, etc., and various amino acid residues such as tyrosine, tryptophan, and phenylalanine. We have shown how Raman spectra can differentiate between bovine serum albumin (BSA) and lysozyme protein based on their difference in sequence and structure (primary, secondary and tertiary). Although it is challenging to elucidate the structure of a protein using a Raman spectrum alone, Raman spectra can be used to differentiate small changes in conformations of proteins such as BSA during melting. Various new advancements in technique and data analyses in Raman spectroscopic studies of proteins have been discussed. The last part of the review focuses on the importance of the ROA spectrum to understand additional features about proteins. The ROA spectrum is rich in information about the protein backbone due to its rigidity compared to its side chains. Furthermore, the ROA spectra of lysozyme and BSA have been presented to show how ROA provides extra information about the solvent properties of proteins.

Protein secondary structure classification revisited: processing DSSP information with PSSC

A first step toward three-dimensional protein structure description is the characterization of secondary structure. The most widely used program for secondary structure assignment remains DSSP, introduced in 1983, with currently more than 400 citations per year. DSSP output is in a one-letter representation, where much of the information on DSSP's internal description is lost. Recently it became evident that DSSP overlooks most π-helical structures, which are more prevalent and important than anticipated before. We introduce an alternative concept, representing the internal structure characterization of DSSP as an eight-character string that is human-interpretable and easy to parse by software. We demonstrate how our protein secondary structure characterization (PSSC) code allows for inspection of complicated structural features. It recognizes ten times more π-helical residues than does the standard DSSP. The plausibility of introduced changes in interpreting DSSP information is demonstrated by better clustering of secondary structures in (φ, ψ) dihedral angle space. With a sliding sequence window (SSW), helical assignments with PSSC remain invariant compared with an assignment based on the complete structure. In contrast, assignment with DSSP can be changed by residues in the neighborhood that are in fact not interacting with the residue under consideration. We demonstrate how one can easily define new secondary structure classification schemes with PSSC and perform the classifications. Our approach works without changing the DSSP source code and allows for more detailed protein characterization.

Design of monomeric water-soluble β-hairpin and β-sheet peptides

Since the first report in 1993 (JACS 115, 5887-5888) of a peptide able to form a monomeric β-hairpin structure in aqueous solution, the design of peptides forming either β-hairpins (two-stranded antiparallel β-sheets) or three-stranded antiparallel β-sheets has become a field of growing interest and activity. These studies have yielded great insights into the principles governing the stability and folding of β-hairpins and antiparallel β-sheets. This chapter provides an overview of the reported β-hairpin/β-sheet peptides focussed on the applied design criteria, reviews briefly the factors contributing to β-hairpin/β-sheet stability, and describes a protocol for the de novo design of β-sheet-forming peptides based on them. Guidelines to select appropriate turn and strand residues and to avoid self-association are provided. The methods employed to check the success of new designed peptides are also summarized. Since NMR is the best technique to that end, NOEs and chemical shifts characteristic of β-hairpins and three-stranded antiparallel β-sheets are given.

β-Bulges: extensive structural analyses of β-sheets irregularities

β-Sheets are quite frequent in protein structures and are stabilized by regular main-chain hydrogen bond patterns. Irregularities in β-sheets, named β-bulges, are distorted regions between two consecutive hydrogen bonds. They disrupt the classical alternation of side chain direction and can alter the directionality of β-strands. They are implicated in protein-protein interactions and are introduced to avoid β-strand aggregation. Five different types of β-bulges are defined. Previous studies on β-bulges were performed on a limited number of protein structures or one specific family. These studies evoked a potential conservation during evolution. In this work, we analyze the β-bulge distribution and conservation in terms of local backbone conformations and amino acid composition. Our dataset consists of 66 times more β-bulges than the last systematic study (Chan et al. Protein Science 1993, 2:1574-1590). Novel amino acid preferences are underlined and local structure conformations are highlighted by the use of a structural alphabet. We observed that β-bulges are preferably localized at the N- and C-termini of β-strands, but contrary to the earlier studies, no significant conservation of β-bulges was observed among structural homologues. Displacement of β-bulges along the sequence was also investigated by Molecular Dynamics simulations.

Identifying the hotspots on the top faces of WD40-repeat proteins from their primary sequences by β-bulges and DHSW tetrads

The analysis of 36 available crystal structures of WD40 repeat proteins reveals widespread existence of a beta-bulge formed at the beginning of strand a and the end of strand b, termed as WD_b–a bulge: among a total of 259 WD40 blades, there are 243 such β-bulges. The R₁ positions in these WD_b–a bulges have fair distributions of Arg, His, Ile, Leu, Lys, Met, Phe, Trp, Tyr and Val residues. These residues protrude on the top face of the WD40 proteins and can serve as hotspots for protein-protein interactions. An analysis of 29 protein complexes formed by 17 WD proteins reveals that these R₁ residues, along with two other residues (R₁-2 and D-1), are indeed widely involved in protein-protein interactions. Interestingly, these WD_b–a bulges can be easily identified by the 4-amino acid sequences of (V, L, I), R₁, R₂, (V, L, I), along with some other significant amino acids. Thus, the hotspots of WD40 proteins on the top face can be readily predicted based on the primary sequences of the proteins. The literature-reported mutagenesis studies for Met30, MDV1, Tup11, COP1 and SPA1, which crystal structures are not available, can be readily understood based on the feature-based method. Applying the method, the twelve potential hotspots on the top face of Tup11 from S. japonicas have been identified. Our ITC measurements confirm seven of them, Tyr382, Arg284, Tyr426, Tyr508, Leu559, Lys575 and Ile601, are essential for recognizing Fep1. The ITC measurements further convinced that the feature-based method provides accurate prediction of hotspots on the top face.

Roles of beta-turns in protein folding: from peptide models to protein engineering

Reverse turns are a major class of protein secondary structure; they represent sites of chain reversal and thus sites where the globular character of a protein is created. It has been speculated for many years that turns may nucleate the formation of structure in protein folding, as their propensity to occur will favor the approximation of their flanking regions and their general tendency to be hydrophilic will favor their disposition at the solvent-accessible surface. Reverse turns are local features, and it is therefore not surprising that their structural properties have been extensively studied using peptide models. In this article, we review research on peptide models of turns to test the hypothesis that the propensities of turns to form in short peptides will relate to the roles of corresponding sequences in protein folding. Turns with significant stability as isolated entities should actively promote the folding of a protein, and by contrast, turn sequences that merely allow the chain to adopt conformations required for chain reversal are predicted to be passive in the folding mechanism. We discuss results of protein engineering studies of the roles of turn residues in folding mechanisms. Factors that correlate with the importance of turns in folding indeed include their intrinsic stability, as well as their topological context and their participation in hydrophobic networks within the protein's structure.

Structural malleability of plasticins: Preorganized conformations in solution and relevance for antimicrobial activity

Plasticins (23 long-residue glycine-leucine-rich dermaseptin-related peptides produced by the skin of South American hylids) have very similar amino acid sequences, hydrophobicities, and amphipathicities, but differ in their membrane-damaging properties and structurations (i.e. destabilized helix states, beta-hairpin, beta-sheet, and disordered states) at anionic and zwitterionic membrane interfaces. Structural malleability of plasticins in aqueous solutions together with parameters that may govern their ability to fold within beta-hairpin like structures were analyzed through circular dichroism and FTIR spectroscopic studies completed by molecular dynamics simulations in polar mimetic media. The goal of this study was to probe to which extent pre-existent peptide conformations, i.e. intrinsic "conformational landscape", may be responsible for variability in bioactive conformation and antimicrobial/hemolytic mechanisms of action of these peptides in relation with their various membrane disturbing properties. All plasticins present a turn region that does not always result in folding into a beta-hairpin shaped conformation. Residue at position 8 plays a major role in initiating the folding, while position 12 is not critical. Conformational stability has no major impact on antimicrobial efficacy. However, preformed beta-hairpin in solution may act as a conformational lock that prevents switch to alpha-helical structure. This lock lowers the antimicrobial efficiency and explains subtle differences in potencies of the most active antimicrobial plasticins.

Natural polypeptide scaffolds: beta-sheets, beta-turns, and beta-hairpins

This paper provides an introduction to fundamental conformational states of polypeptides in the beta-region of phi,psi space, in which the backbone is extended near to its maximal length, and to more complex architectures in which extended segments are linked by turns and loops. There are several variants on these conformations, and they comprise versatile scaffolds for presentation of side chains and backbone amides for molecular recognition and designed catalysts. In addition, the geometry of these fundamental folds can be readily mimicked in peptidomimetics.

Ubiquitin: a small protein folding paradigm

For the past twenty years, the small, 76-residue protein ubiquitin has been used as a model system to study protein structure, stability, folding and dynamics. In this time, ubiquitin has become a paradigm for both the experimental and computational folding communities. The folding energy landscape is now uniquely characterised with a plethora of information available on not only the native and denatured states, but partially structured states, alternatively folded states and locally unfolded states, in addition to the transition state ensemble. This Perspective focuses on the experimental characterisation of ubiquitin using a comprehensive range of biophysical techniques.

Reconstruction of the Conserved β-Bulge in Mammalian Defensins Using d-Amino Acids

Defensins are cationic antimicrobial mini-proteins that play important roles in the innate immune defense against microbial infection. Six invariant Cys residues in each defensin form three structurally indispensable intramolecular disulfide bridges. The only other residue invariant in all known mammalian defensins is a Gly. Structural studies indicate that the invariant Gly residue is located in an atypical, classic-type beta-bulge with the backbone torsion angles (Phi, Psi) disallowed for L-amino acids but permissible for D-enantiomers. We replaced the invariant Gly17 residue in human neutrophil alpha-defensin 2 (HNP2) by L-Ala or one of the D-amino acids Ala, Glu, Phe, Arg, Thr, Val, or Tyr. Although L-Ala17-HNP2 could not be folded, resulting in massive aggregation, all of the D-amino acid-substituted analogs folded with high efficiency. The high resolution x-ray crystal structures of dimeric D-Ala17-HNP2 were determined in three different crystal forms, showing a well preserved beta-bulge identical to those found in other defensins. The seven D-analogs of HNP2 exhibited highly variable bactericidal activity against Gram-positive and Gram-negative test strains, consistent with the premise that interplay between charge and hydrophobicity dictates how amphiphilic defensins kill. Further, the bactericidal activity of these d-amino acid analogs of HNP2 correlated well with their ability to induce leakage from large unilamellar vesicles, supporting membrane permeabilization as the lethal event in microbial killing by HNP2. Our findings identify a conformational prerequisite in the beta-bulge of defensins essential for correct folding and native structure, thereby explaining the molecular basis of the Gly-Xaa-Cys motif conserved in all mammalian defensins.

Measuring the refolding of beta-sheets with different turn sequences on a nanosecond time scale

Whether turns play an active or passive role in protein folding remains a controversial issue at this juncture. Here we use a photolabile cage strategy in combination with laser-flash photolysis and photoacoustic calorimetry to study the effects of different turns on the kinetics of beta-hairpin refolding on a nanosecond time scale. This strategy opens up a temporal window to allow the observation of early kinetic events in the protein refolding process at ambient temperature and pH without interference from any denaturants. Our results provide direct evidence demonstrating that even a one-residue difference in the turn region can change the refolding kinetics of a peptide. This observation suggests an active role for turn formation in directing protein folding.

Discerning the Structure and Energy of Multiple Transition States in Protein Folding using ψ-Analysis

We quantify the degree to which folding occurs along a complex landscape with structurally distinct pathways using psi-analysis in combination with a protein engineering method that identifies native, non-covalent polypeptide interactions and their relative populations at the rate-limiting step. By probing the proximity of two specific partners, this method is extremely well-suited for comparison to theoretical simulations. Using ubiquitin as a model system, we detect individual pathways with site-resolved resolution, demonstrating that the protein folds through a native-like transition state ensemble with a common nucleus that contains heterogeneous features on its periphery. The consensus transition state topology has part of the major helix docked against four properly aligned beta-strands. However, structural heterogeneity exists in the transition state ensemble, wherein peripheral regions are differentially populated according to their relative stability. Pathway diversity reflects the variable order of formation of these peripheral elements, which radiate outward from the common nucleus. These results, which show only moderate agreement with traditional mutational phi-analysis, provide an extraordinarily detailed and quantitative description of protein folding.

Local sequence information in cellular retinoic acid-binding protein I: Specific residue roles in β-turns*

We have recently shown that two of the beta-turns (III and IV) in the ten-stranded, beta-clam protein, cellular retinoic acid-binding protein I (CRABP I), are favored in short peptide fragments, arguing that they are encoded by local interactions (K. S. Rotondi and L. M. Gierasch, Biochemistry, 2003, Vol. 42, pp. 7976-7985). In this paper we examine these turns in greater detail to dissect the specific local interactions responsible for their observed native conformational biases. Conformations of peptides corresponding to the turn III and IV fragments were examined under conditions designed to selectively disrupt stabilizing interactions, using pH variation, chaotrope addition, or mutagenesis to probe specific side-chain influences. We find that steric constraints imposed by excluded volume effects between near neighbor residues (i,i+2), favorable polar (i,i+2) interactions, and steric permissiveness of glycines are the principal factors accounting for the observed native bias in these turns. Longer-range stabilizing interactions across the beta-turns do not appear to play a significant role in turn stability in these short peptides, in contrast to their importance in hairpins. Additionally, our data add to a growing number of examples of the 3:5 type I turn with a beta-bulge as a class of turns with high propensity to form locally defined structure. Current work is directed at the interplay between the local sequence information in the turns and more long-range influences in the mechanism of folding of this predominantly beta-sheet protein.

Expected and unexpected results from combined β-hairpin design elements

A model beta-hairpin dodecapeptide [EFGWVpGKWTIK] was designed by including a favorable D-ProGly Type II' beta-turn sequence and a Trp-zip interaction, while also incorporating a beta-strand unfavorable glycine residue in the N-terminal strand. This peptide is highly folded and monomeric in aqueous solution as determined by combined analysis with circular dichroism and 1H NMR spectroscopy. A peptide representing the folded conformation of the model beta-hairpin [cyclic(EFGWVpGKWTIKpG)] and a linear peptide representing the unfolded conformation [EFGWVPGKWTIK] yield unexpected relative deviations between the CD and 1H NMR spectroscopic results that are attributed to variations in the packing interactions of the aromatic side chains. Mutational analysis of the model beta-hairpin indicates that the Trp-zip interaction favors folding and stability relative to an alternate hydrophobic cluster between Trp and Tyr residues [EFGYVpGKWTIK]. The significance of select diagonal interactions in the model beta-hairpin was tested by rearranging the cross-strand hydrophobic interactions to provide a folded peptide [EWFGIpGKTYWK] displaying evidence of an unusual backbone conformation at the hydrophobic cluster. This unusual conformation does not appear to be a result of the glycine residue in the beta-strand, as replacement with a serine results in a peptide [EWFSIpGKTYWK] with a similar and seemingly characteristic CD spectrum. However, an alternate arrangement of hydrophobic residues with a Trp-zip interaction in a similar position to the parent beta-hairpin [EGFWVpGKWITK] results in a folded beta-hairpin conformation. The differences between side chain packing of these peptides precludes meaningful thermodynamic analysis and illustrates the caution necessary when interpreting beta-hairpin folding thermodynamics that are driven, at least in part, by aromatic cross strand interactions.

Native structural propensity in cellular retinoic acid-binding protein I 64-88: the role of locally encoded structure in the folding of a beta-barrel protein

A central question in protein folding is the relative importance of locally encoded structure and cooperative interactions among residues distant in sequence. We have been exploring this question in a predominantly beta-sheet protein, since beta-structure formation clearly relies on both local and global sequence information. We present evidence that a 24-residue peptide corresponding to two linked hairpins of cellular retinoic acid-binding protein I (CRABP I) adopts significant native structure in aqueous solution. Prior work from our laboratory showed that the two turns contained in this fragment (turns III and IV) had the highest tendency of any of the eight turns in this anti-parallel beta-barrel to fold into native turns. In addition, the primary sequence of these two turns is well conserved throughout the structural family to which CRABP I belongs, and residues in the turns and their associated hairpins participate in a network of conserved long-range interactions. We propose that the strong local-sequence biases within the chain segment comprising turns III and IV favor longer-range interactions that are crucial to the folding and native-state stability of CRABP I, and may play a similar role in related intracellular lipid-binding proteins (iLBPs).