Metadynamics for Perspective Drug Design: Computationally Driven Synthesis of New Protein-Protein Interaction Inhibitors Targeting the EphA2 Receptor

Research status of indole-modified natural products

Indole is a heterocyclic compound formed by the fusion of a benzene ring and pyrrole ring, which has rich biological activity. Many indole-containing compounds have been sold on the market due to their excellent pharmacological activity. For example, vincristine and reserpine have been widely used in clinical practice. The diverse structures and biological activities of natural products provide abundant resources for the development of new drugs. Therefore, this review classifies natural products by structure, and summarizes the research progress of indole-containing natural product derivatives, their biological activities, structure-activity relationship and research mechanism which has been studied in the past 13 years, so as to provide a basis for the development of new drug development.

Palladium-Catalyzed Synthesis of Linked Bis-Heterocycles─Synthesis and Investigation of Photophysical Properties

A palladium-catalyzed domino C-N coupling/Cacchi reaction is reported. Design of photoluminescent bis-heterocycles, aided by density functional theory calculations, was performed with synthetic yields up to 98%. The photophysical properties of the products accessed via this strategy were part of a comprehensive study that led to broad emission spectra and quantum yields of up to 0.59. Mechanistic experiments confirmed bromoalkynes as competent intermediates, and a density functional theory investigation suggests a pathway involving initial oxidative addition into the cis C-Br bond of the gem-dihaloolefin.

Molecular Determinants of EphA2 and EphB2 Antagonism Enable the Design of Ligands with Improved Selectivity

With the aim of identifying novel antagonists selective for the EphA receptor family, a combined experimental and computational approach was taken to investigate the molecular basis of the recognition between a prototypical Eph-ephrin antagonist (UniPR1447) and two representative receptors of the EphA and EphB subfamilies, namely, EphA2 and EphB2 receptors. The conformational free-energy surface (FES) of the binding state of UniPR1447 within the ligand binding domain of EphA2 and EphB2, reconstructed from molecular dynamics (MD) simulations performed on the microsecond time scale, was exploited to drive the design and synthesis of a novel antagonist selective for EphA2 over the EphB2 receptor. The availability of compounds with this pharmacological profile will help discriminate the importance of these two receptors in the insurgence and progression of cancer.

Unbinding Kinetics of Muscarinic M3 Receptor Antagonists Explained by Metadynamics Simulations

The residence time (RT), the time for which a drug remains bound to its biological target, is a critical parameter for drug design. The prediction of this key kinetic property has been proven to be challenging and computationally demanding in the framework of atomistic simulations. In the present work, we setup and applied two distinct metadynamics protocols to estimate the RTs of muscarinic M3 receptor antagonists. In the first method, derived from the conformational flooding approach, the kinetics of unbinding is retrieved from a physics-based parameter known as the acceleration factor α (i.e., the running average over time of the potential deposited in the bound state). Such an approach is expected to recover the absolute RT value for a compound of interest. In the second method, known as the tMETA-D approach, a qualitative estimation of the RT is given by the time of simulation required to drive the ligand from the binding site to the solvent bulk. This approach has been developed to reproduce the change of experimental RTs for compounds targeting the same target. Our analysis shows that both computational protocols are able to rank compounds in agreement with their experimental RTs. Quantitative structure-kinetics relationship (SKR) models can be identified and employed to predict the impact of a chemical modification on the experimental RT once a calibration study has been performed.

Expanding the Accessible Chemical Space of SIRT2 Inhibitors through Exploration of Binding Pocket Dynamics

Considerations of binding pocket dynamics are one of the crucial aspects of the rational design of binders. Identification of alternative conformational states or cryptic subpockets could lead to the discovery of completely novel groups of the ligands. However, experimental characterization of pocket dynamics, besides being expensive, may not be able to elucidate all of the conformational states relevant for drug discovery projects. In this study, we propose the protocol for computational simulations of sirtuin 2 (SIRT2) binding pocket dynamics and its integration into the structure-based virtual screening (SBVS) pipeline. Initially, unbiased molecular dynamics simulations of SIRT2:inhibitor complexes were performed using optimized force field parameters of SIRT2 inhibitors. Time-lagged independent component analysis (tICA) was used to design pocket-related collective variables (CVs) for enhanced sampling of SIRT2 pocket dynamics. Metadynamics simulations in the tICA eigenvector space revealed alternative conformational states of the SIRT2 binding pocket and the existence of a cryptic subpocket. Newly identified SIRT2 conformational states outperformed experimentally resolved states in retrospective SBVS validation. After performing prospective SBVS, compounds from the under-represented portions of the SIRT2 inhibitor chemical space were selected for in vitro evaluation. Two compounds, NDJ18 and NDJ85, were identified as potent and selective SIRT2 inhibitors, which validated the in silico protocol and opened up the possibility for generalization and broadening of its application. The anticancer effects of the most potent compound NDJ18 were examined on the triple-negative breast cancer cell line. Results indicated that NDJ18 represents a promising structure suitable for further evaluation.

