Alpha2-adrenoceptors: Structure and ligand binding properties at the molecular level

Alpha₂-Adrenoceptors: structure and ligand binding properties at the molecular level The mouse is the most frequently used animal model in biomedical research, but the use of zebrafish as a model organism to mimic human diseases is on the increase. Therefore it is considered important to understand their pharmacological differences from humans also at the molecular level. The zebrafish Alpha_2-adrenoceptors were expressed in mammalian cells and the binding affinities of 20 diverse ligands were determined and compared to the corresponding human receptors. The pharmacological properties of the human and zebrafish Alpha_2--adrenoceptors were found to be quite well conserved. Receptor models based on the crystal structures of bovine rhodopsin and the human Beta₂-adrenoceptor revealed that most structural differences between the paralogous and orthologous Alpha_2--adrenoceptors were located within the second extracellular loop (XL2). Reciprocal mutations were generated in the mouse and human Alpha_2--adrenoceptors. Ligand binding experiments revealed that substitutions in XL2 reversed the binding profiles of the human and mouse Alpha_2--adrenoceptors for yohimbine, rauwolscine and RS-79948-197, evidence for a role for XL2 in the determination of species-specific ligand binding. Previous mutagenesis studies had not been able to explain the subtype preference of several large Alpha_2--adrenoceptor antagonists. We prepared chimaeric Alpha_2--adrenoceptors where the first transmembrane (TM1) domain was exchanged between the three human Alpha_2--adrenoceptor subtypes. The binding affinities of spiperone, spiroxatrine and chlorpromazine were observed to be significantly improved by TM1 substitutions of the Alpha_2a--adrenoceptor. Docking simulations indicated that indirect effects, such as allosteric modulation, are more likely to be involved in this phenomenon rather than specific side-chain interactions between ligands and receptors.

Calibration of time domain reflectometry (DTR) for soil moisture measurements in cultivated peat soils

Tiivistelmä: TDR-mittausten kalibrointi viljeltyjen turvemaiden kosteuden mittaamiseen

IPSec in 3G packet switched domain

Diplomityö käsittelee IPSec-protokollan (IP Security Protocol) implementointia UMTS:n (Universal Mobile Telecommunications System) pakettikytkentäiseen verkkoon. Runkoverkkoa käytetään mobiilikäyttäjän datan siirtämiseen sekä verkkoelementtien väliseen ohjausinformaation välitykseen. Koska UMTS:n runkoverkot ovat IP-pakettikytkentäisiä verkkoja, IPSec-protokollaa voidaan käyttää lähetettyjen IP-datasähkeiden suojaamiseen. IPSec- ja IKE-protokollien (Internet Key Exchange) käyttö on koettu monimutkaiseksi kiinteissä verkoissa. Tämän saman ongelman edessä tulevat olemaan myös operaattorit, kun he alkavat rakentaa UMTS-verkkojaan. On kuitenkin muistettava se, että tulevaisuudessa lähes kaikki data mukaanlukien ääni ja video on tarkoitus siirtää IP-protokollan avulla. IP-teknologiaan perustuva tiedonsiirron kasvu lisää IPSec-protokollan merkitystä ei ainoastaan runkoverkossa mutta myös radioliityntäverkoissa sekä SS7-merkinantoverkoissa (Signaling System No. 7). Diplomityö on tehty osaksi diplomi-insinöörin tutkintoa Lappeenrannan teknillisessä yliopistossa. Työ on tehty Nokia Networksin palveluksessa Helsingissä, vuosien 2002 ja 2003 välisenä aikana.

Study of ligand-based virtual screening tools in computer-aided drug design

Virtual screening is a central technique in drug discovery today. Millions of molecules can be tested in silico with the aim to only select the most promising and test them experimentally. The topic of this thesis is ligand-based virtual screening tools which take existing active molecules as starting point for finding new drug candidates. One goal of this thesis was to build a model that gives the probability that two molecules are biologically similar as function of one or more chemical similarity scores. Another important goal was to evaluate how well different ligand-based virtual screening tools are able to distinguish active molecules from inactives. One more criterion set for the virtual screening tools was their applicability in scaffold-hopping, i.e. finding new active chemotypes. In the first part of the work, a link was defined between the abstract chemical similarity score given by a screening tool and the probability that the two molecules are biologically similar. These results help to decide objectively which virtual screening hits to test experimentally. The work also resulted in a new type of data fusion method when using two or more tools. In the second part, five ligand-based virtual screening tools were evaluated and their performance was found to be generally poor. Three reasons for this were proposed: false negatives in the benchmark sets, active molecules that do not share the binding mode, and activity cliffs. In the third part of the study, a novel visualization and quantification method is presented for evaluation of the scaffold-hopping ability of virtual screening tools.

