Probing the Origin of Structural Stability of Single and Double Stapled p53 Peptide Analogs Bound to MDM2

Dissecting the geometric and hydrophobic constraints of stapled peptides

Stapled peptides are a promising class of molecules with potential as highly specific probes of protein-protein interactions and as therapeutics. Hydrocarbon stapling affects the peptide properties through the interplay of two factors: enhancing the overall hydrophobicity and constraining the conformational flexibility. By constructing a series of virtual peptides, we study the role of each factor in modulating the structural properties of a hydrocarbon-stapled peptide PM2, which has been shown to enter cells, engage its target Mouse Double Minute 2 (MDM2), and activate p53. Hamiltonian replica exchange molecular dynamics (HREMD) simulations suggest that hydrocarbon stapling favors helical populations of PM2 through a combination of the geometric constraints and the enhanced hydrophobicity of the peptide. To further understand the conformational landscape of the stapled peptides along the binding pathway, we performed HREMD simulations by restraining the peptide at different distances from MDM2. When the peptide approaches MDM2, the binding pocket undergoes dehydration which appears to be greater in the presence of the stapled peptide compared with the linear peptide. In the binding pocket, the helicity of the stapled peptide is increased due to the favorable interactions between the peptide residues as well as the staple and the microenvironment of the binding pocket, contributing to enhanced affinity. The dissection of the multifaceted mechanism of hydrocarbon stapling into individual factors not only deepens fundamental understanding of peptide stapling, but also provides guidelines for the design of new stapled peptides.

Solution-State Preorganization of Cyclic β-Hairpin Ligands Determines Binding Mechanism and Affinities for MDM2

Understanding mechanisms of protein folding and binding is crucial to designing their molecular function. Molecular dynamics (MD) simulations and Markov state model (MSM) approaches provide a powerful way to understand complex conformational change that occurs over long time scales. Such dynamics are important for the design of therapeutic peptidomimetic ligands, whose affinity and binding mechanism are dictated by a combination of folding and binding. To examine the role of preorganization in peptide binding to protein targets, we performed massively parallel explicit-solvent MD simulations of cyclic β-hairpin ligands designed to mimic the p53 transactivation domain and competitively bind mouse double minute 2 homologue (MDM2). Disrupting the MDM2-p53 interaction is a therapeutic strategy to prevent degradation of the p53 tumor suppressor in cancer cells. MSM analysis of over 3 ms of aggregate trajectory data enabled us to build a detailed mechanistic model of coupled folding and binding of four cyclic peptides which we compare to experimental binding affinities and rates. The results show a striking relationship between the relative preorganization of each ligand in solution and its affinity for MDM2. Specifically, changes in peptide conformational populations predicted by the MSMs suggest that entropy loss upon binding is the main factor influencing affinity. The MSMs also enable detailed examination of non-native interactions which lead to misfolded states and comparison of structural ensembles with experimental NMR measurements. In contrast to an MSM study of p53 transactivation domain (TAD) binding to MDM2, MSMs of cyclic β-hairpin binding show a conformational selection mechanism. Finally, we make progress toward predicting accurate off rates of cyclic peptides using multiensemble Markov models (MEMMs) constructed from unbiased and biased simulated trajectories.

Macrocyclic α helical peptide therapeutic modality: A perspective of learnings and challenges

Macrocyclic α-helical peptides have emerged as a compelling new therapeutic modality to tackle targets confined to the intracellular compartment. Within the scope of hydrocarbon-stapling there has been significant progress to date, including the first stapled α-helical peptide to enter into clinical trials. The principal design concept of stapled α-helical peptides is to mimic a cognate (protein) ligand relative to binding its target via an α-helical interface. However, it was the proclivity of such stapled α-helical peptides to exhibit cell permeability and proteolytic stability that underscored their promise as unique macrocyclic peptide drugs for intracellular targets. This perspective highlights key learnings as well as challenges in basic research with respect to structure-based design, innovative chemistry, cell permeability and proteolytic stability that are essential to fulfill the promise of stapled α-helical peptide drug development.

Molecular Simulations Identify Binding Poses and Approximate Affinities of Stapled α-Helical Peptides to MDM2 and MDMX

Traditionally, computing the binding affinities of proteins to even relatively small and rigid ligands by free-energy methods has been challenging due to large computational costs and significant errors. Here, we apply a new molecular simulation acceleration method called MELD (Modeling by Employing Limited Data) to study the binding of stapled α-helical peptides to the MDM2 and MDMX proteins. We employ free-energy-based molecular dynamics simulations (MELD-MD) to identify binding poses and calculate binding affinities. Even though stapled peptides are larger and more complex than most protein ligands, the MELD-MD simulations can identify relevant binding poses and compute relative binding affinities. MELD-MD appears to be a promising method for computing the binding properties of peptide ligands with proteins.

Investigating Protein-Peptide Interactions Using the Schrödinger Computational Suite

