Thermodynamic and Structural Effects of Conformational Constraints in Protein−Ligand Interactions. Entropic Paradoxy Associated with Ligand Preorganization

A high-resolution structural characterization and physicochemical study of how a peptoid binds to an oncoprotein MDM2

Peptoids are a promising drug modality targeting disease-related proteins, but how a peptoid engages in protein binding is poorly understood. This is primarily due to a lack of high-resolution peptoid-protein complex structures and systematic physicochemical studies. Here, we present the first crystal structure of a peptoid bound to a protein, providing high-resolution structural information about how a peptoid binds to a protein. We previously reported a rigid peptoid, oligo(N-substituted alanine) (oligo-NSA), and developed an oligo-NSA-type peptoid that binds to MDM2. X-ray crystallographic analysis of the peptoid bound to MDM2 showed that the peptoid recognizes the MDM2 surface predominantly through the interaction of the N-substituents, while the main chain acts as a scaffold. Additionally, conformational, thermodynamic, and kinetic analysis of the peptoid and its derivatives with a less rigid main chain revealed that rigidification of the peptoid main chain contributes to improving the protein binding affinity. This improvement is thermodynamically attributed to an increased magnitude of the binding enthalpy change, and kinetically to an increased association rate and decreased dissociation rate. This study provides invaluable insights into the design of protein-targeting peptoids.

Thermodynamic consequences of stapling side-chains on a peptide ligand using a lactam-bridge: A theoretical study on anti-angiogenic peptides targeting VEGF

Peptides are promising therapeutic agents for various biological targets due to their high efficacy and low toxicity, and the design of peptide ligands with high binding affinity to the target of interest is of utmost importance in peptide-based drug design. Introducing a conformational constraint to a flexible peptide ligand using a side-chain lactam-bridge is a convenient and efficient method to improve its binding affinity to the target. However, in general, such a small structural modification to a flexible ligand made with the intent of lowering the configurational entropic penalty for binding may have unintended consequences in different components of the binding enthalpy and entropy, including the configurational entropy component, which are still not clearly understood. Toward probing this, we examine different components of the binding enthalpy and entropy as well as the underlying structure and dynamics, for a side-chain lactam-bridged peptide inhibitor and its flexible analog forming complexes with vascular endothelial growth factor (VEGF), using all-atom molecular dynamics simulations. It is found that introducing a side-chain lactam-bridge constraint into the flexible peptide analog led to a gain in configurational entropy change but losses in solvation entropy, solute internal energy, and solvation energy changes upon binding, pinpointing the opportunities and challenges in drug design. The present study features an interplay between configurational and solvation entropy changes, as well as the one between binding enthalpy and entropy, in ligand-target binding upon imposing a conformational constraint into a flexible ligand.

Calculations of Absolute Free Energies, Enthalpies, and Entropies for Drug Binding

Computer simulation methods can aid in the rational design of drugs aimed at a specific target, typically a protein. The affinity of a drug for its target is given by the free energy of binding. Binding can be further characterized by the enthalpy and entropy changes in the process. Methods exist to determine exact free energies, enthalpies, and entropies that are dependent only on the quality of the potential model and adequate sampling of conformational degrees of freedom. Entropy and enthalpy are roughly an order of magnitude more difficult to calculate than the free energy. This project combines a replica exchange method for enhanced sampling, designed to be efficient for protein-sized systems, with free energy calculations. This approach, replica exchange with dynamical scaling (REDS), uses two conventional simulations at different temperatures so that the entropy can be found from the temperature dependence of the free energy. A third replica is placed between them, with a modified Hamiltonian that allows it to span the temperature range of the conventional replicas. REDS provides temperature-dependent data and aids in sampling. It is applied to the bromodomain-containing protein 4 (BRD4) system. We find that for the force fields used, the free energies are accurate but the entropies and enthalpies are not, with the entropic contribution being too positive. Reproducing the entropy and enthalpy of binding appears to be a more stringent test of the force fields than reproducing the free energy.

