Statistical Methods for Conservation and Alignment Quality in Proteins

Construction of multiple sequence alignments is a fundamental task in Bioinformatics. Multiple sequence alignments are used as a prerequisite in many Bioinformatics methods, and subsequently the quality of such methods can be critically dependent on the quality of the alignment. However, automatic construction of a multiple sequence alignment for a set of remotely related sequences does not always provide biologically relevant alignments.Therefore, there is a need for an objective approach for evaluating the quality of automatically aligned sequences. The profile hidden Markov model is a powerful approach in comparative genomics. In the profile hidden Markov model, the symbol probabilities are estimated at each conserved alignment position. This can increase the dimension of parameter space and cause an overfitting problem. These two research problems are both related to conservation. We have developed statistical measures for quantifying the conservation of multiple sequence alignments. Two types of methods are considered, those identifying conserved residues in an alignment position, and those calculating positional conservation scores. The positional conservation score was exploited in a statistical prediction model for assessing the quality of multiple sequence alignments. The residue conservation score was used as part of the emission probability estimation method proposed for profile hidden Markov models. The results of the predicted alignment quality score highly correlated with the correct alignment quality scores, indicating that our method is reliable for assessing the quality of any multiple sequence alignment. The comparison of the emission probability estimation method with the maximum likelihood method showed that the number of estimated parameters in the model was dramatically decreased, while the same level of accuracy was maintained. To conclude, we have shown that conservation can be successfully used in the statistical model for alignment quality assessment and in the estimation of emission probabilities in the profile hidden Markov models.

Functional studies on bacterial nucleotide-regulated inorganic pyrophosphatases

CBS domains are ~60 amino acid tandemly repeated regulatory modules forming a widely distributed domain superfamily. Found in thousands of proteins from all kingdoms of life, CBS domains have adopted a variety of functions during evolution, one of which is regulation of enzyme activity through binding of adenylate-containing compounds in a hydrophobic cavity. Mutations in human CBS domain-containing proteins cause hereditary diseases. Inorganic pyrophosphatases (PPases) are ubiquitous enzymes, which pull pyrophosphate (PPi) producing reactions forward by hydrolyzing PPi into phosphate. Of the two nonhomologous soluble PPases, dimeric family II PPases, belonging to the DHH family of phosphoesterases, require a transition metal and magnesium for maximal activity. A quarter of the almost 500 family II PPases, found in bacteria and archaea, contain a 120-250 amino acid N-terminal insertion, comprised of two CBS domains separated in sequence by a DRTGG domain. These enzymes are thus named CBS-PPases. The function of the DRTGG domain in proteins is unknown. The aim of this PhD thesis was to elucidate the structural and functional differences of CBS-PPases in comparison to family II PPases lacking the regulatory insert. To this end, we expressed, purified and characterized the CBS-PPases from Clostridium perfringens (cpCBS-PPase) and Moorella thermoacetica (mtCBS-PPase), the latter lacking a DRTGG domain. Both enzymes are homodimers in solution and display maximal activity against PPi in the presence of Co2+ and Mg2+. Uniquely, the DRTGG domain was found to enable tripolyphosphate hydrolysis at rates similar to that of PPi. Additionally, we found that AMP and ADP inhibit, while ATP and AP4A activate CBSPPases, thus enabling regulation in response to changes in cellular energy status. We then observed substrate- and nucleotide-induced conformational transitions in mtCBS-PPase and found that the enzyme exists in two differentially active conformations, interconverted through substrate binding and resulting in a 2.5-fold enzyme activation. AMP binding was shown to produce an alternate conformation, which is reached through a different pathway than the substrate-induced conformation. We solved the structure of the regulatory insert from cpCBS-PPase in complex with AMP and AP4A and proposed that conformational changes in the loops connecting the catalytic and regulatory domains enable activity regulation. We examined the effects of mutations in the CBS domains of mtCBS-PPase on catalytic activity, as well as, nucleotide binding and inhibition.

