Spatial chemical distance based on atomic property fields

A two-stage computational approach to predict novel ligands for a chemosensory receptor

Olfactory receptor (OR) 1A2 is the member of largest superfamily of G protein-coupled receptors (GPCRs). OR1A2 is an ectopically expressed receptor with only 13 known ligands, implicated in reducing hepatocellular carcinoma progression, with enormous therapeutic potential. We have developed a two-stage screening approach to identify novel putative ligands of OR1A2. We first used a pharmacophore model based on atomic property field (APF) to virtually screen a library of 5942 human metabolites. We then carried out structure-based virtual screening (SBVS) for predicting the potential agonists, based on a 3D homology model of OR1A2. This model was developed using a biophysical approach for template selection, based on multiple parameters including hydrophobicity correspondence, applied to the complete set of available GPCR structures to pick the most appropriate template. Finally, the membrane-embedded 3D model was refined by molecular dynamics (MD) simulations in both the apo and holo forms. The refined model in the apo form was selected for SBVS. Four novel small molecules were identified as strong binders to this olfactory receptor on the basis of computed binding energies.

2D- and 3D-QSAR modelling, molecular docking and in vitro evaluation studies on 18β-glycyrrhetinic acid derivatives against triple-negative breast cancer cell line

Triple-negative breast cancers (TNBCs) are one of the most aggressive and complex forms of cancers in women. TNBCs are commonly known for their complex heterogeneity and poor prognosis. The present work aimed to develop a predictive 2D and 3D quantitative structure-activity relationship (QSAR) models against metastatic TNBC cell line. The 2D-QSAR was based on multiple linear regression analysis and validated by Leave-One-Out (LOO) and external test set prediction approach. QSAR model presented regression coefficient values for training set (r²), LOO-based internal regression (q²) and external test set regression (pred_r²) which are 0.84, 0.82 and 0.75, respectively. Five properties, Epsilon4 (electronegativity), ChiV3cluster (valence molecular connectivity index), chi3chain (retention index for three-membered ring), TNN5 (nitrogen atoms separated through 5 bond distance) and nitrogen counts, were identified as important structural features responsible for anticancer activity of MDA-MB-231 inhibitors. Five novel derivatives of glycyrrhetinic acid (GA) named GA-1, GA-2, GA-3, GA-4 and GA-5 were semi-synthesised and screened through the QSAR model. Further, in vitro activities of the derivatives were analysed against human TNBC cell line, MDA-MB-231. The result showed that GA-1 exhibits improved cytotoxic activity to that of parent compound (GA). Further, atomic property field (APF)-based 3D-QSAR and scoring recognise C-30 carboxylic group of GA-1 as major influential factor for its anticancer activity. The significance of C-30 carboxylic group in GA derivatives was also confirmed by molecular docking study against cancer target glyoxalase-I. Finally, the oral bioavailability and toxicity of GA-1 were assessed by computational ADMET studies.Communicated by Ramaswamy H. Sarma.

Hybrid receptor structure/ligand-based docking and activity prediction in ICM: development and evaluation in D3R Grand Challenge 3

In context of D3R Grand Challenge 3 we have investigated several ligand activity prediction protocols that combined elements of a physics-based energy function (ICM VLS score) and the knowledge-based Atomic Property Field 3D QSAR approach. Activity prediction models utilized poses produced by ICM-Dock with ligand bias and 4D receptor conformational ensembles (LigBEnD). Hybrid APF/P (APF/Physics) models were superior to pure physics- or knowledge-based models in our preliminary tests using rigorous three-fold clustered cross-validation and later proved successful in the blind prediction for D3R GC3 sets, consistently performing well across four different targets. The results demonstrate that knowledge-based and physics-based inputs into the machine-learning activity model can be non-redundant and synergistic.

Exploring G Protein-Coupled Receptors (GPCRs) Ligand Space via Cheminformatics Approaches: Impact on Rational Drug Design

The primary goal of rational drug discovery is the identification of selective ligands which act on single or multiple drug targets to achieve the desired clinical outcome through the exploration of total chemical space. To identify such desired compounds, computational approaches are necessary in predicting their drug-like properties. G Protein-Coupled Receptors (GPCRs) represent one of the largest and most important integral membrane protein families. These receptors serve as increasingly attractive drug targets due to their relevance in the treatment of various diseases, such as inflammatory disorders, metabolic imbalances, cardiac disorders, cancer, monogenic disorders, etc. In the last decade, multitudes of three-dimensional (3D) structures were solved for diverse GPCRs, thus referring to this period as the "golden age for GPCR structural biology." Moreover, accumulation of data about the chemical properties of GPCR ligands has garnered much interest toward the exploration of GPCR chemical space. Due to the steady increase in the structural, ligand, and functional data of GPCRs, several cheminformatics approaches have been implemented in its drug discovery pipeline. In this review, we mainly focus on the cheminformatics-based paradigms in GPCR drug discovery. We provide a comprehensive view on the ligand- and structure-based cheminformatics approaches which are best illustrated via GPCR case studies. Furthermore, an appropriate combination of ligand-based knowledge with structure-based ones, i.e., integrated approach, which is emerging as a promising strategy for cheminformatics-based GPCR drug design is also discussed.

