A metabolomic approach to studying mechanisms of polymeric gene delivery to retinal pigment epithelial cells

Tiivistelmä ReferatAbstract Metabolomics is a rapidly growing research field that studies the response of biological systems to environmental factors, disease states and genetic modifications. It aims at measuring the complete set of endogenous metabolites, i.e. the metabolome, in a biological sample such as plasma or cells. Because metabolites are the intermediates and end products of biochemical reactions, metabolite compositions and metabolite levels in biological samples can provide a wealth of information on on-going processes in a living system. Due to the complexity of the metabolome, metabolomic analysis poses a challenge to analytical chemistry. Adequate sample preparation is critical to accurate and reproducible analysis, and the analytical techniques must have high resolution and sensitivity to allow detection of as many metabolites as possible. Furthermore, as the information contained in the metabolome is immense, the data set collected from metabolomic studies is very large. In order to extract the relevant information from such large data sets, efficient data processing and multivariate data analysis methods are needed. In the research presented in this thesis, metabolomics was used to study mechanisms of polymeric gene delivery to retinal pigment epithelial (RPE) cells. The aim of the study was to detect differences in metabolomic fingerprints between transfected cells and non-transfected controls, and thereafter to identify metabolites responsible for the discrimination. The plasmid pCMV-β was introduced into RPE cells using the vector polyethyleneimine (PEI). The samples were analyzed using high performance liquid chromatography (HPLC) and ultra performance liquid chromatography (UPLC) coupled to a triple quadrupole (QqQ) mass spectrometer (MS). The software MZmine was used for raw data processing and principal component analysis (PCA) was used in statistical data analysis. The results revealed differences in metabolomic fingerprints between transfected cells and non-transfected controls. However, reliable fingerprinting data could not be obtained because of low analysis repeatability. Therefore, no attempts were made to identify metabolites responsible for discrimination between sample groups. Repeatability and accuracy of analyses can be influenced by protocol optimization. However, in this study, optimization of analytical methods was hindered by the very small number of samples available for analysis. In conclusion, this study demonstrates that obtaining reliable fingerprinting data is technically demanding, and the protocols need to be thoroughly optimized in order to approach the goals of gaining information on mechanisms of gene delivery.

Afrikkalaisen sikaruton mahdollisia maahantuloreittejä

Afrikkalainen sikarutto on lakisääteisesti vastustettava helposti leviävä sikojen ja villisikojen virustauti, jolla on myös vakavia sosioekonomisia vaikutuksia. Akuutissa tautimuodossa eläimillä esiintyy syömättömyyttä, korkeaa kuumetta, verenpurkaumia iholla, veristä ulostetta sekä mahdollisesti ripulia. Kuolleisuus on lähes 100 % ja tauti johtaa kuolemaan 7-10 vrk tartunnan jälkeen. Tauti ei tartu ihmisiin. Tautia esiintyy sekä kesy- että villisioissa suurimmassa osassa Saharan eteläpuolista Afrikkaa ja Sardiniassa. Vuodesta 2007 lähtien tautia on esiintynyt Kaukasuksen alueella ja vuonna 2011 lähellä Suomen rajaa: Leningradin alueella ja Kuolan niemimaalla. Tässä riskiprofiilissa kartoitetaan reitit ja tapahtumasarjat, jotka voivat johtaa siihen, että afrikkalainen sikarutto tulee Suomeen ensimmäisen kerran. Näistä oleellisimmat ovat: maahantulo infektioalueella matkustaneiden ihmisten mukana, infektoituneen lihan tai lihatuotteen mukana, elävien kesysikojen ja sperman mukana, kontaminoituneiden eläinkuljetusajoneuvojen mukana, kansainvälisen liikenteen ruokajätteen mukana ja Suomeen vaeltavan infektoituneen villisian mukana. Tilatason tautisuojaus sekä tehokas ja kohdennettu tiedottaminen taudin vaaroista ovat avainasemassa, kun halutaan suojata Suomen sikaelinkeinoa afrikkalaiselta sikarutolta.

In vitro expression of belladonna mottle virus genome

