Simulation and graph mining tools for improving gene mapping efficiency

Gene mapping is a systematic search for genes that affect observable characteristics of an organism. In this thesis we offer computational tools to improve the efficiency of (disease) gene-mapping efforts. In the first part of the thesis we propose an efficient simulation procedure for generating realistic genetical data from isolated populations. Simulated data is useful for evaluating hypothesised gene-mapping study designs and computational analysis tools. As an example of such evaluation, we demonstrate how a population-based study design can be a powerful alternative to traditional family-based designs in association-based gene-mapping projects. In the second part of the thesis we consider a prioritisation of a (typically large) set of putative disease-associated genes acquired from an initial gene-mapping analysis. Prioritisation is necessary to be able to focus on the most promising candidates. We show how to harness the current biomedical knowledge for the prioritisation task by integrating various publicly available biological databases into a weighted biological graph. We then demonstrate how to find and evaluate connections between entities, such as genes and diseases, from this unified schema by graph mining techniques. Finally, in the last part of the thesis, we define the concept of reliable subgraph and the corresponding subgraph extraction problem. Reliable subgraphs concisely describe strong and independent connections between two given vertices in a random graph, and hence they are especially useful for visualising such connections. We propose novel algorithms for extracting reliable subgraphs from large random graphs. The efficiency and scalability of the proposed graph mining methods are backed by extensive experiments on real data. While our application focus is in genetics, the concepts and algorithms can be applied to other domains as well. We demonstrate this generality by considering coauthor graphs in addition to biological graphs in the experiments.

Nucleophilic attacks of 1,2-diaminoethane on MeCN ligands: synthesis, x-ray structure, and spectral and electrochemical properties of [Ru(.mu.-O)(.mu.-O2CAr)2[NH2CH2CH2NHC(Me)NH]2(PPh3)2](ClO4)2(Ar = C6H4-p-X; X = H, Me, OMe, Cl)

The diruthenium(III) complex [Ru2O(O2CAr)2(MeCN)4(PPh3)2](ClO4)2 (1), on reaction with 1,2-diaminoethane (en) in MeOH at 25-degrees-C, undergoes nucleophilic attacks at the carbon of two facial MeCN ligands to form [(Ru2O)-O-III(O2CAr)2-{NH2CH2CH2NHC(Me)NH}2(PPh3)2](ClO4)2 (2) (Ar = C6H4-p-X, X = H, Me, OMe, Cl) containing two seven-membered amino-amidine chelating ligands. The molecular structure of 2 with Ar = C6H4-p-OMe was determined by X-ray crystallography. Crystal data are as follows: triclinic, P1BAR, a = 13.942 (5) angstrom, b = 14.528 (2) angstrom, c = 21.758 (6) angstrom, alpha = 109.50 (2)-degrees, beta = 92.52 (3)-degrees, gamma = 112.61 (2)-degrees, V = 3759 (2) angstrom 3, and Z = 2. The complex has an {Ru2(mu-O)(mu-O2CAr2)2(2+)} core. The Ru-Ru and average Ru-O(oxo) distances and the Ru-O-Ru angle are 3.280 (2) angstrom, 1.887 [8] angstrom, and 120.7 (4)-degrees, respectively. The amino group of the chelating ligand is trans to the mu-oxo ligand. The nucleophilic attacks take place on the MeCN ligands cis to the mu-oxo ligand. The visible spectra of 2 in CHCl3 display an absorption band at 565 nm. The H-1 NMR spectra of 2 in CDCl3 are indicative of the formation of an amino-amidine ligand. Complex 2 exhibits metal-centered quasireversible one-electron oxidation and reduction processes in the potential ranges +0.9 to +1.0 V and -0.3 to -0.5 V (vs SCE), respectively, involving the Ru(III)2/Ru(III)Ru(IV) and Ru(III)2/Ru(II)Ru(III) redox couples in CH2Cl2 containing 0.1 M TBAP. The mechanistic aspects of the nucleophilic reaction are discussed.

Towards Understanding Powder Behavior Via Simulation

