In silico screening of LRRK2 WDR domain inhibitors using deep docking and free energy simulations.

The Critical Assessment of Computational Hit-Finding Experiments (CACHE) Challenge series is focused on identifying small molecule inhibitors of protein targets using computational methods. Each challenge contains two phases, hit-finding and follow-up optimization, each of which is followed by experimental validation of the computational predictions. For the CACHE Challenge #1, the Leucine-Rich Repeat Kinase 2 (LRRK2) WD40 Repeat (WDR) domain was selected as the target for in silico hit-finding and optimization. Mutations in LRRK2 are the most common genetic cause of the familial form of Parkinson's disease. The LRRK2 WDR domain is an understudied drug target with no known molecular inhibitors. Herein we detail the first phase of our winning submission to the CACHE Challenge #1. We developed a framework for the high-throughput structure-based virtual screening of a chemically diverse small molecule space. Hit identification was performed using the large-scale Deep Docking (DD) protocol followed by absolute binding free energy (ABFE) simulations. ABFEs were computed using an automated molecular dynamics (MD)-based thermodynamic integration (TI) approach. 4.1 billion ligands from Enamine REAL were screened with DD followed by ABFEs computed by MD TI for 793 ligands. 76 ligands were prioritized for experimental validation, with 59 compounds successfully synthesized and 5 compounds identified as hits, yielding a 8.5% hit rate. Our results demonstrate the efficacy of the combined DD and ABFE approaches for hit identification for a target with no previously known hits. This approach is widely applicable for the efficient screening of ultra-large chemical libraries as well as rigorous protein-ligand binding affinity estimation leveraging modern computational resources.

PROTACable Is an Integrative Computational Pipeline of 3-D Modeling and Deep Learning To Automate the De Novo Design of PROTACs.

Proteolysis-targeting chimeras (PROTACs) that engage two biological targets at once are a promising technology in degrading clinically relevant protein targets. Since factors that influence the biological activities of PROTACs are more complex than those of a small molecule drug, we explored a combination of computational chemistry and deep learning strategies to forecast PROTAC activity and enable automated design. A new method named PROTACable was developed for the de novo design of PROTACs, which includes a robust 3-D modeling workflow to model PROTAC ternary complexes using a library of E3 ligase and linker and an SE(3)-equivariant graph transformer network to predict the activity of newly designed PROTACs. PROTACable is available at https://github.com/giaguaro/PROTACable/.

A novel class of broad-spectrum active-site-directed 3C-like protease inhibitors with nanomolar antiviral activity against highly immune-evasive SARS-CoV-2 Omicron subvariants.

Antivirals with broad coronavirus activity are important for treating high-risk individuals exposed to the constantly evolving SARS-CoV-2 variants of concern (VOCs) as well as emerging drug-resistant variants. We developed and characterized a novel class of active-site-directed 3-chymotrypsin-like protease (3CLpro) inhibitors (C2-C5a). Our lead direct-acting antiviral (DAA), C5a, is a non-covalent, non-peptide with a dissociation constant of 170 nM against recombinant SARS-CoV-2 3CLpro. The compounds C2-C5a exhibit broad-spectrum activity against Omicron subvariants (BA.5, BQ.1.1, and XBB.1.5) and seasonal human coronavirus-229E infection in human cells. Notably, C5a has median effective concentrations of 30-50 nM against BQ.1.1 and XBB.1.5 in two different human cell lines. X-ray crystallography has confirmed the unique binding modes of C2-C5a to the 3CLpro, which can limit virus cross-resistance to emerging Paxlovid-resistant variants. We tested the effect of C5a with two of our newly discovered host-directed antivirals (HDAs): N-0385, a TMPRSS2 inhibitor, and bafilomycin D (BafD), a human vacuolar H+-ATPase [V-ATPase] inhibitor. We demonstrated a synergistic action of C5a in combination with N-0385 and BafD against Omicron BA.5 infection in human Calu-3 lung cells. Our findings underscore that a SARS-CoV-2 multi-targeted treatment for circulating Omicron subvariants based on DAAs (C5a) and HDAs (N-0385 or BafD) can lead to therapeutic benefits by enhancing treatment efficacy. Furthermore, the high-resolution structures of SARS-CoV-2 3CLpro in complex with C2-C5a will facilitate future rational optimization of our novel broad-spectrum active-site-directed 3C-like protease inhibitors.

