Stability of cyclic beta-hairpins: asymmetric contributions from side chains of a hydrogen-bonded cross-strand residue pair

4-(Azolyl)-Benzamidines as a Novel Chemotype for ASIC1a Inhibitors

Acid-sensing ion channels (ASICs) play a key role in the perception and response to extracellular acidification changes. These proton-gated cation channels are critical for neuronal functions, like learning and memory, fear, mechanosensation and internal adjustments like synaptic plasticity. Moreover, they play a key role in neuronal degeneration, ischemic neuronal injury, seizure termination, pain-sensing, etc. Functional ASICs are homo or heterotrimers formed with (ASIC1-ASIC3) homologous subunits. ASIC1a, a major ASIC isoform in the central nervous system (CNS), possesses an acidic pocket in the extracellular region, which is a key regulator of channel gating. Growing data suggest that ASIC1a channels are a potential therapeutic target for treating a variety of neurological disorders, including stroke, epilepsy and pain. Many studies were aimed at identifying allosteric modulators of ASIC channels. However, the regulation of ASICs remains poorly understood. Using all available crystal structures, which correspond to different functional states of ASIC1, and a molecular dynamics simulation (MD) protocol, we analyzed the process of channel inactivation. Then we applied a molecular docking procedure to predict the protein conformation suitable for the amiloride binding. To confirm the effect of its sole active blocker against the ASIC1 state transition route we studied the complex with another MD simulation run. Further experiments evaluated various compounds in the Enamine library that emerge with a detectable ASIC inhibitory activity. We performed a detailed analysis of the structural basis of ASIC1a inhibition by amiloride, using a combination of in silico approaches to visualize its interaction with the ion pore in the open state. An artificial activation (otherwise, expansion of the central pore) causes a complex modification of the channel structure, namely its transmembrane domain. The output protein conformations were used as a set of docking models, suitable for a high-throughput virtual screening of the Enamine chemical library. The outcome of the virtual screening was confirmed by electrophysiological assays with the best results shown for three hit compounds.

Computational design of a cyclic peptide that inhibits the CTLA4 immune checkpoint

Proteins involved in immune checkpoint pathways, such as CTLA4, PD1, and PD-L1, have become important targets for cancer immunotherapy; however, development of small molecule drugs targeting these pathways has proven difficult due to the nature of their protein-protein interfaces. Here, using a hierarchy of computational techniques, we design a cyclic peptide that binds CTLA4 and follow this with experimental verification of binding and biological activity, using bio-layer interferometry, cell culture, and a mouse tumor model. Beginning from a template excised from the X-ray structure of the CTLA4:B7-2 complex, we generate several peptide sequences using flexible docking and modeling steps. These peptides are cyclized head-to-tail to improve structural and proteolytic stability and screened using molecular dynamics simulation and MM-GBSA calculation. The standard binding free energies for shortlisted peptides are then calculated in explicit-solvent simulation using a rigorous multistep technique. The most promising peptide, cyc(EIDTVLTPTGWVAKRYS), yields the standard free energy -6.6 ± 3.5 kcal mol-1, which corresponds to a dissociation constant of ∼15 μmol L-1. The binding affinity of this peptide for CTLA4 is measured experimentally (31 ± 4 μmol L-1) using bio-layer interferometry. Treatment with this peptide inhibited tumor growth in a co-culture of Lewis lung carcinoma (LLC) cells and antigen primed T cells, as well as in mice with an orthotropic Lewis lung carcinoma allograft model.

Diagonal Interactions between Glutamate and Arginine Analogs with Varying Side-Chain Lengths in a β-Hairpin

Cross-strand interactions are important for the stability of β-sheet structures. Accordingly, cross-strand diagonal interactions between glutamate and arginine analogs with varying side-chain lengths were studied in a series of β-hairpin peptides. The peptides were analyzed by homonuclear two-dimensional nuclear magnetic resonance methods. The fraction folded population and folding free energy of the peptides were derived from the chemical shift data. The fraction folded population trends could be rationalized using the strand propensity of the constituting residues, which was not the case for the peptides with lysine analogs, highlighting the difference between the arginine analogs and lysine analogs. Double-mutant cycle analysis was used to derive the diagonal ion-pairing interaction energetics. The most stabilizing diagonal cross-strand interaction was between the shortest residues (i.e., Asp2-Agp9), most likely due to the least side-chain conformational penalty for ion-pair formation. The diagonal interaction energetics in this study involving the arginine analogs appears to be consistent with and extend beyond our understanding of diagonal ion-pairing interactions involving lysine analogs. The results should be useful for designing β-strand-containing molecules to affect biological processes such as amyloid formation and protein-protein interactions.

