Determinants of strand register in antiparallel beta-sheets of proteins

Computational analysis to comprehend the structure-function properties of fibrinolytic enzymes from Bacillus spp for their efficient integration into industrial applications

Background

The fibrinolytic enzymes from Bacillus sp. are proposed as therapeutics in preventing thrombosis. Computational-based analyses of these enzymes' amino acid composition, basic physiological properties, presence of functional domain and motifs, and secondary and tertiary structure analyses can lead to developing a specific enzyme with improved catalytic activity and other properties that may increase their therapeutic potential.

Methods

The nucleotide sequences of fibrinolytic enzymes produced by the genus Bacillus and its corresponding protein sequences were retrieved from the NCBI database and aligned using the PRALINE programme. The varied physiochemical parameters and structural and functional analysis of the enzyme sequences were carried out with the ExPASy-ProtParam tool, MEME server, SOPMA, PDBsum tool, CYS-REC tool, SWISS-MODEL, SAVES servers, TMHMM program, GlobPlot, and peptide cutter software. The assessed in-silico data were compared with the published experimental results for validation.

Results

The alignment of sixty fibrinolytic serine protease enzymes (molecular mass 12-86 kDa) sequences showed 49 enzymes possess a conserved domain with a catalytic triad of Asp196, His242, and Ser569. The predicted instability and aliphatic indexes were 1.94-37.77, and 68.9-93.41, respectively, indicating high thermostability. The random coil means value suggested the predominance of this secondary structure in these proteases. A set of 50 amino acid residues representing motif 3 signifies the Peptidase S8/S53 domain that was invariably observed in 56 sequences. Additionally, 28 sequences have transmembrane helices, with two having the most disordered areas, and they pose 25 enzyme cleavage sites. A comparative analysis of the experimental work with the results of in-silico study put forward the characteristics of the enzyme sequences JF739176.1 and MF677779.1 to be considered when creating a potential mutant enzyme as these sequences are stable at high pH with thermostability and to exhibit αβ-fibrinogenase activity in both experimental and in-silico studies.

β-Strand-mediated Domain-swapping in the Absence of Hydrophobic Core Repacking

Domain swapping is a process wherein a portion of a protein is exchanged with its counterpart in another copy of the molecule, resulting in the formation of homo-oligomers with concomitant repacking of a hydrophobic core. Here, we report domain swapping triggered upon modifying a β-hairpin sequence within a single-layer β-sheet (SLB) of a model protein, OspA that did not involve the formation of a reorganized hydrophobic core. The replacement of two β-hairpin sequences with a Gly-Gly and shorteing of a β-hairpin resulted in a protein that formed two distinct crystal structures under similar conditions: one was monomeric, similar to the parental molecule, whereas the other was a domain-swapped dimer, mediated by an intermolecular β-sheet in the SLB portion. Based on the dimer interface structure, we replaced the Gly-Gly sequence with three-residue sequences that enable the formation of a consecutive intermolecular β-sheet, including the Cys-Thr-Cys sequence that formed a stable disulfide-linked dimer. These results provide new insights into protein folding, evolution, and the designability of protein structure.

Macrocyclic β-Sheets Stabilized by Hydrogen Bond Surrogates

Mimics of protein secondary and tertiary structure offer rationally-designed inhibitors of biomolecular interactions. β-Sheet mimics have a storied history in bioorganic chemistry and are typically designed with synthetic or natural turn segments. We hypothesized that replacement of terminal inter-β-strand hydrogen bonds with hydrogen bond surrogates (HBS) may lead to conformationally-defined macrocyclic β-sheets without the requirement for natural or synthetic β-turns, thereby providing a minimal mimic of a protein β-sheet. To access turn-less antiparallel β-sheet mimics, we developed a facile solid phase synthesis protocol. We surveyed a dataset of protein β-sheets for naturally observed interstrand side chain interactions. This bioinformatics survey highlighted an over-abundance of aromatic-aromatic, cation-π and ionic interactions in β-sheets. In correspondence with natural β-sheets, we find that minimal HBS mimics show robust β-sheet formation when specific amino acid residue pairings are incorporated. In isolated β-sheets, aromatic interactions endow superior conformational stability over ionic or cation-π interactions. Circular dichroism and NMR spectroscopies, along with high-resolution X-ray crystallography, support our design principles.

