A marine analgesic peptide, Contulakin-G, and neurotensin are distinct agonists for neurotensin receptors: uncovering structural determinants of desensitization properties.

Neurotensin receptors have been studied as molecular targets for the treatment of pain, schizophrenia, addiction, or cancer. Neurotensin (NT) and Contulakin-G, a glycopeptide isolated from a predatory cone snail Conus geographus, share a sequence similarity at the C-terminus, which is critical for activation of neurotensin receptors. Both peptides are potent analgesics, although affinity and agonist potency of Contulakin-G toward neurotensin receptors are significantly lower, as compared to those for NT. In this work, we show that the weaker agonist properties of Contulakin-G result in inducing significantly less desensitization of neurotensin receptors and preserving their cell-surface density. Structure-activity relationship (SAR) studies suggested that both glycosylation and charged amino acid residues in Contulakin-G or NT played important roles in desensitizing neurotensin receptors. Computational modeling studies of human neurotensin receptor NTS1 and Contulakin-G confirmed the role of glycosylation in weakening interactions with the receptors. Based on available SAR data, we designed, synthesized, and characterized an analog of Contulakin-G in which the glycosylated amino acid residue, Gal-GalNAc-Thr10, was replaced by memantine-Glu10 residue. This analog exhibited comparable agonist potency and weaker desensitization properties as compared to that of Contulakin-G, while producing analgesia in the animal model of acute pain following systemic administration. We discuss our study in the context of feasibility and safety of developing NT therapeutic agents with improved penetration across the blood-brain barrier. Our work supports engineering peptide-based agonists with diverse abilities to desensitize G-protein coupled receptors and further emphasizes opportunities for conotoxins as novel pharmacological tools and drug candidates.

Combined STAT3 and BCR-ABL1 inhibition induces synthetic lethality in therapy-resistant chronic myeloid leukemia.

Mutations in the BCR-ABL1 kinase domain are an established mechanism of tyrosine kinase inhibitor (TKI) resistance in Philadelphia chromosome-positive leukemia, but fail to explain many cases of clinical TKI failure. In contrast, it is largely unknown why some patients fail TKI therapy despite continued suppression of BCR-ABL1 kinase activity, a situation termed BCR-ABL1 kinase-independent TKI resistance. Here, we identified activation of signal transducer and activator of transcription 3 (STAT3) by extrinsic or intrinsic mechanisms as an essential feature of BCR-ABL1 kinase-independent TKI resistance. By combining synthetic chemistry, in vitro reporter assays, and molecular dynamics-guided rational inhibitor design and high-throughput screening, we discovered BP-5-087, a potent and selective STAT3 SH2 domain inhibitor that reduces STAT3 phosphorylation and nuclear transactivation. Computational simulations, fluorescence polarization assays and hydrogen-deuterium exchange assays establish direct engagement of STAT3 by BP-5-087 and provide a high-resolution view of the STAT3 SH2 domain/BP-5-087 interface. In primary cells from chronic myeloid leukemia (CML) patients with BCR-ABL1 kinase-independent TKI resistance, BP-5-087 (1.0 µM) restored TKI sensitivity to therapy-resistant CML progenitor cells, including leukemic stem cells. Our findings implicate STAT3 as a critical signaling node in BCR-ABL1 kinase-independent TKI resistance, and suggest that BP-5-087 has clinical utility for treating malignancies characterized by STAT3 activation.

BCR-ABL1 compound mutations combining key kinase domain positions confer clinical resistance to ponatinib in Ph chromosome-positive leukemia.

Ponatinib is the only currently approved tyrosine kinase inhibitor (TKI) that suppresses all BCR-ABL1 single mutants in Philadelphia chromosome-positive (Ph(+)) leukemia, including the recalcitrant BCR-ABL1(T315I) mutant. However, emergence of compound mutations in a BCR-ABL1 allele may confer ponatinib resistance. We found that clinically reported BCR-ABL1 compound mutants center on 12 key positions and confer varying resistance to imatinib, nilotinib, dasatinib, ponatinib, rebastinib, and bosutinib. T315I-inclusive compound mutants confer high-level resistance to TKIs, including ponatinib. In vitro resistance profiling was predictive of treatment outcomes in Ph(+) leukemia patients. Structural explanations for compound mutation-based resistance were obtained through molecular dynamics simulations. Our findings demonstrate that BCR-ABL1 compound mutants confer different levels of TKI resistance, necessitating rational treatment selection to optimize clinical outcome.

