Toward structure prediction of cyclic peptides

Accurate Structure Prediction for Cyclic Peptides Containing Proline Residues with High-Temperature Molecular Dynamics

Cyclic peptides (CPs) are emerging as promising drug candidates. Numerous natural CPs and their analogs are effective therapeutics against various diseases. Notably, many of them contain peptidyl cis-prolyl bonds. Due to the high rotational barrier of peptide bonds, conventional molecular dynamics simulations struggle to effectively sample the cis/trans-isomerization of peptide bonds. Previous studies have highlighted the high accuracy of the residue-specific force field (RSFF) and the high sampling efficiency of high-temperature molecular dynamics (high-T MD). Herein, we propose a protocol that combines high-T MD with RSFF2C and a recently developed reweighting method based on probability densities for accurate structure prediction of proline-containing CPs. Our method successfully predicted 19 out of 23 CPs with the backbone rmsd < 1.0 Å compared to X-ray structures. Furthermore, we performed high-T MD and density reweighting on the sunflower trypsin inhibitor (SFTI-1)/trypsin complex to demonstrate its applicability in studying CP-complexes containing cis-prolines. Our results show that the conformation of SFTI-1 in aqueous solution is consistent with its bound conformation, potentially facilitating its binding.

Gaussian accelerated molecular dynamics simulations facilitate prediction of the permeability of cyclic peptides

Despite their widespread use as therapeutics, clinical development of small molecule drugs remains challenging. Among the many parameters that undergo optimization during the drug development process, increasing passive cell permeability (i.e., log(P)) can have some of the largest impact on potency. Cyclic peptides (CPs) have emerged as a viable alternative to small molecules, as they retain many of the advantages of small molecules (oral availability, target specificity) while being highly effective at traversing the plasma membrane. However, the relationship between the dominant conformations that typify CPs in an aqueous versus a membrane environment and cell permeability remain poorly characterized. In this study, we have used Gaussian accelerated molecular dynamics (GaMD) simulations to characterize the effect of solvent on the free energy landscape of lariat peptides, a subset of CPs that have recently shown potential for drug development (Kelly et al., JACS 2021). Differences in the free energy of lariat peptides as a function of solvent can be used to predict permeability of these molecules, and our results show that permeability is most greatly influenced by N-methylation and exposure to solvent. Our approach lays the groundwork for using GaMD as a way to virtually screen large libraries of CPs and drive forward development of CP-based therapeutics.

Comprehensive Evaluation of 10 Docking Programs on a Diverse Set of Protein-Cyclic Peptide Complexes

Cyclic peptides have emerged as a highly promising class of therapeutic molecules owing to their favorable pharmacokinetic properties, including stability and permeability. Currently, many clinically approved cyclic peptides are derived from natural products or their derivatives, and the development of molecular docking techniques for cyclic peptide discovery holds great promise for expanding the applications and potential of this class of molecules. Given the availability of numerous docking programs, there is a pressing need for a systematic evaluation of their performance, specifically on protein-cyclic peptide systems. In this study, we constructed an extensive benchmark data set called CPSet, consisting of 493 protein-cyclic peptide complexes. Based on this data set, we conducted a comprehensive evaluation of 10 docking programs, including Rosetta, AutoDock CrankPep, and eight protein-small molecule docking programs (i.e., AutoDock, AudoDock Vina, Glide, GOLD, LeDock, rDock, MOE, and Surflex). The evaluation encompassed the assessment of the sampling power, docking power, and scoring power of these programs. The results revealed that all of the tested protein-small molecule docking programs successfully sampled the binding conformations when using the crystal conformations as the initial structures. Among them, rDock exhibited outstanding performance, achieving a remarkable 94.3% top-100 sampling success rate. However, few programs achieved successful predictions of the binding conformations using tLEaP-generated conformations as the initial structures. Within this scheme, AutoDock CrankPep yielded the highest top-100 sampling success rate of 29.6%. Rosetta's scoring function outperformed the others in selecting optimal conformations, resulting in an impressive top-1 docking success rate of 87.6%. Nevertheless, all the tested scoring functions displayed limited performance in predicting binding affinity, with MOE@Affinity dG exhibiting the highest Pearson's correlation coefficient of 0.378. It is therefore suggested to use an appropriate combination of different docking programs for given tasks in real applications. We expect that this work will offer valuable insights into selecting the appropriate docking programs for protein-cyclic peptide complexes.

Addressing Challenges of Macrocyclic Conformational Sampling in Polar and Apolar Solvents: Lessons for Chameleonicity

We evaluated a workflow to reliably sample the conformational space of a set of 47 peptidic macrocycles. Starting from SMILES strings, we use accelerated molecular dynamics simulations to overcome high energy barriers, in particular, the cis-trans isomerization of peptide bonds. We find that our approach performs very well in polar solvents like water and dimethyl sulfoxide. Interestingly, the protonation state of a secondary amine in the ring only slightly influences the conformational ensembles of our test systems. For several of the macrocycles, determining the conformational distribution in chloroform turns out to be considerably more challenging. Especially, the choice of partial charges crucially influences the ensembles in chloroform. We address these challenges by modifying initial structures and the choice of partial charges. Our results suggest that special care has to be taken to understand the configurational distribution in apolar solvents, which is a key step toward a reliable prediction of membrane permeation of macrocycles and their chameleonic properties.

