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.

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.

Characterisation and source identification of pollution episodes caused by long-range transported aerosols

Aerosol particles can cause detrimental environmental and health effects. The particles and their precursor gases are emitted from various anthropogenic and natural sources. It is important to know the origin and properties of aerosols to efficiently reduce their harmful effects. The diameter of aerosol particles (Dp) varies between ~0.001 and ~100 μm. Fine particles (PM2.5: Dp < 2.5 μm) are especially interesting because they are the most harmful and can be transported over long distances. The aim of this thesis is to study the impact on air quality by pollution episodes of long-range transported aerosols affecting the composition of the boundary-layer atmosphere in remote and relatively unpolluted regions of the world. The sources and physicochemical properties of aerosols were investigated in detail, based on various measurements (1) in southern Finland during selected long-range transport (LRT) pollution episodes and unpolluted periods and (2) over the Atlantic Ocean between Europe and Antarctica during a voyage. Furthermore, the frequency of LRT pollution episodes of fine particles in southern Finland was investigated over a period of 8 years, using long-term air quality monitoring data. In southern Finland, the annual mean PM2.5 mass concentrations were low but LRT caused high peaks of daily mean concentrations every year. At an urban background site in Helsinki, the updated WHO guideline value (24-h PM2.5 mean 25 μg/m3) was exceeded during 1-7 LRT episodes each year during 1999-2006. The daily mean concentrations varied between 25 and 49 μg/m3 during the episodes, which was 3-6 times higher than the mean concentration in the long term. The in-depth studies of selected LRT episodes in southern Finland revealed that biomass burning in agricultural fields and wildfires, occurring mainly in Eastern Europe, deteriorated air quality on a continental scale. The strongest LRT episodes of fine particles resulted from open biomass-burning fires but the emissions from other anthropogenic sources in Eastern Europe also caused significant LRT episodes. Particle mass and number concentrations increased strongly in the accumulation mode (Dp ~ 0.09-1 μm) during the LRT episodes. However, the concentrations of smaller particles (Dp < 0.09 μm) remained low or even decreased due to the uptake of vapours and molecular clusters by LRT particles. The chemical analysis of individual particles showed that the proportions of several anthropogenic particle types increased (e.g. tar balls, metal oxides/hydroxides, spherical silicate fly ash particles and various calcium-rich particles) in southern Finland during an LRT episode, when aerosols originated from the polluted regions of Eastern Europe and some open biomass-burning smoke was also brought in by LRT. During unpolluted periods when air masses arrived from the north, the proportions of marine aerosols increased. In unpolluted rural regions of southern Finland, both accumulation mode particles and small-sized (Dp ~ 1-3 μm) coarse mode particles originated mostly from LRT. However, the composition of particles was totally different in these size fractions. In both size fractions, strong internal mixing of chemical components was typical for LRT particles. Thus, the aging of particles has significant impacts on their chemical, hygroscopic and optical properties, which can largely alter the environmental and health effects of LRT aerosols. Over the Atlantic Ocean, the individual particle composition of small-sized (Dp ~ 1-3 μm) coarse mode particles was affected by continental aerosol plumes to distances of at least 100-1000 km from the coast (e.g. pollutants from industrialized Europe, desert dust from the Sahara and biomass-burning aerosols near the Gulf of Guinea). The rate of chloride depletion from sea-salt particles was high near the coasts of Europe and Africa when air masses arrived from polluted continental regions. Thus, the LRT of continental aerosols had significant impacts on the composition of the marine boundary-layer atmosphere and seawater. In conclusion, integration of the results obtained using different measurement techniques captured the large spatial and temporal variability of aerosols as observed at terrestrial and marine sites, and assisted in establishing the causal link between land-bound emissions, LRT and air quality.

Observations on the first steps of atmospheric particle formation and growth

