Ligand Specificity of Group I Biotin Protein Ligase of Mycobacterium tuberculosis

Background: Fatty acids are indispensable constituents of mycolic acids that impart toughness & permeability barrier to the cell envelope of M. tuberculosis. Biotin is an essential co-factor for acetyl-CoA carboxylase (ACC) the enzyme involved in the synthesis of malonyl-CoA, a committed precursor, needed for fatty acid synthesis. Biotin carboxyl carrier protein (BCCP) provides the co-factor for catalytic activity of ACC. Methodology/Principal Findings: BPL/BirA (Biotin Protein Ligase), and its substrate, biotin carboxyl carrier protein (BCCP) of Mycobacterium tuberculosis (Mt) were cloned and expressed in E. coli BL21. In contrast to EcBirA and PhBPL, the similar to 29.5 kDa MtBPL exists as a monomer in native, biotin and bio-5'AMP liganded forms. This was confirmed by molecular weigt profiling by gel filtration on Superdex S-200 and Dynamic Light Scattering (DLS). Computational docking of biotin and bio-5'AMP to MtBPL show that adenylation alters the contact residues for biotin. MtBPL forms 11 H-bonds with biotin, relative to 35 with bio-5'AMP. Docking simulations also suggest that bio-5'AMP hydrogen bonds to the conserved `GRGRRG' sequence but not biotin. The enzyme catalyzed transfer of biotin to BCCP was confirmed by incorporation of radioactive biotin and by Avidin blot. The K-m for BCCP was similar to 5.2 mu M and similar to 420 nM for biotin. MtBPL has low affinity (K-b = 1.06 x 10(-6) M) for biotin relative to EcBirA but their K-m are almost comparable suggesting that while the major function of MtBPL is biotinylation of BCCP, tight binding of biotin/bio-5'AMP by EcBirA is channeled for its repressor activity. Conclusions/Significance: These studies thus open up avenues for understanding the unique features of MtBPL and the role it plays in biotin utilization in M. tuberculosis.

Parallel implementation of AutoDock

Computational docking of ligands to protein structures is a key step in structure-based drug design. Currently, the time required for each docking run is high and thus limits the use of docking in a high-throughput manner, warranting parallelization of docking algorithms. AutoDock, a widely used tool, has been chosen for parallelization. Near-linear increases in speed were observed with 96 processors, reducing the time required for docking ligands to HIV-protease from 81 min, as an example, on a single IBM Power-5 processor ( 1.65 GHz), to about 1 min on an IBM cluster, with 96 such processors. This implementation would make it feasible to perform virtual ligand screening using AutoDock.

Predicting the Pose of β-Casomorphin-5 and 7 in the Opioid Receptors

The opioid receptors consist of three main subtypes; μ, δ, and κ. Previous binding studies have shown that fragments of the milk protein, β-casein, known as β-casomorphins are agonists of these receptors which are selective for the μ receptor subtype. Using the crystal structures of these three receptors, computational molecular docking studies were done using the software GOLD to determine the conformation of β-casomorphin-5 and 7 when they bind to these three opioid receptors. GOLD was able to discriminate among the three receptors when docking the rigid ligands co-crystalized with the receptors. However, GOLD could not discriminate among the three receptors for either of the highly flexible β-casomorphins. A per amino acid scoring method was developed to overcome this problem. This method was used to predict the conformation of both β-casomorphin-5 and 7 in the μ receptor and determine that the two amino acid residues, Lys303 and Trp318 of the μ receptor are responsible for discriminating among the three receptor subtypes for binding of the β-casomorphin-5 and 7.

MPO Inhibitors Selected by Virtual Screening

The hemeprotein myeloperoxidase (MPO) participates in innate immune defense through its ability to generate potent microbicidal oxidants. However, these oxidants are also key mediators of the tissue damage associated with many inflammatory diseases. Thus, there is considerable interest in developing therapeutically useful MPO inhibitors. Here, we used structure-based drug design (SBDD) and ligand-based drug design (LBDD) to select for potentially new and selective MPO inhibitors. A pharmacophore model was developed based on the crystal structure of human MPO in complex with salicylhydroxamic acid (SHA), a known inhibitor of the enzyme. The pharmacophore model was used to screen the ZINC database for potential ligands, which were further filtered on the basis of their physical-chemical properties and docking score. The filtered compounds were visually inspected, and nine were purchased for experimental studies. Surprisingly, almost all of the selected compounds belonged to the aromatic hydrazide class, which had been previously described as MPO inhibitors. The compounds selected by virtual screening were shown to inhibit the chlorinating activity of MPO; the top four compounds displayed IC(50) values ranging from 1.0 to 2.8 mM. MPO inactivation by the most effective compound was shown to be irreversible. Overall, our results show that SBDD and LBDD may be useful for the rational development of new MPO inhibitors.

FERM domain interaction with myosin negatively regulates FAK in cardiomyocyte hypertrophy

Focal adhesion kinase (FAK) regulates cellular processes that affect several aspects of development and disease. The FAK N-terminal FERM (4.1 protein-ezrin-radixin-moesin homology) domain, a compact clover-leaf structure, binds partner proteins and mediates intramolecular regulatory interactions. Combined chemical cross-linking coupled to MS, small-angle X-ray scattering, computational docking and mutational analyses showed that the FAK FERM domain has a molecular cleft (similar to 998 angstrom(2)) that interacts with sarcomeric myosin, resulting in FAK inhibition. Accordingly, mutations in a unique short amino acid sequence of the FERM myosin cleft, FP-1, impaired the interaction with myosin and enhanced FAK activity in cardiomyocytes. An FP-1 decoy peptide selectively inhibited myosin interaction and increased FAK activity, promoting cardiomyocyte hypertrophy through activation of the AKT-mammalian target of rapamycin pathway. Our findings uncover an inhibitory interaction between the FAK FERM domain and sarcomeric myosin that presents potential opportunities to modulate the cardiac hypertrophic response through changes in FAK activity.

Optimization of anti-Trypanosoma cruzi oxadiazoles leads to identification of compounds with efficacy in infected mice

We recently showed that oxadiazoles have anti-Trypanosoma cruzi activity at micromolar concentrations. These compounds are easy to synthesize and show a number of clear and interpretable structure-activity relationships (SAR), features that make them attractive to pursue potency enhancement. We present here the structural design, synthesis, and anti-T. cruzi evaluation of new oxadiazoles denoted 5a-h and 6a-h. The design of these compounds was based on a previous model of computational docking of oxadiazoles on the T. cruzi protease cruzain. We tested the ability of these compounds to inhibit catalytic activity of cruzain, but we found no correlation between the enzyme inhibition and the antiparasitic activity of the compounds. However, we found reliable SAR data when we tested these compounds against the whole parasite. While none of these oxadiazoles showed toxicity for mammalian cells, oxadiazoles 6c (fluorine), 6d (chlorine), and 6e (bromine) reduced epimastigote proliferation and were cidal for trypomastigotes of T. cruzi Y strain. Oxadiazoles 6c and 6d have IC50 of 9.5 +/- 2.8 and 3.5 +/- 1.8 mu M for trypomastigotes, while Benznidazole, which is the currently used drug for Chagas disease treatment, showed an IC50 of 11.3 +/- 2.8 mu M. Compounds 6c and 6d impair trypomastigote development and invasion in macrophages, and also induce ultrastructural alterations in trypomastigotes. Finally, compound 6d given orally at 50 mg/kg substantially reduces the parasitemia in T. cruzi-infected BALB/c mice. Our drug design resulted in potency enhancement of oxadiazoles as anti-Chagas disease agents, and culminated with the identification of oxadiazole 6d, a trypanosomicidal compound in an animal model of infection. (C) 2012 Elsevier Ltd. All rights reserved.

Epigallocatechin gallate (EGCG) affects glucose metabolism and enhances fitness and life span in Drosophila melanogaster

In this study, we tested whether a standardized epigallocatechin-3-gallate (EGCG) rich green tea extract (comprising > 90% EGCG) affects fitness and lifespan as well as parameters of glucose metabolism and energy homeostasis in the fruit fly, Drosophila melanogaster. Following the application of the green tea extract a significant increase in the mean lifespan (+ 3.3 days) and the 50% survival (+ 4.3 days) as well as improved fitness was detected. These effects went along an increased expression of Spargel, the homolog of mammalian PGC1α, which has been reported to affect lifespan in flies. Intriguingly, in flies, treatment with the green tea extract decreased glucose concentrations, which were accompanied by an inhibition of α-amylase and α-glucosidase activity. Computational docking analysis proved the potential of EGCG to dock into the substrate binding pocket of α-amylase and to a greater extent into α-glucosidase. Furthermore, we demonstrate that EGCG downregulates insulin-like peptide 5 and phosphoenolpyruvate carboxykinase, major regulators of glucose metabolism, as well as the Drosophila homolog of leptin, unpaired 2. We propose that a decrease in glucose metabolism in connection with an upregulated expression of Spargel contribute to the better fitness and the extended lifespan in EGCG-treated flies.

New cytochromes P450: Tale of two variants

Cytochrome P450s, a superfamily of heme enzymes found in most living organisms. They are responsible for metabolism of many therapeutic drugs, industrial pollutants, carcinogens, and additives to foodstuffs, as well as some endogenous compounds including fatty acids and steroids. First pass drug metabolism studies represent mainly liver and small intestine elimination, and are viewed as the standard to predict therapeutic outcome. However, drug plasma levels determined after administration do not always correlate with therapeutic efficacy of the drug. Therefore, a possible explanation may come by understanding drug metabolism in extrahepatic tissues and/or at the site of drug action. Identification and characterization of novel tissue specific isoforms of P450 generated by alternative splicing of known P450 genes or as yet unidentified genes is essential to predict pharmacological outcome of drugs or the fate of a carcinogen that act at sites remote from liver. ^ Using RT-PCR, brain-specific cytochrome P450s were detected in samples of human autopsy brain. So far, we have identified two human brain variants including P450 2D7 and P450 1A1. We have shown the presence of the P450 1A1 brain specific splice variant in African Americans, Caucasians and Indians albeit different patterns of liver to brain variant ratio were seen distributed throughout each population. Interestingly, the splice variant was detected only in the brain but not in any other tissues from the same individual. Homology modeling was used to compare the variant 3D structure to the liver form structure and differences in the substrate access channels and substrate binding sites were noticed. Automated computational docking was used to predict the metabolic fate of the potent carcinogenic substrate, benzo[a]pyrene. P450 1A1 brain variant showed no binding orientations that could produce the active metabolite, whereas P450 1A1 liver form did reveal orientations capable of generating active carcinogenic product. In vitro P32 labeling studies verified the docking predictions. Therefore, the data support the hypothesis that P450 brain splice variants mediate the metabolism of xenobiotics by mechanisms distinct from the well-studied liver counterparts. ^

Structure-activity relationships for a new family of sulfonylurea herbicides

Acetohydroxyacid synthase (AHAS; EC 2.2.1.6) catalyzes the first common step in branched-chain amino acid biosynthesis. The enzyme is inhibited by several chemical classes of compounds and this inhibition is the basis of action of the sulfonylurea and imidazolinone herbicides. The commercial sulfonylureas contain a pyrimidine or a triazine ring that is substituted at both meta positions, thus obeying the initial rules proposed by Levitt. Here we assess the activity of 69 monosubstituted sulfonylurea analogs and related compounds as inhibitors of pure recombinant Arabidopsis thaliana AHAS and show that disubstitution is not absolutely essential as exemplified by our novel herbicide, monosulfuron (2-nitro-N-(4'-methyl-pyrimidin-2'-yl) phenyl-sulfonylurea), which has a pyrimidine ring with a single meta substituent. A subset of these compounds was tested for herbicidal activity and it was shown that their effect in vivo correlates well with their potency in vitro as AHAS inhibitors. Three-dimensional quantitative structure-activity relationships were developed using comparative molecular field analysis and comparative molecular similarity indices analysis. For the latter, the best result was obtained when steric, electrostatic, hydrophobic and H-bond acceptor factors were taken into consideration. The resulting fields were mapped on to the published crystal structure of the yeast enzyme and it was shown that the steric and hydrophobic fields are in good agreement with sulfonylurea-AHAS interaction geometry.

Unexpected Activity of a Novel Kunitz-Type Inhibtior: Inhibition of Cysteine Proteases but not Serine Proteases

Kunitz-type (KT) protease inhibitors are low molecular weight proteins classically defined as serine protease inhibitors. We identified a novel secreted KT inhibitor associated with the gut and parenchymal tissues of the infective juvenile stage of Fasciola hepatica, a helminth parasite of medical and veterinary importance. Unexpectedly, recombinant KT inhibitor (rFhKT1) exhibited no inhibitory activity towards serine proteases but was a potent inhibitor of the major secreted cathepsin L cysteine proteases of F. hepatica, FhCL1 and FhCL2, and of human cathepsins L and K (Ki = 0.4 nM - 27 nM). FhKT1 prevented the auto-catalytic activation of FhCL1 and FhCL2 and formed stable complexes with the mature enzymes. Pull-down experiments from adult parasite culture medium showed that rFhKT1 interacts specifically with native secreted FhCL1, FhCL2 and FhCL5. Substitution of the unusual P1 Leu15 within the exposed reactive loop of FhKT1 for the more commonly found Arg (FhKT1Leu15/Arg15) had modest adverse effects on the cysteine protease inhibition but conferred potent activity against the serine protease trypsin (Ki = 1.5 nM). Computational docking and sequence analysis provided hypotheses for the exclusive binding of FhKT1 to cysteine proteases, the importance of the Leu15 in anchoring the inhibitor into the S2 active site pocket, and the inhibitor's selectivity towards FhCL1, FhCL2 and human cathepsins L and K. FhKT1 represents a novel evolutionary adaptation of KT protease inhibitors by F. hepatica, with its prime purpose likely in the regulation of the major parasite-secreted proteases and/or cathepsin L-like proteases of its host.

Computational strategies to include protein flexibility in Ligand Docking and Virtual Screening

The dynamic character of proteins strongly influences biomolecular recognition mechanisms. With the development of the main models of ligand recognition (lock-and-key, induced fit, conformational selection theories), the role of protein plasticity has become increasingly relevant. In particular, major structural changes concerning large deviations of protein backbones, and slight movements such as side chain rotations are now carefully considered in drug discovery and development. It is of great interest to identify multiple protein conformations as preliminary step in a screening campaign. Protein flexibility has been widely investigated, in terms of both local and global motions, in two diverse biological systems. On one side, Replica Exchange Molecular Dynamics has been exploited as enhanced sampling method to collect multiple conformations of Lactate Dehydrogenase A (LDHA), an emerging anticancer target. The aim of this project was the development of an Ensemble-based Virtual Screening protocol, in order to find novel potent inhibitors. On the other side, a preliminary study concerning the local flexibility of Opioid Receptors has been carried out through ALiBERO approach, an iterative method based on Elastic Network-Normal Mode Analysis and Monte Carlo sampling. Comparison of the Virtual Screening performances by using single or multiple conformations confirmed that the inclusion of protein flexibility in screening protocols has a positive effect on the probability to early recognize novel or known active compounds.

Computational analyses of the catalytic and heparin-binding sites and their interactions with glycosaminoglycans in glycoside hydrolase family 79 endo-d-glucuronidase (heparanase)

Mammalian heparanase is an endo-β-glucuronidase associated with cell invasion in cancer metastasis, angiogenesis and inflammation. Heparanase cleaves heparan sulfate proteoglycans in the extracellular matrix and basement membrane, releasing heparin/heparan sulfate oligosaccharides of appreciable size. This in turn causes the release of growth factors, which accelerate tumor growth and metastasis. Heparanase has two glycosaminoglycan-binding domains; however, no three-dimensional structure information is available for human heparanase that can provide insights into how the two domains interact to degrade heparin fragments. We have constructed a new homology model of heparanase that takes into account the most recent structural and bioinformatics data available. Heparin analogs and glycosaminoglycan mimetics were computationally docked into the active site with energetically stable ring conformations and their interaction energies were compared. The resulting docked structures were used to propose a model for substrates and conformer selectivity based on the dimensions of the active site. The docking of substrates and inhibitors indicates the existence of a large binding site extending at least two saccharide units beyond the cleavage site (toward the nonreducing end) and at least three saccharides toward the reducing end (toward heparin-binding site 2). The docking of substrates suggests that heparanase recognizes the N-sulfated and O-sulfated glucosamines at subsite +1 and glucuronic acid at the cleavage site, whereas in the absence of 6-O-sulfation in glucosamine, glucuronic acid is docked at subsite +2. These findings will help us to focus on the rational design of heparanase-inhibiting molecules for anticancer drug development by targeting the two heparin/heparan sulfate recognition domains.

Design, Synthesis, and DNA Binding Properties of Photoisomerizable Azobenzene−Distamycin Conjugates: An Experimental and Computational Study

Here, we present the synthesis, photochemical, and DNA binding properties of three photoisomerizable azobenzene−distamycin conjugates in which two distamycin units were linked via electron-rich alkoxy or electron-withdrawing carboxamido moieties with the azobenzene core. Like parent distamycin A, these molecules also demonstrated AT-specific DNA binding. Duplex DNA binding abilities of these conjugates were found to depend upon the nature and length of the spacer, the location of protonatable residues, and the isomeric state of the conjugate. The changes in the duplex DNA binding efficiency of the individual conjugates in the dark and with their respective photoirradiated forms were examined by circular dichroism, thermal denaturation of DNA, and Hoechst displacement assay with poly[d(A-T).d(T-A)] DNA in 150 mM NaCl buffer. Computational structural analyses of the uncomplexed ligands using ab initio HF and MP2 theory and molecular docking studies involving the conjugates with duplex d[(GC(AT)10CG)]2 DNA were performed to rationalize the nature of binding of these conjugates.

Using Correlated Parameters for Improved Ranking of Protein-Protein Docking Decoys

A successful protein-protein docking study culminates in identification of decoys at top ranks with near-native quaternary structures. However, this task remains enigmatic because no generalized scoring functions exist that effectively infer decoys according to the similarity to near-native quaternary structures. Difficulties arise because of the highly irregular nature of the protein surface and the significant variation of the nonbonding and solvation energies based on the chemical composition of the protein-protein interface. In this work, we describe a novel method combining an interface-size filter, a regression model for geometric compatibility (based on two correlated surface and packing parameters), and normalized interaction energy (calculated from correlated nonbonded and solvation energies), to effectively rank decoys from a set of 10,000 decoys. Tests on 30 unbound binary protein-protein complexes show that in 16 cases we can identify at least one decoy in top three ranks having <= 10 angstrom backbone root mean square deviation from true binding geometry. Comparisons with other state-of-art methods confirm the improved ranking power of our method without the use of any experiment-guided restraints, evolutionary information, statistical propensities, or modified interaction energy equations. Tests on 118 less-difficult bound binary protein-protein complexes with <= 35% sequence redundancy at the interface showed that in 77% cases, at least 1 in 10,000 decoys were identified with <= 5 angstrom backbone root mean square deviation from true geometry at first rank. The work will promote the use of new concepts where correlations among parameters provide more robust scoring models. It will facilitate studies involving molecular interactions, including modeling of large macromolecular assemblies and protein structure prediction. (C) 2010 Wiley Periodicals, Inc. J Comput Chem 32: 787-796, 2011.

PRUNE and PROBE-two modular web services for protein-protein docking

The protein-protein docking programs typically perform four major tasks: (i) generation of docking poses, (ii) selecting a subset of poses, (iii) their structural refinement and (iv) scoring, ranking for the final assessment of the true quaternary structure. Although the tasks can be integrated or performed in a serial order, they are by nature modular, allowing an opportunity to substitute one algorithm with another. We have implemented two modular web services, (i) PRUNE: to select a subset of docking poses generated during sampling search (http://pallab.serc.iisc.ernet.in/prune) and (ii) PROBE: to refine, score and rank them (http://pallab.serc.iisc.ernet.in/probe). The former uses a new interface area based edge-scoring function to eliminate > 95% of the poses generated during docking search. In contrast to other multi-parameter-based screening functions, this single parameter based elimination reduces the computational time significantly, in addition to increasing the chances of selecting native-like models in the top rank list. The PROBE server performs ranking of pruned poses, after structure refinement and scoring using a regression model for geometric compatibility, and normalized interaction energy. While web-service similar to PROBE is infrequent, no web-service akin to PRUNE has been described before. Both the servers are publicly accessible and free for use.