Epigenetic molecular recognition: a biomolecular modeling perspective.

The abnormal regulation of epigenetic protein families is associated with the onset and progression of various human diseases. However, epigenetic processes remain relatively obscure at the molecular level, thus preventing the rational design of chemical therapeutics. An array of robust computational and modeling approaches can complement experiments to shed light on the complex mechanisms of epigenetic molecular recognition and can guide medicinal chemists in designing selective and potent drug molecules. Herein we present a review of studies focused on epigenetic molecular recognition from a biomolecular modeling viewpoint. Although the known epigenetic targets are numerous, this review focuses on the more limited protein families on which computational modeling has been successfully applied. Therefore, we review three main topics: 1) histone deacetylases, 2) histone demethylases, and 3) histone tail dynamics. A brief review of the biological background and biomedical relevance is presented for each topic, followed by a detailed discussion of the computational studies and their relevance.

MDWiZ: a platform for the automated translation of molecular dynamics simulations.

A variety of popular molecular dynamics (MD) simulation packages were independently developed in the last decades to reach diverse scientific goals. However, such non-coordinated development of software, force fields, and analysis tools for molecular simulations gave rise to an array of software formats and arbitrary conventions for routine preparation and analysis of simulation input and output data. Different formats and/or parameter definitions are used at each stage of the modeling process despite largely contain redundant information between alternative software tools. Such Babel of languages that cannot be easily and univocally translated one into another poses one of the major technical obstacles to the preparation, translation, and comparison of molecular simulation data that users face on a daily basis. Here, we present the MDWiZ platform, a freely accessed online portal designed to aid the fast and reliable preparation and conversion of file formats that allows researchers to reproduce or generate data from MD simulations using different setups, including force fields and models with different underlying potential forms. The general structure of MDWiZ is presented, the features of version 1.0 are detailed, and an extensive validation based on GROMACS to LAMMPS conversion is presented. We believe that MDWiZ will be largely useful to the molecular dynamics community. Such fast format and force field exchange for a given system allows tailoring the chosen system to a given computer platform and/or taking advantage of a specific capabilities offered by different software engines.

PITOMBA: Parameter Interface for Oligosaccharide Molecules Based on Atoms.

A novel four-bead coarse-grained (CG) model for carbohydrates denoted PITOMBA was devised using a bottom-up approach based on the atomistic GROMOS 53A6GLYC force field and on experimental thermodynamical data. The model was developed to be used in conjunction with the SPC CG water model (J. Chem. Phys. 2011, 134, 084110) and the GROMOS force field functional form. Explicit electrostatic interactions are considered by assigning point charges to each CG bead. Validation of the model is presented to a variety of structural and thermodynamic properties for mono- and oligosaccharides in solution. In addition, the model development philosophy allows for prompt extensions to include hexopyranose chains with diverse glycosidic linkages and branches.

Coarse-Graining of TIP4P/2005, TIP4P-Ew, SPC/E, and TIP3P to Monatomic Anisotropic Water Models Using Relative Entropy Minimization.

