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.

Light-induced geometric isomerization of 1,2-diphenylcyclopropanes included within Y zeolites: Role of cation-guest binding

Through a systematic study of several diphenylcyclopropane derivatives, we have inferred that the cations present within a zeolite control the excited-state chemistry of these systems. In the parent 1,2-diphenylcylopropane, the cation binds to the two phenyl rings in a sandwich-type arrangement, and such a mode of binding prevents cis-to-trans isomerization. Once an ester or amide group is introduced into the system (derivatives of 2beta,3beta-diphenylcyclopropane-1alpha-carboxylic acid), the cation binds to the carbonyl group present in these chromophores and such a binding has no influence on the cis-trans isomerization process. Cation-reactant structures computed at density functional theory level have been very valuable in rationalizing the observed photochemical behavior of diphenylcyclopropane derivatives included in zeolites. While the parent system, 1,2-diphenyleylopropane, has been extensively investigated in the context of chiral induction in solution, owing to its failure to isomerize from cis to trans, the same could not be investigated in zeolites. However, esters of 2beta,3beta-diphenylcyclopropane-1alpha-carboxylic acid could be studied within zeolites in the context of chiral induction. Chiral induction as high 20% ee and 55% de has been obtained with selected systems. These numbers, although low, are much higher than what has been obtained in solution with the same system or with the parent system by other investigators (maximum similar to10% ee).

A molecular dynamics study based post facto free energy analysis of the binding of bovine angiogenin with UMP and CMP ligands

Angiogenin is a protein belonging to the superfamily of RNase A. The RNase activity of this protein is essential for its angiogenic activity. Although members of the RNase A family carry out RNase activity, they differ markedly in their strength and specificity. In this paper, we address the problem of higher specificity of angiogenin towards cytosine against uracil in the first base binding position. We have carried out extensive nano-second level molecular dynamics(MD) computer simulations on the native bovine angiogenin and on the CMP and UMP complexes of this protein in aqueous medium with explicit molecular solvent. The structures thus generated were subjected to a rigorous free energy component analysis to arrive at a plausible molecular thermodynamic explanation for the substrate specificity of angiogenin.

Denaturant mediated unfolding of both native and molten globule states of maltose binding protein are accompanied by large [Delta]Cp's

Maltose binding protein (MBP) is a large, monomeric two domain protein containing 370 amino acids. In the absence of denaturant at neutral pH, the protein is in the native state, while at pH 3.0 it forms a molten globule. The molten globule lacks a tertiary circular dichroism signal but has secondary structure similar to that of the native state. The molten globule binds 8-anilino-1-naphthalene sulfonate (ANS). The unfolding thermodynamics of MBP at both pHs were measured by carrying out a series of isothermal urea melts at temperatures ranging from 274–329 K. At 298 K, values of [Delta]G°, [Delta]Cp, and Cm were 3.1 ± 0.2 kcal mol−1, 5.9 ± 0.8 kcal mol−1 K−1 (15.9 cal (mol-residue)−1 K−1), and 0.8 M, respectively, at pH 3.0 and 14.5 ± 0.4 kcal mol−1, 8.3 ± 0.7 kcal mol−1 K−1 (22.4 kcal (mol-residue)−1 K−1), and 3.3 M, respectively, at pH 7.1. Guanidine hydrochloride denaturation at pH 7.1 gave values of [Delta]G° and [Delta]Cp similar to those obtained with urea. The m values for denaturation are strongly temperature dependent, in contrast to what has been previously observed for small globular proteins. The value of [Delta]Cp per mol-residue for the molten globule is comparable to corresponding values of [Delta]Cp for the unfolding of typical globular proteins and suggests that it is a highly ordered structure, unlike molten globules of many small proteins. The value of [Delta]Cp per mol-residue for the unfolding of the native state is among the highest currently known for any protein.

PocketAlign A Novel Algorithm for Aligning Binding Sites in Protein Structures

