Regulation of Osteoblast Differentiation and Gene Expression by Phosphodiesterases and Bioactive Peptides

Bone is a mineralized tissue that enables multiple mechanical and metabolic functions to be carried out in the skeleton. Bone contains distinct cell types: osteoblasts (bone-forming cells), osteocytes (mature osteoblast that embedded in mineralized bone matrix) and the osteoclasts (bone-resorbing cells). Remodelling of bone begins early in foetal life, and once the skeleton is fully formed in young adults, almost all of the metabolic activity is in this form. Bone is constantly destroyed or resorbed by osteoclasts and then replaced by osteoblasts. Many bone diseases, i.e. osteoporosis, also known as bone loss, typically reflect an imbalance in skeletal turnover. The cyclic adenosine monophosphate (cAMP) and the cyclic guanosine monophosphate (cGMP) are second messengers involved in a variety of cellular responses to such extracellular agents as hormones and neurotransmitters. In the hormonal regulation of bone metabolism, i.e. via parathyroid hormone (PTH), parathyroid hormone-related peptide (PTHrp) and prostaglandin E2 signal via cAMP. cAMP and cGMP are formed by adenylate and guanylate cyclases and are degraded by phosphodiesterases (PDEs). PDEs determine the amplitudes of cyclic nucleotide-mediated hormonal responses and modulate the duration of the signal. The activities of the PDEs are regulated by multiple inputs from other signalling systems and are crucial points of cross-talk between the pathways. Food-derived bioactive peptides are reported to express a variety of functions in vivo. The angiotensin-converting enzymes (ACEs) are involved in the regulation of the specific maturation or degradation of a number of mammalian bioactive peptides. The bioactive peptides offer also a nutriceutical and a nutrigenomic aspect to bone cell biology. The aim of this study was to investigate the influence of PDEs and bioactive peptides on the activation and the differentiation of human osteoblast cells. The profile of PDEs in human osteoblast-like cells and the effect of glucocorticoids on the function of cAMP PDEs, were investigated at the mRNA and enzyme levels. The effects of PDEs on bone formation and osteoblast gene expression were determined with chemical inhibitors and siRNAs (short interfering RNAs). The influence of bioactive peptides on osteoblast gene expression and proliferation was studied at the mRNA and cellular levels. This work provides information on how PDEs are involved in the function and the differentiation of osteoblasts. The findings illustrate that gene-specific silencing with an RNA interference (RNAi) method is useful in inhibiting, the gene expression of specific PDEs and further, PDE7 inhibition upregulates several osteogenic genes and increases bALP activity and mineralization in human mesenchymal stem cells-derived osteoblasts. PDEs appear to be involved in a mechanism by which glucocorticoids affect cAMP signaling. This may provide a potential route in the formation of glucocorticoid-induced bone loss, involving the down-regulation of cAMP-PDE. PDEs may play an important role in the regulation of osteoblastic differentiation. Isoleucine-proline-proline (IPP), a bioactive peptide, possesses the potential to increase osteoblast proliferation, differentiation and signalling.

Neuronal K-Cl Cotransporter : Transcriptional Mechanisms of KCC2 Gene Regulation

K-Cl cotransporter 2 (KCC2) maintains a low intracellular Cl concentration required for fast hyperpolarizing responses of neurons to classical inhibitory neurotransmitters γ-aminobutyric acid (GABA) and glycine. Decreased Cl extrusion observed in genetically modified KCC2-deficient mice leads to depolarizing GABA responses, impaired brain inhibition, and as a consequence to epileptic seizures. Identification of mechanisms regulating activity of the SLC12A5 gene, which encodes the KCC2 cotransporter, in normal and pathological conditions is, thus, of extreme importance. Multiple reports have previously elucidated in details a spatio-temporal pattern of KCC2 expression. Among the characteristic features are an exclusive neuronal specificity, a dramatic upregulation during embryonic and early postnatal development, and a significant downregulation by neuronal trauma. Numerous studies confirmed these expressional features, however transcriptional mechanisms predetermining the SLC12A5 gene behaviour are still unknown. The aim of the presented thesis is to recognize such transcriptional mechanisms and, on their basis, to create a transcriptional model that would explain the established SLC12A5 gene behaviour. Up to recently, only one KCC2 transcript has been thought to exist. A particular novelty of the presented work is the identification of two SLC12A5 gene promoters (SLC12A5-1a and SLC12A5-1b) that produce at least two KCC2 isoforms (KCC2a and KCC2b) differing by their N-terminal parts. Even though a functional 86Rb+ assay reveals no significant difference between transport activities of the isoforms, consensus sites for several protein kinases, found in KCC2a but not in KCC2b, imply a distinct kinetic regulation. As a logical continuation, the current work presents a detailed analysis of the KCC2a and KCC2b expression patterns. This analysis shows an exclusively neuron-specific pattern and similar expression levels for both isoforms during embryonic and neonatal development in rodents. During subsequent postnatal development, the KCC2b expression dramatically increases, while KCC2a expression, depending on central nervous system (CNS) area, either remains at the same level or moderately decreases. In an attempt to explain both the neuronal specificity and the distinct expressional kinetics of the KCC2a and KCC2b isoforms during postnatal development, the corresponding SLC12A5-1a and SLC12A5-1b promoters have been subjected to a comprehensive bioinformatical analysis. Binding sites of several transcription factors (TFs), conserved in the mammalian SLC12A5 gene orthologs, have been identified that might shed light on the observed behaviour of the SLC12A5 gene. Possible roles of these TFs in the regulating of the SLC12A5 gene expression have been elucidated in subsequent experiments and are discussed in the current thesis.

