Ensemble-Based Virtual Screening Led to the Discovery of New Classes of Potent Tyrosinase Inhibitors

Probing the Activation Mechanisms of Agonist DPI-287 to Delta-Opioid Receptor and Novel Agonists Using Ensemble-Based Virtual Screening with Molecular Dynamics Simulations

Pain drugs targeting mu-opioid receptors face major addiction problems that have caused an epidemic. The delta-opioid receptor (DOR) has shown to not cause addictive effects when bound to an agonist. While the active conformation of the DOR in complex with agonist DPI-287 has been recently solved, there are still no FDA-approved agonists targeting it, providing the opportunity for structure-based virtual screening. In this study, the conformational plasticity of the DOR was probed using molecular dynamics (MD) simulations, identifying two representative conformations from clustering analysis. The two MD conformations as well as the crystal conformation of DOR were used to screen novel compounds from the ZINC database (17 million compounds), in which 69 drugs were picked as potential compounds based on their docking scores. Notably, 37 out of the 69 compounds were obtained from the simulated conformations. The binding stability of the 69 compounds was further investigated using MD simulations. Based on the MM-GBSA binding energy and the predicted drug properties, eight compounds were chosen as the most favorable, six of which were from the simulated conformations. Using a dynamic network model, the communication between the crystal agonist and the top eight molecules with the receptor was analyzed to confirm if these novel compounds share a similar activation mechanism to the crystal ligand. Encouragingly, docking of these eight compounds to the other two opioid receptors (kappa and mu) suggests their good selectivity toward DOR.

Recent discovery of tyrosinase inhibitors in traditional Chinese medicines and screening methods

ETHNOPHARMACOLOGICAL RELEVANCE

Tyrosinase, the key rate-limiting enzyme for melanogenesis, is one of the main targets for skin senescence and some pigmented skin diseases, such as albinism and melanoma. Tyrosinase inhibitors are capable of reducing melanin generation and deposition in the skin through blocking the reaction chain of formation. Thus, it has been used for anti-melanoma and showed the potential to be developed into novel skin whitening and spot removing products. With the trend of back-to-nature, natural tyrosinase inhibitors are receiving more and more attention. Traditional Chinese medicines (TCMs) as the promising source of novel chemotypes and pharmacophores, are huge treasures for the discovery of natural tyrosinase inhibitors characterized with green, safe, and highly efficient.

AIM OF THIS REVIEW

This review aims to provide a systematic overview of natural tyrosinase inhibitors and a detailed summary of the types of TCMs from which they originate. In addition, this paper also highlights the screening methods developed for exploring tyrosinase inhibitors in recent years, compares the advantages and disadvantages of various methods under the guidance of different screening principles, and predicts their applications in the future.

MATERIALS AND METHODS

Relevant literature have been obtained using the keywords "tyrosinase inhibitors", "traditional Chinese medicines", "whitening", and "screening" in scientific databases, such as "PubMed", "SciFinder", "Web of Science", "Elsevier", "China Knowledge Resource Integrated databases". Information was also collected from Chinese pharmacopoeia, Chinese herbal classics books, "Google Scholar", "Baidu Scholar", and other literature sources, etc. RESULTS: An overview about the tyrosinase inhibitors derived from TCMs since 2002 has been compiled via the above-mentioned sources. Up to now, 186 components, mainly belonging to flavonoids, lignans, terpenoids, Diels-Alder adducts, simple phenylpropanoids and stilbenes, from 61 kinds of TCMs have been reported to possess tyrosinase inhibitory activity, among which flavonoids are mainly focused on. Furthermore, on the basis of bioactive detection strategies, the screening methods for tyrosinase inhibitors have been classified into bioaffinity-based, intrinsic enzymatic-based, and computer-aided drug design (CADD). Precisely because screening approaches are essential for rapid identification of tyrosinase inhibitors from TCMs, the principles, advantages and disadvantages, and specific applications of each method are presented along with a comparison of applicability.

CONCLUSIONS

The summary of TCMs-derived inhibitors gives a clue on the discovery of candidates with the property to whiten the skin. Meanwhile, the outlook of developed screening methods provides technical references for the efficient exploration of safer and more effective tyrosinase inhibitors from TCMs.

Actives-Based Receptor Selection Strongly Increases the Success Rate in Structure-Based Drug Design and Leads to Identification of 22 Potent Cancer Inhibitors

