Fragment based group QSAR and molecular dynamics mechanistic studies on arylthioindole derivatives targeting the α-β interfacial site of human tubulin

Simultaneous Pharmacokinetics Estimation of Nateglinide and Pioglitazone by RP-HPLC: Computational Study to Unlock the Synergism

Nateglinide (NAT) and Pioglitazone (PIO) are an antidiabetic drugs combination and currently under clinical trial in countries like Japan. In this study, an alternative, a simple, sensitive high-performance liquid chromatography method has been developed (limit of detection: 15 ng/mL and limit of quantification: 50 ng/mL) for simultaneous estimation of this drug combination in rat plasma. Most remarkably, bioavailability of NAT has been increased markedly on coadministration with PIO, than when it was administered alone. Thus, PIO is assumed to retard the catabolism of NAT by inhibiting metabolic liver-microsomal enzyme, especially CYP2C9. Using a Waters Nova-Pak C 18 column (150 × 3.9 mm, 4 μm) and a mobile phase of acetonitrile: 10 mM KH2PO4 (60: 40, V/V (volume by volume)) pH 3.5, the analysis was performed at 210 nm with a flow rate of 1.5 mL/min. In silico docking via molecular dynamics simulation revealed that NAT-CYP2C9 binding affinity may be reduced after PIO attachment, presumably due to the binding site overlapping of the two drugs. Thus, it has been proposed that NAT and PIO may be an efficient synergistic fixed dose combination against diabetes mellitus, and the above method can foster a simple but highly sensitive bioanalytical estimation for routine analysis.

Molecular modeling and lead design of substituted zanamivir derivatives as potent anti-influenza drugs

BACKGROUND

Influenza virus spreads infection by two main surface glycoproteins, namely hemagglutinin (HA) and neuraminidase (NA). NA cleaves the sialic acid receptors eventually releasing newly formed virus particles which then invade new cells. Inhibition of NA could limit the replication of virus to one round which is insufficient to cause the disease.

RESULTS

An experimentally reported series of acylguanidine zanamivir derivatives was used to develop GQSAR model targeting NA in different strains of influenza virus, H1N1 and H3N2. A combinatorial library was developed and their inhibitory activities were predicted using the GQSAR model.

CONCLUSION

The top leads were analyzed by docking which revealed the binding modes of these inhibitors in the active site of NA (150-loop). The top compound (AMA) was selected for carrying out molecular dynamics simulations for 15 ns which provided insights into the time dependent dynamics of the designed leads. AMA possessed a docking score of -8.26 Kcal/mol with H1N1 strain and -7.00 Kcal/mol with H3N2 strain. Ligand-bound complexes of both H1N1 and H3N2 were observed to be stable for 11 ns and 7 ns respectively. ADME descriptors were also calculated to study the pharmacokinetic properties of AMA which revealed its drug-like properties.

Mechanistic insights into mode of action of rice allene oxide synthase on hydroxyperoxides: An intermediate step in herbivory-induced jasmonate pathway

Various types of oxygenated fatty acids termed 'oxylipins' are involved in plant response to herbivory. Oxylipins like jasmonic acid (JA) and green leafy volatiles (GLVs) are formed by the action of enzymes like allene oxide synthase (AOS) and hydroxyperoxide lyase (HPL) respectively. In this study, we focus on AOS of Oryza sativa sb. Japonica, that interact with 9- and 13- hydroxyperoxides to produce intermediates of jasmonate pathway and compare it with rice HPL that yields GLVs. We attempt to elucidate the interaction pattern by computational docking protocols keeping the Arabidopsis AOS system as the reference model system. Both 9-hydroxyperoxide and 13-hydroxyperoxide fit into the active site of AOS completely with Phe347, Phe92, Ile463, Val345, and Asn278 being the common interacting residues. Phe347 and Phe92 were mutated with Leucine and docked again with the hydroxyperoxides. The Phe347→Leu347 mutant showed a different mode of action than AOS-hydroxyperoxide complex with Trp413 in direct bonding with the OOH group of 9-hydroxyperoxide. The loss of Lys88-OOH interaction in 13-hydroxyperoxide and loss-of-interaction of Leu347 indicated the importance of Phe347 residue in hydroxyperoxide catalysis. The second mutant Phe92→Leu92 also shows a very different interaction pattern with 13-hydroxyperoxide but not with 9-hydroxyperoxide.Therefore, it can be concluded that Phe347 is more crucial for AOS functionality than Phe92. The aromatic ring of a Phenylalanine residue is important for catalysis and its mutation affects the binding of the two ligands. Another important residue is Asn278 which is an important part of the AOS catalytic site for maintaining stability and can be compared with the Arabidopsis AOS residue Asn321. Lastly, the interaction of HPL with these two derivatives involves Leu363 residue instead of Phe347 and thus, validating the importance of Phe→Leu substitution to be the reason of different modes of action that result in completely different products from same substrates.