The aim of this study was to investigate powder and tablet behavior at the level of mechanical interactions between single particles. Various aspects of powder packing, mixing, compression, and bond formation were examined with the aid of computer simulations. The packing and mixing simulations were based on spring forces interacting between particles. Packing and breakage simulations included systems in which permanent bonds were formed and broken between particles, based on their interaction strengths. During the process, a new simulation environment based on Newtonian mechanics and elementary interactions between the particles was created, and a new method for evaluating mixing was developed. Powder behavior is a complicated process, and many of its aspects are still unclear. Powders as a whole exhibit some aspects of solids and others of liquids. Therefore, their physics is far from clear. However, using relatively simple models based on particle-particle interaction, many powder properties could be replicated during this work. Simulated packing densities were similar to values reported in the literature. The method developed for describing powder mixing correlated well with previous methods. The new method can be applied to determine mixing in completely homogeneous materials, without dividing them into different components. As such, it can describe the efficiency of the mixing method, regardless of the powder's initial setup. The mixing efficiency at different vibrations was examined, and we found that certain combinations of amplitude, direction, and frequencies resulted in better mixing while using less energy. Simulations using exponential force potentials between particles were able to explain the elementary compression behavior of tablets, and create force distributions that were similar to the pressure distributions reported in the literature. Tablet-breaking simulations resulted in breaking strengths that were similar to measured tablet breaking strengths. In general, many aspects of powder behavior can be explained with mechanical interactions at the particle level, and single particle properties can be reliably linked to powder behavior with accurate simulations.

Preparation, characterization, spectral and thermal analyses of (N2H5)2MCl4·2H2O (M = Fe, Co, Ni and Cu); crystal structure of the iron complex

Hydrazinium metal chlorides, (N2H5)2MCl4·2H2O (where M = Fe, Co, Ni and Cu), have been prepared from the aqueous solutions of the respective metal chlorides and hydrazine hydrochloride (N2H4·HCl or N2H4·2HCl) and investigated by spectral and thermal analyses. The crystal structure of the iron complex has been determined by direct methods and refined by full-matrix least-squares to an R of 0.023 and Rw of 0.031 for 1495 independent reflections. The structure shows ferrous ion in an octahedral environment bonded by two hydrazinium cations, two chloride anions and two water molecules. In the complex cation [Fe(N2H5)2(H2O)2Cl2]2+, the coordinated groups are in trans positions.

Synthesis and spectral studies of polyamide-phosphate esters from phosphine-oxide-containing diols with aryl phosphorodichloridates and their thermal and flammability behaviour

Polyamide-phosphate esters were synthesized by interfacial polycondensation of aryl phosphorodichloridates with the diols of phenoxaphosphine and phosphine oxide in the presence of a phase-transfer catalyst. The polymers were characterized by infra-red and 1H, 13C and 31P nuclear magnetic resonance (n.m.r.) spectroscopy. The molecular weights were determined by end-group analysis using 31P n.m.r. spectral data. The phenoxaphosphine-containing polymers showed superior thermostability and flame retardancy over the phosphine-oxide-containing polymers.

Synthesis, spectral and structural studies on metal complexes of Schiff bases involving vitamin B6 and histamine

Six metal complexes of Schiff bases involving Vitamin B6 and the decarboxylated amino acid histamine have been synthesised and characterized. Crystal structures have been determined for [CuL1(H2O)Br]-NO31(L1= pyridoxylidenehistamine) and [Cu2L22(NO3)2]·6H2O 2(L2= 5′-phosphopyridoxylidenehistaminate). The crystal structure of complex 1[space group P[1 with combining macron], a= 8.161(2), b= 10.368(2), c= 11.110(2)Å, α= 105.18(1), β= 102.12(1), γ= 72.10(1)° and Z= 2; R= 0.072, R′= 0.083] consists of square-pyramidally co-ordinated copper with the tridentate Schiff base in the zwitterionic form, whereas in 2[space group P[1 with combining macron], a= 8.727(1), b= 10.308(1), c= 12.845(2)Å, α= 110.00(1), β= 78.94(1), γ= 114.35(1)° and Z= 1; R= 0.035, R′= 0.034] the copper has the same co-ordination geometry but the tetradentate Schiff-base ligand exists as a monoanion. The conformational parameters deduced from such structures are important for understanding the stereochemical aspects of Vitamin B6-catalysed model reactions involving histidine.

Optimizing simulation model on forest policy.

XVIII IUFRO World Congress, Ljubljana 1986.

Forest sector simulation models as methodological tools in forest policy analysis.

XVIII IUFRO World Congress, Ljubljana 1986.

On the holographic simulation of non-circular zone plates

New methods involving the manipulation of fundamental wavefronts (e.g., plane and spherical) with simple optical components such as pinholes and spherical lenses have been developed for the fabrication of elliptic, hyperbolic and conical holographic zone plates. Also parabolic zone plates by holographic techniques have been obtained for the first time. The performance behaviour of these zone plates has been studied. Further a phenomenological explanation is offered for the observed improved fringe contrast obtained with a spherical reference wave.

Application of an ecological simulation model to Lake Päijänne

Tiivistelmä: Simulointimallin soveltaminen Pohjois-Päijänteellä

Direct Growth of Gold Nanorods on Gold and Indium Tin Oxide Surfaces: Spectral, Electrochemical, and Electrocatalytic Studies

