Computational studies on Begomoviral AC2/C2 proteins.

Geminiviridae is a large family of circular, single stranded DNA viruses, which infects and causes devastating diseases on economically important crops. They are subdivided into nine genera. Members of the genus begomovirus encode a pathogenic protein called AC2/C2 which interacts that inactivates many plant proteins and trans-activates a number of host genes via the C-terminal transactivation domain. Hence, a sequence analysis on C-terminal region of AC2/C2 was completed. Analysis of 124 bipartite and 463 mono partite begomo viral AC2/C2 proteins revealed major differences in protein length, composition and position of acidic, aromatic and hydrophobic residues. Secondary structure analysis of AC2/C2 revealed the possible formation of C-terminal α-helix, which is similar to the acidic activation domain of many transcriptional activator proteins. Previous studies demonstrated that AC2 utilizes conserved late element (CLE) for the transactivation of viral genes and genome-wide mapping of same consensus in A. thaliana yielded 122 promoters with exact CLE consensus sequence. Analysis of protein interaction network for 106 CLE containing genes, 87 AC2 trans activated genes and 10 AC2 interacting proteins revealed a possible regulation of hundreds of host proteins which helps begomoviruses to produce a successful viral infection.

Computational protein design and protein-ligand interaction studies for the improvement of organophosphorus degrading potential of Deinococcus radiodurans.

The organophosphorus hydrolase enzyme is involved in the catalyzing reaction that involve hydrolysis of organophosphate toxic compounds. An enzyme from Deinococcus radiodurans reported as homologous to phosphotriesterase and show activity against organophosphate. In the past activity of this enzyme is low and efforts made to improve the activity by experimental mutation study. However only very few organophosphates tested against very few catalytic site mutations. In order to improve the catalytic power of the organophosphorus hydrolase enzyme, we carried out systematic functional hotspot based protein engineering strategy. The mutants tested against 46 know organophosphate compounds using molecular docking study. Finally, we carried out an extensive molecular docking study to predict the binding of 46 organophosphate compounds to wild-type protein and mutant organophosphorus hydrolase enzyme. At the end we are able to improve the degrading potential of organophosphorus hydrolase enzyme against organophosphate toxic compounds. This preliminary study and the outcome would be useful guide for the experimental scientist involved in the bioremediation of toxic organophosphate compounds.

Fragment-based virtual screening approach and molecular dynamics simulation studies for identification of BACE1 inhibitor leads.

Traditional structure-based virtual screening method to identify drug-like small molecules for BACE1 is so far unsuccessful. Location of BACE1, poor Blood Brain Barrier permeability and P-glycoprotein (Pgp) susceptibility of the inhibitors make it even more difficult. Fragment-based drug design method is suitable for efficient optimization of initial hit molecules for target like BACE1. We have developed a fragment-based virtual screening approach to identify/optimize the fragment molecules as a starting point. This method combines the shape, electrostatic, and pharmacophoric features of known fragment molecules, bound to protein conjugate crystal structure, and aims to identify both chemically and energetically feasible small fragment ligands that bind to BACE1 active site. The two top-ranked fragment hits were subjected for a 53 ns MD simulation. Principle component analysis and free energy landscape analysis reveal that the new ligands show the characteristic features of established BACE1 inhibitors. The potent method employed in this study may serve for the development of potential lead molecules for BACE1-directed Alzheimer's disease therapeutics.

Computational analysis of BACE1-ligand complex crystal structures and linear discriminant analysis for identification of BACE1 inhibitors with anti P-glycoprotein binding property.

More than 100 years of research on Alzheimer's disease didn't yield a potential cure for this dreadful disease. Poor Blood Brain Barrier (BBB) permeability and P-glycoprotein binding of BACE1 inhibitors are the major causes for the failure of these molecules during clinical trials. The design of BACE1 inhibitors with a balance of sufficient affinity to the binding site and little or no interaction with P-glycoproteins is indispensable. Identification and understanding of protein-ligand interactions are essential for ligand optimization process. Structure-based drug design (SBDD) efforts led to a steady accumulation of BACE1-ligand crystal complexes in the PDB. This study focuses on analyses of 153 BACE1-ligand complexes for the direct contacts (hydrogen bonds and weak interactions) observed between protein and ligand and indirect contacts (water-mediated hydrogen bonds), observed in BACE1-ligand complex crystal structures. Intraligand hydrogen bonds were analyzed, with focus on ligand P-glycoprotein efflux. The interactions are dissected specific to subsites in the active site and discussed. The observed protein-ligand and intraligand interactions were used to develop the linear discriminant model for the identification of BACE1 inhibitors with less or no P-glycoprotein binding property. Excellent statistical results and model's ability to correctly predict a new data-set with an accuracy of 92% is achieved. The results are retrospectively analyzed to give input for the design of potential BACE1 inhibitors.

