Interaction of water molecules in non-identical protein structures

The interaction of the protein atoms with the surrounding water oxygen atoms has been computed for 392 protein chains from 369 protein structures belonging to 90% non-homologous high resolution (<= 1.5 angstrom) protein Structures with a crystallographic R-factor <= 20%. The percentage composition of the polar atoms is found to be 36.3%. An average of 82.55% of water oxygen atoms are found to be in the primary hydration shell and 15.12% in the secondary hydration shell. The average Percentage of interactions of water oxygen atoms with the polar atoms of the main chain and side chain are 54% and 46%. respectively. The interaction of the acidic residues, aspartate and glutamate, with the water oxygen atoms is more when compared to that of the other residues.

Introduction of high-content screening method for studying drug-induced mitochondrial dysfunction in hepatocytes

Drug induced liver injury is one of the frequent reasons for the drug removal from the market. During the recent years there has been a pressure to develop more cost efficient, faster and easier ways to investigate drug-induced toxicity in order to recognize hepatotoxic drugs in the earlier phases of drug development. High Content Screening (HCS) instrument is an automated microscope equipped with image analysis software. It makes the image analysis faster and decreases the risk for an error caused by a person by analyzing the images always in the same way. Because the amount of drug and time needed in the analysis are smaller and multiple parameters can be analyzed from the same cells, the method should be more sensitive, effective and cheaper than the conventional assays in cytotoxicity testing. Liver cells are rich in mitochondria and many drugs target their toxicity to hepatocyte mitochondria. Mitochondria produce the majority of the ATP in the cell through oxidative phosphorylation. They maintain biochemical homeostasis in the cell and participate in cell death. Mitochondria is divided into two compartments by inner and outer mitochondrial membranes. The oxidative phosphorylation happens in the inner mitochondrial membrane. A part of the respiratory chain, a protein called cytochrome c, activates caspase cascades when released. This leads to apoptosis. The aim of this study was to implement, optimize and compare mitochondrial toxicity HCS assays in live cells and fixed cells in two cellular models: human HepG2 hepatoma cell line and rat primary hepatocytes. Three different hepato- and mitochondriatoxic drugs (staurosporine, rotenone and tolcapone) were used. Cells were treated with the drugs, incubated with the fluorescent probes and then the images were analyzed using Cellomics ArrayScan VTI reader. Finally the results obtained after optimizing methods were compared to each other and to the results of the conventional cytotoxicity assays, ATP and LDH measurements. After optimization the live cell method and rat primary hepatocytes were selected to be used in the experiments. Staurosporine was the most toxic of the three drugs and caused most damage to the cells most quickly. Rotenone was not that toxic, but the results were more reproducible and thus it would serve as a good positive control in the screening. Tolcapone was the least toxic. So far the conventional analysis of cytotoxicity worked better than the HCS methods. More optimization needs to be done to get the HCS method more sensitive. This was not possible in this study due to time limit.

Empirical torsional potential functions from protein structure data. Phi- and psi-potentials for non-glycyl amino acid residues.

The torsional potential functions Vt(phi) and Vt(psi) around single bonds N--C alpha and C alpha--C, which can be used in conformational studies of oligopeptides, polypeptides and proteins, have been derived, using crystal structure data of 22 globular proteins, fitting the observed distribution in the (phi, psi)-plane with the value of Vtot(phi, psi), using the Boltzmann distribution. The averaged torsional potential functions, obtained from various amino acid residues in L-configuration, are Vt(phi) = 1.0 cos (phi + 60 degrees); Vt(psi) = 0.5 cos (psi + 60 degrees) - 1.0 cos (2 psi + 30 degrees) - 0.5 cos (3 psi + 30 degrees). The dipeptide energy maps Vtot(phi, psi) obtained using these functions, instead of the normally accepted torsional functions, were found to explain various observations, such as the absence of the left-handed alpha helix and the C7 conformation, and the relatively high density of points near the line psi = 0 degrees. These functions derived from observational data on protein structures, will, it is hoped, explain various previously unexplained facts in polypeptide conformation.

Over-expression and purification strategies for recombinant multi-protein oligomers: A case study of Mycobacterium tuberculosis sigma/anti-sigma factor protein complexes