Novel Inhibitors of androgen receptor's DNA binding domain identified using an ultra-large virtual screening.

Androgen receptor (AR) inhibition remains the primary strategy to combat the progression of prostate cancer (PC). However, all clinically used AR inhibitors target the ligand-binding domain (LBD), which is highly susceptible to truncations through splicing or mutations that confer drug resistance. Thus, there exists an urgent need for AR inhibitors with novel modes of action. We thus launched a virtual screening of an ultra-large chemical library to find novel inhibitors of the AR DNA-binding domain (DBD) at two sites: protein-DNA interface (P-box) and dimerization site (D-box). The compounds selected through vigorous computational filtering were then experimentally validated. We identified several novel chemotypes that effectively suppress transcriptional activity of AR and its splice variant V7. The identified compounds represent previously unexplored chemical scaffolds with a mechanism of action that evades the conventional drug resistance manifested through LBD mutations. Additionally, we describe the binding features required to inhibit AR DBD at both P-box and D-box target sites.

Large-Scale Virtual Screening for the Discovery of SARS-CoV-2 Papain-like Protease (PLpro) Non-covalent Inhibitors.

The rapid global spread of the SARS-CoV-2 virus facilitated the development of novel direct-acting antiviral agents (DAAs). The papain-like protease (PLpro) has been proposed as one of the major SARS-CoV-2 targets for DAAs due to its dual role in processing viral proteins and facilitating the host's immune suppression. This dual role makes identifying small molecules that can effectively neutralize SARS-CoV-2 PLpro activity a high-priority task. However, PLpro drug discovery faces a significant challenge due to the high mobility and induced-fit effects in the protease's active site. Herein, we virtually screened the ZINC20 database with Deep Docking (DD) to identify prospective noncovalent PLpro binders and combined ultra-large consensus docking with two pharmacophore (ph4)-filtering strategies. The analysis of active compounds revealed their somewhat-limited diversity, likely attributed to the induced-fit nature of PLpro's active site in the crystal structures, and therefore, the use of rigid docking protocols poses inherited limitations. The top hits were assessed against recombinant viral proteins and live viruses, demonstrating desirable inhibitory activities. The best compound VPC-300195 (IC50: 15 µM) ranks among the top noncovalent PLpro inhibitors discovered through in silico methodologies. In the search for novel SARS-CoV-2 PLpro-specific chemotypes, the identified inhibitors could serve as diverse templates for the development of effective noncovalent PLpro inhibitors.

Targeting SARS-CoV-2 papain-like protease in the postvaccine era.

While vaccines remain at the forefront of global healthcare responses, pioneering therapeutics against SARS-CoV-2 are expected to fill the gaps for waning immunity. Rapid development and approval of orally available direct-acting antivirals targeting crucial SARS-CoV-2 proteins marked the beginning of the era of small-molecule drugs for COVID-19. In that regard, the papain-like protease (PLpro) can be considered a major SARS-CoV-2 therapeutic target due to its dual biological role in suppressing host innate immune responses and in ensuring viral replication. Here, we summarize the challenges of targeting PLpro and innovative early-stage PLpro-specific small molecules. We propose that state-of-the-art computer-aided drug design (CADD) methodologies will play a critical role in the discovery of PLpro compounds as a novel class of COVID-19 drugs.

Structure-Based Study to Overcome Cross-Reactivity of Novel Androgen Receptor Inhibitors.

The mutation-driven transformation of clinical anti-androgen drugs into agonists of the human androgen receptor (AR) represents a major challenge for the treatment of prostate cancer patients. To address this challenge, we have developed a novel class of inhibitors targeting the DNA-binding domain (DBD) of the receptor, which is distanced from the androgen binding site (ABS) targeted by all conventional anti-AR drugs and prone to resistant mutations. While many members of the developed 4-(4-phenylthiazol-2-yl)morpholine series of AR-DBD inhibitors demonstrated the effective suppression of wild-type AR, a few represented by 4-(4-(3-fluoro-2-methoxyphenyl)thiazol-2-yl)morpholine (VPC14368) exhibited a partial agonistic effect toward the mutated T878A form of the receptor, implying their cross-interaction with the AR ABS. To study the molecular basis of the observed cross-reactivity, we co-crystallized the T878A mutated form of the AR ligand binding domain (LBD) with a bound VPC14368 molecule. Computational modelling revealed that helix 12 of AR undergoes a characteristic shift upon VPC14368 binding causing the agonistic behaviour. Based on the obtained structural data we then designed derivatives of VPC14368 to successfully eliminate the cross-reactivity towards the AR ABS, while maintaining significant anti-AR DBD potency.