The surfactant-assisted seed-mediated growth method was used for the formation of gold nanorods (GNRs) directly on gold (Au) and indium tin oxide (ITO) surfaces. Citrate-stabilized similar to 2.6 nm spherical gold nanoparticles (AuNPs) were first self-assembled on ITO or Au surfaces modified with (3-mercaptopropyl)-trimethoxysilane (MPTS) sol-gel film and then immersed in a cationic surfactant growth solution to form GNRs. The growth of GNRs on the MPTS sol gel film modified ITO surface was monitored by UV-visible spectroscopy. The ITO surface with the attached spherical AuNPs shows a surface plasmon resonance band at 550 nm. The intensity of this absorption band increases while increasing the immersion time of the AuNP-modified ITO surface into the growth solution, and after 5 h, an additional shoulder band around 680 nm was observed. The intensity of this shoulder band increases, and it was shifted to longer wavelength as the immersion time of the AuNP-modified ITO surface into the growth solution increases. After 20 h, a predominant wave at 720 nm was observed along with a band at 550 nm. Further immersion of the modified ITO surface into the growth solution did not change the absorption characteristics. The bands observed at 550 and 720 nm were characteristics of GNRs, corresponding to transverse and longitudinal waves, respectively. The AFM images showed the presence of GNRs on the surface of the MPTS sol gel modified ITO surface with a typical length of similar to 100-120 nm and a width of similar to 20-22 nm in addition to a few spherical AuNPs, indicating that seeded spherical AuNPs were not completely involved in the GNRs' formation. Finally, the electrocatalytic activity of the surface-grown GNRs on the MPTS sol gel film modified Au electrode toward the oxidation of ascorbic acid (AA) was studied. Unlike a polycrystalline Au electrode, the surface-grown GNR-modified electrode shows two well-defined voltammetric peaks for AA at 0.01 and 0.35 V in alkaline, neutral, and acidic pHs. The cause for the observed two oxidation peaks for AA was due to the presence of both nanorods and spherical nanoparticles on the electrode surface. The presence of spherical AuNPs on the MPTS sol gel film oxidized AA at more positive potential, whereas the GNRs oxidized AA at less positive potential. The observed 340 mV less positive potential shift in the oxidation of AA suggested that GNRs are better electrocatalysts for the oxidation of AA than the spherical AuNPs.

Computer simulation of multi-elemental fusion reactor materials

Thermonuclear fusion is a sustainable energy solution, in which energy is produced using similar processes as in the sun. In this technology hydrogen isotopes are fused to gain energy and consequently to produce electricity. In a fusion reactor hydrogen isotopes are confined by magnetic fields as ionized gas, the plasma. Since the core plasma is millions of degrees hot, there are special needs for the plasma-facing materials. Moreover, in the plasma the fusion of hydrogen isotopes leads to the production of high energetic neutrons which sets demanding abilities for the structural materials of the reactor. This thesis investigates the irradiation response of materials to be used in future fusion reactors. Interactions of the plasma with the reactor wall leads to the removal of surface atoms, migration of them, and formation of co-deposited layers such as tungsten carbide. Sputtering of tungsten carbide and deuterium trapping in tungsten carbide was investigated in this thesis. As the second topic the primary interaction of the neutrons in the structural material steel was examined. As model materials for steel iron chromium and iron nickel were used. This study was performed theoretically by the means of computer simulations on the atomic level. In contrast to previous studies in the field, in which simulations were limited to pure elements, in this work more complex materials were used, i.e. they were multi-elemental including two or more atom species. The results of this thesis are in the microscale. One of the results is a catalogue of atom species, which were removed from tungsten carbide by the plasma. Another result is e.g. the atomic distributions of defects in iron chromium caused by the energetic neutrons. These microscopic results are used in data bases for multiscale modelling of fusion reactor materials, which has the aim to explain the macroscopic degradation in the materials. This thesis is therefore a relevant contribution to investigate the connection of microscopic and macroscopic radiation effects, which is one objective in fusion reactor materials research.

Large Eddy Simulation

Large eddy simulation (LES) is an emerging technique for obtaining an approximation to turbulent flow fields It is an improvement over the widely prevalent practice of obtaining means of turbulent flows when the flow has large scale, low frequency, unsteadiness An introduction to the method, its general formulation, and the more common modelling for flows without reaction, is discussed Some attempts at extension to flows with combustion have been made Examples from present work for flows with and without combustion are given The final example of the LES of the combustor of a helicopter engine illustrates the state-of-the-art in application of the technique

Analytical and simulation studies of a multiprocessor system for high-speed numerical computations

This paper describes the architecture of a multiprocessor system which we call the Broadcast Cube System (BCS) for solving important computation intensive problems such as systems of linear algebraic equations and Partial Differential Equations (PDEs), and highlights its features. Further, this paper presents an analytical performance study of the BCS, and it describes the main details of the design and implementation of the simulator for the BCS.