Computer-aided drug design, an important component of the early stages of the drug discovery pipeline, routinely identifies large numbers of false positive hits that are subsequently confirmed to be experimentally inactive compounds. We have developed a methodology to improve true positive prediction rates in structure-based drug design and have successfully applied the protocol to twenty target systems and identified the top three performing conformers for each of the targets. Receptor performance was evaluated based on the area under the curve of the receiver operating characteristic curve for two independent sets of known actives. For a subset of five diverse cancer-related disease targets, we validated our approach through experimental testing of the top 50 compounds from a blind screening of a small molecule library containing hundreds of thousands of compounds. Our methods of receptor and compound selection resulted in the identification of 22 novel inhibitors in the low μM-nM range, with the most potent being an EGFR inhibitor with an IC₅₀ value of 7.96 nM. Additionally, for a subset of five independent target systems, we demonstrated the utility of Gaussian accelerated molecular dynamics to thoroughly explore a target system's potential energy surface and generate highly predictive receptor conformations.

Discovery of Tyrosinase Inhibitors: Structure-Based Virtual Screening and Biological Evaluation

Tyrosinase (EC 1.14.18.1) plays an indispensable role in the rate-limiting steps of melanin biosynthesis, and its uncontrolled activity may result in various diseases, such as albinism, melanoma, freckles, etc. The inhibition of tyrosinase activity may provide a useful and efficient strategy to treat hyperpigmentation disorders. However, the widely used tyrosinase inhibitors, like α-arbutin, hydroquinone, and kojic acid, have many shortcomings, such as lower efficacy and much more side effects. Herein, we reported the use of homology modeling and multistep structure-based virtual screening for the discovery of novel tyrosinase inhibitors. In this study, 10 initial potential hits (compounds T1–T10) were evaluated for enzyme inhibition and kinetic study, with kojic acid being used as a control. Among them, the IC50 values of both T1 (11.56 ± 0.98 μmol/L) and T5 (18.36 ± 0.82 μmol/L) were superior to that of kojic acid (23.12 ± 1.26 μmol/L). Moreover, T1 and T5 were also identified as the effective noncompetitive tyrosinase inhibitors by the subsequent kinetic study. Above all, T1 and T5 may represent the promising drug candidates for hyperpigmentation therapy in pharmaceutical fields, as well as the effective whitening agents in cosmetic applications.

Discovery of Kinase and Carbonic Anhydrase Dual Inhibitors by Machine Learning Classification and Experiments

A multi-target small molecule modulator is advantageous for treating complicated diseases such as cancers. However, the strategy and application for discovering a multi-target modulator have been less reported. This study presents the dual inhibitors for kinase and carbonic anhydrase (CA) predicted by machine learning (ML) classifiers, and validated by biochemical and biophysical experiments. ML trained by CA I and CA II inhibitor molecular fingerprints predicted candidates from the protein-specific bioactive molecules approved or under clinical trials. For experimental tests, three sulfonamide-containing kinase inhibitors, 5932, 5946, and 6046, were chosen. The enzyme assays with CA I, CA II, CA IX, and CA XII have allowed the quantitative comparison in the molecules’ inhibitory activities. While 6046 inhibited weakly, 5932 and 5946 exhibited potent inhibitions with 100 nM to 1 μM inhibitory constants. The ML screening was extended for finding CAs inhibitors of all known kinase inhibitors. It found XMU-MP-1 as another potent CA inhibitor with an approximate 30 nM inhibitory constant for CA I, CA II, and CA IX. Differential scanning fluorimetry confirmed the direct interaction between CAs and small molecules. Cheminformatics studies, including docking simulation, suggest that each molecule possesses two separate functional moieties: one for interaction with kinases and the other with CAs.

Hydroxamic Acid as a Potent Metal-Binding Group for Inhibiting Tyrosinase