Critical flux and fouling in ultrafiltration of proteins

In this thesis different parameters influencing critical flux in protein ultrafiltration and membrane foul-ing were studied. Short reviews of proteins, cross-flow ultrafiltration, flux decline and criticalflux and the basic theory of Partial Least Square analysis (PLS) are given at the beginning. The experiments were mainly performed using dilute solutions of globular proteins, commercial polymeric membranes and laboratory scale apparatuses. Fouling was studied by flux, streaming potential and FTIR-ATR measurements. Critical flux was evaluated by different kinds of stepwise procedures and by both con-stant pressure and constant flux methods. The critical flux was affected by transmembrane pressure, flow velocity, protein concentration, mem-brane hydrophobicity and protein and membrane charges. Generally, the lowest critical fluxes were obtained at the isoelectric points of the protein and the highest in the presence of electrostatic repulsion between the membrane surface and the protein molecules. In the laminar flow regime the critical flux increased with flow velocity, but not any more above this region. An increase in concentration de-creased the critical flux. Hydrophobic membranes showed fouling in all charge conditionsand, furthermore, especially at the beginning of the experiment even at very low transmembrane pressures. Fouling of these membranes was thought to be due to protein adsorption by hydrophobic interactions. The hydrophilic membranes used suffered more from reversible fouling and concentration polarisation than from irreversible foul-ing. They became fouled at higher transmembrane pressures becauseof pore blocking. In this thesis some new aspects on critical flux are presented that are important for ultrafiltration and fractionation of proteins.

Neurofilament proteins in the postnatal rat hippocampus. Developmental expression and changes in experimental epilepsy

Neurofilament proteins (NFs) are the major components of the intermediate filaments of the neuronal cytoskeleton. The three different NF proteins; the low (NF-L), medium (NF-M),and dendrites.NF proteins play an important role in neuronal development, and plasticity,and seem to contribute to the pathophysiology of several diseases. However, the detailed expression patterns of NF proteins in the course of postnatal aturation, and in response to seizures in the rat have remained unknown. In this work, I have studied the developmental expression and cellular distribution of the three NF proteins in the rat hippocampus during the postnatal development. The reactivity of NF proteins in response to kainic acid (KA)-induced status epilepticus (SE)was studied in the hippocampus of 9-day-old rats, and using in vitro organotypic hippocampal slices cultures prepared from P6-7 rats. The results showed that NF-L and NF-M proteins are expressed already at the postnatal day 1, while the expression of NF-H mainly occurred during the second postnatal week. The immunoreactivity of NF proteins varied depending on the cell type and sub-cellular location in the hippocampus. In adult rats, KA-induced SE typically results in severe and permanent NF degradation. However, in our P9 rats KA-induced SE resulted in a transient increase in the expression of NF proteins during the first few hours but not degradation. No neuronal death or mossy fiber sprouting was observed at any time after SE. The in vitro studies with OHCs, which mimick the in vivo developing models where a local injection of KA is applied(e.g. intrahippocampal), indicated that NF proteins were rapidly degraded in response to KA treatment, this effect being effectively inhibited by the treatment with the AMPA receptor antagonist CNQX, and calpain inhibitor MDL-28170. These compounds also significantly ameliorated the KA-induced region-specific neuronal damage. The NMDA receptor antagonist and the L-type Ca2+ channel blocker did not have any significant effect. In conclusion, the results indicate that the developmental expression of NF in the rat hippocampus is differentially regulated and targeted in the different hippocampal cell types during the postnatal development. Furthermore, despite SE, the mechanisms leading to NF degradation and neuronal death are not activated in P9 rats unlike in adults. The reason for this remains unknown. The results in organotypic hippocampal cultures confirm the validity of this in vitro model to study development processes, and to perform pharmacological studies. The results also suggest that calpain proteases as interesting pharmacological targets to reduce neuronal damage after acute excitotoxic insults.

Novel Proteins Involved in Dynamics of The Photosynthetic Apparatus