Conformational Regulation of alpha2beta1 Integrin in Ligand Binding

Integrins are heterodimeric cell adhesion receptors involved in cell-cell and cell-extracellular matrix (ECM) interactions. They transmit bidirectional signals across the cell membrane. This results in a wide range of biological events from cell differentiation to apoptosis. alpha2beta1 integrin is an abundant collagen receptor expressed on the surface of several cell types. In addition to ECM ligands, alpha2beta1 integrins are bound by echovirus 1 (EV1) which uses alpha2beta1 as a receptor to initiate its life cycle in the infected cell. The aim of this thesis project was to provide further insight into the mechanisms of alpha2beta1 integrin ligand recognition and receptor activation. Collagen fibrils are the principal tensile elements of the ECM. Yet, the interaction of alpha2beta1 integrin with the fibrillar form of collagen I has received relatively little attention. This research focused on the ability of alpha2beta1 integrin to act as a receptor for type I collagen fibrils. Also the molecular requirements of the EV1 interaction with alpha2beta1 were studied. Conventionally, ligand binding has been suggested to require integrin activation and the binding may further trigger integrin signalling. Another main objective of this study was to elucidate both the inside-out and outside-in signalling mechanisms of alpha2beta1 integrin in adherent cells. The results indicated that alpha2beta1 integrin is the principal integrin-type collagen receptor for type I collagen fibrils, and alpha2beta1 may participate in the regulation of pericellular collagen fibrillogenesis. Furthermore, alpha2beta1 integrin inside-out activation appeared to be synergistically regulated by integrin clustering and conformational activation. The triggering of alpha2beta1 integrin outside-in signalling, however, was shown to require both conformational changes and clustering. In contrast to ECM ligands, EV1 appeared to take advantage of the bent, inactive form of alpha2beta1 integrin in initiating its life cycle in the cell. This research together with other recent studies, has shed light on the molecular mechanisms of integrin activation. It is becoming evident that large ligands are able to bind to the bent form of integrin, which has been previously considered to be physiologically inactive. Consequently, our understanding of the conformational modulation of integrins upon activation is changing.

Utilization of a Time Domain Simulator in the Technical and Economic Analysis of a Wind Turbine Electric Drive Train

The amount of installed wind power has been growing exponentially during the past ten years. As wind turbines have become a significant source of electrical energy, the interactions between the turbines and the electric power network need to be studied more thoroughly than before. Especially, the behavior of the turbines in fault situations is of prime importance; simply disconnecting all wind turbines from the network during a voltage drop is no longer acceptable, since this would contribute to a total network collapse. These requirements have been a contributor to the increased role of simulations in the study and design of the electric drive train of a wind turbine. When planning a wind power investment, the selection of the site and the turbine are crucial for the economic feasibility of the installation. Economic feasibility, on the other hand, is the factor that determines whether or not investment in wind power will continue, contributing to green electricity production and reduction of emissions. In the selection of the installation site and the turbine (siting and site matching), the properties of the electric drive train of the planned turbine have so far been generally not been taken into account. Additionally, although the loss minimization of some of the individual components of the drive train has been studied, the drive train as a whole has received less attention. Furthermore, as a wind turbine will typically operate at a power level lower than the nominal most of the time, efficiency analysis in the nominal operating point is not sufficient. This doctoral dissertation attempts to combine the two aforementioned areas of interest by studying the applicability of time domain simulations in the analysis of the economicfeasibility of a wind turbine. The utilization of a general-purpose time domain simulator, otherwise applied to the study of network interactions and control systems, in the economic analysis of the wind energy conversion system is studied. The main benefits of the simulation-based method over traditional methods based on analytic calculation of losses include the ability to reuse and recombine existing models, the ability to analyze interactions between the components and subsystems in the electric drive train (something which is impossible when considering different subsystems as independent blocks, as is commonly done in theanalytical calculation of efficiencies), the ability to analyze in a rather straightforward manner the effect of selections other than physical components, for example control algorithms, and the ability to verify assumptions of the effects of a particular design change on the efficiency of the whole system. Based on the work, it can be concluded that differences between two configurations can be seen in the economic performance with only minor modifications to the simulation models used in the network interaction and control method study. This eliminates the need ofdeveloping analytic expressions for losses and enables the study of the system as a whole instead of modeling it as series connection of independent blocks with no lossinterdependencies. Three example cases (site matching, component selection, control principle selection) are provided to illustrate the usage of the approach and analyze its performance.

