DEPTH: a web server to compute depth and predict small-molecule binding cavities in proteins

Depth measures the extent of atom/residue burial within a protein. It correlates with properties such as protein stability, hydrogen exchange rate, protein-protein interaction hot spots, post-translational modification sites and sequence variability. Our server, DEPTH, accurately computes depth and solvent-accessible surface area (SASA) values. We show that depth can be used to predict small molecule ligand binding cavities in proteins. Often, some of the residues lining a ligand binding cavity are both deep and solvent exposed. Using the depth-SASA pair values for a residue, its likelihood to form part of a small molecule binding cavity is estimated. The parameters of the method were calibrated over a training set of 900 high-resolution X-ray crystal structures of single-domain proteins bound to small molecules (molecular weight < 1.5 KDa). The prediction accuracy of DEPTH is comparable to that of other geometry-based prediction methods including LIGSITE, SURFNET and Pocket-Finder (all with Matthew's correlation coefficient of similar to 0.4) over a testing set of 225 single and multi-chain protein structures. Users have the option of tuning several parameters to detect cavities of different sizes, for example, geometrically flat binding sites. The input to the server is a protein 3D structure in PDB format. The users have the option of tuning the values of four parameters associated with the computation of residue depth and the prediction of binding cavities. The computed depths, SASA and binding cavity predictions are displayed in 2D plots and mapped onto 3D representations of the protein structure using Jmol. Links are provided to download the outputs. Our server is useful for all structural analysis based on residue depth and SASA, such as guiding site-directed mutagenesis experiments and small molecule docking exercises, in the context of protein functional annotation and drug discovery.

Remarkable photocytotoxicity in hypoxic HeLa cells by a dipyridophenazine copper(II) Schiff base thiolate

Copper(II) complexes Cu(satp)(L)] (1-3) of a Schiff base thiolate (salicylidene-2-aminothiophenol, H(2)satP) and phenanthroline bases (L), viz. 1,10-phenanthroline (phen in 1), dipyrido3,2-d:2',3'-f]quinoxaline (dpq in 2) and dipyrido3,2-a:2',3'-c]phenazine (dppz in 3), were prepared, characterized and their anaerobic DNA photocleavage activity and hypoxic photocytotoxicity studied. The redox active complexes show the Cu(II)-Cu(I) couple near -0.5 V for 1 and near 0.0 V vs. SCE (saturated calomel electrode) for 2 and 3. The one-electron paramagnetic complexes (similar to 1.85 mu(B)) are avid DNA binders giving K(b) values within 1.0 x 10(5) - 8.0 x 10(5) M(-1). Thermal melting and viscosity data along with molecular docking calculations suggest DNA groove and/or partial intercalative binding of the complexes. The complexes show anaerobic DNA cleavage activity in red light under argon via type-I pathway, while DNA photocleavage in air proceeds via hydroxyl radical pathway. The DFT (density functional theory) calculations reveal a thyil radical pathway for the anaerobic DNA photocleavage activity and suggest the possibility of generation of a transient copper(I) species due to bond breakage between the copper and sulfur to generate the thyil radical. An oxidation of the copper(I) species is likely by oxygen in an aerobic medium or by the buffer medium in an anaerobic condition. Complex 3 exhibits significant photocytotoxicity in HeLa cells (IC(50) = 8.3(+/- 1.0) mu M) in visible light, while showing lower dark toxicity (IC(50) = 17.2(+/- 1.0) mu M). A significant reduction in the dark toxicity is observed under hypoxic cellular conditions (IC(50) = 30.0(+/- 1.0) mu M in dark), while retaining its photocytotoxicity (IC(50) = 8.0(+/- 1.0) mu M). (C) 2011 Elsevier Inc. All rights reserved.

The Structure of the Thioredoxin-Triosephosphate Isomerase Complex Provides Insights into the Reversible Glutathione-Mediated Regulation of Triosephosphate Isomerase

