Thermodynamic and Structural Effects of Macrocyclization as a Constraining Method in Protein-Ligand Interactions

Different molecular recognition by three domains of the full-length GRB2 to SOS1 proline-rich motifs and EGFR phosphorylated sites

The adaptor protein human GRB2 plays crucial roles in mediating signal transduction from cell membrane receptors to RAS and its downstream proteins by recruiting SOS1. Recent studies have revealed that GRB2 also serves as a scaffold for liquid-liquid phase separation (LLPS) with SOS1 and transmembrane receptors, which is thought to regulate the magnitude of cell signalling pathways. In this study, we employed solution NMR spectroscopy to investigate the interactions of the full-length GRB2 with proline-rich motifs (PRMs) derived from ten potential GRB2-binding sites in SOS1, as well as a peptide from a phosphorylation site of EGFR. Our findings indicate that the binding affinity of the two SH3 domains of GRB2 for PRMs differs by a factor of ten to twenty, with the N-terminal SH3 domain (NSH3) exhibiting a markedly higher affinity. The interactions of PRMs with the SH3 domains affected not only the regions surrounding the PRM binding sites on the SH3 domains but also the linker area connecting the three domains and parts of the SH2 domain. Analysis of the interaction between the phosphorylated EGFR binding site and the SH2 domain revealed chemical shift perturbations in regions distal from the known binding site of SH2. Moreover, we observed that the inter-domain interactions of the two SH3 domains with the SH2 domain of GRB2 are asymmetric. These findings suggest that the local binding of PRMs and phosphorylated EGFR to GRB2 impacts the overall structure of the GRB2 molecule, including domain orientation and dimerisation, which may contribute to LLPS formation.

Discovery of macrocyclic covalent inhibitors for severe acute respiratory syndrome coronavirus 2 3CL protease

The coronavirus disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been spread worldwide for more than 3 years. Although the hospitalization rate and mortality have decreased dramatically due to wide vaccination effort and improved treatment options, the disease is still a global health issue due to constant viral mutations, causing negative impact on social and economic activities. In addition, long COVID and complications arising from COVID-19 weeks after infection have become a concern for public health experts. Therefore, better treatments for COVID-19 are still needed. Herein, we describe a class of macrocyclic peptidomimetic compounds that are potent inhibitors of SARS-Cov-2 3CL protease (3CLpro). Significantly, some of the compounds showed a higher stability against human liver microsomes (HLM t1/2 > 180 min) and may be suitable for oral administration without the need for a pharmacokinetic (PK) boosting agent such as ritonavir.

Structure and Dynamics of Drk-SH2 Domain and Its Site-Specific Interaction with Sev Receptor Tyrosine Kinase

The Drosophila downstream receptor kinase (Drk), a homologue of human GRB2, participates in the signal transduction from the extracellular to the intracellular environment. Drk receives signals through the interaction of its Src homology 2 (SH2) domain with the phosphorylated tyrosine residue in the receptor tyrosine kinases (RTKs). Here, we present the solution NMR structure of the SH2 domain of Drk (Drk-SH2), which was determined in the presence of a phosphotyrosine (pY)-containing peptide derived from a receptor tyrosine kinase, Sevenless (Sev). The solution structure of Drk-SH2 possess a common SH2 domain architecture, consisting of three β strands imposed between two α helices. Additionally, we interpret the site-specific interactions of the Drk-SH2 domain with the pY-containing peptide through NMR titration experiments. The dynamics of Drk-SH2 were also analysed through NMR-relaxation experiments as well as the molecular dynamic simulation. The docking simulations of the pY-containing peptide onto the protein surface of Drk-SH2 provided the orientation of the peptide, which showed a good agreement with the analysis of the SH2 domain of GRB2.

Non-oxidized and oxidized flaxseed orbitides differently induce HepG2 cell apoptosis: involvement of cellular uptake and membrane death receptor DR4

BACKGROUND

Flaxseed orbitides have health-promoting properties, particularly potent anti-cancer activity. However, flaxseed orbitides containing a methionine structure, such as [1-9-NαC]-linusorb B2 (CLB), are easily oxidized to sulfoxide ([1-9-NαC],[1-Rs,Ss-MetO]-linusorb-B2 (CLC)) and sulfone ([1-9-NαC], [1-MetO]-linusorb B2 (CLK)), with CLC having less anti-cancer ability than CLB. It is unclear why oxidized flaxseed orbitides are less effective against cancer than non-oxidized flaxseed orbitide.

RESULTS

