Molecular dynamics simulation study of PTP1B with allosteric inhibitor and its application in receptor based pharmacophore modeling

Allosteric inhibition of protein tyrosine phosphatase 1B (PTP1B), has paved a new path to design specific inhibitors for PTP1B, which is an important drug target for the treatment of type II diabetes and obesity. The PTP1B_1–282-allosteric inhibitor complex crystal structure lacks α7 (287–298) and moreover there is no available 3D structure of PTP1B_1–298 in open form. As the interaction between α7 and α6–α3 helices plays a crucial role in allosteric inhibition, α7 was modeled to the PTP1B_1–282 in open form complexed with an allosteric inhibitor (compound-2) and a 5 ns MD simulation was performed to investigate the relative orientation of the α7–α6–α3 helices. The simulation conformational space was statistically sampled by clustering analyses. This approach was helpful to reveal certain clues on PTP1B allosteric inhibition. The simulation was also utilized in the generation of receptor based pharmacophore models to include the conformational flexibility of the protein-inhibitor complex. Three cluster representative structures of the highly populated clusters were selected for pharmacophore model generation. The three pharmacophore models were subsequently utilized for screening databases to retrieve molecules containing the features that complement the allosteric site. The retrieved hits were filtered based on certain drug-like properties and molecular docking simulations were performed in two different conformations of protein. Thus, performing MD simulation with α7 to investigate the changes at the allosteric site, then developing receptor based pharmacophore models and finally docking the retrieved hits into two distinct conformations will be a reliable methodology in identifying PTP1B allosteric inhibitors. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Pharmacophore generation,atom-based 3D-QSAR and molecular dynamics simulation analyses of pyridine-3-carboxamide-6-yl-urea analogues as potential gyrase B inhibitors

DNA gyrase subunit B (GyrB) is an attractive drug target for the development of antibacterial agents with therapeutic potential. In the present study, computational studies based on pharmacophore modelling, atom-based QSAR, molecular docking, free binding energy calculation and dynamics simulation were performed on a series of pyridine-3-carboxamide-6-yl-urea derivatives. A pharmacophore model using 49 molecules revealed structural and chemical features necessary for these molecules to inhibit GyrB. The best fitted model AADDR.13 was generated with a coefficient of determination (r²) of 0.918. This model was validated using test set molecules and had a good r² of 0.78. 3D contour maps generated by the 3D atom-based QSAR revealed the key structural features responsible for the GyrB inhibitory activity. Extra precision molecular docking showed hydrogen bond interactions with key amino acid residues of ATP-binding pocket, important for inhibitor binding. Further, binding free energy was calculated by the MM-GBSA rescoring approach to validate the binding affinity. A 10 ns MD simulation of inhibitor #47 showed the stability of the predicted binding conformations. We identified 10 virtual hits by in silico high-throughput screening. A few new molecules were also designed as potent GyrB inhibitors. The information obtained from these methodologies may be helpful to design novel inhibitors of GyrB.     

In Silico Strategy for Targeting the mTOR Kinase at Rapamycin Binding Site by Small Molecules

Computer aided drug-design methods proved to be powerful tools for the identification of new therapeutic agents. We employed a structure-based workflow to identify new inhibitors targeting mTOR kinase at rapamycin binding site. By combining molecular dynamics (MD) simulation and pharmacophore modelling, a simplified structure-based pharmacophore hypothesis was built starting from the FKBP12-rapamycin-FRB ternary complex retrieved from RCSB Protein Data Bank (PDB code 1FAP). Then, the obtained model was used as filter to screen the ZINC biogenic compounds library, containing molecules derived from natural sources or natural-inspired compounds. The resulting hits were clustered according to their similarity; moreover, compounds showing the highest pharmacophore fit-score were chosen from each cluster. The selected molecules were subjected to docking studies to clarify their putative binding mode. The binding free energy of the obtained complexes was calculated by MM/GBSA method and the hits characterized by the lowest ΔG_bind values were identified as potential mTOR inhibitors. Furthermore, the stability of the resulting complexes was studied by means of MD simulation which revealed that the selected compounds were able to form a stable ternary complex with FKBP12 and FRB domain, thus underlining their potential ability to inhibit mTOR with a rapamycin-like mechanism.     

Targeting tribbles homolog 3 (TRIB3) protein against type 2 diabetes for the identification of potential inhibitors by in silico screening