Development of VPC-70619, a Small-Molecule N-Myc Inhibitor as a Potential Therapy for Neuroendocrine Prostate Cancer.

The Myc family of transcription factors are involved in the development and progression of numerous cancers, including prostate cancer (PCa). Under the pressure of androgen receptor (AR)-directed therapies resistance can occur, leading to the lethal form of PCa known as neuroendocrine prostate cancer (NEPC), characterized among other features by N-Myc overexpression. There are no clinically approved treatments for NEPC, translating into poor patient prognosis and survival. Therefore, there is a pressing need to develop novel therapeutic avenues to treat NEPC patients. In this study, we investigate the N-Myc-Max DNA binding domain (DBD) as a potential target for small molecule inhibitors and utilize computer-aided drug design (CADD) approaches to discover prospective hits. Through further exploration and optimization, a compound, VPC-70619, was identified with notable anti-N-Myc potency and strong antiproliferative activity against numerous N-Myc expressing cell lines, including those representing NEPC.

Artificial intelligence-enabled virtual screening of ultra-large chemical libraries with deep docking.

With the recent explosion of chemical libraries beyond a billion molecules, more efficient virtual screening approaches are needed. The Deep Docking (DD) platform enables up to 100-fold acceleration of structure-based virtual screening by docking only a subset of a chemical library, iteratively synchronized with a ligand-based prediction of the remaining docking scores. This method results in hundreds- to thousands-fold virtual hit enrichment (without significant loss of potential drug candidates) and hence enables the screening of billion molecule-sized chemical libraries without using extraordinary computational resources. Herein, we present and discuss the generalized DD protocol that has been proven successful in various computer-aided drug discovery (CADD) campaigns and can be applied in conjunction with any conventional docking program. The protocol encompasses eight consecutive stages: molecular library preparation, receptor preparation, random sampling of a library, ligand preparation, molecular docking, model training, model inference and the residual docking. The standard DD workflow enables iterative application of stages 3-7 with continuous augmentation of the training set, and the number of such iterations can be adjusted by the user. A predefined recall value allows for control of the percentage of top-scoring molecules that are retained by DD and can be adjusted to control the library size reduction. The procedure takes 1-2 weeks (depending on the available resources) and can be completely automated on computing clusters managed by job schedulers. This open-source protocol, at https://github.com/jamesgleave/DD_protocol , can be readily deployed by CADD researchers and can significantly accelerate the effective exploration of ultra-large portions of a chemical space.

Development of 2-(5,6,7-Trifluoro-1H-Indol-3-yl)-quinoline-5-carboxamide as a Potent, Selective, and Orally Available Inhibitor of Human Androgen Receptor Targeting Its Binding Function-3 for the Treatment of Castration-Resistant Prostate Cancer.

Prostate cancer (PCa) patients undergoing androgen deprivation therapy almost invariably develop castration-resistant prostate cancer (CRPC). Targeting the androgen receptor (AR) Binding Function-3 (BF3) site offers a promising option to treat CRPC. However, BF3 inhibitors have been limited by poor potency or inadequate metabolic stability. Through extensive medicinal chemistry, molecular modeling, and biochemistry, we identified 2-(5,6,7-trifluoro-1H-Indol-3-yl)-quinoline-5-carboxamide (VPC-13789), a potent AR BF3 antagonist with markedly improved pharmacokinetic properties. We demonstrate that VPC-13789 suppresses AR-mediated transcription, chromatin binding, and recruitment of coregulatory proteins. This novel AR antagonist selectively reduces the growth of both androgen-dependent and enzalutamide-resistant PCa cell lines. Having demonstrated in vitro efficacy, we developed an orally bioavailable prodrug that reduced PSA production and tumor volume in animal models of CRPC with no observed toxicity. VPC-13789 is a potent, selective, and orally bioavailable antiandrogen with a distinct mode of action that has a potential as novel CRPC therapeutics.

