Development of a Fingerprint-Based Scoring Function for the Prediction of the Binding Mode of Carbonic Anhydrase II Inhibitors

Identification of natural xanthine oxidase inhibitors: Virtual screening, anti-xanthine oxidase activity, and interaction mechanism

Xanthine oxidase (XO) is a crucial target for hyperuricemia treatment(s). Naturally occurred XO inhibitors with minimal toxicity and high efficacy have attracted researchers' attention. With the goal of quickly identifying natural XO inhibitors, an integrated computational screening strategy was constructed by molecular docking and calculating the free energy of binding. Twenty-seven hits were achieved from a database containing 19,377 natural molecules. This includes fourteen known XO inhibitors and four firstly-reported inhibitors (isolicoflavonol, 5,7-dihydroxycoumarin, parvifolol D and clauszoline M, IC50 < 40 μM). Iolicoflavonol (hit 8, IC50 = 8.45 ± 0.68 μM) and 5,7-dihydroxycoumarin (hit 25, IC50 = 10.91 ± 0.71 μM) displayed the great potency as mixed-type inhibitors. Docking study and molecular dynamics simulation revealed that both hits could interact with XO's primarily active site residues ARG880, MOS1328, and ASN768 of XO. Fluorescence spectroscopy studies showed that hit 8 bound to the active cavity region of XO, causing changes in XO's conformation and hydrophobicity. Hits 8 and 25 exhibit favorable Absorption, Distribution, Metabolism, and Excretion (ADME) properties. Additionally, no cytotoxicity against human liver cells was observed at their median inhibition concentrations against XO. Therefore, the present study offers isolicoflavonol and 5,7-dihydroxycoumarin with the potential to be disease-modifying agents for hyperuricemia.

Prediction of activity and selectivity profiles of human Carbonic Anhydrase inhibitors using machine learning classification models

The development of selective inhibitors of the clinically relevant human Carbonic Anhydrase (hCA) isoforms IX and XII has become a major topic in drug research, due to their deregulation in several types of cancer. Indeed, the selective inhibition of these two isoforms, especially with respect to the homeostatic isoform II, holds great promise to develop anticancer drugs with limited side effects. Therefore, the development of in silico models able to predict the activity and selectivity against the desired isoform(s) is of central interest. In this work, we have developed a series of machine learning classification models, trained on high confidence data extracted from ChEMBL, able to predict the activity and selectivity profiles of ligands for human Carbonic Anhydrase isoforms II, IX and XII. The training datasets were built with a procedure that made use of flexible bioactivity thresholds to obtain well-balanced active and inactive classes. We used multiple algorithms and sampling sizes to finally select activity models able to classify active or inactive molecules with excellent performances. Remarkably, the results herein reported turned out to be better than those obtained by models built with the classic approach of selecting an a priori activity threshold. The sequential application of such validated models enables virtual screening to be performed in a fast and more reliable way to predict the activity and selectivity profiles against the investigated isoforms.

Consensus Docking in Drug Discovery

Background:

Molecular docking is probably the most popular and profitable approach in computer-aided drug design, being the staple technique for predicting the binding mode of bioactive compounds and for performing receptor-based virtual screening studies. The growing attention received by docking, as well as the need for improving its reliability in pose prediction and virtual screening performance, has led to the development of a wide plethora of new docking algorithms and scoring functions. Nevertheless, it is unlikely to identify a single procedure outperforming the other ones in terms of reliability and accuracy or demonstrating to be generally suitable for all kinds of protein targets.

Methods:

In this context, consensus docking approaches are taking hold in computer-aided drug design. These computational protocols consist in docking ligands using multiple docking methods and then comparing the binding poses predicted for the same ligand by the different methods. This analysis is usually carried out calculating the root-mean-square deviation among the different docking results obtained for each ligand, in order to identify the number of docking methods producing the same binding pose.

Results:

The consensus docking approaches demonstrated to improve the quality of docking and virtual screening results compared to the single docking methods. From a qualitative point of view, the improvement in pose prediction accuracy was obtained by prioritizing ligand binding poses produced by a high number of docking methods, whereas with regards to virtual screening studies, high hit rates were obtained by prioritizing the compounds showing a high level of pose consensus.

Conclusion:

In this review, we provide an overview of the results obtained from the performance assessment of various consensus docking protocols and we illustrate successful case studies where consensus docking has been applied in virtual screening studies.

Application of MM-PBSA Methods in Virtual Screening

Computer-aided drug design techniques are today largely applied in medicinal chemistry. In particular, receptor-based virtual screening (VS) studies, in which molecular docking represents the gold standard in silico approach, constitute a powerful strategy for identifying novel hit compounds active against the desired target receptor. Nevertheless, the need for improving the ability of docking in discriminating true active ligands from inactive compounds, thus boosting VS hit rates, is still pressing. In this context, the use of binding free energy evaluation approaches can represent a profitable tool for rescoring ligand-protein complexes predicted by docking based on more reliable estimations of ligand-protein binding affinities than those obtained with simple scoring functions. In the present review, we focused our attention on the Molecular Mechanics-Poisson Boltzman Surface Area (MM-PBSA) method for the calculation of binding free energies and its application in VS studies. We provided examples of successful applications of this method in VS campaigns and evaluation studies in which the reliability of this approach has been assessed, thus providing useful guidelines for employing this approach in VS.

Development of a cheminformatics platform for selectivity analyses of carbonic anhydrase inhibitors