Tyrosinase, a metalloenzyme containing a dicopper cofactor, plays a central role in synthesizing melanin from tyrosine. Many studies have aimed to identify small-molecule inhibitors of tyrosinase for pharmaceutical, cosmetic, and agricultural purposes. In this study, we report that hydroxamic acid is a potent metal-binding group for interacting with dicopper atoms, thereby inhibiting tyrosinase. Hydroxamate-containing molecules, including anticancer drugs targeting histone deacetylase, vorinostat and panobinostat, significantly inhibited mushroom tyrosinase, with inhibitory constants in the submicromolar range. Of the tested molecules, benzohydroxamic acid was the most potent. Its inhibitory constant of 7 nM indicates that benzohydroxamic acid is one of the most potent tyrosinase inhibitors. Results from differential scanning fluorimetry revealed that direct binding mediates inhibition. The enzyme kinetics were studied to assess the inhibitory mechanism of the hydroxamate-containing molecules. Experiments with B16F10 cell lysates confirmed that the new inhibitors are inhibitory against mammalian tyrosinase. Docking simulation data revealed intermolecular contacts between hydroxamate-containing molecules and tyrosinase.

New Inhibitors of Laccase and Tyrosinase by Examination of Cross-Inhibition between Copper-Containing Enzymes

Coppers play crucial roles in the maintenance homeostasis in living species. Approximately 20 enzyme families of eukaryotes and prokaryotes are known to utilize copper atoms for catalytic activities. However, small-molecule inhibitors directly targeting catalytic centers are rare, except for those that act against tyrosinase and dopamine-β-hydroxylase (DBH). This study tested whether known tyrosinase inhibitors can inhibit the copper-containing enzymes, ceruloplasmin, DBH, and laccase. While most small molecules minimally reduced the activities of ceruloplasmin and DBH, aside from known inhibitors, 5 of 28 tested molecules significantly inhibited the function of laccase, with the K_i values in the range of 15 to 48 µM. Enzyme inhibitory kinetics classified the molecules as competitive inhibitors, whereas differential scanning fluorimetry and fluorescence quenching supported direct bindings. To the best of our knowledge, this is the first report on organic small-molecule inhibitors for laccase. Comparison of tyrosinase and DBH inhibitors using cheminformatics predicted that the presence of thione moiety would suffice to inhibit tyrosinase. Enzyme assays confirmed this prediction, leading to the discovery of two new dual tyrosinase and DBH inhibitors.

Identification of Phytoconstituents as Potent Inhibitors of Casein Kinase-1 Alpha Using Virtual Screening and Molecular Dynamics Simulations

Casein kinase-1 alpha (CK1α) is a multifunctional protein kinase that belongs to the serine/threonine kinases of the CK1α family. It is involved in various signaling pathways associated with chromosome segregation, cell metabolism, cell cycle progression, apoptosis, autophagy, etc. It has been known to involve in the progression of many diseases, including cancer, neurodegeneration, obesity, and behavioral disorders. The elevated expression of CK1α in diseased conditions facilitates its selective targeting for therapeutic management. Here, we have performed virtual screening of phytoconstituents from the IMPPAT database seeking potential inhibitors of CK1α. First, a cluster of compounds was retrieved based on physicochemical parameters following Lipinski's rules and PAINS filter. Further, high-affinity hits against CK1α were obtained based on their binding affinity score. Furthermore, the ADMET, PAINS, and PASS evaluation was carried out to select more potent hits. Finally, following the interaction analysis, we elucidated three phytoconstituents, Semiglabrinol, Curcusone_A, and Liriodenine, posturing considerable affinity and specificity towards the CK1α binding pocket. The result was further evaluated by molecular dynamics (MD) simulations, dynamical cross-correlation matrix (DCCM), and principal components analysis (PCA), which revealed that binding of the selected compounds, especially Semiglabrinol, stabilizes CK1α and leads to fewer conformational fluctuations. The MM-PBSA analysis suggested an appreciable binding affinity of all three compounds toward CK1α.

Recent advances in the design and discovery of synthetic tyrosinase inhibitors

Tyrosinase is a copper-containing metalloenzyme that is responsible for the rate-limiting catalytic step in the melanin biosynthesis and enzymatic browning. As a promising target, tyrosinase inhibitors can be used as skin whitening agents and food preservatives, thus having broad potential in the fields of food, cosmetics, agriculture and medicine. From 2015 to 2020, numerous synthetic inhibitors of tyrosinase have been developed to overcome the challenges of low efficacy and side effects. This review summarizes the enzyme structure and biological functions of tyrosinase and demonstrates the recent advances of synthetic tyrosinase inhibitors from the perspective of medicinal chemistry, providing a better understanding of the catalytic mechanisms and more effective tyrosinase inhibitors.