Stigmatellin Y - An anti-biofilm compound from Bacillus subtilis BR4 possibly interferes in PQS-PqsR mediated quorum sensing system in Pseudomonas aeruginosa.

Hitherto this is the first report pertaining to production of biofilm inhibitory compound(s) (BIC) from Bacillus subtilis BR4 against Pseudomonas aeruginosa (ATCC 27853) coupled with production optimization. In order to achieve this, combinations of media components were formulated by employing statistical tools such as Plackett-Burman analysis and central composite rotatable design (CCRD). It was evident that at 35mlL-1 glycerol and 3.8gL-1 casamino acid, anti-biofilm activity and production of extracellular protein significantly increased by 1.5-fold and 1.2-fold, respectively. These results corroborate that the combination of glycerol and casamino acid plays a key role in the production of BIC. Further, metabolic profiling of BIC was carried out using liquid chromatography/tandem mass spectrometry (LC-MS/MS) based on m/z value. The presence of Stigmatellin Y was predicted with monoisotopic neutral mass of 484.2825Da. In support of optimization study, higher production of BIC was confirmed in the optimized-media-grown BR4 (OPT-BR4) than in the ideal-media-grown BR4 (ID-BR4) by LC-MS/MS analysis. PqsR in P. aeruginosa is a potential target for anti-virulent therapy. Molecular docking study has revealed that Stigmatellin Y interacts with PqsR in the similar orientation like a cognate signal (PQS) and synthetic inhibitor. In addition, Stigmatellin Y was found to exhibit interaction with four more amino acid residues of PqsR to establish strong affinity. Stigmatellin Y thus might play a role of competitor for PQS to distract PQS-PqsR mediated communication in P. aeruginosa. The present investigation thus paves new avenues to develop anti-Pseudomonas virulent therapy.

Target specific proteochemometric model development for BACE1 - protein flexibility and structural water are critical in virtual screening.

BACE1 is an attractive target in Alzheimer's disease (AD) treatment. A rational drug design effort for the inhibition of BACE1 is actively pursued by researchers in both academic and pharmaceutical industries. This continued effort led to the steady accumulation of BACE1 crystal structures, co-complexed with different classes of inhibitors. This wealth of information is used in this study to develop target specific proteochemometric models and these models are exploited for predicting the prospective BACE1 inhibitors. The models developed in this study have performed excellently in predicting the computationally generated poses, separately obtained from single and ensemble docking approaches. The simple protein-ligand contact (SPLC) model outperforms other sophisticated high end models, in virtual screening performance, developed during this study. In an attempt to account for BACE1 protein active site flexibility information in predictive models, we included the change in the area of solvent accessible surface and the change in the volume of solvent accessible surface in our models. The ensemble and single receptor docking results obtained from this study indicate that the structural water mediated interactions improve the virtual screening results. Also, these waters are essential for recapitulating bioactive conformation during docking study. The proteochemometric models developed in this study can be used for the prediction of BACE1 inhibitors, during the early stage of AD drug discovery.

Neutralization function affected by single amino acid replacement in the HIV-1 antibody targets.