The function of a protein in a cell often involves coordinated interactions with one or several regulatory partners. It is thus imperative to characterize a protein both in isolation as well as in the context of its complex with an interacting partner. High resolution structural information determined by X-ray crystallography and Nuclear Magnetic Resonance offer the best route to characterize protein complexes. These techniques, however, require highly purified and homogenous protein samples at high concentration. This requirement often presents a major hurdle for structural studies. Here we present a strategy based on co-expression and co-purification to obtain recombinant multi-protein complexes in the quantity and concentration range that can enable hitherto intractable structural projects. The feasibility of this strategy was examined using the sigma factor/anti-sigma factor protein complexes from Mycobacterium tuberculosis. The approach was successful across a wide range of sigma factors and their cognate interacting partners. It thus appears likely that the analysis of these complexes based on variations in expression constructs and procedures for the purification and characterization of these recombinant protein samples would be widely applicable for other multi-protein systems. (C) 2010 Elsevier Inc. All rights reserved.

PeptideMine - A webserver for the design of peptides for protein-peptide binding studies derived from protein-protein interactomes

Background: Signal transduction events often involve transient, yet specific, interactions between structurally conserved protein domains and polypeptide sequences in target proteins. The identification and validation of these associating domains is crucial to understand signal transduction pathways that modulate different cellular or developmental processes. Bioinformatics strategies to extract and integrate information from diverse sources have been shown to facilitate the experimental design to understand complex biological events. These methods, primarily based on information from high-throughput experiments, have also led to the identification of new connections thus providing hypothetical models for cellular events. Such models, in turn, provide a framework for directing experimental efforts for validating the predicted molecular rationale for complex cellular processes. In this context, it is envisaged that the rational design of peptides for protein-peptide binding studies could substantially facilitate the experimental strategies to evaluate a predicted interaction. This rational design procedure involves the integration of protein-protein interaction data, gene ontology, physico-chemical calculations, domain-domain interaction data and information on functional sites or critical residues. Results: Here we describe an integrated approach called ``PeptideMine'' for the identification of peptides based on specific functional patterns present in the sequence of an interacting protein. This approach based on sequence searches in the interacting sequence space has been developed into a webserver, which can be used for the identification and analysis of peptides, peptide homologues or functional patterns from the interacting sequence space of a protein. To further facilitate experimental validation, the PeptideMine webserver also provides a list of physico-chemical parameters corresponding to the peptide to determine the feasibility of using the peptide for in vitro biochemical or biophysical studies. Conclusions: The strategy described here involves the integration of data and tools to identify potential interacting partners for a protein and design criteria for peptides based on desired biochemical properties. Alongside the search for interacting protein sequences using three different search programs, the server also provides the biochemical characteristics of candidate peptides to prune peptide sequences based on features that are most suited for a given experiment. The PeptideMine server is available at the URL: http://caps.ncbs.res.in/peptidemine

Chloroquine Inhibits Heme-Dependent Protein Synthesis in Plasmodium falciparum

A cell-free protein-synthesizing system has been reconstituted using the S-30 fraction or ribosomes and the S-100 fraction from Plasmodium falciparum. Addition of heme in vitro stimulates cell-free protein synthesis strikingly. Chloroquine inhibits the heme-dependent protein synthesis in the parasite lysate. The drug has also been found to inhibit parasite protein synthesis in situ at therapeutic concentrations soon after addition to parasite cultures. Ribosomes as well as the S-100 fraction isolated from such chloroquine-treated cultures are defective in protein synthesis. Addition of hemin plus glucose 6-phosphate or high concentrations of GTP, cAMP, and an active preparation of eIF-2 to the parasite cell-free system restores protein synthesis to a significant extent in chloroquine-treated cultures. Under conditions of inhibition of protein synthesis in situ by chloroquine in the culture, the parasite eukaryotic initiation factor 2-alpha- (eIF-2-alpha) is phosphorylated in the parasite lysate to a greater extent than that observed in the control culture. Addition of hemin in vitro suppresses this phosphorylation. eIF-2-alpha kinase activity is present in the parasite lysate and is not a contaminant derived from the human erythrocytes used to culture the parasite. The heme-chloroquine interactive effects can also be demonstrated with purified eIF-2-alpha kinase from rabbit reticulocyte lysate. It is proposed that chloroquine inhibits heme-dependent protein synthesis in the parasite and this is an early event mediating the growth-inhibitory effects of the drug.