Non-oxidized ([1-9-NαC]-linusorb-B3 (CLA) and CLB) and oxidized (CLC and CLK) flaxseed orbitides were found to significantly upregulate the levels of pro-apoptotic proteins, including Bax/Bcl-2, CytoC, caspase-3, and caspase-8, in a dose-dependent manner, with non-oxidized flaxseed orbitides being more effective than oxidized flaxseed orbitides. Mechanically, the cellular absorption of non-oxidized flaxseed orbitides was higher than that of oxidized flaxseed orbitides. Moreover, the significant fluorescence quenching of DR4 protein by flaxseed orbitides (especially non-oxidized orbitides) indicated the formation of a DR4-orbitide complex. Molecular docking demonstrated that non-oxidized orbitides could easily dock into the active cavity of DR4 protein. Further blocking DR4 significantly reduced the ability of non-oxidized flaxseed orbitides to stimulate caspase-3 expression, whereas oxidized flaxseed orbitides retained this ability.

CONCLUSION

Non-oxidized flaxseed orbitides are more effective against cancer than oxidized flaxseed orbitides due to higher cellular uptake and activation of the DR4-mediated death receptor signaling pathway. © 2024 Society of Chemical Industry.

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.

An Intriguing Purview on the Design of Macrocyclic Inhibitors for Unexplored Protein Kinases through Their Binding Site Comparison

Kinases play an important role in regulating various intracellular signaling pathways that control cell proliferation, differentiation, survival, and other cellular processes, and their deregulation causes more than 400 diseases. Consequently, macrocyclization can be considered a noteworthy approach to developing new therapeutic agents for human diseases. Macrocyclization has emerged as an effective drug discovery strategy over the past decade to improve target selectivity and potency of small molecules. Small compounds with linear structures upon macrocyclization can lead to changes in their physicochemical and biological properties by firmly reducing conformational flexibility. A number of distinct protein kinases exhibit similar binding sites. Comparison of protein binding sites provides crucial insights for drug discovery and development. Binding site similarities are helpful in understanding polypharmacology, identifying potential off-targets, and repurposing known drugs. In this review, we focused on comparing the binding sites of those kinases for which macrocyclic inhibitors are available/studied so far. Furthermore, we calculated the volume of the binding site pocket for each targeted kinase and then compared it with the binding site pocket of the kinase for which only acyclic inhibitors were designed to date. Our review and analysis of several explored kinases might be useful in targeting new protein kinases for macrocyclic drug discovery.

Thermodynamic characterization of a macrocyclic Zika virus NS2B/NS3 protease inhibitor and its acyclic analogs

Cyclization of small molecules is a widely applied strategy in drug design for ligand optimization to improve affinity, as it eliminates the putative need for structural preorganization of the ligand before binding, or to improve pharmacokinetic properties. In this work, we provide a deeper insight into the binding thermodynamics of a macrocyclic Zika virus NS2B/NS3 protease inhibitor and its linear analogs. Characterization of the thermodynamic binding profiles by isothermal titration calorimetry experiments revealed an unfavorable entropy of the macrocycle compared to the open linear reference ligands. Molecular dynamic simulations and X-ray crystal structure analysis indicated only minor benefits from macrocyclization to fixate a favorable conformation, while linear ligands retained some flexibility even in the protein-bound complex structure, possibly explaining the initially surprising effect of a higher entropic penalty for the macrocyclic ligand.

Marine Cyclic Peptides: Antimicrobial Activity and Synthetic Strategies

Oceans are a rich source of structurally unique bioactive compounds from the perspective of potential therapeutic agents. Marine peptides are a particularly interesting group of secondary metabolites because of their chemistry and wide range of biological activities. Among them, cyclic peptides exhibit a broad spectrum of antimicrobial activities, including against bacteria, protozoa, fungi, and viruses. Moreover, there are several examples of marine cyclic peptides revealing interesting antimicrobial activities against numerous drug-resistant bacteria and fungi, making these compounds a very promising resource in the search for novel antimicrobial agents to revert multidrug-resistance. This review summarizes 174 marine cyclic peptides with antibacterial, antifungal, antiparasitic, or antiviral properties. These natural products were categorized according to their sources—sponges, mollusks, crustaceans, crabs, marine bacteria, and fungi—and chemical structure—cyclic peptides and depsipeptides. The antimicrobial activities, including against drug-resistant microorganisms, unusual structural characteristics, and hits more advanced in (pre)clinical studies, are highlighted. Nocathiacins I–III (91–93), unnarmicins A (114) and C (115), sclerotides A (160) and B (161), and plitidepsin (174) can be highlighted considering not only their high antimicrobial potency in vitro, but also for their promising in vivo results. Marine cyclic peptides are also interesting models for molecular modifications and/or total synthesis to obtain more potent compounds, with improved properties and in higher quantity. Solid-phase Fmoc- and Boc-protection chemistry is the major synthetic strategy to obtain marine cyclic peptides with antimicrobial properties, and key examples are presented guiding microbiologist and medicinal chemists to the discovery of new antimicrobial drug candidates from marine sources.