Protein-protein interactions are crucial for many biological functions. The redox interactome encompasses numerous weak transient interactions in which thioredoxin plays a central role. Proteomic studies have shown that thioredoxin binds to numerous proteins belonging to various cellular processes, including energy metabolism. Thioredoxin has cross talk with other redox mechanisms involving glutathionylation and has functional overlap with glutaredoxin in deglutathionylation reactions. In this study, we have explored the structural and biochemical interactions of thioredoxin with the glycolytic enzyme, triosephosphate isomerase. Nuclear magnetic resonance chemical shift mapping methods and molecular dynamics-based docking have been applied in deriving a structural model of the thioredoxin-triosephosphate isomerase complex. The spatial proximity of active site cysteine residues of thioredoxin to reactive thiol groups on triosephosphate isomerase provides a direct link to the observed deglutathionylation of cysteine 217 in triosephosphate isomerase, thereby reversing the inhibitory effect of S-glutathionylation of triosephosphate isomerase.

The Structure of the Thioredoxin−Triosephosphate Isomerase Complex Provides Insights into the Reversible Glutathione-Mediate Regulation of Triosephosphate Isomerase

Protein−protein interactions are crucial for many biological functions. The redox interactome encompasses numerous weak transient interactions in which thioredoxin plays a central role. Proteomic studies have shown that thioredoxin binds to numerous proteins belonging to various cellular processes, including energy metabolism. Thioredoxin has cross talk with other redox mechanisms involving glutathionylation and has functional overlap with glutaredoxin in deglutathionylation reactions. In this study, we have explored the structural and biochemical interactions of thioredoxin with the glycolytic enzyme, triosephosphate isomerase. Nuclear magnetic resonance chemical shift mapping methods and molecular dynamics-based docking have been applied in deriving a structural model of the thioredoxin−triosephosphate isomerase complex. The spatial proximity of active site cysteine residues of thioredoxin to reactive thiol groups on triosephosphate isomerase provides a direct link to the observed deglutathionylation of cysteine 217 in triosephosphate isomerase, thereby reversing the inhibitory effect of S-glutathionylation of triosephosphate isomerase.

PAMAM Dendrimer-Drug Interactions: Effect of pH on the Binding and Release Pattern

Understanding the dendrimer-drug interaction is of great importance to design and optimize the dendrimer-based drug delivery system. Using atomistic molecular dynamics (MD) simulations, we have analyzed the release pattern of four ligands (two soluble drugs, namely, salicylic acid (Sal), L-alanine (Ala), and two insoluble drugs, namely, phenylbutazone (Pbz) and primidone (Prim)), which were initially encapsulated inside the ethylenediamine (EDA) cored polyamidoamine (PAMAM) dendrimer using the docking method. We have computed the potential of mean force (PMF) variation with generation 5 (G5)-PAMAM dendrimer complexed with drug molecules using umbrella sampling. From our calculated PMF values, we observe that soluble drugs (Sal and Ala) have lower energy barriers than insoluble drugs (Pbz and Prim). The order of ease of release pattern for these drugs from G5 protonated PAMAM dendrimer was found to be Ala > Sal > Prim > Pbz. In the case of insoluble drugs (Prim and Pbz), because of larger size, we observe much nonpolar contribution, and thus, their larger energy barriers can be reasoned to van der Waals contribution. From the hydrogen bonding analysis of the four PAMAM drug complexes under study, we found intermolecular hydrogen bonding to show less significant contribution to the free energy barrier. Another interesting feature appears while calculating the PMF profile of G5NP (nonprotonated)-PAMAM Pbz and G5NP (nonprotonated)-PAMAM-Sal complex. The PMF was found to be less when the drug is bound to nonprotonated dendrimer compared to the protonated dendrimer. Our results suggest that encapsulation of the drug molecule into the host PAMAM dendrimer should be carried out at higher pH values (near pH 10). When such complex enters the human body, the pH is around 7.4 and at that physiological pH, the dendrimer holds the drug tightly. Hence the release of drug can occur at a controlled rate into the bloodstream. Thus, our findings provide a microscopic picture of the encapsulation and controlled release of drugs in the case of dendrimer-based host-guest systems.

Regioselective synthesis and antimicrobial studies of ester linked 1,4-disubstituted 1,2,3-bistriazoles

A series of 1,4-disubstituted 1,2,3-bistriazoles was synthesized via click chemistry by cycloaddition of various bisalkynes with benzyl/2-phenylethyl azide. Synthesized triazoles were characterized by IR, H-1 NMR, C-13 NMR and mass spectral techniques. All the compounds were evaluated for antibacterial/antifungal activities and found to possess moderate to good antimicrobial activities. Further the docking study for the most active compound against DNA Gyrase was also carried out. (C) 2012 Elsevier Ltd. All rights reserved.