Tribbles homolog 3 (TRIB3) protein is inhibiting the insulin signaling by directly binding to the Akt/PKB leading to insulin resistance in the pancreas causing type 2 diabetes mellitus. Hence, TRIB3 protein is considered as a possible drug target for the new lead identification against type 2 diabetes. In the present study, the homology model of TRIB3 protein was generated to explore its biochemical function and molecular interactions in the new lead identification. The energy minimization of TRIB3 protein was carried out and evaluated by validation protocols for structure reliability. The druggable binding site of TRIB3 protein was identified for the virtual screening and molecular docking studies. The Asinex-fragments library of 22634 small molecules was docked at TRIB3 active site using the Glide module to identify new chemical entities. A total of 9 molecules were identified as final hits from virtual screening and their potency was ranked using Glide score, Glide energies, and residues interactions. The 6 prioritized lead molecules were further optimized using AutoDock, Prime MM/GBSA, and percentage of human oral absorption for the identification of potential leads. The molecules L2, L5, and L6 are identified as lead inhibitors and are showing consistent interactions with key residues Glu194 and Lys196 of TRIB3 protein. The identified potential leads were analyzed by ADME properties for their drug likeness and HergIC50 values are predicted for the prevention of preclinical failures. The present work sheds light on the identification of the best lead molecules against TRIB3 protein and offers a route to design as novel potential drug candidates for T2DM.     

Identification of potential HIV-1 integrase strand transfer inhibitors: In silico virtual screening and QM/MM docking studies

HIV-1 integrase (IN) is a retroviral enzyme that catalyses integration of the reverse-transcribed viral DNA into the host genome, which is necessary for efficient viral replication. In this study, we have performed an in silico virtual screening for the identification of potential HIV-1 IN strand transfer (ST) inhibitors. Pharmacophore modelling and atom-based 3D-QSAR studies were carried out for a series of compounds belonging to 3-Hydroxypyrimidine-2,4-diones. Based on the ligand-based pharmacophore model, we obtained a five-point pharmacophore with two hydrogen bond acceptors (A), one hydrogen bond donor (D), one hydrophobic group (H) and one aromatic ring (R) as pharmacophoric features. The pharmacophore hypothesis AADHR was used as a 3D query in a sequential virtual screening study to filter small molecule databases Maybridge, ChemBridge and Asinex. Hits matching with pharmacophore hypothesis AADHR were retrieved and passed progressively through Lipinski’s rule of five filtering, molecular docking and hierarchical clustering. The five compounds with best hits with novel and diverse chemotypes were subjected to QM/MM docking, which showed improved docking accuracy. We further performed molecular dynamics simulation and found three compounds that form stable interactions with key residues. These compounds could be used as a leads for further drug development and rational design of HIV-1 IN inhibitors.     

New leads for selective GSK-3 inhibition: pharmacophore mapping and virtual screening studies

Summary Glycogen Synthase Kinase-3 is a regulatory serine/threonine kinase, which is being targeted for the treatment of a number of human diseases including type-2 diabetes mellitus, neurodegenerative diseases, cancer and chronic inflammation. Selective GSK-3 inhibition is an important requirement owing to the possibility of side effects arising from other kinases. A pharmacophore mapping strategy is employed in this work to identify new leads for selective GSK-3 inhibition. Ligands known to show selective GSK-3 inhibition were employed in generating a pharmacophore map using distance comparison method (DISCO). The derived pharmacophore map was validated using (i) important interactions involved in selective GSK-3 inhibitions, and (ii) an in-house database containing different classes of GSK-3 selective, non-selective and inactive molecules. New Lead identification was carried out by performing virtual screening using validated pharmacophoric query and three chemical databases namely NCI, Maybridge and Leadquest. Further data reduction was carried out by employing virtual filters based on (i) Lipinski’s rule of 5 (ii) van der Waals bumps and (iii) restricting the number of rotatable bonds to seven. Final screening was carried out using FlexX based molecular docking study.     

Identification of novel nt-MGAM inhibitors for potential treatment of type 2 diabetes: Virtual screening,atom based 3D-QSAR model,docking analysis and ADME study

