Computer-aided drug design of falcipain inhibitors: virtual screening, structure-activity relationships, hydration site thermodynamics, and reactivity analysis

Falcipains (FPs) are hemoglobinases of Plasmodium falciparum that are validated targets for the development of antimalarial chemotherapy. A combined ligand- and structure-based virtual screening of commercial databases was performed to identify structural analogs of virtual screening hits previously discovered in our laboratory. A total of 28 low micromolar inhibitors of FP-2 and FP-3 were identified and the structure-activity relationship (SAR) in each series was elaborated. The SAR of the compounds was unusually steep in some cases and could not be explained by a traditional analysis of the ligand-protein interactions (van der Waals, electrostatics, and hydrogen bonds). To gain further insights, a statistical thermodynamic analysis of explicit solvent in the ligand binding domains of FP-2 and FP-3 was carried out to understand the roles played by water molecules in binding of these inhibitors. Indeed, the energetics associated with the displacement of water molecules upon ligand binding explained some of the complex trends in the SAR. Furthermore, low potency of a subset of FP-2 inhibitors that could not be understood by the water energetics was explained in the context of poor chemical reactivity of the reactive centers of these compounds. The present study highlights the importance of considering energetic contributors to binding beyond traditional ligand-protein interactions.     

Synthesis and antimalarial activity of novel dihydro-artemisinin derivatives

The Plasmodium falciparum cysteine protease falcipain-2, one of the most promising targets for antimalarial drug design, plays a key role in parasite survival as a major peptide hydrolase within the hemoglobin degradation pathway. In this work, a series of novel dihydroartemisinin derivatives based on (thio)semicarbazone scaffold were designed and synthesized as potential falcipain-2 inhibitors. The in vitro biological assay indicated that most of the target compounds showed excellent inhibition activity against P. falciparum falcipain-2, with IC(50) values in the 0.29-10.63 μM range. Molecular docking studies were performed to investigate the binding affinities and interaction modes for the inhibitors. The preliminary SARs were summarized and could serve as a foundation for further investigation in the development of antimalarial drugs.     

Novel anti-Plasmodial hits identified by virtual screening of the ZINC database

Increased resistance of Plasmodium falciparum to most available drugs challenges the control of malaria. Studies with protease inhibitors have suggested important roles for the falcipain family of cysteine proteases. These enzymes act in concert with other proteases to hydrolyze host erythrocyte hemoglobin in the parasite food vacuole. In order to find potential new antimalarial drugs, we screened in silico the ZINC database using two different protocols involving structure- and ligand-based methodologies. Our search identified 19 novel low micromolar inhibitors of cultured chloroquine resistant P. falciparum. The most active compound presented an IC₅₀ value of 0.5 μM against cultured parasites and it also inhibited the cysteine protease falcipain-2 (IC₅₀ = 25.5 μM). These results identify novel classes of antimalarials that are structurally different from those currently in use and which can be further derivatized to deliver leads suitable for optimisation.     

Toward the discovery of inhibitors of babesipain-1, a Babesia bigemina cysteine protease: in vitro evaluation, homology modeling and molecular docking studies

Babesia bigemina is a protozoan parasite that causes babesiosis, a disease with a world-wide distribution in mammals, principally affecting cattle and man. The unveiling of the genome of B. bigemina is a project in active progress that has already revealed a number of new targets with potential interest for the design of anti-babesiosis drugs. In this context, babesipain-1 has been identified as a proteolytically active enzyme whose three-dimensional structure has not been resolved yet, but which is known to be inhibited by cysteine proteases inhibitors such as E64, ALLN, leupeptin, and vinyl sulfones. In this work, we introduce (1) a homology model of babesipain-1; (2) a comparison between babesipain-1 and falcipain-2, a cysteine protease of the malaria parasite Plasmodium falciparum; (3) in vitro data for babesipain-1 inhibition by HEDICINs and HECINs, previously reported as modest inhibitors of falcipain-2; and (4) the docked binding conformations of HEDICINs and HECINs in the model of babesipain-1. HEDICINs presented similar preferred binding conformations for both babesipain-1 and falcipain-2. However, in vitro bioassay shows that HEDICINs and HECINs are better inhibitors of babesipain-1 than of falcipain-2, which could be explained by observed differences between the active pockets of these proteins in silico. Results presented herein provide a valuable contribution to future computer-aided molecular design of new babesipain-1 inhibitors.     

