Pharmacophore modeling for protein tyrosine phosphatase 1B inhibitors

Hyphenated 3D-QSAR statistical model-drug repurposing analysis for the identification of potent neuraminidase inhibitor

The Influenza A virus is one of the principle causes of respiratory illness in human. The surface glycoprotein of the influenza virus, neuraminidase (NA), has a vital role in the release of new viral particle and spreads infection in the respiratory tract. It has been long recognized as a valid drug target for influenza A virus infection. Oseltamivir is used as a standard drug of choice for the treatment of influenza. However, the emergence of mutants with novel mutations has increased the resistance to potent NA inhibitor. In the present investigation, we have employed computer-assisted combinatorial techniques in the screening of 8621 molecules from Drug Bank to find potent NA inhibitors. A three-dimensional pharmacophore model was generated from the previously reported 28 carbocylic influenza NA inhibitors along with oseltamivir using PHASE module of Schrödinger Suite. The model generated consists of one hydrogen bond acceptor (A), one hydrogen bond donors (D), one hydrophobic group (H), and one positively charged group (P), ADHP. The hypothesis was further validated for its integrity and significance using enrichment analysis. Subsequently, an atom-based 3D-QSAR model was built using the common pharmacophore hypothesis (CPH). The developed 3D-QSAR model was found to be statistically significant with R² value of 0.9866 and Q² value of 0.7629. Further screening was accomplished using three-stage docking process using the Glide algorithm. The resultant lead molecules were examined for its drug-like properties using the Qikprop algorithm. Finally, the calculated pIC₅₀ values of the lead compounds were validated by the AutoQSAR algorithm. Overall, the results from our analysis highlights that lisinopril (DB00722) is predicted to bind better with NA than currently approved drug. In addition, it has the best match in binding geometry conformations with the existing NA inhibitor. Note that the antiviral activity of lisinopril is reported in the literature. However, our paper is the first report on lisinopril activity against influenza A virus infection. These results are envisioned to help design the novel NA inhibitors with an increased antiviral efficacy.

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.

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 PTP1B1-282-allosteric inhibitor complex crystal structure lacks alpha7 (287-298) and moreover there is no available 3D structure of PTP1B1-298 in open form. As the interaction between alpha7 and alpha6-alpha3 helices plays a crucial role in allosteric inhibition, alpha7 was modeled to the PTP1B1-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 alpha7-alpha6-alpha3 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 alpha7 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.