An Insight into the Medicinal Chemistry Perspective of Macrocyclic Derivatives with Antitumor Activity: A Systematic Review

The profound pharmacological properties of macrocyclic compounds have led to their development as drugs. In conformationally pre-organized ring structures, the multiple functions and stereochemical complexity provided by the macrocycle result in high affinity and selectivity of protein targets while maintaining sufficient bioavailability to reach intracellular locations. Therefore, the construction of macrocycles is an ideal choice to solve the problem of “undruggable” targets. Inspection of 68 macrocyclic drugs on the market showed that 10 of them were used to treat cancer, but this structural class still has been poorly explored within drug discovery. This perspective considers the macrocyclic compounds used for anti-tumor with different targets, their advantages and disadvantages, and the various synthetic methods of them.

Recent Advances in the Chemistry and Therapeutic Evaluation of Naturally Occurring and Synthetic Withanolides

Natural products are a major source of biologically active compounds that make promising lead molecules for developing efficacious drug-like molecules. Natural withanolides are found in many flora and fauna, including plants, algae, and corals, that traditionally have shown multiple health benefits and are known for their anti-cancer, anti-inflammatory, anti-bacterial, anti-leishmaniasis, and many other medicinal properties. Structures of these withanolides possess a few reactive sites that can be exploited to design and synthesize more potent and safe analogs. In this review, we discuss the literature evidence related to the medicinal implications, particularly anticancer properties of natural withanolides and their synthetic analogs, and provide perspectives on the translational potential of these promising compounds.

Structure-Based Epitope Design: Toward a Greater Antibody-SARS-CoV-2 RBD Affinity

Efficient COVID-19 vaccines are widely acknowledged as the best way to end the global pandemic. SARS-CoV-2 receptor-binding domain (RBD) plays fundamental roles related to cell infection. Antibodies could be developed to target RBD and represent a potential approach for the neutralization of the virus. Epitopes used to produce antibodies are generally linear peptides and thus possess multiple confirmations that do not reflect the actual topology of the targeted part in the native protein. On the other hand, macrocyclic epitopes could constitute closer mimics of the native protein topology and, as such, could generate superior antibodies. In this study, we demonstrated the vital effect of the size and the three-dimensional shape of epitopes on the activity of the developed antibodies against the RBD of SARS-CoV-2. The molecular dynamics studies showed the greater stability of the cyclic epitopes compared with the linear counterparts, which was reflected in the affinity of the produced antibodies. The antibodies developed using macrocyclic epitopes showed superiority with respect to binding to RBD compared to antibodies formed from linear peptides. This study constitutes a roadmap for developing superior antibodies that could be used to inhibit the activity of SARS-CoV-2.

Synthesis and structural characterization of a monocarboxylic inhibitor for GRB2 SH2 domain

A monocarboxylic inhibitor was designed and synthesized to disrupt the protein-protein interaction (PPI) between GRB2 and phosphotyrosine-containing proteins. Biochemical characterizations show compound 7 binds with the Src homology 2 (SH2) domain of GRB2 and is more potent than EGFR₁₀₆₈ phosphopeptide 14-mer. X-ray crystallographic studies demonstrate compound 7 occupies the GRB2 binding site for phosphotyrosine-containing sequences and reveal key structural features for GRB2-inhibitor binding. This compound with a -1 formal charge offers a new direction for structural optimization to generate cell-permeable inhibitors for this key protein target of the aberrant Ras-MAPK signaling cascade.

Conformational Strain of Macrocyclic Peptides in Ligand-Receptor Complexes Based on Advanced Refinement of Bound-State Conformers

Macrocyclic peptides are an important modality in drug discovery, but molecular design is limited due to the complexity of their conformational landscape. To better understand conformational propensities, global strain energies were estimated for 156 protein-macrocyclic peptide cocrystal structures. Unexpectedly large strain energies were observed when the bound-state conformations were modeled with positional restraints. Instead, low-energy conformer ensembles were generated using xGen that fit experimental X-ray electron density maps and gave reasonable strain energy estimates. The ensembles featured significant conformational adjustments while still fitting the electron density as well or better than the original coordinates. Strain estimates suggest the interaction energy in protein-ligand complexes can offset a greater amount of strain for macrocyclic peptides than for small molecules and non-peptidic macrocycles. Across all molecular classes, the approximate upper bound on global strain energies had the same relationship with molecular size, and bound-state ensembles from xGen yielded favorable binding energy estimates.

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.