Design of new plasmepsin inhibitors: a virtual high throughput screening approach on the EGEE grid

Though different species of the genus Plasmodium may be responsible for malaria, the variant caused by P. falciparum is often very dangerous and even fatal if untreated. Hemoglobin degradation is one of the key metabolic processes for the survival of the Plasmodium parasite in its host. Plasmepsins, a family of aspartic proteases encoded by the Plasmodium genome, play a prominent role in host hemoglobin cleavage. In this paper we demonstrate the use of virtual screening, in particular molecular docking, employed at a very large scale to identify novel inhibitors for plasmepsins II and IV. A large grid infrastructure, the EGEE grid, was used to address the problem of large computation resources required for docking hundreds of thousands of chemical compounds on different plasmepsin targets of P. falciparum. A large compound library of about 1 million chemical compounds was docked on 5 different targets of plasmepsins using two different docking software, namely FlexX and AutoDock. Several strategies were employed to analyze the results of this virtual screening approach including docking scores, ideal binding modes, and interactions to key residues of the protein. Three different classes of structures with thiourea, diphenylurea, and guanidino scaffolds were identified to be promising hits. While the identification of diphenylurea compounds is in accordance with the literature and thus provides a sort of "positive control", the identification of novel compounds with a guanidino scaffold proves that high throughput docking can be effectively used to identify novel potential inhibitors of P. falciparum plasmepsins. Thus, with the work presented here, we do not only demonstrate the relevance of computational grids in drug discovery but also identify several promising small molecules which have the potential to serve as candidate inhibitors for P. falciparum plasmepsins. With the use of the EGEE grid infrastructure for the virtual screening campaign against the malaria causing parasite P. falciparum we have demonstrated that resource sharing on an eScience infrastructure such as EGEE provides a new model for doing collaborative research to fight diseases of the poor.     

Substrate activity screening: a fragment-based method for the rapid identification of nonpeptidic protease inhibitors

A new fragment-based method for the rapid development of novel and distinct classes of nonpeptidic protease inhibitors, Substrate Activity Screening (SAS), is described. This method consists of three steps: (1) a library of N-acyl aminocoumarins with diverse, low molecular weight N-acyl groups is screened to identify protease substrates using a simple fluorescence-based assay, (2) the identified N-acyl aminocoumarin substrates are optimized by rapid analogue synthesis and evaluation, and (3) the optimized substrates are converted to inhibitors by direct replacement of the aminocoumarin with known mechanism-based pharmacophores. The SAS method was successfully applied to the cysteine protease cathepsin S, which is implicated in autoimmune diseases. Multiple distinct classes of nonpeptidic substrates were identified upon screening an N-acyl aminocoumarin library. Two of the nonpeptidic substrate classes were optimized to substrates with >8000-fold improvements in cleavage efficiency for each class. Select nonpeptidic substrates were then directly converted to low molecular weight, novel aldehyde inhibitors with nanomolar affinity to cathepsin S. This study demonstrates the unique characteristics and merits of this first substrate-based method for the rapid identification and optimization of weak fragments and provides the framework for the development of completely nonpeptidic inhibitors to many different proteases.     

Tuning and predicting biological affinity: aryl nitriles as cysteine protease inhibitors

A series of aryl nitrile-based ligands were prepared to investigate the effect of their electrophilicity on the affinity against the cysteine proteases rhodesain and human cathepsin L. Density functional theory calculations provided relative reactivities of the nitriles, enabling prediction of their biological affinity and cytotoxicity and a clear structure-activity relationship.     