Genesis and Emplacement of Carbonatites and Lamprophyres in the Svecofennian Domain

A small carbonatite dyke swarm has been identified at Naantali, southwest Finland. Several swarms of shoshonitic lamprophyres are also known along the Archean-Proterozoic boundary in eastern Finland and northwest Russia. These intrusions, along with the carbonatite intrusion at Halpanen, eastern Finland, represent a stage of widespread low-volume mantle-sourced alkaline magmatism in the Svecofennian Domain. Using trace element and isotope geochemistry coupled with precise geochronology from these rocks, a model is presented for the Proterozoic metasomatic evolution of the Fennoscandian subcontinental lithospheric mantle. At ~2.2-2.06 Ga, increased biological production in shallow seas linked to continental rifting, resulted in increased burial rates of organic carbon. Subduction between ~1.93-1.88 Ga returned organic carbon-enriched sediments of mixed Archean and Proterozoic provenance to the mantle. Dehydration reactions supplied water to the mantle wedge, driving arc volcanism, while mica, amphibole and carbonate were brought deeper into the mantle with the subducting slab. The cold subducted slab was heated conductively from the surrounding warm mantle, while pressures continued to gradually increase as a result of crustal thickening. The sediments began to melt in a two stage process, first producing a hydrous alkaline silicate melt, which infiltrated the mantle wedge and crystallised as metasomatic veins. At higher temperatures, carbonatite melt was produced, which preferentially infiltrated the pre-existing metasomatic vein network. At the onset of post-collisional extension, deep fault structures formed, providing conduits for mantle melts to reach the upper crust. Low-volume partial melting of the enriched mantle at depths of at least 110 km led to the formation of first carbonatitic magma and subsequently lamprophyric magma. Carbonatite was emplaced in the upper crust at Naantali at 1795.7 ± 6.8 Ma; lamprophyres along the Archean-Proterozoic boundary were emplaced between 1790.1 ± 3.3 Ma and 1781 ± 20 Ma.

Novel Players in the Integrin Signaling Orchestra: TCPTP and MDGI

Metastases are the major cause of cancer deaths. Tumor cell dissemination from the primary tumor utilizes dysregulated cellular adhesion and upregulated proteolytic degradation of the extracellular matrix for progeny formation in distant organs. Integrins are transmembrane adhesive receptors mediating cellcell and cellmatrix interactions that are crucial for regulating cell migration, invasion, proliferation, and survival. Consequently, increased integrin activity is associated with augmented migration and invasion capacity in several cancer types. Heterodimeric integrins consist of an alpha - and beta-subunit that are held together in a bent conformation when the receptor is inactive, but extension and separation of subdomains is observed during receptor activation. Either inside-out or outside-in activation of receptors is possible through the intracellular molecule binding to an integrin cytoplasmic domain or extracellular ligand association with an integrin ectodomain, respectively. Several regulatory binding partners have been characterized for integrin cytoplasmic beta-domains, but the regulators interacting with the cytoplasmic alpha-domains have remained elusive. In this study, we performed yeast two-hybrid screens to identify novel binding partners for the cytoplasmic integrin alpha-domains. Further examination of two plausible candidates revealed a significant coregulatory role of an integrin alpha-subunit for cellular signaling processes. T-cell protein tyrosine phosphatase (TCPTP) showed a specific interaction with the cytoplasmic tail of integrin alpha1. This association stimulated TCPTP phosphatase activity, leading to negative regulation of epidermal growth factor receptor (EGFR) signaling and diminished anchorage-independent growth. Another candidate, mammary-derived growth inhibitor (MDGI), exhibited binding to several different integrin cytoplasmic alpha-tails through a conserved GFFKR sequence. MDGI overexpression in breast cancer cells altered EGFR trafficking and caused a remarkable accumulation of EGFR in the cytoplasm. We further demonstrated in vivo that MDGI expression induced a novel form of anti-EGFR therapy resistance. Moreover, MDGI binding to α-tails retained integrin in an inactive conformation attenuating integrin-mediated adhesion, migration, and invasion. In agreement with these results, sustained MDGI expression in breast cancer patients correlated with an increased 10-year distant disease-free survival. Taken together, the integrin signaling network is far from a complete view and future work will doubtless broaden our understanding further.

