Molecular modeling study on chemically diverse series of cyclooxygenase-2 selective inhibitors: generation of predictive pharmacophore model using Catalyst

Pharmacophore-based virtual screening of ZINC database, molecular modeling and designing new derivatives as potential HDAC6 inhibitors

To date, many HDAC6 inhibitors have been identified and developed but none is clinically approved as of now. Through this study, we aim to obtain novel HDAC6 selective inhibitors and provide new insights into the detailed structural design of potential HDAC6 inhibitors. A HypoGen-based 3D QSAR HDAC6 pharmacophore was built and used as a query model to screen approximately 8 million ZINC database compounds. First, the ZINC Database was filtered using ADMET, followed by pharmacophore-based library screening. Using fit value and estimated activity cutoffs, a final set of 54 ZINC hits was obtained that were further investigated using molecular docking with the crystal structure of human histone deacetylase 6 catalytic domain 2 in complex with Trichostatin A (PDB ID: 5EDU). Through detailed in silico screening of the ZINC database, we shortlisted three hits as the lead molecules for designing novel HDAC6 inhibitors with better efficacy. Docking with 5EDU, followed by ADMET and TOPKAT analysis of modified ZINC hits provided 9 novel potential HDAC6 inhibitors that possess better docking scores and 2D interactions as compared to the control ZINC hit molecules. Finally, a 50 ns MD analysis run followed by Protein-Ligand Interaction Energy (PLIE) analysis of the top scored hits provided a novel molecule N1 that showed promisingly similar results to that of Ricolinostat (a known HDAC6 inhibitor). The comparable result of the designed hits to established HDAC6 inhibitors suggests that these compounds might prove to be successful HDAC6 inhibitors in future. Designed novel hits that might act as good HDAC6 inhibitors derived from ZINC database using combined molecular docking and modeling approaches.

LSD1-Based Reversible Inhibitors Virtual Screening and Binding Mechanism Computational Study

As one of the crucial targets of epigenetics, histone lysine-specific demethylase 1 (LSD1) is significant in the occurrence and development of various tumors. Although several irreversible covalent LSD1 inhibitors have entered clinical trials, the large size and polarity of the FAD-binding pocket and undesired toxicity have focused interest on developing reversible LSD1 inhibitors. In this study, targeting the substrate-binding pocket of LSD1, structure-based and ligand-based virtual screenings were adopted to expand the potential novel structures with molecular docking and pharmacophore model strategies, respectively. Through drug-likeness evaluation, ADMET screening, molecular dynamics simulations, and binding free energy screening, we screened out one and four hit compounds from the databases of 2,029,554 compounds, respectively. Generally, these hit compounds can be divided into two categories, amide (Lig2 and Comp2) and 1,2,4-triazolo-4,3-α-quinazoline (Comp3, Comp4, Comp7). Among them, Comp4 exhibits the strongest binding affinity. Finally, the binding mechanisms of the hit compounds were further calculated in detail by the residue free energy decomposition. It was found that van der Waals interactions contribute most to the binding, and FAD is also helpful in stabilizing the binding and avoiding off-target effects. We believe this work not only provides a solid theoretical foundation for the design of LSD1 substrate reversible inhibitors, but also expands the diversity of parent nucleus, offering new insights for synthetic chemists.

Application of per-residue energy decomposition to identify the set of amino acids critical for in silico prediction of COX-2 inhibitory activity

The enormous magnitude of scientific research carried out in the field of NSAIDs and cyclooxygenases (COXs) is known. They are crucial in pain management. COX-2 inhibitors have evolved over the years; from traditional NSAIDs to isoform-specific. The present study is aimed to identify a cluster of amino acids in the catalytic site whose energy contribution can better explain COX-2 inhibitory activity accurately than the binding energy of the whole protein. Initially, MD simulations (25 ns) and MM-PBSA calculations were performed for 8 diarylheterocyclic inhibitors. Per-residue energy decomposition studies were carried out to elucidate the energy contribution of each amino acid, and their correlation with COX-2 inhibitory activity was enumerated. A cluster of catalytic amino acids whose free energy sum has a high correlation with biological data was identified. The cluster of Gln178, Ser339, Tyr341, Arg499, Phe504, Val509 and Ala513 showed the correlation of -0.60. Further, the study was extended to a total of 26 COX-2 inhibitors belonging to different classes to validate the applicability of the cluster of amino acids identified. Results clearly suggest that the cluster of amino acids identified provide accurate screening method, and can be applied to predict COX-2 inhibitory activity of small molecules.