The Schrödinger software suite contains a broad array of computational chemistry and molecular modeling tools that can be used to study the interaction of peptides with proteins. These include molecular docking using Glide and Piper, relative binding free energy predictions with FEP+, conformational searches using MacroModel and Desmond, and structural refinement using Prime and PrimeX. In this review we provide a comprehensive overview of these tools and describe their potential application in the identification and optimization of peptide ligands for proteins.

Markov models of the apo-MDM2 lid region reveal diffuse yet two-state binding dynamics and receptor poses for computational docking

MDM2 is a negative regulator of p53 activity and an important target for cancer therapeutics. The N-terminal lid region of MDM2 modulates interactions with p53 via competition for its binding cleft, exchanging slowly between docked and undocked conformations in the absence of p53. To better understand these dynamics, we constructed Markov State Models (MSMs) from large collections of unbiased simulation trajectories of apo-MDM2, and find strong evidence for diffuse, yet two-state folding and binding of the N-terminal region to the p53 receptor site. The MSM also identifies holo-like receptor conformations highly suitable for computational docking, despite initiating trajectories from closed-cleft receptor structures unsuitable for docking. Fixed-anchor docking studies using a test set of high-affinity small molecules and peptides show simulated receptor ensembles achieve docking successes comparable to cross-docking studies using crystal structures of receptors bound by alternative ligands. For p53, the best-scoring receptor structures have the N-terminal region lid region bound in a helical conformation mimicking the bound structure of p53, suggesting lid region association induces receptor conformations suitable for binding. These results suggest that MD + MSM approaches can sample binding-competent receptor conformations suitable for computational peptidomimetic design, and that inclusion of disordered regions may be essential to capturing the correct receptor dynamics.

Bridged Analogues for p53-Dependent Cancer Therapy Obtained by S-Alkylation

A small library of anticancer, cell-permeating, stapled peptides based on potent dual-specific antagonist of p53-MDM2/MDMX interactions, PMI-N8A, was synthesized, characterized and screened for anticancer activity against human colorectal cancer cell line, HCT-116. Employed synthetic modifications included: S-alkylation-based stapling, point mutations increasing hydrophobicity in key residues as well as improvement of cell-permeability by introduction of polycationic sequence(s) that were woven into the sequence of parental peptide. Selected analogue, ArB14Co, was also tested in vivo and exhibited potent anticancer bioactivity at the low dose (3.0 mg/kg). Collectively, our findings suggest that application of stapling in combination with rational design of polycationic short analogues may be a suitable approach in the development of physiologically active p53-MDM2/MDMX peptide inhibitors.

Hydrocarbon stapled peptides as modulators of biological function

Peptide-based drug discovery has experienced a significant upturn within the past decade since the introduction of chemical modifications and unnatural amino acids has allowed for overcoming some of the drawbacks associated with peptide therapeutics. Strengthened by such features, modified peptides become capable of occupying a niche that emerges between the two major classes of today's therapeutics-small molecules (<500 Da) and biologics (>5000 Da). Stabilized α-helices have proven particularly successful at impairing disease-relevant PPIs previously considered "undruggable." Among those, hydrocarbon stapled α-helical peptides have emerged as a novel class of potential peptide therapeutics. This review provides a comprehensive overview of the development and applications of hydrocarbon stapled peptides discussing the benefits and limitations of this technique.

Macrocyclic α-Helical Peptide Drug Discovery

Macrocyclic α-helical peptides have emerged as a promising new drug class and within the scope of hydrocarbon-stapled peptides such molecules have advanced into the clinic. The overarching concept of designing proteomimetics of an α-helical ‘ligand’ which binds its cognate ‘target’ relative to α-helical interfacing protein-protein interactions has been well-validated and expanded through numerous investigations for a plethora of therapeutic targets oftentimes referred to as “undruggable” with respect to other modalities (e.g., small-molecule or proteins). This chapter highlights the evolution of macrocyclic α-helical peptides in terms of target space, biophysical and computational chemistry, structural diversity and synthesis, drug design and chemical biology. It is noteworthy that hydrocarbon-stapled peptides have successfully risen to the summit of such drug discovery campaigns.

Heterogeneous Hydration of p53/MDM2 Complex

Water-mediated interactions play critical roles in biomolecular recognition processes. Explicit solvent molecular dynamics (MD) simulations and the variational implicit-solvent model (VISM) are used to study those hydration properties during binding for the biologically important p53/MDM2 complex. Unlike simple model solutes, in such a realistic and heterogeneous solute–solvent system with both geometrical and chemical complexity, the local water distribution sensitively depends on nearby amino acid properties and the geometric shape of the protein. We show that the VISM can accurately describe the locations of high and low density solvation shells identified by the MD simulations and can explain them by a local coupling balance of solvent–solute interaction potentials and curvature. In particular, capillary transitions between local dry and wet hydration states in the binding pocket are captured for interdomain distance between 4 to 6 Å, right at the onset of binding. The underlying physical connection between geometry and polarity is illustrated and quantified. Our study offers a microscopic and physical insight into the heterogeneous hydration behavior of the biologically highly relevant p53/MDM2 system and demonstrates the fundamental importance of hydrophobic effects for biological binding processes. We hope our study can help to establish new design rules for drugs and medical substances.