Isolation and characterization of thiamin-binding protein from chicken egg white

A thiamin-binding protein was isolated and characterized from chicken egg white by affinity chromatography on thiamin pyrophosphate coupled to aminoethyl-Sepharose. The high specificity of interaction between the thiamin-binding protein and the riboflavin-binding protein of the egg white, with a protein/protein molar ratio of 1.0, led to the development of an alternative procedure that used the riboflavin-binding protein immobilized on CNBr-activated Sepharose as the affinity matrix. The thiamin-binding protein thus isolated was homogeneous by the criteria of polyacrylamide-gel disc electrophoresis, double immunodiffusion and sodium dodecyl sulphate/polyacrylamide-gel electrophoresis, had a mol.wt. of 38,000 +/- 2000 and was not a glycoprotein. The protein bound [14C]thiamin was a molar ratio of 1.0, with dissociation constant (Kd) 0.3 micrometer.

A Deterministic Optimization Approach to Protein Sequence Design Using Continuous Models

Determining the sequence of amino acid residues in a heteropolymer chain of a protein with a given conformation is a discrete combinatorial problem that is not generally amenable for gradient-based continuous optimization algorithms. In this paper we present a new approach to this problem using continuous models. In this modeling, continuous "state functions" are proposed to designate the type of each residue in the chain. Such a continuous model helps define a continuous sequence space in which a chosen criterion is optimized to find the most appropriate sequence. Searching a continuous sequence space using a deterministic optimization algorithm makes it possible to find the optimal sequences with much less computation than many other approaches. The computational efficiency of this method is further improved by combining it with a graph spectral method, which explicitly takes into account the topology of the desired conformation and also helps make the combined method more robust. The continuous modeling used here appears to have additional advantages in mimicking the folding pathways and in creating the energy landscapes that help find sequences with high stability and kinetic accessibility. To illustrate the new approach, a widely used simplifying assumption is made by considering only two types of residues: hydrophobic (H) and polar (P). Self-avoiding compact lattice models are used to validate the method with known results in the literature and data that can be practically obtained by exhaustive enumeration on a desktop computer. We also present examples of sequence design for the HP models of some real proteins, which are solved in less than five minutes on a single-processor desktop computer Some open issues and future extensions are noted.

Homeostatic regulation of free retinol-binding protein and free thyroxine pools of plasma by their plasma carrier proteins in chicken

The mechanism underlying homeostatic regulation of the plasma levels of free retinol-binding protein and free thyroxine, the systemic distribution of which is of great importance, has been investigated. A simple method has been developed to determine the rate of dissociation of a ligand from the binding protein. Analysis of the dissociation process of retinol-binding protein from prealbumin-2 reveals that the free retinol-binding protein pool undergoes massive flux, and the prealbumin-2 participates in homeostatic regulation of the free retinol-binding protein pool. Studies on the dissociation process of thyroxine from its plasma carrier proteins show that the various plasma carrier proteins share two roles. Of the two types of protein, the thyroxine-binding globulin (the high affinity binding protein) contributes only 27% of the free thyroxine in a rapid transition process, despite its being the major binding protein. But prealbumin-2, which has lower affinity towards thyroxine, participates mainly in a rapid flux of the free thyroxine pool. Thus thyroxine-binding globulin acts predominantly as a plasma reservoir of thyroxine, and also probably in the �buffering� action on plasma free thyroxine level, in the long term, while prealbumin-2 participates mainly in the maintainance of constancy of free thyroxine levels even in the short term. The existence of these two types of binding protein facilitates compensation for the metabolic flux of the free ligand and maintenance of the thyroxine pool within a very narrow range.

High-quality annotation of promoter regions for 913 bacterial genomes

