Guidelines for protein design: the energetics of beta sheet side chain interactions

In silico protein dynamics in the human cytoplasm: Partial folding, misfolding, fold switching, and non-native interactions

We examine the influence of cellular interactions in all-atom models of a section of the Homo sapiens cytoplasm on the early folding events of the three-helix bundle protein B (PB). While genetically engineered PB is known to fold in dilute water box simulations in three microseconds, the three initially unfolded PB copies in our two cytoplasm models using a similar force field did not reach the native state during 30-microsecond simulations. We did however capture the formation of all three helices in a compact native-like topology. Folding in vivo is delayed because intramolecular contact formation within PB is in direct competition with intermolecular contacts between PB and surrounding macromolecules. In extreme cases, intermolecular beta-sheets are formed. Interactions with other macromolecules are also observed to promote structure formation, for example when a PB helix in our simulations is shielded from solvent by macromolecular crowding. Sticking and crowding in our models initiate sampling of helix/sheet structural plasticity of PB. Relatedly, in past in vitro experiments, similar GA domains were shown to switch between two different folds. Finally, we also observed that stickiness between PB and the cellular environment can be modulated in our simulations through the reduction in protein hydrophobicity when we reversed PB back to the wild-type sequence. This study demonstrates that even fast-folding proteins can get stuck in non-native states in the cell, making them useful models for protein-chaperone interactions and early stages of aggregate formation relevant to cellular disease.

Anticancer peptides mechanisms, simple and complex

Peptide therapy has started since 1920s with the advent of insulin application, and now it has emerged as a new approach in treatment of diseases including cancer. Using anti-cancer peptides (ACPs) is a promising way of cancer therapy as ACPs are continuing to be approved and arrived at major pharmaceutical markets. Traditional cancer treatments face different problems like intensive adverse effects to patient's body, cell resistance to conventional chemical drugs and in some worse cases the occurrence of cell multidrug resistance (MDR) of cancerous tissues against chemotherapy. On the other hand, there are some benefits conceived for peptides usage in treatment of diseases specifically cancer, as these compounds present favorable characteristics such as smaller size, high activity, low immunogenicity, good biocompatibility in vivo, convenient and rapid way of synthesis, amenable to sequence modification and revision and there is no limitation for the type of cargo they carry. It is possible to achieve an optimum molecular and functional structure of peptides based on previous experience and bank of peptide motif data which may result in novel peptide design. Bioactive peptides are able to form pores in cell membrane and induce necrosis or apoptosis of abnormal cells. Moreover, recent researches have focused on the tumor recognizing peptide motifs with the ability to permeate to cancerous cells with the aim of cancer treatment at earlier stages. In this strategy the most important factors for addressing cancer are choosing peptides with easy accessibility to tumor cell without cytotoxicity effect towards normal cells. The peptides must also meet acceptable pharmacokinetic requirements. In this review, the characteristics of peptides and cancer cells are discussed. The various mechanisms of peptides' action proposed against cancer cells make the next part of discussion. It will be followed by giving information on peptides application, various methods of peptide designing along with introducing various databases. Future aspects of peptides for employing in area of cancer treatment come as conclusion at the end.

Single-crystal-to-single-crystal translation of a helical supramolecular polymer to a helical covalent polymer

Polymers possessing helical conformation in the solid state are in high demand. We report a helical peptide-polymer via the topochemical ene-azide cycloaddition (TEAC) polymerization. The molecules of the designed Gly-Phe-based dipeptide, decorated with ene and azide, assemble in its crystals as β-sheets and as supramolecular helices in two mutually perpendicular directions. While the NH…O H-bonding facilitates β-sheet-like stacking along one direction, weak CH…N H-bonding between the azide-nitrogen and vinylic-hydrogen of molecules belonging to the adjacent stacks arranges them in a head-to-tail manner as supramolecular helices. In the crystal lattice, the azide and alkene of adjacent molecules in the supramolecular helix are suitably preorganized for their TEAC reaction. The dipeptide underwent regio- and stereospecific polymerization upon mild heating in a single-crystal-to-single-crystal fashion, yielding a triazoline-linked helical covalent polymer that could be characterized by single-crystal X-ray diffraction studies. Upon heating, the triazoline-linked polymer undergoes denitrogenation to aziridine-linked polymer, as evidenced by differential scanning calorimetry, thermogravimetric analysis, and solid-state NMR analyses.

The Effects of Charged Amino Acid Side-Chain Length on Diagonal Cross-Strand Interactions between Carboxylate- and Ammonium-Containing Residues in a β-Hairpin

The β-sheet is one of the common protein secondary structures, and the aberrant aggregation of β-sheets is implicated in various neurodegenerative diseases. Cross-strand interactions are an important determinant of β-sheet stability. Accordingly, both diagonal and lateral cross-strand interactions have been studied. Surprisingly, diagonal cross-strand ion-pairing interactions have yet to be investigated. Herein, we present a systematic study on the effects of charged amino acid side-chain length on a diagonal ion-pairing interaction between carboxylate- and ammonium-containing residues in a β-hairpin. To this end, 2D-NMR was used to investigate the conformation of the peptides. The fraction folded population and the folding free energy were derived from the chemical shift data. The fraction folded population for these peptides with potential diagonal ion pairs was mostly lower compared to the corresponding peptide with a potential lateral ion pair. The diagonal ion-pairing interaction energy was derived using double mutant cycle analysis. The Asp2-Dab9 (Asp: one methylene; Dab: two methylenes) interaction was the most stabilizing (−0.79 ± 0.14 kcal/mol), most likely representing an optimal balance between the entropic penalty to enable the ion-pairing interaction and the number of side-chain conformations that can accommodate the interaction. These results should be useful for designing β-sheet containing molecular entities for various applications.

Swapping the Positions in a Cross-Strand Lateral Ion-Pairing Interaction between Ammonium- and Carboxylate-Containing Residues in a β-Hairpin