Computational Prediction of Cyclic Peptide Structural Ensembles and Application to the Design of Keap1 Binders

The Nrf2 transcription factor is a master regulator of the cellular response to oxidative stress, and Keap1 is its primary negative regulator. Activating Nrf2 by inhibiting the Nrf2-Keap1 protein-protein interaction has shown promise for treating cancer and inflammatory diseases. A loop derived from Nrf2 has been shown to inhibit Keap1 selectively, especially when cyclized, but there are no reliable design methods for predicting an optimal macrocyclization strategy. In this work, we employed all-atom, explicit-solvent molecular dynamics simulations with enhanced sampling methods to predict the relative degree of preorganization for a series of peptides cyclized with a set of bis-thioether "staples". We then correlated these predictions to experimentally measured binding affinities for Keap1 and crystal structures of the cyclic peptides bound to Keap1. This work showcases a computational method for designing cyclic peptides by simulating and comparing their entire solution-phase ensembles, providing key insights into designing cyclic peptides as selective inhibitors of protein-protein interactions.

Ensemble Generation for Linear and Cyclic Peptides Using a Reservoir Replica Exchange Molecular Dynamics Implementation in GROMACS

The profile of shapes presented by a cyclic peptide modulates its therapeutic efficacy and is represented by the ensemble of its sampled conformations. Although some algorithms excel at creating a diverse ensemble of cyclic peptide conformations, they seldom address the entropic contribution of flexible conformations and often have significant practical difficulty producing an ensemble with converged and reliable thermodynamic properties. In this study, an accelerated molecular dynamics (MD) method, namely, reservoir replica exchange MD (R-REMD or Res-REMD), was implemented in GROMACS ver. 4.6.7 and benchmarked on two small cyclic peptide model systems: a cyclized furin cleavage site of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike (cyclo-(CGPRRARSG)) and oxytocin (disulfide-bonded CYIQNCPLG). Additionally, we also benchmarked Res-REMD on alanine dipeptide and Trpzip2 to demonstrate its validity and efficiency over REMD. For Trpzip2, Res-REMD coupled with an umbrella-sampling-derived reservoir generated similar folded fractions as regular REMD but on a much faster time scale. For cyclic peptides, Res-REMD appeared to be marginally faster than REMD in ensemble generation. Finally, Res-REMD was more effective in sampling rare events such as trans to cis peptide bond isomerization. We provide a GitHub page with the modified GROMACS source code for running Res-REMD at https://github.com/PlotkinLab/Reservoir-REMD.

Docking cyclic peptides formed by a disulfide bond through a hierarchical strategy

MOTIVATION

Cyclization is a common strategy to enhance the therapeutic potential of peptides. Many cyclic peptide drugs have been approved for clinical use, in which the disulfide-driven cyclic peptide is one of the most prevalent categories. Molecular docking is a powerful computational method to predict the binding modes of molecules. For protein-cyclic peptide docking, a big challenge is considering the flexibility of peptides with conformers constrained by cyclization.

RESULTS

Integrating our efficient peptide 3D conformation sampling algorithm MODPEP2.0 and knowledge-based scoring function ITScorePP, we have proposed an extended version of our hierarchical peptide docking algorithm, named HPEPDOCK2.0, to predict the binding modes of the peptide cyclized through a disulfide against a protein. Our HPEPDOCK2.0 approach was extensively evaluated on diverse test sets and compared with the state-of-the-art cyclic peptide docking program AutoDock CrankPep (ADCP). On a benchmark dataset of 18 cyclic peptide-protein complexes, HPEPDOCK2.0 obtained a native contact fraction of above 0.5 for 61% of the cases when the top prediction was considered, compared with 39% for ADCP. On a larger test set of 25 cyclic peptide-protein complexes, HPEPDOCK2.0 yielded a success rate of 44% for the top prediction, compared with 20% for ADCP. In addition, HPEPDOCK2.0 was also validated on two other test sets of 10 and 11 complexes with apo and predicted receptor structures, respectively. HPEPDOCK2.0 is computationally efficient and the average running time for docking a cyclic peptide is about 34 min on a single CPU core, compared with 496 min for ADCP. HPEPDOCK2.0 will facilitate the study of the interaction between cyclic peptides and proteins and the development of therapeutic cyclic peptide drugs.

AVAILABILITY AND IMPLEMENTATION

http://huanglab.phys.hust.edu.cn/hpepdock/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Evaluating the Conformations and Dynamics of Peptoid Macrocycles

Peptoid macrocycles are versatile and chemically diverse peptidomimetic oligomers. However, the conformations and dynamics of these macrocycles have not been evaluated comprehensively and require extensive further investigation. Recent studies indicate that two degrees of freedom, and four distinct conformations, adequately describe the behavior of each monomer backbone unit in most peptoid oligomers. On the basis of this insight, we conducted molecular dynamics simulations of model macrocycles using an exhaustive set of idealized possible starting conformations. Simulations of various sizes of peptoid macrocycles yielded a limited set of populated conformations. In addition to reproducing all relevant experimentally determined conformations, the simulations accurately predicted a cyclo-octamer conformation for which we now present the first experimental observation. Sets of three adjacent dihedral angles (ϕ_i, ψ_i, ω_i+1) exhibited correlated crankshaft motions over the course of simulation for peptoid macrocycles of six residues and larger. These correlated motions may occur in the form of an inversion of one amide bond and the concerted rotation of the preceding ϕ and ψ angles to their mirror-image conformation, a variation on "crankshaft flip" motions studied in polymers and peptides. The energy landscape of these peptoid macrocycles can be described as a network of conformations interconnected by transformations of individual crankshaft flips. For macrocycles of up to eight residues, our mapping of the landscape is essentially complete.