Coarse-grained models are becoming a competitive alternative for modeling processes that occur over time and length scales beyond the reach of fully atomistic molecular simulations. Ideally, coarse-grained models should not only achieve high computational efficiency but also provide accurate predictions and fundamental insight into the role of molecular interactions, the characteristic behavior, and properties of the system they model. In this work we derive a series of monatomic coarse-grained water models mX(REM) from the most popular atomistic water models X = TIP3P, SPC/E, TIP4P-Ew, and TIP4P/2005, using the relative entropy minimization (REM) method. Each coarse-grained water molecule is represented by a single particle that interacts through short-ranged anisotropic interactions that encourage the formation of "hydrogen-bonded" structures. We systematically investigate the features of the coarse-grained models in reproducing over 20 structural, dynamic, and thermodynamic properties of the reference atomistic water models-including the existence and locus of the characteristic density anomaly. The mX(REM) coarse-grained models reproduce quite faithfully the radial and angular distribution function of water, produce a temperature of maximum density (TMD), and stabilize the ice I crystal. Moreover, the ratio between the TMD and the melting temperature of the crystal in the mX(REM) models and liquid-ice equilibrium properties show reasonable agreement with the results of the corresponding atomistic models. The mX(REM) models, however, severely underestimate the cohesive energy of the condensed water phases. We investigate which specific limitations of the coarse-grained models arise from the REM methodology, from the monatomic nature of the models, and from the Stillinger-Weber interaction potential form. Our analysis indicates that a small compromise in the accuracy of structural properties can result in a significant increase of the overall accuracy and representability of the coarse-grained water models. We evaluate the accuracy of the atomistic and the monatomic anisotropic coarse-grained water models, including the mW water model, in reproducing experimental water properties. We find that mW and mTIP4P/2005(REM) score closer to experiment than widely used atomistic water models. We conclude that monatomic models of water with short-range, anisotropic "hydrogen-bonding" three-body interactions can be competitive in accuracy with fully atomistic models for the study of a wide range of properties and phenomena at less than 1/100th of the computational cost.

Molecular dynamics simulations indicate an induced-fit mechanism for LSD1/CoREST-H3-histone molecular recognition.

BACKGROUND: Lysine Specific Demethylase (LSD1 or KDM1A) in complex with its co-repressor protein CoREST catalyzes the demethylation of the H3 histone N-terminal tail and is currently one of the most promising epigenetic targets for drug discovery against cancer and neurodegenerative diseases. Models of non-covalent binding, such as lock and key, induced-fit, and conformational selection could help explaining the molecular mechanism of LSD1/CoREST-H3-histone association, thus guiding drug discovery and design efforts. Here, we quantify the extent to which LSD1/CoREST substrate binding is consistent with these hypothetical models using LSD1/CoREST conformational ensembles obtained through extensive explicit solvent molecular dynamics (MD) simulations. RESULTS: We find that an induced-fit model is the most representative of LSD1/CoREST-H3-histone non-covalent binding and accounts for the local conformational changes occurring in the H3-histone binding site. We also show that conformational selection - despite in principle not ruled out by this finding - is minimal, and only relevant when global properties are considered, e.g. the nanoscale motion of the LSD1/CoREST clamp. CONCLUSION: The induced-fit mechanism revealed by our MD simulation study will aid the inclusion of protein dynamics for the discovery and design of LSD1 inhibitors targeting the H3-histone binding region. On a general basis, our study indicates the importance of using multiple metrics or selection schemes when testing alternative hypothetical mechanistic models of non-covalent binding.

Variational Implicit-Solvent Modeling of Host-Guest Binding: A Case Study on Cucurbit[7]uril|

The synthetic host cucurbit[7]uril (CB[7]) binds aromatic guests or metal complexes with ultrahigh affinity compared with that typically displayed in protein-ligand binding. Due to its small size, CB[7] serves as an ideal receptor-ligand system for developing computational methods for molecular recognition. Here, we apply the recently developed variational implicit-solvent model (VISM), numerically evaluated by the level-set method, to study hydration effects in the high-affinity binding of the B2 bicyclo[2.2.2]octane derivative to CB[7]. For the unbound host, we find that the host cavity favors the hydrated state over the dry state due to electrostatic effects. For the guest binding, we find reasonable agreement to experimental binding affinities. Dissection of the individual VISM free-energy contributions shows that the major driving forces are water-mediated hydrophobic interactions and the intrinsic (vacuum) host-guest van der Waals interactions. These findings are in line with recent experiments and molecular dynamics simulations with explicit solvent. It is expected that the level-set VISM, with further refinement on the electrostatic descriptions, can efficiently predict molecular binding and recognition in a wide range of future applications.