Docking-based virtual screening studies aiming at the covalent inhibition of SARS-CoV-2 MPro by targeting the cysteine 145

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the COVID-19 which has infected millions of people worldwide. The main protease of SARS-CoV-2 (M^Pro) has been recognized as a key target for the development of antiviral compounds. Taking advantage of the X-ray crystal complex with reversible covalent inhibitors interacting with the catalytic cysteine 145 (Cys145), we explored flexible docking studies to select alternative compounds able to target this residue as covalent inhibitors. First, docking studies of three known electrophilic compounds led to results consistent with co-crystallized data validating the method for SARS-CoV-2 M^Pro covalent inhibition. Then, libraries of soft electrophiles (overall 41 757 compounds) were submitted to docking-based virtual screening resulting in the identification of 17 molecules having their electrophilic group close to the Cys145 residue. We also investigated flexible docking studies of a focused approved covalent drugs library including 32 compounds with various electrophilic functional groups. Among them, the calculations resulted in the identification of four compounds, namely dimethylfumarate, fosfomycin, ibrutinib and saxagliptin, able first, to bind to the active site of the protein and second, to form a covalent bond with the catalytic cysteine.     

Synthesis of a diverse library of mechanism-based cysteine protease inhibitors

We report improvements of our method for the solid-phase synthesis of mechanism-based mercaptomethyl ketone inhibitors of cysteine proteases (Lee, A.; Huang, L.; Ellman, J. A. J. Am. Chem. Soc. 1999, 121, 9907-9914). Specifically, Fmoc-protected chloromethyl ketones were used, rather than the Alloc-protected counterparts. In addition, we further demonstrated that diverse polar functionality can be incorporated at the R1', R1, and R2 sites, in contrast to our previous efforts, where primarily hydrophobic groups were incorporated at these positions. On the basis of these results, a 2016-membered library of potential mercaptomethyl ketone inhibitors was prepared that incorporated diverse functionality. The library was screened against cathepsin B, which is implicated in cancer, resulting in the identification of single-digit nanomolar inhibitors. Because of the diverse functionality incorporated in this library, it should be a rich source of potent inhibitors against many other cysteine proteases.     

Deep Learning and Structure-Based Virtual Screening for Drug Discovery against NEK7: A Novel Target for the Treatment of Cancer

NIMA-related kinase7 (NEK7) plays a multifunctional role in cell division and NLRP3 inflammasone activation. A typical expression or any mutation in the genetic makeup of NEK7 leads to the development of cancer malignancies and fatal inflammatory disease, i.e., breast cancer, non-small cell lung cancer, gout, rheumatoid arthritis, and liver cirrhosis. Therefore, NEK7 is a promising target for drug development against various cancer malignancies. The combination of drug repurposing and structure-based virtual screening of large libraries of compounds has dramatically improved the development of anticancer drugs. The current study focused on the virtual screening of 1200 benzene sulphonamide derivatives retrieved from the PubChem database by selecting and docking validation of the crystal structure of NEK7 protein (PDB ID: 2WQN). The compounds library was subjected to virtual screening using Auto Dock Vina. The binding energies of screened compounds were compared to standard Dabrafenib. In particular, compound 762 exhibited excellent binding energy of −42.67 kJ/mol, better than Dabrafenib (−33.89 kJ/mol). Selected drug candidates showed a reactive profile that was comparable to standard Dabrafenib. To characterize the stability of protein–ligand complexes, molecular dynamic simulations were performed, providing insight into the molecular interactions. The NEK7–Dabrafenib complex showed stability throughout the simulated trajectory. In addition, binding affinities, pIC50, and ADMET profiles of drug candidates were predicted using deep learning models. Deep learning models predicted the binding affinity of compound 762 best among all derivatives, which supports the findings of virtual screening. These findings suggest that top hits can serve as potential inhibitors of NEK7. Moreover, it is recommended to explore the inhibitory potential of identified hits compounds through in-vitro and in-vivo approaches.     

Design, synthesis, and evaluation of in vivo potency and selectivity of epoxysuccinyl-based inhibitors of papain-family cysteine proteases

The papain-family cathepsins are cysteine proteases that are emerging as promising therapeutic targets for a number of human disease conditions ranging from osteoporosis to cancer. Relatively few selective inhibitors for this family exist, and the in vivo selectivity of most existing compounds is unclear. We present here the synthesis of focused libraries of epoxysuccinyl-based inhibitors and their screening in crude tissue extracts. We identified a number of potent inhibitors that display selectivity for endogenous cathepsin targets both in vitro and in vivo. Importantly, the selectivity patterns observed in crude extracts were generally retained in vivo, as assessed by active-site labeling of tissues from treated animals. Overall, this study identifies several important compound classes and highlights the use of activity-based probes to assess pharmacodynamic properties of small-molecule inhibitors in vivo.     