Cyclization and Docking Protocol for Cyclic Peptide-Protein Modeling Using HADDOCK2.4

An emerging class of therapeutic molecules are cyclic peptides with over 40 cyclic peptide drugs currently in clinical use. Their mode of action is, however, not fully understood, impeding rational drug design. Computational techniques could positively impact their design, but modeling them and their interactions remains challenging due to their cyclic nature and their flexibility. This study presents a step-by-step protocol for generating cyclic peptide conformations and docking them to their protein target using HADDOCK2.4. A dataset of 30 cyclic peptide-protein complexes was used to optimize both cyclization and docking protocols. It supports peptides cyclized via an N- and C-terminus peptide bond and/or a disulfide bond. An ensemble of cyclic peptide conformations is then used in HADDOCK to dock them onto their target protein using knowledge of the binding site on the protein side to drive the modeling. The presented protocol predicts at least one acceptable model according to the critical assessment of prediction of interaction criteria for each complex of the dataset when the top 10 HADDOCK-ranked single structures are considered (100% success rate top 10) both in the bound and unbound docking scenarios. Moreover, its performance in both bound and fully unbound docking is similar to the state-of-the-art software in the field, Autodock CrankPep. The presented cyclization and docking protocol should make HADDOCK a valuable tool for rational cyclic peptide-based drug design and high-throughput screening.

Efficient 3D conformer generation of cyclic peptides formed by a disulfide bond

Cyclic peptides formed by disulfide bonds have been one large group of common drug candidates in drug development. Structural information of a peptide is essential to understand its interaction with its target. However, due to the high flexibility of peptides, it is difficult to sample the near-native conformations of a peptide. Here, we have developed an extended version of our MODPEP approach, named MODPEP2.0, to fast generate the conformations of cyclic peptides formed by a disulfide bond. MODPEP2.0 builds the three-dimensional (3D) structures of a cyclic peptide from scratch by assembling amino acids one by one onto the cyclic fragment based on the constructed rotamer and cyclic backbone libraries. Being tested on a data set of 193 diverse cyclic peptides, MODPEP2.0 obtained a considerable advantage in both accuracy and computational efficiency, compared with other sampling algorithms including PEP-FOLD, ETKDG, and modified ETKDG (mETKDG). MODPEP2.0 achieved a high sampling accuracy with an average C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document}α RMSD of 2.20 Å and 1.66 Å when 10 and 100 conformations were considered, respectively, compared with 3.41 Å and 2.62 Å for PEP-FOLD, 3.44 Å and 3.16 Å for ETKDG, 3.09 Å and 2.72 Å for mETKDG. MODPEP2.0 also reproduced experimental peptide structures for 81.35% of the test cases when an ensemble of 100 conformations were considered, compared with 54.95%, 37.50% and 50.00% for PEP-FOLD, ETKDG, and mETKDG. MODPEP2.0 is computationally efficient and can generate 100 peptide conformations in one second. MODPEP2.0 will be useful in sampling cyclic peptide structures and modeling related protein-peptide interactions, facilitating the development of cyclic peptide drugs.

Computational Tools and Strategies to Develop Peptide-Based Inhibitors of Protein-Protein Interactions

Protein-protein interactions play crucial and subtle roles in many biological processes and modifications of their fine mechanisms generally result in severe diseases. Peptide derivatives are very promising therapeutic agents for modulating protein-protein associations with sizes and specificities between those of small compounds and antibodies. For the same reasons, rational design of peptide-based inhibitors naturally borrows and combines computational methods from both protein-ligand and protein-protein research fields. In this chapter, we aim to provide an overview of computational tools and approaches used for identifying and optimizing peptides that target protein-protein interfaces with high affinity and specificity. We hope that this review will help to implement appropriate in silico strategies for peptide-based drug design that builds on available information for the systems of interest.

Molecular Dynamics Simulations of Transmembrane Cyclic Peptide Nanotubes Using Classical Force Fields, Hydrogen Mass Repartitioning, and Hydrogen Isotope Exchange Methods: A Critical Comparison

Self-assembled cyclic peptide nanotubes with alternating D- and L-amino acid residues in the sequence of each subunit have attracted a great deal of attention due to their potential for new nanotechnology and biomedical applications, mainly in the field of antimicrobial peptides. Molecular dynamics simulations can be used to characterize these systems with atomic resolution at different time scales, providing information that is difficult to obtain via wet lab experiments. However, the performance of classical force fields typically employed in the simulation of biomolecules has not yet been extensively tested with this kind of highly constrained peptide. Four different classical force fields (AMBER, CHARMM, OPLS, and GROMOS), using a nanotube formed by eight D,L-α-cyclic peptides inserted into a lipid bilayer as a model system, were employed here to fill this gap. Significant differences in the pseudo-cylindrical cavities formed by the nanotubes were observed, the most important being the diameter of the nanopores, the number and location of confined water molecules, and the density distribution of the solvent molecules. Furthermore, several modifications were performed on GROMOS54a7, aiming to explore acceleration strategies of the MD simulations. The hydrogen mass repartitioning (HMR) and hydrogen isotope exchange (HIE) methods were tested to slow down the fastest degrees of freedom. These approaches allowed a significant increase in the time step employed in the equation of the motion integration algorithm, from 2 fs up to 5–7 fs, with no serious changes in the structural and dynamical properties of the nanopores. Subtle differences with respect to the simulations with the unmodified force fields were observed in the concerted movements of the cyclic peptides, as well as in the lifetime of several H-bonds. All together, these results are expected to contribute to better understanding of the behavior of self-assembled cyclic peptide nanotubes, as well as to support the methods tested to speed up general MD simulations; additionally, they do provide a number of quantitative descriptors that are expected to be used as a reference to design new experiments intended to validate and complement computational studies of antimicrobial cyclic peptides.