In this study, a virtual screening procedure was applied to identify new potential nt-MGAM inhibitors as a possible medication for type 2 diabetes. To this aim, a series of salacinol analogues were first investigated by docking analysis for their binding to the X-ray structure of the biological target nt-MGAM. Key interactions for ligand binding into the receptor active site were identified which shared common features to those found for other known inhibitors, which strengthen the results of this study. 3D QSAR model was then built and showed to be statistically significant and with a good predictive power for the training (R² = 0.99, SD = 0.17, F = 555.3 and N = 27) and test set (Q² = 0.81, Pearson(r) = 0.92, RMSE = 0.52, N = 08). The model was then used to virtually screen the ZINC database with the aim of identifying novel chemical scaffolds as potential nt-MGAM inhibitors. Further, in silico predicted ADME properties were investigated for the most promising molecules. The outcome of this investigation sheds light on the molecular characteristics of the binding of salacinol analogues to nt-MGAM enzyme and identifies new possible inhibitors which have the potential to be developed into drugs, thus significantly contributing to the design and optimization of therapeutic strategies against type 2 diabetes.     

Pharmacophore modelling and electronic feature analysis of hydroxamic acid derivatives,the HIV integrase inhibitors

Hydroxamic acid derivatives with metal ion binding properties were collected from the literature to generate a pharmacophore and 3D-QSAR model for HIV strand transfer inhibition. The derived pharmacophore model (AAAHRR) recognizes both metal ion binding site and hydrophobic group. The QSAR model generated using this hypothesis expressed statistical significance (r ² = 0.971 for the training set and q ² = 0.913 for the test set). The ability of this pharmacophore model to retrieve other metal ion binding inhibitors was examined by screening the ChemBank database (ligandinfo) incorporated with 10 known strand transfer inhibitors. The studied favourable and unfavourable contours of chemical features (H-bond donor, acceptor and hydrophobic sites) revealed the role of hydrophobic substitution at the fluorobenzene ring and cyclization of the metal ion binding hydroxamic acid in effective integrase inhibition. Analysis of the frontier orbitals, HOMO and LUMO revealed that the nucleophilic / electrophilic interactions depend on the significant overlapping observed at the azaindole and hydroxamic acid groups. In essence, the generated pharmacophore model is competent enough to disclose the essential site-specific interactions involved in the inhibition of HIV integrase, and hence can be used in virtual screening to identify novel scaffolds as leads with increased anti-viral potency.     

Homology modeling,docking and structure-based pharmacophore of inhibitors of DNA methyltransferase

DNA methyltransferase 1 (DNMT1) is an emerging epigenetic target for the treatment of cancer and other diseases. To date, several inhibitors from different structural classes have been published. In this work, we report a comprehensive molecular modeling study of 14 established DNTM1 inhibitors with a herein developed homology model of the catalytic domain of human DNTM1. The geometry of the homology model was in agreement with the proposed mechanism of DNA methylation. Docking results revealed that all inhibitors studied in this work have hydrogen bond interactions with a glutamic acid and arginine residues that play a central role in the mechanism of cytosine DNA methylation. The binding models of compounds such as curcumin and parthenolide suggest that these natural products are covalent blockers of the catalytic site. A pharmacophore model was also developed for all DNMT1 inhibitors considered in this work using the most favorable binding conformations and energetic terms of the docked poses. To the best of our knowledge, this is the first pharmacophore model proposed for compounds with inhibitory activity of DNMT1. The results presented in this work represent a conceptual advance for understanding the protein–ligand interactions and mechanism of action of DNMT1 inhibitors. The insights obtained in this work can be used for the structure-based design and virtual screening for novel inhibitors targeting DNMT1.     

Development of predictive pharmacophore model for in silico screening,and 3D QSAR CoMFA and CoMSIA studies for lead optimization,for designing of potent tumor necrosis factor alpha converting enzyme inhibitors

A chemical feature-based pharmacophore model was developed for Tumor Necrosis Factor-α converting enzyme (TACE) inhibitors. A five point pharmacophore model having two hydrogen bond acceptors (A), one hydrogen bond donor (D) and two aromatic rings (R) with discrete geometries as pharmacophoric features was developed. The pharmacophore model so generated was then utilized for in silico screening of a database. The pharmacophore model so developed was validated by using four compounds having proven TACE inhibitory activity which were grafted into the database. These compounds mapped well onto the five listed pharmacophoric features. This validated pharmacophore model was also used for alignment of molecules in CoMFA and CoMSIA analysis. The contour maps of the CoMFA/CoMSIA models were utilized to provide structural insight for activity improvement of potential novel TACE inhibitors. The pharmacophore model so developed could be used for in silico screening of any commercial/in house database for identification of TACE inhibiting lead compounds, and the leads so identified could be optimized using the developed CoMSIA model. The present work highlights the tremendous potential of the two mutually complementary ligand-based drug designing techniques (i.e. pharmacophore mapping and 3D-QSAR analysis) using TACE inhibitors as prototype biologically active molecules.     