A fundamental task in bioinformatics involves a transfer of knowledge from one protein molecule onto another by way of recognizing similarities. Such similarities are obtained at different levels, that of sequence, whole fold, or important substructures. Comparison of binding sites is important to understand functional similarities among the proteins and also to understand drug cross-reactivities. Current methods in literature have their own merits and demerits, warranting exploration of newer concepts and algorithms, especially for large-scale comparisons and for obtaining accurate residue-wise mappings. Here, we report the development of a new algorithm, PocketAlign, for obtaining structural superpositions of binding sites. The software is available as a web-service at http://proline.physicslisc.emetin/pocketalign/. The algorithm encodes shape descriptors in the form of geometric perspectives, supplemented by chemical group classification. The shape descriptor considers several perspectives with each residue as the focus and captures relative distribution of residues around it in a given site. Residue-wise pairings are computed by comparing the set of perspectives of the first site with that of the second, followed by a greedy approach that incrementally combines residue pairings into a mapping. The mappings in different frames are then evaluated by different metrics encoding the extent of alignment of individual geometric perspectives. Different initial seed alignments are computed, each subsequently extended by detecting consequential atomic alignments in a three-dimensional grid, and the best 500 stored in a database. Alignments are then ranked, and the top scoring alignments reported, which are then streamed into Pymol for visualization and analyses. The method is validated for accuracy and sensitivity and benchmarked against existing methods. An advantage of PocketAlign, as compared to some of the existing tools available for binding site comparison in literature, is that it explores different schemes for identifying an alignment thus has a better potential to capture similarities in ligand recognition abilities. PocketAlign, by finding a detailed alignment of a pair of sites, provides insights as to why two sites are similar and which set of residues and atoms contribute to the similarity.

Insulin Growth Factor-2 Binding Protein 3 (IGF2BP3) Is a Glioblastoma-specific Marker That Activates Phosphatidylinositol 3-Kinase/Mitogen-activated Protein Kinase (PI3K/MAPK) Pathways by Modulating IGF-2

Glioblastoma is the most common and malignant form of primary astrocytoma. Upon investigation of the insulin-like growth factor (IGF) pathway, we found the IGF2BP3/IMP3 transcript and protein to be up-regulated in GBMs but not in lower grade astrocytomas (p<0.0001). IMP3 is an RNA binding protein known to bind to the 5'-untranslated region of IGF-2 mRNA, thereby activating its translation. Overexpression-and knockdown-based studies establish a role for IMP3 in promoting proliferation, anchorage-independent growth, invasion, and chemoresistance. IMP3 overexpressing B16F10 cells also showed increased tumor growth, angiogenesis, and metastasis, resulting in poor survival in a mouse model. Additionally, the infiltrating front, perivascular, and subpial regions in a majority of the GBMs stained positive for IMP3. Furthermore, two different murine glioma models were used to substantiate the above findings. In agreement with the translation activation functions of IMP3, we also found increased IGF-2 protein in the GBM tumor samples without a corresponding increase in its transcript levels. Also, in vitro IMP3 overexpression/knockdown modulated the IGF-2 protein levels without altering its transcript levels. Additionally, IGF-2 neutralization and supplementation studies established that the proproliferative effects of IMP3 were indeed mediated through IGF-2. Concordantly, PI3K and MAPK, the downstream effectors of IGF-2, are activated by IMP3 and are found to be essential for IMP3-induced cell proliferation. Thus, we have identified IMP3 as a GBM-specific proproliferative and proinvasive marker acting through IGF-2 resulting in the activation of oncogenic PI3K and MAPK pathways.

Ligand-modulated Parallel Mechanical Unfolding Pathways of Maltose-binding Proteins

Protein folding and unfolding are complex phenomena, and it is accepted that multidomain proteins generally follow multiple pathways. Maltose-binding protein (MBP) is a large (a two-domain, 370-amino acid residue) bacterial periplasmic protein involved in maltose uptake. Despite the large size, it has been shown to exhibit an apparent two-state equilibrium unfolding in bulk experiments. Single-molecule studies can uncover rare events that are masked by averaging in bulk studies. Here, we use single-molecule force spectroscopy to study the mechanical unfolding pathways of MBP and its precursor protein (preMBP) in the presence and absence of ligands. Our results show that MBP exhibits kinetic partitioning on mechanical stretching and unfolds via two parallel pathways: one of them involves a mechanically stable intermediate (path I) whereas the other is devoid of it (path II). The apoMBP unfolds via path I in 62% of the mechanical unfolding events, and the remaining 38% follow path II. In the case of maltose-bound MBP, the protein unfolds via the intermediate in 79% of the cases, the remaining 21% via path II. Similarly, on binding to maltotriose, a ligand whose binding strength with the polyprotein is similar to that of maltose, the occurrence of the intermediate is comparable (82% via path I) with that of maltose. The precursor protein preMBP also shows a similar behavior upon mechanical unfolding. The percentages of molecules unfolding via path I are 53% in the apo form and 68% and 72% upon binding to maltose and maltotriose, respectively, for preMBP. These observations demonstrate that ligand binding can modulate the mechanical unfolding pathways of proteins by a kinetic partitioning mechanism. This could be a general mechanism in the unfolding of other large two-domain ligand-binding proteins of the bacterial periplasmic space.