Cyclosporin A: Conformational Complexity and Chameleonicity

The chameleonic behavior of cyclosporin A (CsA) was investigated through conformational ensembles employing multicanonical molecular dynamics simulations that could sample the cis and trans isomers of N-methylated amino acids; these assessments were conducted in explicit water, dimethyl sulfoxide, acetonitrile, methanol, chloroform, cyclohexane (CHX), and n-hexane (HEX) using AMBER ff03, AMBER10:EHT, AMBER12:EHT, and AMBER14:EHT force fields. The conformational details were discussed employing the free-energy landscapes (FELs) at T = 300 K; it was observed that the experimentally determined structures of CsA were only a part of the conformational space. Comparing the ROESY measurements in CHX-d12 and HEX-d14, the major conformations in those apolar solvents were essentially the same as that in CDCl₃ except for the observation of some sidechain rotamers. The effects of the metal ions on the conformations, including the cis/trans isomerization, were also investigated. Based on the analysis of FELs, it was concluded that the AMBER ff03 force field best described the experimentally derived conformations, indicating that CsA intrinsically formed membrane-permeable conformations and that the metal ions might be the key to the cis/trans isomerization of N-methylated amino acids before binding a partner protein.

Quantum Chemistry Calculations for Metabolomics

A primary goal of metabolomics studies is to fully characterize the small-molecule composition of complex biological and environmental samples. However, despite advances in analytical technologies over the past two decades, the majority of small molecules in complex samples are not readily identifiable due to the immense structural and chemical diversity present within the metabolome. Current gold-standard identification methods rely on reference libraries built using authentic chemical materials ("standards"), which are not available for most molecules. Computational quantum chemistry methods, which can be used to calculate chemical properties that are then measured by analytical platforms, offer an alternative route for building reference libraries, i.e., in silico libraries for "standards-free" identification. In this review, we cover the major roadblocks currently facing metabolomics and discuss applications where quantum chemistry calculations offer a solution. Several successful examples for nuclear magnetic resonance spectroscopy, ion mobility spectrometry, infrared spectroscopy, and mass spectrometry methods are reviewed. Finally, we consider current best practices, sources of error, and provide an outlook for quantum chemistry calculations in metabolomics studies. We expect this review will inspire researchers in the field of small-molecule identification to accelerate adoption of in silico methods for generation of reference libraries and to add quantum chemistry calculations as another tool at their disposal to characterize complex samples.

Photoinduced reconfiguration to control the protein-binding affinity of azobenzene-cyclized peptides

The impact of next-generation biorecognition elements (ligands) will be determined by the ability to remotely control their binding activity for a target biomolecule in complex environments. Compared to conventional mechanisms for regulating binding affinity (pH, ionic strength, or chaotropic agents), light provides higher accuracy and rapidity, and is particularly suited for labile targets. In this study, we demonstrate a general method to develop azobenzene-cyclized peptide ligands with light-controlled affinity for target proteins. Light triggers a cis/trans isomerization of the azobenzene, which results in a major structural rearrangement of the cyclic peptide from a non-binding to a binding configuration. Critical to this goal are the ability to achieve efficient photo-isomerization under low light dosage and the temporal stability of both cis and trans isomers. We demonstrated our method by designing photo-switchable peptides targeting vascular cell adhesion marker 1 (VCAM1), a cell marker implicated in stem cell function. Starting from a known VCAM1-binding linear peptide, an ensemble of azobenzene-cyclized variants with selective light-controlled binding were identified by combining in silico design with experimental characterization via spectroscopy and surface plasmon resonance. Variant cycloAZOB[G-VHAKQHRN-K] featured rapid, light-controlled binding of VCAM1 (KD,trans/KD,cis ∼ 130). Biotin-cycloAZOB[G-VHAKQHRN-K] was utilized to label brain microvascular endothelial cells (BMECs), showing co-localization with anti-VCAM1 antibodies in cis configuration and negligible binding in trans configuration.

Elucidating Solution Structures of Cyclic Peptides Using Molecular Dynamics Simulations