During the past few years, a considerable number of research articles have been published relating to the structure and function of the major photosynthetic protein complexes, photosystem (PS) I, PSII, cytochrome (Cyt) b6f, and adenosine triphosphate (ATP) synthase. Sequencing of the Arabidopsis thaliana (Arabidopsis) genome together with several high-quality proteomics studies has, however, revealed that the thylakoid membrane network of plant chloroplasts still contains a number of functionally unknown proteins. These proteins may have a role as auxiliary proteins guiding the assembly, maintenance, and turnover of the thylakoid protein complexes, or they may be as yet unknown subunits of the photosynthetic complexes. Novel subunits are most likely to be found in the NAD(P)H dehydrogenase (NDH) complex, the structure and function of which have remained obscure in the absence of detailed crystallographic data, thus making this thylakoid protein complex a particularly interesting target of investigation. In this thesis, several novel thylakoid-associated proteins were identified by proteomics-based methods. The major goal of characterization of the stroma thylakoid associated polysome-nascent chain complexes was to determine the proteins that guide the dynamic life cycle of PSII. In addition, a large protein complex of ≥ 1,000 kDa, residing in the stroma thylakoid, was characterized in greater depth and it was found to be a supercomplex composed of the PSI and NDH complexes. A set of newly identified proteins from Arabidopsis thylakoids was subjected to detailed characterization using the reverse genetics approach and extensive biochemical and biophysical analysis. The role of the novel proteins, either as auxiliary proteins or subunits of the photosynthetic protein complexes, was revealed. Two novel thylakoid lumen proteins, TLP18.3 and AtCYP38, function as auxiliary proteins assisting specific steps of the assembly/repair of PSII. The role of the 10-kDa thylakoid lumen protein PsbR is related to the optimization of oxygen evolution of PSII by assisting the assembly of the PsbP protein. Two integral thylakoid membrane proteins, NDH45 and NDH48, are novel subunits of the chloroplast NDH complex. Finally, the thylakoid lumen immunophilin AtCYP20-2 is suggested to interact with the NDH complex, instead of PSII as was hypothesized earlier.

Novel CO2 Regulated Proteins in Synechocystis PCC 6803

The large biodiversity of cyanobacteria together with the increasing genomics and proteomics metadata provide novel information for finding new commercially valuable metabolites. With the advent of global warming, there is growing interest in the processes that results in efficient CO2 capture through the use of photosynthetic microorganisms such as cyanobacteria. This requires a detailed knowledge of how cyanobacteria respond to the ambient CO2. My study was aimed at understanding the changes in the protein profile of the model organism, Synechocystis PCC 6803 towards the varying CO₂ level. In order to achieve this goal I have employed modern proteomics tools such as iTRAQ and DIGE, recombinant DNA techniques to construct different mutants in cyanobacteria and biophysical methods to study the photosynthetic properties. The proteomics study revealed several novel proteins, apart from the well characterized proteins involved in carbon concentrating mechanisms (CCMs), that were upregulated upon shift of the cells from high CO₂ concentration (3%) to that in air level (0.039%). The unknown proteins, Slr0006 and flavodiiron proteins (FDPs) Sll0217-Flv4 and Sll0219-Flv2, were selected for further characterization. Although slr0006 was substantially upregulated under C_i limiting conditions, inactivation of the gene did not result in any visual phenotype under various environmental conditions indicating that this protein is not essential for cell survival. However, quantitative proteomics showed the induction of novel plasmid and chromosome encoded proteins in deltaslr0006 under air level CO₂ conditions. The expression of the slr0006 gene was found to be strictly dependent on active photosynthetic electron transfer. Slr0006 contains conserved dsRNA binding domain that belongs to the Sua5/YrdC/YciO protein family. Structural modelling of Slr0006 showed an alpha/beta twisted open-sheet structure and a positively charged cavity, indicating a possible binding site for RNA. The 3D model and the co-localization of Slr0006 with ribosomal subunits suggest that it might play a role in translation or ribosome biogenesis. On the other hand, deletions in the sll0217-sll218- sll0219 operon resulted in enhanced photodamage of PSII and distorted energy transfer from phycobilisome (PBS) to PSII, suggesting a dynamic photoprotection role of the operon. Constructed homology models also suggest efficient electron transfer in heterodimeric Flv2/Flv4, apparently involved in PSII photoprotection. Both Slr0006 and FDPs exhibited several common features, including negative regulation by NdhR and ambiguous cellular localization when subjected to different concentrations of divalent ions. This strong association with the membranes remained undisturbed even in the presence of detergent or high salt. My finding brings ample information on three novel proteins and their functions towards carbon limitation. Nevertheless, many pathways and related proteins remain unexplored. The comprehensive understanding of the acclimation processes in cyanobacteria towards varying environmental CO₂ levels will help to uncover adaptive mechanisms in other organisms, including higher plants.

