Sanjeev Kumar Singh

Extra precision docking, free energy calculation and molecular dynamics simulation studies of CDK2 inhibitors.

Molecular docking, free energy calculation and molecular dynamics (MD) simulation studies have been performed, to explore the putative binding modes of 3,5-diaminoindazoles, imidazo(1,2-b)pyridazines and triazolo(1,5-a) pyridazines series of Cyclin-dependent kinase (CDK2) inhibitors. To evaluate the effectiveness of docking protocol in flexible docking, we have selected crystallographic bound compound to validate our docking procedure as evident from root mean square deviations (RMSDs). We found different binding sites namely catalytic, inhibitory phosphorylation, cyclin binding and CKS-binding site of the CDK2 contributing towards the binding of these compounds. Moreover, correlation between free energy of binding and biological activity yielded a statistically significant correlation coefficient. Finally, three representative protein-ligand complexes were subjected to molecular dynamics simulation to determine the stability of the predicted conformations. The low value of the RMSDs between the initial complex structure and the energy minimized final average complex structure suggests that the derived docked complexes are close to equilibrium. We suggest that the phenylacetyl type of substituents and cyclohexyl moiety make the favorable interactions with a number of residues in the active site, and show better inhibitory activity to improve the pharmacokinetic profile of compounds against CDK2. The structure-based drug design strategy described in this study will be highly useful for the development of new inhibitors with high potency and selectivity.

In silico screening of indinavir-based compounds targeting proteolytic activity in HIV PR: binding pocket fit approach

The intense research on small molecule inhibitors of Human immunodeficiency virus (HIV)-protease (PR) has produced a diverse class of chemical scaffolds which includes clinically available HIV PR inhibitors (PRI). Till now, these inhibitors are insignificant for targeting proteolytic activity and few drug molecules on alterations can enhance the inhibition of PR enzyme. Here, we developed a method for screening of new hits from Cambridge structural database, based on binding mode of indinavir interaction participating atoms. Knowledge-based ligand screening technique approximately informs that new hits are also having same binding mode-like indinavir interaction patterns. Considering the importance of ligand fitting in binding pocket, we developed induced-fit models for each compound and we obtained accurate energy values in terms of binding and interaction energy. We found that newly search molecules are interacting better than known drug—indinavir and these new compounds are comparatively having better drug-like property. Finally, we demonstrated that pocket specific docking, energy utilization, interactions, and ADME for screened compounds are showing new hit compounds of indinavir are better HIV PRI and these new compounds can also show better activity in in vivo and in vitro conditions.

Validation of potential inhibitors for SrtA against Bacillus anthracis by combined approach of ligand-based and molecular dynamics simulation

The development of SrtA inhibitors targeting the biothreat organism namely Bacillus anthracis was achieved by the combined approach of pharmacophore modeling, binding interactions, electron transferring capacity, ADME, and Molecular dynamics studies. In this study, experimentally reported Ba-SrtA inhibitors (pyridazinone and pyrazolethione derivatives) were considered for the development of enhanced pharmacophoric model. The obtained AAAHR hypothesis was a pure theoretical concept that accounts for common molecular interaction network present in experimentally active pyridazinone and pyrazolethione derivatives. Pharmacophore-based screening of AAAHR hypothesis provides several new compounds, and those compounds were treated with four phases of docking protocols with combined Glide-QPLD docking approach. In this approach, scoring and charge accuracy variations were seen to be dominated by QM/MM approach through the allocation of partial charges. Finally, we reported the best compounds from binding db, Chembridge db, and Toslab based on scoring values, energy parameters, electron transfer reaction, ADME, and cell adhesion inhibition activity. The dynamic state of interaction and binding energy assess that new compounds are more active inside the binding pocket and these compounds on experimental validations will survive as better inhibitors for targeting the cell adhesion mechanism of Ba-SrtA.

Exploring the selectivity of a ligand complex with CDK2/CDK1: a molecular dynamics simulation approach