Cross-strand lateral ion-pairing interactions are important for antiparallel β-sheet stability. Statistical studies suggested that swapping the position of cross-strand lateral residues should not significantly affect the interaction. Herein, we swapped the position of ammonium- and carboxylate-containing residues with different side-chain lengths in a cross-strand lateral ion-pairing interaction in a β-hairpin. The peptides were analyzed by 2D-NMR. The fraction folded population and folding free energy were derived from the chemical shift data. The ion-pairing interaction energy was derived using double mutant cycle analysis. The general trends for the fraction folded population and interaction energetics remained similar upon swapping the position of the interacting charged residues. The most stabilizing cross-strand interactions were between short residues, similar to the unswapped study. However, the fraction folded populations for most of the swapped peptides were higher compared to the corresponding unswapped peptides. Furthermore, subtle differences in the ion-pairing interaction energy upon swapping were observed, most likely due to the “unleveled” relative positioning of the interacting residues created by the inherent right-handed twist of the structure. These results should be useful for developing functional peptides that rely on lateral ion-pairing interactions across antiparallel β-strands.

Engineered Reader Proteins for Enhanced Detection of Methylated Lysine on Histones

Histone post-translational modifications (PTMs) are crucial for many cellular processes including mitosis, transcription, and DNA repair. The cellular readout of histone PTMs is dependent on both the chemical modification and histone site, and the array of histone PTMs on chromatin is dynamic throughout the eukaryotic life cycle. Accordingly, methods that report on the presence of PTMs are essential tools for resolving open questions about epigenetic processes and for developing therapeutic diagnostics. Reader domains that recognize histone PTMs have shown potential as advantageous substitutes for anti-PTM antibodies, and engineering efforts aimed at enhancing reader domain affinities would advance their efficacy as antibody alternatives. Here we describe engineered chromodomains from Drosophila melanogaster and humans that bind more tightly to H3K9 methylation (H3K9me) marks and result in the tightest reported reader-H3K9me interaction to date. Point mutations near the binding interface of the HP1 chromodomain were screened in a combinatorial fashion, and a triple mutant was found that binds 20-fold tighter than the native scaffold without any loss in PTM-site selectivity. The beneficial mutations were then translated to a human homologue, CBX1, resulting in an even tighter interaction with H3K9me3. Furthermore, we show that these engineered readers (eReaders) increase detection of H3K9me marks in several biochemical assays and outperform a commercial anti-H3K9me antibody in detecting H3K9me-containing nucleosomes in vitro, demonstrating the utility of eReaders to complement antibodies in epigenetics research.

The relative stability of trpzip1 and its mutants determined by computation and experiment

The single-point mutations of tprzip1 are indicated at left, and their relative energetics are compared at right.

Prions and Prion-like assemblies in neurodegeneration and immunity: The emergence of universal mechanisms across health and disease

Prion-like behaviour is an abrupt process, an "all-or-nothing" transition between a monomeric species and an "infinite" fibrillated form. Once a nucleation point is formed, the process is unstoppable as fibrils self-propagate by recruiting and converting all monomers into the amyloid fold. After the "mad cow" episode, prion diseases have made the headlines, but more and more prion-like behaviours have emerged in neurodegenerative diseases, where formation of fibrils and large conglomerates of proteins deeply disrupt the cell homeostasis. More interestingly, in the last decade, examples emerged to suggest that prion-like conversion can be used as a positive gain of function, for memory storage or structural scaffolding. More recent experiments show that we are only seeing the tip of the iceberg and that, for example, prion-like amplification is found in many pathways of the immune response. In innate immunity, receptors on the cellular surface or within the cells 'sense' danger and propagate this information as signal, through protein-protein interactions (PPIs) between 'receptor', 'adaptor' and 'effector' proteins. In innate immunity, the smallest signal of a foreign element or pathogen needs to trigger a macroscopic signal output, and it was found that adaptor polymerize to create an extreme signal amplification. Interestingly, our body uses multiple structural motifs to create large signalling platform; a few innate proteins use amyloid scaffolds but most of the polymers discovered are composed by self-assembly in helical filaments. Some of these helical assemblies even have intercellular "contamination" in a "true" prion action, as demonstrated for ASC specks and MyD88 filaments. Here, we will describe the current knowledge in neurodegenerative diseases and innate immunity and show how these two very different fields can cross-seed discoveries.

Self-assembly of photochromic diarylethene-peptide conjugates stabilized by β-sheet formation at the liquid/graphite interface

Two-dimensional (2-D) self-assembly of diarylethene (DAE)-peptide conjugates at the octanoic acid/graphite interface was investigated by scanning tunnelling microscopy (STM). The open-ring isomer of a DAE-peptide conjugate formed a stable 2-D molecular assembly with an antiparallel β-sheet structure. Quantitative analysis of surface coverage depending on concentration revealed a stronger stabilization effect of the oligopeptide than that of the alkyl group with a similar side chain length.

In vivo localization of human acetylcholinesterase-derived species in a β-sheet conformation at the core of senile plaques in Alzheimer's disease

Many neurodegenerative diseases are characterized by amyloid deposition. In Alzheimer's disease (AD), β-amyloid (Aβ) peptides accumulate extracellularly in senile plaques. The AD amyloid cascade hypothesis proposes that Aβ production or reduced clearance leads to toxicity. In contrast, the cholinergic hypothesis argues for a specific pathology of brain cholinergic pathways. However, neither hypothesis in isolation explains the pattern of AD pathogenesis. Evidence suggests that a connection exists between these two scenarios: the synaptic form of human acetylcholinesterase (hAChE-S) associates with plaques in AD brains; among hAChE variants, only hAChE-S enhances Aβ fibrillization in vitro and Aβ deposition and toxicity in vivo Only hAChE-S contains an amphiphilic C-terminal domain (T40, AChE_575-614), with AChE_586-599 homologous to Aβ and forming amyloid fibrils, which implicates T40 in AD pathology. We previously showed that the amyloid scavenger, insulin-degrading enzyme (IDE), generates T40-derived amyloidogenic species that, as a peptide mixture, seed Aβ fibrillization. Here, we characterized 11 peptides from a T40-IDE digest for β-sheet conformation, surfactant activity, fibrillization, and seeding capability. We identified residues important for amyloidogenicity and raised polyclonal antibodies against the most amyloidogenic peptide. These new antisera, alongside other specific antibodies, labeled sections from control, hAChE-S, hAPPswe, and hAChE-S/hAPPswe transgenic mice. We observed that hAChE-S β-sheet species co-localized with Aβ in mature plaque cores, surrounded by hAChE-S α-helical species. This observation provides the first in vivo evidence of the conformation of hAChE-S species within plaques. Our results may explain the role of hAChE-S in Aβ deposition and aggregation, as amyloidogenic hAChE-S β-sheet species might seed Aβ aggregation.