Development of an Androgen Receptor Inhibitor Targeting the N-Terminal Domain of Androgen Receptor for Treatment of Castration Resistant Prostate Cancer.

Prostate cancer patients undergoing androgen deprivation therapy almost invariably develop castration-resistant prostate cancer. Resistance can occur when mutations in the androgen receptor (AR) render anti-androgen drugs ineffective or through the expression of constitutively active splice variants lacking the androgen binding domain entirely (e.g., ARV7). In this study, we are reporting the discovery of a novel AR-NTD covalent inhibitor 1-chloro-3-[(5-([(2S)-3-chloro-2-hydroxypropyl]amino)naphthalen-1-yl)amino]propan-2-ol (VPC-220010) targeting the AR-N-terminal Domain (AR-NTD). VPC-220010 inhibits AR-mediated transcription of full length and truncated variant ARV7, downregulates AR response genes, and selectively reduces the growth of both full-length AR- and truncated AR-dependent prostate cancer cell lines. We show that VPC-220010 disrupts interactions between AR and known coactivators and coregulatory proteins, such as CHD4, FOXA1, ZMIZ1, and several SWI/SNF complex proteins. Taken together, our data suggest that VPC-220010 is a promising small molecule that can be further optimized into effective AR-NTD inhibitor for the treatment of CRPC.

Drugging the 'undruggable'. Therapeutic targeting of protein-DNA interactions with the use of computer-aided drug discovery methods.

Transcription factors (TFs) act as major oncodrivers in many cancers and are frequently regarded as high-value therapeutic targets. The functionality of TFs relies on direct protein-DNA interactions, which are notoriously difficult to target with small molecules. However, this prior view of the 'undruggability' of protein-DNA interfaces has shifted substantially in recent years, in part because of significant advances in computer-aided drug discovery (CADD). In this review, we highlight recent examples of successful CADD campaigns resulting in drug candidates that directly interfere with protein-DNA interactions of several key cancer TFs, including androgen receptor (AR), ETS-related gene (ERG), MYC, thymocyte selection-associated high mobility group box protein (TOX), topoisomerase II (TOP2), and signal transducer and activator of transcription 3 (STAT3). Importantly, these findings open novel and compelling avenues for therapeutic targeting of over 1600 human TFs implicated in many conditions including and beyond cancer.

Development of Novel Inhibitors Targeting the D-Box of the DNA Binding Domain of Androgen Receptor.

The inhibition of the androgen receptor (AR) is an established strategy in prostate cancer (PCa) treatment until drug resistance develops either through mutations in the ligand-binding domain (LBD) portion of the receptor or its deletion. We previously identified a druggable pocket on the DNA binding domain (DBD) dimerization surface of the AR and reported several potent inhibitors that effectively disrupted DBD-DBD interactions and consequently demonstrated certain antineoplastic activity. Here we describe further development of small molecule inhibitors of AR DBD dimerization and provide their broad biological characterization. The developed compounds demonstrate improved activity in the mammalian two-hybrid assay, enhanced inhibition of AR-V7 transcriptional activity, and improved microsomal stability. These findings position us for the development of AR inhibitors with entirely novel mechanisms of action that would bypass most forms of PCa treatment resistance, including the truncation of the LBD of the AR.

Automated discovery of noncovalent inhibitors of SARS-CoV-2 main protease by consensus Deep Docking of 40 billion small molecules.

Recent explosive growth of 'make-on-demand' chemical libraries brought unprecedented opportunities but also significant challenges to the field of computer-aided drug discovery. To address this expansion of the accessible chemical universe, molecular docking needs to accurately rank billions of chemical structures, calling for the development of automated hit-selecting protocols to minimize human intervention and error. Herein, we report the development of an artificial intelligence-driven virtual screening pipeline that utilizes Deep Docking with Autodock GPU, Glide SP, FRED, ICM and QuickVina2 programs to screen 40 billion molecules against SARS-CoV-2 main protease (Mpro). This campaign returned a significant number of experimentally confirmed inhibitors of Mpro enzyme, and also enabled to benchmark the performance of twenty-eight hit-selecting strategies of various degrees of stringency and automation. These findings provide new starting scaffolds for hit-to-lead optimization campaigns against Mpro and encourage the development of fully automated end-to-end drug discovery protocols integrating machine learning and human expertise.