Photochemical Methods for Peptide Macrocyclisation

Photochemical reactions have been the subject of renewed interest over the last two decades, leading to the development of many new, diverse and powerful chemical transformations. More recently, these developments have been expanded to enable the photochemical macrocyclisation of peptides and small proteins. These constructs benefit from increased stability, structural rigidity and biological potency over their linear counterparts, providing opportunities for improved therapeutic agents. In this review, an overview of both the established and emerging methods for photochemical peptide macrocyclisation is presented, highlighting both the limitations and opportunities for further innovation in the field.

Binding Interactions of Peptide Aptamers

Peptide aptamers are short amino acid chains that are capable of binding specifically to ligands in the same way as their much larger counterparts, antibodies. Ligands of therapeutic interest that can be targeted are other peptide chains or loops located on the surface of protein receptors (e.g., GCPR), which take part in cell-to-cell communications either directly or via the intermediary of hormones or signalling molecules. To confer on aptamers the same sort of conformational rigidity that characterises an antibody binding site, aptamers are often constructed in the form of cyclic peptides, on the assumption that this will encourage stronger binding interactions than would occur if the aptamers were simply linear chains. However, no formal studies have been conducted to confirm the hypothesis that linear peptides will engage in stronger binding interactions with cyclic peptides than with other linear peptides. In this study, the interaction of a model cyclic decamer with a series of linear peptide constructs was compared with that of a linear peptide with the same sequence, showing that the cyclic configuration does confer benefits by increasing the strength of binding.

Some thermodynamic effects of varying nonpolar surfaces in protein-ligand interactions

Understanding how making structural changes in small molecules affects their binding affinities for targeted proteins is central to improving strategies for rational drug design. To assess the effects of varying the nature of nonpolar groups upon binding entropies and enthalpies, we designed and prepared a set of Grb2-SH2 domain ligands, Ac-pTyr-Ac6c-Asn-(CH2)n-R, in which the size and electrostatic nature of R groups at the pTyr+3 site were varied. The complexes of these ligands with the Grb2-SH2 domain were evaluated in a series of studies in which the binding thermodynamics were determined using isothermal titration calorimetry, and binding interactions were examined in crystallographic studies of two different complexes. Notably, adding nonpolar groups to the pTyr+3 site leads to higher binding affinities, but the magnitude and energetic origins of these effects vary with the nature of the R substituent. For example, enhancements to binding affinities using aliphatic R groups are driven by more favorable changes in binding entropies, whereas aryl R groups improve binding free energies through a combination of more favorable changes in binding enthalpies and entropies. However, enthalpy/entropy compensation plays a significant role in these associations and mitigates against any significant variation in binding free energies, which vary by only 0.8 kcal•mol-1, with changes in the electrostatic nature and size of the R group. Crystallographic studies show that differences in ΔG° or ΔH° correlate with buried nonpolar surface area, but they do not correlate with the total number of polar or van der Waals contacts. The relative number of ordered water molecules and relative order in the side chains at pTyr+3 correlate with differences in -TΔS°. Overall, these studies show that burial of nonpolar surface can lead to enhanced binding affinities arising from dominating entropy- or enthalpy-driven hydrophobic effects, depending upon the electrostatic nature of the apolar R group.

Design and Development of a Macrocyclic Series Targeting Phosphoinositide 3-Kinase δ

A macrocyclization approach has been explored on a series of benzoxazine phosphoinositide 3-kinase δ inhibitors, resulting in compounds with improved potency, permeability, and in vivo clearance while maintaining good solubility. The thermodynamics of binding was explored via surface plasmon resonance, and the binding of lead macrocycle 19 was found to be almost exclusively entropically driven compared with progenitor 18, which demonstrated both enthalpic and entropic contributions. The pharmacokinetics of macrocycle 19 was also explored in vivo, where it showed reduced clearance when compared with the progenitor 18. This work adds to the growing body of evidence that macrocyclization could provide an alternative and complementary approach to the design of small-molecule inhibitors, with the potential to deliver differentiated properties.

Structure-thermodynamics-relationships of hepatitis C viral NS3 protease inhibitors