Principles of Inter-Amino-Acid Recognition Revealed by Binding Energies between Homogeneous Oligopeptides

We have determined the interaction strengths of the common naturally occurring amino acids using a complete binding affinity matrix of 20 × 20 pairs of homo-octapeptides consisting of the 20 common amino acids between stationary and mobile states. We used a bead-based fluorescence assay for these measurements. The results provide a basis for analyzing specificity, polymorphisms, and selectivity of inter-amino-acid interactions. Comparative analyses of the binding energies, i.e., the free energies of association (ΔG A), reveal contributions assignable to both main-chain-related and side-chain-related interactions originating from the chemical structures of these 20 common amino acids. Side-chain-side-chain and side-chain-main-chain interactions are found to be pronounced in an identified set of amino acid pairs that determine the basis of inter-amino-acid recognition.

BetaDL: A protein beta-sheet predictor utilizing a deep learning model and independent set solution

The sequence-based prediction of beta-residue contacts and beta-sheet structures contain key information for protein structure prediction. However, the determination of beta-sheet structures poses numerous challenges due to long-range beta-residue interactions and the huge number of possible beta-sheet structures. Recently gaining attention has been the prediction of residue contacts based on deep learning models whose results have led to improvement in protein structure prediction. In addition, to reduce the computational complexity of determining beta-sheet structures, it has been suggested that this problem be transformed into graph-based solutions. Consequently, the current work proposes BetaDL, a combination of a deep learning and a graph-based beta-sheet structure predictor. BetaDL adopts deep learning models to capture beta-residue contacts and improve beta-sheet structure predictions. In addition, a graph-based approach is presented to model the beta-sheets conformational space and a new score function is introduced to evaluate beta-sheets. Furthermore, the present study demonstrates that the beta-sheet structure can be predicted within an acceptable computational time by the utilization of a heuristic maximum weight independent set solution. When compared to state-of-the-art methods, experimental results from BetaSheet916 and BetaSheet1452 datasets indicate that BetaDL improves the accuracy of beta-residue contact and beta-sheet structure prediction. Using BetaDL, beta-sheet structures are predicted with a 4% and 6% improvement in the F1-score at the residue and strand levels, respectively. BetaDL's source code and data are available at http://kerg.um.ac.ir/index.php/datasets/#BetaDL.

A threonine zipper that mediates protein-protein interactions: Structure and prediction

We present the structure of an engineered protein-protein interface between two beta barrel proteins, which is mediated by interactions between threonine (Thr) residues. This Thr zipper structure suggests that the protein interface is stabilized by close-packing of the Thr residues, with only one intermonomer hydrogen bond (H-bond) between two of the Thr residues. This Thr-rich interface provides a unique opportunity to study the behavior of Thr in the context of many other Thr residues. In previous work, we have shown that the side chain (χ1 ) dihedral angles of interface and core Thr residues can be predicted with high accuracy using a hard sphere plus stereochemical constraint (HS) model. Here, we demonstrate that in the Thr-rich local environment of the Thr zipper structure, we are able to predict the χ1 dihedral angles of most of the Thr residues. Some, however, are not well predicted by the HS model. We therefore employed explicitly solvated molecular dynamics (MD) simulations to further investigate the side chain conformations of these residues. The MD simulations illustrate the role that transient H-bonding to water, in combination with steric constraints, plays in determining the behavior of these Thr side chains. Broader Audience Statement: Protein-protein interactions are critical to life and the search for ways to disrupt adverse protein-protein interactions involved in disease is an ongoing area of drug discovery. We must better understand protein-protein interfaces, both to be able to disrupt existing ones and to engineer new ones for a variety of biotechnological applications. We have discovered and characterized an artificial Thr-rich protein-protein interface. This novel interface demonstrates a heretofore unknown property of Thr-rich surfaces: mediating protein-protein interactions.

Packing in protein cores

Proteins are biological polymers that underlie all cellular functions. The first high-resolution protein structures were determined by x-ray crystallography in the 1960s. Since then, there has been continued interest in understanding and predicting protein structure and stability. It is well-established that a large contribution to protein stability originates from the sequestration from solvent of hydrophobic residues in the protein core. How are such hydrophobic residues arranged in the core; how can one best model the packing of these residues, and are residues loosely packed with multiple allowed side chain conformations or densely packed with a single allowed side chain conformation? Here we show that to properly model the packing of residues in protein cores it is essential that amino acids are represented by appropriately calibrated atom sizes, and that hydrogen atoms are explicitly included. We show that protein cores possess a packing fraction of [Formula: see text], which is significantly less than the typically quoted value of 0.74 obtained using the extended atom representation. We also compare the results for the packing of amino acids in protein cores to results obtained for jammed packings from discrete element simulations of spheres, elongated particles, and composite particles with bumpy surfaces. We show that amino acids in protein cores pack as densely as disordered jammed packings of particles with similar values for the aspect ratio and bumpiness as found for amino acids. Knowing the structural properties of protein cores is of both fundamental and practical importance. Practically, it enables the assessment of changes in the structure and stability of proteins arising from amino acid mutations (such as those identified as a result of the massive human genome sequencing efforts) and the design of new folded, stable proteins and protein-protein interactions with tunable specificity and affinity.

A cell is more than the sum of its (dilute) parts: A brief history of quinary structure

