Use of 1,2,3-trisubstituted cyclopropanes as conformationally constrained peptide mimics in SH2 antagonists

Insights into Modeling Approaches in Chemistry: Assessing Ligand-Protein Binding Thermodynamics Based on Rigid-Flexible Model Molecules

In the field of chemistry, model compounds find extensive use for investigating complex objects. One prime example of such object is the protein-ligand supramolecular interaction. Prediction the enthalpic and entropic contribution to the free energy associated with this process, as well as the structural and dynamic characteristics of protein-ligand complexes poses considerable challenges. This review exemplifies modeling approaches used to study protein-ligand binding (PLB) thermodynamics by employing pairs of conformationally constrained/flexible model molecules. Strategically designing the model molecules can reduce the number of variables that influence thermodynamic parameters. This enables scientists to gain deeper insights into the enthalpy and entropy of PLB, which is relevant for medicinal chemistry and drug design. The model studies reviewed here demonstrate that rigidifying ligands may induce compensating changes in the enthalpy and entropy of binding. Some "rules of thumb" have started to emerge on how to minimize entropy-enthalpy compensation and design efficient rigidified or flexible ligands.

Natural Products and Their Mimics as Targets of Opportunity for Discovery

Diverse structural types of natural products and their mimics have served as targets of opportunity in our laboratory to inspire the discovery and development of new methods and strategies to assemble polyfunctional and polycyclic molecular architectures. Furthermore, our efforts toward identifying novel compounds having useful biological properties led to the creation of new targets, many of which posed synthetic challenges that required the invention of new methodology. In this Perspective, selected examples of how we have exploited a diverse range of natural products and their mimics to create, explore, and solve a variety of problems in chemistry and biology will be discussed. The journey was not without its twists and turns, but the unexpected often led to new revelations and insights. Indeed, in our recent excursion into applications of synthetic organic chemistry to neuroscience, avoiding the more-traveled paths was richly rewarding.

Relative Binding Enthalpies from Molecular Dynamics Simulations Using a Direct Method

The potential for reliably predicting relative binding enthalpies, ΔΔE, from a direct method utilizing molecular dynamics is examined for a system of three phosphotyrosyl peptides binding to a protein receptor, the Src SH2 domain. The binding enthalpies were calculated from the potential energy differences between the bound and the unbound end-states of each peptide from equilibrium simulations in explicit water. The statistical uncertainties in the ensemble-mean energy values from multiple, independent simulations were obtained using a bootstrap method. Simulations were initiated with different starting coordinates as well as different velocities. Statistical uncertainties in ΔΔE are 2 to 3 kcal/mol based on calculations from 40, 10 ns trajectories for each system (three SH2–peptide complexes or unbound peptides). Uncertainties in relative component energies, comprising solute–solute, solute–solvent and solvent–solvent interactions, are considerably larger. Energy values were estimated from an unweighted ensemble averaging of multiple trajectories with the a priori assumption that all trajectories are equally likely. Distributions in energy–rmsd space indicate that the trajectories sample the same basin and the difference in mean energy values between trajectories is due to sampling of alternative local regions of this superbasin. The direct estimate of relative binding enthalpies is concluded to be a reasonable approach for well-ordered systems with ΔΔE values greater than ∼3 kcal/mol, although the approach would benefit from future work to determine properly distributed starting points that would enable efficient sampling of conformational space using multiple trajectories.

Correlating structure and energetics in protein-ligand interactions: paradigms and paradoxes

Predicting protein-binding affinities of small molecules, even closely related ones, is a formidable challenge in biomolecular recognition and medicinal chemistry. A thermodynamic approach to optimizing affinity in protein-ligand interactions requires knowledge and understanding of how altering the structure of a small molecule will be manifested in protein-binding enthalpy and entropy changes; however, there is a relative paucity of such detailed information. In this review, we examine two strategies commonly used to increase ligand potency. The first of these involves introducing a cyclic constraint to preorganize a small molecule in its biologically active conformation, and the second entails adding nonpolar groups to a molecule to increase the amount of hydrophobic surface that is buried upon binding. Both of these approaches are motivated by paradigms suggesting that protein-binding entropy changes should become more favorable, but paradoxes can emerge that defy conventional wisdom.