Combining molecular docking and molecular dynamics studies for modelling Staphylococcus aureus MurD inhibitory activity

The ATP-dependent bacterial MurD enzyme catalyses the formation of the peptide bond between cytoplasmic intermediate UDP-N-acetylmuramoyl-L-alanine and D-glutamic acid. This is essential for bacterial cell wall peptidoglycan synthesis in both Gram-positive and Gram-negative bacteria. MurD is recognized as an important target for the development of new antibacterial agents. In the present study we prepared the 3D-stucture of the catalytic pocket of the Staphylococcus aureus MurD enzyme by homology modelling. Extra-precision docking, binding free energy calculation by the MM–GBSA approach and a 40 ns molecular dynamics (MD) simulation of 2-thioxothiazolidin-4-one based inhibitor $1 was carried out to elucidate its inhibition potential for the S. aureus MurD enzyme. Molecular docking results showed that Lys19, Gly147, Tyr148, Lys328, Thr330 and Phe431 residues are responsible for the inhibitor–protein complex stabilization. Binding free energy calculation revealed electrostatic solvation and van der Waals energy components as major contributors for the inhibitor binding. The inhibitor-modelled S. aureus protein complex had a stable conformation in response to the atomic flexibility and interaction, when subjected to MD simulation at 40 ns in aqueous solution. We designed some molecules as potent inhibitors of S. aureus MurD, and to validate the stability of the designed molecule D1-modelled protein complex we performed a 20 ns MD simulation. Results obtained from this study can be utilized for the design of potent S. aureus MurD inhibitors.     

Comparison of a homology model and the crystallographic structure of human 11β-hydroxysteroid dehydrogenase type 1 (11βHSD1) in a structure-based identification of inhibitors

Summary Human 11β-hydroxysteroid dehydrogenase type 1 (11βHSD1) catalyzes the interconversion of cortisone into active cortisol. 11βHSD1 inhibition is a tempting target for the treatment of a host of human disorders that might benefit from blockade of glucocorticoid action, such as obesity, metabolic syndrome, and diabetes type 2. Here, we report an in silico screening study aimed at identifying new selective inhibitors of human 11βHSD1 enzyme. In the first step, homology modeling was employed to build the 3D structure of 11βHSD1. Further, molecular docking was used to validate the predicted model by showing that it was able to discriminate between known 11βHSD1 inhibitors or substrates and non-inhibitors. The homology model was found to reproduce closely the crystal structure that became publicly available in the final stages of this work. Finally, we carried out structure-based virtual screening experiments on both the homology model and the crystallographic structure with a database of 114’000 natural molecules. Among these, 15 molecules were consistently selected as inhibitors based on both the model and crystal structures of the enzyme, implying a good quality for the homology model. Among these putative 11βHSD1 inhibitors, two were flavonone derivatives that have already been shown to be potent inhibitors of the enzyme.     

Pharmacophore Screening on Piperidinecarboxamides Derivatives Based on GALAHAD and CoMFA Models

In order to discover the novel anticonvulsant drugs, pharmacophore screening of the anticonvulsant inhibitors was enforced. Genetic Algorithm with Linear Assignment for Hypermolecular Alignment of Datasets (GALAHAD) and Comparative Molecular Field Analysis (CoMFA) studies were combined to implement our research. Firstly, multiple models were generated using GALAHAG based on high active molecules. Secondly, several of them were validated using the CoMFA study. Finally, a good values of q² from training set and promising predictive power from test set were obtained based on one model simutaneously. One model had been selected as the most reasonable pharmacophore model. The results of the CoMFA study based on the model 1 suggested that both steric and electrostatic interactions played important roles.     

Design of new CD38 inhibitors based on CoMFA modelling and molecular docking analysis of 4‑amino-8-quinoline carboxamides and 2,4-diamino-8-quinazoline carboxamides