Most knowledge of protein structure and function is derived from experiments performed with purified protein resuspended in dilute, buffered solutions. However, proteins function in the crowded, complex cellular environment. Although the first four levels of protein structure provide important information, a complete understanding requires consideration of quinary structure. Quinary structure comprises the transient interactions between macromolecules that provides organization and compartmentalization inside cells. We review the history of quinary structure in the context of several metabolic pathways, and the technological advances that have yielded recent insight into protein behavior in living cells. The evidence demonstrates that protein behavior in isolated solutions deviates from behavior in the physiological environment.

Miniaturizing VEGF: Peptides mimicking the discontinuous VEGF receptor-binding site modulate the angiogenic response

The angiogenic properties of VEGF are mediated through the binding of VEGF to its receptor VEGFR2. The VEGF/VEGFR interface is constituted by a discontinuous binding region distributed on both VEGF monomers. We attempted to reproduce this discontinuous binding site by covalently linking into a single molecular entity two VEGF segments involved in receptor recognition. We designed and synthesized by chemical ligation a set of peptides differing in length and flexibility of the molecular linker joining the two VEGF segments. The biological activity of the peptides was characterized in vitro and in vivo showing a VEGF-like activity. The most biologically active mini-VEGF was further analyzed by NMR to determine the atomic details of its interaction with the receptor.

Protein design: Past, present, and future

Building on the pioneering work of Ho and DeGrado (J Am Chem Soc 1987, 109, 6751-6758) in the late 1980s, protein design approaches have revealed many fundamental features of protein structure and stability. We are now in the era that the early work presaged - the design of new proteins with practical applications and uses. Here we briefly survey some past milestones in protein design, in addition to highlighting recent progress and future aspirations.

Boosting protein stability with the computational design of β-sheet surfaces

β-sheets often have one face packed against the core of the protein and the other facing solvent. Mutational studies have indicated that the solvent-facing residues can contribute significantly to protein stability, and that the preferred amino acid at each sequence position is dependent on the precise structure of the protein backbone and the identity of the neighboring amino acids. This suggests that the most advantageous methods for designing β-sheet surfaces will be approaches that take into account the multiple energetic factors at play including side chain rotamer preferences, van der Waals forces, electrostatics, and desolvation effects. Here, we show that the protein design software Rosetta, which models these energetic factors, can be used to dramatically increase protein stability by optimizing interactions on the surfaces of small β-sheet proteins. Two design variants of the β-sandwich protein from tenascin were made with 7 and 14 mutations respectively on its β-sheet surfaces. These changes raised the thermal midpoint for unfolding from 45°C to 64°C and 74°C. Additionally, we tested an empirical approach based on increasing the number of potential salt bridges on the surfaces of the β-sheets. This was not a robust strategy for increasing stability, as three of the four variants tested were unfolded.

Physicochemical Properties of Ion Pairs of Biological Macromolecules

Ion pairs (also known as salt bridges) of electrostatically interacting cationic and anionic moieties are important for proteins and nucleic acids to perform their function. Although numerous three-dimensional structures show ion pairs at functionally important sites of biological macromolecules and their complexes, the physicochemical properties of the ion pairs are not well understood. Crystal structures typically show a single state for each ion pair. However, recent studies have revealed the dynamic nature of the ion pairs of the biological macromolecules. Biomolecular ion pairs undergo dynamic transitions between distinct states in which the charged moieties are either in direct contact or separated by water. This dynamic behavior is reasonable in light of the fundamental concepts that were established for small ions over the last century. In this review, we introduce the physicochemical concepts relevant to the ion pairs and provide an overview of the recent advancement in biophysical research on the ion pairs of biological macromolecules.

Aromatic cluster mutations produce focal modulations of β-sheet structure

Site-directed mutagenesis is a powerful tool for altering the structure and function of proteins in a focused manner. Here, we examined how a model β-sheet protein could be tuned by mutation of numerous surface-exposed residues to aromatic amino acids. We designed these aromatic side chain "clusters" at highly solvent-exposed positions in the flat, single-layer β-sheet of Borrelia outer surface protein A (OspA). This unusual β-sheet scaffold allows us to interrogate the effects of these mutations in the context of well-defined structure but in the absence of the strong scaffolding effects of globular protein architecture. We anticipated that the introduction of a cluster of aromatic amino acid residues on the β-sheet surface would result in large conformational changes and/or stabilization and thereby provide new means of controlling the properties of β-sheets. Surprisingly, X-ray crystal structures revealed that the introduction of aromatic clusters produced only subtle conformational changes in the OspA β-sheet. Additionally, despite burying a large degree of hydrophobic surface area, the aromatic cluster mutants were slightly less stable than the wild-type scaffold. These results thereby demonstrate that the introduction of aromatic cluster mutations can serve as a means for subtly modulating β-sheet conformation in protein design.

Effect of charged amino acid side chain length on lateral cross-strand interactions between carboxylate- and guanidinium-containing residues in a β-hairpin

β-Sheet is one of the major protein secondary structures. Oppositely charged residues are frequently observed across neighboring strands in antiparallel sheets, suggesting the importance of cross-strand ion pairing interactions. The charged amino acids Asp, Glu, Arg, and Lys have different numbers of hydrophobic methylenes linking the charged functionality to the backbone. To investigate the effect of side chain length of guanidinium- and carboxylate-containing residues on lateral cross-strand ion pairing interactions at non-hydrogen-bonded positions, β-hairpin peptides containing Zbb-Agx (Zbb = Asp, Glu, Aad in increasing length; Agx = Agh, Arg, Agb, Agp in decreasing length) sequence patterns were studied by NMR methods. The fraction folded population and folding energy were derived from the chemical shift deviation data. Peptides with high fraction folded populations involved charged residue side chain lengths that supported high strand propensity. Double mutant cycle analysis was used to determine the interaction energy for the potential lateral ion pairs. Minimal interaction was observed between residues with short side chains, most likely due to the diffused positive charge on the guanidinium group, which weakened cross-strand electrostatic interactions with the carboxylate side chain. Only the Aad-Arg/Agh interactions with long side chains clearly exhibited stabilizing energetics, possibly relying on hydrophobics. A survey of a non-redundant protein structure database revealed that the statistical sheet pair propensity followed the trend Asp-Arg < Glu-Arg, implying the need for matching long side chains. This suggested the need for long side chains on both guanidinium-bearing and carboxylate-bearing residues to stabilize the β-hairpin motif.