Organic azide inhibitors of cysteine proteases

Cysteine proteases are crucial regulatory enzymes in human physiology and disease. Inhibitors are usually designed with reactive electrophiles to covalently bond to the catalytic cysteinyl sulfur, and consequently they also indiscriminately interact with biological thiolates and other nucleophiles, leading to toxic side effects in vivo. Here we describe an alternative to using reactive electrophiles, demonstrating the use of a much less reactive azidomethylene substituent (-CH2-N3) that confers potent inhibition of cysteine proteases. This new approach resulted in potent, reversible, competitive inhibitors of caspase-1 (IC50 < 10 nM), with significant advantages over aldehydes such as high stability in vitro to thiols (10 mM dithiothreitol (pH 7.2), 20 mM glutathione (pH 7.2, 9, 11)) and aqueous media, as well as some highly desirable druglike features. It was also demonstrated that azides can be incorporated into inhibitors of other caspases (e.g. 3, 8) and cathepsins (e.g. K, S, B), indicating the versatility of this valuable new approach to cysteine protease inhibition.     

A peptide aldehyde microarray for high-throughput profiling of cellular events

Microarrays provide exciting opportunities in the field of large-scale proteomics. With the aim to elucidate enzymatic activity and profiles within native biological samples, we developed a microarray comprising a focused positional-scanning library of enzyme inhibitors. The library was diversified across P(1)-P(4) positions, creating 270 different inhibitor sublibraries which were immobilized onto avidin slides. The peptide aldehyde-based small-molecule microarray (SMM) specifically targeted cysteine proteases, thereby enabling large-scale functional assessment of this subgroup of proteases, within fluorescently labeled samples, including pure proteins, cellular lysates, and infected samples. The arrays were shown to elicit binding fingerprints consistent with those of model proteins, specifically caspases and purified cysteine proteases from parasites (rhodesein and cruzain). When tested against lysates from apoptotic Hela and red blood cells infected with Plasmodium falciparum, clear signatures were obtained that were readily attributable to the activity of constituent proteases within these samples. Characteristic binding profiles were further able to distinguish various stages of the parasite infection in erythrocyte lysates. By converting one of our brightest microarray hits into a probe, putative protein markers were identified and pulled down from within apoptotic Hela lysates, demonstrating the potential of target validation and discovery. Taken together, these results demonstrate the utility of targeted SMMs in dissecting cellular biology in complex proteomic samples.     

Identification of novel PfDHODH inhibitors as antimalarial agents via pharmacophore-based virtual screening followed by molecular docking and in vivo antimalarial activity

Plasmodium falciparum dihydroorotate dehydrogenase (PfDHODH) catalyses the fourth reaction of de novo pyrimidine biosynthesis in parasites, and represents an important target for the treatment of malaria. In this study, we describe pharmacophore-based virtual screening combined with docking study and biological evaluation as a rational strategy for identification of novel hits as antimalarial agents. Pharmacophore models were established from known PfDHODH inhibitors using the GALAHAD module with IC₅₀ values ranging from 0.033 μM to 142 μM. The best pharmacophore model consisted of three hydrogen bond acceptor, one hydrogen bond donor and one hydrophobic features. The pharmacophore models were validated through receiver operating characteristic and Günere–Henry scoring methods. The best pharmacophore model as a 3D search query was searched against the IBS database. Several compounds with different structures (scaffolds) were retrieved as hit molecules. Among these compounds, those with a Q_FIT value of more than 81 were docked in the PfDHODH enzyme to further explore the binding modes of these compounds. In silico pharmacokinetic and toxicities were predicted for the best docked molecules. Finally, the identified hits were evaluated in vivo for their antimalarial activity in a parasite inhibition assay. The hits reported here showed good potential to become novel antimalarial agents.     

The role of in vitro ADME assays in antimalarial drug discovery and development