Design of Peptides that Fold and Self-Assemble on Graphite

The graphite-water interface provides a unique environment for polypeptides that generally favors ordered structures more than in solution. Therefore, systems consisting of designed peptides and graphitic carbon might serve as a convenient medium for controlled self-assembly of functional materials. Here, we computationally designed cyclic peptides that spontaneously fold into a β-sheet-like conformation at the graphite-water interface and self-assemble, and we subsequently observed evidence of such assembly by atomic force microscopy. Using a novel protocol, we screened nearly 2000 sequences, optimizing for formation of a unique folded conformation while discouraging unfolded or misfolded conformations. A head-to-tail cyclic peptide with the sequence GTGSGTGGPGGGCGTGTGSGPG showed the greatest apparent propensity to fold spontaneously, and this optimized sequence was selected for larger scale molecular dynamics simulations, rigorous free-energy calculations, and experimental validation. In simulations ranging from hundreds of nanoseconds to a few microseconds, we observed spontaneous folding of this peptide at the graphite-water interface under many different conditions, including multiple temperatures (295 and 370 K), with different initial orientations relative to the graphite surface, and using different molecular dynamics force fields (CHARMM and Amber). The thermodynamic stability of the folded conformation on graphite over a range of temperatures was verified by replica-exchange simulations and free-energy calculations. On the other hand, in free solution, the folded conformation was found to be unstable, unfolding in tens of picoseconds. Intermolecular hydrogen bonds promoted self-assembly of the folded peptides into linear arrangements where the peptide backbone exhibited a tendency to align along one of the six zigzag directions of the graphite basal plane. For the optimized peptide, atomic force microscopy revealed growth of single-molecule-thick linear patterns of 6-fold symmetry, consistent with the simulations, while no such patterns were observed for a control peptide with the same amino acid composition but a scrambled sequence.

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.

Comparative Analysis of Sulfonium-π, Ammonium-π, and Sulfur-π Interactions and Relevance to SAM-Dependent Methyltransferases

We report the measurement and analysis of sulfonium-π, thioether-π, and ammonium-π interactions in a β-hairpin peptide model system, coupled with computational investigation and PDB analysis. These studies indicated that the sulfonium-π interaction is the strongest and that polarizability contributes to the stronger interaction with sulfonium relative to ammonium. Computational studies demonstrate that differences in solvation of the trimethylsulfonium versus the trimethylammonium group also contribute to the stronger sulfonium-π interaction. In comparing sulfonium-π versus sulfur-π interactions in proteins, analysis of SAM- and SAH-bound enzymes in the PDB suggests that aromatic residues are enriched in close proximity to the sulfur of both SAM and SAH, but the populations of aromatic interactions of the two cofactors are not significantly different, with the exception of the Me-π interactions in SAM, which are the most prevalent interaction in SAM but are not possible for SAH. This suggests that the weaker interaction energies due to loss of the cation-π interaction in going from SAM to SAH may contribute to turnover of the cofactor.

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.

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.

Cyclic peptide engineered from phytocystatin inhibitory hairpin loop as an effective modulator of falcipains and potent antimalarial

Cystatins are classical competitive inhibitors of C1 family cysteine proteases (papain family). Phytocystatin superfamily shares high sequence homology and typical tertiary structure with conserved glutamine-valine-glycine (Q-X-V-X-G) loop blocking the active site of C1 proteases. Here, we develop a cysteine-bounded cyclic peptide (CYS-cIHL) and linear peptide (CYS-IHL), using the conserved inhibitory hairpin loop amino acid sequence. Using an in silico approach based on modeling, protein-peptide docking, molecular dynamics simulations and calculation of free energy of binding, we designed and validated inhibitory peptides against falcipain-2 (FP-2) and -3 (FP-3), cysteine proteases from the malarial parasite Plasmodium falciparum. Falcipains are critical hemoglobinases of P. falciparum that are validated targets for the development of antimalarial therapies. CYS-cIHL was able to bind with micromolar affinity to FP-2 and modulate its binding with its substrate, hemoglobin in in vitro and in vivo assays. CYS-cIHL could effectively block parasite growth and displayed antimalarial activity in culture assays with no cytotoxicity towards human cells. These results indicated that cyclization can substantially increase the peptide affinity to the target. Furthermore, this can be applied as an effective strategy for engineering peptide inhibitory potency against proteases.

