Combining molecular dynamics simulation and ligand-receptor contacts analysis as a new approach for pharmacophore modeling: beta-secretase 1 and check point kinase 1 as case studies

Computational studies for the development of extracellular vesicle-based biosensors

Cancer affects millions of people, and early detection and efficient treatment are two strong levers to hurdle this disease. Recent studies on exosomes, a subset of extracellular vesicles, have deliberately shown the potential to function as a biomarker or treatment tool, thereby attracting the attention of researchers who work on developing biosensors. Due to the ability of computational methods to predict of the behavior of biomolecules, the combination of experimental and computational methods would enhance the analytical performance of the biosensor, including sensitivity, accuracy, and specificity, even detecting such vesicles from bodily fluids. In this regard, the role of computational methods such as molecular docking, molecular dynamics simulation, and density functional theory is overviewed in the development of biosensors. This review highlights the investigations and studies that have been reported using these methods to design exosome-based biosensors. This review concludes with the role of the quantum mechanics/molecular mechanics method in the investigation of chemical processes of biomolecular systems and the deficiencies in using this approach to develop exosome-based biosensors. In addition, the artificial intelligence theory is explained briefly to show its importance in the study of these biosensors.

Computational discovery of novel PI3KC2α inhibitors using structure-based pharmacophore modeling, machine learning and molecular dynamic simulation

PI3KC2α is a lipid kinase associated with cancer metastasis and thrombosis. In this study, we present a novel computational workflow integrating structure-based pharmacophore modeling, machine learning (ML), and molecular dynamics (MD) simulations to discover new PI3KC2α inhibitors. Key innovations include the generation of diverse pharmacophores from both crystallographic and docking-derived complexes, coupled with data augmentation via ligand conformational sampling to enhance ML robustness. The optimal model, developed using XGBoost with genetic function algorithm (GFA) and Shapley additive explanations (SHAP), identified four critical pharmacophores and three descriptors governing bioactivity. Virtual screening of the NCI database using these pharmacophores yielded three hits, with H_1 (NCI: 725847) demonstrating MD-derived binding stability and affinity comparable to the potent inhibitor PITCOIN1 (IC50 = 95 nM). This study represents the first application of a conformation-augmented ML framework to PI3KC2α inhibition, offering a blueprint for targeting underexplored kinases with limited structural data.

Predictive modelling of peroxisome proliferator-activated receptor gamma (PPARγ) IC50 inhibition by emerging pollutants using light gradient boosting machine

Peroxisome proliferator-activated receptor gamma (PPARγ), a critical nuclear receptor, plays a pivotal role in regulating metabolic and inflammatory processes. However, various environmental contaminants can disrupt PPARγ function, leading to adverse health effects. This study introduces a novel approach to predict the inhibitory activity (IC50 values) of 140 chemical compounds across 13 categories, including pesticides, organochlorines, dioxins, detergents, flame retardants, and preservatives, on PPARγ. The predictive model, based on the light-gradient boosting machine (LightGBM) algorithm, was trained on a dataset of 1804 molecules showed r2 values of 0.82 and 0.59, Mean Absolute Error (MAE) of 0.38 and 0.58, and Root Mean Square Error (RMSE) of 0.54 and 0.76 for the training and test sets, respectively. This study provides novel insights into the interactions between emerging contaminants and PPARγ, highlighting the potential hazards and risks these chemicals may pose to public health and the environment. The ability to predict PPARγ inhibition by these hazardous contaminants demonstrates the value of this approach in guiding enhanced environmental toxicology research and risk assessment.

Molecular mechanism underlying effect of D93 and D289 protonation states on inhibitor-BACE1 binding: exploration from multiple independent Gaussian accelerated molecular dynamics and deep learning