What to Make of Zero: Resolving the Statistical Noise from Conformational Reorganization in Alchemical Binding Free Energy Estimates with Metadynamics Sampling

We introduce the self-relative binding free energy (self-RBFE) approach to evaluate the intrinsic statistical variance of dual-topology alchemical binding free energy estimators. The self-RBFE is the relative binding free energy between a ligand and a copy of the same ligand, and its true value is zero. Nevertheless, because the two copies of the ligand move independently, the self-RBFE value produced by a finite-length simulation fluctuates and can be used to measure the variance of the model. The results of this validation provide evidence that a significant fraction of the errors observed in benchmark studies reflect the statistical fluctuations of unconverged estimates rather than the models' accuracy. Furthermore, we find that ligand reorganization is a significant contributing factor to the statistical variance of binding free energy estimates and that metadynamics-accelerated conformational sampling of the torsional degrees of freedom of the ligand can drastically reduce the time to convergence.

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.

Biological evaluation of carbohydrate-based aprepitant analogs for neuroblastoma treatment

Different studies using Aprepitant, a NK1R antagonist currently used as a clinical drug for treating chemotherapy-related nausea and vomiting, have demonstrated that pharmacological inhibition of NK1R effectively reduces the growth of several tumor types such as neuroblastoma (NB). In a previous work, we demonstrated that a series of carbohydrate-based Aprepitant analogs, derived from either d-galactose or l-arabinose, have shown high affinity and NK1R antagonistic activity with a broad-spectrum anticancer activity and an important selectivity. In this new study, we explore the selective cytotoxic effects of these derivatives for the treatment of NB. Furthermore, we describe the design and stereoselective synthesis of a new generation of d-glucose derivatives as Aprepitant analogs, supported by docking studies. This approach showed that most of our carbohydrate-based analogs are significantly more selective than Aprepitant. The galactosyl derivative 2α, has demonstrated a marked in vitro selective cytotoxic activity against NB, with IC50 values in the same range as those of Aprepitant and its prodrug Fosaprepitant. Interestingly, the derivative 2α has shown similar apoptotic effect to that of Aprepitant. Moreover, we can select the glucosyl amino derivative 10α as an interesting hit exhibiting higher in vitro cytotoxic activity against NB than Aprepitant, being 1.2 times more selective.

Small Cyclic Diglycolamides: Tautomerism, Solvent Extraction and Coordination with f-Elements: One Strain to Rule Them All

The search for new effective extractants is an important task for the management of high-level liquid waste (HLW) generated during the reprocessing of spent nuclear fuel. Here, we synthesized a series of diglycolamides with cyclic substituents for the first time. We disclosed their coordination with f-element nitrates [La(NO3)3 and UO2(NO3)2] by SC-XRD study and complexation properties toward Am(III), Ln(III), and U(VI) during solvent extraction from nitric acid solutions. Using dynamic nuclear magnetic resonance (NMR) and density functional theory (DFT) calculations, the importance of tautomerism in the extraction properties of diglycolamides was shown.

MicroED as a Powerful Tool for Structure Determination of Macrocyclic Drug Compounds Directly from Their Powder Formulations

Macrocycles are important drug leads with many advantages including the ability to target flat and featureless binding sites as well as to act as molecular chameleons and thereby reach intracellular targets. However, due to their complex structures and inherent flexibility, macrocycles are difficult to study structurally, and there are limited structural data available. Herein, we use the cryo-EM method MicroED to determine the novel atomic structures of several macrocycles that have previously resisted structural determination. We show that structures of similar complexity can now be obtained rapidly from nanograms of material and that different conformations of flexible compounds can be derived from the same experiment. These results will have an impact on contemporary drug discovery as well as natural product exploration.

Computed Protein-Protein Enthalpy Signatures as a Tool for Identifying Conformation Sampling Problems