Peptidyl inhibition of Spt4-Spt5: Protein-protein inhibitors for targeting the transcriptional pathway related to C9orf72 expansion repeats

Spt4/Spt5 is an useful target as it is likely a transcription factor that has implications for long non-coding RNA repeats related to frontotemporal dementia (FTD) found in the C9orf72 disease pathology. Inhibitors for Spt4/Spt5 using peptides as a starting point for assays as a means for developing small molecules, which could likely lead to therapeutic development for inhibition for Spt4/Spt5 with CNS characteristics. To elucidate the specific steps of identification and modification of key interacting residues from Spt4/Spt5 complex with further effect prediction, a set of different computational methods was applied. Newly characterized, theoretically derived peptides docked on Spt4/Spt5 models, based on X-ray crystallography sources, allowed us to complete molecular dynamics simulations and docking studies for peptide libraries that give us high confident set of peptides for use to screen for Spt4/Spt5 inhibition. Several peptides with increased specificity to the Spt4/Spt5 interface were found and can be screened in cell-based assays and enzymatic assays for peptide screens that lead to small molecule campaigns. Spt4/Spt5 comprises an attractive target for neurological diseases, and applying these peptides into a screening campaign will promote the goal of therapeutic searches for FTD drug discovery.

Defining the mobility range of a hinge-type connection using molecular dynamics and metadynamics

A designed disulfide-rich β-hairpin peptide that dimerizes spontaneously served as a hinge-type connection between proteins. Here, we analyze the range of dynamics of this hinge dimer with the aim of proposing new applications for the DNA-encodable peptide and establishing guidelines for the computational analysis of other disulfide hinges. A recent structural analysis based on nuclear magnetic resonance spectroscopy and ion mobility spectrometry revealed an averaged conformation in the hinge region which motivated us to investigate the dynamic behavior using a combination of molecular dynamics simulation, metadynamics and free energy surface analysis to characterize the conformational space available to the hinge. Principal component analysis uncovered two slow modes of the peptide, namely, the opening and closing motion and twisting of the two β-hairpins assembling the hinge. Applying a collective variable (CV) that mimics the first dominating mode, led to a major expansion of the conformational space. The description of the dynamics could be achieved by analysis of the opening angle and the twisting of the β-hairpins and, thus, offers a methodology that can also be transferred to other derivatives. It has been demonstrated that the hinge peptide’s lowest energy conformation consists of a large opening angle and strong twist but is separated by small energy barriers and can, thus, adopt a closed and untwisted structure. With the aim of proposing further applications for the hinge peptide, we simulated its behavior in the sterically congested environment of a four-helix bundle. Preliminary investigations show that one helix is pushed out and a three-helix bundle forms. The insights gained into the dynamics of the tetra-disulfide peptide and analytical guidelines developed in this study may contribute to the understanding of the structure and function of more complex hinge-type proteins, such as the IgG antibody family.

A set of conformationally well-defined L/D-peptide epitopes provides a serological bar code for autoantibody subtypes

Which conformational parameters lead to an antibody-affine peptide antigen? And in how many different conformations can we actually present the respective conformational epitope? To provide answers from a chemical point of view, we direct the bending and tethering of peptide backbones by the utilisation of a hydrophobic cluster, disulfides, and d-amino acids. Each mutation is employed pairwise on directly opposite sides of a β-hairpin. In combination, these synthetic modules guide the formation of complementary β-sheet-like structures, whereby the oppositely configured (l/d-)bi-disulfide pairs form with high regioselectivity. The conformational properties of the peptides are assessed by NMR spectroscopy and correlated with their antibody affinity in ELISA. From a pool of thus designed peptide antigens with distinctive complementary affinities against known rheumatoid arthritis (RA) autoantibodies, we select a set of epitopes for an immunoassay with sera of RA patients. We want to put emphasis on the idea, that the different conformational properties of the chosen antigens, containing the same epitope sequence, are mirrored in the distribution of autoantibody subtypes (or of the antibody polyclonality, respectively). Such directly comparable information can only be delivered by a set of peptides, rather than a single one. The hairpin-restriction technology of l/d-configured bi-disulfide amino acid pairs is not limited to RA but applicable to other shape-persistent hairpin motifs which are supposed to identify subgroups of protein receptors.