bbcontacts: prediction of β-strand pairing from direct coupling patterns

MOTIVATION

It has recently become possible to build reliable de novo models of proteins if a multiple sequence alignment (MSA) of at least 1000 homologous sequences can be built. Methods of global statistical network analysis can explain the observed correlations between columns in the MSA by a small set of directly coupled pairs of columns. Strong couplings are indicative of residue-residue contacts, and from the predicted contacts a structure can be computed. Here, we exploit the structural regularity of paired β-strands that leads to characteristic patterns in the noisy matrices of couplings. The β-β contacts should be detected more reliably than single contacts, reducing the required number of sequences in the MSAs.

RESULTS

bbcontacts predicts β-β contacts by detecting these characteristic patterns in the 2D map of coupling scores using two hidden Markov models (HMMs), one for parallel and one for antiparallel contacts. β-bulges are modelled as indel states. In contrast to existing methods, bbcontacts uses predicted instead of true secondary structure. On a standard set of 916 test proteins, 34% of which have MSAs with < 1000 sequences, bbcontacts achieves 50% precision for contacting β-β residue pairs at 50% recall using predicted secondary structure and 64% precision at 64% recall using true secondary structure, while existing tools achieve around 45% precision at 45% recall using true secondary structure.

AVAILABILITY AND IMPLEMENTATION

bbcontacts is open source software (GNU Affero GPL v3) available at https://bitbucket.org/soedinglab/bbcontacts .

A peptide from human semenogelin I self-assembles into a pH-responsive hydrogel

The peptide GSFSIQYTYHV derived from human semenogelin I forms a transparent hydrogel through spontaneous self-assembly in water at neutral pH. Linear rheology measurements demonstrate that the gel shows a dominating elastic response over a large frequency interval. CD, fluorescence and FTIR spectroscopy and cryo-TEM studies imply long fibrillar aggregates of extended β-sheet. Dynamic light scattering data indicate that the fibril lengths are of the order of micrometers. Time-dependent thioflavin T fluorescence shows that fibril formation by GSFSIQYTYHV is a nucleated reaction. The peptide may serve as basis for development of smart biomaterials of low immunogenicity suitable for biomedical applications, including drug delivery and wound healing.

Anatomy of β-strands at protein-protein interfaces

The development of inhibitors for protein-protein interactions frequently involves the mimicry of secondary structure motifs. While helical protein-protein interactions have been heavily targeted, a similar level of success for the inhibition of β-strand and β-sheet rich interfaces has been elusive. We describe an assessment of the full range of β-strand interfaces whose high-resolution structures are available in the Protein Data Bank. This analysis identifies complexes where a β-stand or β-sheet contributes significantly to binding. The results highlight the molecular recognition complexity in strand-mediated interactions relative to helical interfaces and offer guidelines for the construction of β-strand and β-sheet mimics as ligands for protein receptors. The online data set will potentially serve as an entry-point to new classes of protein-protein interaction inhibitors.

Intrinsic α-helical and β-sheet conformational preferences: a computational case study of alanine

A fundamental question in protein science is what is the intrinsic propensity for an amino acid to be in an α-helix, β-sheet, or other backbone dihedral angle ( ϕ-ψ) conformation. This question has been hotly debated for many years because including all protein crystal structures from the protein database, increases the probabilities for α-helical structures, while experiments on small peptides observe that β-sheet-like conformations predominate. We perform molecular dynamics (MD) simulations of a hard-sphere model for Ala dipeptide mimetics that includes steric interactions between nonbonded atoms and bond length and angle constraints with the goal of evaluating the role of steric interactions in determining protein backbone conformational preferences. We find four key results. For the hard-sphere MD simulations, we show that (1) β-sheet structures are roughly three and half times more probable than α-helical structures, (2) transitions between α-helix and β-sheet structures only occur when the backbone bond angle τ (NCα C) is greater than 110°, and (3) the probability distribution of τ for Ala conformations in the "bridge" region of ϕ-ψ space is shifted to larger angles compared to other regions. In contrast, (4) the distributions obtained from Amber and CHARMM MD simulations in the bridge regions are broader and have increased τ compared to those for hard sphere simulations and from high-resolution protein crystal structures. Our results emphasize the importance of hard-sphere interactions and local stereochemical constraints that yield strong correlations between ϕ-ψ conformations and τ.

Characterization of intermolecular structure of β(2)-microglobulin core fragments in amyloid fibrils by vacuum-ultraviolet circular dichroism spectroscopy and circular dichroism theory

Intermolecular structures are important factors for understanding the conformational properties of amyloid fibrils. In this study, vacuum-ultraviolet circular dichroism (VUVCD) spectroscopy and circular dichroism (CD) theory were used for characterizing the intermolecular structures of β2-microglobulin (β2m) core fragments in the amyloid fibrils. The VUVCD spectra of β2m20-41, β2m21-31, and β2m21-29 fragments in the amyloid fibrils exhibited characteristic features, but they were affected not only by the backbone conformations but also by the aromatic side-chain conformations. To estimate the contributions of aromatic side-chains to the spectra, the theoretical spectra were calculated from the simulated structures of β2m21-29 amyloid fibrils with various types of β-sheet stacking (parallel or antiparallel) using CD theory. We found that the experimental spectrum of β2m21-29 fibrils is largely affected by aromatic-backbone couplings, which are induced by the interaction between transitions within the aromatic and backbone chromophores, and these couplings are sensitive to the type of stacking among the β-sheets of the fibrils. Further theoretical analyses of simulated structures incorporating mutated aromatic residues suggested that the β2m21-29 fibrils are composed of amyloid accumulations in which the parallel β-sheets stack in an antiparallel manner and that the characteristic Phe-Tyr interactions among the β-sheet stacks affect the aromatic-backbone coupling. These findings indicate that the coupling components, which depend on the characteristic intermolecular structures, induce the spectral differences among three fragments in the amyloid fibrils. These advanced spectral analyses using CD theory provide a useful method for characterizing the intermolecular structures of protein and peptide fragment complexes.