Atmospheric aerosol particles have a significant impact on air quality, human health and global climate. The climatic effects of secondary aerosol are currently among the largest uncertainties limiting the scientific understanding of future and past climate changes. To better estimate the climatic importance of secondary aerosol particles, detailed information on atmospheric particle formation mechanisms and the vapours forming the aerosol is required. In this thesis we studied these issues by applying novel instrumentation in a boreal forest to obtain direct information on the very first steps of atmospheric nucleation and particle growth. Additionally, we used detailed laboratory experiments and process modelling to determine condensational growth properties, such as saturation vapour pressures, of dicarboxylic acids, which are organic acids often found in atmospheric samples. Based on our studies, we came to four main conclusions: 1) In the boreal forest region, both sulphurous compounds and organics are needed for secondary particle formation, the previous contributing mainly to particle formation and latter to growth; 2) A persistent pool of molecular clusters, both neutral and charged, is present and participates in atmospheric nucleation processes in boreal forests; 3) Neutral particle formation seems to dominate over ion-mediated mechanisms, at least in the boreal forest boundary layer; 4) The subcooled liquid phase saturation vapour pressures of C3-C9 dicarboxylic acids are of the order of 1e-5 1e-3 Pa at atmospheric temperatures, indicating that a mixed pre-existing particulate phase is required for their condensation in atmospheric conditions. The work presented in this thesis gives tools to better quantify the aerosol source provided by secondary aerosol formation. The results are particularly useful when estimating, for instance, anthropogenic versus biogenic influences and the fractions of secondary aerosol formation explained by neutral or ion-mediated nucleation mechanisms, at least in environments where the average particle formation rates are of the order of some tens of particles per cubic centimeter or lower. However, as the factors driving secondary particle formation are likely to vary depending on the environment, measurements on atmospheric nucleation and particle growth are needed from around the world to be able to better describe the secondary particle formation, and assess its climatic effects on a global scale.

Quantum Chemical Studies on Tropospheric Nucleation Mechanisms Involving Sulfuric Acid

Nucleation is the first step of the process by which gas molecules in the atmosphere condense to form liquid or solid particles. Despite the importance of atmospheric new-particle formation for both climate and health-related issues, little information exists on its precise molecular-level mechanisms. In this thesis, potential nucleation mechanisms involving sulfuric acid together with either water and ammonia or reactive biogenic molecules are studied using quantum chemical methods. Quantum chemistry calculations are based on the numerical solution of Schrödinger's equation for a system of atoms and electrons subject to various sets of approximations, the precise details of which give rise to a large number of model chemistries. A comparison of several different model chemistries indicates that the computational method must be chosen with care if accurate results for sulfuric acid - water - ammonia clusters are desired. Specifically, binding energies are incorrectly predicted by some popular density functionals, and vibrational anharmonicity must be accounted for if quantitatively reliable formation free energies are desired. The calculations reported in this thesis show that a combination of different high-level energy corrections and advanced thermochemical analysis can quantitatively replicate experimental results concerning the hydration of sulfuric acid. The role of ammonia in sulfuric acid - water nucleation was revealed by a series of calculations on molecular clusters of increasing size with respect to all three co-ordinates; sulfuric acid, water and ammonia. As indicated by experimental measurements, ammonia significantly assists the growth of clusters in the sulfuric acid - co-ordinate. The calculations presented in this thesis predict that in atmospheric conditions, this effect becomes important as the number of acid molecules increases from two to three. On the other hand, small molecular clusters are unlikely to contain more than one ammonia molecule per sulfuric acid. This implies that the average NH3:H2SO4 mole ratio of small molecular clusters in atmospheric conditions is likely to be between 1:3 and 1:1. Calculations on charged clusters confirm the experimental result that the HSO4- ion is much more strongly hydrated than neutral sulfuric acid. Preliminary calculations on HSO4- NH3 clusters indicate that ammonia is likely to play at most a minor role in ion-induced nucleation in the sulfuric acid - water system. Calculations of thermodynamic and kinetic parameters for the reaction of stabilized Criegee Intermediates with sulfuric acid demonstrate that quantum chemistry is a powerful tool for investigating chemically complicated nucleation mechanisms. The calculations indicate that if the biogenic Criegee Intermediates have sufficiently long lifetimes in atmospheric conditions, the studied reaction may be an important source of nucleation precursors.