A variant peptide of buffalo colostrum beta-lactoglobulin inhibits angiotensin I-converting enzyme activity

beta-lactoglobulin is a rich source of bioactive peptides. The LC-MS separated tryptic peptides of buffalo colostrum beta-lactoglobulin (BLG-col) were computed based on MS-MS fragmentation for de novo sequencing. Among the selected peptides (P1-P8), a variant was detected with methionine at position 74 instead of glutamate. The sequences of two peptides were identical to hypocholesterolemic peptides whereas the remaining peptides were in accordance with buffalo milk beta-lactoglobulin. Comparative sequence analysis of BLG-col to milk beta-lactoglobulin was carried out using CLUSTALW2 and a molecular model for BLG-col was constructed (PMDB ID-PM0076812). The synthesized variant pentapeptide (IIAMK, m/z-576 Da) was found to inhibit angiotensin I-converting enzyme (ACE) with an IC50 of 498 +/- 2 mu M, which was rationalized through docking simulations using Molgrow virtual docker. (C) 2012 Elsevier Masson SAS. All rights reserved.

A Rhodanine Derivative CCR-11 Inhibits Bacterial Proliferation by Inhibiting the Assembly and GTPase Activity of FtsZ

A perturbation of FtsZ assembly dynamics has been shown to inhibit bacterial cytokinesis. In this study, the antibacterial activity of 151 rhodanine compounds was assayed using Bacillus subtilis cells. Of 151 compounds, eight strongly inhibited bacterial proliferation at 2 mu M. Subsequently, we used the elongation of B. subtilis cells as a secondary screen to identify potential FtsZ-targeted antibacterial agents. We found that three compounds significantly increased bacterial cell length. One of the three compounds, namely, CCR-11 (E)-2-thioxo-5-({3-(trifluoromethyl)phenyl]furan-2-yl}methylene) thiazolidin-4-one], inhibited the assembly and GTPase activity of FtsZ in vitro. CCR-11 bound to FtsZ with a dissociation constant of 1.5 +/- 0.3 mu M. A docking analysis indicated that CCR-11 may bind to FtsZ in a cavity adjacent to the T7 loop and that short halogen oxygen, H-bonding, and hydrophobic interactions might be important for the binding of CCR-11 with FtsZ. CCR-11 inhibited the proliferation of B. subtilis cells with a half-maximal inhibitory concentration (IC50) of 1.2 +/- 0.2 mu M and a minimal inhibitory concentration of 3 mu M. It also potently inhibited proliferation of Mycobacterium smegmatis cells. Further, CCR-11 perturbed Z-ring formation in B. subtilis cells; however, it neither visibly affected nucleoid segregation nor altered the membrane integrity of the cells. CCR-11 inhibited HeLa cell proliferation with an IC50 value of 18.1 +/- 0.2,mu M (similar to 15 x IC50 of B. subtilis cell proliferation). The results suggested that CCR-11 inhibits bacterial cytokinesis by inhibiting FtsZ assembly, and it can be used as a lead molecule to develop FtsZ-targeted antibacterial agents.

Differential Activities of the Two Closely Related Withanolides, Withaferin A and Withanone: Bioinformatics and Experimental Evidences

Background and Purpose: Withanolides are naturally occurring chemical compounds. They are secondary metabolites produced via oxidation of steroids and structurally consist of a steroid-backbone bound to a lactone or its derivatives. They are known to protect plants against herbivores and have medicinal value including anti-inflammation, anti-cancer, adaptogenic and anti-oxidant effects. Withaferin A (Wi-A) and Withanone (Wi-N) are two structurally similar withanolides isolated from Withania somnifera, also known as Ashwagandha in Indian Ayurvedic medicine. Ashwagandha alcoholic leaf extract (i-Extract), rich in Wi-N, was shown to kill cancer cells selectively. Furthermore, the two closely related purified phytochemicals, Wi-A and Wi-N, showed differential activity in normal and cancer human cells in vitro and in vivo. We had earlier identified several genes involved in cytotoxicity of i-Extract in human cancer cells by loss-of-function assays using either siRNA or randomized ribozyme library. Methodology/Principal Findings: In the present study, we have employed bioinformatics tools on four genes, i.e., mortalin, p53, p21 and Nrf2, identified by loss-of-function screenings. We examined the docking efficacy of Wi-N and Wi-A to each of the four targets and found that the two closely related phytochemicals have differential binding properties to the selected cellular targets that can potentially instigate differential molecular effects. We validated these findings by undertaking parallel experiments on specific gene responses to either Wi-N or Wi-A in human normal and cancer cells. We demonstrate that Wi-A that binds strongly to the selected targets acts as a strong cytotoxic agent both for normal and cancer cells. Wi-N, on the other hand, has a weak binding to the targets; it showed milder cytotoxicity towards cancer cells and was safe for normal cells. The present molecular docking analyses and experimental evidence revealed important insights to the use of Wi-A and Wi-N for cancer treatment and development of new anti-cancer phytochemical cocktails.