The viral envelope glycoproteins are essential for entry into their host cells and studied extensively for designing vaccines. We hypothesize that the glycosylation on the HIV-1 viral envelope glycoprotein 41(gp41) at critical residues offers viral escape from the specific immune surveillant neutralizing antibodies Z13, 4E10 and 10E8 targeted to their linear epitopes in the Membrane Proximal External Region (MPER). The glycosylation occurring on the 50th residue (Asparagine) contained in the target (NWFNIT) can mask itself to be inaccessible for these neutralizing antibodies. The glycosylation rate of the epitopes which are shared by the Z13, 4E10 and 10E8 neutralizing antibodies of HIV-1 were predicited in silico. We analyzed the reliable frequency of glycosylation on the HIV-1 envelope gp41 using prediction tools to unravel the plausibility of the glycosylation by a mannose at 50th residue in the 59 amino acid long HIV-gp41 trimer (PDBID: 2M7W and 2LP7). It is evident that the glycosylation by a mannose that masks these targets is possible only when the 50th amino-acid is N (Asparagine, Asn) which is not possible when N is mutated to D (Aspartatic acid, Asp). The additive advantage for the retrovirus is its error-prone reverse transcriptase which can choose to copy these survivable mutants with Asn N-50 that can be glycosylated as explained by the Copy-choice model. So the glycan shields varying in their intensity and patterns have to be essentially studied to understand the viral escape strategies that will give a way forward towards a successful vaccine that can elicit a neutralizing antibody response to confer protection.

Modeling the closed and open state conformations of the GABA(A) ion channel--plausible structural insights for channel gating.

Recent disclosure of high resolution crystal structures of Gloeobacter violaceus (GLIC) in open state and Erwinia chrysanthemii (ELIC) in closed state provides newer avenues to advance our knowledge and understanding of the physiologically and pharmacologically important ionotropic GABA(A) ion channel. The present modeling study envisions understanding the complex molecular transitions involved in ionic conductance, which were not evident in earlier disclosed homology models. In particular, emphasis was put on understanding the structural basis of gating, gating transition from the closed to the open state on an atomic scale. Homology modeling of two different physiological states of GABA(A) was carried out using their respective templates. The ability of induced fit docking in breaking the critical inter residue salt bridge (Glu155ß(2) and Arg207ß(2)) upon endogenous GABA docking reflects the perceived side chain rearrangements that occur at the orthosteric site and consolidate the quality of the model. Biophysical calculations like electrostatic mapping, pore radius calculation, ion solvation profile, and normal-mode analysis (NMA) were undertaken to address pertinent questions like the following: How the change in state of the ion channel alters the electrostatic environment across the lumen; How accessible is the Cl(-) ion in the open state and closed state; What structural changes regulate channel gating. A "Twist to Turn" global motion evinced at the quaternary level accompanied by tilting and rotation of the M2 helices along the membrane normal rationalizes the structural transition involved in gating. This perceived global motion hints toward a conserved gating mechanism among pLGIC. To paraphrase, this modeling study proves to be a reliable framework for understanding the structure function relationship of the hitherto unresolved GABA(A) ion channel. The modeled structures presented herein not only reveal the structurally distinct conformational states of the GABA(A) ion channel but also explain the biophysical difference between the respective states.

Rationalizing lead optimization by consensus 2D- CoMFA CoMSIA GRIND (3D) QSAR guided fragment hopping in search of γ-secretase inhibitors.

γ-Secretase (Gamma Secretase) is a potential drug target in Alzheimer's disease therapeutics. A sequel lead design study was undertaken on a series of bicyclononanes with an aim of identifying potent isofunctional chemotypes. Fragment-based bioisosteric replacement, which considers shape, chemistry, and electrostatics was carried out to mine over four million medicinally relevant fragments of Brood database. The resulting subset, thus, obtained was further mined using consensus QSAR developed from 2D and CoMFA, CoMSIA, GRIND (3D) QSAR predicted endpoints with superior statistical results. The employed consensus prediction and the predicted endpoint values were found to be in good agreement with the experimental values. The predictive ability of the generated model was validated using different statistical metrics, and similarity-based coverage estimation was carried out to define applicability boundaries. Few analogs designed, using the concept of bioisosterism, were found to be promising and could be considered for synthesis and subsequent screening.

Discovery of safe and orally effective 4-aminoquinaldine analogues as apoptotic inducers with activity against experimental visceral leishmaniasis.

