Protein‐Ligand Docking

Combining Stochastic Deformation/Relaxation and Intermolecular Contacts Analysis for Extracting Pharmacophores from Ligand-Receptor Complexes

We previously combined molecular dynamics (classical or simulated annealing) with ligand-receptor contacts analysis as a means to extract valid pharmacophore model(s) from single ligand-receptor complexes. However, molecular dynamics methods are computationally expensive and time-consuming. Here we describe a novel method for extracting valid pharmacophore model(s) from a single crystallographic structure within a reasonable time scale. The new method is based on ligand-receptor contacts analysis following energy relaxation of a predetermined set of randomly deformed complexes generated from the targeted crystallographic structure. Ligand-receptor contacts maintained across many deformed/relaxed structures are assumed to be critical and used to guide pharmacophore development. This methodology was implemented to develop valid pharmacophore models for PI3K-γ, RENIN, and JAK1. The resulting pharmacophore models were validated by receiver operating characteristic (ROC) analysis against inhibitors extracted from the CHEMBL database. Additionally, we implemented pharmacophores extracted from PI3K-γ to search for new inhibitors from the National Cancer Institute list of compounds. The process culminated in new PI3K-γ/mTOR inhibitory leads of low micromolar IC₅₀s.

Simulated annealing molecular dynamics and ligand-receptor contacts analysis for pharmacophore modeling

AIM

Ligand-based pharmacophore modeling requires long list of inhibitors, while pharmacophores based on single ligand-receptor crystallographic structure can be too restricted or promiscuous.

METHODOLOGY

This prompted us to combine simulated annealing molecular dynamics (SAMD) with ligand-receptor contacts analysis as means to construct pharmacophore model(s) from single ligand-receptor complex. Ligand-receptor contacts that survive numerous heating-cooling SAMD cycles are considered critical and are used to guide pharmacophore development.

RESULTS

This methodology was implemented to develop pharmacophores for acetylcholinesterase and protein kinase C-θ. The resulting models were validated by receiver-operating characteristic analysis and in vitro bioassay. Assay identified four new protein kinase C-θ inhibitors among captured hits, two of which exhibited nanomolar potencies.

CONCLUSION

The results illustrate the ability of the new method to extract valid pharmacophores from single ligand-protein complex.

The Comparison of Docking Search Algorithms and Scoring Functions

This chapter, composed of two parts, firstly provides molecular docking overview and secondly two molecular docking case studies are described. In overview, basic information about molecular docking are presented such as description of search algorithms and scoring functions applied in various docking programs. Brief description of methods utilized in some of the most popular docking programs also applied in our experimental work is provided. AutoDock, CDOCKER, GOLD, FlexX and FRED were used for docking studies of the DC-SIGN protein, while GOLD was further used for docking studies of cathepsin K protein. Methods and results of our studies with their contribution to science and medicine are presented. Content of the chapter is therefore appropriate for public of Medicinal and Organic Chemistry as an overview of docking studies, and also for Computational Chemists at the beginning of their work as the introduction to application of molecular docking programs.

The use of docking-based comparative intermolecular contacts analysis to identify optimal docking conditions within glucokinase and to discover of new GK activators

Glucokinase (GK) is involved in normal glucose homeostasis and therefore it is a valid target for drug design and discovery efforts. GK activators (GKAs) have excellent potential as treatments of hyperglycemia and diabetes. The combined recent interest in GKAs, together with docking limitations and shortages of docking validation methods prompted us to use our new 3D-QSAR analysis, namely, docking-based comparative intermolecular contacts analysis (dbCICA), to validate docking configurations performed on a group of GKAs within GK binding site. dbCICA assesses the consistency of docking by assessing the correlation between ligands' affinities and their contacts with binding site spots. Optimal dbCICA models were validated by receiver operating characteristic curve analysis and comparative molecular field analysis. dbCICA models were also converted into valid pharmacophores that were used as search queries to mine 3D structural databases for new GKAs. The search yielded several potent bioactivators that experimentally increased GK bioactivity up to 7.5-folds at 10 μM.

Application of docking-based comparative intermolecular contacts analysis to validate Hsp90α docking studies and subsequent in silico screening for inhibitors

Heat shock protein (Hsp90α) has been recently implicated in cancer, prompting several attempts to discover and optimize new Hsp90α inhibitors. Towards this end, we docked 83 diverse Hsp90α inhibitors into the ATP-binding site of this chaperone using several docking-scoring settings. Subsequently, we applied our newly developed computational tool--docking-based comparative intramolecular contacts analysis (dbCICA)--to assess the different docking conditions and select the best settings. dbCICA is based on the number and quality of contacts between docked ligands and amino acid residues within the binding pocket. It assesses a particular docking configuration based on its ability to align a set of ligands within a corresponding binding pocket in such a way that potent ligands come into contact with binding site spots distinct from those approached by low-affinity ligands, and vice versa. The optimal dbCICA models were translated into valid pharmacophore models that were used as 3D search queries to mine the National Cancer Institute's structural database for new inhibitors of Hsp90α that could potentially be used as anticancer agents. The process culminated in 15 micromolar Hsp90α ATPase inhibitors.