The high level of attrition of drugs in clinical development has led pharmaceutical companies to increase the efficiency of their lead identification and development through techniques such as combinatorial chemistry and high-throughput (HTP) screening. Since the major reasons for clinical drug candidate failure other than efficacy are pharmacokinetics and toxicity, attention has been focused on assessing properties such as metabolic stability, drug-drug interactions (DDI), and absorption earlier in the drug discovery process. Animal studies are simply too labor-intensive and expensive to use for evaluating every hit, so it has been necessary to develop and implement higher throughput in vitro ADME screens to manage the large number of compounds of interest. The antimalarial drug development program at the Walter Reed Army Institute of Research, Division of Experimental Therapeutics (WRAIR/ET) has adopted this paradigm in its search for a long-term prophylactic for the prevention of malaria. The overarching goal of this program is to develop new, long half-life, orally bioavailable compounds with potent intrinsic activity against liver- and blood-stage parasites. From the WRAIR HTP antimalarial screen, numerous compounds are regularly identified with potent activity. These hits are now immediately evaluated using a panel of in vitro ADME screens to identify and predict compounds that will meet our specific treatment criteria. In this review, the WRAIR ADME screening program for antimalarial drugs is described as well as how we have implemented it to predict the ADME properties of small molecule for the identification of promising drug candidates.     

In silico identification of noncompetitive inhibitors targeting an uncharacterized allosteric site of falcipain-2

Falcipain-2 (FP-2) is a Plasmodium falciparum hemoglobinase widely targeted in the search for antimalarials. FP-2 can be allosterically modulated by various noncompetitive inhibitors that have been serendipitously identified. Moreover, the crystal structures of two inhibitors bound to an allosteric site, termed site 6, of the homolog enzyme human cathepsin K (hCatK) suggest that the equivalent region in FP-2 might play a similar role. Here, we conduct the rational identification of FP-2 inhibitors through virtual screenings (VS) of compounds into several pocket-like conformations of site 6, sampled during molecular dynamics (MD) simulations of the free enzyme. Two noncompetitive inhibitors, ZINC03225317 and ZINC72290660, were confirmed using in vitro enzymatic assays and their poses into site 6 led to calculated binding free energies matching the experimental ones. Our results provide strong evidence about the allosteric inhibition of FP-2 through binding of small molecules to site 6, thus opening the way toward the discovery of new inhibitors against this enzyme.

     

Inhibitors of SARS-3CLpro: virtual screening, biological evaluation, and molecular dynamics simulation studies

SARS-CoV from the coronaviridae family has been identified as the etiological agent of Severe Acute Respiratory Syndrome (SARS), a highly contagious upper respiratory disease that reached epidemic status in 2002. SARS-3CL(pro), a cysteine protease indispensible to the viral life cycle, has been identified as one of the key therapeutic targets against SARS. A combined ligand and structure-based virtual screening was carried out against the Asinex Platinum collection. Multiple low micromolar inhibitors of the enzyme were identified through this search, one of which also showed activity against SARS-CoV in a whole cell CPE assay. Furthermore, multinanosecond explicit solvent simulations were carried out using the docking poses of the identified hits to study the overall stability of the binding site interactions as well as identify important changes in the interaction profile that were not apparent from the docking study. Cumulative analysis of the evaluated compounds and the simulation studies led to the identification of certain protein-ligand interaction patterns which would be useful in further structure based design efforts.     

Identification of potential Tpx inhibitors against pathogen-host interactions

Yersinia organisms cause many infectious diseases by invading human cells and delivering their virulence factors via the type three secretion system (T3SS). One alternative strategy in the fight against these pathogenic organisms is to interfere with their T3SS. Previous studies demonstrated that thiol peroxidase, Tpx is functional in the assembly of T3SS and its inhibition by salicylidene acylhydrazides prevents the secretion of pathogenic effectors. In this study, the aim was to identify potential inhibitors of Tpx using an integrated approach starting with high throughput virtual screening and ending with molecular dynamics simulations of selected ligands. Virtual screening of ZINC database of 500,000 compounds via ligand-based and structure-based pharmacophore models retrieved 10,000 hits. The structure-based pharmacophore model was validated using high-throughput virtual screening (HTVS). After multistep docking (SP and XP), common scaffolds were used to find common substructures and the ligand binding poses were optimized using induced fit docking. The stability of the protein–ligand complex was examined with molecular dynamics simulations and the binding free energy of the complex was calculated. As a final outcome eight compounds with different chemotypes were proposed as potential inhibitors for Tpx. The eight ligands identified by a detailed virtual screening protocol can serve as leads in future drug design efforts against the destructive actions of pathogenic bacteria.