Rohit Farmer

A conserved motif flags acyl carrier proteins for β-branching in polyketide synthesis.

Type I PKSs often utilise programmed β-branching, via enzymes of an “HMG-CoA synthase (HCS) cassette”, to incorporate various side chains at the second carbon from the terminal carboxylic acid of growing polyketide backbones. We identified a strong sequence motif in Acyl Carrier Proteins (ACPs) where β-branching is known. Substituting ACPs confirmed a correlation of ACP type with β-branching specificity. While these ACPs often occur in tandem, NMR analysis of tandem β-branching ACPs indicated no ACP-ACP synergistic effects and revealed that the conserved sequence motif forms an internal core rather than an exposed patch. Modelling and mutagenesis identified ACP Helix III as a probable anchor point of the ACP-HCS complex whose position is determined by the core. Mutating the core affects ACP functionality while ACP-HCS interface substitutions modulate system specificity. Our method for predicting β-carbon branching expands the potential for engineering novel polyketides and lays a basis for determining specificity rules.

Binding efficiencies of carbohydrate ligands with different genotypes of cholera toxin B: molecular modeling, dynamics and docking simulation studies.

Vibrio cholerae produces cholera toxin (CT) that consists of two subunits, A and B, and is encoded by a filamentous phage CTXΦ. The A subunit carries enzymatic activity that ribosylates ADP, whereas the B subunit binds to monosialoganglioside (GM1) receptor in epithelial cells. Molecular analysis of toxigenic V. cholerae strains indicated the presence of multiple ctxB genotypes. In this study, we employed a comparative modeling approach to define the structural features of all known variants of ctxB found in O139 serogroup V. cholerae. Modeling, molecular dynamics and docking simulations studies suggested subtle variations in the binding ability of ctxB variants to carbohydrate ligands of GM1 (galactose, sialic acid and N-acetyl galactosamine). These findings throw light on the molecular efficiencies of pathogenic isolates of V. cholerae harboring natural variants of ctxB in causing the disease, thus suggesting the need to consider ctxB variations when designing vaccines against cholera.

In silico Prediction and Characterization of 3D Structure and Binding Properties of Catalase from the Commercially Important Crab, Scylla serrata

The enzyme catalase breaks down H2O2, a potentially harmful oxidant, to H2O and O2. Besides oxidase activity, the enzyme also exhibits peroxidase activity. Therefore, it plays an important role in maintaining health and regulating pathophysiology of the organisms. However, 3D structure of this important enzyme in invertebrates particularly in crabs is not yet available. Therefore, an attempt has been made to predict the structure of the crab catalase and to envisage its catalytic interaction with H2O2. A three dimensional model of crab catalase was constructed using the NADPH binding site on Beef Liver catalase from Bos taurus (PDBID: 7CAT) as template by comparative modeling approach. Backbone conformation of the modeled structure by PROCHECK revealed that more than 98% of the residues fell in the allowed regions, ERRAT results confirmed good quality of modeled structure and VERIFY3D profile was satisfying. Molecular docking has been used to know the binding modes of hydrogen peroxide with the crab catalase protein. The receptor structures used for docking were derived from molecular dynamics (MD) simulations of homology modeled structure. The docking results showed that the three important determinant residues Arg68, Val70 and Arg108 in catalase were binding with H2O2 as they had strong hydrogen bonding contacts with the substrate. Our analysis provides insight into the structural properties of crab catalase and defines its active sites for binding with substrate. These data are important for further studies of catalase of invertebrates in general and that of crabs in particular.

Metabolic pathway analysis and molecular docking analysis for identification of putative drug targets in Toxoplasma gondii: novel approach.

Toxoplasma gondii is an obligate intracellular apicomplexan parasite that can infect a wide range of warm-blooded animals including humans. In humans and other intermediate hosts, toxoplasma develops into chronic infection that cannot be eliminated by host’s immune response or by currently used drugs. In most cases, chronic infections are largely asymptomatic unless the host becomes immune compromised. Thus, toxoplasma is a global health problem and the situation has become more precarious due to the advent of HIV infections and poor toleration of drugs used to treat toxoplasma infection, having severe side effects and also resistance have been developed to the current generation of drugs. The emergence of these drug resistant varieties of T. gondii has led to a search for novel drug targets. We have performed a comparative analysis of metabolic pathways of the host Homo sapiens and the pathogen T. gondii. The enzymes in the unique pathways of T. gondii, which do not show similarity to any protein from the host, represent attractive potential drug targets. We have listed out 11 such potential drug targets which are playing some important work in more than one pathway. Out of these, one important target is Glutamate dehydrogenase enzyme; it plays crucial part in oxidation reduction, metabolic process and amino acid metabolic process. As this is also present in the targets of tropical diseases of TDR (Tropical disease related Drug) target database and no PDB and MODBASE 3D structural model is available, homology models for Glutamate dehydrogenase enzyme were generated using MODELLER9v6. The model was further explored for the molecular dynamics simulation study with GROMACS, virtual screening and docking studies with suitable inhibitors against the NCI diversity subset molecules from ZINC database, by using AutoDock-Vina. The best ten docking solutions were selected (ZINC01690699, ZINC17465979, ZINC17465983, ZINC18141294_03, ZINC05462670, ZINC01572309, ZINC18055497_01, ZINC18141294, ZINC05462674 and ZINC13152284_01). Further the Complexes were analyzed through LIGPLOT. On the basis of Complex scoring and binding ability it is deciphered that these NCI diversity set II compounds, specifically ZINC01690699 (as it has minimum energy score and one of the highest number of interactions with the active site residue), could be promising inhibitors for T. gondii using Glutamate dehydrogenase as Drug target.