Motivation: The number of bacterial genomes being sequenced is increasing very rapidly and hence, it is crucial to have procedures for rapid and reliable annotation of their functional elements such as promoter regions, which control the expression of each gene or each transcription unit of the genome. The present work addresses this requirement and presents a generic method applicable across organisms. Results: Relative stability of the DNA double helical sequences has been used to discriminate promoter regions from non-promoter regions. Based on the difference in stability between neighboring regions, an algorithm has been implemented to predict promoter regions on a large scale over 913 microbial genome sequences. The average free energy values for the promoter regions as well as their downstream regions are found to differ, depending on their GC content. Threshold values to identify promoter regions have been derived using sequences flanking a subset of translation start sites from all microbial genomes and then used to predict promoters over the complete genome sequences. An average recall value of 72% (which indicates the percentage of protein and RNA coding genes with predicted promoter regions assigned to them) and precision of 56% is achieved over the 913 microbial genome dataset.

Microsomal redox systems in brown adipose tissue: high lipid peroxidation, low cholesterol biosynthesis and no detectable cytochrome P-450

The presence of redox systems in microsomes of brown adipose tissue (BAT) in cold exposed rats was investigated and compared with liver. BAT microsomes showed high activity of lipid peroxidation measured both by the formation of malondialdehyde (MDA) and by oxygen uptake. NADH and NADPH dependent cytochrome c reductase activity were present in both BAT and liver microsomes. Aminopyrine demethylase and aniline hydroxylase activities, the characteristic detoxification enzymes in liver microsomes could not be detected in BAT microsomes. BAT minces showed very poor incorporation of [1-14C]acetate and [2-14C]-mevalonate in unsaponifiable lipid fraction compared to liver. Biosynthesis of cholesterol and ubiquinone, but not fatty acids, and the activity of 3-hydroxy-3-methyl glutaryl CoA reductase appear to be very low in BAT. Examination of difference spectra showed the presence of only cytochrome b 5 in BAT microsomes. In addition to the inability to detect the enzyme activities dependent on cytochrome P-450, a protein with the characteristic spectrum, molecular size in SDS-PAGE and interaction with antibodies in double diffusion test, also could not be detected in BAT microsomes. The high activity of lipid peroxidation in microsomes, being associated with large oxygen uptake and oxidation of NADPH, will also contribute to the energy dissipation as heat in BAT, considered important in thermogenesis.

Characterization and immune response of a protein produced by a cDNA clone of foot and mouth disease virus, type Asia 1 63/72

A 0.9 kb double stranded cDNA of foot and mouth disease virus (FMDV) Type Asia 1, 63/72 was cloned in an expression vector, pUR222. A protein of 38 kd was produced by the clone which reacted with the antibodies raised against the virus. A 20 kd protein which may be derived from the 38 kd protein contained the antigenic epitopes of the protein VP1 of the virus. Injection of 10-20 micrograms of the partially purified 38 and 20 kd proteins or a lysate of cells containing 240 micrograms of the proteins elicited high titers of FMDV specific antibodies in guinea pigs and cattle respectively. Also, at these concentrations, the proteins protected 5 of 8 guinea pigs and 3 of 8 cattle when challenged with a virulent virus.

The variants of the protein toxins abrin and ricin A useful guide to understanding the processing events in the toxin transport

Kinetic data on inhibition of protein synthesis in thymocyte by three abrins and ricin have been obtained. The intrinsic efficiencies of A chains of four toxins to inactivate ribosomes, as analyzed by k1-versus-concentration plots were abrin II, III > ricin > abrin I. The lag times were 90, 66, 75 and 105 min at a 0.0744 nM concentration of each of abrin I, II, III and ricin, respectively. To account for the observed differences in the dose-dependent lag time, functional and structural variables of toxins such as binding efficiency of B chains to receptors and low-pH-induced structural alterations have been analyzed. The association constants obtained by stopped flow studies showed that abrin-I (4.13 × 105 M−1 s−1) association with putative receptor (4-methylumbelliferyl-α-D-galactoside) is nearly two times more often than abrin III (2.6 × 105 M−1 s−1) at 20°C. Equillibrium binding constants of abrin I and II to thymocyte at 37°C were 2.26 × 107 M−1 and 2.8 × 107 M−1 respectively. pH-induced structural alterations as studied by a parallel enhancement in 8-anilino-L-naphthalene sulfonate fluorescence revealed a high degree of qualitative similarity. These results taken with a nearly identical concentration-independent lag time (minimum lag of 41–42 min) indicated that the binding efficiencies and internalization efficiencies of these toxins are the same and that the observed difference in the dose-dependent lag time is causally related to the proposed processing event. The rates of reduction of inter-subunit disulfide bond, an obligatory step in the intoxication process, have been measured and compared under a variety of conditions. Intersubunit disulfide reduction of abrin I is fourfold faster than that of abrin II at pH 7.2. The rate of disulfide reduction in abrin I could be decreased 1 I-fold by adding lactose, compared to that without lactose. The observed differences in the efficiencies of A chains, the dose-dependent lag period, the modulating effect of lactose on the rates of disulfide reduction and similarity in binding properties make the variants a valuable tool to probe the processing events in toxin transport in detail.