AssignSLP_GUI, a software tool exploiting AI for NMR resonance assignment of sparsely labeled proteins

Not all proteins are amenable to uniform isotopic labeling with 13C and 15N, something needed for the widely used, and largely deductive, triple resonance assignment process. Among them are proteins expressed in mammalian cell culture where native glycosylation can be maintained, and proper formation of disulfide bonds facilitated. Uniform labeling in mammalian cells is prohibitively expensive, but sparse labeling with one or a few isotopically enriched amino acid types is an option for these proteins. However, assignment then relies on accessing the best match between a variety of measured NMR parameters and predictions based on 3D structure, often from X-ray crystallography. Finding this match is a challenging process that has benefitted from many computational tools, including trained neural nets for chemical shift prediction, genetic algorithms for searches through a myriad of assignment possibilities, and now AI-based prediction of high-quality structures for protein targets. AssignSLP_GUI, a new version of a software package for assignment of resonances from sparsely-labeled proteins, uses many of these tools. These tools and new additions to the package are highlighted in an application to a sparsely-labeled domain from a glycoprotein, CEACAM1.

NMR analysis suggests the terminal domains of Robo1 remain extended but are rigidified in the presence of heparan sulfate

Human roundabout 1 (hRobo1) is an extracellular receptor glycoprotein that plays important roles in angiogenesis, organ development, and tumor progression. Interaction between hRobo1 and heparan sulfate (HS) has been shown to be essential for its biological activity. To better understand the effect of HS binding we engineered a lanthanide-binding peptide sequence (Loop) into the Ig2 domain of hRobo1. Native mass spectrometry was used to verify that loop introduction did not inhibit HS binding or conformational changes previously suggested by gas phase ion mobility measurements. NMR experiments measuring long-range pseudocontact shifts were then performed on 13C-methyl labeled hRobo1-Ig1-2-Loop in HS-bound and unbound forms. The magnitude of most PCSs for methyl groups in the Ig1 domain increase in the bound state confirming a change in the distribution of interdomain geometries. A grid search over Ig1 orientations to optimize the fit of data to a single conformer for both forms produced two similar structures, both of which differ from existing X-ray crystal structures and structures inferred from gas-phase ion mobility measurements. The structures and degree of fit suggest that the hRobo1-Ig1-2 structure changes slightly and becomes more rigid on HS binding. This may have implications for Robo-Slit signaling.

Modeling protein structure as a stable static equilibrium

We present evidence that the protein structure can be modeled as a stable static equilibrium, determined mainly by compressive supports in the nonpolar interior. That is, protein structures derive their structural strength through the same mechanical principles as do conventional structures like bridges and buildings. This is based on the observation that the experimentally elucidated structural determinants, the interior nonpolar side chains, are engaged in strong compressions in static terms. At the same time, major substructures in proteins, helices and h-bonded strands, because of their geometry, inherently leave gaps in the space they occupy. Under the compressive force, nonpolar side chains from one substructure can protrude into the gaps of another neighboring substructure and block its motion. As a result, interlocking of substructures can form, which builds up the nonpolar core assembly. The native structure then is the one with the structurally most stable core assembly. While intuitively appealing, this is a radical departure from the prevailing thinking that protein native structure is determined by global energy minimum, which is founded on thermodynamic hypothesis. Furthermore, to develop an effective model for analyzing protein structure with conventional tools, a proper mechanical representation must be established. By proving that the stability of the equilibrium in compressive interactions is conditioned on a form of mechanical energy minimum, we show that our notion of native structure can be equally consistent with the thermodynamic hypothesis. By mathematically treating the blocking action, an interaction, as a bar, a physical object, we succeed in representing and analyzing the core assembly as truss, a conventional structure. In this paper we define and expound step-by-step increasingly integrated interlocking patterns. We then analyze the core assemblies of a large set of diverse protein database structures. A native structure can be distinguished from decoys by comparing the composition and strength of their core assemblies. We show the results for two sets of native structures vs decoys.

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.

Computational Identification and Design of Complementary β-Strand Sequences