BACE1 has been regarded as an essential drug design target for treating Alzheimer's disease (AD). Multiple independent Gaussian accelerated molecular dynamics simulations (GaMD), deep learning (DL), and molecular mechanics general Born surface area (MM-GBSA) method are integrated to elucidate the molecular mechanism underlying the effect of D93 and D289 protonation on binding of inhibitors OV6 and 4B2 to BACE1. The GaMD trajectory-based DL successfully identifies significant function domains. Dynamic analysis shows that the protonation of D93 and D289 strongly affects the structural flexibility and dynamic behaviour of BACE1. Free energy landscapes indicate that inhibitor-bound BACE1s have more conformational states in the protonated states than the wild-type (WT) BACE1, and show more binding poses of inhibitors. Binding affinities calculated using the MM-GBSA method indicate that the protonation of D93 and D289 highly disturbs the binding ability of inhibitors to BACE1. In addition, the protonation of two residues significantly affects the hydrogen bonding interactions (HBIs) of OV6 and 4B2 with BACE1, altering their binding activity to BACE1. The binding hot spots of BACE1 recognized by residue-based free energy estimations provide rational targeting sites for drug design towards BACE1. This study is anticipated to provide theoretical aids for drug development towards treatment of AD.

Advancing 3D printed microfluidics with computational methods for sweat analysis

The intricate tapestry of biomarkers, including proteins, lipids, carbohydrates, vesicles, and nucleic acids within sweat, exhibits a profound correlation with the ones in the bloodstream. The facile extraction of samples from sweat glands has recently positioned sweat sampling at the forefront of non-invasive health monitoring and diagnostics. While extant platforms for sweat analysis exist, the imperative for portability, cost-effectiveness, ease of manufacture, and expeditious turnaround underscores the necessity for parameters that transcend conventional considerations. In this regard, 3D printed microfluidic devices emerge as promising systems, offering a harmonious fusion of attributes such as multifunctional integration, flexibility, biocompatibility, a controlled closed environment, and a minimal requisite analyte volume-features that leverage their prominence in the realm of sweat analysis. However, formidable challenges, including high throughput demands, chemical interactions intrinsic to the printing materials, size constraints, and durability concerns, beset the landscape of 3D printed microfluidic devices. Within this paradigm, we expound upon the foundational aspects of 3D printed microfluidic devices and proffer a distinctive perspective by delving into the computational study of printing materials utilizing density functional theory (DFT) and molecular dynamics (MD) methodologies. This multifaceted approach serves manifold purposes: (i) understanding the complexity of microfluidic systems, (ii) facilitating comprehensive analyses, (iii) saving both cost and time, (iv) improving design optimization, and (v) augmenting resolution. In a nutshell, the allure of 3D printing lies in its capacity for affordable and expeditious production, offering seamless integration of diverse components into microfluidic devices-a testament to their inherent utility in the domain of sweat analysis. The synergistic fusion of computational assessment methodologies with materials science not only optimizes analysis and production processes, but also expedites their widespread accessibility, ensuring continuous biomarker monitoring from sweat for end-users.

Computational workflow for discovering small molecular binders for shallow binding sites by integrating molecular dynamics simulation, pharmacophore modeling, and machine learning: STAT3 as case study

STAT3 belongs to a family of seven transcription factors. It plays an important role in activating the transcription of various genes involved in a variety of cellular processes. High levels of STAT3 are detected in several types of cancer. Hence, STAT3 inhibition is considered a promising therapeutic anti-cancer strategy. However, since STAT3 inhibitors bind to the shallow SH2 domain of the protein, it is expected that hydration water molecules play significant role in ligand-binding complicating the discovery of potent binders. To remedy this issue, we herein propose to extract pharmacophores from molecular dynamics (MD) frames of a potent co-crystallized ligand complexed within STAT3 SH2 domain. Subsequently, we employ genetic function algorithm coupled with machine learning (GFA-ML) to explore the optimal combination of MD-derived pharmacophores that can account for the variations in bioactivity among a list of inhibitors. To enhance the dataset, the training and testing lists were augmented nearly a 100-fold by considering multiple conformers of the ligands. A single significant pharmacophore emerged after 188 ns of MD simulation to represent STAT3-ligand binding. Screening the National Cancer Institute (NCI) database with this model identified one low micromolar inhibitor most likely binds to the SH2 domain of STAT3 and inhibits this pathway.