A Bifunctional Enzyme That Has Both Monoacylglycerol Acyltransferase and Acyl Hydrolase Activities

Monoacylglycerol acyltransferase (MGAT) catalyzes the synthesis of diacylglycerol, the precursor of triacylglycerol biosynthesis and an important signaling molecule. Here, we describe the isolation and characterization of the peanut (Arachis hypogaea) MGAT gene. The soluble enzyme utilizes invariant histidine-62 and aspartate-67 residues of the acyltransferase motif for its MGAT activity. A sequence analysis revealed the presence of a hydrolase (GXSXG) motif, and enzyme assays revealed the presence of monoacylglycerol (MAG) and lysophosphatidylcholine (LPC) hydrolytic activities, indicating the bifunctional nature of the enzyme. The overexpression of the MGAT gene in yeast (Saccharomyces cerevisiae) caused an increase in triacylglycerol accumulation. Similar to the peanut MGAT, the Arabidopsis (Arabidopsis thaliana) homolog (At1g52760) also exhibited both acyltransferase and hydrolase activities. Interestingly, the yeast homolog lacks the conserved HX4D motif, and it is deficient in the acyltransferase function but exhibits MAG and LPC hydrolase activities. This study demonstrates the presence of a soluble MGAT/hydrolase in plants. The predicted three-dimensional homology modeling and substrate docking suggested the presence of two separate substrate (MAG and LPC)-binding sites in a single polypeptide. Our study describes a soluble bifunctional enzyme that has both MGAT and hydrolase functions.

Crystal structures and binding studies of atovaquone and its derivatives with cytochrome bc(1): a molecular basis for drug design

Crystal structure of trans-atovaquone (antimalarial drug), its polymorph and its stereoisomer (cis) along with five other derivatives with different functional groups have been analyzed. Based on the conformational features of these compounds and the characteristics of the nature of intermolecular interactions, valuable insights into the atomistic details of protein-inhibitor interactions have been derived by docking studies. Atovaquone and its derivatives pack in the crystal lattice using intermolecular O-H center dot center dot center dot O hydrogen bond dimer motifs supported by surrogate weak interactions including C-H center dot center dot center dot O and C-H center dot center dot center dot Cl hydrogen bonds. The docking results of these molecules with cytochrome bc(1) show preferences to form N-H center dot center dot center dot O, O-H center dot center dot center dot O and O-H center dot center dot center dot Cl hydrogen bonds. The involvement of halogen atoms in the binding pocket appears to be significant and is contrary to the theoretically predicted mechanism of protein-ligand docking reported earlier based on mimicking experimental binding results of stigmatellin with cytochrome bc(1). The significance of subtle energy factors controlled by weak intermolecular interactions appears to play a major role in drug binding.

CXI-benzo-84 reversibly binds to tubulin at colchicine site and induces apoptosis in cancer cells

Here, we have discovered CXI-benzo-84 as a potential anticancer agent from a library of benzimidazole derivatives using cell based screening strategy. CXI-benzo-84 inhibited cell cycle progression in metaphase stage of mitosis and accumulated spindle assembly checkpoint proteins Mad2 and BubR1 on kinetochores, which subsequently activated apoptotic cell death in cancer cells. CXI-benzo-84 depolymerized both interphase and mitotic microtubules, perturbed EB1 binding to microtubules and inhibited the assembly and GTPase activity of tubulin in vitro. CXI-benzo-84 bound to tubulin at a single binding site with a dissociation constant of 1.2 +/- 0.2 mu M. Competition experiments and molecular docking suggested that CXI-benzo-84 binds to tubulin at the colchicine-site. Further, computational analysis provided a significant insight on the binding site of CXI-benzo-84 on tubulin. In addition to its potential use in cancer chemotherapy, CXI-benzo-84 may also be useful to screen colchicine-site agents and to understand the colchicine binding site on tubulin. (C) 2013 Elsevier Inc. All rights reserved.