Effect of charged amino acid side chain length on lateral cross-strand interactions between carboxylate-containing residues and lysine analogues in a β-hairpin

β-Sheets are one of the fundamental three-dimensional building blocks for protein structures. Oppositely charged amino acids are frequently observed directly across one another in antiparallel sheet structures, suggesting the importance of cross-strand ion pairing interactions. Despite the apparent electrostatic nature of ion pairing interactions, the charged amino acids Asp, Glu, Arg, Lys have different numbers of hydrophobic methylenes linking the charged functionality to the backbone. Accordingly, the effect of charged amino acid side chain length on cross-strand ion pairing interactions at lateral non-hydrogen bonded positions was investigated in a β-hairpin motif. The negatively charged residues with a carboxylate (Asp, Glu, Aad in increasing length) were incorporated at position 4, and the positively charged residues with an ammonium (Dap, Dab, Orn, Lys in increasing length) were incorporated at position 9. The fraction folded population and folding free energy were derived from the chemical shift deviation data. Double mutant cycle analysis was used to determine the interaction energy for the potential lateral ion pairs. Only the Asp/Glu-Dap interactions with shorter side chains and the Aad-Orn/Lys interactions with longer side chains exhibited stabilizing energetics, mostly relying on electrostatics and hydrophobics, respectively. This suggested the need for length matching of the interacting residues to stabilize the β-hairpin motif. A survey of a nonredundant protein structure database revealed that the statistical sheet pair propensity followed the trend Asp-Lys < Glu-Lys, also implying the need for length matching of the oppositely charged residues.

Effect of charged amino acid side chain length at non-hydrogen bonded strand positions on β-hairpin stability

β-Sheets have been implicated in various neurological disorders, and ∼20% of protein residues adopt a sheet conformation. Therefore, studies on the structural origin of sheet stability can provide fundamental knowledge with potential biomedical applications. Oppositely charged amino acids are frequently observed across one another in antiparallel β-sheets. Interestingly, the side chains of natural charged amino acids Asp, Glu, Arg, Lys have different numbers of hydrophobic methylenes linking the backbone to the hydrophilic charged functionalities. To explore the inherent effect of charged amino acid side chain length on antiparallel sheets, the stability of a designed hairpin motif containing charged amino acids with varying side chain lengths at non-hydrogen bonded positions was studied. Peptides with the guest position on the N-terminal strand and the C-terminal strand were investigated by NMR methods. The charged amino acids (Xaa) included negatively charged residues with a carboxylate group (Asp, Glu, Aad in increasing length), positively charged residues with an ammonium group (Dap, Dab, Orn, Lys in increasing length), and positively charged residues with a guanidinium group (Agp, Agb, Arg, Agh in increasing length). The fraction folded and folding free energy for each peptide were derived from the chemical shift deviation data. The stability of the peptides with the charged residues at the N-terminal guest position followed the trends: Asp > Glu > Aad, Dap < Dab < Orn ∼ Lys, and Agb < Arg < Agh < Agp. The stability of the peptides with the charged residues at the C-terminal guest position followed the trends: Asp < Glu < Aad, Dap ∼ Dab < Orn ∼ Lys, and Agb < Arg ∼ Agp < Agh. These trends were rationalized by thermodynamic sheet propensity and cross-strand interactions.

The supramolecular chemistry of β-sheets

Interactions among β-sheets occur widely in protein quaternary structure, protein-protein interaction, and protein aggregation and are central in Alzheimer's and other amyloid-related diseases. This Perspective looks at the structural biology of these important yet under-appreciated interactions from a supramolecular chemist's point of view. Common themes in the supramolecular interactions of β-sheets are identified and richly illustrated though examples from proteins, amyloids, and chemical model systems. β-Sheets interact through edge-to-edge hydrogen bonding to form extended layers and through face-to-face hydrophobic or van der Waals interactions to form layered sandwich-like structures. Side chains from adjacent layers can fit together through simple hydrophobic contacts or can participate in complementary interdigitation or knob-hole interactions. The layers can be aligned, offset, or rotated. The right-handed twist of β-sheets provides additional opportunities for stabilization of edge-to-edge contacts and rotated layered structures.

Predicting protein β-sheet contacts using a maximum entropy-based correlated mutation measure

MOTIVATION

The problem of ab initio protein folding is one of the most difficult in modern computational biology. The prediction of residue contacts within a protein provides a more tractable immediate step. Recently introduced maximum entropy-based correlated mutation measures (CMMs), such as direct information, have been successful in predicting residue contacts. However, most correlated mutation studies focus on proteins that have large good-quality multiple sequence alignments (MSA) because the power of correlated mutation analysis falls as the size of the MSA decreases. However, even with small autogenerated MSAs, maximum entropy-based CMMs contain information. To make use of this information, in this article, we focus not on general residue contacts but contacts between residues in β-sheets. The strong constraints and prior knowledge associated with β-contacts are ideally suited for prediction using a method that incorporates an often noisy CMM.

RESULTS

Using contrastive divergence, a statistical machine learning technique, we have calculated a maximum entropy-based CMM. We have integrated this measure with a new probabilistic model for β-contact prediction, which is used to predict both residue- and strand-level contacts. Using our model on a standard non-redundant dataset, we significantly outperform a 2D recurrent neural network architecture, achieving a 5% improvement in true positives at the 5% false-positive rate at the residue level. At the strand level, our approach is competitive with the state-of-the-art single methods achieving precision of 61.0% and recall of 55.4%, while not requiring residue solvent accessibility as an input.

AVAILABILITY

http://www2.warwick.ac.uk/fac/sci/systemsbiology/research/software/

Cross-strand interactions of fluorinated amino acids in β-hairpin constructs