In vitro translation of belladonna mottle virus BDMV(I) genomic RNA in a rabbit reticulocyte lysate system produced proteins of Mr 210,000, 150,000 and 78,000 which form the non-structural proteins. The coat protein, on the other hand, was expressed from a subgenomic RNA which was found to be encapsidated in the empty capsids forming the top component viral particles. The implications of subgenomic RNA encapsidation in viral replication and assembly are discussed.

Functions of human papillomavirus E5 oncogene in epithelial cells and in the onset of cervical cancer

Human papillomaviruses (HPVs) are the causal agents of cervical cancer, which is the second most common cancer among women worldwide. Cellular transformation and carcinogenesis depend on the activities of viral E5, E6 and E7 proteins. Alterations in cell-cell contacts and in communication between epithelial cells take place during cervical carcinogenesis, leading to changes in cell morphology, increased cell motility and finally invasion. The aim of this thesis was to study genome-wide effects of the HPV type 16 (HPV-16) E5 protein on the expression of host cell messenger RNAs (mRNAs) and microRNAs by applying microarray technology. The results showed that the HPV-16 E5 protein alters several cellular pathways involved in cellular adhesion, motility and proliferation as well as in the extracellular matrix. The E5 protein was observed to enhance wound healing of epithelial cell monolayers by increasing cell motility in vivo. HPV-16 E5-induced alterations in the expression of cellular microRNAs and their target genes seem to favour increased proliferation and tumorigenesis. E5 was also shown to affect the expression of adherens junction proteins in HaCaT epithelial keratinocytes. In addition, a study of a membrane cytoskeletal cross-linker protein, ezrin, revealed that when activated, it localizes to adherens junctions. The results suggest that ezrin distribution to forming adherens junctions is due to Rac1 activity in epithelial cells. These studies reveal for the first time the holistic effects of HPV-16 E5 protein in promoting precancerous events in epithelial cells. The results contribute to identifyinging novel markers for cervical precancerous stages and to predicting disease behaviour.

Substance misuse problems during pregnancy with special emphasis on buprenorphine

Maternal drug abuse during pregnancy endangers the future health and wellbeing of the infant and growing child. On the other hand, via maternal abstinence, these problems would never occur; so the problems would be totally preventable. Buprenorphine is widely used in opioid maintenance treatment as a substitute medication. In Finland, during 2000 s buprenorphine misuse has steadily increased. In 2009 almost one third of clientele of substance treatment units were in treatment because of buprenorphine dependence. At Helsinki Women s Clinic the first child with prenatal buprenorphine exposure was born in 2001. During 1992-2001 in the three capital area maternity hospitals (Women s clinic, Maternity hospital, Jorvi hospital) 524 women were followed at special antenatal clinics due to substance abuse problems. Three control women were drawn from birth register to each case woman and matched for parity and same place and date of the index birth. According to register data mortality rate was 38-fold higher among cases than controls within 6-15 years after index birth. Especially, the risk for violent or accidental death was increased. The women with substance misuse problems had also elevated risk for viral hepatitis and psychiatric morbidity. They were more often reimbursed for psychopharmaceuticals. Disability pensions and rehabilitation allowances were more often granted to cases than controls. In total 626 children were born from these pregnancies. According to register data 38% of these children were placed in out-of-home care as part of child protection services by the age of two years, and half of them by the age of 12 years, the median follow-up time was 5.8 years. The risk for out-of-home care was associated with factors identifiable during the pre- and perinatal period. In 2002-2005 67 pregnant women with buprenorphine dependence were followed up at the Helsinki University Hospital, Department of Obstetrics and Gynecology. Their pregnancies were uneventful. The prematurity rate was similar and there were no more major anomalies compared to the national statistics. The neonates were lighter compared to the national statistics. They were also born in good condition, with no perinatal hypoxia as defined by standard clinical parameters or certain biochemical markers in the cord blood: erythropoietin, S100 and cardiac troponin-t. Almost 80% of newborns developed neonatal abstinence syndrome (NAS) and two third of them needed morphine medication for it. Maternal smoking over ten cigarettes per day aggravated and benzodiazepine use attenuated NAS. An infant s highest urinary norbuprenorphine concentration during their first 3 days of life correlated with the duration of morphine treatment. The average length of infant s hospital stay was 25 days.