How the development of condensation particle counters is reforming our view on atmospheric nucleation

Aerosol particles deteriorate air quality, atmospheric visibility and our health. They affect the Earth s climate by absorbing and scattering sunlight, forming clouds, and also via several feed-back mechanisms. The net effect on the radiative balance is negative, i.e. cooling, which means that particles counteract the effect of greenhouse gases. However, particles are one of the poorly known pieces in the climate puzzle. Some of the airborne particles are natural, some anthropogenic; some enter the atmosphere in particle form, while others form by gas-to-particle conversion. Unless the sources and dynamical processes shaping the particle population are quantified, they cannot be incorporated into climate models. The molecular level understanding of new particle formation is still inadequate, mainly due to the lack of suitable measurement techniques to detect the smallest particles and their precursors. This thesis has contributed to our ability to measure newly formed particles. Three new condensation particle counter applications for measuring the concentration of nano-particles were developed. The suitability of the methods for detecting both charged and electrically neutral particles and molecular clusters as small as 1 nm in diameter was thoroughly tested both in laboratory and field conditions. It was shown that condensation particle counting has reached the size scale of individual molecules, and besides measuring the concentration they can be used for getting size information. In addition to atmospheric research, the particle counters could have various applications in other fields, especially in nanotechnology. Using the new instruments, the first continuous time series of neutral sub-3 nm particle concentrations were measured at two field sites, which represent two different kinds of environments: the boreal forest and the Atlantic coastline, both of which are known to be hot-spots for new particle formation. The contribution of ions to the total concentrations in this size range was estimated, and it could be concluded that the fraction of ions was usually minor, especially in boreal forest conditions. Since the ionization rate is connected to the amount of cosmic rays entering the atmosphere, the relative contribution of neutral to charged nucleation mechanisms extends beyond academic interest, and links the research directly to current climate debate.

Molecular pathogenesis of human parechovirus 1 infection

The Parechoviruses (HPEV) belong to the family Picornaviridae of positive-stranded RNA viruses. Although the parechovirus genome shares the general properties of other picornaviruses, the genus has several unique features when compared to other family members. We found that HPEV1 attaches to αv integrins on the cell surface and is internalized through the clathrin-mediated endocytic pathway. During he course of the infection, the Golgi was found to disintegrate and the ER membranes to swell and loose their ribosomes. The replication of HPEV1 was found to take place on small clusters of vesicles which contained the trans-Golgi marker GalT as well as the viral non-structural 2C protein. 2C was additionally found on stretches of modified ER-membranes, seemingly not involved in RNA replication. The viral non-structural 2A and 2C proteins were studied in further detail and were found to display several interesting features. The 2A protein was found to be a RNA-binding protein that preferably binds to positive sense 3 UTR RNA. It was found to bind also duplex RNA containing 3 UTR(+)-3 UTR(-), but not other dsRNA molecules studied. Mutagenesis revealed that the N-terminal basic-rich region as well as the C-terminus, are important for RNA-binding. The 2C protein on the other hand, was found to have both ATP-diphosphohydrolase and AMP kinase activities. Neither dATP nor other NTP:s were suitable substrates. Furthermore, we found that as a result of theses activities the protein is autophosphorylated. The intracellular changes brought about by the individual HPEV1 non-structural proteins were studied through the expression of fusion proteins. None of the proteins expressed were able to induce membrane changes similar to those seen during HPEV1 infection. However, the 2C protein, which could be found on the surface of lipid droplets but also on diverse intracellular membranes, was partly relocated to viral replication complexes in transfected, superinfected cells. Although Golgi to ER traffic was arrested in HPEV1-infected cells, none of the individually expressed non-structural proteins had any visible effect on the anterograde membrane traffic. Our results suggest that the HPEV1 replication strategy is different from that of many other picornaviruses. Furthermore, this study shows how relatively small differences in genome sequence result in very different intracellular pathology.