Inhibition of dipeptidyl peptidase IV (DPP IV) is one of the mechanisms explaining the hypoglycemic effect of berberine

Berberine was investigated as an inhibitor of human dipeptidyl peptidase IV (DPP IV) in an attempt to explain its anti-hyperglycemic activities. The investigation included simulated docking experiments to fit berberine within the binding pocket of DPP IV. Berberine was found to readily fit within the binding pocket of DPP IV in a low energy orientation characterized with optimal electrostatic attractive interactions bridging the isoquinolinium positively charged nitrogen atom (berberine) and the negatively charged acidic residue of glutamic acid-205 (GLU205) of DPP IV. Experimentally, berberine was found to inhibit human recombinant DPP IV in vitro with IC(50) = 13.3 microM. Our findings suggest that DPP IV inhibition is, at least, one of the mechanisms that explain the anti-hyperglycemic activity of berberine. The fact that berberine was recently reported to potently inhibit the pro-diabetic target human protein tyrosine phosphatase 1B (h-PTP 1B) discloses a novel dual natural h-PTP 1B/DPP IV inhibitor.

Berberine potently inhibits protein tyrosine phosphatase 1B: Investigation by docking simulation and experimental validation

Berberine was investigated as an inhibitor of human protein tyrosine phosphatase 1B (h-PTP 1B) in an attempt to explain its anti-hyperglycemic activitiy. The investigation included simulated docking experiments to fit berberine within the binding pocket of h-PTP 1B. Berberine was found to readily fit within the binding pocket of h-PTP 1B in a low energy orientation characterized with optimal electrostatic attractive interactions bridging the isoquinolinium positively charged nitrogen atom of berberine and the negatively charged acidic residue of ASP 48 of h-PTP 1B. Experimentally, berberine was found to potently competitively inhibit recombinant h-PTP 1B in vitro (Ki value = 91.3 nM). Our findings strongly suggest that h-PTP 1B inhibition is at least one of the reasons for the reported anti-hyperglycemic activities of berberine.

Enhanced Efficiency of Direct-Space Structure Solution from Powder X-ray Diffraction Data in the Case of Conformationally Flexible Molecules

A strategy is reported for assessing the feasibility of molecular conformations within direct-space structure-solution calculations of organic molecular crystal structures from powder X-ray diffraction data, focusing in particular on the genetic algorithm technique for structure solution in which fitness is defined as a function of the whole-profile figure-of-merit Rwp. The strategy employs a readily computed distance-based function to assess the feasibility of the molecular conformation in each trial structure generated in the genetic algorithm calculation, and structures considered to have low-feasibility conformations are penalized within the evolutionary process. The strategy is shown to lead to significant improvements in the success rate of structure-solution calculations in the case of flexible molecules with a significant number of conformational degrees of freedom.

Combining docking, scoring and molecular field analyses to probe influenza neuraminidase-ligand interactions

In this project, several docking conditions, scoring functions and corresponding protein-aligned molecular field analysis (CoMFA) models were evaluated for a diverse set of neuraminidase (NA) inhibitors. To this end, a group of inhibitors were docked into the active site of NA. The docked structures were utilized to construct a corresponding protein-aligned CoMFA models by employing probe-based (H+, OH, CH3) energy grids and genetic partial least squares (G/PLS) statistical analysis. A total of 16 different docking configurations were evaluated, of which some succeeded in producing self-consistent and predictive CoMFA models. However, the best model coincided with docking the ionized ligands into the hydrated form of the binding site via PLP1 scoring function (r2LOO=0.735, r2PRESS against 24 test compounds=0.828). The highest-ranking CoMFA models were employed to probe NA-ligand interactions. Further validation by comparison with a co-crystallized ligand-NA crystallographic structure was performed. This combination of docking/scoring/CoMFA modeling provided interesting insights into the binding of different NA inhibitors.

Evaluation of library ranking efficacy in virtual screening

We present the results of a comprehensive study in which we explored how the docking procedure affects the performance of a virtual screening approach. We used four docking engines and applied 10 scoring functions to the top-ranked docking solutions of seeded databases against six target proteins. The scores of the experimental poses were placed within the total set to assess whether the scoring function required an accurate pose to provide the appropriate rank for the seeded compounds. This method allows a direct comparison of library ranking efficacy. Our results indicate that the LigandFit/Ligscore1 and LigandFit/GOLD docking/scoring combinations, and to a lesser degree FlexX/FlexX, Glide/Ligscore1, DOCK/PMF (Tripos implementation), LigandFit1/Ligscore2 and LigandFit/PMF (Tripos implementation) were able to retrieve the highest number of actives at a 10% fraction of the database when all targets were looked upon collectively. We also show that the scoring functions rank the observed binding modes higher than the inaccurate poses provided that the experimental poses are available. This finding stresses the discriminatory ability of the scoring algorithms, when better poses are available, and suggests that the number of false positives can be lowered with conformers closer to bioactive ones.