Protein-protein interactions are vital to biological processes, but the shape and size of their interfaces make them hard to target using small molecules. Cyclic peptides have shown promise as protein-protein interaction modulators, as they can bind protein surfaces with high affinity and specificity. Dozens of cyclic peptides are already FDA approved, and many more are in various stages of development as immunosuppressants, antibiotics, antivirals, or anticancer drugs. However, most cyclic peptide drugs so far have been natural products or derivatives thereof, with de novo design having proven challenging. A key obstacle is structural characterization: cyclic peptides frequently adopt multiple conformations in solution, which are difficult to resolve using techniques like NMR spectroscopy. The lack of solution structural information prevents a thorough understanding of cyclic peptides' sequence-structure-function relationship. Here we review recent development and application of molecular dynamics simulations with enhanced sampling to studying the solution structures of cyclic peptides. We describe novel computational methods capable of sampling cyclic peptides' conformational space and provide examples of computational studies that relate peptides' sequence and structure to biological activity. We demonstrate that molecular dynamics simulations have grown from an explanatory technique to a full-fledged tool for systematic studies at the forefront of cyclic peptide therapeutic design.

Screening of Yeast Display Libraries of Enzymatically Treated Peptides to Discover Macrocyclic Peptide Ligands

We present the construction and screening of yeast display libraries of post-translationally modified peptides wherein site-selective enzymatic treatment of linear peptides is achieved using bacterial transglutaminase. To this end, we developed two alternative routes, namely (i) yeast display of linear peptides followed by treatment with recombinant transglutaminase in solution; or (ii) intracellular co-expression of linear peptides and transglutaminase to achieve peptide modification in the endoplasmic reticulum prior to yeast surface display. The efficiency of peptide modification was evaluated via orthogonal detection of epitope tags integrated in the yeast-displayed peptides by flow cytometry, and via comparative cleavage of putative cyclic vs. linear peptides by tobacco etch virus (TEV) protease. Subsequently, yeast display libraries of transglutaminase-treated peptides were screened to isolate binders to the N-terminal region of the Yes-Associated Protein (YAP) and its WW domains using magnetic selection and fluorescence activated cell sorting (FACS). The identified peptide cyclo[E-LYLAYPAH-K] featured a K_D of 1.75 μM for YAP and 0.68 μM for the WW domains of YAP as well as high binding selectivity against albumin and lysozyme. These results demonstrate the usefulness of enzyme-mediated cyclization in screening combinatorial libraries to identify cyclic peptide binders.

Discovery of Membrane-Permeating Cyclic Peptides via mRNA Display

Small synthetic peptides capable of crossing biological membranes represent valuable tools in cell biology and drug delivery. While several cell-penetrating peptides (CPPs) of natural or synthetic origin have been reported, no peptide is currently known to cross both cytoplasmic and outer embryonic membranes. Here, we describe a method to engineer membrane-permeating cyclic peptides (MPPs) with broad permeation activity by screening mRNA display libraries of cyclic peptides against embryos at different developmental stages. The proposed method was demonstrated by identifying peptides capable of permeating Drosophila melanogaster (fruit fly) embryos and mammalian cells. The selected peptide cyclo[Glut-MRKRHASRRE-K*] showed a strong permeation activity of embryos exposed to minimal permeabilization pretreatment, as well as human embryonic stem cells and a murine fibroblast cell line. Notably, in both embryos and mammalian cells, the cyclic peptide outperformed its linear counterpart and the control MPPs. Confocal microscopy and single cell flow cytometry analysis were utilized to assess the degree of permeation both qualitatively and quantitatively. These MPPs have potential application in studying and nondisruptively controlling intracellular or intraembryonic processes.

Affibody-Binding Ligands

While antibodies remain established therapeutic and diagnostic tools, other protein scaffolds are emerging as effective and safer alternatives. Affibodies in particular are a new class of small proteins marketed as bio-analytic reagents. They feature tailorable binding affinity, low immunogenicity, high tissue permeation, and high expression titer in bacterial hosts. This work presents the development of affibody-binding peptides to be utilized as ligands for their purification from bacterial lysates. Affibody-binding candidates were identified by screening a peptide library simultaneously against two model affibodies (anti-immunoglobulin G (IgG) and anti-albumin) with the aim of selecting peptides targeting the conserved domain of affibodies. An ensemble of homologous sequences identified from screening was synthesized on Toyopearl^® resin and evaluated via binding studies to select sequences that afford high product binding and recovery. The affibody-peptide interaction was also evaluated by in silico docking, which corroborated the targeting of the conserved domain. Ligand IGKQRI was validated through purification of an anti-ErbB2 affibody from an Escherichia coli lysate. The values of binding capacity (~5 mg affibody per mL of resin), affinity (K_D ~1 μM), recovery and purity (64-71% and 86-91%), and resin lifetime (100 cycles) demonstrate that IGKQRI can be employed as ligand in affibody purification processes.

Computational Methods for Studying Conformational Behaviors of Cyclic Peptides

Molecular dynamics (MD) simulations play more and more important roles in studying conformations of cyclic peptides in solution. Here we describe how to use replica-exchange molecular dynamics (REMD) implemented in Gromacs software package to simulate peptides with backbone cyclization and stapled peptides with side-chain linkages. Some of our methods for trajectory analyses and our residue-specific force fields are also described.

Docking Flexible Cyclic Peptides with AutoDock CrankPep

While a new therapeutic cyclic peptide is approved nearly every year, docking large macrocycles has remained challenging. Here, we present a new version of our peptide docking software AutoDock CrankPep (ADCP), extended to dock peptides cyclized through their backbone and/or side chain disulfide bonds. We show that within the top 10 solutions, ADCP identifies the proper interactions for 71% of a data set of 38 complexes, thus making it a useful tool for rational peptide-based drug design.