We describe herein the design, synthesis, and thermodynamic characterization of fluorinated β-hairpin constructs. Introduction of hexafluoroleucine (Hfl) did not perturb β-hairpin formation, as judged by (1)H NMR structures of four peptides determined to <1 Å backbone RMSDs, allowing direct comparison of thermodynamic stabilities of fluorinated peptides to their hydrocarbon counterparts. Judicious fluorination of peptides often results in increased thermal and chemical stability of the resultant folded structures. However, we found that when cross-strand residue partners were varied, the side-chain interaction energies followed the order Leu-Leu > Hfl-Leu > Hfl-Hfl. All peptides were more structured in 90% MeOH than in aqueous buffers. The peptides with Hfl-Leu or Hfl-Hfl cross-strand partners showed increased interaction energies in this solvent compared to those in water, in contrast to the insignificant effect on Leu-Leu. Our results inform the binding and assembly of peptides containing Hfl in the context of β-sheet structures and may be useful in interpreting binding of fluorinated ligands and peptides to biological targets.

Allosteric peptides bind a caspase zymogen and mediate caspase tetramerization

The caspases are a family of cytosolic proteases with essential roles in inflammation and apoptosis. Drug discovery efforts have focused on developing molecules directed against the active sites of caspases, but this approach has proved challenging and has not yielded any approved therapeutics. Here we describe a new strategy for generating inhibitors of caspase-6, a potential therapeutic target in neurodegenerative disorders, by screening against its zymogen form. Using phage display to discover molecules that bind the zymogen, we report the identification of a peptide that specifically impairs the function of caspase-6 in vitro and in neuronal cells. Remarkably, the peptide binds at a tetramerization interface that is uniquely present in zymogen caspase-6, rather than binding into the active site, and acts via a new allosteric mechanism that promotes caspase tetramerization. Our data illustrate that screening against the zymogen holds promise as an approach for targeting caspases in drug discovery.

Contributions of aromatic pairs to the folding and stability of long-lived human γD-crystallin

Human γD-crystallin (HγD-Crys) is a highly stable protein that remains folded in the eye lens for the majority of an individual's lifetime. HγD-Crys exhibits two homologous crystallin domains, each containing two Greek key motifs with eight β-strands. Six aromatic pairs (four Tyr/Tyr, one Tyr/Phe and one Phe/Phe) are present in the β-hairpin sequences of the Greek keys. Ultraviolet damage to the aromatic residues in lens crystallins may contribute to the genesis of cataract. Mutant proteins with these aromatic residues substituted with alanines were constructed and expressed in E. coli. All mutant proteins except F115A and F117A had lower thermal stability than the WT protein. In equilibrium experiments in guanidine hydrochloride (GuHCl), all mutant proteins had lower thermodynamic stability than the WT protein. N-terminal domain (N-td) substitutions shifted the N-td transition to lower GuHCl concentration, but the C-terminal domain (C-td) transition remained unaffected. C-td substitutions led to a more cooperative unfolding/refolding process, with both the N-td and C-td transitions shifted to lower GuHCl concentration. The aromatic pairs conserved for each Greek key motif (Greek key pairs) had larger contributions to both thermal stability and thermodynamic stability than the other pairs. Aromatic-aromatic interaction was estimated as 1.5-2.0 kcal/mol. In kinetic experiments, N-td substitutions accelerated the early phase of unfolding, while C-td substitutions accelerated the late phase, suggesting independent domain unfolding. Only substitutions of the second Greek key pair of each crystallin domain slowed refolding. The second Greek keys may provide nucleation sites during the folding of the double-Greek-key crystallin domains.

Giant macrolactams based on β-sheet peptides

This paper reports the use of natural amino acids, the tripeptide β-strand mimic Hao, and the β-turn mimic δ-linked ornithine to generate water-soluble 54-, 78-, and 102-membered-ring macrolactams. These giant macrocycles were efficiently prepared by synthesis of the corresponding protected linear peptides, followed by solution-phase cyclization and deprotection. The protected linear peptide precursors were synthesized on 2-chlorotrityl chloride resin by conventional Fmoc-based solid-phase peptide synthesis. Macrocyclization was typically performed using HCTU and N,N-diisopropylethylamine in DMF at ca. 0.5 mM concentration. The macrocycles were isolated in 13-45% overall yield after HPLC purification and lyophilization. 1D, 2D TOCSY, and 2D ROESY (1)H NMR studies of the 54- and 78-membered-ring macrolactams establish that these compounds fold to form β-sheet structures in aqueous solutions.

Tryptophan residues: scarce in proteins but strong stabilizers of β-hairpin peptides

Tryptophan plays important roles in protein stability and recognition despite its scarcity in proteins. Except as fluorescent groups, they have been used rarely in peptide design. Nevertheless, Trp residues were crucial for the stability of some designed minimal proteins. In 2000, Trp-Trp pairs were shown to contribute more than any other hydrophobic interaction to the stability of β-hairpin peptides. Since then, Trp-Trp pairs have emerged as a paradigm for the design of stable β-hairpins, such as the Trpzip peptides. Here, we analyze the nature of the stabilizing capacity of Trp-Trp pairs by reviewing the β-hairpin peptides containing Trp-Trp pairs described up to now, the spectroscopic features and geometry of the Trp-Trp pairs, and their use as binding sites in β-hairpin peptides. To complete the overview, we briefly go through the other relevant β-hairpin stabilizing Trp-non-Trp interactions and illustrate the use of Trp in the design of short peptides adopting α-helical and mixed α/β motifs. This review is of interest in the field of rational design of proteins, peptides, peptidomimetics, and biomaterials.

De novo designed peptides for biological applications

In recent years our ability to design and assemble peptide-based materials and objects de novo (i.e. from first principles) has improved considerably. This brings us to a point where the resulting assemblies are quite sophisticated and amenable to engineering in new functions. Whilst such systems could be used in a variety of ways, biological applications are of particular interest because of the demand for biocompatible, readily produced systems with potential as drug-delivery agents, components of biosensors and scaffolds for 3D cell and tissue culture. This tutorial review describes the building blocks (or tectons) that are being used in peptide assembly, highlights a range of materials and objects that have been produced, notably hydrogels and virus-like particles, and introduces a number of potential applications for the designs.

Bayesian models and algorithms for protein β-sheet prediction