Molecular epidemiology of human rhinoviruses

The first part of this work investigates the molecular epidemiology of a human enterovirus (HEV), echovirus 30 (E-30). This project is part of a series of studies performed in our research team analyzing the molecular epidemiology of HEV-B viruses. A total of 129 virus strains had been isolated in different parts of Europe. The sequence analysis was performed in three different genomic regions: 420 nucleotides (nt) in the VP4/VP2 capsid protein coding region, the entire VP1 capsid protein coding gene of 876 nt, and 150 nt in the VP1/2A junction region. The analysis revealed a succession of dominant sublineages within a major genotype. The temporally earlier genotypes had been replaced by a genetically homogenous lineage that has been circulating in Europe since the late 1970s. The same genotype was found by other research groups in North America and Australia. Globally, other cocirculating genetic lineages also exist. The prevalence of a dominant genotype makes E-30 different from other previously studied HEVs, such as polioviruses and coxsackieviruses B4 and B5, for which several coexisting genetic lineages have been reported. The second part of this work deals with molecular epidemiology of human rhinoviruses (HRVs). A total of 61 field isolates were studied in the 420-nt stretch in the capsid coding region of VP4/VP2. The isolates were collected from children under two years of age in Tampere, Finland. Sequences from the clinical isolates clustered in the two previously known phylogenetic clades. Seasonal clustering was found. Also, several distinct serotype-like clusters were found to co-circulate during the same epidemic season. Reappearance of a cluster after disappearing for a season was observed. The molecular epidemiology of the analyzed strains turned out to be complex, and we decided to continue our studies of HRV. Only five previously published complete genome sequences of HRV prototype strains were available for analysis. Therefore, all designated HRV prototype strains (n=102) were sequenced in the VP4/VP2 region, and the possibility of genetic typing of HRV was evaluated. Seventy-six of the 102 prototype strains clustered in HRV genetic group A (HRV-A) and 25 in group B (HRV-B). Serotype 87 clustered separately from other HRVs with HEV species D. The field strains of HRV represented as many as 19 different genotypes, as judged with an approximate demarcation of a 20% nt difference in the VP4/VP2 region. The interserotypic differences of HRV were generally similar to those reported between different HEV serotypes (i.e. about 20%), but smaller differences, less than 10%, were also observed. Because some HRV serotypes are genetically so closely related, we suggest that the genetic typing be performed using the criterion "the closest prototype strain". This study is the first systematic genetic characterization of all known HRV prototype strains, providing a further taxonomic proposal for classification of HRV. We proposed to divide the genus Human rhinoviruses into HRV-A and HRV-B. The final part of the work comprises a phylogenetic analysis of a subset (48) of HRV prototype strains and field isolates (12) in the nonstructural part of the genome coding for the RNA-dependent RNA polymerase (3D). The proposed division of the HRV strains in the species HRV-A and HRV-B was also supported by 3D region. HRV-B clustered closer to HEV species B, C, and also to polioviruses than to HRV-A. Intraspecies variation within both HRV-A and HRV-B was greater in the 3D coding region than in the VP4/VP2 coding region, in contrast to HEV. Moreover, the diversity of HRV in 3D exceeded that of HEV. One group of HRV-A, designated HRV-A', formed a separate cluster outside other HRV-A in the 3D region. It formed a cluster also in the capsid region, but located within HRV-A. This may reflect a different evolutionary history of distinct genomic regions among HRV-A. Furthermore, the tree topology within HRV-A in the 3D region differed from that in the VP4/VP2, suggesting possible recombination events in the evolution of the strains. No conflicting phylogenies were observed in any of the 12 field isolates. Possible recombination was further studied using the Similarity and Bootscanning analyses of the complete genome sequences of HRV available in public databases. Evidence for recombination among HRV-A was found, as HRV2 and HRV39 showed higher similarity in the nonstructural part of the genome. Whether HRV2 and HRV39 strains - and perhaps also some other HRV-A strains not yet completely sequenced - are recombinants remains to be determined.