A comprehensive review on tyrosinase inhibitors

Tyrosinase is a multi-copper enzyme which is widely distributed in different organisms and plays an important role in the melanogenesis and enzymatic browning. Therefore, its inhibitors can be attractive in cosmetics and medicinal industries as depigmentation agents and also in food and agriculture industries as antibrowning compounds. For this purpose, many natural, semi-synthetic and synthetic inhibitors have been developed by different screening methods to date. This review has focused on the tyrosinase inhibitors discovered from all sources and biochemically characterised in the last four decades.

The compromise of virtual screening and its impact on drug discovery

Introduction: Docking and structure-based virtual screening (VS) have been standard approaches in structure-based design for over two decades. However, our understanding of the limitations, potential, and strength of these techniques has enhanced, raising expectations. Areas covered: Based on a survey of reports in the past five years, we assess whether VS: (1) predicts binding poses in agreement with crystallographic data (when available); (2) is a superior screening tool, as often claimed; (3) is successful in identifying chemical scaffolds that can be starting points for subsequent lead optimization cycles. Data shows that knowledge of the target and its chemotypes in postprocessing lead to viable hits in early drug discovery endeavors. Expert opinion: VS is capable of accurate placements in the pocket for the most part, but does not consistently score screening collections accurately. What matters is capitalization on available resources to get closer to a viable lead or optimizable series. Integration of approaches, subjective hit selection guided by knowledge of the receptor or endogenous ligand, libraries driven by experimental guides, validation studies to identify the best docking/scoring that reproduces experimental findings, constraints regarding receptor-ligand interactions, thoroughly designed methodologies, and predefined cutoff scoring criteria strengthen VS's position in pharmaceutical research.

Predictive models for tyrosinase inhibitors: Challenges from heterogeneous activity data determined by different experimental protocols

Quantitative Structure-Activity Relationship (QSAR) models of tyrosinase inhibitors were built using Random Forest (RF) algorithm and evaluated by the out-of-bag estimation (R²_OOB) and 10-fold cross validation (Q²_CV). We found that the performances of QSAR models were closely correlated with the systematic errors of inhibitory activities of tyrosinase inhibitors arising from the different measuring protocols. By defining ERR_sys, outliers with larger errors can be efficiently identified and removed from heterogeneous activity data. A reasonable QSAR model (R²_OOB of 0.74 and Q²_CV of 0.80) was obtained by the exclusion of 13 outliers with larger systematic errors. It is a clear example of the challenge for QSAR model that can overwhelm heterogeneous data from different experimental protocols.

Structural ensemble-based docking simulation and biophysical studies discovered new inhibitors of Hsp90 N-terminal domain

Heat shock protein 90 (Hsp90) is one of the most abundant cellular proteins and plays a substantial role in the folding of client proteins. The inhibition of Hsp90 has been regarded as an attractive therapeutic strategy for treating cancer because many oncogenic kinases are Hsp90 client proteins. In this study, we report new inhibitors that directly bind to N-terminal ATP-binding pocket of Hsp90. Optimized structure-based virtual screening predicted candidate molecules, which was followed by confirmation using biophysical and cell-based assays. Among the reported crystal structures, we chose the two structures that show the most favourable early enrichments of true-positives in the receiver operating characteristic curve. Four molecules showed significant changes in the signals of 2D [¹H, ¹⁵N] correlation NMR spectroscopy. Differential scanning calorimetry analysis supported the results indicating direct binding. Quantified dissociation constant values of the molecules, determined by a series of 2D NMR experiments, lie in the range of 0.1–33 μM. Growth inhibition assay with breast and lung cancer cells confirmed the cellular activities of the molecules. Cheminformatics revealed that the molecules share limited chemical similarities with known inhibitors. Molecular dynamics simulations detailed the putative binding modes of the inhibitors.

Understanding the inhibitory mechanism of tea polyphenols against tyrosinase using fluorescence spectroscopy, cyclic voltammetry, oximetry, and molecular simulations