The ß-sheet is a regular secondary structure element which consists of linear segments called ß-strands. They are involved in many important biological processes, and some are known to be related to serious diseases such as neurologic disorders and amyloidosis. The self-assembly of ß-sheet peptides also has practical applications in material sciences since they can be building blocks of repeated nanostructures. Therefore, computational algorithms for identification of ß-sheet formation can offer useful insight into the mechanism of disease-prone protein segments and the construction of biocompatible nanomaterials. Despite the recent advances in structure-based methods for the assessment of atomic interactions, identifying amyloidogenic peptides has proven to be extremely difficult since they are structurally very flexible. Thus, an alternative strategy is required to describe ß-sheet formation. It has been hypothesized and observed that there are certain amino acid propensities between ß-strand pairs. Based on this hypothesis, a database search algorithm, B-SIDER, is developed for the identification and design of ß-sheet forming sequences. Given a target sequence, the algorithm identifies exact or partial matches from the structure database and constructs a position-specific score matrix. The score matrix can be utilized to design novel sequences that can form a ß-sheet specifically with the target.

Folding in Place: Design of β-Strap Motifs to Stabilize the Folding of Hairpins with Long Loops

Despite their pivotal role in defining antibody affinity and protein function, β-hairpins harboring long noncanonical loops remain synthetically challenging because of the large entropic penalty associated with their conformational folding. Little is known about the contribution and impact of stabilizing motifs on the folding of β-hairpins with loops of variable length and plasticity. Here, we report a design of minimalist β-straps (strap = strand + cap) that offset the entropic cost of long-loop folding. The judicious positioning of noncovalent interactions (hydrophobic cluster and salt-bridge) within the novel 8-mer β-strap design RW(V/H)W···WVWE stabilizes hairpins with up to 10-residue loops of varying degrees of plasticity (Tm up to 52 °C; 88 ± 1% folded at 18 °C). This "hyper" thermostable β-strap outperforms the previous gold-standard technology of β-strand-β-cap (16-mer) and provides a foundation for producing new classes of long hairpins as a viable and practical alternative to macrocyclic peptides.

Inhibitor and substrate cooperate to inhibit amyloid fibril elongation of α-synuclein

In amyloid fibril elongation, soluble growth substrate binds to the fibril-end and converts into the fibril conformation. This process is targeted by inhibitors that block fibril-ends. Here, we investigated how the elongation of α-synuclein (αS) fibrils, which are associated with Parkinson's disease and other synucleinopathies, is inhibited by αS variants with a preformed hairpin in the critical N-terminal region comprising residues 36-57. The inhibitory efficiency is strongly dependent on the specific position of the hairpin. We find that the inhibitor and substrate concentration dependencies can be analyzed with models of competitive enzyme inhibition. Remarkably, the growth substrate, i.e., wild-type αS, supports inhibition by stabilizing the elongation-incompetent blocked state. This observation allowed us to create inhibitor-substrate fusions that achieved inhibition at low nanomolar concentration. We conclude that inhibitor-substrate cooperativity can be exploited for the design of fibril growth inhibitors.

Propensities of Amino Acid Pairings in Secondary Structure of Globular Proteins

A class of secondary structure prediction algorithms use the information from the statistics of the residue pairs found in secondary structural elements. Because the protein folding process is dominated by backbone hydrogen bonding, an approach based on backbone hydrogen-bonded residue pairings would improve the predicting capabilities of these class algorithms. The reliability of the prediction algorithms depends on the quality of the statistics, therefore, of the data set. In this study, it was aimed to determine the propensities of the backbone hydrogen-bonded residue pairings for secondary structural elements of α-helix and β-sheet in globular proteins using a new and comprehensive data set created from the peptides deposited in Worldwide Protein Data Bank. A master data set including 4882 globular peptide chains with resolution better than 2.5 Å, sequence identity smaller than 25% and length of no shorter than 100 residues were created. Separate data sub sets also were created for helix and sheet structures from master set and each sub set includes 4594 and 4483 chains, respectively. Backbone hydrogen-bonded residue pairings in helices and sheets were detected and the propensities of them were represented as odds ratios (observed/[random or expected]) in matrices. Propensities assigned by this study to the residue pairings in secondary structural elements (as helix, overall strands, parallel strands and antiparallel strands) differ from the previous studies by 19 to 34%. These dissimilarities are important and they would cause further improvements in secondary structure prediction algorithms.

ROSETTA-informed design of structurally stabilized cyclic anti-amyloid peptides