Substrate selection influences molecular recognition in a screen for lymphoid tyrosine phosphatase inhibitors.

Assay design is an important variable that influences the outcome of an inhibitor screen. Here, we have investigated the hypothesis that protein tyrosine phosphatase inhibitors with improved biological activity could be identified from a screen by using a biologically relevant peptide substrate, rather than traditional phosphotyrosine mimetic substrates. A 2000-member library of drugs and drug-like compounds was screened for inhibitors of lymphoid tyrosine phosphatase (LYP) by using both a peptide substrate (Ac-ARLIEDNE-pCAP-TAREG-NH2, peptide 1) and a small-molecule phosphotyrosine mimetic substrate (difluoromethyl umbelliferyl phosphate, DiFMUP). The results demonstrate that compounds that inhibited enzyme activity on the peptide substrate had greater biological activity than compounds that only inhibited enzyme activity on DiFMUP. Finally, epigallocatechin-3,5-digallate was identified as the most potent inhibitor of lymphoid tyrosine phosphatase activity to date, with an IC50 of 50 nM and significant activity in T-cells. Molecular docking simulations provided a first model for binding of this potent inhibitor to LYP; this will constitute the platform for ongoing lead optimization efforts.

Thiuram disulfides as pseudo-irreversible inhibitors of lymphoid tyrosine phosphatase.

We screened a small library of thiuram disulfides for inhibition of lymphoid tyrosine phosphatase (LYP) activity. The parent thiuram disulfide, disulfiram, inhibited LYP activity in vitro and in Jurkat T cells, whereas diethyldithiocarbamate failed to inhibit LYP at the concentrations tested. Compound 13, an N-(2-thioxothiazolidin-4-one) analogue, was found to be the most potent LYP inhibitor in this series, with an IC50 value of 3 µM. Compound 13 inhibits LYP pseudo-irreversibly, as evidenced by the time-dependence of inhibition, with a K(i) value of 1.1 µM and a k(inact) value of 0.004 s⁻¹. The inhibition of LYP by compound 13 could not be reversed significantly by incubation with glutathione or by prolonged dialysis, but could be partially reversed by incubation with dithiothreitol. Compound 13 also inhibited LYP activity in Jurkat T cells.

Expanding the druggable space of the LSD1/CoREST epigenetic target: new potential binding regions for drug-like molecules, peptides, protein partners, and chromatin.

Lysine specific demethylase-1 (LSD1/KDM1A) in complex with its corepressor protein CoREST is a promising target for epigenetic drugs. No therapeutic that targets LSD1/CoREST, however, has been reported to date. Recently, extended molecular dynamics (MD) simulations indicated that LSD1/CoREST nanoscale clamp dynamics is regulated by substrate binding and highlighted key hinge points of this large-scale motion as well as the relevance of local residue dynamics. Prompted by the urgent need for new molecular probes and inhibitors to understand LSD1/CoREST interactions with small-molecules, peptides, protein partners, and chromatin, we undertake here a configurational ensemble approach to expand LSD1/CoREST druggability. The independent algorithms FTMap and SiteMap and our newly developed Druggable Site Visualizer (DSV) software tool were used to predict and inspect favorable binding sites. We find that the hinge points revealed by MD simulations at the SANT2/Tower interface, at the SWIRM/AOD interface, and at the AOD/Tower interface are new targets for the discovery of molecular probes to block association of LSD1/CoREST with chromatin or protein partners. A fourth region was also predicted from simulated configurational ensembles and was experimentally validated to have strong binding propensity. The observation that this prediction would be prevented when using only the X-ray structures available (including the X-ray structure bound to the same peptide) underscores the relevance of protein dynamics in protein interactions. A fifth region was highlighted corresponding to a small pocket on the AOD domain. This study sets the basis for future virtual screening campaigns targeting the five novel regions reported herein and for the design of LSD1/CoREST mutants to probe LSD1/CoREST binding with chromatin and various protein partners.