Positional cloning and pathway analysis of the asthma susceptibility gene, NPSR1

In the present study, we identified a novel asthma susceptibility gene, NPSR1 (neuropeptide S receptor 1) on chromosome 7p14.3 by the positional cloning strategy. An earlier significant linkage mapping result among Finnish Kainuu asthma families was confirmed in two independent cohorts: in asthma families from Quebec, Canada and in allergy families from North Karelia, Finland. The linkage region was narrowed down to a 133-kb segment by a hierarchial genotyping method. The observed 77-kb haplotype block showed 7 haplotypes and a similar risk and nonrisk pattern in all three populations studied. All seven haplotypes occur in all three populations at frequences > 2%. Significant elevated relative risks were detected for elevated total IgE (immunoglobulin E) or asthma. Risk effects of the gene variants varied from 1.4 to 2.5. NPSR1 belongs to the G protein-coupled receptor (GPCR) family with a topology of seven transmembrane domains. NPSR1 has 9 exons, with the two main transcripts, A and B, encoding proteins of 371 and 377 amino acids, respectively. We detected a low but ubiquitous expression level of NPSR1-B in various tissues and endogenous cell lines while NPSR1-A has a more restricted expression pattern. Both isoforms were expressed in the lung epithelium. We observed aberrant expression levels of NPSR1-B in smooth muscle in asthmatic bronchi as compared to healthy. In an experimental mouse model, the induced lung inflammation resulted in elevated Npsr1 levels. Furthermore, we demonstrated that the activation of NPSR1 with its endogenous agonist, neuropeptide S (NPS), resulted in a significant inhibition of the growth of NPSR1-A overexpressing stable cell lines (NPSR1-A cells). To determine which target genes were regulated by the NPS-NPSR1 pathway, NPSR1-A cells were stimulated with NPS, and differentially expressed genes were identified using the Affymetrix HGU133Plus2 GeneChip. A total of 104 genes were found significantly up-regulated and 42 down-regulated 6 h after NPS administration. The up-regulated genes included many neuronal genes and some putative susceptibility genes for respiratory disorders. By Gene Ontology enrichment analysis, the biological process terms, cell proliferation, morphogenesis and immune response were among the most altered. The expression of four up-regulated genes, matrix metallopeptidase 10 (MMP10), INHBA (activin A), interleukin 8 (IL8) and EPH receptor A2 (EPHA2), were verified and confirmed by quantitative reverse-transcriptase-PCR. In conclusion, we identified a novel asthma susceptibility gene, NPSR1, on chromosome 7p14.3. NPS-NPSR1 represents a novel pathway that regulates cell proliferation and immune responses, and thus may have functional relevance in the pathogenesis of asthma.

Characterization of LRRTM and NGR gene families: Expression and functions