Protein-Protein Interaction Inhibitors Targeting the Eph-Ephrin System with a Focus on Amino Acid Conjugates of Bile Acids

The role of the Eph-ephrin system in the etiology of pathological conditions has been consolidated throughout the years. In this context, approaches directed against this signaling system, intended to modulate its activity, can be strategic therapeutic opportunities. Currently, the most promising class of compounds able to interfere with the Eph receptor-ephrin protein interaction is composed of synthetic derivatives of bile acids. In the present review, we summarize the progresses achieved, in terms of chemical expansions and structure-activity relationships, both in the steroidal core and the terminal carboxylic acid group, along with the pharmacological characterization for the most promising Eph-ephrin antagonists in in vivo settings.

Metabolic Soft Spot and Pharmacokinetics: Functionalization of C-3 Position of an Eph-Ephrin Antagonist Featuring a Bile Acid Core as an Effective Strategy to Obtain Oral Bioavailability in Mice

UniPR129, an L-β-homotryptophan conjugate of the secondary bile acid lithocholic acid (LCA), acts as an effective protein-protein interaction (PPI) inhibitor of the Eph-ephrin system but suffers from a poor oral bioavailability in mice. To improve UniPR129 bioavailability, a metabolic soft spot, i.e., the 3α-hydroxyl group on the LCA steroidal ring, was functionalized to 3-hydroxyimine. In vitro metabolism of UniPR129 and 3-hydroxyimine derivative UniPR500 was compared in mouse liver subcellular fractions, and main metabolites were profiled by high resolution (HR-MS) and tandem (MS/MS) mass spectrometry. In mouse liver microsomes (MLM), UniPR129 was converted into several metabolites: M1 derived from the oxidation of the 3-hydroxy group to 3-oxo, M2-M7, mono-hydroxylated metabolites, M8-M10, di-hydroxylated metabolites, and M11, a mono-hydroxylated metabolite of M1. Phase II reactions were only minor routes of in vitro biotransformation. UniPR500 shared several metabolic pathways with parent UniPR129, but it showed higher stability in MLM, with a half-life (t_1/2) of 60.4 min, if compared to a t_1/2 = 16.8 min for UniPR129. When orally administered to mice at the same dose, UniPR500 showed an increased systemic exposure, maintaining an in vitro valuable pharmacological profile as an EphA2 receptor antagonist and an overall improvement in its physico-chemical profile (solubility, lipophilicity), if compared to UniPR129. The present work highlights an effective strategy for the pharmacokinetic optimization of aminoacid conjugates of bile acids as small molecule Eph-ephrin antagonists.

A higher flexibility at the SARS-CoV-2 main protease active site compared to SARS-CoV and its potentialities for new inhibitor virtual screening targeting multi-conformers

The main-protease (Mpro) catalyzes a crucial step for the SARS-CoV-2 life cycle. The recent SARS-CoV-2 presents the main protease (MCoV2pro) with 12 mutations compared to SARS-CoV (MCoV1pro). Recent studies point out that these subtle differences lead to mobility variances at the active site loops with functional implications. We use metadynamics simulations and a sort of computational analysis to probe the dynamic, pharmacophoric and catalytic environment differences between the monomers of both enzymes. So, we verify how much intrinsic distinctions are preserved in the functional dimer of MCoV2pro, as well as its implications for ligand accessibility and optimized drug screening. We find a significantly higher accessibility to open binding conformers in the MCoV2pro monomer compared to MCoV1pro. A higher hydration propensity for the MCoV2pro S2 loop with the A46S substitution seems to exercise a key role. Quantum calculations suggest that the wider conformations for MCoV2pro are less catalytically active in the monomer. However, the statistics for contacts involving the N-finger suggest higher maintenance of this activity at the dimer. Docking analyses suggest that the ability to vary the active site width can be important to improve the access of the ligand to the active site in different ways. So, we carry out a multiconformational virtual screening with different ligand bases. The results point to the importance of taking into account the protein conformational multiplicity for new promissors anti MCoV2pro ligands. We hope these results will be useful in prospecting, repurposing and/or designing new anti SARS-CoV-2 drugs.Communicated by Ramaswamy H. Sarma.