Selective AKR1C3 inhibitors do not recapitulate the anti-leukaemic activities of the pan-AKR1C inhibitor medroxyprogesterone acetate

Background:We and others have identified the aldo-keto reductase AKR1C3 as a potential drug target in prostate cancer, breast cancer and leukaemia. As a consequence, significant effort is being invested in the development of AKR1C3-selective inhibitors.Methods:We report the screening of an in-house drug library to identify known drugs that selectively inhibit AKR1C3 over the closely related isoforms AKR1C1, 1C2 and 1C4. This screen initially identified tetracycline as a potential AKR1C3-selective inhibitor. However, mass spectrometry and nuclear magnetic resonance studies identified that the active agent was a novel breakdown product (4-methyl(de-dimethylamine)-tetracycline (4-MDDT)).Results:We demonstrate that, although 4-MDDT enters AML cells and inhibits their AKR1C3 activity, it does not recapitulate the anti-leukaemic actions of the pan-AKR1C inhibitor medroxyprogesterone acetate (MPA). Screens of the NCI diversity set and an independently curated small-molecule library identified several additional AKR1C3-selective inhibitors, none of which had the expected anti-leukaemic activity. However, a pan AKR1C, also identified in the NCI diversity set faithfully recapitulated the actions of MPA.Conclusions:In summary, we have identified a novel tetracycline-derived product that provides an excellent lead structure with proven drug-like qualities for the development of AKR1C3 inhibitors. However, our findings suggest that, at least in leukaemia, selective inhibition of AKR1C3 is insufficient to elicit an anticancer effect and that multiple AKR1C inhibition may be required.

In silico prediction of 3D structure of Mn superoxide dismutase of Scylla serrata and its binding properties with inhibitors

In the present study, we used computational methods to model crab and rat MnSOD using the crystal structure of MnSOD from Homo sapiens (PDB code: 1MSD) as template by comparative modeling approach. We performed molecular dynamics simulations to study dynamic behavior of the crab MnSOD. The modeled proteins were validated and subjected to molecular docking analyses. Molecular docking tool was used to elucidate a comparative binding mode of the crab and rat SOD with potent inhibitors of SOD such as hydrogen peroxide (H2O2), potassium cyanide (KCN) and sodium dodecyl sulphate (SDS). The predicted valid structure of crab MnSOD did not show any interaction with KCN but close interaction with H2O2 and SDS. A possible inhibitory mechanism of SDS and H2O2 due to their interaction with the amino acids present in the active site of the MnSOD of the above two animals are elucidated. This allowed us to predict the binding modes of the proteins to elucidate probable mode of action and sites of interference.

Prediction and analysis of paralogous proteins in Trichomonas vaginalis genome.

Trichomonas vaginalis causes trichomoniasis, second most sexually transmitted disease. The genome sequence draft of T. vaginalis was published by The Institute of Genomic Research reveals an abnormally large genome size of 160 Mb. It was speculated that a significant portion of the proteome contains paralogous proteins. The present study was aimed at identification and analysis of the paralogous proteins. The all against all search approach is used to identify the paralogous proteins. The dataset of proteins was retrieved from TIGR and TrichDB FTP server. The BLAST-P program performed all against all database searches against the protein database of Trichomonas vaginalis available at NCBI genome database. In the present study about 50,000 proteins were searched where 2,700 proteins were found to be paralogous under the rigid selection criteria. The Pfam database search has identified significant number of paralogous proteins which were further categorized among different 1496 paralogous protein in pfam families, 1027 paralogous protein contains domain, 60 proteins were having different repeats and 1092 paralogous protein sequences of clans. Such identification and functional annotation of paralogous proteins will also help in removing paralogous proteins from possible drug targets in future. Presence of huge number of paralogous proteins across wide range of gene families and domains may be one of the possible mechanisms involved in the T. vaginalis genome expansion and evolution.