Formin proteins in normal tissues and cancer

The actin cytoskeleton is a dynamic structure that determines cell shape. Actin turnover is mandatory for migration in normal and malignant cells. In epithelial cancers invasion is frequently accompanied by epithelial to mesenchymal transition (EMT). In EMT, cancer cells acquire a migratory phenotype through transcriptional reprogramming. EMT requires substantial re-organization of actin. During the past decade, new actin regulating proteins have been discovered. Among these are members of the formin family. To study formin expression in tissues and cells, antibodies for detection of formin proteins FMNL1 (Formin-like protein 1), FMNL2 (Formin-like protein 2) and FHOD1 (Formin homology 2 domain containing protein 1) were used. The expression of formins was characterized in normal tissues and selected cancers using immunohistochemistry. The functional roles of formins were studied in cancer cell lines. We found that FMNL2 is widely expressed. It is a filopodial component in cultured melanoma cells. In clinical melanoma, FMNL2 expression has prognostic significance. FHOD1 is a formin expressed in mesenchymal cell types. FHOD1 expression is increased in oral squamous cell carcinoma (SCC) EMT. Importantly, FHOD1 participates in invasion of cultured oral SCC cells. FMNL1 expression is low in normal epithelia, but high in leukocytes and smooth muscle cells. Expression of FMNL1 can be found in carcinoma; we detected FMNL1 expressing cells in basal type of breast cancer. Our results indicate that formins are differentially expressed in normal tissues and that their expression may shift in cancer. Functionally FMNL2 and FHOD1 participate in processes related to cancer progression. Studying formins is increasingly important since they are potential drug targets.

Life without Phd? – Phd Oxygen Sensor Proteins in Hypoxic Carcinoma Cell Survival

The microenvironment within the tumor plays a central role in cellular signaling. Rapidly proliferating cancer cells need building blocks for structures as well as nutrients and oxygen for energy production. In normal tissue, the vasculature effectively transports oxygen, nutrient and waste products, and maintains physiological pH. Within a tumor however, the vasculature is rarely sufficient for the needs of tumor cells. This causes the tumor to suffer from lack of oxygen (hypoxia) and nutrients as well as acidification, as the glycolytic end product lactate is accumulated. Cancer cells harbor mutations enabling survival in the rough microenvironment. One of the best characterized mutations is the inactivation of the von Hippel-Lindau protein (pVHL) in clear cell renal cell carcinoma (ccRCC). Inactivation causes constitutive activation of hypoxia-inducible factor HIF which is an important survival factor regulating glycolysis, neovascularization and apoptosis. HIFs are normally regulated by HIF prolyl hydroxylases (PHDs), which in the presence of oxygen target HIF α-subunit to ubiquitination by pVHL and degradation by proteasomes. In my thesis work, I studied the role of PHDs in the survival of carcinoma cells in hypoxia. My work revealed an essential role of PHD1 and PHD3 in cell cycle regulation through two cyclin-dependent kinase inhibitors (CKIs) p21 and p27. Depletion of PHD1 or PHD3 caused a cell cycle arrest and subjected the carcinoma cells to stress and impaired the survival.