Probing the effect of conformational constraint on phosphorylated ligand binding to an SH2 domain using polarizable force field simulations

Preorganizing a ligand in the conformation it adopts upon binding to a protein has long been considered to be an effective way to improve affinity by making the binding entropy more favorable. However, recent thermodynamic studies of a series of complexes of the Grb2 SH2 domain with peptide analogues having constrained and flexible replacements for a phosphotyrosine residue revealed that less favorable binding entropies may result from constraining ligands in their biologically active conformations. Toward probing the origin of this unexpected finding, we examined the complexes of four phosphotyrosine-derived analogues with the Grb2 SH2 domain using molecular dynamics simulations with a polarizable force field. Significantly, the computed values for the relative binding free energies, entropies, and enthalpies of two pairs of constrained and unconstrained ligands reproduced the trends that were determined experimentally, although the relative differences were overestimated. These calculations also revealed that a large fraction of the ligands lacking the constraining element exist in solution as compact, macrocyclic-like structures that are stabilized by interactions between the phosphate groups and the amide moieties of the C-terminal pY+2 residues. In contrast, the three-membered ring in the constrained ligands prevents the formation of such macrocyclic structures, leading instead to globally extended, less ordered conformations. Quasiharmonic analysis of these conformational ensembles suggests that the unconstrained ligands possess significantly lower entropies in solution, a finding that is consistent with the experimental observation that the binding entropies for the unconstrained ligands are more favorable than for their constrained counterparts. This study suggests that introducing local constraints in flexible molecules may have unexpected consequences, and a detailed understanding of the conformational preferences of ligands in their unbound states is a critical prerequisite to correlating changes in their chemical structure with protein binding entropies and enthalpies.

Constraining binding hot spots: NMR and molecular dynamics simulations provide a structural explanation for enthalpy-entropy compensation in SH2-ligand binding

NMR spectroscopy and molecular dynamics (MD) simulations were used to probe the structure and dynamics of complexes of three phosphotyrosine-derived peptides with the Src SH2 domain in an effort to uncover a structural explanation for enthalpy-entropy compensation observed in the binding thermodynamics. The series of phosphotyrosine peptide derivatives comprises the natural pYEEI Src SH2 ligand, a constrained mimic, in which the phosphotyrosine (pY) residue is preorganized in the bound conformation for the purpose of gaining an entropic advantage to binding, and a flexible analogue of the constrained mimic. The expected gain in binding entropy of the constrained mimic was realized; however, a balancing loss in binding enthalpy was also observed that could not be rationalized from the crystallographic structures. We examined protein dynamics to evaluate whether the observed enthalpic penalty might be the result of effects arising from altered motions in the complex. (15)N-relaxation studies and positional fluctuations from molecular dynamics indicate that the main-chain dynamics of the protein show little variation among the three complexes. Root mean squared (rms) coordinate deviations vary by less than 1.5 A for all non-hydrogen atoms for the crystal structures and in the ensemble average structures calculated from the simulations. In contrast to this striking similarity in the structures and dynamics, there are a number of large chemical shift differences from residues across the binding interface, but particularly from key Src SH2 residues that interact with pY, the "hot spot" residue, which contributes about one-half of the binding free energy. Rank-order correlations between chemical shifts and ligand binding enthalpy for several pY-binding residues, coupled with available mutagenesis and calorimetric data, suggest that subtle structural perturbations (<1 A) from the conformational constraint of the pY residue sufficiently alter the geometry of enthalpically critical interactions in the binding pocket to cause the loss of binding enthalpy, leading to the observed enthalpy-entropy compensation. We find no evidence to support the premise that enthalpy-entropy compensation is an inherent property and conclude that preorganization of Src SH2 ligand residues involved in binding hot spots may eventuate in suboptimal interactions with the domain. We propose that introducing constraints elsewhere in the ligand could minimize enthalpy-entropy compensation effects. The results illustrate the utility of the NMR chemical shift to highlight small, but energetically significant, perturbations in structure that might otherwise go unnoticed in an apparently rigid protein.