In vitro and in silico biological studies of novel thiazolo3,2-a]pyrimidine-6-carboxylate derivatives

In the present study, we report the synthesis, characterization of new series of thiazolo3,2-a]pyrimidine-6-carboxylate derivatives 3a-f and 4a-f. The newly synthesized compounds were screened for in vitro antimicrobial and antiviral activities. The probable mode of action of these active compounds was determined through in silico docking study by docking the receptor methionyl-tRNA synthetase and human inosine-5'-monophosphate dehydrogenase (IMPDH) for antibacterial and antiviral activities, respectively. Among the compounds, 4c exhibited excellent in vitro antimicrobial activity against all tested strains with binding and docking energies -35.6 and -12.4 kcal/mol, respectively. The antiviral studies were carried out for the selected compounds in which 4a exhibited 73.69 and 54.42 % of inhibition of buffalopox and camelpox viruses, respectively. Furthermore, compound 4a showed minimum docking and binding energy along with the maximum hydrogen/hydrophobic interaction with IMPDH. The study contributes towards identification and screening of potential antimicrobial and antiviral agent's against the pathogens.

Identification and characterization of novel indole based small molecules as anticancer agents through SIRT1 inhibition

In our pursuit to develop new potential anticancer leads, we designed a combination of structural units of indole and substituted triazole; and a library of 1-{1-methyl-2-4-phenyl-5-(propan-2-ylsulfanyl)-4H-1,2,4-triazol-3-yl ]-1H-indol-3-yl}methanamine derivatives was synthesized and characterized. Cytotoxic evaluations of these molecules over a panel of three human cancer cell lines were carried out. Few molecules exhibited potent growth inhibitory action against the treated cancer cell lines at lower micro molar concentration. An in vitro assay investigation of these active compounds using recombinant human SIRT1 enzyme showed that one of the compounds (IT-14) inhibited the deacetylation activity of the enzyme. The in vivo study of IT-14 exemplified its promising action by reducing the prostate weight to the body weight ratio in prostate hyperplasia animal models. A remarkable decrease in the disruption of histoarchitecture of the prostate tissues isolated from IT-14 treated animal compared to that of the positive control was observed. The molecular interactions with SIRT1 enzyme were also supported by molecular docking simulations. Hence this compound can act as a lead molecule to treat prostatic hyperplasia. (C) 2013 Elsevier Masson SAS. All rights reserved.

PLIC: protein-ligand interaction clusters

Most of the biological processes are governed through specific protein-ligand interactions. Discerning different components that contribute toward a favorable protein-ligand interaction could contribute significantly toward better understanding protein function, rationalizing drug design and obtaining design principles for protein engineering. The Protein Data Bank (PDB) currently hosts the structure of similar to 68 000 protein-ligand complexes. Although several databases exist that classify proteins according to sequence and structure, a mere handful of them annotate and classify protein-ligand interactions and provide information on different attributes of molecular recognition. In this study, an exhaustive comparison of all the biologically relevant ligand-binding sites (84 846 sites) has been conducted using PocketMatch: a rapid, parallel, in-house algorithm. PocketMatch quantifies the similarity between binding sites based on structural descriptors and residue attributes. A similarity network was constructed using binding sites whose PocketMatch scores exceeded a high similarity threshold (0.80). The binding site similarity network was clustered into discrete sets of similar sites using the Markov clustering (MCL) algorithm. Furthermore, various computational tools have been used to study different attributes of interactions within the individual clusters. The attributes can be roughly divided into (i) binding site characteristics including pocket shape, nature of residues and interaction profiles with different kinds of atomic probes, (ii) atomic contacts consisting of various types of polar, hydrophobic and aromatic contacts along with binding site water molecules that could play crucial roles in protein-ligand interactions and (iii) binding energetics involved in interactions derived from scoring functions developed for docking. For each ligand-binding site in each protein in the PDB, site similarity information, clusters they belong to and description of site attributes are provided as a relational database-protein-ligand interaction clusters (PLIC).