Effective formation of stable and versatile double-stranded β-sheets templated by a hydrogen-bonded duplex

An effective approach to construct stable and versatile double-stranded β-sheets composed of tetra- and penta-peptides through a hydrogen-bonded duplex template has been explored.

Structure-based design of ferritin nanoparticle immunogens displaying antigenic loops of Neisseria gonorrhoeae

Effective vaccines are urgently needed to combat gonorrhea, a common sexually transmitted bacterial infection, for which treatment options are diminishing due to rapid emergence of antibiotic resistance. We have used a rational approach to the development of gonorrhea vaccines, and genetically engineered nanoparticles to present antigenic peptides of Neisseria gonorrhoeae, the causative agent of gonorrhea. We hypothesized that the ferritin nanocage could be used as a platform to display an ordered array of N. gonorrhoeae antigenic peptides on its surface. MtrE, the outer membrane channel of the highly conserved gonococcal MtrCDE active efflux pump, is an attractive vaccine target due to its importance in protecting N. gonorrhoeae from host innate effectors and antibiotic resistance. Using computational approaches, we designed constructs that expressed chimeric proteins of the Helicobacter pylori ferritin and antigenic peptides that correspond to the two surface-exposed loops of N. gonorrhoeae MtrE. The peptides were inserted at the N terminus or in a surface-exposed ferritin loop between helices αA and αB. Crystal structures of the chimeric proteins revealed that the proteins assembled correctly into a 24-mer nanocage structure. Although the inserted N. gonorrhoeae peptides were disordered, it was clear that they were displayed on the nanocage surface, but with multiple conformations. Our results confirmed that the ferritin nanoparticle is a robust platform to present antigenic peptides and therefore an ideal system for rational design of immunogens.

Eight at one stroke – a synthetic tetra-disulfide peptide epitope

A tetra-disulfide peptide dimer, representing an antiparallel hinge, is synthesised without the need for orthogonal cysteine protecting groups.

Two Opposing d-Amino Acids Give Zigzag Hairpin Epitopes an Additional Kink to Create Antibody-Selective Peptide Antigens

We have developed peptides that are able to distinguish between subgroups of polyclonal antibodies. These β-hairpin peptides act as conformational epitopes with specific shape and flexibility; they have been analyzed by NMR and CD spectroscopy, and have been shown to identify known disease markers. As a standalone mini β-sheet, a hairpin is stabilized by alternating pairs of hydrogen-bonded and non-bonded amino acids on its two opposing peptide strands. A single d mutation disrupts this secondary structure, the correlated double-d mutation of two opposing amino acids compensates for this destabilizing effect. The designed kink was introduced into both hydrogen-bonded and -non-bonded positions of an all-l hairpin that is a known conformational epitope in molecular recognition. Our peptides enabled the discrimination of different human rheumatoid arthritis autoantibodies in an ELISA assay.

Synthesis of a biological active β-hairpin peptide by addition of two structural motifs

The idea of privileged scaffolds - that there seem to be more bioactive compounds found around some structures than others - is well established for small drug molecules, but has little significance for standalone peptide secondary structures whose adaptable shapes escape the definition of a 3D motif in the absence of a protein scaffold. Here, we joined two independent biological functions in a single highly restricted peptide to support the hypothesis that the β-hairpin shape is the common basis of two otherwise unrelated biological recognition processes. To achieve this, the hydrophobic cluster HWX₄LV from the decapeptide cyclic hairpin model peptide C¹-C¹⁰cyclo-CHWEGNKLVC was included in the bicyclic peptide 2. The designed β-hairpin peptide C⁴-C¹⁷, C⁸-C¹³bicyclo-KHQCHWECTZGRCRLVCGRSGS (2, Z=citrulline), serves, on the one hand, as a specific epitope for rheumatoid autoantibodies and, on the other hand, shows a not negligible antibiotic effect against the bacterial strain E. coli AS19.

Modifying the OPLS-AA force field to improve hydration free energies for several amino acid side chains using new atomic charges and an off-plane charge model for aromatic residues