Preorganization in biological systems: Are conformational constraints worth the energy?

It is generally assumed that preorganizing a flexible ligand in the three-dimensional shape it adopts when bound to a macromolecular receptor will provide a derivative having an increased binding affinity, primarily because the rigidified molecule is expected to benefit from a lesser entropic penalty during complexation. We now provide the first experimental evidence that demonstrates this common belief is not universally true. Indeed, we find that ligand preorganization may be accompanied by an unfavorable entropy of binding, even when the constrained ligand exhibits a higher binding affinity than its flexible control. Thus, the effects that ligand preorganization have upon energetics and structure in protein-ligand interactions must be reevaluated.

Synthesis and properties of cyclopropane-derived peptidomimetics

Small peptides exhibit a wide range of biological activities, but although there are some notable exceptions, they are not generally useful as drugs. This has spurred widespread interest in designing peptidomimetics and introducing them as replacements of portions of native peptides to enhance their biological properties. Special attention has been focused upon rigid replacements because of their potential to preorganize the resulting pseudopeptide in a conformation corresponding to its bound structure. Toward this goal, we invented trisubstituted cyclopropanes as novel peptidomimetics, anticipating that the cyclopropane ring would locally orient the backbone and the corresponding amino acid side chain in the biologically active conformation. Selected aspects of the syntheses and applications of these cyclopropane-derived peptidomimetics are presented in this Account.

Phosphotyrosyl mimetics in the development of signal transduction inhibitors

Phosphotyrosyl (pTyr) residues play important roles in cellular signal transduction by facilitating recognition and binding necessary for critical protein-protein interactions, and for this reason pTyr motifs represent attractive starting points in the development of signaling antagonists. Although the pTyr phosphoryl moiety is central in these phenomena, its incorporation into signaling inhibitors is contraindicated due to enzymatic lability and limited bioavailability associated with phosphate esters. To address these limitations, an entire field of study has arisen devoted to the design and utilization of pTyr mimetics. This Account provides a perspective on the roles of pTyr residues in signal transduction and approaches to pTyr mimetic development.

Role of solution conformation and flexibility of short peptide ligands that bind to the p56(lck) SH2 domain

A general approach in drug design is making ligands more rigid in order to avoid loss in conformational entropy (deltaS(conf)) upon receptor binding. We hypothesized that in the high affinity binding of pYEEI peptide ligands to the p56(lck) SH2 domain this loss in deltaS(conf) might be diminished due to preorganization of the fourfold negatively charged pYEEI peptide in the bound, extended, conformation. A thermodynamic analysis was performed on the peptides Ac-pYEEI-NH(2), Ac-pYAAI-NH(2) and Ac-pYGGI-NH(2) using surface plasmon resonance (SPR) competition experiments to assay affinity constants at different temperatures. To study the effect of solution conformation and flexibility a computational conformation analysis was performed from which low energy conformations in solution were calculated, and S(conf) estimated. It was found that the calculated low energy conformations for especially the pYE moiety in solution resemble that in the bound state. In the calculated minimum energy conformation in solution isoleucine is bent towards the pY aromatic ring, the occurrence of such conformation is experimentally confirmed by NMR. The estimated values for S(conf) of the EE- and AA-peptide were similar, suggesting no predominant role of preorganization of the solution conformation due to electrostatic repulsion. Apparently the thermodynamics obey the same entropy-enthalpy compensation relationship, which also was found to hold for other peptides and peptidomimetics binding to p60(src) family SH2 domains. The implications of the results for drug design are discussed.

Protein-protein interactions: mechanisms and modification by drugs

Protein-protein interactions form the proteinaceous network, which plays a central role in numerous processes in the cell. This review highlights the main structures, properties of contact surfaces, and forces involved in protein-protein interactions. The properties of protein contact surfaces depend on their functions. The characteristics of contact surfaces of short-lived protein complexes share some similarities with the active sites of enzymes. The contact surfaces of permanent complexes resemble domain contacts or the protein core. It is reasonable to consider protein-protein complex formation as a continuation of protein folding. The contact surfaces of the protein complexes have unique structure and properties, so they represent prospective targets for a new generation of drugs. During the last decade, numerous investigations have been undertaken to find or design small molecules that block protein dimerization or protein(peptide)-receptor interaction, or on the other hand, induce protein dimerization.