Thermodynamic parameters were determined for structurally-related inhibitors of HCV NS3 protease to assess how binding entropies and enthalpies vary with incremental changes at the P2 and P3 inhibitor subsites. Changing the heterocyclic substituent at P2 from a pyridyl to a 7-methoxy-2-phenyl-4-quinolyl group leads to a 710-fold increase in affinity. Annelating a benzene ring onto a pyridine ring leads to quinoline-derived inhibitors having higher affinities, but the individual enthalpy and entropy contributions are markedly different for each ligand pair. Introducing a phenyl group at C2 of the heterocyclic ring at P2 uniformly leads to higher affinity analogs with more favorable binding entropies, while adding a methoxy group at C7 of the quinoline ring at P2 provides derivatives with more favorable binding enthalpies. Significant enthalpy/entropy compensation is observed for structural changes made to inhibitors lacking a 2-phenyl substituent, whereas favorable changes in both binding enthalpies and entropies accompany structural modifications when a 2-phenyl group is present. Overall, binding energetics of inhibitors having a 2-phenyl-4-quinolyl group at P2 are dominated by entropic effects, whereas binding of the corresponding norphenyl analogs are primarily enthalpy driven. Notably, the reversal from an entropy driven association to an enthalpy driven one for this set of inhibitors also correlates with alternate binding modes. When the steric bulk of the side chain at P3 is increased from a hydrogen atom to a tert-butyl group, there is a 770-fold improvement in affinity. The 30-fold increase resulting from the first methyl group is solely the consequence of a more favorable change in entropy, whereas subsequent additions of methyl groups leads to modest increases in affinity that arise primarily from incremental improvements in binding enthalpies accompanied with smaller favorable entropic contributions.

Strategies for Late-Stage Optimization: Profiling Thermodynamics by Preorganization and Salt Bridge Shielding

Structural fixation of a ligand in its bioactive conformation may, due to entropic reasons, improve affinity. We present a congeneric series of thrombin ligands with a variety of functional groups triggering preorganization prior to binding. Fixation in solution and complex formation have been characterized by crystallography, isothermal titration calorimetry (ITC), and molecular dynamics (MD) simulations. First, we show why these preorganizing modifications do not affect the overall binding mode and how key interactions are preserved. Next, we demonstrate how preorganization thermodynamics can be largely dominated by enthalpy rather than entropy because of the significant population of low-energy conformations. Furthermore, a salt bridge is shielded by actively reducing its surface exposure, thus leading to an enhanced enthalpic binding profile. Our results suggest that the consideration of the ligand solution ensemble by MD simulation is necessary to predict preorganizing modifications that enhance the binding behavior of already promising binders.

Structural and functional properties of Grb2 SH2 dimer in CD28 binding

Growth factor receptor-bound protein 2 (Grb2) is an adaptor protein that plays a critical role in cellular signal transduction. It contains a central Src homology 2 (SH2) domain flanked by two Src homology 3 (SH3) domains. Binding of Grb2 SH2 to the cytoplasmic region of CD28, phosphorylated Tyr (pY) containing the peptide motif pY-X-N-X, is required for costimulatory signaling in T cells. In this study, we purified the dimer and monomer forms of Grb2 SH2, respectively, and analyzed their structural and functional properties. Size exclusion chromatography analysis showed that both dimer and monomer exist as stable states. Thermal stability analysis using circular dichroism showed that the dimer mostly dissociates into the monomer around 50°C. CD28 binding experiments showed that the affinity of the dimer to the phosphopeptide was about three fold higher than that of the monomer, possibly due to the avidity effect. The present crystal structure analysis of Grb2 SH2 showed two forms; one is monomer at 1.15 Å resolution, which is currently the highest resolution analysis, and another is dimer at 2.00 Å resolution. In the dimer structure, the C-terminal region, comprising residues 123–152, was extended towards the adjacent molecule, in which Trp121 was the hinge residue. The stable dimer purified using size exclusion chromatography would be due to the C-terminal helix “swapping”. In cases where a mutation caused Trp121 to be replaced by Ser in Grb2 SH2, this protein still formed dimers, but lost the ability to bind CD28.

Protein-ligand binding enthalpies from near-millisecond simulations: Analysis of a preorganization paradox