In this study, based on molecular docking analysis and comparative molecular field analysis (CoMFA) modelling of a series of 71 CD38 inhibitors including 4?amino-8-quinoline carboxamides and 2,4-diamino-8-quinazoline carboxamides, new CD38 inhibitors were designed. The interactions of the molecules with the greatest and the lowest activities with the nicotinamide mononucleotide (NMN) binding site were investigated by molecular docking analysis. A CoMFA model with four partial least squares regression (PLSR) components was developed to predict the CD38 inhibitory activity of the molecules. The r² values for the training and test sets were 0.89 and 0.82, respectively. The Q² values for leave-one-out cross-validation (LOO-CV) and leave-many-out cross-validation (LMO-CV) tests on the training set were 0.65 and 0.64, respectively. The CoMFA model was validated by calculating several statistical parameters. CoMFA contour maps were interpreted, and structural features that influence the CD38 inhibitory activity of molecules were determined. Finally, seven new CD38 inhibitors with greater activity with respect to the greatest active molecules were designed.     

Binding mode of inhibitors and Cryptosporidium parvum IMP dehydrogenase: A combined ligand- and receptor-based study

A combined ligand- and target-based approach was used to analyse the interaction models of Cryptosporidium parvum inosine 5’-monophosphate dehydrogenase (CpIMPDH) with selective inhibitors. First, a ligand-based pharmacophore model was generated from 20 NAD⁺ competitive CpIMPDH inhibitors with the HipHop module. The characteristic of the NAD⁺ binding site of CpIMPDH was then described, and the binding modes of the representative inhibitors were studied by molecular docking. The combination of the pharmacophore model and the docking results allowed us to evaluate the pharmacophore features and structural information of the NAD⁺ binding site of CpIMPDH. This research supports the proposal of an interaction model inside the NAD⁺ binding site of CpIMPDH, consisting of four key interaction points: two hydrophobic-aromatic groups, a hydrophobic-aliphatic group and a hydrogen bond donor. This study also provides guidance for the design of more potent CpIMPDH inhibitors for the treatment of Cryptosporidium infections.     

Identification of new dual spleen tyrosine kinase (Syk) and phosphoionositide-3-kinase δ (PI3Kδ) inhibitors using ligand and structure-based integrated ideal pharmacophore models

Owing to the complex pathophysiology of autoimmune disorders, it is very challenging to develop successful treatment strategies. Single-target agents are not desired therapeutics for such multi-factorial disorders. Considering the current need for the treatment of complex autoimmune disorders, dual inhibitors of Syk and PI3Kδ have been designed using ligand and structure-based molecular modelling strategies. In the present work, structure and ligand-based pharmacophore modelling was implemented for a varied set of Syk and PI3Kδ inhibitors. Ligand-based pharmacophore models (LBPMs) were developed for two kinases: ADPR.14 (r²_train = 0.809) for Syk, comprising one hydrogen bond acceptor, one hydrogen bond donor, one positive ionisable and one ring aromatic feature, and for PI3Kδ: AAARR.45 (r²_train = 0.942) consisting of three hydrogen bond acceptor and two ring aromatic features. The generated e-pharmacophore models revealed an additional ring aromatic and hydrophobic feature important for Syk and PI3Kδ inhibition, respectively. Subsequently, LBPMs were modified resulting in APDRR.14 hypothesis for Syk inhibitors and AAAHRR.45 hypothesis for PI3Kδ inhibitors employed for virtual screening. Thus, the combination of ligand and structure-based pharmacophore modelling helped in developing ideal pharmacophore models that may be an efficient tool for the designing of novel dual inhibitors of Syk and PI3Kδ.     

Pharmacophore Mode lingand Virtual Screening to Design the Potential Influenza Virus Endonuclease Inhibitors

Influenza virus endonuclease is an attractive target for antiviral therapy in the treatment of influenza infection. The purpos e of this study is to design a novel antiviral agent with improved biological activities against the influenza virus endonuclease. In this study, chemical feature‐based 3D pharmacophore models were developed from 41 known influenza virus endonuclease inhibitors. The best quantitative pharmacohore model (Hypo 1), which consists of two hydrogen‐bond acceptors and two hydrophobic features, yields the highest correlation coefficient (R = 0.886). Hypo 1 was further validated by the cross validation method and the test set compounds. The application of this model for predicting the activities of 11 known influenza virus endonuclease inhibitors in the test set shows great success. The correlation coefficient of 0.942 and a cross validation of 95;% confidence level prove that this model is reliable in identifying structurally diverse compounds for influenza virus endonuclease inhibition. The most active compound (compound 1) from the training set was docked into the active site of the influenza virus endonuclease as an additional verification that the pharmacophore model is accurate. The docked conformation showed important hydrogen bond interactions between the compound and two amino acids, Lys 134 and Lys 137. After validation, this model was used to screen the NCI chemical database to identify new influenza virus endonuclease inhibitors. Our study shows that the to pranking compound out of the 10 newly identified compounds using fit value ranking has an estimated activity of 0.049 μM. These newly identified lead compounds can be further experimentally validated using in vitro techniques.     