Understanding the thermodynamic signature of protein-peptide binding events is a major challenge in computational chemistry. The complexity generated by both components possessing many degrees of freedom poses a significant issue for methods that attempt to directly compute the enthalpic contribution to binding. Indeed, the prevailing assumption has been that the errors associated with such approaches would be too large for them to be meaningful. Nevertheless, we currently have no indication of how well the present methods would perform in terms of predicting the enthalpy of binding for protein-peptide complexes. To that end, we carefully assembled and curated a set of 11 protein-peptide complexes where there is structural and isothermal titration calorimetry data available and then computed the absolute enthalpy of binding. The initial "out of the box" calculations were, as expected, very modest in terms of agreement with the experiment. However, careful inspection of the outliers allows for the identification of key sampling problems such as distinct conformations of peptide termini not being sampled or suboptimal cofactor parameters. Additional simulations guided by these aspects can lead to a respectable correlation with isothermal titration calorimetry (ITC) experiments (R2 of 0.88 and an RMSE of 1.48 kcal/mol overall). Although one cannot know prospectively whether computed ITC values will be correct or not, this work shows that if experimental ITC data are available, then this in conjunction with computed ITC, can be used as a tool to know if the ensemble being simulated is representative of the true ensemble or not. That is important for allowing the correct interpretation of the detailed dynamics of the system with respect to the measured enthalpy. The results also suggest that computational calorimetry is becoming increasingly feasible. We provide the data set as a resource for the community, which could be used as a benchmark to help further progress in this area.

MicroED as a powerful tool for structure determination of macrocyclic drug compounds directly from their powder formulations

Macrocycles are important drug leads with many advantages including the ability to target flat and featureless binding sites as well as act as molecular chameleons and thereby reach intracellular targets. However, due to their complex structures and inherent flexibility, macrocycles are difficult to study structurally and there are limited structural data available. Herein, we use the cryo-EM method MicroED to determine the novel atomic structures of several macrocycles which have previously resisted structural determination. We show that structures of similar complexity can now be obtained rapidly from nanograms of material, and that different conformations of flexible compounds can be derived from the same experiment. These results will have impact on contemporary drug discovery as well as natural product exploration.

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.

Boronic Acid: A Novel Pharmacophore Targeting Src Homology 2 (SH2) Domain of STAT3

SH2 domains have been recognized as promising targets for various human diseases. However, targeting SH2 domains with phosphopeptides or small-molecule inhibitors derived from bioisosteres of the phosphate group is still challenging. Identifying novel bioisosteres of the phosphate group to achieve favorable in vivo potency is urgently needed. Here, we report the feasibility of targeting the STAT3-SH2 domain with a boronic acid group and the identification of a highly potent inhibitor compound 7 by replacing the carboxylic acid of compound 4 with a boronic acid. Compound 7 shows higher binding affinity, better cellular potency, more favorable PK profiles, and higher in vivo antitumor activity than 4. The stronger anticancer effect of 7 partially stems from its covalent binding mode with the SH2 domain, verified by the washout experiments. The relatively high level of sequence conservation among SH2 domains makes the results presented here of general significance.

Forces Driving a Magic Bullet to Its Target: Revisiting the Role of Thermodynamics in Drug Design, Development, and Optimization

Drug discovery strategies have advanced significantly towards prioritizing target selectivity to achieve the longstanding goal of identifying “magic bullets” amongst thousands of chemical molecules screened for therapeutic efficacy. A myriad of emerging and existing health threats, including the SARS-CoV-2 pandemic, alarming increase in bacterial resistance, and potentially fatal chronic ailments, such as cancer, cardiovascular disease, and neurodegeneration, have incentivized the discovery of novel therapeutics in treatment regimens. The design, development, and optimization of lead compounds represent an arduous and time-consuming process that necessitates the assessment of specific criteria and metrics derived via multidisciplinary approaches incorporating functional, structural, and energetic properties. The present review focuses on specific methodologies and technologies aimed at advancing drug development with particular emphasis on the role of thermodynamics in elucidating the underlying forces governing ligand–target interaction selectivity and specificity. In the pursuit of novel therapeutics, isothermal titration calorimetry (ITC) has been utilized extensively over the past two decades to bolster drug discovery efforts, yielding information-rich thermodynamic binding signatures. A wealth of studies recognizes the need for mining thermodynamic databases to critically examine and evaluate prospective drug candidates on the basis of available metrics. The ultimate power and utility of thermodynamics within drug discovery strategies reside in the characterization and comparison of intrinsic binding signatures that facilitate the elucidation of structural–energetic correlations which assist in lead compound identification and optimization to improve overall therapeutic efficacy.

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.