The hydration free energies of amino acid side chains are an important determinant of processes that involve partitioning between different environments, including protein folding, protein complex formation, and protein-membrane interactions. Several recent papers have shown that calculated hydration free energies for polar and aromatic residues (Trp, His, Tyr, Asn, Gln, Asp, Glu) in several common molecular dynamics force fields differ significantly from experimentally measured values. We have attempted to improve the hydration energies for these residues by modifying the partial charges of the OPLS-AA force field based on natural population analysis of density functional theory calculations. The resulting differences between calculated hydration free energies and experimental results for the seven side chain analogs are less than 0.1 kcal/mol. Simulations of the synthetic Trp-rich peptide Trpzip2 show that the new charges lead to significantly improved geometries for interacting Trp-side chains. We also investigated an off-plane charge model for aromatic rings that more closely mimics their electronic configuration. This model results in an improved free energy of hydration for Trp and a somewhat altered benzene-sodium potential of mean force with a more favorable energy for direct benzene-sodium contact.

Disulfide-Mediated β-Strand Dimers: Hyperstable β-Sheets Lacking Tertiary Interactions and Turns

Disulfide bonds between cysteine residues are essential to the structure and folding of many proteins. Yet their role in the design of structured peptides and proteins has frequently been limited to use as intrachain covalent staples that reinforce existing structure or induce knot-like conformations. In β-hairpins, their placement at non-H-bonding positions across antiparallel strands has proven useful for achieving fully folded positive controls. Here we report a new class of designed β-sheet peptide dimers with strand-central disulfides as a key element. We have found that the mere presence of a disulfide bond near the middle of a short peptide chain is sufficient to nucleate some antiparallel β-sheet structure; addition of β-capping units and other favorable cross-strand interactions yield hyperstable sheets. Strand-central cystines were found to be superior to the best designed reversing turns in terms of nucleating β-sheet structure formation. We have explored the limitations and possibilities of this technique (the use of disulfides as sheet nucleators), and we provide a set of rules and rationales for the application and further design of disulfide-tethered "turnless" β-sheets.

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.

Tightening up the structure, lighting up the pathway: Application of molecular constraints and light to manipulate protein folding, self-assembly and function

Chemical cross-linking provides an effective avenue to reduce the conformational entropy of polypeptide chains and hence has become a popular method to induce or force structural formation in peptides and proteins. Recently, other types of molecular constraints, especially photoresponsive linkers and functional groups, have also found increased use in a wide variety of applications. Herein, we provide a concise review of using various forms of molecular strategies to constrain proteins, thereby stabilizing their native states, gaining insight into their folding mechanisms, and/or providing a handle to trigger a conformational process of interest with light. The applications discussed here cover a wide range of topics, ranging from delineating the details of the protein folding energy landscape to controlling protein assembly and function.

Asymmetric contribution of aromatic interactions stems from spatial positioning of the interacting aryl pairs in β-hairpins

Isolated aromatic interactions in designed octapeptide β-hairpin scaffolds display a near-universal T-shaped face-to-edge geometry in all positional permutations, with the exception of aryl-Trp interactions. The heterogeneous asymmetric indole ring of Trp competes for a "shielding" face geometry, which lowers the scaffold stability in FtE aryl-Trp pairs. Assessment of the contributions of aryl pairs (in the absence of solvent-driven interactions) to the overall β-hairpin structure reveals the superiority of Trp-Phe and Trp-Tyr contributions over the well-established scaffold stabilization by Trp-Trp.

Evolution of the first genetic cells and the universal genetic code: a hypothesis based on macromolecular coevolution of RNA and proteins