Crystallographic structure of wild-type SARS-CoV-2 main protease acyl-enzyme intermediate with physiological C-terminal autoprocessing site.

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), the pathogen that causes the disease COVID-19, produces replicase polyproteins 1a and 1ab that contain, respectively, 11 or 16 nonstructural proteins (nsp). Nsp5 is the main protease (Mpro) responsible for cleavage at eleven positions along these polyproteins, including at its own N- and C-terminal boundaries, representing essential processing events for subsequent viral assembly and maturation. We have determined X-ray crystallographic structures of this cysteine protease in its wild-type free active site state at 1.8 Å resolution, in its acyl-enzyme intermediate state with the native C-terminal autocleavage sequence at 1.95 Å resolution and in its product bound state at 2.0 Å resolution by employing an active site mutation (C145A). We characterize the stereochemical features of the acyl-enzyme intermediate including critical hydrogen bonding distances underlying catalysis in the Cys/His dyad and oxyanion hole. We also identify a highly ordered water molecule in a position compatible for a role as the deacylating nucleophile in the catalytic mechanism and characterize the binding groove conformational changes and dimerization interface that occur upon formation of the acyl-enzyme. Collectively, these crystallographic snapshots provide valuable mechanistic and structural insights for future antiviral therapeutic development including revised molecular docking strategies based on Mpro inhibition.

Bioanalytical Assay Development and Validation for the Pharmacokinetic Study of GMC1, a Novel FKBP52 Co-chaperone Inhibitor for Castration Resistant Prostate Cancer.

BACKGROUND: GMC1 (2-(1H-benzimidazol-2-ylsulfanyl)-N-[(Z)-(4-methoxyphenyl) methylideneamino] acetamide) effectively inhibits androgen receptor function by binding directly to FKBP52. This is a novel mechanism for the treatment of castration resistant prostate cancer (CRPC). METHODS: an LC-MS/MS method was developed and validated to quantify GMC1 in plasma and urine from pharmacokinetics studies in rats. An ultra-high-performance liquid chromatography (UHPLC) system equipped with a Waters XTerra MS C18 column was used for chromatographic separation by gradient elution with 0.1% (v/v) formic acid in water and methanol. A Sciex 4000 QTRAP® mass spectrometer was used for analysis by multiple reaction monitoring (MRM) in positive mode; the specific ions [M+H]+m/z 340.995 → m/z 191.000 and [M+H]+ m/z 266.013 → m/z 234.000 were monitored for GMC1 and internal standard (albendazole), respectively. RESULTS: GMC1 and albendazole had retention times of 1.68 and 1.66 min, respectively. The calibration curves for the determination of GMC1 in rat plasma and urine were linear from 1-1000 ng/mL. The LC-MS/MS method was validated with intra- and inter-day accuracy and precision within the 15% acceptance limit. The extraction recovery values of GMC1 from rat plasma and urine were greater than 95.0 ± 2.1% and 97.6 ± 4.6%, respectively, with no significant interfering matrix effect. GMC1 is stable under expected sample handling, storage, preparation and LC-MS/MS analysis conditions. CONCLUSIONS: Pharmacokinetic evaluation of GMC1 revealed that the molecule has a biexponential disposition in rats, is distributed rapidly and extensively, has a long elimination half-life, and appears to be eliminated primarily by first order kinetics.

Dual-Inhibitors of N-Myc and AURKA as Potential Therapy for Neuroendocrine Prostate Cancer.

Resistance to androgen-receptor (AR) directed therapies is, among other factors, associated with Myc transcription factors that are involved in development and progression of many cancers. Overexpression of N-Myc protein in prostate cancer (PCa) leads to its transformation to advanced neuroendocrine prostate cancer (NEPC) that currently has no approved treatments. N-Myc has a short half-life but acts as an NEPC stimulator when it is stabilized by forming a protective complex with Aurora A kinase (AURKA). Therefore, dual-inhibition of N-Myc and AURKA would be an attractive therapeutic avenue for NEPC. Following our computer-aided drug discovery approach, compounds exhibiting potent N-Myc specific inhibition and strong anti-proliferative activity against several N-Myc driven cell lines, were identified. Thereafter, we have developed dual inhibitors of N-Myc and AURKA through structure-based drug design approach by merging our novel N-Myc specific chemical scaffolds with fragments of known AURKA inhibitors. Favorable binding modes of the designed compounds to both N-Myc and AURKA target sites have been predicted by docking. A promising lead compound, 70812, demonstrated low-micromolar potency against both N-Myc and AURKA in vitro assays and effectively suppressed NEPC cell growth.