Transition Metal Catalysis Controlled by Hydrogen Bonding in the Second Coordination Sphere

Transition metal catalysis is of utmost importance for the development of sustainable processes in academia and industry. The activity and selectivity of metal complexes are typically the result of the interplay between ligand and metal properties. As the ligand can be chemically altered, a large research focus has been on ligand development. More recently, it has been recognized that further control over activity and selectivity can be achieved by using the "second coordination sphere", which can be seen as the region beyond the direct coordination sphere of the metal center. Hydrogen bonds appear to be very useful interactions in this context as they typically have sufficient strength and directionality to exert control of the second coordination sphere, yet hydrogen bonds are typically very dynamic, allowing fast turnover. In this review we have highlighted several key features of hydrogen bonding interactions and have summarized the use of hydrogen bonding to program the second coordination sphere. Such control can be achieved by bridging two ligands that are coordinated to a metal center to effectively lead to supramolecular bidentate ligands. In addition, hydrogen bonding can be used to preorganize a substrate that is coordinated to the metal center. Both strategies lead to catalysts with superior properties in a variety of metal catalyzed transformations, including (asymmetric) hydrogenation, hydroformylation, C-H activation, oxidation, radical-type transformations, and photochemical reactions.

Photomodulation of ultrastable host–guest complexes in water and their application in light-controlled steroid release

Host–guest complexation of dithienylethene photoswitches with cucurbit[8]uril leads to photoresponsive binding pairs with picomolar affinity in water.

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.

Ligand Binding Introduces Significant Allosteric Shifts in the Locations of Protein Fluctuations

Allostery is usually considered to be a mechanism for transmission of signals associated with physical or dynamic changes in some part of a protein. Here, we investigate the changes in fluctuations across the protein upon ligand binding based on the fluctuations computed with elastic network models. These results suggest that binding reduces the fluctuations at the binding site but increases fluctuations at remote sites, but not to fully compensating extents. If there were complete conservation of entropy, then only the enthalpies of binding would matter and not the entropies; however this does not appear to be the case. Experimental evidence also suggests that energies and entropies of binding can compensate but that the extent of compensation varies widely from case to case. Our results do however always show transmission of an allosteric signal to distant locations where the fluctuations are increased. These fluctuations could be used to compute entropies to improve evaluations of the thermodynamics of binding. We also show the allosteric relationship between peptide binding in the GroEL trans-ring that leads directly to the release of GroES from the GroEL-GroES cis-ring. This finding provides an example of how calculating these changes to protein dynamics induced by the binding of an allosteric ligand can regulate protein function and mechanism.

Anchor extension: a structure-guided approach to design cyclic peptides targeting enzyme active sites

Despite recent success in computational design of structured cyclic peptides, de novo design of cyclic peptides that bind to any protein functional site remains difficult. To address this challenge, we develop a computational "anchor extension" methodology for targeting protein interfaces by extending a peptide chain around a non-canonical amino acid residue anchor. To test our approach using a well characterized model system, we design cyclic peptides that inhibit histone deacetylases 2 and 6 (HDAC2 and HDAC6) with enhanced potency compared to the original anchor (IC₅₀ values of 9.1 and 4.4 nM for the best binders compared to 5.4 and 0.6 µM for the anchor, respectively). The HDAC6 inhibitor is among the most potent reported so far. These results highlight the potential for de novo design of high-affinity protein-peptide interfaces, as well as the challenges that remain.

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-Ac₆c-Asn-(CH₂)_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.

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.

Molecular Docking Studies of a Cyclic Octapeptide-Cyclosaplin from Sandalwood