Identification Mechanism of BACE1 on Inhibitors Probed by Using Multiple Separate Molecular Dynamics Simulations and Comparative Calculations of Binding Free Energies

β-amyloid cleaving enzyme 1 (BACE1) is regarded as an important target of drug design toward the treatment of Alzheimer's disease (AD). In this study, three separate molecular dynamics (MD) simulations and calculations of binding free energies were carried out to comparatively determine the identification mechanism of BACE1 for three inhibitors, 60W, 954 and 60X. The analyses of MD trajectories indicated that the presence of three inhibitors influences the structural stability, flexibility and internal dynamics of BACE1. Binding free energies calculated by using solvated interaction energy (SIE) and molecular mechanics generalized Born surface area (MM-GBSA) methods reveal that the hydrophobic interactions provide decisive forces for inhibitor-BACE1 binding. The calculations of residue-based free energy decomposition suggest that the sidechains of residues L91, D93, S96, V130, Q134, W137, F169 and I179 play key roles in inhibitor-BACE1 binding, which provides a direction for future drug design toward the treatment of AD.

Augmenting bioactivity by docking-generated multiple ligand poses to enhance machine learning and pharmacophore modelling: discovery of new TTK inhibitors as case study

Dual specificity protein kinase threonine/Tyrosine kinase (TTK) is one of the mitotic kinases. High levels of TTK are detected in several types of cancer. Hence, TTK inhibition is considered a promising therapeutic anti-cancer strategy. In this work, we used multiple docked poses of TTK inhibitors to augment training data for machine learning QSAR modeling. Ligand-Receptor Contacts Fingerprints and docking scoring values were used as descriptor variables. Escalating docking-scoring consensus levels were scanned against orthogonal machine learners, and the best learners (Random Forests and XGBoost) were coupled with genetic algorithm and Shapley additive explanations (SHAP) to determine critical descriptors for predicting anti-TTK bioactivity and for pharmacophore generation. Three successful pharmacophores were deduced and subsequently used for in silico screening against the NCI database. A total of 14 hits were evaluated in vitro for their anti-TTK bioactivities. One hit of novel chemotype showed reasonable dose-response curve with experimental IC50 of 1.0 μM. The presented work indicates the validity of data augmentation using multiple docked poses for building successful machine learning models and pharmacophore hypotheses.

Multidimensional virtual screening approaches combined with drug repurposing to identify potential covalent inhibitors of SARS-CoV-2 3CL protease

The outbreak of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has caused an unprecedented global pandemic, and new cases are still on the rise due to the absence of effective medicines. However, developing new drugs within a short time is extremely difficult. Repurposing the existing drugs provides a fast and effective strategy to identify promising inhibitors. Here we focus on the screening of drugs database for discovering potential covalent inhibitors that target 3-chymotrypsin-like protease (3CLpro), an essential enzyme mediating viral replication and transcription. Firstly, we constructed a receptor-ligand pharmacophore model and verified it through decoy set. The importance of pharmacophore features was evaluated by combining molecular dynamics simulation with interaction analyses. Then, covalent docking was used to perform further screening. According to docking score and Prime/Molecular Mechanics Generalized Born Surface Area (MM-GBSA) score, total ten compounds obtained good scores and successfully established covalent bonds with the catalytic Cys145 residue. They also formed favorable interactions with key residues in active sites and closely integrated with 3CLpro with binding modes similar to known 3CLpro inhibitor. Finally, the top four hits DB08732, DB04653, DB01871 and DB07299 were further subjected to 100 ns molecular dynamics (MD) simulation and MM-GBSA binding free energy calculations. The results suggest that the four candidates show good binding affinities for 3CLpro, which warrants further evaluation for their in-vitro/in-vivo activities. Overall, our research methods provide a valuable reference for discovering promising inhibitors against SARS-CoV-2 and help to fight against the epidemic.Communicated by Ramaswamy H. Sarma.

Discovery of new STAT3 inhibitors as anticancer agents using ligand-receptor contact fingerprints and docking-augmented machine learning

