Dramatic down-regulation of oxidoreductases in human hepatocellular carcinoma hepG2 cells: proteomics and gene ontology unveiling new frontiers in cancer enzymology

Background

Oxidoreductases are enzymes that catalyze many redox reactions in normal and neoplastic cells. Their actions include catalysis of the transformation of free, neutral oxygen gas into oxygen free radicals, superoxide, hydroperoxide, singlet oxygen and hydrogen peroxide. These activated forms of oxygen contribute to oxidative stress that modifies lipids, proteins, DNA and carbohydrates. On the other hand, oxidoreductases constitute one of the most important free radical scavenger systems typified by catalase, superoxide dismutase and glutathione peroxidase. In this work, proteomics, Gene Ontology mapping and Directed Acyclic Graphs (DAG) are employed to detect and quantify differential oxidoreductase enzyme expressions between HepG2 cells and normal human liver tissues.

Results

For the set of bioinformatics calculations whose BLAST searches are performed using the BLAST program BLASTP 2.2.13 [Nov-27-2005], DAG of the Gene Ontology's Molecular Function annotations show that oxidoreductase activity parent node of the liver proteome contains 331 annotated protein sequences, 7 child nodes and an annotation score of 188.9, whereas that of HepG2 cells has 188 annotated protein sequences, 3 child nodes and an annotation score of only 91.9. Overwhelming preponderance of oxidoreductases in the liver is additionally supported by the isomerase DAGs: nearly all the reactions described in the normal liver isomerase DAG are oxidoreductase isomerization reactions, whereas only one of the three child nodes in the HepG2 isomerase DAG is oxidoreductase. Upon normalization of the annotation scores to the parent Molecular Function nodes, oxidoreductases are down-regulated in HepG2 cells by 58%. Similarly, for the set of bioinformatics calculations whose BLAST searches are carried out using BLASTP 2.2.15 [Oct-15-2006], oxidoreductases are down-regulated in HepG2 cells by 56%.

Conclusion

Proteomics and Gene Ontology reveal, for the first time, differential enzyme activities between HepG2 cells and normal human liver tissues, which may be a promising new prognostic marker of Hepatocellular carcinoma. Two independent sets of bioinformatics calculations that employ two BLAST program versions, and searched different databases, arrived at essentially the same conclusion: oxidoreductases are down-regulated in HepG2 cells by approximately 57%, when compared to normal human liver tissues. Down-regulation of oxidoreductases in hepatoma is additionally supported by Gene Ontology analysis of isomerises.     

Sample prep for proteomics of breast cancer: proteomics and gene ontology reveal dramatic differences in protein solubilization preferences of radioimmunoprecipitation assay and urea lysis buffers

Background

An important step in the proteomics of solid tumors, including breast cancer, consists of efficiently extracting most of proteins in the tumor specimen. For this purpose, Radio-Immunoprecipitation Assay (RIPA) buffer is widely employed. RIPA buffer's rapid and highly efficient cell lysis and good solubilization of a wide range of proteins is further augmented by its compatibility with protease and phosphatase inhibitors, ability to minimize non-specific protein binding leading to a lower background in immunoprecipitation, and its suitability for protein quantitation.

Results

In this work, the insoluble matter left after RIPA buffer extraction of proteins from breast tumors are subjected to another extraction step, using a urea-based buffer. It is shown that RIPA and urea lysis buffers fractionate breast tissue proteins primarily on the basis of molecular weights. The average molecular weight of proteins that dissolve exclusively in urea buffer is up to 60% higher than in RIPA. Gene Ontology (GO) and Directed Acyclic Graphs (DAG) are used to map the collective biological and biophysical attributes of the RIPA and urea proteomes. The Cellular Component and Molecular Function annotations reveal protein solubilization preferences of the buffers, especially the compartmentalization and functional distributions. It is shown that nearly all extracellular matrix proteins (ECM) in the breast tumors and matched normal tissues are found, nearly exclusively, in the urea fraction, while they are mostly insoluble in RIPA buffer. Additionally, it is demonstrated that cytoskeletal and extracellular region proteins are more soluble in urea than in RIPA, whereas for nuclear, cytoplasmic and mitochondrial proteins, RIPA buffer is preferred. Extracellular matrix proteins are highly implicated in cancer, including their proteinase-mediated degradation and remodelling, tumor development, progression, adhesion and metastasis. Thus, if they are not efficiently extracted by RIPA buffer, important information may be missed in cancer research.

Conclusion

For proteomics of solid tumors, a two-step extraction process is recommended. First, proteins in the tumor specimen should be extracted with RIPA buffer. Second, the RIPA-insoluble material should be extracted with the urea-based buffer employed in this work.     

A comparison of alternative methods to compute conditional genotype probabilities for genetic evaluation with finite locus models