Natural products from plants, such as chemopreventive agents, attract huge attention because of their low toxicity and high specificity. The rational drug design in combination with structure-based modeling and rapid screening methods offer significant potential for identifying and developing lead anticancer molecules. Thus, the molecular docking method plays an important role in screening a large set of molecules based on their free binding energies and proposes structural hypotheses of how the molecules can inhibit the target. Several peptide-based therapeutics have been developed to combat several health disorders, including cancers, metabolic disorders, heart-related diseases, and infectious diseases. Despite the discovery of hundreds of such therapeutic peptides however, only few peptide-based drugs have made it to the market. Moreover, the in silico activities of cyclic peptides towards molecular targets, such as protein kinases, proteases, and apoptosis related proteins have not been extensively investigated. In this study, we explored the in silico kinase and protease inhibitor potentials of cyclosaplin, and studied the interactions of cyclosaplin with other apoptosis-related proteins. Previously, the structure of cyclosaplin was elucidated by molecular modeling associated with dynamics that were used in the current study as well. Docking studies showed strong affinity of cyclosaplin towards cancer-related proteins. The binding affinity closer to 10 kcal/mol indicated efficient binding. Cyclosaplin showed strong binding affinities towards protein kinases such as EGFR, VEGFR2, PKB, and p38, indicating its potential role in protein kinase inhibition. Moreover, it displayed strong binding affinity to apoptosis-related proteins and revealed the possible role of cyclosaplin in apoptotic cell death. The protein–ligand interactions using LigPlot displayed some similar interactions between cyclosaplin and peptide-based ligands, especially in case of protein kinases and a few apoptosis related proteins. Thus, the in silico analyses gave the insights of cyclosaplin being a potential apoptosis inducer and protein kinase inhibitor.

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.

Broad-scale analysis of thermodynamic signatures in medicinal chemistry: are enthalpy-favored binders the better development option?

Thermodynamic profiles of ligand binding, particularly enthalpically favored binding signatures, have been suggested as a criterion to support the decision-making process around which compounds to select for further optimization in drug development. The concept was enthusiastically taken up, but turned out to be too superficial, either because many aspects determining thermodynamic profiles are insufficiently appreciated or because it is difficult to compare such data on a global scale. The impact of water, changes in protonation states, along with buffer dependencies and incompatible measurement conditions that are far from standard conditions hamper such broad-scale comparisons. However, thermodynamic signatures can make us aware of the impact of these aspects and provide important hints for improving our understanding of the binding process and defining criteria for drug optimization.

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.

Intramolecular H-bonds govern the recognition of a flexible peptide by an antibody

Molecular recognition is a fundamental event at the core of essentially every biological process. In particular, intermolecular H-bonds have been recognized as key stabilizing forces in antibody-antigen interactions resulting in exquisite specificity and high affinity. Although equally abundant, the role of intramolecular H-bonds is far less clear and not universally acknowledged. Herein, we have carried out a molecular-level study to dissect the contribution of intramolecular H-bonds in a flexible peptide for the recognition by an antibody. We show that intramolecular H-bonds may have a profound, multifaceted and favorable effect on the binding affinity by up to 2 kcal mol-1 of free energy. Collectively, our results suggest that antibodies are fine tuned to recognize transiently stabilized structures of flexible peptides in solution, for which intramolecular H-bonds play a key role.

CDNs-STING Interaction Mechanism Investigations and Instructions on Design of CDN-Derivatives

Cyclic dinucleotides (CDNs) present thousand-fold differences of dissociation constants to STING, a pivotal protein in cytosolic dsDNA immunity. To understand how subtle chemical changes in CDNs lead to these substantial variances, a precise ranking of binding affinity is needed. However, the large size and flexibility of CDNs elevate the entropic effect and pose a challenge for this precise prediction. Therefore, in this paper, we developed a new protocol, a combination of selective-integrated tempering sampling of ligands and molecular docking, to take into account the entropic effects originating from extensive ligand configurational space and solvation on binding affinity evaluations. The calculated ranking orders of CDNs and CDN-derivatives to wild type STING and R232H mutant are in agreement with experimental measurements. Further molecular dynamics analysis revealed that the interaction between phosphonate groups and 232R differentiates the binding affinities. The 2'-5' linked phosphonate groups have a larger tendency to form hydrogen bonds with 232R than those with 3'-5' linkages. Moreover, the new protocol identified structural features that enhanced CDNs-STING binding, such as anti-glycosidic bonds and large pro-R distances, which explains the high binding affinity of dithio-RpRp-2'3'-CDA to STING and is expected to provide valuable guidance in the lead-drug optimization.