A qualitative hypothesis based on coevolution of protein and nucleic acid macromolecules was developed to explain the evolution of the first genetic cells, from the likely organic chemical-rich environment of early earth, through to the Last Universal Common Ancestor (LUCA). The evolution of the first genetic cell was divided into three phases, proto-genetic cells I, II and III, and the transition to each milestone is described, based on development of chemical cross-catalysis, bio-cross-catalysis, and the universal genetic code, respectively. Selection of macromolecular properties of both peptides and nucleic acids, in response to environmental factors, was likely to be a key aspect of early evolution. The development of hereditable nucleic acids with various key functions; translation, transcription and replication, is described. These functions are envisaged to have coevolved with protein enzymes, from simple organic precursors. Genetically heritable nucleotides may have developed after the local earth environment had cooled below 63 °C. Around this temperature G-C bases would have been preferentially utilized for nucleotide synthesis. Under these conditions RNA type nucleotides were then likely selected from a range of different types of nucleotide backbones through template-based synthesis. Initial development of the genetic coding system was simplified by the availability of proto-messenger RNA sequences that contained only G and C bases, and the need to encode only four amino acids. The step-wise addition of further amino acids to the code was predicted to parallel the growing metabolic complexity of the proto-genetic cell. On completion of this evolutionary process the proto-genetic cell is envisaged to have become the LUCA, the last common ancestor of bacteria, eukaryote and archaea domains. Key issues addressed by the model include: (a) the transition from non-hereditable random sequences of peptides and nucleic acids to specific proteins coded by hereditable nucleotide sequences, (b) the origin of homochiral amino acids and sugars, and (c) the mutation limits on the sizes of early nucleic acid genomes. The first genome was limited to a size of about 200 base pairs.

Computational screening and selection of cyclic peptide hairpin mimetics by molecular simulation and kinetic network models

Designing peptidomimetic compounds to have a preorganized structure in solution is highly nontrivial. To show how simulation-based approaches can help speed this process, we performed an extensive simulation study of designed cyclic peptide mimics of a β-hairpin from bacterial protein LapD involved in a protein-protein interaction (PPI) pertinent to bacterial biofilm formation. We used replica exchange molecular dynamics (REMD) simulation to screen 20 covalently cross-linked designs with varying stereochemistry and selected the most favorable of these for massively parallel simulation on Folding@home in explicit solvent. Markov state models (MSMs) built from the trajectory data reveal how subtle chemical modifications can have a significant effect on conformational populations, leading to the overall stabilization of the target structure. In particular, we identify a key steric interaction between a methyl substituent and a valine side chain that acts to allosterically shift population between native and near-native states, which could be exploited in future designs. Visualization of this mechanism is aided considerably by the tICA method, which identifies degrees of freedom most important in slow conformational transitions. The combination of quantitative detail and human comprehension provided by MSMs suggests such approaches will be increasingly useful for design.

How quickly can a β-hairpin fold from its transition state?

Understanding the structural nature of the free energy bottleneck(s) encountered in protein folding is essential to elucidating the underlying dynamics and mechanism. For this reason, several techniques, including Φ-value analysis, have previously been developed to infer the structural characteristics of such high free-energy or transition states. Herein we propose that one (or few) appropriately placed backbone and/or side chain cross-linkers, such as disulfides, could be used to populate a thermodynamically accessible conformational state that mimics the folding transition state. Specifically, we test this hypothesis on a model β-hairpin, Trpzip4, as its folding mechanism has been extensively studied and is well understood. Our results show that cross-linking the two β-strands near the turn region increases the folding rate by an order of magnitude, to about (500 ns)⁻¹, whereas cross-linking the termini results in a hyperstable β-hairpin that has essentially the same folding rate as the uncross-linked peptide. Taken together, these findings suggest that cross-linking is not only a useful strategy to manipulate folding free energy barriers, as shown in other studies, but also, in some cases, it can be used to stabilize a folding transition state analogue and allow for direct assessment of the folding process on the downhill side of the free energy barrier. The calculated free energy landscape of the cross-linked Trpzip4 also supports this picture. An empirical analysis further suggests, when folding of β-hairpins does not involve a significant free energy barrier, the folding time (τ) follows a power law dependence on the number of hydrogen bonds to be formed (n_H), namely, τ = τ₀n_H^α, with τ₀ = 20 ns and α = 2.3.

Comparative analysis of cross strand aromatic–Phe interactions in designed peptide β-hairpins

Examination of the preferential interaction geometries of the aromatic amino acids Phe, Tyr and Trp with the benzyl ring of Phe in designed octapeptide hairpin scaffolds reveals stabilizing contributions of a Trp–Phe pair, even in amphipathic solvents.

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.

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.

Investigating the function of single nucleotide polymorphisms in the CTSB gene: a computational approach