Design, synthesis, and evaluation of matrix metalloprotease inhibitors bearing cyclopropane-derived peptidomimetics as P1' and P2' replacements

We have previously used trisubstituted cyclopropanes as peptide replacements to induce conformational constraints in known pseudopeptide inhibitors of a number of important enzymes. Cyclopropane-derived peptide mimics are novel in that they are among the few replacements that locally orient the peptide backbone and the amino acid side chain in a predefined manner. Although these dipeptide isosteres have been employed to orient amino acid side chains mimicking the gauche(-) conformation of chi(1)-space, their ability to project the side chains into an anti orientation has not been evaluated. As a first step toward this goal, the conformationally constrained pseudopeptides 8 and 10 and their corresponding flexible analogues 9 and 11 were prepared and tested as inhibitors of matrix metalloproteinases (MMPs). These compounds are analogues of 4 and 5, which were known to be potent MMP inhibitors. The anti orientations of the isopropyl side chain in 8 and the aromatic ring in 10 relative to the peptide backbone substituents on the cyclopropane were predicted to correspond to the known orientations of the P1' and P2' side chains of 5 when bound to MMPs. Hence, 8 and 10 were designed explicitly to probe topological features of the S1' or the S2' binding pockets of the MMPs. They were also designed to explore the importance of the P1'-P2' amide group, which is known to form highly conserved hydrogen bonds in several MMP-inhibitor complexes, and the viability of introducing a retro amide linkage between P2' and P3'. Pseudopeptides 8 and 9 were found to be weak competitive inhibitors of a series of MMPs. Any entropically favorable conformational constraints that were induced by the cyclopropane in 8 were thus overwhelmed by the loss of the hydrogen bonding capability associated with the P1'-P2' amide group. On the other hand, compounds 10 and 11, which contain a P2'-P3' retro amide group, were modest competitive inhibitors of a series of MMPs. The results obtained for 10 and 11 suggest that there may be a loss of hydrogen bonding capability associated with introducing the P2'-P3' retro amide group. However, because the conformationally constrained pseudopeptide 10 was significantly more potent than its flexible analogue 11, trisubstituted cyclopropanes related to 3 may serve as useful rigid dipeptide replacements in some biologically active pseudopeptides.

Calorimetric and structural studies of 1,2,3-trisubstituted cyclopropanes as conformationally constrained peptide inhibitors of Src SH2 domain binding

Isothermal titration calorimetry and X-ray crystallography have been used to determine the structural and thermodynamic consequences associated with constraining the pTyr residue of the pYEEI ligand for the Src Homology 2 domain of the Src kinase (Src SH2 domain). The conformationally constrained peptide mimics that were used are cyclopropane-derived isosteres whereby a cyclopropane ring substitutes to the N-Calpha-Cbeta atoms of the phosphotyrosine. Comparison of the thermodynamic data for the binding of the conformationally constrained peptide mimics relative to their equivalent flexible analogues as well as a native tetrapeptide revealed an entropic advantage of 5-9 cal mol(-1) K(-1) for the binding of the conformationally constrained ligands. However, an unexpected drop in enthalpy for the binding of the conformationally constrained ligands relative to their flexible analogues was also observed. To evaluate whether these differences reflected conformational variations in peptide binding modes, we have determined the crystal structure of a complex of the Src SH2 domain bound to one of the conformationally constrained peptide mimics. Comparison of this new structure with that of the Src SH2 domain bound to a natural 11-mer peptide (Waksman et al. Cell 1993, 72, 779-790) revealed only very small differences. Hence, cyclopropane-derived peptides are excellent mimics of the bound state of their flexible analogues. However, a rigorous analysis of the structures and of the surface areas at the binding interface, and subsequent computational derivation of the energetic binding parameters, failed to predict the observed differences between the binding thermodynamics of the rigidified and flexible ligands, suggesting that the drop in enthalpy observed with the conformationally constrained peptide mimic arises from sources other than changes in buried surface areas, though the exact origin of the differences remains unclear.