Identification of Novel TRPC5 Inhibitors by Pharmacophore-Based and Structure-Based Approaches

Canonical transient receptor potential-5 (TRPC5), which belongs to the subfamily of transient receptor potential (TRP) channels, is a non-selective cation channel mainly expressed in the central nervous system and shows more restricted expression in the periphery. TRPC5 plays a crucial role in human physiology and pathology, for instance, anxiety, depression, epilepsy, pain, memory and chronic kidney disease (CKD). However, due to lack of the effective and selective inhibitors, its physiological and pathological mechanism remains so far unknown. It is therefore pivotal to identify potential TRPC5 inhibitors. We have applied ligand-based virtual screening (LBVS) and structure-based virtual screening (SBVS) methods. The pharmacophore models of TRPC5 antagonists generated by using the HypoGen and HipHop algorithms were used as a query model for the screening of potential inhibitors against the Specs database. The resultant hits from LBVS were further screened by SBVS. SBVS was carried out based on the homology model generation of human TRPC5, binding site identification, molecular dynamics optimization and molecular docking studies. In our systematic screening approaches, we have identified 7 hits compounds with comparable dock score after Lipinski and Veber rules, ADMET, PAINS analysis, cluster analysis, and similarity analysis. In conclusion, the current research provides novel backbones for the new-generation of TRPC5 inhibitors.

Discovery of a potent cyclooxygenase-2 inhibitor, S4, through docking-based pharmacophore screening, in vivo and in vitro estimations

Cyclooxygenase (COX; EC: 1.14.99.1), the key enzyme in prostaglandin production in the human body, is a major pharmacological target for developing anti-inflammatory agents.

The effect of various atomic partial charge schemes to elucidate consensus activity-correlating molecular regions: a test case of diverse QSAR models

The estimation of atomic partial charges of the small molecules to calculate molecular interaction fields (MIFs) is an important process in field-based quantitative structure-activity relationship (QSAR). Several studies showed the influence of partial charge schemes that drastically affects the prediction accuracy of the QSAR model and focused on the selection of appropriate charge models that provide highest cross-validated correlation coefficient ([Formula: see text] or q(2)) to explain the variation in chemical structures against biological endpoints. This study shift this focus in a direction to understand the molecular regions deemed to explain SAR in various charge models and recognize a consensus picture of activity-correlating molecular regions. We selected eleven diverse dataset and developed MIF-based QSAR models using various charge schemes including Gasteiger-Marsili, Del Re, Merck Molecular Force Field, Hückel, Gasteiger-Hückel, and Pullman. The generalized resultant QSAR models were then compared with Open3DQSAR model to interpret the MIF descriptors decisively. We suggest the regions of activity contribution or optimization can be effectively determined by studying various charge-based models to understand SAR precisely.

Novel design strategy for checkpoint kinase 2 inhibitors using pharmacophore modeling, combinatorial fusion, and virtual screening

Checkpoint kinase 2 (Chk2) has a great effect on DNA-damage and plays an important role in response to DNA double-strand breaks and related lesions. In this study, we will concentrate on Chk2 and the purpose is to find the potential inhibitors by the pharmacophore hypotheses (PhModels), combinatorial fusion, and virtual screening techniques. Applying combinatorial fusion into PhModels and virtual screening techniques is a novel design strategy for drug design. We used combinatorial fusion to analyze the prediction results and then obtained the best correlation coefficient of the testing set (r test) with the value 0.816 by combining the Best(train)Best(test) and Fast(train)Fast(test) prediction results. The potential inhibitors were selected from NCI database by screening according to Best(train)Best(test) + Fast(train)Fast(test) prediction results and molecular docking with CDOCKER docking program. Finally, the selected compounds have high interaction energy between a ligand and a receptor. Through these approaches, 23 potential inhibitors for Chk2 are retrieved for further study.