Inhibiting the activity of tyrosinase is a very effective and safe way to prevent enzymatic browning in food and to resist pests in agriculture. Tea polyphenols (TPs), regarded as safe and non-toxic food additives, have been reported due to their potential inhibitory capability against tyrosinase, but their ambiguous inhibitory mechanisms have severely limited their application. In the present work, fluorescence spectroscopy, cyclic voltammetry (CV), oximetry and molecular simulation approaches were employed to shed light on the underlying inhibitory mechanisms of TPs with different structures including (+)-catechin, (−)-epicatechin gallate (ECG) and (−)-epigallocatechin gallate (EGCG) against tyrosinase. Fluorescence spectra show that the three TPs are capable of binding tyrosinase with a molar proportion of 1 : 1. The analysis of CV curves and oxygen utilization suggests that these three TPs can be oxidized by tyrosinase, revealing that these three TPs are suicide inhibitors of tyrosinase. Furthermore, ECG and catechin make tyrosinase irreversibly inactivated due to their catechol group (ring B) being catalyzed by tyrosinase through a cresolase-like pathway, while EGCG inhibits the activity of tyrosinase by competing with or delaying the oxidation of substrate. Molecular simulations further confirm that ring B of ECG and catechin makes a significant contribution to tyrosinase inhibitory activities, and has a direct interaction with the coupled binuclear copper ions in the optimal orientation required by the cresolase-like pathway.

Revealing the inactivation mechanism of TPs on inhibiting the activity of tyrosinase.

Thiopurine Drugs Repositioned as Tyrosinase Inhibitors

Drug repositioning is the application of the existing drugs to new uses and has the potential to reduce the time and cost required for the typical drug discovery process. In this study, we repositioned thiopurine drugs used for the treatment of acute leukaemia as new tyrosinase inhibitors. Tyrosinase catalyses two successive oxidations in melanin biosynthesis: the conversions of tyrosine to dihydroxyphenylalanine (DOPA) and DOPA to dopaquinone. Continuous efforts are underway to discover small molecule inhibitors of tyrosinase for therapeutic and cosmetic purposes. Structure-based virtual screening predicted inhibitor candidates from the US Food and Drug Administration (FDA)-approved drugs. Enzyme assays confirmed the thiopurine leukaemia drug, thioguanine, as a tyrosinase inhibitor with the inhibitory constant of 52 μM. Two other thiopurine drugs, mercaptopurine and azathioprine, were also evaluated for their tyrosinase inhibition; mercaptopurine caused stronger inhibition than thioguanine did, whereas azathioprine was a poor inhibitor. The inhibitory constant of mercaptopurine (16 μM) was comparable to that of the well-known inhibitor kojic acid (13 μM). The cell-based assay using B16F10 melanoma cells confirmed that the compounds inhibit mammalian tyrosinase. Particularly, 50 μM thioguanine reduced the melanin content by 57%, without apparent cytotoxicity. Cheminformatics showed that the thiopurine drugs shared little chemical similarity with the known tyrosinase inhibitors.

Skin whitening agents: medicinal chemistry perspective of tyrosinase inhibitors

Melanogenesis is a process to synthesize melanin, which is a primary responsible for the pigmentation of human skin, eye and hair. Although numerous enzymatic catalyzed and chemical reactions are involved in melanogenesis process, the enzymes such as tyrosinase and tyrosinase-related protein-1 (TRP-1) and TRP-2 played a major role in melanin synthesis. Specifically, tyrosinase is a key enzyme, which catalyzes a rate-limiting step of the melanin synthesis, and the downregulation of tyrosinase is the most prominent approach for the development of melanogenesis inhibitors. Therefore, numerous inhibitors that target tyrosinase have been developed in recent years. The review focuses on the recent discovery of tyrosinase inhibitors that are directly involved in the inhibition of tyrosinase catalytic activity and functionality from all sources, including laboratory synthetic methods, natural products, virtual screening and structure-based molecular docking studies.

Antithrombotic properties of JJ1, a potent and novel thrombin inhibitor