The selectivity for a specific human Carbonic Anhydrase (hCA) isoform is an important property a hCA inhibitor (CAI) should be endowed with, in order to constitute a valuable therapeutic tool for the treatment of a desired pathology. In this context, we developed a chemoinformatic platform that allows the analysis of the structure and selectivity profile of known CAIs reported in literature, with the aim of identifying structural motifs connected to ligand selectivity, thus providing useful guidelines for the design of novel ligands selective for the desired hCA isoform. The platform is able to perform ultrafast structure and selectivity analyses through ligand fingerprint similarity, with no need of structural information about the target receptor and ligands' binding mode. It is easily accessible to the non-expert user through the implementation of a KNIME Analytic Platform workflow and could be extended to analyze the selectivity profile of known ligands of different target proteins.

Virtual screening identifies a PIN1 inhibitor with possible antiovarian cancer effects

Peptidyl-prolyl cis-trans isomerase, NIMA-interacting 1 (PIN1) is a peptidyl-prolyl isomerase that binds phospho-Ser/Thr-Pro motifs in proteins and catalyzes the cis-trans isomerization of proline peptide bonds. PIN1 is overexpressed in several cancers including high-grade serous ovarian cancer. Since few therapies are effective against this cancer, PIN1 could be a therapeutic target but effective PIN1 inhibitors are lacking. To identify molecules with in vivo inhibitory effects on PIN1, we used consensus docking to model existing PIN1-ligand X-ray structures and to screen a chemical database for candidate inhibitors. Ten molecules were selected and tested in cellular assays, leading to the identification of VS10 that bound and inhibited PIN1. VS10 treatment reduced the viability of ovarian cancer cell lines by inducing proteasomal PIN1 degradation, without effects on PIN1 transcription, and also reduced the levels of downstream targets β-catenin, cyclin D1, and pSer473-Akt. VS10 is a selective PIN1 inhibitor that may offer new opportunities for treating PIN1-overexpressing tumors.

An Overview of Scoring Functions Used for Protein-Ligand Interactions in Molecular Docking

Currently, molecular docking is becoming a key tool in drug discovery and molecular modeling applications. The reliability of molecular docking depends on the accuracy of the adopted scoring function, which can guide and determine the ligand poses when thousands of possible poses of ligand are generated. The scoring function can be used to determine the binding mode and site of a ligand, predict binding affinity and identify the potential drug leads for a given protein target. Despite intensive research over the years, accurate and rapid prediction of protein-ligand interactions is still a challenge in molecular docking. For this reason, this study reviews four basic types of scoring functions, physics-based, empirical, knowledge-based, and machine learning-based scoring functions, based on an up-to-date classification scheme. We not only discuss the foundations of the four types scoring functions, suitable application areas and shortcomings, but also discuss challenges and potential future study directions.

Novel 8-Substituted Coumarins That Selectively Inhibit Human Carbonic Anhydrase IX and XII

A novel series of 8-substituted coumarin-based compounds, characterized by the presence of alkylpiperazine and arylpiperazine chains, were synthesized and tested for their inhibitory activity against four human carbonic anhydrase (hCA) isoforms. All compounds displayed nanomolar potency against the cancer-related hCA IX and hCA XII; moreover, they were shown to be devoid of any inhibitory activity toward the cytosolic hCA I and hCA II up to 10 µM concentration in the assay system. Therefore, the synthesized coumarin ligands demonstrated to be potent and selective hCA IX/XII inhibitors, and were shown to be as potent as the reference inhibitor acetazolamide against hCA XII, with single-digit nanomolar K_i values. Molecular modeling studies provided a rationale for explaining the selectivity profile of these non-classic hCA inhibitors and their interactions with the enzymes, according to their specific mechanism of action, thus paving the way for future structure-based lead optimization studies.

Computationally driven discovery of phenyl(piperazin-1-yl)methanone derivatives as reversible monoacylglycerol lipase (MAGL) inhibitors

Monoacylglycerol lipase (MAGL) is an attractive therapeutic target for many pathologies, including neurodegenerative diseases, cancer as well as chronic pain and inflammatory pathologies. The identification of reversible MAGL inhibitors, devoid of the side effects associated to prolonged MAGL inactivation, is a hot topic in medicinal chemistry. In this study, a novel phenyl(piperazin-1-yl)methanone inhibitor of MAGL was identified through a virtual screening protocol based on a fingerprint-driven consensus docking (CD) approach. Molecular modeling and preliminary structure-based hit optimization studies allowed the discovery of derivative 4, which showed an efficient reversible MAGL inhibition (IC₅₀ = 6.1 µM) and a promising antiproliferative activity on breast and ovarian cancer cell lines (IC₅₀ of 31-72 µM), thus representing a lead for the development of new and more potent reversible MAGL inhibitors. Moreover, the obtained results confirmed the reliability of the fingerprint-driven CD approach herein developed.

Target-based drug discovery through inversion of quantitative structure-drug-property relationships and molecular simulation: CA IX-sulphonamide complexes

In this work, a target-based drug screening method is proposed exploiting the synergy effect of ligand-based and structure-based computer-assisted drug design. The new method provides great flexibility in drug design and drug candidates with considerably lower risk in an efficient manner. As a model system, 45 sulphonamides (33 training, 12 testing ligands) in complex with carbonic anhydrase IX were used for development of quantitative structure-activity-lipophilicity (property)-relationships (QSPRs). For each ligand, nearly 5,000 molecular descriptors were calculated, while lipophilicity (logk_w) and inhibitory activity (logK_i) were used as drug properties. Genetic algorithm-partial least squares (GA-PLS) provided a QSPR model with high prediction capability employing only seven molecular descriptors. As a proof-of-concept, optimal drug structure was obtained by inverting the model with respect to reference drug properties. 3509 ligands were ranked accordingly. Top 10 ligands were further validated through molecular docking. Large-scale MD simulations were performed to test the stability of structures of selected ligands obtained through docking complemented with biophysical experiments.