Physical mapping of the genomic DNA of the Oryctes rhinoceros baculovirus, KI

A non-occluded baculovirus, OBV-KI has been isolated from the insect pest, Oryctes rhinoceros. The viral genome is estimated to be 123 kb, with a G + C content of 43 mol% and no detectible methylated bases. A restriction map of the OBV-KI genome for BamHI, EcoRI, HindIII, PstI, SalI and XbaI has been constructed.

Structure-function relationships of icosahedral plant viruses

X-ray diffraction studies on single crystals of a few viruses have led to the elucidation of their three dimensional structure at near atomic resolution. Both the tertiary structure of the coat protein subunit and the quaternary morganization of the icosahedral capsid in these viruses are remarkably similar. These studies have led to a critical re-examination of the structural principles in the architecture of isometric viruses and suggestions of alternative mechanisms of assembly. Apart from their role in the assembly of the virus particle, the coat proteins of certian viruses have been shown to inhibit the replication of the cognate RNA leading to cross-protection. The coat protein amino acid sequence and the genomic sequence of several spherical plant RNA viruses have been determined in the last decade. Experimental data on the mechanisms of uncoating, gene expression and replication of several classes of viruses have also become available. The function of the non-structural proteins of some viruses have been determined. This rapid progress has provided a wealth of information on several key steps in the life cycle of RNA viruses. The function of the viral coat protein, capsid architecture, assembly and disassembly and replication of isometric RNA plant viruses are discussed in the light of this accumulated knowledge.

Glycoprotein Interactions in the Assembly of Hantaviruses

Hantaviruses are one of the five genera of the vector-borne virus family Bunyaviridae. While other members of the family are transmitted via arthropods, hantaviruses are carried and transmitted by rodents and insectivores. Occasional transmission to humans occurs via inhalation of aerosolized rodent excreta. When transmitted to man hantaviruses cause hemorrhagic fever with renal syndrome (HFRS, in Eurasia, mortality ~10%) and hantavirus cardiopulmonary syndrome (HCPS, in the Americas, mortality ~40%). The single-stranded, negative-sense RNA genome of hantaviruses is in segments S, M and L that respectively encode for nucleocapsid (N), glycoproteins Gn and Gc, and RNA-dependent RNA-polymerase (RdRp or L protein). The genome segments, encapsidated by N protein to form ribonucleoprotein (RNP), are enclosed inside a lipid envelope decorated by spikes formed of Gn and Gc. The focus of this study was to understand the mechanisms and interactions through which the virion is formed and maintained. We observed that when extracted from virions both Gn and Gc favor homo- over hetero-oligomerization. The minimal glycoprotein complexes extracted from virion by detergent were observed, by using ultracentrifugation and gel filtration, to be tetrameric Gn and homodimeric Gc. These results led us to suggest a model where tetrameric Gn complexes are interconnected through homodimeric Gc units to form the grid-like surface architecture described for hantaviruses. This model was found to correlate with the three-dimensional (3D) reconstruction of virion surface created using cryo-electron tomography (cryo-ET). The 3D-density map showed the spike complex formed of Gn and Gc to be 10 nm high and to display a four-fold symmetry with dimensions of 15 nm times 15 nm. This unique square-shaped complex on a roughly round virion creates a hitch for the assembly, since a sphere cannot be broken into rectangles. Thus additional interactions are likely required for the virion assembly. In cryo-ET we observed that the RNP makes occasional contacts to the viral membrane, suggesting an interaction between the spike and RNP. We were able to demonstrate this interaction using various techniques, and showed that both Gn and Gc contribute to the interaction. This led us to suggest that in addition to the interactions between Gn and Gc, also the interaction between spike and RNP is required for assembly. We found galectin-3 binding protein (referred to as 90K) to co-purify with the virions and showed an interaction between 90K and the virion. Analysis of plasma samples taken from patients hospitalized for Puumala virus infection showed increased concentrations of 90K in the acute phase and the increased 90K level was found to correlate with several parameters that reflect the severity of acute HFRS. The results of these studies confirmed, but also challenged some of the dogmas on the structure and assembly of hantaviruses. We confirmed that Gn and RNP do interact, as long assumed. On the other hand we demonstrated that the glycoproteins Gn and Gc exist as homo-oligomers or appear in large hetero-oligomeric complexes, rather than form primarily heterodimers as was previously assumed. This work provided new insight into the structure and assembly of hantaviruses.