Design checkpoint kinase 2 inhibitors by pharmacophore modeling and virtual screening techniques

Damage to DNA is caused by ionizing radiation, genotoxic chemicals or collapsed replication forks. When DNA is damaged or cells fail to respond, a mutation that is associated with breast or ovarian cancer may occur. Mammalian cells control and stabilize the genome using a cell cycle checkpoint to prevent damage to DNA or to repair damaged DNA. Checkpoint kinase 2 (Chk2) is one of the important kinases, which strongly affects DNA-damage and plays an important role in the response to the breakage of DNA double-strands and related lesions. Therefore, this study concerns Chk2. Its purpose is to find potential inhibitors using the pharmacophore hypotheses (PhModels) and virtual screening techniques. PhModels can identify inhibitors with high biological activities and virtual screening techniques are used to screen the database of the National Cancer Institute (NCI) to retrieve compounds that exhibit all of the pharmacophoric features of potential inhibitors with high interaction energy. Ten PhModels were generated using the HypoGen best algorithm. The established PhModel, Hypo01, was evaluated by performing a cost function analysis of its correlation coefficient (r), root mean square deviation (RMSD), cost difference, and configuration cost, with the values 0.955, 1.28, 192.51, and 16.07, respectively. The result of Fischer's cross-validation test for the Hypo01 model yielded a 95% confidence level, and the correlation coefficient of the testing set (rtest) had a best value of 0.81. The potential inhibitors were then chosen from the NCI database by Hypo01 model screening and molecular docking using the cdocker docking program. Finally, the selected compounds exhibited the identified pharmacophoric features and had a high interaction energy between the ligand and the receptor. Eighty-three potential inhibitors for Chk2 are retrieved for further study.

Theoretical evaluation of flotation performance of carboxyl hydroxamic acids with different number of polar groups on the surfaces of diaspore (010) and kaolinite (001)

The adsorption behaviors of three carboxyl hydroxamic acids on diaspore (010) and kaolinite (001) have been studied by density functional theory (DFT) and molecular dynamics (MD) method. The results indicated that carboxyl hydroxamic acids could adsorb on diaspore surface by ionic bonds and hydrogen bonds, and adsorb on kaolinite surface by hydrogen bonds. The models of carboxyl hydroxamic acids adsorbed on diaspore and kaolinite surfaces are proposed.

Identification of novel HIV 1--protease inhibitors: application of ligand and structure based pharmacophore mapping and virtual screening

A combined ligand and structure-based drug design approach provides a synergistic advantage over either methods performed individually. Present work bestows a good assembly of ligand and structure-based pharmacophore generation concept. Ligand-oriented study was accomplished by employing the HypoGen module of Catalyst in which we have translated the experimental findings into 3-D pharmacophore models by identifying key features (four point pharmacophore) necessary for interaction of the inhibitors with the active site of HIV-1 protease enzyme using a training set of 33 compounds belonging to the cyclic cyanoguanidines and cyclic urea derivatives. The most predictive pharmacophore model (hypothesis 1), consisting of four features, namely, two hydrogen bond acceptors and two hydrophobic, showed a correlation (r) of 0.90 and a root mean square of 0.71 and cost difference of 56.59 bits between null cost and fixed cost. The model was validated using CatScramble technique, internal and external test set prediction. In the second phase of our study, a structure-based five feature pharmacophore hypothesis was generated which signifies the importance of hydrogen bond donor, hydrogen bond acceptors and hydrophobic interaction between the HIV-1 protease enzyme and its inhibitors. This work has taken a significant step towards the full integration of ligand and structure-based drug design methodologies as pharmacophoric features retrieved from structure-based strategy complemented the features from ligand-based study hence proving the accuracy of the developed models. The ligand-based pharmacophore model was used in virtual screening of Maybridge and NCI compound database resulting in the identification of four structurally diverse druggable compounds with nM activities.