Aim: Recent genome-wide association studies have revealed large numbers of single nucleotide polymorphisms (SNPs) related to Alzheimer’s disease. Here, we have investigated the gene CTSB, which plays a crucial role in encoding CTSB, a lysosomal cysteine proteinase protein. CTSB is also involved in the proteolytic processing of amyloid precursor protein (APP), which is believed to be a causative factor in Alzheimer’s disease. Materials & methods: Several bioinformatics algorithms such as, Sorting Intolerant from Tolerant (SIFT), Polymorphism Phenotyping (PolyPhen) and CUPSAT could identify the synonymous SNPs and nonsynonymous SNPs (nsSNPs), which are predicted to be deleterious and nondeleterious, respectively. Similar tools were used to predict the impact of single amino acid substitutions on CTSB protein activity. The FASTSNP server and UTRscan were used to predict the influence on splicing regulations. The stability and solvent-accessible surface area of modeled mutated proteins were analyzed using PBEQ solver and NetASA view. Furthermore, the DSP program was used to determine the secondary structures of the modeled protein. Results: A total of 999 SNPs in CTSB were retrieved from the SNP database; 55 nsSNPs, 35 synonymous SNPs, 165 mRNA were found in the 3´untranslated region SNPs, 12 SNPs were found in the 5´untranslated region in addition to 732 intronic SNPs. Potential functions of SNPs in the CTSB gene were identified using different web servers. For example, SIFT, PolyPhen and CUPSAT servers predicted ten nsSNPs to be intolerant, three nsSNPs to be damaging and eight nsSNPs to have the potential to destabilize protein structure. The FASTSNP server predicted 12 SNPs to influence splicing regulation, whereas two SNPs could predict a risk in the range of 3–4 (medium to high). Furthermore, mutant proteins were modeled and the total energy values were compared with the native CTSB protein. It was observed that on the surface of the protein, a mutation from threonine to serine at position 235 (rs17573) caused the greatest impact on stability. Conclusion: The genome-wide association studies database has already found rs7003814 of the CTSB gene reported against Alzheimer’s disease. Our study demonstrates the presence of other deleterious nsSNPs, which may play a crucial role in predicting Alzheimer’s disease risk.

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.

Three types of induced tryptophan optical activity compared in model dipeptides: theory and experiment

The tryptophan (Trp) aromatic residue in chiral matrices often exhibits a large optical activity and thus provides valuable structural information. However, it can also obscure spectral contributions from other peptide parts. To better understand the induced chirality, electronic circular dichroism (ECD), vibrational circular dichroism (VCD), and Raman optical activity (ROA) spectra of Trp-containing cyclic dipeptides c-(Trp-X) (where X = Gly, Ala, Trp, Leu, nLeu, and Pro) are analyzed on the basis of experimental spectra and density functional theory (DFT) computations. The results provide valuable insight into the molecular conformational and spectroscopic behavior of Trp. Whereas the ECD is dominated by Trp π-π* transitions, VCD is dominated by the amide modes, well separated from minor Trp contributions. The ROA signal is the most complex. However, an ROA marker band at 1554 cm(-1) indicates the local χ(2) angle value in this residue, in accordance with previous theoretical predictions. The spectra and computations also indicate that the peptide ring is nonplanar, with a shallow potential so that the nonplanarity is primarily induced by the side chains. Dispersion-corrected DFT calculations provide better results than plain DFT, but comparison with experiment suggests that they overestimate the stability of the folded conformers. Molecular dynamics simulations and NMR results also confirm a limited accuracy of the dispersion-DFT model in nonaqueous solvents. Combination of chiral spectroscopies with theoretical analysis thus significantly enhances the information that can be obtained from the induced chirality of the Trp aromatic residue.

Using thioamides to site-specifically interrogate the dynamics of hydrogen bond formation in β-sheet folding

Thioamides are sterically almost identical to their oxoamide counterparts, but they are weaker hydrogen bond acceptors. Therefore, thioamide amino acids are excellent candidates for perturbing the energetics of backbone-backbone H-bonds in proteins and hence should be useful in elucidating protein folding mechanisms in a site-specific manner. Herein, we validate this approach by applying it to probe the dynamic role of interstrand H-bond formation in the folding kinetics of a well-studied β-hairpin, tryptophan zipper. Our results show that reducing the strength of the peptide's backbone-backbone H-bonds, except the one directly next to the β-turn, does not change the folding rate, suggesting that most native interstrand H-bonds in β-hairpins are formed only after the folding transition state.