Applications of Molecular Dynamics Simulation in Structure Prediction of Peptides and Proteins

Compared with rapid accumulation of protein sequences from high-throughput DNA sequencing, obtaining experimental 3D structures of proteins is still much more difficult, making protein structure prediction (PSP) potentially very useful. Currently, a vast majority of PSP efforts are based on data mining of known sequences, structures and their relationships (informatics-based). However, if closely related template is not available, these methods are usually much less reliable than experiments. They may also be problematic in predicting the structures of naturally occurring or designed peptides. On the other hand, physics-based methods including molecular dynamics (MD) can utilize our understanding of detailed atomic interactions determining biomolecular structures. In this mini-review, we show that all-atom MD can predict structures of cyclic peptides and other peptide foldamers with accuracy similar to experiments. Then, some notable successes in reproducing experimental 3D structures of small proteins through MD simulations (some with replica-exchange) of the folding were summarized. We also describe advancements of MD-based refinement of structure models, and the integration of limited experimental or bioinformatics data into MD-based structure modeling.

Conformation and Permeability: Cyclic Hexapeptide Diastereomers

Conformational ensembles of eight cyclic hexapeptide diastereomers in explicit cyclohexane, chloroform, and water were analyzed by multicanonical molecular dynamics (McMD) simulations. Free-energy landscapes (FELs) for each compound and solvent were obtained from the molecular shapes and principal component analysis at T = 300 K; detailed analysis of the conformational ensembles and flexibility of the FELs revealed that permeable compounds have different structural profiles even for a single stereoisomeric change. The average solvent-accessible surface area (SASA) in cyclohexane showed excellent correlation with the cell permeability, whereas this correlation was weaker in chloroform. The average SASA in water correlated with the aqueous solubility. The average polar surface area did not correlate with cell permeability in these solvents. A possible strategy for designing permeable cyclic peptides from FELs obtained from McMD simulations is proposed.

β-Branched Amino Acids Stabilize Specific Conformations of Cyclic Hexapeptides

Cyclic peptides (CPs) are a promising class of molecules for drug development, particularly as inhibitors of protein-protein interactions. Predicting low-energy structures and global structural ensembles of individual CPs is critical for the design of bioactive molecules, but these are challenging to predict and difficult to verify experimentally. In our previous work, we used explicit-solvent molecular dynamics simulations with enhanced sampling methods to predict the global structural ensembles of cyclic hexapeptides containing different permutations of glycine, alanine, and valine. One peptide, cyclo-(VVGGVG) or P7, was predicted to be unusually well structured. In this work, we synthesized P7, along with a less well-structured control peptide, cyclo-(VVGVGG) or P6, and characterized their global structural ensembles in water using NMR spectroscopy. The NMR data revealed a structural ensemble similar to the prediction for P7 and showed that P6 was indeed much less well-structured than P7. We then simulated and experimentally characterized the global structural ensembles of several P7 analogs and discovered that β-branching at one critical position within P7 is important for overall structural stability. The simulations allowed deconvolution of thermodynamic factors that underlie this structural stabilization. Overall, the excellent correlation between simulation and experimental data indicates that our simulation platform will be a promising approach for designing well-structured CPs and also for understanding the complex interactions that control the conformations of constrained peptides and other macrocycles.

Understanding and designing head-to-tail cyclic peptides

Cyclic peptides (CPs) are an exciting class of molecules with a variety of applications. However, design strategies for CP therapeutics, for example, are generally limited by a poor understanding of their sequence-structure relationships. This knowledge gap often leads to a trial-and-error approach for designing CPs for a specific purpose, which is both costly and time-consuming. Herein, we describe the current experimental and computational efforts in understanding and designing head-to-tail CPs along with their respective challenges. In addition, we provide several future directions in the field of computational CP design to improve its accuracy, efficiency and applicability. These advances, combined with experimental techniques, shall ultimately provide a better understanding of these interesting molecules and a reliable working platform to rationally design CPs with desired characteristics.

A critical assessment of force field accuracy against NMR data for cyclic peptides containing β-amino acids

Hybrid cyclic α/β-peptides, in which one or more β-amino acids are incorporated into the backbone, are gaining increasing interest as potential therapeutics, thanks to their ability to achieve enhanced binding affinities for a biological target through pre-organization in solution. The in silico prediction of their three dimensional structure through strategies such as MD simulations would substantially advance the rational design process. However, whether the molecular mechanics force fields are accurate in sampling highly constrained cyclopeptides containing β-amino acids remains to be verified. Here, we present a systematic assessment of the ability of 8 widely used force fields to reproduce 79 NMR observables (including chemical shifts and 3J scalar couplings) on five cyclic α/β-peptides that contain the integrin recognition motif isoDGR. Most of the investigated force fields, which include force fields from AMBER, OPLS, CHARMM and GROMOS families, display very good agreement with experimental 3J(HN,Hα), suggesting that MD simulations could be an appropriate tool in the rational design of therapeutic cyclic α-peptides. However, for NMR observables directly related to β-amino acids, we observed a poor agreement with experiments and a remarkable dependence of our evaluation on the choice of Karplus parameters. The force field weaknesses herein unveiled might constitute a source of inspiration for further force field optimization.

Succinimide-Based Conjugates Improve IsoDGR Cyclopeptide Affinity to αvβ3 without Promoting Integrin Allosteric Activation

