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.

The zinc finger proteins Mxr1p and repressor of phosphoenolpyruvate carboxykinase (ROP) have the same DNA binding specificity but regulate Methanol Metabolism antagonistically in pichia pastoris

The methanol-inducible alcohol oxidase I (AOXI) promoter of the methylotrophic yeast, Pichia pastoris, is used widely for the production of recombinant proteins. AOXI transcription is regulated by the zinc finger protein Mxr1p (methanol expression regulator 1). ROP (repressor of phosphoenolpyruvate carboxykinase, PEPCK) is a methanol- and biotin starvation-inducible zinc finger protein that acts as a negative regulator of PEPCK in P. pastoris cultured in biotin-deficient, glucose-ammonium medium. The function of ROP during methanol metabolism is not known. In this study, we demonstrate that ROP represses methanol-inducible expression of AOXI when P. pastoris is cultured in a nutrient-rich medium containing yeast extract, peptone, and methanol (YPM). Deletion of the gene encoding ROP results in enhanced expression of AOXI and growth promotion whereas overexpression of ROP results in repression of AOXI and growth retardation of P. pastoris cultured in YPM medium. Surprisingly, deletion or overexpression of ROP has no effect on AOXI gene expression and growth of P. pastoris cultured in a minimal medium containing yeast nitrogen base and methanol (YNBM). Subcellular localization studies indicate that ROP translocates from cytosol to nucleus of cells cultured in YPM but not YNBM. In vitro DNA binding studies indicate that AOXI promoter sequences containing 5' CYCCNY 3' motifs serve as binding sites for Mxr1p as well as ROP. Thus, Mxr1p and ROP exhibit the same DNA binding specificity but regulate methanol metabolism antagonistically in P. pastoris. This is the first report on the identification of a transcriptional repressor of methanol metabolism in any yeast species.

Depth: a web server to compute depth, cavity sizes, detect potential small-molecule ligand-binding cavities and predict the pK(a) of ionizable residues in proteins