β-amyloid oligomers are thought to be the most toxic species formed en route to fibril deposition in Alzheimer's disease. Transthyretin is a natural sequestering agent of β-amyloid oligomers: the binding site to β-amyloid has been traced to strands G/H of the inner β-sheet of transthyretin. A linear peptide, with the same primary sequence as the β-amyloid binding domain on transthyretin, was moderately effective at inhibiting β-amyloid fibril growth. Insertion of a β-turn template and cyclization greatly increased stability against proteolysis and improved efficacy as an amyloid inhibitor. However, the cyclic peptide still contained a significant amount of disorder. Using the Simple Cyclic Peptide Application within ROSETTA as an in silico predictor of cyclic peptide conformation and stability, we investigated putative structural enhancements, including stabilization by disulfide linkages and insertion of a second β-turn template. Several candidates were synthesized and tested for secondary structure and ability to inhibit β-amyloid aggregation. The results demonstrate that cyclization, β-sheet structure and conformational homogeneity are all preferable design features, whereas disulfide bond formation across the two β-strands is not preferable.

B-SIDER: Computational Algorithm for the Design of Complementary β-Sheet Sequences

The β-sheet is an element of protein secondary structure, and intra-/intermolecular β-sheet interactions play pivotal roles in biological regulatory processes including scaffolding, transporting, and oligomerization. In nature, a β-sheet formation is tightly regulated because dysregulated β-stacking often leads to severe diseases such as Alzheimer's, Parkinson's, systemic amyloidosis, or diabetes. Thus, the identification of intrinsic β-sheet-forming propensities can provide valuable insight into protein designs for the development of novel therapeutics. However, structure-based design methods may not be generally applicable to such amyloidogenic peptides mainly owing to high structural plasticity and complexity. Therefore, an alternative design strategy based on complementary sequence information is of significant importance. Herein, we developed a database search method called β-Stacking Interaction DEsign for Reciprocity (B-SIDER) for the design of complementary β-strands. This method makes use of the structural database information and generates target-specific score matrices. The discriminatory power of the B-SIDER score function was tested on representative amyloidogenic peptide substructures against a sequence-based score matrix (PASTA 2.0) and two popular ab initio protein design score functions (Rosetta and FoldX). B-SIDER is able to distinguish wild-type amyloidogenic β-strands as favored interactions in a more consistent manner than other methods. B-SIDER was prospectively applied to the design of complementary β-strands for a splitGFP scaffold. Three variants were identified to have stronger interactions than the original sequence selected through a directed evolution, emitting higher fluorescence intensities. Our results indicate that B-SIDER can be applicable to the design of other β-strands, assisting in the development of therapeutics against disease-related amyloidogenic peptides.

The Membrane-Integrated Steric Chaperone Lif Facilitates Active Site Opening of Pseudomonas aeruginosa Lipase A

Lipases are essential and widely used biocatalysts. Hence, the production of lipases requires a detailed understanding of the molecular mechanism of its folding and secretion. Lipase A from Pseudomonas aeruginosa, PaLipA, constitutes a prominent example that has additional relevance because of its role as a virulence factor in many diseases. PaLipA requires the assistance of a membrane-integrated steric chaperone, the lipase-specific foldase Lif, to achieve its enzymatically active state. However, the molecular mechanism of how Lif activates its cognate lipase has remained elusive. Here, we show by molecular dynamics simulations at the atomistic level and potential of mean force computations that Lif catalyzes the activation process of PaLipA by structurally stabilizing an intermediate PaLipA conformation, particularly a β-sheet in the region of residues 17-30, such that the opening of PaLipA's lid domain is facilitated. This opening allows substrate access to PaLipA's catalytic site. A surprising and so far not fully understood aspect of our study is that the open state of PaLipA is unstable compared to the closed one according to our computational and in vitro biochemical results. We thus speculate that further interactions of PaLipA with the Xcp secretion machinery and/or components of the extracellular matrix contribute to the remaining activity of secreted PaLipA. © 2019 Wiley Periodicals, Inc.

Cross-strand disulfides in the hydrogen bonding site of antiparallel β-sheet (aCSDhs): Forbidden disulfides that are highly strained, easily broken