Molecular dynamics-based identification of novel natural mortalin-p53 abrogators as anticancer agents

INTRODUCTION

Cancer is one of the leading causes of mortality worldwide that requires attention in terms of extensive study and research. Eradication of mortalin-p53 interaction that leads to the inhibition of transcriptional activation or blocking of p53 from functioning as a suppressor and induction of nuclear translocation of p53 can prove to be one of the useful approaches for cancer management.

RESULTS

In this study, we used structure-based approach to target the p53-binding domain of mortalin in order to prevent mortalin-p53 complex formation. We screened compounds from ZINC database against the modeled mortalin protein using Glide virtual screening. The top two compounds, DTOM (ZINC 28639308) and TTOM (ZINC 38143676) with Glide score of -12.27 and -12.16, respectively, were identified with the potential to abrogate mortalin-p53 interaction. Finally, molecular dynamics simulations were used to analyze the dynamic stability of the ligand-bound complex and it was observed that residues Tyr196, Asn198, Val264 and Thr267 were involved in intermolecular interactions in both the simulated ligand-bound complexes, and thus, these residues may have a paramount role in stabilizing the binding of the ligands with the protein.

CONCLUSION

These detailed insights can further facilitate the development of potent inhibitors against mortalin-p53 complex.

Fragment based G-QSAR and molecular dynamics based mechanistic simulations into hydroxamic-based HDAC inhibitors against spinocerebellar ataxia

Expansion of polyglutamine (CAG) triplets within the coding gene ataxin 2 results in transcriptional repression, forming the molecular basis of the neurodegenerative disorder named spinocerebellar ataxia type-2 (SCA2). HDAC inhibitors (HDACi) have been elements of great interest in polyglutamine disorders such as Huntington's and Ataxia's. In this study, we have selected hydroxamic acid derivatives as HDACi and performed fragment-based G-QSAR, molecular docking studies and molecular dynamics simulations for elucidating the dynamic mode of action of HDACi with His-Asp catalytic dyad of HDAC4. The model was statistically validated to establish its predictive robustness. The model was statistically significant with r(2) value of .6297, cross-validated co-relation coefficient q(2) value of .5905 and pred_r(2) (predicted square co-relation coefficient) value of .85. An F-test value of 56.11 confirms absolute robustness of the model. Two combinatorial libraries comprising of 3180 compounds were created with hydroxamate moiety as the template and their pIC50 activities were predicted based on the G-QSAR model. The combinatorial library created was screened on the basis of predicted activity (pIC50), with two resultant top scoring compounds, HIC and DHC. The interaction of the compounds with His-Asp dyad in terms of H-bond interactions with His802, Asp840, Pro942, and Gly975 residues of HDAC4 was evaluated by docking and 20 ns long molecular dynamics simulations. This study provides valuable leads for structural substitutions required for hydroxamate moiety to exhibit enhanced inhibitory activity against HDAC4. The reported compounds demonstrated good binding and thus can be considered as potent therapeutic leads against ataxia.

Molecular principles behind Boceprevir resistance due to mutations in hepatitis C NS3/4A protease