Cyclin‐dependent kinases (CDKs) are core components of the cell cycle machinery that govern the transition between phases during cell cycle progression. Abnormalities in CDKs activity and regulation are common features of cancer, making CDK family members attractive targets for the development of anticancer drugs. Their inhibitors have entered in clinical trials to treat cancer. Very recently, Heathcote et al. (J. Med. Chem. 2010, 53:8508–8522) have found a ligand BS194 that has a high affinity with CDK2 (IC50 = 3 nm) but shows low affinity with CDK1 (IC50 = 30 nm). To understand the selectivity, we used homology modeling, molecular docking, molecular dynamics, and free‐energy calculation to analyze the interactions. A rational three‐dimensional model of the CDK1/BS194 complex is built. We found that Leu83 is a key residue that recognizes BS194 more effectively with CDK2 with good binding free energies rather than CDK1. Energetic analysis reveals that van der Waals interaction and non‐polar contributions to solvent are favorable in the formation of complexes and amine group of the ligand, which plays a crucial role for binding selectivity between CDK2 and CDK1. Copyright

Pharmacophore modelling and atom-based 3D-QSAR studies on N-methyl pyrimidones as HIV-1 integrase inhibitors.

Pharmacophore modelling and atom-based 3D-QSAR studies were carried out for a series of compounds belonging to N-methyl pyrimidones as HIV-1 integrase inhibitors. Based on the ligand-based pharmacophore model, we got 5-point pharmacophore model AADDR, with two hydrogen bond acceptors (A), two hydrogen bond donors (D) and one aromatic ring (R). The generated pharmacophore-based alignment was used to derive a predictive atom-based 3D-QSAR model for the training set (r2 = 0.92, SD = 0.16, F = 84.8, N = 40) and for test set (Q2 = 0.71, RMSE = 0.06, Pearson R = 0.90, N = 10). From these results, AADDR pharmacophore feature was selected as best common pharmacophore hypothesis, and atom-based 3D-QSAR results also support the outcome by means of favourable and unfavourable regions of hydrophobic and electron-withdrawing groups for the most potent compound 30. These results can be useful for further design of new and potent HIV-1 IN inhibitors.

Exploration of fluoroquinolone resistance in Streptococcus pyogenes: comparative structure analysis of wild-type and mutant DNA gyrase.

Quinolone resistance‐determining region is known to be the druggability site of the target protein that undergoes frequent mutation and thus renders quinolone resistance. In the present study, ligands were tested for their inhibitory activity against DNA gyrase of Streptococcus pyogenes involved in DNA replication. In silico mutational analysis on modelled gyrase A revealed that GLU85 had the most possible interactions with all the ligands used for the study. The amino acid residue GLU85 had also been predicted with an essential role of maintaining the three‐dimensional structure of the protein. When introduced with a mutation (GLU 85 LYS) on this particular residue, it had readily denatured the whole α‐helix (from 80 to 90 amino acids). This was confirmed through the molecular dynamics simulation and revealed that this single mutation can cause many functional and structural changes. Furthermore, LYS85 mutation has altered the original secondary structure of the protein, which in turn led to the steric hindrance during the ligand–receptor interaction. The results based on the G‐score revealed that ligands have reduced interaction with the mutant protein. The semisynthetic fluoroquinolone 6d, which is an exception, forms a strong interaction with the mutant protein and was experimentally verified using the antimicrobial test. Hence, the present study unravels the fact that mutation at the drug binding site is the major cause for different level of resistance by the S. pyogenes when exposed against the varying concentrations of the fluoroquinolones. Furthermore, a comparative assessment of quinolone derivative with the older generation fluoroquinolones will be of great impact for S. pyogenes–related infections. Copyright

In silico and in vitro studies of cinnamaldehyde and their derivatives against LuxS in Streptococcus pyogenes: effects on biofilm and virulence genes.

The LuxS‐based signalling pathway has an important role in physiological and pathogenic functions that are capable of causing different infections. In the present study, cinnamaldehyde (CN) and their derivatives were evaluated for their inhibitory efficiency against LuxS by molecular modelling, docking, dynamics and free‐energy calculations. Sequence and structure‐similarity analysis of LuxS protein, five different amino acids were found to be highly conserved, of which GLY128 was identified as the key residue involved in the effective binding of the ligands. Quantum‐polarized ligand docking protocol showed that 2nitro and 4nitro CN has a higher binding efficiency than CN, which very well corroborates with the in vitro studies. COMSTAT analysis for the microscopic images of the S. pyogenes biofilm showed that the ligands have antibiofilm potential. In addition, the results of quantitative polymerase chain reaction (qPCR) analysis revealed that the transcripts treated with the compounds showed decrease in luxS expression, which directly reflects with the reduction in expression of speB. No substantial effect was observed on the virulence regulator (srv) transcript. These results confirm that speB is controlled by the regulation of luxS. The decreased rate of S. pyogenes survival in the presence of these ligands envisaged the fact that the compounds could readily enhance opsonophagocytosis with the reduction of virulence factor secretion. Thus, the overall data supports the use of CN derivatives against quorum sensing‐mediated infections caused by S. pyogenes. Copyright