Some disulfide bonds perform important structural roles in proteins, but another group has functional roles via redox reactions. Forbidden disulfides are stressed disulfides found in recognizable protein contexts, which currently constitute more than 10% of all disulfides in the PDB. They likely have functional redox roles and constitute a major subset of all redox-active disulfides. The torsional strain of forbidden disulfides is typically higher than for structural disulfides, but not so high as to render them immediately susceptible to reduction under physionormal conditions. Previously we characterized the most abundant forbidden disulfide in the Protein Data Bank, the aCSDn: a canonical motif in which disulfide-bonded cysteine residues are positioned directly opposite each other on adjacent anti-parallel β-strands such that the backbone hydrogen-bonded moieties are directed away from each other. Here we perform a similar analysis for the aCSDh, a less common motif in which the opposed cysteine residues are backbone hydrogen bonded. Oxidation of two Cys in this context places significant strain on the protein system, with the β-chains tilting toward each other to allow disulfide formation. Only left-handed aCSDh conformations are compatible with the inherent right-handed twist of β-sheets. aCSDhs tend to be more highly strained than aCSDns, particularly when both hydrogen bonds are formed. We discuss characterized roles of aCSDh motifs in proteins of the dataset, which include catalytic disulfides in ribonucleotide reductase and ahpC peroxidase as well as a redox-active disulfide in P1 lysozyme, involved in a major conformation change. The dataset also includes many binding proteins.

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.

Ribosomal RACK1:Protein Kinase C βII Modulates Intramolecular Interactions between Unstructured Regions of Eukaryotic Initiation Factor 4G (eIF4G) That Control eIF4E and eIF3 Binding