Molecular profiling of malignant pleural mesothelioma and pulmonary fibrosis

Malignant mesothelioma (MM) is a rare, usually incurable, disease mainly caused by former exposure to asbestos. Even though MM has a strong etiological link, genetic factors may play a role, since not all cases can be linked to former asbestos exposure. This thesis focuses on lung diseases, mainly malignant mesothelioma (MM), and idiopathic pulmonary fibrosis (IPF), which resembles asbestosis. The specific asbestos-related pathways associated with malignant as well as non-malignant lung diseases, still need to be clarified. Since most patients diagnosed with MM or asbestosis/fibrosis have a dismal prognosis and few therapeutic options are available, early diagnosis and better understanding of the disease pathogenesis are of the utmost importance. The first objective of this thesis was to identify asbestos specific differentially expressed genes. This was approached by using high-resolution gene expression arrays, and three different human lung cell lines, as well as with three different bioinformatics approaches. Since the first study aimed to elucidate potential early changes, the second study was used to screen DNA copy number changes in MM tumour samples. This was performed using genome wide microarrays for identification of DNA copy number changes characterstic for MM. Study III focused on the role of gremlin in the regulation of bone morphogenetic protein (BMPs) in IPF. Further studies were conducted in asbestos-exposed cell cultures as well as in an asbestos-induced mouse model. Furthermore, GATA-6 was studied in MM and metastatic pleural adenocarcinoma. The GATA transcription factors are important during embryonic development, but their role in cancer is still unclear. GATA-6 is a co-factor/target of thyroid transcription factor 1 (TTF-1), which is used in differential diagnostics of pleural MM and adenocarcinoma. Bioinformatics probed the genes and biological processes ordered in terms of significance, clusters, and highly enriched chromosomal regions. The study revealed several already identified targets, produced new ideas about genes which are central for asbestos exposure, as well as provided supplementary data for researchers to check their own novel findings or ideas. The analysis revealed DNA copy number changes characteristic for MM tumors. The most common regions of loss were detected in 1p, 3p, 6q, 9p, 13, 14, and 22, and gains at 17q. The histological features in asbestosis and IPF are very similar, wherefore IPF can be studied in asbestos models. The BMP antagonist gremlin was up-regulated by asbestos exposure in human epithelial cell lines, which was also observed in Study I. The transforming growth factor (TGF) -β and BMP expression and signaling activities were measured from murine and human fibrotic lungs. BMP-7 signaling was down-regulated in response to up-regulation of gremlin, and restoration of BMP-7 signaling prevented progression of fibrosis in mice. Therefore, the study suggests that the restoration of BMP-7 signaling in fibrotic lung could potentially aid in the treatment of IPF patients. Study IV revealed that GATA-6 was strongly expressed in the majority of the MM cases, and correlated statistically significant with longer survival in subgroups of MM.

Foodborne Yersinia : Identification and Molecular Epidemiology of Isolates from Human Infections