Calorimetric studies of protein-ligand binding sometimes yield thermodynamic data that are difficult to understand. Today, molecular simulations can be used to seek insight into such calorimetric puzzles, and, when simulations and experiments diverge, the results can usefully motivate further improvements in computational methods. Here, we apply near-millisecond duration simulations to estimate the relative binding enthalpies of four peptidic ligands with the Grb2 SH2 domain. The ligands fall into matched pairs, where one member of each pair has an added bond that preorganizes the ligand for binding and thus may be expected to favor binding entropically, due to a smaller loss in configurational entropy. Calorimetric studies have shown that the constrained ligands do in fact bind the SH2 domain more tightly than the flexible ones, but, paradoxically, the improvement in affinity for the constrained ligands is enthalpic, rather than entropic. The present enthalpy calculations yield the opposite trend, as they suggest that the flexible ligands bind more exothermically. Additionally, the small relative binding enthalpies are found to be balances of large differences in the energies of structural components such as ligand and the binding site residues. As a consequence, the deviations from experiment in the relative binding enthalpies represent small differences between these large numbers and hence may be particularly susceptible to error, due, for example, to approximations in the force field. We also computed first-order estimates of changes in configurational entropy on binding. These too are, arguably, paradoxical, as they tend to favor binding of the flexible ligands. The paradox is explained in part by the fact that the more rigid constrained ligands reduce the entropy of binding site residues more than their flexible analogs do, at least in the simulations. This result offers a rather general counterargument to the expectation that preorganized ligands should be associated with more favorable binding entropies, other things being equal.

Present drug-likeness filters in medicinal chemistry during the hit and lead optimization process: how far can they be simplified?

During the past decade, decreasing the attrition rate of drug development candidates reaching the market has become one of the major challenges in pharmaceutical research and drug development (R&D). To facilitate the decision-making process, and to increase the probability of rapidly finding and developing high-quality compounds, a variety of multiparametric guidelines, also known as rules and ligand efficiency (LE) metrics, have been developed. However, what are the 'best' descriptors and how far can we simplify these drug-likeness prediction tools in terms of the numerous, complex properties that they relate to?

Flexibility is important for inhibition of the MDM2/p53 protein-protein interaction by cyclic β-hairpins

Protein-protein interactions that have large, flat and featureless binding sites are difficult drug targets. In the development of their modulators conventional drug discovery strategies are often unsuccessful. Gaining a detailed understanding of the binding mode of protein-protein interaction inhibitors is therefore of vast importance for their future pharmaceutical use. The MDM2/p53 protein pair is a highly promising target for cancer treatment. Disruption of the protein complex using p53 α-helix mimetics has been shown to be a successful strategy to control p53 activity. To gain further insight into the binding of inhibitors to MDM2, the flexibility of four cyclic β-hairpins that act as α-helical mimetics and potential MDM2/p53 interaction inhibitors was investigated in relation to their inhibitory activity. MDM2-binding of the mimetics was determined using fluorescence polarization and surface plasmon resonance assays, whereas their conformation and dynamics in solution was described by the combined experimental and computational NAMFIS analysis. Molecular flexibility was shown to be important for the activity of the cyclic β-hairpin based MDM2 inhibitors.

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.

Flexibility vs rigidity of amphipathic peptide conjugates when interacting with lipid bilayers

For the first time, the photoisomerization of a diarylethene moiety (DAET) in peptide conjugates was used to probe the effects of molecular rigidity/flexibility on the structure and behavior of model peptides bound to lipid membranes. The DAET unit was incorporated into the backbones of linear peptide-based constructs, connecting two amphipathic sequences (derived from the β-stranded peptide (KIGAKI)₃ and/or the α-helical peptide BP100). A β-strand-DAET-α-helix and an α-helix-DAET-α-helix models were synthesized and studied in phospholipid membranes. Light-induced photoisomerization of the linker allowed the generation of two forms of each conjugate, which differed in the conformational mobility of the junction between the α-helical and/or the β-stranded part of these peptidomimetic molecules. A detailed study of their structural, orientational and conformational behavior, both in isotropic solution and in phospholipid model membranes, was carried out using circular dichroism and solid-state ¹⁹F-NMR spectroscopy. The study showed that the rigid and flexible forms of the two conjugates had appreciably different structures only when embedded in an anisotropic lipid environment and only in the gel phase. The influence of the rigidity/flexibility of the studied conjugates on the lipid thermotropic phase transition was also investigated by differential scanning calorimetry. Both models were found to destabilize the lamellar gel phases.

•

DAET building blocks can be used to study rigidity/flexibility effects in supramolecular model systems.

•

Photoswitchable DAET linkers perturb only up to 3–4 adjacent amino acid residues.

•

Membrane-bound amphiphilic secondary structure elements exert a negligible influence on each other when linked by DAET.

•

The rigidity of peptide conjugates affected their structural behavior only in the lipid gel phase.

Increased Conformational Flexibility of a Macrocycle-Receptor Complex Contributes to Reduced Dissociation Rates