Ligand-biased ensemble receptor docking (LigBEnD): a hybrid ligand/receptor structure-based approach

Ligand docking to flexible protein molecules can be efficiently carried out through ensemble docking to multiple protein conformations, either from experimental X-ray structures or from in silico simulations. The success of ensemble docking often requires the careful selection of complementary protein conformations, through docking and scoring of known co-crystallized ligands. False positives, in which a ligand in a wrong pose achieves a better docking score than that of native pose, arise as additional protein conformations are added. In the current study, we developed a new ligand-biased ensemble receptor docking method and composite scoring function which combine the use of ligand-based atomic property field (APF) method with receptor structure-based docking. This method helps us to correctly dock 30 out of 36 ligands presented by the D3R docking challenge. For the six mis-docked ligands, the cognate receptor structures prove to be too different from the 40 available experimental Pocketome conformations used for docking and could be identified only by receptor sampling beyond experimentally explored conformational subspace.

Docking to multiple pockets or ligand fields for screening, activity prediction and scaffold hopping

BACKGROUND

Two recent technological advances dramatically reducing the rate of false-negatives in activity prediction by docking flexible 3D models of compounds include multi-conformational docking (mPockDock) and the docking of candidates to atomic property fields derived by co-crystallized ligands (mApfDock).

RESULTS

The mApfDock and mPockDock provide the AUC of 90.4 and 83.8%, respectively. The mApfDock gave better performance when compounds required large induced-fit pocket changes unseen in crystallography, whereas the mPockDock is superior when the co-crystallized ligands do not represent sufficient chemical and binding location diversity.

CONCLUSION

Both approaches proved to be efficient for scaffold hopping; they are complementary when the coverage of the co-crystallized complexes is poor but become convergent when the complexes are diverse enough.

Re-engineering aryl methylcarbamates to confer high selectivity for inhibition of Anopheles gambiae versus human acetylcholinesterase

To identify potential human-safe insecticides against the malaria mosquito we undertook an investigation of the structure-activity relationship of aryl methylcarbamates inhibitors of acetylcholinesterase (AChE). Compounds bearing a β-branched 2-alkoxy or 2-thioalkyl group were found to possess good selectivity for inhibition of Anopheles gambiae AChE over human AChE; up to 530-fold selectivity was achieved with carbamate 11d. A 3D QSAR model is presented that is reasonably consistent with log inhibition selectivity of 34 carbamates. Toxicity of these compounds to live Anopheles gambiae was demonstrated using both tarsal contact (filter paper) and topical application protocols.

Identifying mechanism-of-action targets for drugs and probes

Notwithstanding their key roles in therapy and as biological probes, 7% of approved drugs are purported to have no known primary target, and up to 18% lack a well-defined mechanism of action. Using a chemoinformatics approach, we sought to "de-orphanize" drugs that lack primary targets. Surprisingly, targets could be easily predicted for many: Whereas these targets were not known to us nor to the common databases, most could be confirmed by literature search, leaving only 13 Food and Drug Administration-approved drugs with unknown targets; the number of drugs without molecular targets likely is far fewer than reported. The number of worldwide drugs without reasonable molecular targets similarly dropped, from 352 (25%) to 44 (4%). Nevertheless, there remained at least seven drugs for which reasonable mechanism-of-action targets were unknown but could be predicted, including the antitussives clemastine, cloperastine, and nepinalone; the antiemetic benzquinamide; the muscle relaxant cyclobenzaprine; the analgesic nefopam; and the immunomodulator lobenzarit. For each, predicted targets were confirmed experimentally, with affinities within their physiological concentration ranges. Turning this question on its head, we next asked which drugs were specific enough to act as chemical probes. Over 100 drugs met the standard criteria for probes, and 40 did so by more stringent criteria. A chemical information approach to drug-target association can guide therapeutic development and reveal applications to probe biology, a focus of much current interest.

An alphabetic code based atomic level molecular similarity search in databases

Atomic level molecular similarity and diversity studies have gained considerable importance through their wide application in Bioinformatics and Chemo-informatics for drug design. The availability of large volumes of data on chemical compounds requires new methodologies for efficient and effective searching of its archives in less time with optimal computational power. We describe an alphabetic algorithm for similarity searching based on atom-atom bonding preference for ligands. We represented 170 cyclindependent kinase 2 inhibitors using strings of pre-defined alphabets for searching using known protein sequence alignment tools. Thus, a common pattern was extracted using this set of compounds for database searching to retrieve similar active compounds. Area under the receiver operating characteristic (ROC) curve was used for the discrimination of similar and dissimilar compounds in the databases. An average retrieval rate of about 60% is obtained in cross-validation using the home-grown dataset and the directory of useful decoys (DUD, formally known as the ZINC database) data. This will help in the effective retrieval of similar compounds using database search.

Ligand binding site superposition and comparison based on Atomic Property Fields: identification of distant homologues, convergent evolution and PDB-wide clustering of binding sites

A new binding site comparison algorithm using optimal superposition of the continuous pharmacophoric property distributions is reported. The method demonstrates high sensitivity in discovering both, distantly homologous and convergent binding sites. Good quality of superposition is also observed on multiple examples. Using the new approach, a measure of site similarity is derived and applied to clustering of ligand binding pockets in PDB.