Chimeric Antigen Receptor-modified T cells targeting EphA2 for the immunotherapy of paediatric bone tumours

Chimeric Antigen Receptor (CAR) T-cell therapy, as an approved treatment option for patients with B cell malignancies, demonstrates that genetic modification of autologous immune cells is an effective anti-cancer regimen. Erythropoietin-producing Hepatocellular receptor tyrosine kinase class A2 (EphA2) is a tumour associated antigen expressed on a range of sarcomas, including paediatric osteosarcoma (OS) and Ewing sarcoma (ES). We tested human EphA2 directed CAR T cells for their capacity to target and kill human OS and ES tumour cells using in vitro and in vivo assays, demonstrating that EphA2 CAR T cells have potent anti-tumour efficacy in vitro and can eliminate established OS and ES tumours in vivo in a dose and delivery route dependent manner. Next, in an aggressive metastatic OS model we demonstrated that systemically infused EphA2 CAR T cells can traffic to and eradicate tumour deposits in murine livers and lungs. These results support further pre-clinical evaluation of EphA2 CAR T cells to inform the design of early phase clinical trial protocols to test the feasibility and safety of this immune cell therapy in paediatric bone sarcoma patients.

Evaluation of the Anti-Tumor Activity of Small Molecules Targeting Eph/Ephrins in APC min/J Mice

The Eph receptors are the largest receptors tyrosine kinases (RTKs) family in humans and together with ephrin ligands constitute a complex cellular communication system often dysregulated in many tumors. The role of the Eph-ephrin system in colorectal cancer (CRC) has been investigated and different expression of Eph receptors have been associated with tumor development and progression. In light of this evidence, we investigated if a pharmacological approach aimed at inhibiting Eph/ephrin interaction through small molecules could prevent tumor growth in APC ^min/J mice. The 8-week treatment with the Eph-ephrin antagonist UniPR129 significantly reduced the number of adenomas in the ileum and decreased the diameter of adenomas in the same region. Overall our data suggested as UniPR129 could be able to slow down the tumor development in APC ^min/J mice. These results further confirm literature data about Eph kinases as a new valuable target in the intestinal cancer and for the first time showed the feasibility of the Eph-ephrin inhibition as a useful pharmacological approach against the intestinal tumorigenesis. In conclusion this work paves the way for further studies with Eph-ephrin inhibitors in order to confirm the Eph antagonism as innovative pharmacological approach with preventive benefit in the intestinal tumor development.

Ephrin or not? Six tough questions on Eph targeting

Introduction: The Eph-ephrin is a cell-cell communication system generating a forward signal in cell expressing Eph receptors and a reverse signal in ephrin-ligand expressing cells. While clearly involved in the insurgence and progression of cancer, the understanding of the molecular mechanisms regulated by this system needs development; this is a hurdle to the development of therapeutic strategies that can target the Eph receptors and/or their ephrin ligands.Areas covered: We have taken the opportunity to share some key questions on the most effective strategies to target the Eph-ephrin system. This article is based on our experience of the field and therefore is a Perspective and not comprehensive examination of the literature.Expert opinion: Targeting of the Eph-ephrin system has emerged as a potentially valuable approach for cancer therapy. Pharmacological tools have been reported in the last 15 years and these include forward signaling blockers such as kinases inhibitors and antagonists of forward and reverse signaling. Also, biologics including antibodies and recombinant proteins have been developed and some have reached early clinical stages. Data deem the Eph-ephrin system as a signaling axis that is an elusive target. A better understanding of the basic pharmacology behind the activity of available agents and a comprehensive knowledge of the ephrin biology are necessary. We are looking forward to knowing the opinion of the readers.

Exploring Ligand Stability in Protein Crystal Structures Using Binding Pose Metadynamics