STAT3 belongs to a family of seven vital transcription factors. High levels of STAT3 are detected in several types of cancer. Hence, STAT3 inhibition is considered a promising therapeutic anti-cancer strategy. In this work, we used multiple docked poses of STAT3 inhibitors to augment training data for machine learning QSAR modeling. Ligand-Receptor Contact Fingerprints and scoring values were implemented as descriptor variables. Escalating docking-scoring consensus levels were scanned against orthogonal machine learners, and the best learners (Random Forests and XGBoost) were coupled with genetic algorithm and Shapley additive explanations (SHAP) to identify critical descriptors that determine anti-STAT3 bioactivity to be translated into pharmacophore model(s). Two successful pharmacophores were deduced and subsequently used for in silico screening against the National Cancer Institute (NCI) database. A total of 26 hits were evaluated in vitro for their anti-STAT3 bioactivities. Out of which, three hits of novel chemotypes, showed cytotoxic IC50 values in the nanomolar range (35 nM to 6.7 μM). However, two are potent dihydrofolate reductase (DHFR) inhibitors and therefore should have significant indirect STAT3 inhibitory effects. The third hit (cytotoxic IC50 = 0.44 μM) is purely direct STAT3 inhibitor (devoid of DHFR activity) and caused, at its cytotoxic IC50, more than two-fold reduction in the expression of STAT3 downstream genes (c-Myc and Bcl-xL). The presented work indicates that the concept of data augmentation using multiple docked poses is a promising strategy for generating valid machine learning models capable of discriminating active from inactive compounds.

Pharmacophore Modeling of Targets Infested with Activity Cliffs via Molecular Dynamics Simulation Coupled with QSAR and Comparison with other Pharmacophore Generation Methods: KDR as Case Study

Activity cliffs (ACs) are defined as pairs of structurally similar compounds with large difference in their potencies against certain biotarget. We recently proposed that potent AC members induce significant entropically-driven conformational modifications of the target that unveil additional binding interactions, while their weakly-potent counterparts are enthalpically-driven binders with little influence on the protein target. We herein propose to extract pharmacophores for ACs-infested target(s) from molecular dynamics (MD) frames of purely "enthalpic" potent binder(s) complexed within the particular target. Genetic function algorithm/machine learning (GFA/ML) can then be employed to search for the best possible combination of MD pharmacophore(s) capable of explaining bioactivity variations within a list of inhibitors. We compared the performance of this approach with established ligand-based and structure-based methods. Kinase inserts domain receptor (KDR) was used as a case study. KDR plays a crucial role in angiogenic signalling and its inhibitors have been approved in cancer treatment. Interestingly, GFA/ML selected, MD-based, pharmacophores were of comparable performances to ligand-based and structure-based pharmacophores. The resulting pharmacophores and QSAR models were used to capture hits from the national cancer institute list of compounds. The most active hit showed anti-KDR IC50 of 2.76 μM.

Discovery of new Cdc2-like kinase 4 (CLK4) inhibitors via pharmacophore exploration combined with flexible docking-based ligand/receptor contact fingerprints and machine learning

Cdc2-like kinase 4 (CLK4) inhibitors are of potential therapeutic value in many diseases particularly cancer. In this study, we combined extensive ligand-based pharmacophore exploration, ligand-receptor contact fingerprints generated by flexible docking, physicochemical descriptors and machine learning-quantitative structure-activity relationship (ML-QSAR) analysis to investigate the pharmacophoric/binding requirements for potent CLK4 antagonists. Several ML methods were attempted to tie these properties with anti-CLK4 bioactivities including multiple linear regression (MLR), random forests (RF), extreme gradient boosting (XGBoost), probabilistic neural network (PNN), and support vector regression (SVR). A genetic function algorithm (GFA) was combined with each method for feature selection. Eventually, GFA-SVR was found to produce the best self-consistent and predictive model. The model selected three pharmacophores, three ligand-receptor contacts and two physicochemical descriptors. The GFA-SVR model and associated pharmacophore models were used to screen the National Cancer Institute (NCI) structural database for novel CLK4 antagonists. Three potent hits were identified with the best one showing an anti-CLK4 IC50 value of 57 nM.