Protein recognition by short peptide reversible inhibitors of the chromatin-modifying LSD1/CoREST lysine demethylase.

The combinatorial assembly of protein complexes is at the heart of chromatin biology. Lysine demethylase LSD1(KDM1A)/CoREST beautifully exemplifies this concept. The active site of the enzyme tightly associates to the N-terminal domain of transcription factors of the SNAIL1 family, which therefore can competitively inhibit the binding of the N-terminal tail of the histone substrate. Our enzymatic, crystallographic, spectroscopic, and computational studies reveal that LSD1/CoREST can bind to a hexapeptide derived from the SNAIL sequence through recognition of a positively charged α-helical turn that forms upon binding to the enzyme. Variations in sequence and length of this six amino acid ligand modulate affinities enabling the same binding site to differentially interact with proteins that exert distinct biological functions. The discovered short peptide inhibitors exhibit antiproliferative activities and lay the foundation for the development of peptidomimetic small molecule inhibitors of LSD1.

Molecular recognition and ligand association.

We review recent developments in our understanding of molecular recognition and ligand association, focusing on two major viewpoints: (a) studies that highlight new physical insight into the molecular recognition process and the driving forces determining thermodynamic signatures of binding and (b) recent methodological advances in applications to protein-ligand binding. In particular, we highlight the challenges posed by compensating enthalpic and entropic terms, competing solute and solvent contributions, and the relevance of complex configurational ensembles comprising multiple protein, ligand, and solvent intermediate states. As more complete physics is taken into account, computational approaches increase their ability to complement experimental measurements, by providing a microscopic, dynamic view of ensemble-averaged experimental observables. Physics-based approaches are increasingly expanding their power in pharmacology applications.

Solvent fluctuations in hydrophobic cavity-ligand binding kinetics.

Water plays a crucial part in virtually all protein-ligand binding processes in and out of equilibrium. Here, we investigate the role of water in the binding kinetics of a ligand to a prototypical hydrophobic pocket by explicit-water molecular dynamics (MD) simulations and implicit diffusional approaches. The concave pocket in the unbound state exhibits wet/dry hydration oscillations whose magnitude and time scale are significantly amplified by the approaching ligand. In turn, the ligand's stochastic motion intimately couples to the slow hydration fluctuations, leading to a sixfold-enhanced friction in the vicinity of the pocket entrance. The increased friction considerably decelerates association in the otherwise barrierless system, indicating the importance of molecular-scale hydrodynamic effects in cavity-ligand binding arising due to capillary fluctuations. We derive and analyze the diffusivity profile and show that the mean first passage time distribution from the MD simulation can be accurately reproduced by a standard Brownian dynamics simulation if the appropriate position-dependent friction profile is included. However, long-time decays in the water-ligand (random) force autocorrelation demonstrate violation of the Markovian assumption, challenging standard diffusive approaches for rate prediction. Remarkably, the static friction profile derived from the force correlations strongly resembles the profile derived on the Markovian assumption apart from a simple shift in space, which can be rationalized by a time-space retardation in the ligand's downhill dynamics toward the pocket. The observed spatiotemporal hydrodynamic coupling may be of biological importance providing the time needed for conformational receptor-ligand adjustments, typical of the induced-fit paradigm.

On the Role of Dewetting Transitions in Host-Guest Binding Free Energy Calculations.

We use thermodynamic integration (TI) and explicit solvent molecular dynamics (MD) simulation to estimate the absolute free energy of host-guest binding. In the unbound state, water molecules visit all of the internally accessible volume of the host, which is fully hydrated on all sides. Upon binding of an apolar guest, the toroidal host cavity is fully dehydrated; thus, during the intermediate λ stages along the integration, the hydration of the host fluctuates between hydrated and dehydrated states. Estimating free energies by TI can be especially challenging when there is a considerable difference in hydration between the two states of interest. We investigate these aspects using the popular TIP3P and TIP4P water models. TI free energy estimates through MD largely depend on water-related interactions, and water dynamics significantly affect the convergence of binding free energy calculations. Our results indicate that wetting/dewetting transitions play a major role in slowing the convergence of free energy estimation. We employ two alternative approaches-one analytical and the other empirically based on actual MD sampling-to correct for the standard state free energy. This correction is sizable (up to 4 kcal/mol), and the two approaches provide corrections that differ by about 1 kcal/mol. For the system considered here, the TIP4P water model combined with an analytical correction for the standard state free energy provides higher overall accuracy. This observation might be transferable to other systems in which water-related contributions dominate the binding process.