Identification of correct protein–ligand binding poses is important in structure-based drug design and crucial for the evaluation of protein–ligand binding affinity. Protein–ligand coordinates are commonly obtained from crystallography experiments that provide a static model of an ensemble of conformations. Binding pose metadynamics (BPMD) is an enhanced sampling method that allows for an efficient assessment of ligand stability in solution. Ligand poses that are unstable under the bias of the metadynamics simulation are expected to be infrequently occupied in the energy landscape, thus making minimal contributions to the binding affinity. Here, the robustness of the method is studied using crystal structures with ligands known to be incorrectly modeled, as well as 63 structurally diverse crystal structures with ligand fit to electron density from the Twilight database. Results show that BPMD can successfully differentiate compounds whose binding pose is not supported by the electron density from those with well-defined electron density.

Optimization of EphA2 antagonists based on a lithocholic acid core led to the identification of UniPR505, a new 3α-carbamoyloxy derivative with antiangiogenetic properties

The EphA2 receptor has been validated in animal models as new target for treating tumors depending on angiogenesis and vasculogenic mimicry. In the present work, we extended our current knowledge on structure-activity relationship (SAR) data of two related classes of antagonists of the EphA2 receptor, namely 5β-cholan-24-oic acids and 5β-cholan-24-oyl l-β-homotryptophan conjugates, with the aim to develop new antiangiogenic compounds able to efficiently prevent the formation of blood vessels. As a result of our exploration, we identified UniPR505, N-[3α-(Ethylcarbamoyl)oxy-5β-cholan-24-oyl]-l-β-homo-tryptophan (compound 14), as a submicromolar antagonist of the EphA2 receptor capable to block EphA2 phosphorylation and to inhibit neovascularization in a chorioallantoic membrane (CAM) assay.

Drug-gut microbiota metabolic interactions: the case of UniPR1331, selective antagonist of the Eph-ephrin system, in mice

The interest on the role of gut microbiota in the biotransformation of drugs and xenobiotics has grown over the last decades and a deeper understanding of the mutual interactions is expected to help future improvements in the fields of drug development, toxicological risk assessment and precision medicine. In this paper, a microbiome drug metabolism case is presented, involving a lipophilic small molecule, N-(3β-hydroxy-Δ⁵-cholen-24-oyl)-l-tryptophan, UniPR1331, active as antagonist of the Eph-ephrin system and effective in vivo in a murine orthotopic model of glioblastoma multiforme (GBM). Following the administration of a single 30 mg/kg dose (p.o.) to mice, maximal plasma levels were reached 30 min after dosing and rapidly declined thereafter. To explain the observed in vivo behaviour, in vitro phase I and II metabolism assays were conducted employing mouse and human liver subcellular fractions and profiling main metabolites by means of tandem (HPLC-ESI-MS/MS) and high resolution mass spectrometry (HPLC-ESI-HR-MS). In the presence of in vitro mouse liver fractions, UniPR1331 showed a low phase I metabolic clearance, despite the identification of a 3-oxo and several hydroxylated metabolites. Conversely, after oral administration of UniPR1331 to mice, a novel isobaric metabolite was detected that (i) was subjected, as parent UniPR1331, to enterohepatic circulation (ii) had not been previously identified in vitro in mouse liver microsomes and (iii) was not observed forming after intraperitoneal (i.p.) administration of UniPR1331. An in vitro faecal fermentation assay produced the same chemical entity supporting a major role of gut microbiota in the in vivo clearance of UniPR1331.

Investigating the stability of dengue virus envelope protein dimer using well-tempered metadynamics simulations

We explored the stability of the dengue virus envelope (E) protein dimer since it is widely assumed that the E protein dimer is stabilized by drug ligands or antibodies in an acidic environment, neutralizing the virus's ability to fuse with human cells. During this process, a large conformational change of the E protein dimer is required. We performed Molecular Dynamics simulations to mimic the conformational change and stability of the dimer in neutral and acidic conditions with the well-tempered metadynamics method. Furthermore, as a few neutralizing antibodies discovered from dengue patients were reported, we used the same simulation method to examine the influence of a selected antibody on the dimer stability in both neutral and acidic conditions. We also investigated the antibody's influence on a point-mutated E protein that had been reported to interrupt the protein-antibody interaction and result in more than 95% loss of the antibody's binding ability. Our simulation results are highly consistent with the experimental conclusion that binding of the antibody to the E protein dimer neutralizes the virus, especially in a low pH condition, while the mutation of W101A or N153A significantly reduces the antibody's ability in stabilizing the E protein dimer. We demonstrate that well-tempered metadynamics can be used to accurately explore the antibody's interaction on large protein complexes such as the E protein dimer, and the computational approach in this work is promising in future antibody development.