Interleukin-10 and Small Molecule SHIP1 Allosteric Regulators Trigger Anti-inflammatory Effects through SHIP1/STAT3 Complexes.

The anti-inflammatory actions of interleukin-10 (IL10) are thought to be mediated primarily by the STAT3 transcription factor, but pro-inflammatory cytokines such as interleukin-6 (IL6) also act through STAT3. We now report that IL10, but not IL6 signaling, induces formation of a complex between STAT3 and the inositol polyphosphate-5-phosphatase SHIP1 in macrophages. Both SHIP1 and STAT3 translocate to the nucleus in macrophages. Remarkably, sesquiterpenes of the Pelorol family, which we previously described as allosteric activators of SHIP1 phosphatase activity, could induce SHIP1/STAT3 complex formation in cells and mimic the anti-inflammatory action of IL10 in a mouse model of colitis. Using crystallography and docking studies we identified a drug-binding pocket in SHIP1. Our studies reveal new mechanisms of action for both STAT3 and SHIP1 and provide a rationale for use of allosteric SHIP1-activating compounds, which mimic the beneficial anti-inflammatory actions of IL10. VIDEO ABSTRACT.

Deep Docking: A Deep Learning Platform for Augmentation of Structure Based Drug Discovery.

Drug discovery is a rigorous process that requires billion dollars of investments and decades of research to bring a molecule "from bench to a bedside". While virtual docking can significantly accelerate the process of drug discovery, it ultimately lags the current rate of expansion of chemical databases that already exceed billions of molecular records. This recent surge of small molecules availability presents great drug discovery opportunities, but also demands much faster screening protocols. In order to address this challenge, we herein introduce Deep Docking (DD), a novel deep learning platform that is suitable for docking billions of molecular structures in a rapid, yet accurate fashion. The DD approach utilizes quantitative structure-activity relationship (QSAR) deep models trained on docking scores of subsets of a chemical library to approximate the docking outcome for yet unprocessed entries and, therefore, to remove unfavorable molecules in an iterative manner. The use of DD methodology in conjunction with the FRED docking program allowed rapid and accurate calculation of docking scores for 1.36 billion molecules from the ZINC15 library against 12 prominent target proteins and demonstrated up to 100-fold data reduction and 6000-fold enrichment of high scoring molecules (without notable loss of favorably docked entities). The DD protocol can readily be used in conjunction with any docking program and was made publicly available.

Computer-Aided Ligand Discovery for Estrogen Receptor Alpha.

Breast cancer (BCa) is one of the most predominantly diagnosed cancers in women. Notably, 70% of BCa diagnoses are Estrogen Receptor α positive (ERα+) making it a critical therapeutic target. With that, the two subtypes of ER, ERα and ERß, have contrasting effects on BCa cells. While ERα promotes cancerous activities, ERß isoform exhibits inhibitory effects on the same. ER-directed small molecule drug discovery for BCa has provided the FDA approved drugs tamoxifen, toremifene, raloxifene and fulvestrant that all bind to the estrogen binding site of the receptor. These ER-directed inhibitors are non-selective in nature and may eventually induce resistance in BCa cells as well as increase the risk of endometrial cancer development. Thus, there is an urgent need to develop novel drugs with alternative ERα targeting mechanisms that can overcome the limitations of conventional anti-ERα therapies. Several functional sites on ERα, such as Activation Function-2 (AF2), DNA binding domain (DBD), and F-domain, have been recently considered as potential targets in the context of drug research and discovery. In this review, we summarize methods of computer-aided drug design (CADD) that have been employed to analyze and explore potential targetable sites on ERα, discuss recent advancement of ERα inhibitor development, and highlight the potential opportunities and challenges of future ERα-directed drug discovery.