Exploiting activity cliffs for building pharmacophore models and comparison with other pharmacophore generation methods: sphingosine kinase 1 as case study

Activity cliffs (ACs) are defined as closely analogous compounds of significant affinity discrepancies against certain biotarget. In this paper we propose to use AC pair(s) for extracting valid binding pharmacophores through exposing corresponding protein complexes to stochastic deformation/relaxation followed by applying genetic algorithm/machine learning (GA-ML) for selecting optimal pharmacophore(s) that best classify a long list of inhibitors. We compared the performances of ligand-based and structure-based pharmacophores with counterparts generated by this newly introduced technique. Sphingosine kinase 1 (SPHK-1) was used as case study. SPHK-1 is a lipid kinase that plays pivotal role in the regulation of a variety of biological processes including, cell growth, apoptosis, and inflammation. The new approach proved to yield pharmacophore and ML models of comparable accuracies to established ligand-based and structure-based pharmacophores. The resulting pharmacophores and ML models were used to capture hits from the national cancer institute list of compounds and predict their bioactivity categories. Two hits of novel chemotypes showed selective and low micromolar inhibitory IC50 values against SPHK-1.

1,4-Naphthoquinone Is a Potent Inhibitor of IRAK1 Kinases and the Production of Inflammatory Cytokines in THP-1 Differentiated Macrophages

Quinones are a class of cyclic organic compounds that are widely distributed in nature and have been shown to exhibit anti-inflammatory, antioxidant, and anticancerous activities. However, the molecular mechanisms/signaling by which these molecules exert their effect are still not fully understood. In this study, a group of quinone-derived compounds were examined for their potential inhibitory effect against human IRAK1 and IRAK4 kinases in vitro. We have identified five compounds: 1,4-naphthoquinone, emodin, shikonin, plumbagin, and menadione (vitamin K3) as active and selective inhibitors of human IRAK1 enzyme in vitro. The biochemical binding and molecular interactions between the active compounds and IRAK1's catalytic site were demonstrated in silico using structural-based docking and dynamic simulation analysis. Also, 1,4-naphthoquinone was found to effectively inhibit the growth of cancer cell lines overexpressing IRAK1. Furthermore, 1,4-naphthoquinone potently suppressed the production and secretion of key proinflammatory cytokine proteins IL-8, IL-1β, IL-10, TNF-α, and IL-6 in LPS-stimulated PMA-induced human THP-1 macrophages. In conclusion, 1,4-naphthoquinone is an effective inhibitor of IRAK1 kinases and their mediated inflammatory cytokines production in LPS-stimulated PMA-induced human THP-1 macrophages.

Docking-generated multiple ligand poses for bootstrapping bioactivity classifying Machine Learning: Repurposing covalent inhibitors for COVID-19-related TMPRSS2 as case study

In the present work we introduce the use of multiple docked poses for bootstrapping machine learning-based QSAR modelling. Ligand-receptor contact fingerprints are implemented as descriptor variables. We implemented this method for the discovery of potential inhibitors of the serine protease enzyme TMPRSS2 involved the infectivity of coronaviruses. Several machine learners were scanned, however, Xgboost, support vector machines (SVM) and random forests (RF) were the best with testing set accuracies reaching 90%. Three potential hits were identified upon using the method to scan known untested FDA approved drugs against TMPRSS2. Subsequent molecular dynamics simulation and covalent docking supported the results of the new computational approach.

Dynamic-shared Pharmacophore Approach as Tool to Design New Allosteric PRC2 Inhibitors, Targeting EED Binding Pocket

The Polycomb Repressive complex 2 (PRC2) maintains a repressive chromatin state and silences many genes, acting as methylase on histone tails. This enzyme was found overexpressed in many types of cancer. In this work, we have set up a Computer-Aided Drug Design approach based on the allosteric modulation of PRC2. In order to minimize the possible bias derived from using a single set of coordinates within the protein-ligand complex, a dynamic workflow was developed. In details, molecular dynamic was used as tool to identify the most significant ligand-protein interactions from several crystallized protein structures. The identified features were used for the creation of dynamic pharmacophore models and docking grid constraints for the design of new PRC2 allosteric modulators. Our protocol was retrospectively validated using a dataset of active and inactive compounds, and the results were compared to the classic approaches, through ROC curves and enrichment factor. Our approach suggested some important interaction features to be adopted for virtual screening performance improvement.