Development of energetic pharmacophore for the designing of 1,2,3,4-tetrahydropyrimidine derivatives as selective cyclooxygenase-2 inhibitors

We present here the Energetic pharmacophore model representing complementary features of the 1,2,3,4-tetrahydropyrimidine for selective cyclooxygenase-2 (COX-2) inhibition. For the development of pharmacophore hypothesis, a total of 43 previously reported compounds were docked on active site of COX-2 enzyme. The generated pharmacophore features were ranked using energetic terms of Glide XP docking for 1,2,3,4-tetrahydropyrimidine scaffold to optimize its structure requirement for COX-2 inhibition. The thirty new 4,5,6-triphenyl-1,2,3,4-tetrahydropyrimidine derivatives were synthesized and assessed for selective COX-2 inhibitory activity. Two compounds 4B1 and 4B11 were found to be potent and selective COX-2 inhibitors. The molecular docking studies revealed that the newly synthesized compounds can be docked into COX-2 binding site and also provide the molecular basis for their activity.

Design of novel FLT-3 inhibitors based on dual-layer 3D-QSAR model and fragment-based compounds in silico

FMS-like tyrosine kinase 3 (FLT-3) is strongly correlated with acute myeloid leukemia, but no FLT-3-inhibitor cocomplex structure is available to assist the design of therapeutic inhibitors. Hence, we propose a dual-layer 3D-QSAR model for FLT-3 that integrates the pharmacophore, CoMFA, and CoMSIA. We then coupled the model with the fragment-based design strategy to identify novel FLT-3 inhibitors. In the first layer, the previously established model, Hypo02, was evaluated in terms of its correlation coefficient (r), RMS, cost difference, and configuration cost, with values of 0.930, 1.24, 106.45, and 16.44, respectively. Moreover, Fischer's cross-validation test of data generated by Hypo02 yielded a 98% confidence level, and the validation of the testing set yielded a best r value of 0.87. The features of Hypo02 were separated into two parts and then used to screen the MiniMaybridge fragment compound database. Nine novel FLT-3 inhibitors were generated in this layer. In the second layer, Hypo02 was subjected to an alignment rule to generate CoMFA- and CoMSIA-based models, for which the partial least-squares validation method was utilized. The values of q(2), r(2), and predictive r(2) were 0.58, 0.98, and 0.76, respectively, derived from the CoMFA model with steric and electrostatic fields. The CoMSIA model with five different fields yielded values of 0.54, 0.97, and 0.76 for q(2), r(2), and predictive r(2), respectively. The CoMFA and CoMSIA models were used to constrain 3D structures of the nine novel FLT-3 inhibitors. This dual-layer 3D-QSAR model constitutes a valuable tool to easily and quickly screen and optimize novel potential FLT-3 inhibitors for the treatment of acute myeloid leukemia.

Pharmacophore modeling and virtual screening studies of checkpoint kinase 1 inhibitors

In this study, chemical feature-based 3-dimensional (3D) pharmacophore models of Checkpoint kinase 1 (Chk1) inhibitors were developed based on the known inhibitors of Chk1. The best pharmacophore model Hypo1 was characterized by the best correlation coefficient (0.9577), and the lowest root mean square deviation (0.8871). Hypo1 consists of one hydrogen-bond acceptor, one hydrogen-bond donor, and two hydrophobic features, as well as one excluded volume. This pharmacophore model was further validated by both test set and cross validation methods. A comparison analysis of Hypo1 with chemical features in the active site of Chk1 indicates that the pharmacophore model Hypo1 can correctly reflect the interactions between Chk1 and its ligands. Then Hypo1 was used to screen chemical databases, including Specs and Chinese Nature Product Database (CNPD) for potential lead compounds. The hit compounds were subsequently subjected to filtering by Lipinski's rule of five and docking study to refine the retrieved hits. Finally some of the most potent (estimated) compounds were selected from the final refined hits and suggested for further experimental investigation.