Yersinia enterocolitica and Yersinia pseudotuberculosis are among the major enteropathogenic bacteria causing infections in humans in many industrialized countries. In Finland, Y. pseudotuberculosis has caused 10 outbreaks among humans during 1997-2008. Some of these outbreaks have been very extensive involving over 400 cases; mainly children attending schools and day-care. Y. enterocolitica, on the contrary, has caused mainly a large number of sporadic human infections in Finland. Y. pseudotuberculosis is widespread in nature, causing infections in a variety of domestic and wild animals. Foodborne transmission of human infections has long been suspected, however, attempts to trace the pathogen have been unsuccessful before this study that epidemiologically linked Y. pseudotuberculosis to a specific food item. Furthermore, due to modern food distribution systems, foodborne outbreaks usually involve many geographically separate infection clusters difficult to identify as part of the same outbreak. Among pathogenic Y. enterocolitica, the global predominance of one genetically homogeneous type (bioserotype 4/O:3) is a challenge to the development of genetic typing methods discriminatory enough for epidemiological purposes, for example, for tracing back to the sources of infections. Furthermore, the diagnostics of Y. enterocolitica infections is hampered because clinical laboratories easily misidentify some other members of the Yersinia species (Y. enterocolitica–like species) as Y. enterocolitica. This results in misleading information on the prevalence and clinical significance of various Yersinia isolates. The aim of this study was to develop and optimize molecular typing methods to be used in epidemiological investigations of Y. enterocolitica and Y. pseudotuberculosis, particularly in active surveillance and outbreak investigations of Y. pseudotuberculosis isolates. The aim was also to develop a simplified set of phenotypic tests that could be used in routine diagnostic laboratories for the correct identification of Y. enterocolitica and Y. enterocolitica –like species. A PFGE method designed here for typing of Y. pseudotuberculosis was efficient in linking the geographically dispersed and apparently unrelated Y. pseudotuberculosis infections as parts of the same outbreak. It proved to be useful in active laboratory-based surveillance of Y. pseudotuberculosis outbreaks. Throughout the study period, information about the diversity of genotypes among outbreak and non-outbreak related strains of human origin was obtained. Also, to our knowledge, this was the first study to epidemiologically link a Y. pseudotuberculosis outbreak of human illnesses to a specific food item, iceberg lettuce. A novel epidemiological typing method based on the use of a repeated genomic region (YeO:3RS) as a probe was developed for the detection and differentiation between strains of Y. enterocolitica subspecies palearctica. This method was able to increase the discrimination in a set of 106 previously PFGE typed Finnish Y. enterocolitica bioserotype 4/O:3 strains among which two main PFGE genotypes had prevailed. The developed simplified method was a more reliable tool than the commercially available biochemical test kits for differentiation between Y. enterocolitica and Y. enterocolitica –like species. In Finland, the methods developed for Y. enterocolitica and Y. pseudotuberculosis have been used to improve the identification protocols and in subsequent outbreak investigations.

Computational Studies on Bimetallic Clusters

The chemical and physical properties of bimetallic clusters have attracted considerable attention due to the potential technological applications of mixed-metal systems. It is of fundamental interests to study clusters because they are the link between atomic surface and bulk properties. More information of metal-metal bond in small clusters can be hence released. The studies in my thesis mainly focus on the two different kinds of bimetallic clusters: the clusters consisting of extraordinary shaped all metal four-membered rings and a series of sodium auride clusters. As described in most general organic chemistry books nowadays, a group of compounds are classified as aromatic compounds because of their remarkable stabilities, particular geometrical and energetic properties and so on. The notation of aromaticity is essentially qualitative. More recently, the connection has been made between aromaticity and energetic and magnetic properties. Also, the discussions of the aromatic nature of molecular rings are no longer limited to organic compounds obeying the Hückel’s rule. In our research, we mainly applied the GIMIC method to several bimetallic clusters at the CCSD level, and compared the results with those obtained by using chemical shift based methods. The magnetically induced ring currents can be generated easily by employing GIMIC method, and the nature of aromaticity for each system can be therefore clarified. We performed intensive quantum chemical calculations to explore the characters of the anionic sodium auride clusters and the corresponding neutral clusters since it has been fascinating in investigating molecules with gold atom involved due to its distinctive physical and chemical properties. As small gold clusters, the sodium auride clusters seem to form planar structures. With the addition of a negative charge, the gold atom in anionic clusters prefers to carry the charge and orients itself away from other gold atoms. As a result, the energetically lowest isomer for an anionic cluster is distinguished from the one for the corresponding neutral cluster. Mostly importantly, we presented a comprehensive strategy of ab initio applications to computationally implement the experimental photoelectron spectra.

Discovering hidden structures in molecular data using a Bayesian partition model approach