Diversity in Functional Organization of Class I and Class II Biotin Protein Ligase

The cell envelope of Mycobacterium tuberculosis (M. tuberculosis) is composed of a variety of lipids including mycolic acids, sulpholipids, lipoarabinomannans, etc., which impart rigidity crucial for its survival and pathogenesis. Acyl CoA carboxylase (ACC) provides malonyl-CoA and methylmalonyl-CoA, committed precursors for fatty acid and essential for mycolic acid synthesis respectively. Biotin Protein Ligase (BPL/BirA) activates apo-biotin carboxyl carrier protein (BCCP) by biotinylating it to an active holo-BCCP. A minimal peptide (Schatz), an efficient substrate for Escherichia coli BirA, failed to serve as substrate for M. tuberculosis Biotin Protein Ligase (MtBPL). MtBPL specifically biotinylates homologous BCCP domain, MtBCCP87, but not EcBCCP87. This is a unique feature of MtBPL as EcBirA lacks such a stringent substrate specificity. This feature is also reflected in the lack of self/promiscuous biotinylation by MtBPL. The N-terminus/HTH domain of EcBirA has the selfbiotinable lysine residue that is inhibited in the presence of Schatz peptide, a peptide designed to act as a universal acceptor for EcBirA. This suggests that when biotin is limiting, EcBirA preferentially catalyzes, biotinylation of BCCP over selfbiotinylation. R118G mutant of EcBirA showed enhanced self and promiscuous biotinylation but its homologue, R69A MtBPL did not exhibit these properties. The catalytic domain of MtBPL was characterized further by limited proteolysis. Holo-MtBPL is protected from proteolysis by biotinyl-59 AMP, an intermediate of MtBPL catalyzed reaction. In contrast, apo-MtBPL is completely digested by trypsin within 20 min of co-incubation. Substrate selectivity and inability to promote self biotinylation are exquisite features of MtBPL and are a consequence of the unique molecular mechanism of an enzyme adapted for the high turnover of fatty acid biosynthesis.

Identification of Local Conformational Similarity in Structurally Variable Regions of Homologous Proteins Using Protein Blocks

Structure comparison tools can be used to align related protein structures to identify structurally conserved and variable regions and to infer functional and evolutionary relationships. While the conserved regions often superimpose well, the variable regions appear non superimposable. Differences in homologous protein structures are thought to be due to evolutionary plasticity to accommodate diverged sequences during evolution. One of the kinds of differences between 3-D structures of homologous proteins is rigid body displacement. A glaring example is not well superimposed equivalent regions of homologous proteins corresponding to a-helical conformation with different spatial orientations. In a rigid body superimposition, these regions would appear variable although they may contain local similarity. Also, due to high spatial deviation in the variable region, one-to-one correspondence at the residue level cannot be determined accurately. Another kind of difference is conformational variability and the most common example is topologically equivalent loops of two homologues but with different conformations. In the current study, we present a refined view of the ``structurally variable'' regions which may contain local similarity obscured in global alignment of homologous protein structures. As structural alphabet is able to describe local structures of proteins precisely through Protein Blocks approach, conformational similarity has been identified in a substantial number of `variable' regions in a large data set of protein structural alignments; optimal residue-residue equivalences could be achieved on the basis of Protein Blocks which led to improved local alignments. Also, through an example, we have demonstrated how the additional information on local backbone structures through protein blocks can aid in comparative modeling of a loop region. In addition, understanding on sequence-structure relationships can be enhanced through our approach. This has been illustrated through examples where the equivalent regions in homologous protein structures share sequence similarity to varied extent but do not preserve local structure.