Hantavirus infection: Insights into entry, assembly and pathogenesis

Hantaviruses (family Bunyaviridae, genus Hantavirus) are enveloped viruses incorporating a segmented, negative-sense RNA genome. Each hantavirus is carried by its specific host, either a rodent or an insectivore (shrew), in which the infection is asymptomatic and persistent. In humans, hantaviruses cause Hemorrhagic fever with renal syndrome (HFRS) in Eurasia and Hantavirus cardiopulmonary syndrome (HCPS) in the Americas. In Finland, Puumala virus (genus Hantavirus) is the causative agent of NE, a mild form of HFRS. The HFRS-type diseases are often associated with renal failure and proteinuria that might be mechanistically explained by infected kidney tubular cell degeneration in patients. Previously, it has been shown that non-pathogenic hantavirus, Tula virus (TULV), could cause programmed cell death, apoptosis, in cell cultures. This suggested that the infected kidney tubular degeneration could be caused directly by virus replication. In the first paper of this thesis the molecular mechanisms involved in TULV-induced apoptosis was further elucidated. A virus replication-dependent down-regulation of ERK1/2, concomitantly with the induced apoptosis, was identified. In addition, this phenomenon was not restricted to TULV or to non-pathogenic hantaviruses in general since also a pathogenic hantavirus, Seoul virus, could inhibit ERK1/2 activity. Hantaviruses consist of membrane-spanning glycoproteins Gn and Gc, RNA-dependent RNA polymerase (L protein) and nucleocapsid protein N, which encapsidates the viral genome, and thus forms the ribonucleoprotein (RNP). Interaction between the cytoplasmic tails of viral glycoproteins and RNP is assumed to be the only means how viral genetic material is incorporated into infectious virions. In the second paper of this thesis, it was shown by immunoprecipitation that viral glycoproteins and RNP interact in the purified virions. It was further shown that peptides derived from the cytoplasmic tails (CTs) of both Gn and Gc could bind RNP and recombinant N protein. In the fourth paper the cytoplamic tail of Gn but not Gc was shown to interact with genomic RNA. This interaction was probably rather unspecific since binding of Gn-CT with unrelated RNA and even single-stranded DNA were also observed. However, since the RNP consists of both N protein and N protein-encapsidated genomic RNA, it is possible that the viral genome plays a role in packaging of RNPs into virions. On the other hand, the nucleic acid-binding activity of Gn may have importance in the synthesis of viral RNA. Binding sites of Gn-CT with N protein or nucleic acids were also determined by peptide arrays, and they were largely found to overlap. The Gn-CT of hantaviruses contain a conserved zinc finger (ZF) domain with an unknown function. Some viruses need ZFs in entry or post-entry steps of the viral life cycle. Cysteine residues are required for the folding of ZFs by coordinating zinc-ions, and alkylation of these residues can affect virus infectivity. In the third paper, it was shown that purified hantavirions could be inactivated by treatment with cysteine-alkylating reagents, especially N-ethyl maleimide. However, the effect could not be pin-pointed to the ZF of Gn-CT since also other viral proteins reacted with maleimides, and it was, therefore, impossible to exclude the possibility that other cysteines besides those that were essential in the formation of ZF are required for hantavirus infectivity.

Two forms of Pf1 Inovirus: X-ray diffraction studies on a structural phase transition and a calculated libration normal mode of the asymmetric unit