The isoDGR sequence is an integrin-binding motif that has been successfully employed as a tumor-vasculature-homing molecule or for the targeted delivery of drugs and diagnostic agents to tumors. In this context, we previously demonstrated that cyclopeptide 2, the product of the conjugation of c(CGisoDGRG) (1) to 4-( N-maleimidomethyl)cyclohexane-1-carboxamide, can be successfully used as a tumor-homing ligand for nanodrug delivery to neoplastic tissues. Here, combining NMR, computational, and biochemical methods, we show that the succinimide ring contained in 2 contributes to stabilizing interactions with α_vβ₃, an integrin overexpressed in the tumor vasculature. Furthermore, we demonstrate that various cyclopeptides containing the isoDGR sequence embedded in different molecular scaffolds do not induce α_vβ₃ allosteric activation and work as pure integrin antagonists. These results could be profitably exploited for the rational design of novel isoDGR-based ligands and tumor-targeting molecules with improved α_vβ₃-binding properties and devoid of adverse integrin-activating effects.

Designing Well-Structured Cyclic Pentapeptides Based on Sequence-Structure Relationships

Cyclic peptides are a promising class of molecules for unique applications. Unfortunately, cyclic peptide design is severely limited by the difficulty in predicting the conformations they will adopt in solution. In this work, we use explicit-solvent molecular dynamics simulations to design well-structured cyclic peptides by studying their sequence-structure relationships. Critical to our approach is an enhanced sampling method that exploits the essential transitional motions of cyclic peptides to efficiently sample their conformational space. We simulated a range of cyclic pentapeptides from all-glycine to a library of cyclo-(X1X2AAA) peptides to map their conformational space and determine cooperative effects of neighboring residues. By combining the results from all cyclo-(X1X2AAA) peptides, we developed a scoring function to predict the structural preferences for X1-X2 residues within cyclic pentapeptides. Using this scoring function, we designed a cyclic pentapeptide, cyclo-(GNSRV), predicted to be well structured in aqueous solution. Subsequent circular dichroism and NMR spectroscopy revealed that this cyclic pentapeptide is indeed well structured in water, with a nuclear Overhauser effect and J-coupling values consistent with the predicted structure.

Mapping the sequence-structure relationships of simple cyclic hexapeptides

Cyclic peptides are promising protein-protein interaction modulators with high binding affinities and specificities, as well as enhanced stabilities and oral availabilities over linear analogs. Despite their relatively small size and cyclic architecture, it is currently difficult to predict the favored conformation(s) of most classes of cyclic peptides. An improved understanding of the sequence-structure relationships for cyclic peptides will offer an avenue for the rational design of cyclic peptides as possible therapeutics. In this work, we systematically explored the sequence-structure relationships for two cyclic hexapeptide systems using molecular dynamics simulation techniques. Starting with an all-glycine cyclic hexapeptide, cyclo-G₆, we systematically replaced glycine residues with alanines and characterized the structural ensembles of different variants. The same process was repeated with valines to investigate the effects of larger side chains. An analysis of the origin of structure preferences was performed using thermodynamics decomposition and several general observations are reported.

Toward accurately modeling N-methylated cyclic peptides

Cyclic peptides have unique properties and can target protein surfaces specifically and potently. N-Methylation provides a promising way to further optimize the pharmacokinetic and structural profiles of cyclic peptides. The capability to accurately model structures adopted by N-methylated cyclic peptides would facilitate rational design of this interesting and useful class of molecules. We apply molecular dynamics simulations with advanced enhanced sampling methods to efficiently characterize the structural ensembles of N-methylated cyclic peptides, while simultaneously evaluating the overall performance of several simulation force fields. We find that one of the residue-specific force fields, RSFF2, is able to recapitulate experimental structures of the N-methylated cyclic peptide benchmarks tested here when the correct amide isomers are used as initial configurations and enforced during the simulations. Thus, using our simulation approach, it is possible to accurately and efficiently predict the structures of N-methylated cyclic peptides if sufficient information is available to determine the correct amide cis/trans configuration. Moreover, our results suggest that, upon further optimization of RSFF2 to more reliably predict cis/trans isomers, molecular dynamics simulations will be able to de novo predict N-methylated cyclic peptides in the near future, strongly motivating such continued optimization.

Analysis of Loops that Mediate Protein-Protein Interactions and Translation into Submicromolar Inhibitors

Effective strategies for mimicking α-helix and β-strand epitopes have been developed, producing valuable inhibitors for some classes of protein-protein interactions (PPIs). However, there are no general strategies for translating loop epitopes into useful PPI inhibitors. In this work, we use the LoopFinder program to identify diverse sets of "hot loops," which are loop epitopes that mediate PPIs. These include loops that are well-suited to mimicry with macrocyclic compounds, and loops that are most similar to variable loops on antibodies and ankyrin repeat proteins. We present data-driven criteria for scoring loop-mediated PPIs, uncovering a trove of potentially druggable interactions. We also use unbiased clustering to identify common structures among the hot loops. To translate these insights into real-world inhibitors, we describe a robust, diversity-oriented strategy for the rapid production and evaluation of cyclized loops. This method is applied to a computationally identified loop in the PPI between stonin2 and Eps15, producing submicromolar inhibitors. The most potent inhibitor is well-structured in water and successfully mimics the native epitope. Overall, these computational and experimental strategies provide new opportunities to design inhibitors for an otherwise intractable set of PPIs.

