Predicting Order and Disorder for β-Peptide Foldamers in Water

Quest for New Generation Biocompatible Materials: Tailoring β-Peptide Structure and Interactions via Synergy of Experiments and Modelling

Peptide-based self-assembly has been used to produce a wide range of nanostructures. While most of these systems involve self-assembly of α-peptides, more recently β-peptides have also been shown to undergo supramolecular self-assembly, and have been used to produce materials for applications in tissue engineering, cell culture and drug delivery. In order to engineer new materials with specific structure and function, theoretical molecular modelling can provide significant insights into the collective balance of non-covalent interactions that drive the self-assembly and determine the structure of the resultant supramolecular materials under different conditions. However, this approach has only recently become feasible for peptide-based self-assembled nanomaterials, particularly those that incorporate non α-amino acids. This perspective provides an overview of the challenges associated with computational modelling of the self-assembly of β-peptides and the recent success using a combination of experimental and computational techniques to provide insights into the self-assembly mechanisms and fully atomistic models of these new biocompatible materials.

Conformational landscapes of artificial peptides predicted by various force fields: are we ready to simulate β-amino acids?

With the introduction of artificial peptides as antimicrobial agents and organic catalysts, numerous efforts have been made to design foldamers with desirable structures and functions. Computational tools are a helpful proxy for revealing the dynamic structures at atomic resolution and understanding foldamer's complex structure-function relationships. However, the performance of conventional force fields in predicting the structures of artificial peptides has not been systematically evaluated. In this study, we critically assessed three popular force fields, AMBER ff14SB, CHARMM36m, and OPLS-AA/L, in predicting conformational propensities of a β-peptide foldamer at monomer and hexamer levels. Simulation results were compared to those obtained from quantum chemistry calculations and experimental data. We also utilised replica exchange molecular dynamics simulations to investigate the energy landscape of each force field and assess the similarities and differences between force fields. We compared different solvent systems in the AMBER ff14SB and CHARMM36m frameworks and confirmed the unanimous role of hydrogen bonds in shaping energy landscapes. We anticipate that our data will pave the way for further improvements to force fields and for understanding the role of solvents in peptide folding, crystallisation, and engineering.

Folding Coarse-Grained Oligomer Models with PyRosetta

Non-biological foldamers are a promising class of macromolecules that share similarities to classical biopolymers such as proteins and nucleic acids. Currently, designing novel foldamers is a non-trivial process, often involving many iterations of trial synthesis and characterization until folded structures are observed. In this work, we aim to tackle these foldamer design challenges using computational modeling techniques. We developed CG PyRosetta, an extension to the popular protein folding python package, PyRosetta, which introduces coarse-grained (CG) residues into PyRosetta, enabling the folding of toy CG foldamer models. Although these models are simplified, they can help explore overarching physical hypotheses about how oligomers can form. Through systematic variation of CG parameters in these models, we can investigate various folding hypotheses at the CG scale to inform the design process of new foldamer chemistries. In this study, we demonstrate CG PyRosetta's ability to identify minimum energy structures with a diverse structural search over a range of simple models, as well as two hypothesis-driven parameter scans investigating the effects of side-chain size and internal backbone angle on secondary structures. We are able to identify several types of secondary structures from single- and double-helices to sheet-like and knot-like structures. We show how side-chain size and backbone bond angle both play an important role in the structure and energetics of these toy models. Optimal side-chain sizes promote favorable packing of side chains, while specific backbone bond angles influence the specific helix type found in folded structures.

Controlling the conformational stability of coiled-coil peptides with a single stereogenic center of a peripheral β-amino acid residue

The key issue in the research on foldamers remains the understanding of the relationship between the monomers structure and conformational properties at the oligomer level. In peptidomimetic foldamers, the main goal of which is to mimic the structure of proteins, a main challenge is still better understanding of the folding of peptides and the factors that influence their conformational stability. We probed the impact of the modification of the peptide periphery with trans- and cis-2-aminocyclopentanecarboxylic acid (ACPC) on the structure and stability of the model coiled-coil using circular dichroism (CD), analytical ultracentrifugation (AUC) and two-dimensional nuclear magnetic resonance spectroscopy (2D NMR). Although, trans-ACPC and cis-ACPC-containing mutants differ by only one peripheral stereogenic center, their conformational stability is strikingly different.

Improved Modeling of Peptidic Foldamers Using a Quantum Chemical Parametrization Based on Torsional Minimum Energy Path Matching