Prediction of the 3D structure greatly benefits from the information related to secondary structure, solvent accessibility, and nonlocal contacts that stabilize a protein's structure. We address the problem of \beta-sheet prediction defined as the prediction of \beta--strand pairings, interaction types (parallel or antiparallel), and \beta-residue interactions (or contact maps). We introduce a Bayesian approach for proteins with six or less \beta-strands in which we model the conformational features in a probabilistic framework by combining the amino acid pairing potentials with a priori knowledge of \beta-strand organizations. To select the optimum \beta-sheet architecture, we significantly reduce the search space by heuristics that enforce the amino acid pairs with strong interaction potentials. In addition, we find the optimum pairwise alignment between \beta-strands using dynamic programming in which we allow any number of gaps in an alignment to model \beta-bulges more effectively. For proteins with more than six \beta-strands, we first compute \beta-strand pairings using the BetaPro method. Then, we compute gapped alignments of the paired \beta-strands and choose the interaction types and \beta--residue pairings with maximum alignment scores. We performed a 10-fold cross-validation experiment on the BetaSheet916 set and obtained significant improvements in the prediction accuracy.

Control of proton and electron transfer in de novo designed, biomimetic β hairpins

Tyrosine side chains are involved in proton coupled electron transfer reactions (PCET) in many complex proteins, including photosystem II (PSII) and ribonucleotide reductase. For example, PSII contains two redox-active tyrosines, TyrD (Y160D2) and TyrZ (Y161D1), which have different protein environments, midpoint potentials, and roles in catalysis. TyrD has a midpoint potential lower than that of TyrZ, and its protein environment is distinguished by potential π-cation interactions with arginine residues. Designed biomimetic peptides provide a system that can be used to investigate how the protein matrix controls PCET reactions. As a model for the redox-active tyrosines in PSII, we are employing a designed, 18 amino acid β hairpin peptide in which PCET reactions occur between a tyrosine (Tyr5) and a cross-strand histidine (His14). In this peptide, the single tyrosine is hydrogen-bonded to an arginine residue, Arg16, and a second arginine, Arg12, has a π-cation interaction with Tyr5. In this report, the effect of these hydrogen bonding and electrostatic interactions on the PCET reactions is investigated. Electrochemical titrations show that histidine substitutions change the nature of PCET reactions, and optical titrations show that Arg16 substitution changes the pK of Tyr5. Removal of Arg16 or Arg12 increases the midpoint potential for tyrosine oxidation. The effects of Arg12 substitution are consistent with the midpoint potential difference, which is observed for the PSII redox-active tyrosine residues. Our results demonstrate that a π-cation interaction, hydrogen bonding, and PCET reactions alter redox-active tyrosine function. These interactions can contribute equally to the control of midpoint potential and reaction rate.

The external pore loop interacts with S6 and S3-S4 linker in domain 4 to assume an essential role in gating control and anticonvulsant action in the Na(+) channel

Carbamazepine, phenytoin, and lamotrigine are widely prescribed anticonvulsants in neurological clinics. These drugs bind to the same receptor site, probably with the diphenyl motif in their structure, to inhibit the Na⁺ channel. However, the location of the drug receptor remains controversial. In this study, we demonstrate close proximity and potential interaction between an external aromatic residue (W1716 in the external pore loop) and an internal aromatic residue (F1764 in the pore-lining part of the sixth transmembrane segment, S6) of domain 4 (D4), both being closely related to anticonvulsant and/or local anesthetic binding to the Na⁺ channel. Double-mutant cycle analysis reveals significant cooperativity between the two phenyl residues for anticonvulsant binding. Concomitant F1764C mutation evidently decreases the susceptibility of W1716C to external Cd²⁺ and membrane-impermeable methanethiosulfonate reagents. Also, the W1716E/F1764R and G1715E/F1764R double mutations significantly alter the selectivity for Na⁺ over K⁺ and markedly shift the activation curve, respectively. W1716 and F1764 therefore very likely form a link connecting the outer and inner compartments of the Na⁺ channel pore (in addition to the selectivity filter). Anticonvulsants and local anesthetics may well traverse this “S6 recess” without trespassing on the selectivity filter. Furthermore, we found that Y1618K, a point mutation in the S3-4 linker (the extracellular extension of D4S4), significantly alters the consequences of carbamazepine binding to the Na⁺ channel. The effect of Y1618K mutation, however, is abolished by concomitant point mutations in the vicinity of Y1618, but not by those in the internally located inactivation machinery, supporting a direct local rather than a long-range allosteric action. Moreover, Y1618 could interact with D4 pore residues W1716 and L1719 to have a profound effect on both channel gating and anticonvulsant action. We conclude that there are direct interactions among the external S3-4 linker, the external pore loop, and the internal S6 segment in D4, making the external pore loop a pivotal point critically coordinating ion permeation, gating, and anticonvulsant binding in the Na⁺ channel.

Prodepth: predict residue depth by support vector regression approach from protein sequences only

Residue depth (RD) is a solvent exposure measure that complements the information provided by conventional accessible surface area (ASA) and describes to what extent a residue is buried in the protein structure space. Previous studies have established that RD is correlated with several protein properties, such as protein stability, residue conservation and amino acid types. Accurate prediction of RD has many potentially important applications in the field of structural bioinformatics, for example, facilitating the identification of functionally important residues, or residues in the folding nucleus, or enzyme active sites from sequence information. In this work, we introduce an efficient approach that uses support vector regression to quantify the relationship between RD and protein sequence. We systematically investigated eight different sequence encoding schemes including both local and global sequence characteristics and examined their respective prediction performances. For the objective evaluation of our approach, we used 5-fold cross-validation to assess the prediction accuracies and showed that the overall best performance could be achieved with a correlation coefficient (CC) of 0.71 between the observed and predicted RD values and a root mean square error (RMSE) of 1.74, after incorporating the relevant multiple sequence features. The results suggest that residue depth could be reliably predicted solely from protein primary sequences: local sequence environments are the major determinants, while global sequence features could influence the prediction performance marginally. We highlight two examples as a comparison in order to illustrate the applicability of this approach. We also discuss the potential implications of this new structural parameter in the field of protein structure prediction and homology modeling. This method might prove to be a powerful tool for sequence analysis.