Synthesis and DNA binding studies of Ni(II), Co(II), Cu(II) and Zn(II) metal complexes of N(1),N(5)-bis[pyridine-2-methylene]-thiocarbohydrazone Schiff-base ligand

The thiocarbohydrazone Schiff-base ligand with a nitrogen and sulphur donor was synthesized through condensation of pyridine-2-carbaldehyde and thiocarbohydrazide. Schiff-base ligands have the ability to conjugate with metal salts. A series of metal complexes with a general formula [MCl(2)(H(2)L)]center dot nH(2)O (M=Ni, Co, Cu and Zn) were synthesized by forming complexes of the N(1),N5-bis[pyridine-2-methylene]thiocarbohydrazone (H2L) Schiff-base ligand. These metal complexes and ligand were characterized by using ultraviolet-visible (UV-Vis), Fourier Transform Infrared (FT-IR), (1)H and (13)C NMR spectroscopy and mass spectroscopy, physicochemical characterization, CHNS and conductivity. The biological activity of the synthesized ligand was investigated by using Escherichia coli DNA as target. The DNA interaction of the synthesized ligand and complexes on E. coli plasmid DNA was investigated in the aqueous medium by UV-Vis spectroscopy and the binding constant (K(b)) was calculated. The DNA binding studies showed that the metal complexes had an improved interaction due to trans-geometrical isomers of the complexes than ligand isomers in cis-positions. (C) 2011 Elsevier B.V. All rights reserved.

Conformational changes of the chaperone SecB upon binding to a model substrate - bovine pancreatic trypsin inhibitor (BPTI)

SecB is a homotetrameric cytosolic chaperone that forms part of the protein translocation machinery in E. coli. Due to SecB, nascent polypeptides are maintained in an unfolded translocation-competent state devoid of tertiary structure and thus are guided to the translocon. In vitro SecB rapidly binds to a variety of ligands in a non-native state. We have previously investigated the bound state conformation of the model substrate bovine pancreatic trypsin inhibitor (BPTI) as well as the conformation of SecB itself by using proximity relationships based on site-directed spin labeling and pyrene fluorescence methods. It was shown that SecB undergoes a conformational change during the process of substrate binding. Here, we generated SecB mutants containing but a single cysteine per subunit or an exposed highly reactive new cysteine after removal of the nearby intrinsic cysteines. Quantitative spin labeling was achieved with the methanethiosulfonate spin label (MTS) at positions C97 or E90C, respectively. Highfield (W-band) electron paramagnetic resonance (EPR) measurements revealed that with BPTI present the spin labels are exposed to a more polar/hydrophilic environment. Nanoscale distance measurements with double electron-electron resonance (DEER) were in excellent agreement with distances obtained by molecular modeling. Binding of BPTI also led to a slight change in distances between labels at C97 but not at E90C. While the shorter distance in the tetramer increased, the larger diagonal distance decreased. These findings can be explained by a widening of the tetrameric structure upon substrate binding much like the opening of two pairs of scissors.

Synergy between the N-terminal and C-terminal domains of Mycobacterium tuberculosis HupB is essential for high-affinity binding, DNA supercoiling and inhibition of RecA-promoted strand exchange