Residue depth accurately measures burial and parameterizes local protein environment. Depth is the distance of any atom/residue to the closest bulk water. We consider the non-bulk waters to occupy cavities, whose volumes are determined using a Voronoi procedure. Our estimation of cavity sizes is statistically superior to estimates made by CASTp and VOIDOO, and on par with McVol over a data set of 40 cavities. Our calculated cavity volumes correlated best with the experimentally determined destabilization of 34 mutants from five proteins. Some of the cavities identified are capable of binding small molecule ligands. In this study, we have enhanced our depth-based predictions of binding sites by including evolutionary information. We have demonstrated that on a database (LigASite) of similar to 200 proteins, we perform on par with ConCavity and better than MetaPocket 2.0. Our predictions, while less sensitive, are more specific and precise. Finally, we use depth (and other features) to predict pK(a)s of GLU, ASP, LYS and HIS residues. Our results produce an average error of just <1 pH unit over 60 predictions. Our simple empirical method is statistically on par with two and superior to three other methods while inferior to only one. The DEPTH server (http://mspc.bii.a-star.edu.sg/depth/) is an ideal tool for rapid yet accurate structural analyses of protein structures.

Electron transfer amongst flavo- and hemo-proteins: diffusible species effect the relay processes, not protein-protein binding

Hitherto, electron transfer (ET) between redox proteins has been deemed to occur via donor-acceptor binding, and diffusible reactive species are considered as deleterious side-products in such systems. Herein, ET from cytochrome P450 reductase (CPR, an animal membrane flavoprotein) and horseradish peroxidase (HRP, a plant hemoprotein) to cytochrome c (Cyt c, a soluble animal hemoprotein) was probed under diverse conditions, using standard assays. ET in the CPR-Cyt c system was critically inhibited by cyanide and sub-equivalent levels of polar one-electron cyclers like copper ions, vitamin C/Trolox and superoxide dismutase. In the presence of lipids, inhibition was also afforded by amphipathic molecules vitamin E, palmitoyl-vitamin C and the membrane hemoprotein, cytochrome b(5). Such nonspecific inhibition (by diverse agents in both aqueous and lipid phases) indicated that electron transfer/relay was effected by small diffusible agents, whose lifetimes are shortened by the diverse radical scavengers. When CPR was retained in a dialysis membrane and Cyt c presented outside in free solution, ET was still observed. Further, HRP (taken at nM levels) catalyzed oxidation of a phenolic substrate was significantly inhibited upon the incorporation of sub-nM levels of Cyt c. The findings imply that CPR-Cyt c or HRP-Cyt c binding is not crucial for ET. Further, fundamental quantitative arguments (based on diffusion/collision) challenge the erstwhile protein-protein binding-assisted ET hypothesis. It is proven beyond reasonable doubt that mobile and diffusible electron carriers (ions and radicals) serve as ``redox-relay agents'' in the biological ET models/setup studied.

Orchestration of Haemophilus influenzae RecJ Exonuclease by Interaction with Single-Stranded DNA-Binding Protein

RecJ exonuclease plays crucial roles in several DNA repair and recombination pathways, and its ubiquity in bacterial species points to its ancient origin and vital cellular function. RecJ exonuclease from Haemophilus influenzae is a 575-amino-acid protein that harbors the characteristic motifs conserved among RecJ homologs. The purified protein exhibits a process 5'-3' single-stranded-DNA-specific exonuclease activity. The exonuclease activity of H. influenzae RecJ (HiRecJ) was supported by Mg2+ or Mn2+ and inhibited by Cd2+ suggesting a different mode of metal binding in HiRecJ as compared to Escherichia coli RecJ (EcoRecJ). Site-directed mutagenesis of highly conserved residues in HiRecJ abolished enzymatic activity. Interestingly, substitution of alanine for aspartate 77 resulted in a catalytically inactive enzyme that bound to DNA with a significantly higher affinity as compared to the wild-type enzyme. Noticeably, steady-state kinetic studies showed that H. influenzae single-stranded DNA-binding protein (HiSSB) increased the affinity of HiRecJ for single-stranded DNA and stimulated its exonuclease activity. HiSSB, whose C-terminal tail had been deleted, failed to enhance RecJ exonuclease activity. More importantly, HiRecJ was found to directly associate with its cognate single-stranded DNA-binding protein (SSB), as demonstrated by various in vitro assays, Interaction studies carried out with the truncated variants of HiRecJ and HiSSB revealed that the two proteins interact via the C-terminus of SSB protein and the core-catalytic domain of RecJ. Taken together, these results emphasize direct interactio between RecJ and SSB, which confers functional cooperativity to these two proteins. In addition, these results implicate SSB as being involved in the recruitment of RecJ to DNA and provide insights into the interplay between these proteins in repair and recombination pathways.

Structural mobility and transformations in globular proteins

Although globular proteins are endowed with well defined three-dimensional structures, they exhibit substantial mobility within the framework of the given threedimensional structure. The different types of mobility found in proteins by and large correspond to the different levels of organisational hierarchy in protein architecture. They are of considerable structural and functional significance, and can be broadly classified into(a) thermal and conformational fluctuations, (b) segmental mobility, (c) interdomain mobility and (d) intersubunit mobility. Protein crystallographic studies has provided a wealth of information on all of them. The temperature factors derived from X-ray diffraction studies provide a measure of atomic displacements caused by thermal and conformational fluctuations. The variation of displacement along the polypeptide chain have provided functionally significant information on the flexibility of different regions of the molecule in proteins such as myoglobin, lysozyme and prealbumin. Segmental mobility often involves the movement of a region or a segment of a molecule with respect to the rest, as in the transition between the apo and the holo structures of lactate dehydrogenase. It may also involve rigidification of a disordered region of the molecule as in the activation of the zymogens of serine proteases. Transitions between the apo and the holo structures of alcohol dehydrogenase,and between the free and the sugar bound forms of hexokinase, are good examples of interdomain mobility caused by hinge-bending. The capability of different domains to move semi-independently contributes greatly to the versatility of immunoglobulin molecules. Interdomain mobility in citrate synthase appears to be more complex and its study has led to an alternative description of domain closure. The classical and the most thoroughly studied case of intersubunit mobility is that in haemoglobin. The stereochemical mechanism of the action of this allosteric protein clearly brings out the functional subtilities that could be achieved through intersubunit movements. In addition to ligand binding and activation,environmental changes also often cause structural transformations. The reversible transformation between 2 Zn insulin and 4 Zn insulin is caused by changes in the ionic strength of the medium. Adenylate Kinase provides a good example for functionally significant reversible conformational transitions induced by variation in pH. Available evidences indicate that reversible structural transformations in proteins could also be caused by changes in the aqueous environment, including those in the amount of water surrounding protein molecules.

Crystallographic and modelling studies on Mycobacterium tuberculosis RuvA: Additional role of RuvB-binding domain and inter species variability

RuvA, along with RuvB, is involved in branch migration of heteroduplex DNA in homologous recombination. The structures of three new crystal forms of RuvA from Mycobacterium tuberculosis (MtRuvA) have been determined. The RuvB-binding domain is cleaved off in one of them. Detailed models of the complexes of octameric RuvA from different species with the Holliday junction have also been constructed. A thorough examination of the structures presented here and those reported earlier brings to light the hitherto unappreciated role of the RuvB-binding domain in determining inter-domain orientation and oligomerization. These structures also permit an exploration of the interspecies variability of structural features such as oligomerization and the conformation of the loop that carries the acidic pin, in terms of amino acid substitutions. These models emphasize the additional role of the RuvB-binding domain in Holliday junction binding. This role along with its role in oligomerization could have important biological implications.

Estrogen modulation of retinol-binding protein in immature chicks: Comparison with riboflavin carrier protein

The kinetics of estrogen (E) modulation of retinol-binding protein (RBP) production in the liver of immature chicks were compared with those governing de novo induction of riboflavin carrier protein (RCP) in the same tissue. A single dose of E markedly enhanced the plasma levels of RBP without any detectable lag period to reach peak value by 24 h and this was followed by a decline to attain the baseline by 4 days. There was no amplification of the response during secondary stimulation unlike the case with RCP induction. With multiple E administration, the 4-fold increased plasma RBP concentrations were sustained at a steady state during both primary and secondary stimulations, whereas concomitant RCP concentration progressively increased with each hormone administration and this response was further amplified during secondary stimulation. Unlike RCP induction, enhanced RBP accumulation was not strictly E dose dependent although a minimal threshold level of the steroid was required to elicit measurable response. Progesterone (P) could neither modulate nor substitute for E in enhancing plasma levels of either of the 2 proteins while the anti-estrogens, en- and zuclomifene citrate severely suppressed the production of both the proteins. RCP induction was completely inhibited by both α-amanitin and cycloheximide for prolonged periods while E-stimulated RBP production was affected only partially by α-amanitin. Likewise, cycloheximide inhibition of RBP accumulation followed a pattern similar to that of hepatic general protein synthesis.

Effect of Crowding Agents, Signal Peptide, and Chaperone SecB on the Folding and Aggregation of E. coli Maltose Binding Protein

SecB, a soluble cytosolic chaperone component of the Secexport pathway, binds to newly synthesized precursor proteins and prevents their premature aggregation and folding and subsequently targets them to the translocation machinery on the membrane. PreMBP, the precursor form of maltose binding protein, has a 26-residue signal sequence attached to the N-terminus of MBP and is a physiological substrate of SecB. We examine the effect of macromolecular crowding and SecB on the stability and refolding of denatured preMBP and MBP. PreMBP was less stable than MBP (ΔTm =7( 0.5 K) in both crowded and uncrowded solutions. Crowding did not cause any substantial changes in the thermal stability ofMBP(ΔTm=1(0.4 K) or preMBP (ΔTm=0(0.6 K), as observed in spectroscopically monitored thermal unfolding experiments. However, both MBP and preMBP were prone to aggregation while refolding under crowded conditions. In contrast to MBP aggregates, which were amorphous, preMBP aggregates form amyloid fibrils.Under uncrowded conditions, a molar excess of SecB was able to completely prevent aggregation and promote disaggregation of preformed aggregates of MBP. When a complex of the denatured protein and SecB was preformed, SecB could completely prevent aggregation and promote folding of MBP and preMBP even in crowded solution. Thus, in addition to maintaining substrates in an unfolded, export-competent conformation, SecB also suppresses the aggregation of its substrates in the crowded intracellular environment. SecB is also able to promote passive disaggregation of macroscopic aggregates of MBP in the absence of an energy source such as ATP or additional cofactors. These experiments also demonstrate that signal peptide can reatly influence protein stability and aggregation propensity.

Estrogen Induction Of Biotin-Binding Protein In Immature Chicks - Kinetics, Hormonal Specificity And Modulation

Employing a specific radioimmunoassay for quantification, the kinetics of estrogen-induced elevation in the plasma concentration of biotin-binding protein (BBP) in immature male chicks was investigated. A single injection of the steroid hormone enhanced the plasma BBP content several-fold at 6 h, reaching peak levels around 48 h and declining thereafter. A 2-fold amplification of the response was evident during secondary stimulation with the hormone. The magnitude of the response was hormonal dose-dependent while the initial lag phase and the time of peak protein accumulation were unaltered within the hormonal doses tested. The circulatory half-life of the specific protein in normal and estrogenized birds was 10 h. Hyperthyroidism markedly decreased the hormonal response while the opposite effect was seen during hypothyroidism. The antiestrogens E- and Z-clomiphene citrate effectively blocked the protein induction whereas progesterone, either alone or in combination with estrogen, was ineffective in modulating the induction. Cycloheximide administration drastically inhibited the inductive response. The above observations clearly suggest that the genes corresponding to the two isofunctional proteins of chicken egg, viz. BBP and avidin, are differentially regulated.

Analysis of proteins with the `hot dog' fold: Prediction of function and identification of catalytic residues of hypothetical proteins

Background: The hot dog fold has been found in more than sixty proteins since the first report of its existence about a decade ago. The fold appears to have a strong association with fatty acid biosynthesis, its regulation and metabolism, as the proteins with this fold are predominantly coenzyme A-binding enzymes with a variety of substrates located at their active sites. Results: We have analyzed the structural features and sequences of proteins having the hot dog fold. This study reveals that though the basic architecture of the fold is well conserved in these proteins, significant differences exist in their sequence, nature of substrate and oligomerization. Segments with certain conserved sequence motifs seem to play crucial structural and functional roles in various classes of these proteins. Conclusion: The analysis led to predictions regarding the functional classification and identification of possible catalytic residues of a number of hot dog fold-containing hypothetical proteins whose structures were determined in high throughput structural genomics projects.

CSSI-PRO: a method for secondary structure type editing, assignment and estimation in proteins using linear combination of backbone chemical shifts

Estimation of secondary structure in polypeptides is important for studying their structure, folding and dynamics. In NMR spectroscopy, such information is generally obtained after sequence specific resonance assignments are completed. We present here a new methodology for assignment of secondary structure type to spin systems in proteins directly from NMR spectra, without prior knowledge of resonance assignments. The methodology, named Combination of Shifts for Secondary Structure Identification in Proteins (CSSI-PRO), involves detection of specific linear combination of backbone H-1(alpha) and C-13' chemical shifts in a two-dimensional (2D) NMR experiment based on G-matrix Fourier transform (GFT) NMR spectroscopy. Such linear combinations of shifts facilitate editing of residues belonging to alpha-helical/beta-strand regions into distinct spectral regions nearly independent of the amino acid type, thereby allowing the estimation of overall secondary structure content of the protein. Comparison of the predicted secondary structure content with those estimated based on their respective 3D structures and/or the method of Chemical Shift Index for 237 proteins gives a correlation of more than 90% and an overall rmsd of 7.0%, which is comparable to other biophysical techniques used for structural characterization of proteins. Taken together, this methodology has a wide range of applications in NMR spectroscopy such as rapid protein structure determination, monitoring conformational changes in protein-folding/ligand-binding studies and automated resonance assignment.