Hierarchical virtual screening: identification of potential high-affinity and selective β(3)-adrenergic receptor agonists

The hierarchical virtual screening (HVS) study, consisting of pharmacophore modelling, docking and VS of the generated focussed virtual library, has been carried out to identify novel high-affinity and selective β(3)-adrenergic receptor (β-AR) agonists. The best pharmacophore model, comprising one H-bond donor, two hydrophobes, one positive ionizable and one negative ionizable feature, was developed based on a training set of 51 β(3)-AR agonists using the pharmacophore generation protocol implemented in Discovery Studio. The model was further validated with the test set, external set and ability of the pharmacophoric features to complement the active site amino acids of the homology modelled β(3)-AR developed using MODELLER software. The focussed virtual library was generated using the structure-based insights gained from our earlier reported comprehensive study focussing on the structural basis of β-AR subtype selectivity of representative agonists and antagonists. The HVS with the sequential use of the best pharmacophore model and homology modelled β(3)-AR in the screening of the generated focussed library has led to the identification of potential virtual leads as novel high-affinity and selective β(3)-AR agonists.     

A structure-based design approach for the identification of novel inhibitors: application to an alanine racemase

We report a new structure-based strategy for the identification of novel inhibitors. This approach has been applied to Bacillus stearothermophilus alanine racemase (AlaR), an enzyme implicated in the biosynthesis of the bacterial cell wall. The enzyme catalyzes the racemization of l- and d-alanine using pyridoxal 5-phosphate (PLP) as a cofactor. The restriction of AlaR to bacteria and some fungi and the absolute requirement for d-alanine in peptidoglycan biosynthesis make alanine racemase a suitable target for drug design. Unfortunately, known inhibitors of alanine racemase are not specific and inhibit the activity of other PLP-dependent enzymes, leading to neurological and other side effects.This article describes the development of a receptor-based pharmacophore model for AlaR, taking into account receptor flexibility (i.e. a `dynamic' pharmacophore model). In order to accomplish this, molecular dynamics (MD) simulations were performed on the full AlaR dimer from Bacillus stearothermophilus (PDB entry, 1sft) with a d-alanine molecule in one active site and the non-covalent inhibitor, propionate, in the second active site of this homodimer. The basic strategy followed in this study was to utilize conformations of the protein obtained during MD simulations to generate a dynamic pharmacophore model using the property mapping capability of the LigBuilder program. Compounds from the Available Chemicals Directory that fit the pharmacophore model were identified and have been submitted for experimental testing.The approach described here can be used as a valuable tool for the design of novel inhibitors of other biomolecular targets.     

Pharmacophore identification, in silico screening, and virtual library design for inhibitors of the human factor Xa

Factor Xa inhibitors are innovative anticoagulant agents that provide a better safety/efficacy profile compared to other anticoagulative drugs. A chemical feature-based modeling approach was applied to identify crucial pharmacophore patterns from 3D crystal structures of inhibitors bound to human factor Xa (Pdb entries 1fjs, 1kns, 1eqz) using the software LIGANDSCOUT and CATALYST. The complex structures were selected regarding the criteria of high inhibitory potency (i.e. all ligands show K(i) values against factor Xa in the subnanomolar range) and good resolution (i.e. at least 2.2 A) in order to generate selective and high quality pharmacophore models. The resulting chemical-feature based hypotheses were used for virtual screening of commercial molecular databases such as the WDI database. Furthermore, a ligand-based molecular modeling approach was performed to obtain common-feature hypotheses that represent the relevant chemical interactions between 10 bioactive factor Xa inhibitors and the protein, respectively. In a next step a virtual combinatorial library was designed in order to generate new compounds with similar chemical and spatial properties as known inhibitors. The software tool ILIB DIVERSE was used for this procedure in order to provide new scaffolds of this group of anticoagulants. Finally we present the combination of these two techniques, hence virtual screening was performed with selective pharmacophore models in a focused virtual combinatorial database. De novo derived molecular scaffolds that were able to adequately satisfy the pharmacophore criteria are revealed and are promising templates for candidates for further development.