Ligand-Binding Calculations with Metadynamics

All-atom molecular dynamics simulations can capture the dynamic degrees of freedom that characterize molecular recognition, the knowledge of which constitutes the cornerstone of rational approaches to drug design and optimization. In particular, enhanced sampling algorithms, such as metadynamics, are powerful tools to dramatically reduce the computational cost required for a mechanistic description of the binding process. Here, we describe the essential details characterizing these simulation strategies, focusing on the critical step of identifying suitable reaction coordinates, as well as on the different analysis algorithms to estimate binding affinity and residence times. We conclude with a survey of published applications that provides explicit examples of successful simulations for several targets.

Interplay of cysteine exposure and global protein dynamics in small-molecule recognition by a regulator of G-protein signaling protein

Regulator of G protein signaling (RGS) proteins play a pivotal role in regulation of G protein-coupled receptor (GPCR) signaling and are therefore becoming an increasingly important therapeutic target. Recently discovered thiadiazolidinone (TDZD) compounds that target cysteine residues have shown different levels of specificities and potencies for the RGS4 protein, thereby suggesting intrinsic differences in dynamics of this protein upon binding of these compounds. In this work, we investigated using atomistic molecular dynamics (MD) simulations the effect of binding of several small-molecule inhibitors on perturbations and dynamical motions in RGS4. Specifically, we studied two conformational models of RGS4 in which a buried cysteine residue is solvent-exposed due to side-chain motions or due to flexibility in neighboring helices. We found that TDZD compounds with aromatic functional groups perturb the RGS4 structure more than compounds with aliphatic functional groups. Moreover, small-molecules with aromatic functional groups but lacking sulfur atoms only transiently reside within the protein and spontaneously dissociate to the solvent. We further measured inhibitory effects of TDZD compounds using a protein-protein interaction assay on a single-cysteine RGS4 protein showing trends in potencies of compounds consistent with our simulation studies. Thermodynamic analyses of RGS4 conformations in the apo-state and on binding to TDZD compounds revealed links between both conformational models of RGS4. The exposure of cysteine side-chains appears to facilitate initial binding of TDZD compounds followed by migration of the compound into a bundle of four helices, thereby causing allosteric perturbations in the RGS/Gα protein-protein interface.

Pharmacological evaluation of new bioavailable small molecules targeting Eph/ephrin interaction

Eph/ephrin system is an emerging target for cancer therapy but the lack of potent, stable and orally bioavailable compounds is impairing the development of the field. Since 2009 our research group has been devoted to the discovery and development of small molecules targeting Eph/ephrin system and our research culminated with the synthesis of UniPR129, a potent but problematic Eph/ephrin antagonist. Herein, we describe the in vitro pharmacological properties of two derivatives (UniPR139 and UniPR502) stemmed from structure of UniPR129. These two compounds acted as competitive and reversible antagonists of all Eph receptors reducing both ephrin-A1 and -B1 binding to EphAs and EphBs receptors in the low micromolar range. The compounds acted as antagonists inhibiting ephrin-A1-dependent EphA2 activation and UniPR139 exerted an anti-angiogenic effect, inhibiting HUVEC tube formation in vitro and VEGF-induced vessel formation in the chick chorioallantoic membrane assay. Finally, the oral bioavailability of UniPR139 represents a step forward in the search of molecules targeting the Eph/ephrin system and offers a new pharmacological tool useful for future in vivo studies.

Targeting Eph/ephrin system in cancer therapy

It is well established that the Eph/ephrin system plays a central role in the embryonic development, with minor implications in the physiology of the adult. However, it is overexpressed and deregulated in a variety of tumors, with a primary involvement in tumorigenesis, tumor angiogenesis, metastasis development, and cancer stem cell regeneration. Targeting the Eph/ephrin system with biologicals, including antibodies and recombinant proteins, reduces tumor growth in animal models of hematological malignancies, breast, prostate, colon, head and neck cancers and glioblastoma. Currently, some of these biopharmaceutical agents are under investigations in phase I or phase II clinical trials. Peptides and small molecules targeting protein-protein-interaction (PPI) are in the late preclinical phase where they are showing promising activity in models of glioblastoma, ovarian and lung cancer. The present review summarizes the most critical findings proposing the Eph/ephrin signaling system as a new target in molecularly targeted oncology.