A leishmaniasis study: Structure-based screening and molecular dynamics mechanistic analysis for discovering potent inhibitors of spermidine synthase☆

Protozoa Leishmania donovani (Ld) is the main cause of the endemic disease leishmaniasis. Spermidine synthase (SS), an important enzyme in the synthetic pathway of polyamines in Ld, is an essential element for the survival of this protozoan. Targeting SS may provide an important aid for the development of drugs against Ld. However, absence of tertiary structure of spermidine synthase of Leishmania donovani (LSS) limits the possibilities of structure based drug designing. Presence of the same enzyme in the host itself further challenges the drug development process. We modeled the tertiary structure of LSS using homology modeling approach making use of homologous X-ray crystallographic structure of spermidine synthase of Trypanosoma cruzi (TSS) (2.5Å resolution). The modeled structure was stabilized using Molecular Dynamics simulations. Based on active site structural differences between LSS and human spermidine synthase (HSS), we screened a large dataset of compounds against modeled protein using Glide virtual screen docking and selected two best inhibitors based on their docking scores (-10.04 and -13.11 respectively) with LSS and having least/no binding with the human enzyme. Finally Molecular Dynamics simulations were used to assess the dynamic stability of the ligand bound structures and to elaborate on the binding modes. This article is part of a Special Issue entitled: Computational Methods for Protein Interaction and Structural Prediction.

3D-QSAR CoMFA and CoMSIA Study on Benzodipyrazoles as Cyclin Dependent Kinase 2 Inhibitors

Cyclin dependent kinase 2 (CDK2) has appeared as an important drug target over the years with a multitude of therapeutic potentials. With the intention of designing compounds with enhanced inhibitory potencies against CDK2, the 3D-QSAR CoMFA and CoMSIA study on benzodipyrazoles series is presented here. The developed models showed a strong correlative and predictive capability having a cross validated correlation co-efficient of (r(2)(cv)) 0.699 for CoMFA and 0.794 for CoMSIA models. A very good conventional and predicted correlation co-efficients were also obtained: CoMFA (r(2)(ncv), r(2)(pred): 0.883, 0.754), CoMSIA (0.937, 0.815). The models were found to be statistically robust and are expected to be of an aid to design and/or prioritize drug likes for synthesis.

Assessment of dual inhibition property of newly discovered inhibitors against PCAF and GCN5 through in silico screening, molecular dynamics simulation and DFT approach.

Abstract p300/CBP-associated factor (PCAF) is one among the histone acetyltransferase (HAT) family enzymes. It is involved in the regulation of transcription by modifying the chromatin structure indirectly through the acetylation of histones. It has been emerged as a promising drug target for various types of cancer. A four-point pharmacophore with two hydrogen bond acceptor, one aromatic ring and one hydrophobic feature, was generated for six highly active isothiazolone derivatives as PCAF inhibitors in order to elucidate their anticancer activity. The generated pharmacophore was used for screening three different databases such as Maybridge, Life Chemicals and Chembridge databases. The screened compounds were further filtered through docking studies. Then the compounds were further carried for ADME prediction. The best three compounds BTB09406, F1418-0051 and F1880-1727 were docked to GCN5 to explore the dual inhibitory properties. The conformational stability of the protein–ligand complexes were analyzed through molecular dynamics simulation. Three best compounds were finally went through electronic structure analysis using density functional theory (DFT) at B3LYP/6-31**G level to understand their molecular reactivity. The results obtained from this study exploit that the three best compounds (BTB09406, F1418-0051 and F1880-1727) were found to have more potent and dual inhibitory properties.