Discovery of new JNK3 inhibitory chemotypes via QSAR-Guided selection of docking-based pharmacophores and comparison with other structure-based pharmacophore modeling methods

The kinase c-Jun N-terminal Kinase 3 (JNK3) plays an important role in neurodegenerative diseases. JNK3 inhibitors have shown promising results in treating Alzheimer's and Parkinson's diseases. This prompted us to model this interesting target via three established structure-based computational workflows; namely, docking-based Comparative Intermolecular Contacts Analysis (db-CICA), pharmacophore modeling via molecular-dynamics based Ligand-Receptor Contact Analysis (md-LRCA), and QSAR-guided selection of crystallographic pharmacophores. Moreover, we compared the performances of resulting pharmacophores against binding models generated via a newly introduced technique, namely, QSAR-guided selection of docking-based pharmacophores. The resulting pharmacophores were validated by receiver operating characteristic (ROC) curve analysis and used as virtual search queries to screen the National Cancer Institute (NCI) database for promising anti-JNK3 hits of novel chemotypes. Eleven nanomolar and low micromolar hits were identified, three of which were captured by QSAR-guided docking-based pharmacophores.

Pharmacophore modeling coupled with molecular dynamic simulation approach to identify new leads for meprin-β metalloprotease

Human meprin beta metalloprotease, a small subgroup of the astacin family, is a potent drug target for the treatment of several disorders such as fibrosis, neurodegenerative disease in particular Alzheimer and inflammatory bowel diseases. In this study, a ligand-based pharmacophore approach has been used for the selection of potentially active compounds to understand the inhibitory activities of meprin-β by using the sulfonamide scaffold based inhibitors. Using this dataset, a pharmacophore model (Hypo1) was selected on the basis of a highest correlation coefficient (0.959), lowest total cost (105.89) and lowest root mean square deviation (1.31 Å) values. All the pharmacophore hypotheses generated from the candidate inhibitors comprised four features: two hydrogen-bond acceptor, one hydrogen-bond donor and one zinc binder feature. The best validated pharmacophore model (Hypo1) was used for virtual screening of compounds from several databases. The selective hit compounds were filtered by drug likeness property, acceptable ADMET profile, molecular docking and DFT study. Molecular dynamic simulations with the final 10 hit compounds revealed that a large number of non-covalent interactions were formed with the active site and specificity sub-pockets of the meprin beta metalloprotease. This study assists in the development of the new lead molecules as well as gives a better understanding of their interaction with meprin-β.

Discovery of New Phosphoinositide 3-kinase Delta (PI3Kδ) Inhibitors via Virtual Screening using Crystallography-derived Pharmacophore Modelling and QSAR Analysis

BACKGROUND

PI3Kδ is predominantly expressed in hematopoietic cells and participates in the activation of leukocytes. PI3Kδ inhibition is a promising approach for treating inflammatory diseases and leukocyte malignancies. Accordingly, we decided to model PI3Kδ binding.

METHODS

Seventeen PI3Kδ crystallographic complexes were used to extract 94 pharmacophore models. QSAR modelling was subsequently used to select the superior pharmacophore(s) that best explain bioactivity variation within a list of 79 diverse inhibitors (i.e., upon combination with other physicochemical descriptors).

RESULTS

The best QSAR model (r2 = 0.71, r2 LOO = 0.70, r2 press against external testing list of 15 compounds = 0.80) included a single crystallographic pharmacophore of optimal explanatory qualities. The resulting pharmacophore and QSAR model were used to screen the National Cancer Institute (NCI) database for new PI3Kδ inhibitors. Two hits showed low micromolar IC50 values.

CONCLUSION

Crystallography-based pharmacophores were successfully combined with QSAR analysis for the identification of novel PI3Kδ inhibitors.

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 IC50s.

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.