Constraining a peptide in its bioactive conformation by macrocyclization represents a powerful strategy to design modulators of challenging biomolecular targets. This holds particularly true for the development of inhibitors of protein‐protein interactions which often involve interfaces lacking defined binding pockets. Such flat surfaces are demanding targets for traditional small molecules rendering macrocyclic peptides promising scaffolds for novel therapeutics. However, the contribution of peptide dynamics to binding kinetics is barely understood which impedes the design process. Herein, we report unexpected trends in the binding kinetics of two closely related macrocyclic peptides that bind their receptor protein with high affinity. Isothermal titration calorimetry, ¹⁹F NMR experiments and molecular dynamics simulations reveal that increased conformational flexibility of the macrocycle–receptor complex reduces dissociation rates and contributes to complex stability. This observation has impact on macrocycle design strategies that have so far mainly focused on the stabilization of bioactive ligand conformations.

Ligand Strain and Entropic Effects on the Binding of Macrocyclic and Linear Inhibitors: Molecular Modeling of Penicillopepsin Complexes

Using extensive molecular dynamics simulations, we investigate the structure and dynamics of the complexes formed between penicillopepsin and two peptidomimetic inhibitors: a linear compound, isovaleryl(P₄)-valine(P₃)-asparagine(P₂)-leucine(P₁)-phosphonate-phenylalanine(P₁'), and its macrocylic analog that includes a methylene bridge between the Asn(P₂) and Phe(P₁') side chains. The macrocyclic inhibitor, which has a 420-fold stronger affinity than that of the acyclic one, has been considered to lower the entropic penalty for binding. To better understand this binding preference, the solution structure of the inhibitors is studied by molecular dynamics simulations. Subsequently, we assess the influence of the enzyme/inhibitor contacts, the enzyme-induced inhibitor strain, the variation of the ligand configurational entropy and the enzyme reorganization by combining molecular-mechanics Poisson-Boltzmann surface area and normal mode calculations with conformational entropy calculations. We find that there is no relevant entropic stabilization on the binding of the cyclic inhibitor with respect to the acyclic analog because the methylene bridge does not reduce appreciably the conformational flexibility of the free inhibitor. The most important factors explaining the stronger affinity of the macrocyclic inhibitor are the conformational filtering and the lower ligand strain induced by the methylene bridge.

Elucidating the energetic contributions to the binding free energy

New exact equations are derived for the terms contributing to the binding free energy (ΔG⁰) of a ligand-receptor pair using our recently introduced formalism which we here call perturbation-divergence formalism (PDF). Specifically, ΔG⁰ equals the sum of the average of the perturbation (pertaining to new interactions) and additional dissipative terms. The average of the perturbation includes the sum of the average receptor-ligand interactions and the average of the change of solvation energies upon association. The Kullback-Leibler (KL) divergence quantifies the energetically dissipative terms, which are due to the configurational changes and, using the chain rule of KL divergence, can be decomposed into (i) dissipation due to limiting the external liberation (translation and rotation) of the ligand relative to the receptor and (ii) dissipation due to conformational (internal) changes inside the receptor and the ligand. We also identify all exactly canceling energetic terms which do not contribute to ΔG⁰. Furthermore, the PDF provides a new approach towards dimensionality reduction in the representation of the association process and towards relating the dynamic (high dimensional) with the thermodynamic (one-dimensional) changes.

A novel complexity-to-diversity strategy for the diversity-oriented synthesis of structurally diverse and complex macrocycles from quinine

Recent years have witnessed a global decline in the productivity and advancement of the pharmaceutical industry. A major contributing factor to this is the downturn in drug discovery successes. This can be attributed to the lack of structural (particularly scaffold) diversity and structural complexity exhibited by current small molecule screening collections. Macrocycles have been shown to exhibit a diverse range of biological properties, with over 100 natural product-derived examples currently marketed as FDA-approved drugs. Despite this, synthetic macrocycles are widely considered to be a poorly explored structural class within drug discovery, which can be attributed to their synthetic intractability. Herein we describe a novel complexity-to-diversity strategy for the diversity-oriented synthesis of novel, structurally complex and diverse macrocyclic scaffolds from natural product starting materials. This approach exploits the inherent structural (including functional) and stereochemical complexity of natural products in order to rapidly generate diversity and complexity. Readily-accessible natural product-derived intermediates serve as structural templates which can be divergently functionalized with different building blocks to generate a diverse range of acyclic precursors. Subsequent macrocyclisation then furnishes compounds that are each based around a distinct molecular scaffold. Thus, high levels of library scaffold diversity can be rapidly achieved. In this proof-of-concept study, the natural product quinine was used as the foundation for library synthesis, and six novel structurally diverse, highly complex and functionalized macrocycles were generated.

How Molecular Conformational Changes Affect Changes in Free Energy