Advancements in the analysis techniques have led to a rapid accumulation of biological data in databases. Such data often are in the form of sequences of observations, examples including DNA sequences and amino acid sequences of proteins. The scale and quality of the data give promises of answering various biologically relevant questions in more detail than what has been possible before. For example, one may wish to identify areas in an amino acid sequence, which are important for the function of the corresponding protein, or investigate how characteristics on the level of DNA sequence affect the adaptation of a bacterial species to its environment. Many of the interesting questions are intimately associated with the understanding of the evolutionary relationships among the items under consideration. The aim of this work is to develop novel statistical models and computational techniques to meet with the challenge of deriving meaning from the increasing amounts of data. Our main concern is on modeling the evolutionary relationships based on the observed molecular data. We operate within a Bayesian statistical framework, which allows a probabilistic quantification of the uncertainties related to a particular solution. As the basis of our modeling approach we utilize a partition model, which is used to describe the structure of data by appropriately dividing the data items into clusters of related items. Generalizations and modifications of the partition model are developed and applied to various problems. Large-scale data sets provide also a computational challenge. The models used to describe the data must be realistic enough to capture the essential features of the current modeling task but, at the same time, simple enough to make it possible to carry out the inference in practice. The partition model fulfills these two requirements. The problem-specific features can be taken into account by modifying the prior probability distributions of the model parameters. The computational efficiency stems from the ability to integrate out the parameters of the partition model analytically, which enables the use of efficient stochastic search algorithms.

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.

Molecular dynamics studies of nanoparticle impacts

This thesis concerns the dynamics of nanoparticle impacts on solid surfaces. These impacts occur, for instance, in space, where micro- and nanometeoroids hit surfaces of planets, moons, and spacecraft. On Earth, materials are bombarded with nanoparticles in cluster ion beam devices, in order to clean or smooth their surfaces, or to analyse their elemental composition. In both cases, the result depends on the combined effects of countless single impacts. However, the dynamics of single impacts must be understood before the overall effects of nanoparticle radiation can be modelled. In addition to applications, nanoparticle impacts are also important to basic research in the nanoscience field, because the impacts provide an excellent case to test the applicability of atomic-level interaction models to very dynamic conditions. In this thesis, the stopping of nanoparticles in matter is explored using classical molecular dynamics computer simulations. The materials investigated are gold, silicon, and silica. Impacts on silicon through a native oxide layer and formation of complex craters are also simulated. Nanoparticles up to a diameter of 20 nm (315000 atoms) were used as projectiles. The molecular dynamics method and interatomic potentials for silicon and gold are examined in this thesis. It is shown that the displacement cascade expansionmechanism and crater crown formation are very sensitive to the choice of atomic interaction model. However, the best of the current interatomic models can be utilized in nanoparticle impact simulation, if caution is exercised. The stopping of monatomic ions in matter is understood very well nowadays. However, interactions become very complex when several atoms impact on a surface simultaneously and within a short distance, as happens in a nanoparticle impact. A high energy density is deposited in a relatively small volume, which induces ejection of material and formation of a crater. Very high yields of excavated material are observed experimentally. In addition, the yields scale nonlinearly with the cluster size and impact energy at small cluster sizes, whereas in macroscopic hypervelocity impacts, the scaling 2 is linear. The aim of this thesis is to explore the atomistic mechanisms behind the nonlinear scaling at small cluster sizes. It is shown here that the nonlinear scaling of ejected material yield disappears at large impactor sizes because the stopping mechanism of nanoparticles gradually changes to the same mechanism as in macroscopic hypervelocity impacts. The high yields at small impactor size are due to the early escape of energetic atoms from the hot region. In addition, the sputtering yield is shown to depend very much on the spatial initial energy and momentum distributions that the nanoparticle induces in the material in the first phase of the impact. At the later phases, the ejection of material occurs by several mechanisms. The most important mechanism at high energies or at large cluster sizes is atomic cluster ejection from the transient liquid crown that surrounds the crater. The cluster impact dynamics detected in the simulations are in agreement with several recent experimental results. In addition, it is shown that relatively weak impacts can induce modifications on the surface of an amorphous target over a larger area than was previously expected. This is a probable explanation for the formation of the complex crater shapes observed on these surfaces with atomic force microscopy. Clusters that consist of hundreds of thousands of atoms induce long-range modifications in crystalline gold.