Molecular dynamics studies of β-hairpin folding with the presence of the sodium ion

Metal ions are ubiquitous in protein systems and play a significant role during their folding processes. Nineteen independent structures were determined for the Na(+)/β-hairpin interacting systems, and their folding pathways are different. (i) For Na(S47), the turn is rapidly shaped with the help of Na(+) and acts as the folding nucleus for the rest regions. Two intermediate states are observed and the resulted structure is the most folded. (ii) For Na(B41), Na(B52), Na(B54), Na(S55) and Na(B56), the inclusive Na(+) ions are anchored by β-strands. The local structures around the Na(+) ions and the turn regions fold simultaneously and serve as two independent folding nuclei. (iii) The other systems have no folding nuclei and correspond to low-folded structures. Long-range electrostatic interactions contribute a lot to the folding, especially from the four negatively charged residues (Glu42, Asp46, Asp47 and Glu56). The initial positions of the Na(+) ions are largely responsible for the different folding behaviors. The interactions with sidechain- rather than backbone-O atoms generally lead to more compact structures. Another factor affecting the folding is whether the O atoms are associated with native H-bonds, and those involved show decreased affinities to metal ions. The addition of water solvent does not induce obvious folding and conformational transitions to the Na(+)/β-hairpin interacting systems.

Impact ofβ-Turn Sequence onβ-Hairpin Dynamics Studied with Infrared-Detected Temperature Jump

Folding dynamics forβ-structure loss and disordered structure gain were studied in a modelβ-hairpin peptide based on Cochran’s tryptophan zipper peptide Trpzip2, but with an altered Thr-Gly (TG) turn sequence, that is, SWTWETGKWTWK, using laser-induced temperature-jump (T-jump) kinetics with IR detection. As has been shown previously, the TG turn sequence reduces the thermodynamicβ-hairpin stability as compared to the Asn-Gly sequence used in Trpzip2 (TZ2-NG). In this study, we found that the TG-turn slows down the overall relaxation dynamics as compared to TZ2-NG, which were studied at higher temperatures where the time constants show little difference between relaxation of theβ-strand and the disordered conformation. These time constants become equivalent at lower temperatures for TZ2-TG than was seen for TZ2-NG. The correlation of thermodynamic stability and rates of relaxation suggests that the change from NG to TG turn results in a slowing of folding, lowerkf, with less change of the unfolding rate,ku, assuming two state behavior at higher temperatures.

β-hairpin peptide that targets vascular endothelial growth factor (VEGF) receptors: design, NMR characterization, and biological activity

VEGF receptors have been the target of intense research aimed to develop molecules able to inhibit or stimulate angiogenesis. Based on the x-ray structure of the complex placental growth factor-VEGF receptor 1(D2), we designed a VEGF receptor-binding peptide reproducing the placental growth factor β-hairpin region Gln(87)-Val(100) that is involved in receptor recognition. A conformational analysis showed that the designed peptide adopts the expected fold in pure water. Moreover, a combination of NMR interaction analysis and cell binding studies were used to demonstrate that the peptide targets VEGF receptors. The VEGF receptor 1(D2)-interacting residues were characterized at the molecular level, and they correspond to the residues recognizing the placental growth factor sequence Gln(87)-Val(100). Finally, the peptide biological activity was characterized in vitro and in vivo, and it showed a VEGF-like behavior. Indeed, the peptide activated VEGF-dependent intracellular pathways, induced endothelial cell proliferation and rescue from apoptosis, and promoted angiogenesis in vivo. This compound is one of the few peptides known with proangiogenic activity, which makes it a candidate for the development of a novel peptide-based drug for medical applications in therapeutic angiogenesis.

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.

β-Sheet 13C structuring shifts appear only at the H-bonded sites of hairpins

The (13)C chemical shifts measured for designed β-hairpins indicate that the structuring shifts (upfield for Cα and C', downfield for Cβ) previously reported as diagnostic for β-structuring in proteins appear only at the H-bonded strand residues. The resulting periodicity of structuring shift magnitudes is not, however, a consequence of H-bonding status; rather, it reflects a previously unrecognized alternation in the backbone torsion angles of β-strands. This feature of hairpins is also likely to be present in proteins. The study provides reference values for the expectation shifts for (13)C sites in β-structures that should prove useful in the characterization of the folding equilibria of β-sheet models.