Some leucine-rich repeat (LRR) -containing membrane proteins are known regulators of neuronal growth and synapse formation. In this work I characterize two gene families encoding neuronal LRR membrane proteins, namely the LRRTM (leucine-rich repeat, transmembrane neuronal) and NGR (Nogo-66 receptor) families. I studied LRRTM and NGR family member's mRNA tissue distribution by RT-PCR and by in situ hybridization. Subcellular localization of LRRTM1 protein was studied in neurons and in non-neuronal cells. I discovered that LRRTM and NGR family mRNAs are predominantly expressed in the nervous system, and that each gene possesses a specific expression pattern. I also established that LRRTM and NGR family mRNAs are expressed by neurons, and not by glial cells. Within neurons, LRRTM1 protein is not transported to the plasma membrane; rather it localizes to endoplasmic reticulum. Nogo-A (RTN4), MAG, and OMgp are myelin-associated proteins that bind to NgR1 to limit axonal regeneration after central nervous system injury. To better understand the functions of NgR2 and NgR3, and to explore the possible redundancy in the signaling of myelin inhibitors of neurite growth, I mapped the interactions between NgR family and the known and candidate NgR1 ligands. I identified high-affinity interactions between RTN2-66, RTN3-66 and NgR1. I also demonstrate that Rtn3 mRNA is expressed in the same glial cell population of mouse spinal cord white matter as Nogo-A mRNA, and thus it could have a role in myelin inhibition of axonal growth. To understand how NgR1 interacts with multiple structurally divergent ligands, I aimed first to map in more detail the nature of Nogo-A:NgR1 interactions, and then to systematically map the binding sites of multiple myelin ligands in NgR1 by using a library of NgR1 expression constructs encoding proteins with one or multiple surface residues mutated to alanine. My analysis of the Nogo-A:NgR1 -interactions revealed a novel interaction site between the proteins, suggesting a trivalent Nogo-A:NgR1-interaction. Our analysis also defined a central binding region on the concave side of NgR1's LRR domain that is required for the binding of all known ligands, and a surrounding region critical for binding MAG and OMgp. To better understand the biological role of LRRTMs, I generated Lrrtm1 and Lrrtm3 knock out mice. I show here that reporter genes expressed from the targeted loci can be used for maping the neuronal connections of Lrrtm1 and Lrrtm3 expressing neurons in finer detail. With regard to LRRTM1's role in humans, we found a strong association between a 70 kb-spanning haplotype in the proposed promoter region of LRRTM1 gene and two possibly related phenotypes: left-handedness and schizophrenia. Interestingly, the responsible haplotype was linked to phenotypic variability only when paternally inherited. In summary, I identified two families of neuronal receptor-like proteins, and mapped their expression and certain protein-protein interactions. The identification of a central binding region in NgR1 shared by multiple ligands may facilitate the design and development of small molecule therapeutics blocking binding of all NgR1 ligands. Additionally, the genetic association data suggests that allelic variation upstream of LRRTM1 may play a role in the development of left-right brain asymmetry in humans. Lrrtm1 and Lrrtm3 knock out mice developed as a part of this study will likely be useful for schizophrenia and Alzheimer s disease research.

Approaches for improving the safety and efficacy of adenoviral gene therapy

Cancer is a devastating disease with poor prognosis and no curative treatment, when widely metastatic. Conventional therapies, such as chemotherapy and radiotherapy, have efficacy but are not curative and systemic toxicity can be considerable. Almost all cancers are caused due to changes in the genetic material of the transformed cells. Cancer gene therapy has emerged as a new treatment option, and past decades brought new insights in developing new therapeutic drugs for curing cancer. Oncolytic viruses constitute a novel therapeutic approach given their capacity to replicate in and kill specifically tumor cells as well as reaching tumor distant metastasis. Adenoviral gene therapy has been suggested to cause liver toxicity. This study shows that new developed adenoviruses, in particular Ad5/19p-HIT, can be redirected towards kidney while adenovirus uptake by liver is minimal. Moreover, low liver transduction resulted in a favorable tumor to liver ratio of virus load. Further, we established a new immunocompetent animal model Syrian hamsters. Wild type adenovirus 5 was found to replicate in Hap-T1 hamster tumors and normal tissues. There are no antiviral drugs available to inhibit adenovirus replication. In our study, chlorpromazine and cidofovir efficiently abrogated virus replication in vitro and showed significant reduction in vivo in tumors and liver. Once safety concerns were addressed together with the new given antiviral treatment options, we further improved oncolytic adenoviruses for better tumor penetration, local amplification and host system modulation. Further, we created Ad5/3-9HIF-Δ24-VEGFR-1-Ig, oncolytic adenovirus for improved infectivity and antiangiogenic effect for treatment of renal cancer. This virus exhibited increased anti-tumor effect and specific replication in kidney cancer cells. The key player for good efficacy of oncolytic virotherapy is the host immune response. Thus, we engineered a triple targeted adenovirus Ad5/3-hTERT-E1A-hCD40L, which would lead to tumor elimination due to tumor-specific oncolysis and apoptosis together with an anti-tumor immune response prompted by the immunomodulatory molecule. In conclusion, the results presented in this thesis constitute advances in our understanding of oncolytic virotherapy by successful tumor targeting, antiviral treatment options as a safety switch in case of replication associated side-effects, and modulation of the host immune system towards tumor elimination.  

Ion channel gene variants predisposing to severe human ventricular arrhythmias