Measuring Entropy in Molecular Recognition by Proteins

Molecular recognition by proteins is fundamental to the molecular basis of biology. Dissection of the thermodynamic landscape governing protein-ligand interactions has proven difficult because determination of various entropic contributions is quite challenging. Nuclear magnetic resonance relaxation measurements, theory, and simulations suggest that conformational entropy can be accessed through a dynamical proxy. Here, we review the relationship between measures of fast side-chain motion and the underlying conformational entropy. The dynamical proxy reveals that the contribution of conformational entropy can range from highly favorable to highly unfavorable and demonstrates the potential of this key thermodynamic variable to modulate protein-ligand interactions. The dynamical so-called entropy meter also refines the role of solvent entropy and directly determines the loss in rotational-translational entropy that occurs upon formation of high-affinity complexes. The ability to quantify the roles of entropy through an entropy meter based on measurable dynamical properties promises to highlight its role in protein function.

Enhanced Permeability and Binding Activity of Isobutylene-Grafted Peptides

We present a new peptide-macrocyclization strategy with an isobutylene graft. The reaction is mild and proceeds rapidly and efficiently both for linear and cyclic peptides. The resulting isobutylene-grafted peptides possess improved passive membrane permeability due to the shielding of the polar backbone of the amides, as demonstrated by NMR spectroscopy and molecular dynamics simulations. The isobutylene-stapled structures are fully stable in human plasma and in the presence of glutathione. This strategy can be applied to bioactive cyclic peptides such as somatostatin. Importantly, we found that structural preorganization forced by the isobutylene graft leads to a significant improvement in binding. The combined advantages of directness, selectivity, and smallness could allow application to peptide macrocyclization based on this attachment of the isobutylene graft.

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.

Preorganization and Cooperation for Highly Efficient and Reversible Capture of Low-Concentration CO2 by Ionic Liquids

A novel strategy based on the concept of preorganization and cooperation has been designed for a superior capacity to capture low-concentration CO₂ by imide-based ionic liquids. By using this strategy, for the first time, an extremely high gravimetric CO₂ capacity of up to 22 wt % (1.65 mol mol^-1 ) and excellent reversibility (16 cycles) have been achieved from 10 vol. % of CO₂ in N₂ when using an ionic liquid having a preorganized anion. Through a combination of quantum-chemical calculations and spectroscopic investigations, it is suggested that cooperative interactions between CO₂ and multiple active sites in the preorganized anion are the driving force for the superior CO₂ capacity and excellent reversibility.

From Peptides to Peptidomimetics: A Strategy Based on the Structural Features of Cyclopropane

Peptidomimetics, non-natural mimicries of bioactive peptides, comprise an important class of drug molecules. The essence of the peptidomimetic design is to mimic the key conformation assumed by the bioactive peptides upon binding to their targets. Regulation of the conformation of peptidomimetics is important not only to enhance target binding affinity and selectivity, but also to confer cell-membrane permeability for targeting protein-protein interactions in cells. The rational design of peptidomimetics with suitable three-dimensional structures is challenging, however, due to the inherent flexibility of peptides and their dynamic conformational changes upon binding to the target biomolecules. In this Minireview, a three-dimensional structural diversity-oriented strategy based on the characteristic structural features of cyclopropane to address this challenging issue in peptidomimetic chemistry is described.

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.

Development of a multipurpose scaffold for the display of peptide loops