A simple quantitative relationship between the molecular conformational changes and the corresponding changes in the free energy is presented. The change in free energy is the sum of that part of the enthalpic change that is due to the externally applied work (perturbation) and of that part of the entropic change, termed dissipative entropy, that is related to the conformational changes. The dissipative entropy is equivalent to the relative entropy, a concept from information theory, between the distributions of the conformations in the initial and the final states. The remaining change in entropy (nondissipative) cancels exactly with the remaining enthalpic change. The calculation of the dissipative entropy is demonstrated to pose the main difficulty in free energy computation. The straightforward decomposition of the dissipative entropy into contributions from different parts of the system promises to improve the understanding of the role of conformational changes in biochemical reactions.

Counting on natural products for drug design

Natural products and their molecular frameworks have a long tradition as valuable starting points for medicinal chemistry and drug discovery. Recently, there has been a revitalization of interest in the inclusion of these chemotypes in compound collections for screening and achieving selective target modulation. Here we discuss natural-product-inspired drug discovery with a focus on recent advances in the design of synthetically tractable small molecules that mimic nature's chemistry. We highlight the potential of innovative computational tools in processing structurally complex natural products to predict their macromolecular targets and attempt to forecast the role that natural-product-derived fragments and fragment-like natural products will play in next-generation drug discovery.

The Role of Conformational Changes in Molecular Recognition

Conformational changes of molecules are crucial elements in many biochemical processes, and also in molecular recognition. Here, we present a novel exact mathematical equation for the binding free energy of a receptor-ligand pair. It shows that the energetic contribution due to conformational changes upon molecular recognition is defined by the so-called Kullback-Leibler (KL) divergence between the probability distributions of the conformational ensemble of the biomolecule in the bound and free states. We show that conformational changes always contribute positively to the change in free energy and therefore disfavor the association process. Using the example of ligands binding to a flexible cavity of T4 lysozyme, we illustrate that, due to enthalpy-entropy compensation, the conformational entropy is a misleading quantity for assessing the conformational contribution to the binding free energy, in contrast to the KL divergence, which is the correct quantity to use in this context.

Polyketide synthase and non-ribosomal peptide synthetase thioesterase selectivity: logic gate or a victim of fate?

Thioesterases (TEs) are product offloading enzymes from FAS, PKS, and NRPS complexes. We review the diversity, structure, and mechanism of PKS and NRPS TEs and analyze TE loading and release steps as possible logic gates with a view to predicting TE function in new pathways.

First total synthesis of cyclodepsipeptides clavatustide A and B and their enantiomers

The enantiopure synthesis of clavatustides A (1) and B (3) were accomplished by a seven step synthetic protocol starting from commercially available (R)-phenyllactic acid.

Structural and biophysical investigation of the interaction of a mutant Grb2 SH2 domain (W121G) with its cognate phosphopeptide

The adaptor protein Grb2 is a key element of mitogenetically important signaling pathways. With its SH2 domain it binds to upstream targets while its SH3 domains bind to downstream proteins thereby relaying signals from the cell membranes to the nucleus. The Grb2 SH2 domain binds to its targets by recognizing a phosphotyrosine (pY) in a pYxNx peptide motif, requiring an Asn at the +2 position C-terminal to the pY with the residue either side of this Asn being hydrophobic. Structural analysis of the Grb2 SH2 domain in complex with its cognate peptide has shown that the peptide adopts a unique β-turn conformation, unlike the extended conformation that phosphopeptides adopt when bound to other SH2 domains. TrpEF1 (W121) is believed to force the peptide into this unusual conformation conferring this unique specificity to the Grb2 SH2 domain. Using X-ray crystallography, electron paramagnetic resonance (EPR) spectroscopy, and isothermal titration calorimetry (ITC), we describe here a series of experiments that explore the role of TrpEF1 in determining the specificity of the Grb2 SH2 domain. Our results demonstrate that the ligand does not adopt a pre-organized structure before binding to the SH2 domain, rather it is the interaction between the two that imposes the hairpin loop to the peptide. Furthermore, we find that the peptide adopts a similar structure when bound to both the wild-type Grb2 SH2 domain and a TrpEF1Gly mutant. This suggests that TrpEF1 is not the determining factor for the conformation of the phosphopeptide.

Solid-Phase Parallel Synthesis of Functionalised Medium-to-Large Cyclic Peptidomimetics through Three-Component Coupling Driven by Aziridine Aldehyde Dimers

The first solid-phase parallel synthesis of macrocyclic peptides using three-component coupling driven by aziridine aldehyde dimers is described. The method supports the synthesis of 9- to 18-membered aziridine-containing macrocycles, which are then functionalized by nucleophilic opening of the aziridine ring. This constitutes a robust approach for the rapid parallel synthesis of macrocyclic peptides.