An increased availability of genotypes at marker loci has prompted the development of models that include the effect of individual genes. Selection based on these models is known as marker-assisted selection (MAS). MAS is known to be efficient especially for traits that have low heritability and non-additive gene action. BLUP methodology under non-additive gene action is not feasible for large inbred or crossbred pedigrees. It is easy to incorporate non-additive gene action in a finite locus model. Under such a model, the unobservable genotypic values can be predicted using the conditional mean of the genotypic values given the data. To compute this conditional mean, conditional genotype probabilities must be computed. In this study these probabilities were computed using iterative peeling, and three Markov chain Monte Carlo (MCMC) methods – scalar Gibbs, blocking Gibbs, and a sampler that combines the Elston Stewart algorithm with iterative peeling (ESIP). The performance of these four methods was assessed using simulated data. For pedigrees with loops, iterative peeling fails to provide accurate genotype probability estimates for some pedigree members. Also, computing time is exponentially related to the number of loci in the model. For MCMC methods, a linear relationship can be maintained by sampling genotypes one locus at a time. Out of the three MCMC methods considered, ESIP, performed the best while scalar Gibbs performed the worst.     

Effect of a mutation on the structure and dynamics of an alpha-helical antifreeze protein in water and ice

One strategy of psychrophilic organisms to survive subzero temperature is to produce antifreeze protein (AFPs), which inhibit the growth of macromolecular ice. To better understand the binding mechanism, the structure and dynamics of several AFPs have been studied by nuclear magnetic resonance (NMR) and X-ray crystallography. The results have shown that different organisms can use diverse structures (alpha-helix, beta-helix, or different globular folds) to achieve the same function. A number of studies have focused on understanding the relationship between the alpha-helical structure of fish type I AFP and its function as an inhibitor of ice growth. The results have not explained whether the 90% activity loss caused by the conservative mutation of two threonines to serines (Thr13Ser/Thr24Ser) is attributable to a change in protein structure in solution or in ice. We examine here the structure and dynamics of the winter flounder type I AFP and the mutant Thr13Ser/Thr24Ser in both solution and solid states using a wide range of NMR approaches. Both proteins remain fully alpha-helical at all temperatures and in ice, demonstrating that the activity change must therefore not be attributable to changes in the protein fold or dynamics but differences in surface properties.     

Flying lemurs – The 'flying tree shrews'? Molecular cytogenetic evidence for a Scandentia-Dermoptera sister clade

Background  

Flying lemurs or Colugos (order Dermoptera) represent an ancient mammalian lineage that contains only two extant species. Although molecular evidence strongly supports that the orders Dermoptera, Scandentia, Lagomorpha, Rodentia and Primates form a superordinal clade called Supraprimates (or Euarchontoglires), the phylogenetic placement of Dermoptera within Supraprimates remains ambiguous.     

Ligand-dependent differences in estrogen receptor beta-interacting proteins identified in lung adenocarcinoma cells corresponds to estrogenic responses

Background

A recent epidemiological study demonstrated a reduced risk of lung cancer mortality in breast cancer patients using antiestrogens. These and other data implicate a role for estrogens in lung cancer, particularly nonsmall cell lung cancer (NSCLC). Approximately 61% of human NSCLC tumors express nuclear estrogen receptor β (ERβ); however, the role of ERβ and estrogens in NSCLC is likely to be multifactorial. Here we tested the hypothesis that proteins interacting with ERβ in human lung adenocarcinoma cells that respond proliferatively to estradiol (E₂) are distinct from those in non-E₂-responsive cells.

Methods

FLAG affinity purification of FLAG-ERβ-interacting proteins was used to isolate ERβ-interacting proteins in whole cell extracts from E₂ proliferative H1793 and non-E₂-proliferative A549 lung adenocarcinoma cell lines. Following trypsin digestion, proteins were identified using liquid chromatography electrospray ionization tandem mass spectrometry (LC-MS/MS). Proteomic data were analyzed using Ingenuity Pathway Analysis. Select results were confirmed by coimmunoprecipitation.

Results

LC-MS/MS identified 27 non-redundant ERβ-interacting proteins. ERβ-interacting proteins included hsp70, hsp60, vimentin, histones and calmodulin. Ingenuity Pathway Analysis of the ERβ-interacting proteins revealed differences in molecular and functional networks between H1793 and A549 lung adenocarcinoma cells. Coimmunoprecipitation experiments in these and other lung adenocarcinoma cells confirmed that ERβ and EGFR interact in a gender-dependent manner and in response to E₂ or EGF. BRCA1 interacted with ERβ in A549 cell lines and in human lung adenocarcinoma tumors, but not normal lung tissue.

Conclusion