Molecular Dynamics Simulations of Homogeneous Nucleation

Nucleation is the first step in a phase transition where small nuclei of the new phase start appearing in the metastable old phase, such as the appearance of small liquid clusters in a supersaturated vapor. Nucleation is important in various industrial and natural processes, including atmospheric new particle formation: between 20 % to 80 % of atmospheric particle concentration is due to nucleation. These atmospheric aerosol particles have a significant effect both on climate and human health. Different simulation methods are often applied when studying things that are difficult or even impossible to measure, or when trying to distinguish between the merits of various theoretical approaches. Such simulation methods include, among others, molecular dynamics and Monte Carlo simulations. In this work molecular dynamics simulations of the homogeneous nucleation of Lennard-Jones argon have been performed. Homogeneous means that the nucleation does not occur on a pre-existing surface. The simulations include runs where the starting configuration is a supersaturated vapor and the nucleation event is observed during the simulation (direct simulations), as well as simulations of a cluster in equilibrium with a surrounding vapor (indirect simulations). The latter type are a necessity when the conditions prevent the occurrence of a nucleation event in a reasonable timeframe in the direct simulations. The effect of various temperature control schemes on the nucleation rate (the rate of appearance of clusters that are equally able to grow to macroscopic sizes and to evaporate) was studied and found to be relatively small. The method to extract the nucleation rate was also found to be of minor importance. The cluster sizes from direct and indirect simulations were used in conjunction with the nucleation theorem to calculate formation free energies for the clusters in the indirect simulations. The results agreed with density functional theory, but were higher than values from Monte Carlo simulations. The formation energies were also used to calculate surface tension for the clusters. The sizes of the clusters in the direct and indirect simulations were compared, showing that the direct simulation clusters have more atoms between the liquid-like core of the cluster and the surrounding vapor. Finally, the performance of various nucleation theories in predicting simulated nucleation rates was investigated, and the results among other things highlighted once again the inadequacy of the classical nucleation theory that is commonly employed in nucleation studies.

Cathepsin D Deficiency - Molecular and Cellular Mechanisms of Neurodegeneration

Cathepsin D (CTSD) is a lysosomal protease, the deficiency of which is fatal and associated with neurodegeneration. CTSD knock-out mice, which die at the age of four weeks, show intestinal necrosis, loss of lymphoid cells and moderate pathological changes in the brain. An active-site mutation in the CTSD gene underlies a neurodegenerative disease in newborn sheep, characterized by brain atrophy without any changes to visceral tissues. The CTSD deficiences belong to the group of neuronal ceroid-lipofuscinoses (NCLs), severe neurodegenerative lysosomal storage disorders. The aim of this thesis was to examine the molecular and cellular mechanisms behind neurodegeneration in CTSD deficiency. We found the developmental expression pattern of CTSD to resemble that of synaptophysin and the increasing expression of CTSD to coincide with the active period of myelination in the rat brain, suggesting a role for CTSD in early rat brain development. An active-site mutation underlying the congenital ovine NCL not only affected enzymatic activity, but also changed the stability, processing and transport of the mutant protein, possibly contributing to the disease pathogenesis. We also provide CTSD deficiency as a first molecular explanation for human congenital NCL, a lysosomal storage disorder, characterized by neuronal loss and demyelination in the central nervous system. Finally, we show the first evidence for synaptic abnormalities and thalamocortical changes in CTSD-deficient mice at the molecular and ultrastructural levels. Keywords: cathepsin D, congenital, cortex, lysosomal storage disorder, lysosome, mutation, neurodegeneration, neuronal ceroid-lipofuscinosis, overexpression, synapse, thalamus