Inovirus is a helical array of alpha-helical protein asymmetric units surrounding a DNA core. X-ray fibre diffraction studies show that the Pf1 species of Inovirus can undergo a reversible temperature-induced transition between two similar structural forms having slightly different virion helix parameters. Molecular models of the two forms show no evidence for altered interactions between the protein and either the solvent or the viral DNA; but there are significant differences in the shape and orientation of the protein asymmetric unit, related to the changes in the virion parameters. Normal modes involving libration of whole asymmetric units are in a frequency range with appreciable entropy of libration, and the structural transition may be related to changes in libration.

Two forms of PF1 INOVIRUS: X-ray diffraction studies on a structural phase transition and a calculated libration normal mode of the assymetric unit

Inovirus is a helical array of agr-helical protein asymmetric units surrounding a DNA core. X-ray fibre diffraction studies show that the Pf1 species of Inovirus can undergo a reversible temperature-induced transition between two similar structural forms having slightly different virion helix parameters. Molecular models of the two forms show no evidence for altered interactions between the protein and either the solvent or the viral DNA; but there are significant differences in the shape and orientation of the protein asymmetric unit, related to the changes in the virion parameters. Normal modes involving libration of whole asymmetric units are in a frequency range with appreciable entropy of libration, and the structural transition may be related to changes in libration.

Taking Multiple Infections of Cells and Recombination into Account Leads to Small Within-Host Effective-Population-Size Estimates of HIV-1

Whether HIV-1 evolution in infected individuals is dominated by deterministic or stochastic effects remains unclear because current estimates of the effective population size of HIV-1 in vivo, N-e, are widely varying. Models assuming HIV-1 evolution to be neutral estimate N-e similar to 10(2)-10(4), smaller than the inverse mutation rate of HIV-1 (similar to 10(5)), implying the predominance of stochastic forces. In contrast, a model that includes selection estimates N-e>10(5), suggesting that deterministic forces would hold sway. The consequent uncertainty in the nature of HIV-1 evolution compromises our ability to describe disease progression and outcomes of therapy. We perform detailed bit-string simulations of viral evolution that consider large genome lengths and incorporate the key evolutionary processes underlying the genomic diversification of HIV-1 in infected individuals, namely, mutation, multiple infections of cells, recombination, selection, and epistatic interactions between multiple loci. Our simulations describe quantitatively the evolution of HIV-1 diversity and divergence in patients. From comparisons of our simulations with patient data, we estimate N-e similar to 10(3)-10(4), implying predominantly stochastic evolution. Interestingly, we find that N-e and the viral generation time are correlated with the disease progression time, presenting a route to a priori prediction of disease progression in patients. Further, we show that the previous estimate of N-e>10(5) reduces as the frequencies of multiple infections of cells and recombination assumed increase. Our simulations with N-e similar to 10(3)-10(4) may be employed to estimate markers of disease progression and outcomes of therapy that depend on the evolution of viral diversity and divergence.

Interaction of Sesbania Mosaic Virus Movement Protein with VPg and P10: Implication to Specificity of Genome Recognition

Sesbania mosaic virus (SeMV) is a single strand positive-sense RNA plant virus that belongs to the genus Sobemovirus. The mechanism of cell-to-cell movement in sobemoviruses has not been well studied. With a view to identify the viral encoded ancillary proteins of SeMV that may assist in cell-to-cell movement of the virus, all the proteins encoded by SeMV genome were cloned into yeast Matchmaker system 3 and interaction studies were performed. Two proteins namely, viral protein genome linked (VPg) and a 10-kDa protein (P10) c v gft encoded by OFR 2a, were identified as possible interacting partners in addition to the viral coat protein (CP). Further characterization of these interactions revealed that the movement protein (MP) recognizes cognate RNA through interaction with VPg, which is covalently linked to the 59 end of the RNA. Analysis of the deletion mutants delineated the domains of MP involved in the interaction with VPg and P10. This study implicates for the first time that VPg might play an important role in specific recognition of viral genome by MP in SeMV and shed light on the possible role of P10 in the viral movement.

A Shapley Value-Based Approach to Discover Influential Nodes in Social Networks