Congenital long QT syndrome (LQTS) with an estimated prevalence of 1:2000-1:10 000 manifests with prolonged QT interval on electrocardiogram and risk for ventricular arrhythmias and sudden death. Several ion channel genes and hundreds of mutations in these genes have been identified to underlie the disorder. In Finland, four LQTS founder mutations of potassium channel genes account for up to 40-70% of genetic spectrum of LQTS. Acquired LQTS has similar clinical manifestations, but often arises from usage of QT-prolonging medication or electrolyte disturbances. A prolonged QT interval is associated with increased morbidity and mortality not only in clinical LQTS but also in patients with ischemic heart disease and in the general population. The principal aim of this study was to estimate the actual prevalence of LQTS founder mutations in Finland and to calculate their effect on QT interval in the Finnish background population. Using a large population-based sample of over 6000 Finnish individuals from the Health 2000 Survey, we identified LQTS founder mutations KCNQ1 G589D (n=8), KCNQ1 IVS7-2A>G (n=1), KCNH2 L552S (n=2), and KCNH2 R176W (n=16) in 27 study participants. This resulted in a weighted prevalence estimate of 0.4% for LQTS in Finland. Using a linear regression model, the founder mutations resulted in a 22- to 50-ms prolongation of the age-, sex-, and heart rate-adjusted QT interval. Collectively, these data suggest that one of 250 individuals in Finland may be genetically predisposed to ventricular arrhythmias arising from the four LQTS founder mutations. A KCNE1 D85N minor allele with a frequency of 1.4% was associated with a 10-ms prolongation in adjusted QT interval and could thus identify individuals at increased risk of ventricular arrhythmias at the population level. In addition, the previously reported associations of KCNH2 K897T, KCNH2 rs3807375, and NOS1AP rs2880058 with QT interval duration were confirmed in the present study. In a separate study, LQTS founder mutations were identified in a subgroup of acquired LQTS, providing further evidence that congenital LQTS gene mutations may underlie acquired LQTS. Catecholaminergic polymorphic ventricular tachycardia (CPVT) is characterized by exercise-induced ventricular arrhythmias in a structurally normal heart and results from defects in the cardiac Ca2+ signaling proteins, mainly ryanodine receptor type 2 (RyR2). In a patient population of typical CPVT, RyR2 mutations were identifiable in 25% (4/16) of patients, implying that noncoding variants or other genes are involved in CPVT pathogenesis. A 1.1 kb RyR2 exon 3 deletion was identified in two patients independently, suggesting that this region may provide a new target for RyR2-related molecular genetic studies. Two novel RyR2 mutations showing a gain-of-function defect in vitro were identified in three victims of sudden cardiac death. Extended pedigree analyses revealed some surviving mutation carriers with mild structural abnormalities of the heart and resting ventricular arrhythmias suggesting that not all RyR2 mutations lead to a typical CPVT phenotype, underscoring the relevance of tailored risk stratification of a RyR2 mutation carrier.

Monte Carlo Simulations of Molecular Clusters in Nucleation

A better understanding of the limiting step in a first order phase transition, the nucleation process, is of major importance to a variety of scientific fields ranging from atmospheric sciences to nanotechnology and even to cosmology. This is due to the fact that in most phase transitions the new phase is separated from the mother phase by a free energy barrier. This barrier is crossed in a process called nucleation. Nowadays it is considered that a significant fraction of all atmospheric particles is produced by vapor-to liquid nucleation. In atmospheric sciences, as well as in other scientific fields, the theoretical treatment of nucleation is mostly based on a theory known as the Classical Nucleation Theory. However, the Classical Nucleation Theory is known to have only a limited success in predicting the rate at which vapor-to-liquid nucleation takes place at given conditions. This thesis studies the unary homogeneous vapor-to-liquid nucleation from a statistical mechanics viewpoint. We apply Monte Carlo simulations of molecular clusters to calculate the free energy barrier separating the vapor and liquid phases and compare our results against the laboratory measurements and Classical Nucleation Theory predictions. According to our results, the work of adding a monomer to a cluster in equilibrium vapour is accurately described by the liquid drop model applied by the Classical Nucleation Theory, once the clusters are larger than some threshold size. The threshold cluster sizes contain only a few or some tens of molecules depending on the interaction potential and temperature. However, the error made in modeling the smallest of clusters as liquid drops results in an erroneous absolute value for the cluster work of formation throughout the size range, as predicted by the McGraw-Laaksonen scaling law. By calculating correction factors to Classical Nucleation Theory predictions for the nucleation barriers of argon and water, we show that the corrected predictions produce nucleation rates that are in good comparison with experiments. For the smallest clusters, the deviation between the simulation results and the liquid drop values are accurately modelled by the low order virial coefficients at modest temperatures and vapour densities, or in other words, in the validity range of the non-interacting cluster theory by Frenkel, Band and Bilj. Our results do not indicate a need for a size dependent replacement free energy correction. The results also indicate that Classical Nucleation Theory predicts the size of the critical cluster correctly. We also presents a new method for the calculation of the equilibrium vapour density, surface tension size dependence and planar surface tension directly from cluster simulations. We also show how the size dependence of the cluster surface tension in equimolar surface is a function of virial coefficients, a result confirmed by our cluster simulations.