The occurrence of DNA architectural proteins containing two functional domains derived from two different architectural proteins is an interesting emerging research theme in the field of nucleoid structure and function. Mycobacterium tuberculosis HupB, unlike Escherichia coli HU, is a two-domain protein that, in the N-terminal region, shows broad sequence homology with bacterial HU. The long C-terminal extension, on the other hand, contains seven PAKK/KAAK motifs, which are characteristic of the histone H1/H5 family of proteins. In this article, we describe several aspects of HupB function, in comparison with its truncated derivatives lacking either the C-terminus or N-terminus. We found that HupB binds a variety of DNA repair and replication intermediates with K(d) values in the nanomolar range. By contrast, the N-terminal fragment of M. tuberculosis HupB (HupB(MtbN)) showed diminished DNA-binding activity, with K(d) values in the micromolar range, and the C-terminal domain was completely devoid of DNA-binding activity. Unlike HupB(MtbN), HupB was able to constrain DNA in negative supercoils and introduce negative superhelical turns into relaxed DNA. Similarly, HupB exerted a robust inhibitory effect on DNA strand exchange promoted by cognate and noncognate RecA proteins, whereas HupB(MtbN), even at a 50-fold molar excess, had no inhibitory effect. Considered together, these results suggest that synergy between the N-terminal and C-terminal domains of HupB is essential for its DNA-binding ability, and to modulate the topological features of DNA, which has implications for processes such as DNA compaction, gene regulation, homologous recombination, and DNA repair.

Investigations on the RNA binding and phosphorylation of groundnut bud necrosis virus nucleocapsid protein

Groundnut bud necrosis virus belongs to the genus Tospovirus, infects a wide range of crop plants and causes severe losses. To understand the role of the nucleocapsid protein in the viral life cycle, the protein was overexpressed in E. coli and purified by Ni-NTA chromatography. The purified N protein was well folded and was predominantly alpha-helical. Deletion analysis revealed that the C-terminal unfolded region of the N protein was involved in RNA binding. Furthermore, the N protein could be phosphorylated in vitro by Nicotiana benthamiana plant sap and by purified recombinant kinases such as protein kinase CK2 and calcium-dependent protein kinase. This is the first report of phoshphorylation of a nucleocapsid protein in the family Bunyaviridae. The possible implications of the present findings for the viral life cycle are discussed.

Chimeras of Escherichia coli and Mycobacterium tuberculosis Single-Stranded DNA Binding Proteins: Characterization and Function in Escherichia coli

Single stranded DNA binding proteins (SSBs) are vital for the survival of organisms. Studies on SSBs from the prototype, Escherichia coli (EcoSSB) and, an important human pathogen, Mycobacterium tuberculosis (MtuSSB) had shown that despite significant variations in their quaternary structures, the DNA binding and oligomerization properties of the two are similar. Here, we used the X-ray crystal structure data of the two SSBs to design a series of chimeric proteins (m beta 1, m beta 1'beta 2, m beta 1-beta 5, m beta 1-beta 6 and m beta 4-beta 5) by transplanting beta 1, beta 1'beta 2, beta 1-beta 5, beta 1-beta 6 and beta 4-beta 5 regions, respectively of the N-terminal (DNA binding) domain of MtuSSB for the corresponding sequences in EcoSSB. In addition, m beta 1'beta 2(ESWR) SSB was generated by mutating the MtuSSB specific `PRIY' sequence in the beta 2 strand of m beta 1'beta 2 SSB to EcoSSB specific `ESWR' sequence. Biochemical characterization revealed that except for m beta 1 SSB, all chimeras and a control construct lacking the C-terminal domain (Delta C SSB) bound DNA in modes corresponding to limited and unlimited modes of binding. However, the DNA on MtuSSB may follow a different path than the EcoSSB. Structural probing by protease digestion revealed that unlike other SSBs used, m beta 1 SSB was also hypersensitive to chymotrypsin treatment. Further, to check for their biological activities, we developed a sensitive assay, and observed that m beta 1-beta 6, MtuSSB, m beta 1'beta 2 and m beta 1-beta 5 SSBs complemented E. coli Delta ssb in a dose dependent manner. Complementation by the m beta 1-beta 5 SSB was poor. In contrast, m beta 1'beta 2(ESWR) SSB complemented E. coli as well as EcoSSB. The inefficiently functioning SSBs resulted in an elongated cell/filamentation phenotype of E. coli. Taken together, our observations suggest that specific interactions within the DNA binding domain of the homotetrameric SSBs are crucial for their biological function.