Our study concerns an important current problem, that of diffusion of information in social networks. This problem has received significant attention from the Internet research community in the recent times, driven by many potential applications such as viral marketing and sales promotions. In this paper, we focus on the target set selection problem, which involves discovering a small subset of influential players in a given social network, to perform a certain task of information diffusion. The target set selection problem manifests in two forms: 1) top-k nodes problem and 2) lambda-coverage problem. In the top-k nodes problem, we are required to find a set of k key nodes that would maximize the number of nodes being influenced in the network. The lambda-coverage problem is concerned with finding a set of k key nodes having minimal size that can influence a given percentage lambda of the nodes in the entire network. We propose a new way of solving these problems using the concept of Shapley value which is a well known solution concept in cooperative game theory. Our approach leads to algorithms which we call the ShaPley value-based Influential Nodes (SPINs) algorithms for solving the top-k nodes problem and the lambda-coverage problem. We compare the performance of the proposed SPIN algorithms with well known algorithms in the literature. Through extensive experimentation on four synthetically generated random graphs and six real-world data sets (Celegans, Jazz, NIPS coauthorship data set, Netscience data set, High-Energy Physics data set, and Political Books data set), we show that the proposed SPIN approach is more powerful and computationally efficient. Note to Practitioners-In recent times, social networks have received a high level of attention due to their proven ability in improving the performance of web search, recommendations in collaborative filtering systems, spreading a technology in the market using viral marketing techniques, etc. It is well known that the interpersonal relationships (or ties or links) between individuals cause change or improvement in the social system because the decisions made by individuals are influenced heavily by the behavior of their neighbors. An interesting and key problem in social networks is to discover the most influential nodes in the social network which can influence other nodes in the social network in a strong and deep way. This problem is called the target set selection problem and has two variants: 1) the top-k nodes problem, where we are required to identify a set of k influential nodes that maximize the number of nodes being influenced in the network and 2) the lambda-coverage problem which involves finding a set of influential nodes having minimum size that can influence a given percentage lambda of the nodes in the entire network. There are many existing algorithms in the literature for solving these problems. In this paper, we propose a new algorithm which is based on a novel interpretation of information diffusion in a social network as a cooperative game. Using this analogy, we develop an algorithm based on the Shapley value of the underlying cooperative game. The proposed algorithm outperforms the existing algorithms in terms of generality or computational complexity or both. Our results are validated through extensive experimentation on both synthetically generated and real-world data sets.

Ribavirin-Induced Anemia in Hepatitis C Virus Patients Undergoing Combination Therapy

The current standard of care for hepatitis C virus (HCV) infection - combination therapy with pegylated interferon and ribavirin - elicits sustained responses in only similar to 50% of the patients treated. No alternatives exist for patients who do not respond to combination therapy. Addition of ribavirin substantially improves response rates to interferon and lowers relapse rates following the cessation of therapy, suggesting that increasing ribavirin exposure may further improve treatment response. A key limitation, however, is the toxic side-effect of ribavirin, hemolytic anemia, which often necessitates a reduction of ribavirin dosage and compromises treatment response. Maximizing treatment response thus requires striking a balance between the antiviral and hemolytic activities of ribavirin. Current models of viral kinetics describe the enhancement of treatment response due to ribavirin. Ribavirin-induced anemia, however, remains poorly understood and precludes rational optimization of combination therapy. Here, we develop a new mathematical model of the population dynamics of erythrocytes that quantitatively describes ribavirin-induced anemia in HCV patients. Based on the assumption that ribavirin accumulation decreases erythrocyte lifespan in a dose-dependent manner, model predictions capture several independent experimental observations of the accumulation of ribavirin in erythrocytes and the resulting decline of hemoglobin in HCV patients undergoing combination therapy, estimate the reduced erythrocyte lifespan during therapy, and describe inter-patient variations in the severity of ribavirin-induced anemia. Further, model predictions estimate the threshold ribavirin exposure beyond which anemia becomes intolerable and suggest guidelines for the usage of growth hormones, such as erythropoietin, that stimulate erythrocyte production and avert the reduction of ribavirin dosage, thereby improving treatment response. Our model thus facilitates, in conjunction with models of viral kinetics, the rational identification of treatment protocols that maximize treatment response while curtailing side effects.