The development of new anticoagulants is an important goal for the improvement of thrombosis treatment. Recent studies have suggested the importance of thrombin inhibitors in the modulation of thromboembolic disorders. The aim of this study was to discover a new small-molecule thrombin inhibitor. In this study, the compound JJ1, which has a novel scaffold, was selected by structure-based docking simulation to determine its potential inhibitory activity against thrombin. JJ1 was shown to inhibit the catalytic activity of human α-thrombin with a K _i of 0.019 μM by direct binding to the active site and with at least 10,000-fold selectivity relative to that reported for the inhibition of other biologically important serine proteases. JJ1 prolonged clotting times (activated partial thromboplastin time and prothrombin time) and inhibited the activity and production of thrombin. Furthermore, it inhibited thrombin-catalyzed fibrin polymerization and platelet aggregation. Similar to its in vitro antithrombotic activities, JJ1 showed enhanced antithrombotic effects in an in vivo pulmonary embolism and arterial thrombosis model. It also exhibited anticoagulant effects in mice. Collectively, these results demonstrated that JJ1 was a potent, direct, and selective thrombin inhibitor that may be useful in the management of various thrombotic disorders.

Determination of the Bridging Ligand in the Active Site of Tyrosinase

Tyrosinase is a type-3 copper enzyme that is widely distributed in plants, fungi, insects, and mammals. Developing high potent inhibitors against tyrosinase is of great interest in diverse fields including tobacco curing, food processing, bio-insecticides development, cosmetic development, and human healthcare-related research. In the crystal structure of Agaricus bisporus mushroom tyrosinase, there is an oxygen atom bridging the two copper ions in the active site. It is unclear whether the identity of this bridging oxygen is a water molecule or a hydroxide anion. In the present study, we theoretically determine the identity of this critical bridging oxygen by performing first-principles hybrid quantum mechanics/molecular mechanics/Poisson-Boltzmann-surface area (QM/MM-PBSA) calculations along with a thermodynamic cycle that aim to improve the accuracy. Our results show that the binding with water molecule is energy favored and the QM/MM-optimized structure is very close to the crystal structure, whereas the binding with hydroxide anions causes the increase of energy and significant structural changes of the active site, indicating that the identity of the bridging oxygen must be a water molecule rather than a hydroxide anion. The different binding behavior between water and hydroxide anions may explain why molecules with a carboxyl group or too many negative charges have lower inhibitory activity. In light of this, the design of high potent active inhibitors against tyrosinase should satisfy both the affinity to the copper ions and the charge neutrality of the entire molecule.

Identification of potent inhibitors of DNA methyltransferase 1 (DNMT1) through a pharmacophore-based virtual screening approach

DNA methylation is an epigenetic change that results in the addition of a methyl group at the carbon-5 position of cytosine residues. DNA methyltransferase (DNMT) inhibitors can suppress tumour growth and have significant therapeutic value. However, the established inhibitors are limited in their application due to their substantial cytotoxicity. Additionally, the standard drugs for DNMT inhibition are non-selective cytosine analogues with considerable cytotoxic side-effects. In the present study, we have designed a workflow by integrating various ligand-based and structure-based approaches to discover new agents active against DNMT1. We have derived a pharmacophore model with the help of available DNMT1 inhibitors. Utilising this model, we performed the virtual screening of Maybridge chemical library and the identified hits were then subsequently filtered based on the Naïve Bayesian classification model. The molecules that have returned from this classification model were subjected to ensemble based docking. We have selected 10 molecules for the biological assay by inspecting the interactions portrayed by these molecules. Three out of the ten tested compounds have shown DNMT1 inhibitory activity. These compounds were also found to demonstrate potential inhibition of cellular proliferation in human breast cancer MDA-MB-231 cells. In the present study, we have utilized a multi-step virtual screening protocol to identify inhibitors of DNMT1, which offers a starting point to develop more potent DNMT1 inhibitors as anti-cancer agents.

The unravelling of the complex pattern of tyrosinase inhibition

Tyrosinases are responsible for melanin formation in all life domains. Tyrosinase inhibitors are used for the prevention of severe skin diseases, in skin-whitening creams and to avoid fruit browning, however continued use of many such inhibitors is considered unsafe. In this study we provide conclusive evidence of the inhibition mechanism of two well studied tyrosinase inhibitors, KA (kojic acid) and HQ (hydroquinone), which are extensively used in hyperpigmentation treatment. KA is reported in the literature with contradicting inhibition mechanisms, while HQ is described as both a tyrosinase inhibitor and a substrate. By visualization of KA and HQ in the active site of TyrBm crystals, together with molecular modeling, binding constant analysis and kinetic experiments, we have elucidated their mechanisms of inhibition, which was ambiguous for both inhibitors. We confirm that while KA acts as a mixed inhibitor, HQ can act both as a TyrBm substrate and as an inhibitor.

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

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