Water structure, dynamics, and spectral signatures: changes upon model cavity-ligand recognition.

It is becoming increasingly evident that water plays an active role in noncovalent receptor-ligand association. In this study, hydrophobic cavity-ligand association in a model system is characterized through the analysis of the structure, dynamics, and corresponding spectral signatures of water at different stages of the binding process. Molecular dynamics simulations reveal that the reorientation of the water molecules around the ligand becomes faster as the receptor-ligand distance reduces, which is correlated with the decrease in number of water-water hydrogen bonds within the ligand hydration shells. Prompted by the need for calculating physical quantities that can be amenable to experimental validation, the changes in the spectroscopic features upon cavity-ligand binding are investigated. The analysis of both linear and nonlinear infrared spectra allows direct insight into the evolution of water structure and dynamics around the ligand. In particular, characteristic spectroscopic features emerge at key stages of the binding process, which are related to changes in the hydrogen-bond topology of water around the ligand. This study demonstrates that computer simulations and vibrational spectroscopy could be integrated to facilitate the direct study of solvent effects in biomolecular association.

LSD1/CoREST is an allosteric nanoscale clamp regulated by H3-histone-tail molecular recognition.

The complex of lysine-specific demethylase-1 (LSD1/KDM1A) with its corepressor protein CoREST is an exceptionally relevant target for epigenetic drugs. Here, we provide insight into the local and global changes of LSD1/CoREST conformational dynamics that occur upon H3 binding on the basis of a total cumulative time of one microsecond molecular dynamics simulation. The LSD1/CoREST complex functions as an allosteric nanoscale-binding clamp, which is regulated by substrate binding. In the unbound state, LSD1/CoREST reversibly visits clamp states that are more open or significantly more closed compared with the available X-ray crystal structures. The Lys triad of residues Lys355, Lys357, and Lys359 gates the entrance of the H3 pocket. H3 binding shifts the pocket breathing dynamics toward open, higher-volume states while reducing the overall flexibility of the LSD1/CoREST nanoscale clamp. We show that the H3 pocket is an allosteric site for the regulation of the rotation of the amino oxidase domain with respect to the Tower domain. The allosteric mechanism relies on the specific reduction of nanoscale domain rotation upon local H3-tail binding. Instead, clamp opening/closing motions that do not involve domain rotation only reduce in amplitude yet are dominant in the bound state. Overall, our data suggest that the H3 binding pocket is a central target site to (i) switch off LSD1 amino oxidase activity, thus H3-tail demethylation; (ii) block the competitive binding of transcription factors; and (iii) prevent chromatin anchoring to LSD1/CoREST. This study underscores the importance of receptor flexibility for future epigenetic drug discovery.

LSD1/CoREST reversible opening-closing dynamics: discovery of a nanoscale clamp for chromatin and protein binding.

LSD1 associated with its corepressor protein CoREST is an exceptionally relevant target for epigenetic drugs. Hypotheses for the role of LSD1/CoREST as a multidocking site for chromatin and protein binding would require significant molecular flexibility, and LSD1/CoREST large-amplitude conformational dynamics is currently unknown. Here, molecular dynamics simulation reveals that the LSD1/CoREST complex in solution functions as a reversible nanoscale binding clamp. We show that the H3 histone tail binding pocket is a potential allosteric site for regulation of the rotation of SWIRM/SANT2 domains around the Tower domain. Thus, targeting this site and including receptor flexibility are crucial strategies for future drug discovery.