Boundary identification in the domain of a parabolic partial differential equation

Bakgrunden och inspirationen till föreliggande studie är tidigare forskning i tillämpningar på randidentifiering i metallindustrin. Effektiv randidentifiering möjliggör mindre säkerhetsmarginaler och längre serviceintervall för apparaturen i industriella högtemperaturprocesser, utan ökad risk för materielhaverier. I idealfallet vore en metod för randidentifiering baserad på uppföljning av någon indirekt variabel som kan mätas rutinmässigt eller till en ringa kostnad. En dylik variabel för smältugnar är temperaturen i olika positioner i väggen. Denna kan utnyttjas som insignal till en randidentifieringsmetod för att övervaka ugnens väggtjocklek. Vi ger en bakgrund och motivering till valet av den geometriskt endimensionella dynamiska modellen för randidentifiering, som diskuteras i arbetets senare del, framom en flerdimensionell geometrisk beskrivning. I de aktuella industriella tillämpningarna är dynamiken samt fördelarna med en enkel modellstruktur viktigare än exakt geometrisk beskrivning. Lösningsmetoder för den s.k. sidledes värmeledningsekvationen har många saker gemensamt med randidentifiering. Därför studerar vi egenskaper hos lösningarna till denna ekvation, inverkan av mätfel och något som brukar kallas förorening av mätbrus, regularisering och allmännare följder av icke-välställdheten hos sidledes värmeledningsekvationen. Vi studerar en uppsättning av tre olika metoder för randidentifiering, av vilka de två första är utvecklade från en strikt matematisk och den tredje från en mera tillämpad utgångspunkt. Metoderna har olika egenskaper med specifika fördelar och nackdelar. De rent matematiskt baserade metoderna karakteriseras av god noggrannhet och låg numerisk kostnad, dock till priset av låg flexibilitet i formuleringen av den modellbeskrivande partiella differentialekvationen. Den tredje, mera tillämpade, metoden kännetecknas av en sämre noggrannhet förorsakad av en högre grad av icke-välställdhet hos den mera flexibla modellen. För denna gjordes även en ansats till feluppskattning, som senare kunde observeras överensstämma med praktiska beräkningar med metoden. Studien kan anses vara en god startpunkt och matematisk bas för utveckling av industriella tillämpningar av randidentifiering, speciellt mot hantering av olinjära och diskontinuerliga materialegenskaper och plötsliga förändringar orsakade av “nedfallande” väggmaterial. Med de behandlade metoderna förefaller det möjligt att uppnå en robust, snabb och tillräckligt noggrann metod av begränsad komplexitet för randidentifiering.

Applications of graph transformation in tools for domain-specific modeling languages