The increasing interest in novel foldamer constructs demands an accurate computational treatment on an extensive timescale. However, it is still a challenge to derive a force field (FF) that can reproduce the experimentally known fold while also allowing the spontaneous exploration of other structures. Here, aiming at a realistic reproduction of backbone torsional barriers, the relevant proper dihedrals of acyclic β2-, β3- and β2,3-amino acids were added to the CHARMM FF and optimized using a novel, self-consistent iterative procedure based on quantum chemical relaxed scans. The new FF was validated by molecular dynamics simulations on three acyclic peptides. While they resided most of the time in their preferred fold (>80 % in helices and >50 % in hairpin), they also visited other conformations. Owing to the CHARMM36m-consistent parametrization, the proposed extension is suitable for exploring new foldamer structures and assemblies, and their interactions with diverse biomolecules.

Identifying the Coiled-Coil Triple Helix Structure of β-Peptide Nanofibers at Atomic Resolution

Peptide self-assembly represents a powerful bottom-up approach to the fabrication of nanomaterials. β³-Peptides are non-natural peptides composed entirely of β-amino acids, which have an extra methylene in the backbone, and we reported fibers derived from the self-assembly of β³-peptides that adopt 14-helical structures. β³-Peptide assemblies represent a class of stable nanomaterials that can be used to generate bio- and magneto-responsive materials with proteolytic stability. However, the three-dimensional structure of many of these materials remains unknown. To develop structure-based criteria for the design of β³-peptide-based biomaterials with tailored function, we investigated the structure of a tri-β³-peptide nanoassembly by molecular dynamics simulations and X-ray fiber diffraction analysis. Diffraction data was collected from aligned fibrils formed by Ac-β³[LIA] in water and used to inform and validate the model structure. Models with 3-fold radial symmetry resulted in stable fibers with a triple-helical coiled-coil motif and measurable helical pitch and periodicity. The fiber models revealed a hydrophobic core and twist along the fiber axis arising from a maximization of contacts between hydrophobic groups of adjacent tripeptides on the solvent-exposed fiber surface. These atomic structures of macroscale fibers derived from β³-peptide-based materials provide valuable insight into the effects of the geometric placement of the side chains and the influence of solvent on the core fiber structure which is perpetuated in the superstructure morphology.

Helix-loop-helix peptide foldamers and their use in the construction of hydrolase mimetics

De novo designed helix-loop-helix peptide foldamers containing cis-2-aminocyclopentanecarboxylic acid residues were evaluated for their conformational stability and possible use in enzyme mimetic development. The correlation between hydrogen bond network size and conformational stability was demonstrated through CD and NMR spectroscopies. Molecules incorporating a Cys/His/Glu triad exhibited enzyme-like hydrolytic activity.

Therapeutic Potential of Foldamers: From Chemical Biology Tools To Drug Candidates?

Over the past decade, foldamers have progressively emerged as useful architectures to mimic secondary structures of proteins. Peptidic foldamers, consisting of various amino acid based backbones, have been the most studied from a therapeutic perspective, while polyaromatic foldamers have barely evolved from their nascency and remain perplexing for medicinal chemists due to their poor drug-like nature. Despite these limitations, this compound class may still offer opportunities to study challenging targets or provide chemical biology tools. The potential of foldamer drug candidates reaching the clinic is still a stretch. Nevertheless, advances in the field have demonstrated their potential for the discovery of next generation therapeutics. In this perspective, the current knowledge of foldamers is reviewed in a drug discovery context. Recent advances in the early phases of drug discovery including hit finding, target validation, and optimization and molecular modeling are discussed. In addition, challenges and focus areas are debated and gaps highlighted.

Induced folding of protein-sized foldameric β-sandwich models with core β-amino acid residues

The mimicry of protein-sized β-sheet structures with unnatural peptidic sequences (foldamers) is a considerable challenge. In this work, the de novo designed betabellin-14 β-sheet has been used as a template, and α→β residue mutations were carried out in the hydrophobic core (positions 12 and 19). β-Residues with diverse structural properties were utilized: Homologous β(3) -amino acids, (1R,2S)-2-aminocyclopentanecarboxylic acid (ACPC), (1R,2S)-2-aminocyclohexanecarboxylic acid (ACHC), (1R,2S)-2-aminocyclohex-3-enecarboxylic acid (ACEC), and (1S,2S,3R,5S)-2-amino-6,6-dimethylbicyclo[3.1.1]heptane-3-carboxylic acid (ABHC). Six α/β-peptidic chains were constructed in both monomeric and disulfide-linked dimeric forms. Structural studies based on circular dichroism spectroscopy, the analysis of NMR chemical shifts, and molecular dynamics simulations revealed that dimerization induced β-sheet formation in the 64-residue foldameric systems. Core replacement with (1R,2S)-ACHC was found to be unique among the β-amino acid building blocks studied because it was simultaneously able to maintain the interstrand hydrogen-bonding network and to fit sterically into the hydrophobic interior of the β-sandwich. The novel β-sandwich model containing 25 % unnatural building blocks afforded protein-like thermal denaturation behavior.