Protein-protein interactions (PPIs) determine a wide range of biological processes and analysis of these dynamic networks is increasingly becoming a mandatory tool for studying protein function. Using the globular ATPase domain of recombinase RadA as a scaffold, we have developed a peptide display system (RAD display), which allows for the presentation of target peptides, protein domains or full-length proteins and their rapid recombinant production in bacteria. The design of the RAD display system includes differently tagged versions of the scaffold, which allows for flexibility in the protein purification method, and chemical coupling for small molecule labeling or surface immobilization. When combined with the significant thermal stability of the RadA protein, these features create a versatile multipurpose scaffold system. Using various orthogonal biophysical techniques, we show that peptides displayed on the scaffold bind to their natural targets in a fashion similar to linear parent peptides. We use the examples of CK2β/CK2α kinase and TPX2/Aurora A kinase protein complexes to demonstrate that the peptide displayed by the RAD scaffold can be used in PPI studies with the same binding efficacy but at lower costs compared with their linear synthetic counterparts.

Cyclopropane-Based Peptidomimetics Mimicking Wide-Ranging Secondary Structures of Peptides: Conformational Analysis and Their Use in Rational Ligand Optimization

Detailed conformational analyses of our previously reported cyclopropane-based peptidomimetics and conformational analysis-driven ligand optimization are described. Computational calculations and X-ray crystallography showed that the characteristic features of cyclopropane function effectively to constrain the molecular conformation in a three-dimensionally diverse manner. Subsequent principal component analysis revealed that the diversity covers the broad chemical space filled by peptide secondary structures in terms of both main-chain and side-chain conformations. Based on these analyses, a lead stereoisomer targeting melanocortin receptors was identified, and its potency and subtype selectivity were improved by further derivatization. The presented strategy is effective not only for designing non-peptidic ligands from a peptide ligand but also for the rational optimization of these ligands based on the plausible target-binding conformation without requiring the three- dimensional structural information of the target and its peptide ligands.

Synthesis and conformational analysis of peptides embodying 2,3-methanopipecolic acids

When 2,3-methanopipecolic acids replace a proline in peptides, a marked preference (42–92%) for thecisgeometry around the pipecolic amide bond is observed in both water and organic solvents.

Characterization of Promiscuous Binding of Phosphor Ligands to Breast-Cancer-Gene 1 (BRCA1) C-Terminal (BRCT): Molecular Dynamics, Free Energy, Entropy and Inhibitor Design

Inhibition of the protein-protein interaction (PPI) mediated by breast-cancer-gene 1 C-terminal (BRCT) is an attractive strategy to sensitize breast and ovarian cancers to chemotherapeutic agents that induce DNA damage. Such inhibitors could also be used for studies to understand the role of this PPI in DNA damage response. However, design of BRCT inhibitors is challenging because of the inherent flexibility associated with this domain. Several studies identified short phosphopeptides as tight BRCT binders. Here we investigated the thermodynamic properties of 18 phosphopeptides or peptide with phosphate mimic and three compounds with phosphate groups binding to BRCT to understand promiscuous molecular recognition and guide inhibitor design. We performed molecular dynamics (MD) simulations to investigate the interactions between inhibitors and BRCT and their dynamic behavior in the free and bound states. MD simulations revealed the key role of loops in altering the shape and size of the binding site to fit various ligands. The mining minima (M2) method was used for calculating binding free energy to explore the driving forces and the fine balance between configuration entropy loss and enthalpy gain. We designed a rigidified ligand, which showed unfavorable experimental binding affinity due to weakened enthalpy. This was because it lacked the ability to rearrange itself upon binding. Investigation of another phosphate group containing compound, C1, suggested that the entropy loss can be reduced by preventing significant narrowing of the energy well and introducing multiple new compound conformations in the bound states. From our computations, we designed an analog of C1 that introduced new intermolecular interactions to strengthen attractions while maintaining small entropic penalty. This study shows that flexible compounds do not always encounter larger entropy penalty, compared with other more rigid binders, and highlights a new strategy for inhibitor design.