The hepatitis C virus (HCV) infection is a primary cause of chronic hepatitis which eventually progresses to cirrhosis and in some instances might advance to hepatocellular carcinoma. According to the WHO report, HCV infects 130-150 million people globally and every year 350,000 to 500,000 people die from hepatitis C virus infection. Great achievement has been made in viral treatment evolution, after the development of HCV NS3/4A protease inhibitor (Boceprevir). However, efficacy of Boceprevir is compromised by the emergence of drug resistant variants. The molecular principle behind drug resistance of the protease mutants such as (V36M, T54S and R155K) is still poorly understood. Therefore in this study, we employed a series of computational strategies to analyze the binding of antiviral drug, Boceprevir to HCV NS3/4A protease mutants. Our results clearly demonstrate that the point mutations (V36M, T54S and R155K) in protease are associated with lowering of its binding affinity with Boceprevir. Exhaustive analysis of the simulated Boceprevir-bound wild and mutant complexes revealed variations in hydrophobic interactions, hydrogen bond occupancy and salt bridge interactions. Also, substrate envelope analysis scrutinized that the studied mutations reside outside the substrate envelope which may affect the Boceprevir affinity towards HCV protease but not the protease enzymatic activity. Furthermore, structural analyses of the binding site volume and flexibility show impairment in flexibility and stability of the binding site residues in mutant structures. In order to combat Boceprevir resistance, renovation of binding interactions between the drug and protease may be valuable. The structural insight from this study reveals the mechanism of the Boceprevir resistance and the results can be valuable for the design of new PIs with improved efficiency.

Structural investigations of T854A mutation in EGFR and identification of novel inhibitors using structure activity relationships

BACKGROUND

The epidermal growth factor receptor (EGFR) is a member of the ErbB family that is involved in a number of processes responsible for cancer development and progression such as angiogenesis, apoptosis, cell proliferation and metastatic spread. Malfunction in activation of protein tyrosine kinases has been shown to result in uncontrolled cell growth. The EGFR TK domain has been identified as suitable target in cancer therapy and tyrosine kinase inhibitors such as erlotinib have been used for treatment of cancer. Mutations in the region of the EGFR gene encoding the tyrosine kinase (TK) domain causes altered responses to EGFR TK inhibitors (TKI). In this paper we perform molecular dynamics simulations and PCA analysis on wild-type and mutant (T854A) structures to gain insight into the structural changes observed in the target protein upon mutation. We also report two novel inhibitors identified by combined approach of QSAR model development.

RESULTS

The wild-type and mutant structure was observed to be stable for 26 ns and 24 ns respectively. In PCA analysis, the mutant structure proved to be more flexible than wild-type. We developed a 3D-QSAR model using 38 thiazolyl-pyrazoline compounds which was later used for prediction of inhibitory activity of natural compounds of ZINC library. The 3D-QSAR model was proved to be robust by the statistical parameters such as r2 (0.9751), q2(0.9491) and pred_r2(0.9525).

CONCLUSION

Analysis of molecular dynamics simulations results indicate stability loss and increased flexibility in the mutant structure. This flexibility results in structural changes which render the mutant protein drug resistant against erlotinib. We report two novel compounds having high predicted inhibitory activity to EGFR TK domain with both wild-type and mutant structure.

InCoB2014: mining biological data from genomics for transforming industry and health

The 13th International Conference on Bioinformatics (InCoB2014) was held for the first time in Australia, at Sydney, July 31-2 August, 2014. InCoB is the annual scientific gathering of the Asia-Pacific Bioinformatics Network (APBioNet), hosted since 2002 in the Asia-Pacific region. Of 106 full papers submitted to the BMC track of InCoB2014, 50 (47.2%) were accepted in BMC Bioinformatics, BMC Genomics and BMC Systems Biology supplements, with three papers in a new BMC Medical Genomics supplement. While the majority of presenters and authors were from Asia and Australia, the increasing number of US and European conference attendees augurs well for the international flavour of InCoB. Next year's InCoB will be held jointly with the Genome Informatics Workshop (GIW), September 9-11, 2015 in Tokyo, Japan, with a view to integrate bioinformatics communities in the region.