Virtual screening of AmpC/β‐lactamase as target for antimicrobial resistance in Pseudomonas aeruginosa

AmpC is a group I, class C ‐lactamase present in most Enterobacteriaceae and in Pseudomonas aeruginosa and other nonfermenting gram-negative bacilli. The β‐lactam class of antibiotics is one of the most important structural classes of antibacterial compounds and act by inhibiting the bacterial D ,D - transpeptidases that are responsible for the final step of peptidoglycan cross-linking. Our main aim in the study is to screen possible inhibitors against AmpC / β ‐ lactamase (an enzyme responsible for antimicrobial activity in Pseudomonas aeruginosa), through virtual screening of 1364 NCI (National Cancer Institute) diversity set II compounds. Homology Model of AmpC / β ‐ lactamase was constructed using MODELLER and the Model was validated using PROCHECK and Verify 3D programs to obtain a stable structure, which was further used for virtual screening of NCI (National Cancer Institute) diversity set II compounds through molecular Docking studies using Autodock. The amino acid sequence of the β ‐ lactamase was also subjected to ScanProsite web server to find any pattern present in the sequence. After the prediction of 3-dimensional model of AmpC/ β‐lactamase, the possible Active sites ofβ ‐ lactamase were determined using LIGSITEcsc and CastP web servers simultaneously. The Docked complexes were validated and Enumerated based on the Autodock Scoring function to pick out the best inhibitor based on Autodock energy score. Thus from the entire 1364 NCI diversity set II compounds which were Docked, the best four docking solutions were selected (ZINC12670903, ZINC17465965, ZINC11681166 and ZINC13099024). Further the Complexes were analyzed through LIGPLOT for their interaction for the 4 best docked NCI diversity set II compounds. Thus from the Complex scoring and binding ability it is deciphered that these NCI diversity set II compounds could be promising inhibitors for Pseudomonas aeruginosa using AmpC /β ‐ lactamase as Drug target yet pharmacological studies have to confirm it.

Molecular Dynamic Simulation and Inhibitor Prediction of Cysteine Synthase Structured Model as a Potential Drug Target for Trichomoniasis

In our presented research, we made an attempt to predict the 3D model for cysteine synthase (A2GMG5_TRIVA) using homology-modeling approaches. To investigate deeper into the predicted structure, we further performed a molecular dynamics simulation for 10 ns and calculated several supporting analysis for structural properties such as RMSF, radius of gyration, and the total energy calculation to support the predicted structured model of cysteine synthase. The present findings led us to conclude that the proposed model is stereochemically stable. The overall PROCHECK G factor for the homology-modeled structure was −0.04. On the basis of the virtual screening for cysteine synthase against the NCI subset II molecule, we present the molecule 1-N, 4-N-bis [3-(1H-benzimidazol-2-yl) phenyl] benzene-1,4-dicarboxamide (ZINC01690699) having the minimum energy score (−13.0 Kcal/Mol) and a log P value of 6 as a potential inhibitory molecule used to inhibit the growth of T. vaginalis infection.

Structural modelling and molecular dynamics of a multi-stress responsive WRKY TF-DNA complex towards elucidating its role in stress signalling mechanisms in chickpea

Chickpea is a premier food legume crop with high nutritional quality and attains prime importance in the current era of 795 million people being undernourished worldwide. Chickpea production encounters setbacks due to various stresses and understanding the role of key transcription factors (TFs) involved in multiple stresses becomes inevitable. We have recently identified a multi-stress responsive WRKY TF in chickpea. The present study was conducted to predict the structure of WRKY TF to identify the DNA-interacting residues and decipher DNA-protein interactions. Comparative modelling approach produced 3D model of the WRKY TF with good stereochemistry, local/global quality and further revealed W19, R20, K21, and Y22 motifs within a vicinity of 5 Å to the DNA amongst R18, G23, Q24, K25, Y36, Y37, R38 and K47 and these positions were equivalent to the 2LEX WRKY domain of Arabidopsis. Molecular simulations analysis of reference protein -PDB ID 2LEX, along with Car-WRKY TF modelled structure with the DNA coordinates derived from PDB ID 2LEX and docked using HADDOCK were executed. Root Mean Square (RMS) Deviation and RMS Fluctuation values yielded consistently stable trajectories over 50 ns simulation. Strengthening the obtained results, neither radius of gyration, distance and total energy showed any signs of DNA-WRKY complex falling apart nor any significant dissociation event over 50 ns run. Therefore, the study provides first insights into the structural properties of multi-stress responsive WRKY TF-DNA complex in chickpea, enabling genome wide identification of TF binding sites and thereby deciphers their gene regulatory networks.