Our results identify specific differences in ERβ-interacting proteins in lung adenocarcinoma cells corresponding to ligand-dependent differences in estrogenic responses.

     

Human CC chemokine I-309, structural consequences of the additional disulfide bond

I-309 is a member of the CC subclass of chemokines and is one of only three human chemokines known to contain an additional, third disulfide bond. The three-dimensional solution structure of I-309 was determined by (1)H nuclear magnetic resonance spectroscopy and dynamic simulated annealing. The structure of I-309, which remains monomeric at high concentrations, was determined on the basis of 978 experimental restraints. The N-terminal region of I-309 was disordered, as has been previously observed for the CC chemokine eotaxin but not others such as MCP-1 and RANTES. This was followed in I-309 by a well-ordered region between residues 13 and 69 that consisted of a 3(10)-helix, a triple-stranded antiparallel beta-sheet, and finally a C-terminal alpha-helix. Root-mean-square deviations of 0.61 and 1.16 were observed for the backbone and heavy atoms, respectively. A comparison of I-309 to eotaxin and HCC-2 revealed a significant structural change in the C-terminal region of the protein. The alpha-helix normally present in chemokines was terminated early and was followed by a short section of extended strand. These changes were a direct result of the additional disulfide bond present in this protein. An examination of the I-309 structure will aid in an understanding of the specificity of this protein with its receptor, CCR8.     

BAC libraries construction from the ancestral diploid genomes of the allotetraploid cultivated peanut

Background

Cultivated peanut, Arachis hypogaea is an allotetraploid of recent origin, with an AABB genome. In common with many other polyploids, it seems that a severe genetic bottle-neck was imposed at the species origin, via hybridisation of two wild species and spontaneous chromosome duplication. Therefore, the study of the genome of peanut is hampered both by the crop's low genetic diversity and its polyploidy. In contrast to cultivated peanut, most wild Arachis species are diploid with high genetic diversity. The study of diploid Arachis genomes is therefore attractive, both to simplify the construction of genetic and physical maps, and for the isolation and characterization of wild alleles. The most probable wild ancestors of cultivated peanut are A. duranensis and A. ipaënsis with genome types AA and BB respectively.

Results

We constructed and characterized two large-insert libraries in Bacterial Artificial Chromosome (BAC) vector, one for each of the diploid ancestral species. The libraries (AA and BB) are respectively c. 7.4 and c. 5.3 genome equivalents with low organelle contamination and average insert sizes of 110 and 100 kb. Both libraries were used for the isolation of clones containing genetically mapped legume anchor markers (single copy genes), and resistance gene analogues.

Conclusion

These diploid BAC libraries are important tools for the isolation of wild alleles conferring resistances to biotic stresses, comparisons of orthologous regions of the AA and BB genomes with each other and with other legume species, and will facilitate the construction of a physical map.

     

Differential responsiveness to immunoablative therapy in refractory rheumatoid arthritis is associated with level and avidity of anti-cyclic citrullinated protein autoantibodies: a case study

In order to identify pathogenic correlates of refractory rheumatoid arthritis (RA), antibodies against anti-cyclic citrullinated protein (ACPAs) were investigated in RA patients in whom the dysregulated immune system had been ablated by high-dose chemotherapy (HDC) and autologous haematopoietic stem cell transplantation (HSCT). Six patients with refractory RA were extensively characterized in terms of levels of total immunoglobulins, RA-specific autoantibodies (ACPAs and rheumatoid factor) and antibodies against rubella, tetanus toxoid (TT) and phosphorylcholine before and after HDC plus HSCT. Additionally, the avidity of ACPAs was measured before and after treatment and compared with the avidity of TT antibodies following repeated immunizations. Synovial biopsies were obtained by arthroscopy before HDC plus HSCT, and analyzed by immunohistochemistry. In the three patients with clinically long-lasting responses to HDC plus HSCT (median 423 days), significant reductions in ACPA-IgG levels after therapy were observed (median level dropped from 215 to 34 arbitrary units/ml; P = 0.05). In contrast, stable ACPA-IgG levels were observed in three patients who relapsed shortly after HDC plus HSCT (median of 67 days). Clinical responders had ACPA-IgG of lower avidity (r = 0.75; P = 0.08) and higher degree of inflammation histologically (r = 0.73; P = 0.09). Relapse (after 38 to 530 days) in all patients was preceded by rising levels of low avidity ACPA-IgG (after 30 to 388 days), in contrast to the stable titres of high avidity TT antibodies. In conclusion, humoral autoimmune responses were differentially modulated by immunoablative therapy in patients with synovial inflammation and low avidity ACPA-IgG autoantibodies as compared with patients with high levels of high avidity ACPA-IgG. The distinct clinical disease course after immunoablative therapy based on levels and avidity of ACPA-IgG indicates that refractory RA is not a single disease entity.