Computational methods to design cyclic peptides

Cyclic peptides (CPs) are promising modulators of protein-protein interactions (PPIs), but their application remains challenging. It is currently difficult to predict the structures and bioavailability of CPs. The ability to design CPs using computer modeling would greatly facilitate the development of CPs as potent PPI modulators for fundamental studies and as potential therapeutics. Herein, we describe computational methods to generate CP libraries for virtual screening, as well as current efforts to accurately predict the conformations adopted by CPs. These advances are making it possible to envision robust computational design of active CPs. However, unique properties of CPs pose significant challenges associated with sampling CP conformational space and accurately describing CP energetics. These major obstacles to structure prediction likely must be solved before robust design of active CPs can be reliably achieved.

Conformational Restriction of Peptides Using Dithiol Bis-Alkylation

Macrocyclic peptides are highly promising as inhibitors of protein-protein interactions. While many bond-forming reactions can be used to make cyclic peptides, most have limitations that make this chemical space challenging to access. Recently, a variety of cysteine alkylation reactions have been used in rational design and library approaches for cyclic peptide discovery and development. We and others have found that this chemistry is versatile and robust enough to produce a large variety of conformationally constrained cyclic peptides. In this chapter, we describe applications, methods, mechanistic insights, and troubleshooting for dithiol bis-alkylation reactions for the production of cyclic peptides. This method for efficient solution-phase macrocyclization is highly useful for the rapid production and screening of loop-based inhibitors of protein-protein interactions.

Accurate Structure Prediction and Conformational Analysis of Cyclic Peptides with Residue-Specific Force Fields

Cyclic peptides (CPs) are promising candidates for drugs, chemical biology tools, and self-assembling nanomaterials. However, the development of reliable and accurate computational methods for their structure prediction has been challenging. Here, 20 all-trans CPs of 5-12 residues selected from Cambridge Structure Database have been simulated using replica-exchange molecular dynamics with four different force fields. Our recently developed residue-specific force fields RSFF1 and RSFF2 can correctly identify the crystal-like conformations of more than half CPs as the most populated conformation. The RSFF2 performs the best, which consistently predicts the crystal structures of 17 out of 20 CPs with rmsd < 1.1 Å. We also compared the backbone (ϕ, ψ) sampling of residues in CPs with those in short linear peptides and in globular proteins. In general, unlike linear peptides, CPs have local conformational free energies and entropies quite similar to globular proteins.

Reference spectral signature selection using density-based cluster for automatic oil spill detection in hyperspectral images

Reference spectral signature selection is a fundamental work for automatic oil spill detection. To address this issue, a new approach is proposed here, which employs the density-based cluster to select a specific spectral signature from a hyperspectral image. This paper first introduces the framework of oil spill detection from hyperspectral images, indicating that detecting oil spill requires a reference spectral signature of oil spill, parameters of background, and a target detection algorithm. Based on the framework, we give the new reference spectral signature selection approach in details. Then, we demonstrate the estimation of background parameters according to the reflectance of seawater in the infrared bands. Next, the conventional adaptive cosine estimator (ACE) algorithm is employed to achieve oil spill detection. Finally, the proposed approach is tested via several practical hyperspectral images that are collected during the Horizon Deep water oil spill. The experimental results show that this new approach can automatically select the reference spectral signature of oil spills from hyperspectral images and has high detection performance.

Going Out on a Limb: Delineating The Effects of β-Branching, N-Methylation, and Side Chain Size on the Passive Permeability, Solubility, and Flexibility of Sanguinamide A Analogues

It is well established that intramolecular hydrogen bonding and N-methylation play important roles in the passive permeability of cyclic peptides, but other structural features have been explored less intensively. Recent studies on the oral bioavailability of the cyclic heptapeptide sanguinamide A have raised the question of whether steric occlusion of polar groups via β-branching is an effective, yet untapped, tool in cyclic peptide permeability optimization. We report the structures of 17 sanguinamide A analogues designed to test the relative contributions of β-branching, N-methylation, and side chain size to passive membrane permeability and aqueous solubility. We demonstrate that β-branching has little effect on permeability compared to the effects of aliphatic carbon count and N-methylation of exposed NH groups. We highlight a new N-methylated analogue of sanguinamide A with a Leu substitution at position 2 that exhibits solvent-dependent flexibility and improved permeability over that of the natural product.

Molecular Simulation of Conformational Pre-Organization in Cyclic RGD Peptides

To test the ability of molecular simulations to accurately predict the solution-state conformational properties of peptidomimetics, we examined a test set of 18 cyclic RGD peptides selected from the literature, including the anticancer drug candidate cilengitide, whose favorable binding affinity to integrin has been ascribed to its pre-organization in solution. For each design, we performed all-atom replica-exchange molecular dynamics simulations over several microseconds and compared the results to extensive published NMR data. We find excellent agreement with experimental NOE distance restraints, suggesting that molecular simulation can be a useful tool for the computational design of pre-organized solution-state structure. Moreover, our analysis of conformational populations estimates that, despite the potential for increased flexibility due to backbone amide isomerizaton, N-methylation provides about 0.5 kcal/mol of reduced conformational entropy to cyclic RGD peptides. The combination of pre-organization and binding-site compatibility explains the strong binding affinity of cilengitide to integrin.