Promiscuous proteins are commonly observed in biological systems, such as modular domains that recognize phosphopeptides during signal transduction. The use of phosphopeptides and compounds with phosphate groups as inhibitors to protein–protein interactions have attracted increasing interest for years. By using atomistic molecular dynamics simulations, we are able to perform detailed analyses of the dihedral space to explore protein fluctuation upon ligand binding to better understand promiscuous molecular recognition. Free energy calculation can further provide insights into the mechanism of binding, including both enthalpic and entropic contributions for molecular recognition, which assist in inhibitor design. Our calculation results show that pre-rigidifying a ligand is not always advantageous, suggesting the challenge in retaining optimized intermolecular interactions in pre-rigidified ligand. Instead, certain flexible ligands with multiple binding conformations can reduce entropic penalty, and therefore improves binding affinity. According to our computations, we can introduce new intermolecular interactions to flexible ligand to strengthen attractions while maintaining small entropic penalty by retaining its plasticity in the bound conformation. The study might cast light on a new general strategy for designing inhibitors targeting promiscuous modular domains and protein–protein interactions.

Experimental Binding Energies in Supramolecular Complexes

On the basis of many literature measurements, a critical overview is given on essential noncovalent interactions in synthetic supramolecular complexes, accompanied by analyses with selected proteins. The methods, which can be applied to derive binding increments for single noncovalent interactions, start with the evaluation of consistency and additivity with a sufficiently large number of different host-guest complexes by applying linear free energy relations. Other strategies involve the use of double mutant cycles, of molecular balances, of dynamic combinatorial libraries, and of crystal structures. Promises and limitations of these strategies are discussed. Most of the analyses stem from solution studies, but a few also from gas phase. The empirically derived interactions are then presented on the basis of selected complexes with respect to ion pairing, hydrogen bonding, electrostatic contributions, halogen bonding, π-π-stacking, dispersive forces, cation-π and anion-π interactions, and contributions from the hydrophobic effect. Cooperativity in host-guest complexes as well as in self-assembly, and entropy factors are briefly highlighted. Tables with typical values for single noncovalent free energies and polarity parameters are in the Supporting Information.

Dynamic Docking of Conformationally Constrained Macrocycles: Methods and Applications

Many natural products consist of large and flexible macrocycles that engage their targets via multiple contact points. This combination of contained flexibility and large contact area often allows natural products to bind at target surfaces rather than deep pockets, making them attractive scaffolds for inhibiting protein-protein interactions and other challenging therapeutic targets. The increasing ability to manipulate such compounds either biosynthetically or via semisynthetic modification means that these compounds can now be considered as starting points for medchem campaigns rather than solely as ends. Modern medchem benefits substantially from rational improvements made on the basis of molecular docking. As such, docking methods have been enhanced in recent years to deal with the complicated binding modalities and flexible scaffolds of macrocyclic natural products and natural product-like structures. Here, we comprehensively review methods for treating and docking these large macrocyclic scaffolds and discuss some of the resulting advances in medicinal chemistry.

Boosting Affinity by Correct Ligand Preorganization for the S2 Pocket of Thrombin: A Study by Isothermal Titration Calorimetry, Molecular Dynamics, and High-Resolution Crystal Structures

Structural preorganization to fix bioactive conformations at protein binding sites is a popular strategy to enhance binding affinity during late-stage optimization. The rationale for this enhancement relates to entropic advantages assigned to rigidified versus flexible ligands. We analyzed a narrow series of peptidomimetics binding to thrombin. The individual ligands exhibit at P2 a conformationally flexible glycine, more restricted alanine, N-methylglycine, N-methylhomoalanine, and largely rigidified proline moiety. Overall, affinity was found to increase by a factor of 1000, explained partly by an entropic advantage. All ligands adopt the same binding mode with small deviations. The residual mobility of the bound ligands is decreased across the series, and a protein side chain differs in its order/disorder behavior along with changes in the surface-water network pattern established across the newly generated protein-ligand surfaces. The enthalpy/entropy inventory displays a rather complex picture and emphasizes that thermodynamics can only be compared in terms of relative differences within a structurally similar ligand series.