The receptor for activated C kinase (RACK1), a conserved constituent of eukaryotic ribosomes, mediates phosphorylation of eukaryotic initiation factor 4G1(S1093) [eIF4G1(S1093)] and eIF3a(S1364) by protein kinase C βII (PKCβII) (M. I. Dobrikov, E. Y. Dobrikova, and M. Gromeier, Mol Cell Biol 38:e00304-18, 2018, https://doi.org/10.1128/MCB.00304-18). RACK1:PKCβII activation drives a phorbol ester-induced surge of global protein synthesis and template-specific translation induction of PKC-Raf-extracellular signal-regulated kinase 1/2 (ERK1/2)-responsive genes. For unraveling mechanisms of RACK1:PKCβII-mediated translation stimulation, we used sequentially truncated eIF4G1 in coimmunoprecipitation analyses to delineate a set of autoinhibitory elements in the N-terminal unstructured region (surrounding the eIF4E-binding motif) and the interdomain linker (within the eIF3-binding site) of eIF4G1. Computer-based predictions of secondary structure, mutational analyses, and fluorescent titration with the β-sheet dye thioflavin T suggest that eIF4G1(S1093) modulates a 4-stranded β-sheet composed of antiparallel β-hairpins formed by the autoinhibitory elements in eIF4G1's unstructured regions. The intact β-sheet "locks" the eIF4G configuration, preventing assembly with eIF3/40S ribosomal subunits. Upon PKC stimulation, activated RACK1:PKCβII phosphorylates eIF4G(S1093) in the tight 48S initiation complex, possibly facilitating dissociation/recycling of eIF4F.

JP. A multiphase transitioning peptide hydrogel for suturing ultrasmall vessels

Many surgeries are complicated by the need to anastomose, or reconnect, micrometre-scale vessels. Although suturing remains the gold standard for anastomosing vessels, it is difficult to place sutures correctly through collapsed lumen, making the procedure prone to failure. Here, we report a multiphase transitioning peptide hydrogel that can be injected into the lumen of vessels to facilitate suturing. The peptide, which contains a photocaged glutamic acid, forms a solid-like gel in a syringe and can be shear-thin delivered to the lumen of collapsed vessels (where it distends the vessel) and the space between two vessels (where it is used to approximate the vessel ends). Suturing is performed directly through the gel. Light is used to initiate the final gel-sol phase transition that disrupts the hydrogel network, allowing the gel to be removed and blood flow to resume. This gel adds a new tool to the armamentarium for micro- and supermicrosurgical procedures.

Helix-helix interactions: is the medium the message?

In this issue of Structure, Zhang and colleagues compare the helix-helix interaction spaces of an extensive database of soluble and membrane proteins. Intriguingly, the resultant clusters show similar helix interaction geometries between the protein classes, differing in detail only by patterns of local interactions and inter-helical distances.

Separation of peptide fragments of a protein kinase C substrate fused to a β-hairpin by capillary electrophoresis

Synthetic peptides incorporating well-folded β-hairpin peptides possess advantages in a variety of cell biology applications by virtue of increased resistance to proteolytic degradation. In this study, the WKpG β-hairpin peptide fused to a protein kinase C (PKC) substrate was synthesized, and capillary-electrophoretic separation conditions for this peptide and its proteolytic fragments were developed. Fragments of WKpG-PKC were generated by enzymatic treatment with trypsin and Pronase E to produce standards for identification of degradation fragments in a cellular lysate. A simple buffer system of 250 mM H3PO4, pH 1.5 enabled separation of WKpG-PKC and its fragments by capillary electrophoresis in less than 16 min. Using a cellular lysate produced from Ba/F3 cells, the β-hairpin-conjugated substrate and its PKCα-phosphorylated product could be detected and separated from peptidase-generated fragments produced in a cell lysate. The method has potential application for identification and quantification of WKpG-PKC and its fragments in complex biological systems when the peptide is used as a reporter to assay PKC activity.

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.

Cross-strand disulfides in the non-hydrogen bonding site of antiparallel β-sheet (aCSDns): poised for biological switching

aCSDns are forbidden disulfides with protein redox-activity. Within the aCSDn structural motif, a cognate substrate of Trx-like enzymes, the disulfide bonds are strained and metastable, facilitating their role as redox-regulated protein switches.

Neighbor preferences of amino acids and context-dependent effects of amino acid substitutions in human, mouse, and dog

Amino acids show apparent propensities toward their neighbors. In addition to preferences of amino acids for their neighborhood context, amino acid substitutions are also considered to be context-dependent. However, context-dependence patterns of amino acid substitutions still remain poorly understood. Using relative entropy, we investigated the neighbor preferences of 20 amino acids and the context-dependent effects of amino acid substitutions with protein sequences in human, mouse, and dog. For 20 amino acids, the highest relative entropy was mostly observed at the nearest adjacent site of either N- or C-terminus except C and G. C showed the highest relative entropy at the third flanking site and periodic pattern was detected at G flanking sites. Furthermore, neighbor preference patterns of amino acids varied greatly in different secondary structures. We then comprehensively investigated the context-dependent effects of amino acid substitutions. Our results showed that nearly half of 380 substitution types were evidently context dependent, and the context-dependent patterns relied on protein secondary structures. Among 20 amino acids, P elicited the greatest effect on amino acid substitutions. The underlying mechanisms of context-dependent effects of amino acid substitutions were possibly mutation bias at a DNA level and natural selection. Our findings may improve secondary structure prediction algorithms and protein design; moreover, this study provided useful information to develop empirical models of protein evolution that consider dependence between residues.

NXO beta structure mimicry: an ultrashort turn/hairpin mimic that folds in water

An NXO building block derived tetrapeptide mimic emulates a natural proline-glycine β-turn/hairpin in polar media, including water at room temperature.

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.

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.

β-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.

Statistical Analysis of Terminal Extensions of Protein β-Strand Pairs

The long-range interactions, required to the accurate predictions of tertiary structures of β-sheet-containing proteins, are still difficult to simulate. To remedy this problem and to facilitate β-sheet structure predictions, many efforts have been made by computational methods. However, known efforts on β-sheets mainly focus on interresidue contacts or amino acid partners. In this study, to go one step further, we studied β-sheets on the strand level, in which a statistical analysis was made on the terminal extensions of paired β-strands. In most cases, the two paired β-strands have different lengths, and terminal extensions exist. The terminal extensions are the extended part of the paired strands besides the common paired part. However, we found that the best pairing required a terminal alignment, and β-strands tend to pair to make bigger common parts. As a result, 96.97%  of β-strand pairs have a ratio of 25% of the paired common part to the whole length. Also 94.26% and 95.98%  of β-strand pairs have a ratio of 40% of the paired common part to the length of the two β-strands, respectively. Interstrand register predictions by searching interacting β-strands from several alternative offsets should comply with this rule to reduce the computational searching space to improve the performances of algorithms.

A hierarchical order within protein structures underlies large separations between strands in β-sheets

Protein β-sheets often involve nonlocal interactions between parts of the polypeptide chain that are separated by hundreds of residues, raising the question of how these nonlocal contacts form. A recent study of the smallest β-sheets found that their formation was not driven by signals hidden in the primary sequence. Instead, the strands in these sheets were either local in sequence, or, when separated by large sequential distances, the intervening residues were found to fold into compact modules that anchored distant parts of the chain in close spatial proximity. Here, we examine larger β-sheets to investigate the extensibility of this principle. From an analysis of the β-sheets in a nonredundant protein dataset, we find that a highly ordered hierarchical relationship exists in the intervening structure between nonlocal β-strands. This observation is almost universal: virtually all β-sheets, no matter their complexity, appear to adopt an antiparallel model to manage the nonlocal aspects of their assembly, one where the chain, having left the vicinity of an unfinished β-sheet, retraces its steps via the same route to complete the initial sheet. Exceptions typically involve unstructured regions at chain termini. Moreover, an analysis of the residues involved in nonlocal crossstrand interactions did not produce any evidence of a signal hidden in the sequence that might direct long-range interactions. These results build on those reported for the smallest sheets, suggesting that sheet formation is either local in sequence or local in space following prior folding events that anchor disparate parts of the chain in close proximity.

Chemical strategies for controlling protein folding and elucidating the molecular mechanisms of amyloid formation and toxicity

It has been more than a century since the first evidence linking the process of amyloid formation to the pathogenesis of Alzheimer's disease. During the last three decades in particular, increasing evidence from various sources (pathology, genetics, cell culture studies, biochemistry, and biophysics) continues to point to a central role for the pathogenesis of several incurable neurodegenerative and systemic diseases. This is in part driven by our improved understanding of the molecular mechanisms of protein misfolding and aggregation and the structural properties of the different aggregates in the amyloid pathway and the emergence of new tools and experimental approaches that permit better characterization of amyloid formation in vivo. Despite these advances, detailed mechanistic understanding of protein aggregation and amyloid formation in vitro and in vivo presents several challenges that remain to be addressed and several fundamental questions about the molecular and structural determinants of amyloid formation and toxicity and the mechanisms of amyloid-induced toxicity remain unanswered. To address this knowledge gap and technical challenges, there is a critical need for developing novel tools and experimental approaches that will not only permit the detection and monitoring of molecular events that underlie this process but also allow for the manipulation of these events in a spatial and temporal fashion both in and out of the cell. This review is primarily dedicated in highlighting recent results that illustrate how advances in chemistry and chemical biology have been and can be used to address some of the questions and technical challenges mentioned above. We believe that combining recent advances in the development of new fluorescent probes, imaging tools that enabled the visualization and tracking of molecular events with advances in organic synthesis, and novel approaches for protein synthesis and engineering provide unique opportunities to gain a molecular-level understanding of the process of amyloid formation. We hope that this review will stimulate further research in this area and catalyze increased collaboration at the interface of chemistry and biology to decipher the mechanisms and roles of protein folding, misfolding, and aggregation in health and disease.

Parallel β-sheet secondary structure is stabilized and terminated by interstrand disulfide cross-linking

Disulfide bonds between Cys residues in adjacent strands of parallel β-sheets are rare among proteins, which suggests that parallel β-sheet structure is not stabilized by such disulfide cross-links. We report experimental results that show, surprisingly, that an interstrand disulfide bond can stabilize parallel β-sheets formed by an autonomously folding peptide in aqueous solution. NMR analysis reveals that parallel β-sheet structure is terminated beyond the disulfide bond, which causes deviation from the extended backbone conformation at one of the Cys residues.

Structural and functional analysis of the tandem β-zipper interaction of a Streptococcal protein with human fibronectin

Bacterial fibronectin-binding proteins (FnBPs) contain a large intrinsically disordered region (IDR) that mediates adhesion of bacteria to host tissues, and invasion of host cells, through binding to fibronectin (Fn). These FnBP IDRs consist of Fn-binding repeats (FnBRs) that form a highly extended tandem β-zipper interaction on binding to the N-terminal domain of Fn. Several FnBR residues are highly conserved across bacterial species, and here we investigate their contribution to the interaction. Mutation of these residues to alanine in SfbI-5 (a disordered FnBR from the human pathogen Streptococcus pyogenes) reduced binding, but for each residue the change in free energy of binding was <2 kcal/mol. The structure of an SfbI-5 peptide in complex with the second and third F1 modules from Fn confirms that the conserved FnBR residues play equivalent functional roles across bacterial species. Thus, in SfbI-5, the binding energy for the tandem β-zipper interaction with Fn is distributed across the interface rather than concentrated in a small number of "hot spot" residues that are frequently observed in the interactions of folded proteins. We propose that this might be a common feature of the interactions of IDRs and is likely to pose a challenge for the development of small molecule inhibitors of FnBP-mediated adhesion to and invasion of host cells.

Structural role of the terminal disulfide bond in the sweetness of brazzein

Brazzein, a 54 residue sweet-tasting protein, is thought to participate in a multipoint binding interaction with the sweet taste receptor. Proposed sites for interaction with the receptor include 2 surface loops and the disulfide bond that connects the N- and C-termini. However, the importance of each site is not well understood. To characterize the structural role of the termini in the sweetness of brazzein, the position of the disulfide bond connecting the N- and C-termini was shifted by substituting K3-C4-K5 with C3-K4-R5. The apparent affinity and V(max) of the C3-K4-R5-brazzein (CKR-brazzein) variant were only modestly decreased compared with the wild-type (WT) brazzein. We determined a high-resolution structure of CKR-brazzein by nuclear magnetic resonance spectroscopy (backbone root mean square deviation of 0.39 Å). Comparing the structure of CKR-brazzein with that of WT-brazzein revealed that the terminal β-strands of the variant display extended β-structure and increased dynamics relative to WT-brazzein. These results support previous mutagenesis studies and further suggest that, whereas interactions involving the termini are necessary for full function of brazzein, the termini do not constitute the primary site of interaction between brazzein and the sweet taste receptor.

Computational studies of membrane proteins: models and predictions for biological understanding

We discuss recent progresses in computational studies of membrane proteins based on physical models with parameters derived from bioinformatics analysis. We describe computational identification of membrane proteins and prediction of their topology from sequence, discovery of sequence and spatial motifs, and implications of these discoveries. The detection of evolutionary signal for understanding the substitution pattern of residues in the TM segments and for sequence alignment is also discussed. We further discuss empirical potential functions for energetics of inserting residues in the TM domain, for interactions between TM helices or strands, and their applications in predicting lipid-facing surfaces of the TM domain. Recent progresses in structure predictions of membrane proteins are also reviewed, with further discussions on calculation of ensemble properties such as melting temperature based on simplified state space model. Additional topics include prediction of oligomerization state of membrane proteins, identification of the interfaces for protein-protein interactions, and design of membrane proteins. This article is part of a Special Issue entitled: Protein Folding in Membranes.

Hydrophobic interactions in the formation of secondary structures in small peptides

Effects of the attractive and repulsive parts of hydrophobic interactions on α helices and β sheets in small peptides are investigated using a simple atomic potential. Typically, a physical spatial range of attraction tends to favor β sheets, but α helices would be favored if the attractive range were more extended. We also found that desolvation barriers favor β sheets in collapsed conformations of polyalanine, polyvaline, polyleucine, and three fragments of amyloid peptides tested in this study. Our results provide insight into the multifaceted role of hydrophobicity in secondary structure formation, including the α to β transitions in certain amyloid peptides.

Studies on the rules of β-strand alignment in a protein β-sheet structure

To further disclose the underlying mechanisms of protein β-sheet formation, studies were made on the rules of β-strands alignment forming β-sheet structure using statistical and machine learning approaches. Firstly, statistical analysis was performed on the sum of β-strands between each β-strand pairs in protein sequences. The results showed a propensity of near-neighbor pairing (or called "first come first pair") in the β-strand pairs. Secondly, based on the same dataset, the pairwise cross-combinations of real β-strand pairs and four pseudo-β-strand contained pairs were classified by support vector machine (SVM). A novel feature extracting approach was designed for classification using the average amino acid pairing encoding matrix (APEM). Analytical results of the classification indicated that a segment of β-strand had the ability to distinguish β-strands from segments of α-helix and coil. However, the result also showed that a β-strand was not strongly conserved to choose its real partner from all the alternative β-strand partners, which was corresponding with the ordination results of the statistical analysis each other. Thus, the rules of "first come first pair" propensity and the non-conservative ability to choose real partner, were possible important factors affecting the β-strands alignment forming β-sheet structures.