The use of domain-specific languages (DSLs) has been proposed as an approach to cost-e ectively develop families of software systems in a restricted application domain. Domain-specific languages in combination with the accumulated knowledge and experience of previous implementations, can in turn be used to generate new applications with unique sets of requirements. For this reason, DSLs are considered to be an important approach for software reuse. However, the toolset supporting a particular domain-specific language is also domain-specific and is per definition not reusable. Therefore, creating and maintaining a DSL requires additional resources that could be even larger than the savings associated with using them. As a solution, di erent tool frameworks have been proposed to simplify and reduce the cost of developments of DSLs. Developers of tool support for DSLs need to instantiate, customize or configure the framework for a particular DSL. There are di erent approaches for this. An approach is to use an application programming interface (API) and to extend the basic framework using an imperative programming language. An example of a tools which is based on this approach is Eclipse GEF. Another approach is to configure the framework using declarative languages that are independent of the underlying framework implementation. We believe this second approach can bring important benefits as this brings focus to specifying what should the tool be like instead of writing a program specifying how the tool achieves this functionality. In this thesis we explore this second approach. We use graph transformation as the basic approach to customize a domain-specific modeling (DSM) tool framework. The contributions of this thesis includes a comparison of di erent approaches for defining, representing and interchanging software modeling languages and models and a tool architecture for an open domain-specific modeling framework that e ciently integrates several model transformation components and visual editors. We also present several specific algorithms and tool components for DSM framework. These include an approach for graph query based on region operators and the star operator and an approach for reconciling models and diagrams after executing model transformation programs. We exemplify our approach with two case studies MICAS and EFCO. In these studies we show how our experimental modeling tool framework has been used to define tool environments for domain-specific languages.

Frequency-Domain and Wide-Pulse Time-Domain Measurements of Lanthanide Luminescence and Lanthanide-Based Resonance Energy Transfer

Resonance energy transfer (RET) is a non-radiative transfer of the excitation energy from the initially excited luminescent donor to an acceptor. The requirements for the resonance energy transfer are: i) the spectral overlap between the donor emission spectrum and the acceptor absorption spectrum, ii) the close proximity of the donor and the acceptor, and iii) the suitable relative orientations of the donor emission and the acceptor absorption transition dipoles. As a result of the RET process the donor luminescence intensity and the donor lifetime are decreased. If the acceptor is luminescent, a sensitized acceptor emission appears. The rate of RET depends strongly on the donor–acceptor distance (r) and is inversely proportional to r6. The distance dependence of RET is utilized in binding assays. The proximity requirement and the selective detection of the RET-modified emission signal allow homogeneous separation free assays. The term lanthanide-based RET is used when luminescent lanthanide compounds are used as donors. The long luminescence lifetimes, the large Stokes’ shifts and the intense, sharply-spiked emission spectra of the lanthanide donors offer advantages over the conventional organic donor molecules. Both the organic lanthanide chelates and the inorganic up-converting phosphor (UCP) particles have been used as donor labels in the RET based binding assays. In the present work lanthanide luminescence and lanthanide-based resonance energy transfer phenomena were studied. Luminescence lifetime measurements had an essential role in the research. Modular frequency-domain and time-domain luminometers were assembled and used successfully in the lifetime measurements. The frequency-domain luminometer operated in the low frequency domain ( 100 kHz) and utilized a novel dual-phase lock-in detection of the luminescence. One of the studied phenomena was the recently discovered non-overlapping fluorescence resonance energy transfer (nFRET). The studied properties were the distance and temperature dependences of nFRET. The distance dependence was found to deviate from the Förster theory and a clear temperature dependence was observed whereas conventional RET was completely independent of the temperature. Based on the experimental results two thermally activated mechanisms were proposed for the nFRET process. The work with the UCP particles involved the measurement of the luminescence properties of the UCP particles synthesized in our laboratory. The goal of the UCP particle research is to develop UCP donor labels for binding assays. In the present work the effect of the dopant concentrations and the core–shell structure on the total up-conversion luminescence intensity, the red–green emission ratio, and the luminescence lifetime was studied. Also the non-radiative nature of the energy transfer from the UCP particle donors to organic acceptors was demonstrated for the first time in aqueous environment and with a controlled donor–acceptor distance.

The Structure and Function of αI Domains in Collagen Receptor and Leukocyte Integrins