Novel antileishmanials are urgently required to overcome emergence of drug resistance, cytotoxic effects, and difficulties in oral delivery. Toward this, we investigated a series of novel 4-aminoquinaldine derivatives, a new class of molecules, as potential antileishmanials. 4-Aminoquinaldine derivatives presented inhibitory effects on L. donovani promastigotes and amastigotes (50% inhibitory concentration range, 0.94 to 127 µM). Of these, PP-9 and PP-10 were the most effective in vitro and demonstrated strong efficacies in vivo through the intraperitoneal route. They were also found to be effective against both sodium antimony gluconate-sensitive and -resistant Leishmania donovani strains in BALB/c mice when treated orally, resulting in more than 95% protection. Investigation of their mode of action revealed that killing by PP-10 involved moderate inhibition of dihydrofolate reductase and elicitation of the apoptotic cascade. Our studies implicate that PP-10 augments reactive oxygen species generation, evidenced from decreased glutathione levels and increased lipid peroxidation. Subsequent disruption of Leishmania promastigote mitochondrial membrane potential and activation of cytosolic proteases initiated the apoptotic pathway, resulting in DNA fragmentation and parasite death. Our results demonstrate that PP-9 and PP-10 are promising lead compounds with the potential for treating visceral leishmaniasis (VL) through the oral route.

Rationalizing fragment based drug discovery for BACE1: insights from FB-QSAR, FB-QSSR, multi objective (MO-QSPR) and MIF studies.

The ability to identify fragments that interact with a biological target is a key step in FBDD. To date, the concept of fragment based drug design (FBDD) is increasingly driven by bio-physical methods. To expand the boundaries of QSAR paradigm, and to rationalize FBDD using In silico approach, we propose a fragment based QSAR methodology referred here in as FB-QSAR. The FB-QSAR methodology was validated on a dataset consisting of 52 Hydroxy ethylamine (HEA) inhibitors, disclosed by GlaxoSmithKline Pharmaceuticals as potential anti-Alzheimer agents. To address the issue of target selectivity, a major confounding factor in the development of selective BACE1 inhibitors, FB-QSSR models were developed using the reported off target activity values. A heat map constructed, based on the activity and selectivity profile of the individual R-group fragments, and was in turn used to identify superior R-group fragments. Further, simultaneous optimization of multiple properties, an issue encountered in real-world drug discovery scenario, and often overlooked in QSAR approaches, was addressed using a Multi Objective (MO-QSPR) method that balances properties, based on the defined objectives. MO-QSPR was implemented using Derringer and Suich desirability algorithm to identify the optimal level of independent variables (X) that could confer a trade-off between selectivity and activity. The results obtained from FB-QSAR were further substantiated using MIF (Molecular Interaction Fields) studies. To exemplify the potentials of FB-QSAR and MO-QSPR in a pragmatic fashion, the insights gleaned from the MO-QSPR study was reverse engineered using Inverse-QSAR in a combinatorial fashion to enumerate some prospective novel, potent and selective BACE1 inhibitors.

Probing the structure of Leishmania donovani chagasi DHFR-TS: comparative protein modeling and protein-ligand interaction studies.

Dihydrofolate reductase (DHFR) has been used successfully as a drug target in the area of anti-bacterial, anti-cancer and anti-malarial therapy. It also acts as a drug target for Leishmaniasis. Inhibition of DHFR leads to cell death through lack of thymine (nucleotide metabolism). Although the crystal structures of Leishmania major and Trypanosoma cruzi DHFR-thymidylate synthase (TS) have been resolved, to date there is no three-dimensional (3D)-structural information on DHFR-TS of Leishmania donovani chagasi, which causes visceral leishmaniasis. Our aim in this study was to model the 3D structure of L. donovani chagasi DHFR-TS, and to investigate the structural requirements for its inhibition. In this paper we describe a highly refined homology model of L. donovani chagasi DHFR-TS based on available crystallographic structures by using the Homology module of Insight II. Structural refinement and minimization of the generated L. donovani chagasi DHFR-TS model employed the Discover 3 module of Insight II and molecular dynamic simulations. The model was further validated through use of the PROCHECK, Verify_3D, PROSA, PSQS and ERRAT programs, which confirm that the model is reliable. Superimposition of the model structure with the templates L. major A chain, L. major B chain And T. cruzi A chain showed root mean square deviations of 0.69 A, 0.71 A and 1.11 A, respectively. Docking analysis of the L. donovani chagasi DHFR-TS model with methotrexate enabled us to identify specific residues, viz. Val156, Val30, Lys95, Lys75 and Arg97, within the L. donovani chagasi DHFR-TS binding pocket, that play an important role in ligand or substrate binding. Docking studies clearly indicated that these five residues are important determinants for binding as they have strong hydrogen bonding interactions with the ligand.