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.

A backbone-dependent rotamer library with high (ϕ, ψ) coverage using metadynamics simulations.

Backbone-dependent rotamer libraries are commonly used to assign the side chain dihedral angles of amino acids when modeling protein structures. Most rotamer libraries are created by curating protein crystal structure data and using various methods to extrapolate the existing data to cover all possible backbone conformations. However, these rotamer libraries may not be suitable for modeling the structures of cyclic peptides and other constrained peptides because these molecules frequently sample backbone conformations rarely seen in the crystal structures of linear proteins. To provide backbone-dependent side chain information beyond the α-helix, ß-sheet, and PPII regions, we used explicit-solvent metadynamics simulations of model dipeptides to create a new rotamer library that has high coverage in the (ϕ, ψ) space. Furthermore, this approach can be applied to build high-coverage rotamer libraries for noncanonical amino acids. The resulting Metadynamics of Dipeptides for Rotamer Distribution (MEDFORD) rotamer library predicts the side chain conformations of high-resolution protein crystal structures with similar accuracy (~80%) to a state-of-the-art rotamer library. Our ability to test the accuracy of MEDFORD at predicting the side chain dihedral angles of amino acids in noncanonical backbone conformation is restricted by the limited structural data available for cyclic peptides. For the cyclic peptide data that are currently available, MEDFORD and the state-of-the-art rotamer library perform comparably. However, the two rotamer libraries indeed make different rotamer predictions in noncanonical (ϕ, ψ) regions. For noncanonical amino acids, the MEDFORD rotamer library predicts the χ1 values with approximately 75% accuracy.

CATBOSS: Cluster Analysis of Trajectories Based on Segment Splitting.

Molecular dynamics (MD) simulations are an exceedingly and increasingly potent tool for molecular behavior prediction and analysis. However, the enormous wealth of data generated by these simulations can be difficult to process and render in a human-readable fashion. Cluster analysis is a commonly used way to partition data into structurally distinct states. We present a method that improves on the state of the art by taking advantage of the temporal information of MD trajectories to enable more accurate clustering at a lower memory cost. To date, cluster analysis of MD simulations has generally treated simulation snapshots as a mere collection of independent data points and attempted to separate them into different clusters based on structural similarity. This new method, cluster analysis of trajectories based on segment splitting (CATBOSS), applies density-peak-based clustering to classify trajectory segments learned by change detection. Applying the method to a synthetic toy model as well as four real-life data sets-trajectories of MD simulations of alanine dipeptide and valine dipeptide as well as two fast-folding proteins-we find CATBOSS to be robust and highly performant, yielding natural-looking cluster boundaries and greatly improving clustering resolution. As the classification of points into segments emphasizes density gaps in the data by grouping them close to the state means, CATBOSS applied to the valine dipeptide system is even able to account for a degree of freedom deliberately omitted from the input data set. We also demonstrate the potential utility of CATBOSS in distinguishing metastable states from transition segments as well as promising application to cases where there is little or no advance knowledge of intrinsic coordinates, making for a highly versatile analysis tool.

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.

Cyclic peptides: backbone rigidification and capability of mimicking motifs at protein-protein interfaces.

Cyclization is commonly employed in efforts to improve the target binding affinity of peptide-based probes and therapeutics. Many structural motifs have been identified at protein-protein interfaces and provide promising targets for inhibitor design using cyclic peptides. Cyclized peptides are generally assumed to be rigidified relative to their linear counterparts. This rigidification potentially pre-organizes the molecules to interact properly with their targets. However, the actual impact of cyclization on, for example, peptide configurational entropy, is currently poorly understood in terms of both its magnitude and molecular-level origins. Moreover, even with thousands of desired structural motifs at hand, it is currently not possible to a priori identify the ones that are most promising to mimic using cyclic peptides nor to select the ideal linker length. Instead, labor-intensive chemical synthesis and experimental characterization of various cyclic peptide designs are required, in hopes of finding one with improved target affinity. Herein, using molecular dynamics simulations of polyglycines, we elucidated how head-to-tail cyclization impacts peptide backbone dihedral entropy and developed a simple strategy to rapidly screen for structures that can be reliably mimicked by preorganized cyclic peptides. As expected, cyclization generally led to a reduction in backbone dihedral entropy; notably, however, this effect was minimal when the length of polyglycines was >9 residues. We also found that the reduction in backbone dihedral entropy upon cyclization of small polyglycine peptides does not result from more restricted distributions of the dihedrals; rather, it was the correlations between specific dihedrals that caused the decrease in configurational entropy in the cyclic peptides. Using our comprehensive cyclo-Gn structural ensembles, we obtained a holistic picture of what conformations are accessible to cyclic peptides. Using "hot loops" recently identified at protein-protein interfaces as an example, we provide clear guidelines for choosing the "easiest" hot loops for cyclic peptides to mimic and for identifying appropriate cyclic peptide lengths. In conclusion, our results provide an understanding of the thermodynamics and structures of this interesting class of molecules. This information should prove particularly useful for designing cyclic peptide inhibitors of protein-protein interactions.