Integrins are heterodimeric, signaling transmembrane adhesion receptors that connect the intracellular actin microfilaments to the extracellular matrix composed of collagens and other matrix molecules. Bidirectional signaling is mediated via drastic conformational changes in integrins. These changes also occur in the integrin αI domains, which are responsible for ligand binding by collagen receptor and leukocyte specific integrins. Like intact integrins, soluble αI domains exist in the closed, low affinity form and in the open, high affinity form, and so it is possible to use isolated αI domains to study the factors and mechanisms involved in integrin activation/deactivation. Integrins are found in all mammalian tissues and cells, where they play crucial roles in growth, migration, defense mechanisms and apoptosis. Integrins are involved in many human diseases, such as inflammatory, cardiovascular and metastatic diseases, and so plenty of effort has been invested into developing integrin specific drugs. Humans have 24 different integrins, four of which are collagen receptor (α1β1, α2β1, α10β1, α11β1) and five leukocyte specific integrins (αLβ2, αMβ2, αXβ2, αDβ2, αEβ7). These two integrin groups are quite unselective having both primary and secondary ligands. This work presents the first systematic studies performed on these integrin groups to find out how integrin activation affects ligand binding and selectivity. These kinds of studies are important not only for understanding the partially overlapping functions of integrins, but also for drug development. In general, our results indicated that selectivity in ligand recognition is greatly reduced upon integrin activation. Interestingly, in some cases the ligand binding properties of integrins have been shown to be cell type specific. The reason for this is not known, but our observations suggest that cell types with a higher integrin activation state have lower ligand selectivity, and vice versa. Furthermore, we solved the three-dimensional structure for the activated form of the collagen receptor α1I domain. This structure revealed a novel intermediate conformation not previously seen with any other integrin αI domain. This is the first 3D structure for an activated collagen receptor αI domain without ligand. Based on the differences between the open and closed conformation of the αI domain we set structural criteria for a search for effective collagen receptor drugs. By docking a large number of molecules into the closed conformation of the α2I domain we discovered two polyketides, which best fulfilled the set structural criteria, and by cell adhesion studies we showed them to be specific inhibitors of the collagen receptor integrins.

Prediction of 3D structure and ligand-binding properties : ABC transporters and flavodiiron proteins from Synechocystis sp. strain PCC 6803

Cyanobacteria are unicellular, non-nitrogen-fixing prokaryotes, which perform photosynthesis similarly as higher plants. The cyanobacterium Synechocystis sp. strain PCC 6803 is used as a model organism in photosynthesis research. My research described herein aims at understanding the function of the photosynthetic machinery and how it responds to changes in the environment. Detailed knowledge of the regulation of photosynthesis in cyanobacteria can be utilized for biotechnological purposes, for example in the harnessing of solar energy for biofuel production. In photosynthesis, iron participates in electron transfer. Here, we focused on iron transport in Synechocystis sp. strain PCC 6803 and particularly on the environmental regulation of the genes encoding the FutA2BC ferric iron transporter, which belongs to the ABC transporter family. A homology model built for the ATP-binding subunit FutC indicates that it has a functional ATPbinding site as well as conserved interactions with the channel-forming subunit FutB in the transporter complex. Polyamines are important for the cell proliferation, differentiation and apoptosis in prokaryotic and eukaryotic cells. In plants, polyamines have special roles in stress response and in plant survival. The polyamine metabolism in cyanobacteria in response to environmental stress is of interest in research on stress tolerance of higher plants. In this thesis, the potd gene encoding an polyamine transporter subunit from Synechocystis sp. strain PCC 6803 was characterized for the first time. A homology model built for PotD protein indicated that it has capability of binding polyamines, with the preference for spermidine. Furthermore, in order to investigate the structural features of the substrate specificity, polyamines were docked into the binding site. Spermidine was positioned very similarly in Synechocystis PotD as in the template structure and had most favorable interactions of the docked polyamines. Based on the homology model, experimental work was conducted, which confirmed the binding preference. Flavodiiron proteins (Flv) are enzymes, which protect the cell against toxicity of oxygen and/or nitric oxide by reduction. In this thesis, we present a novel type of photoprotection mechanism in cyanobacteria by the heterodimer of Flv2/Flv4. The constructed homology model of Flv2/Flv4 suggests a functional heterodimer capable of rapid electron transfer. The unknown protein sll0218, encoded by the flv2-flv4 operon, is assumed to facilitate the interaction of the Flv2/Flv4 heterodimer and energy transfer between the phycobilisome and PSII. Flv2/Flv4 provides an alternative electron transfer pathway and functions as an electron sink in PSII electron transfer.

Probabilistic contact preferences in protein-ligand and protein-protein complexes

Modeller för intermolekulär växelvärkan utnyttjas brett inom biologin. Analys av kontakter mellan proteiner och läkemedelsforskning representerar typiska tillämpningsområden för dylika modeller. En modell som beskriver sådana molekylära växelverkningar kan utformas med hjälp av biofysisk teori, vilket tenderar att resultera i ytterst tung beräkningsbörda även för enkla tillämpningar. Ett alternativt sätt att formulera modeller är att utnyttja stora databaser som innehåller strukturmätningar gjorda med hjälp av till exempel röntgendiffraktion. Då man använder sig av empiriska mätdata direkt, möjliggör en statistisk modell att osäkerheten och inexaktheten i datat tas till hänsyn på ett adekvat sätt, samtidigt som beräkningsbördan håller sig på en rimligare nivå jämfört med kvantmekaniska metoder som i princip borde ge de optimala resultaten. I avhandlingen utvecklades en 3D modell för numerisk undersökning av intermolekulär växelverkan baserad på Bayesiansk statistik. Modellens syfte är att åstadkomma prognoser för det hurdana eller vilka molekylstrukturer prefereras i en given kontext, d.v.s. är mer sannolika inom ramen för interaktion. Modellen testades i essentiella molekyläromgivningar - en liten molekyl vid sin bindningsplats hos ett protein och en gränsyta mellan proteinerna i ett komplex. De erhållna numeriska resultaten motsvarar väl experimentella resultat som tidigare rapporterats i litteraturen, exempelvis kvalitativa bindningsaffiniteter och kemisk kännedom av vissa aminosyrors rumsliga förmågor att utgöra bindningar. I avhandlingen gjordes ytterligare preliminära tester av den statistiska ansatsen för modellering av den centrala molekylära strukturella anpassningsbarheten. I praktiken är den utvecklade modellen ämnad som ett led i en mer omfattande analysmetod, så som en s.k. farmakofor modell. Molekyylivuorovaikutusten mallintamista hyödynnetään laajasti biologisten kysymysten tarkastelussa. Tyypillisiä esimerkkejä sovelluskohteista ovat proteiinien väliset kontaktit ja lääkesuunnittelu. Vuorovaikutuksia kuvaavan mallin lähtökohta voi olla molekyyleihin liittyvä teoria, jolloin soveltamiseen liittyvä laskenta saattaa olla erityisen raskasta, tai suuri havaintojoukko joka on saatu aikaan esimerkiksi mittaamalla rakenteita röntgendiffraktio menetelmällä. Tilastollinen malli mahdollistaa havaintoaineistossa olevan epätarkkuuden ja epävarmuuden huomioimisen, samalla pitäen laskennallisen kuorman pienempänä verrattuna periaatteessa parhaan tuloksen antavaan kvanttimekaaniseen mallinnukseen. Väitöstyössä kehitettiin bayesiläiseen tilastotieteeseen perustuva 3D malli molekyylien välisten vuorovaikutusten laskennalliseen tarkasteluun. Mallin tehtävä on tuottaa ennusteita sen suhteen, minkä tai millaisten molekyylirakenteiden väliset kompleksit ovat etusijalla, toisin sanoen todennäköisempiä, vuorovaikutustilanteessa. Työssä kehitetyn menetelmän toimivuutta testattiin käyttötarkoituksen suhteen olennaisissa molekyyliympäristöissä - pieni molekyyli sitoutumiskohdassaan proteiinissa sekä rajapinta kahden proteiinin välilllä proteiinikompleksissa. Saadut laskennalliset tulokset vastasivat hyvin vertailuun käytettyjä kirjallisuudesta saatuja kokeellisia tuloksia, kuten laadullisia sitoutumisaffiniteetteja, sekä kemiallista tietoa esimerkiksi tiettyjen aminohappojen avaruudellisesta sidoksenmuodostuksesta. Väitöstyössä myös alustavasti testattiin tilastollista lähestymistapaa tärkeän molekyylien rakenteellisen mukautuvuuden mallintamiseen. Käytännössä malli on tarkoitettu osaksi jotakin laajempaa analyysimenetelmää, kuten farmakoforimallia.