EADock : design of a new molecular docking algorithm and some of its applications

3 Summary 3. 1 English The pharmaceutical industry has been facing several challenges during the last years, and the optimization of their drug discovery pipeline is believed to be the only viable solution. High-throughput techniques do participate actively to this optimization, especially when complemented by computational approaches aiming at rationalizing the enormous amount of information that they can produce. In siiico techniques, such as virtual screening or rational drug design, are now routinely used to guide drug discovery. Both heavily rely on the prediction of the molecular interaction (docking) occurring between drug-like molecules and a therapeutically relevant target. Several softwares are available to this end, but despite the very promising picture drawn in most benchmarks, they still hold several hidden weaknesses. As pointed out in several recent reviews, the docking problem is far from being solved, and there is now a need for methods able to identify binding modes with a high accuracy, which is essential to reliably compute the binding free energy of the ligand. This quantity is directly linked to its affinity and can be related to its biological activity. Accurate docking algorithms are thus critical for both the discovery and the rational optimization of new drugs. In this thesis, a new docking software aiming at this goal is presented, EADock. It uses a hybrid evolutionary algorithm with two fitness functions, in combination with a sophisticated management of the diversity. EADock is interfaced with .the CHARMM package for energy calculations and coordinate handling. A validation was carried out on 37 crystallized protein-ligand complexes featuring 11 different proteins. The search space was defined as a sphere of 15 R around the center of mass of the ligand position in the crystal structure, and conversely to other benchmarks, our algorithms was fed with optimized ligand positions up to 10 A root mean square deviation 2MSD) from the crystal structure. This validation illustrates the efficiency of our sampling heuristic, as correct binding modes, defined by a RMSD to the crystal structure lower than 2 A, were identified and ranked first for 68% of the complexes. The success rate increases to 78% when considering the five best-ranked clusters, and 92% when all clusters present in the last generation are taken into account. Most failures in this benchmark could be explained by the presence of crystal contacts in the experimental structure. EADock has been used to understand molecular interactions involved in the regulation of the Na,K ATPase, and in the activation of the nuclear hormone peroxisome proliferatoractivated receptors a (PPARa). It also helped to understand the action of common pollutants (phthalates) on PPARy, and the impact of biotransformations of the anticancer drug Imatinib (Gleevec®) on its binding mode to the Bcr-Abl tyrosine kinase. Finally, a fragment-based rational drug design approach using EADock was developed, and led to the successful design of new peptidic ligands for the a5ß1 integrin, and for the human PPARa. In both cases, the designed peptides presented activities comparable to that of well-established ligands such as the anticancer drug Cilengitide and Wy14,643, respectively. 3.2 French Les récentes difficultés de l'industrie pharmaceutique ne semblent pouvoir se résoudre que par l'optimisation de leur processus de développement de médicaments. Cette dernière implique de plus en plus. de techniques dites "haut-débit", particulièrement efficaces lorsqu'elles sont couplées aux outils informatiques permettant de gérer la masse de données produite. Désormais, les approches in silico telles que le criblage virtuel ou la conception rationnelle de nouvelles molécules sont utilisées couramment. Toutes deux reposent sur la capacité à prédire les détails de l'interaction moléculaire entre une molécule ressemblant à un principe actif (PA) et une protéine cible ayant un intérêt thérapeutique. Les comparatifs de logiciels s'attaquant à cette prédiction sont flatteurs, mais plusieurs problèmes subsistent. La littérature récente tend à remettre en cause leur fiabilité, affirmant l'émergence .d'un besoin pour des approches plus précises du mode d'interaction. Cette précision est essentielle au calcul de l'énergie libre de liaison, qui est directement liée à l'affinité du PA potentiel pour la protéine cible, et indirectement liée à son activité biologique. Une prédiction précise est d'une importance toute particulière pour la découverte et l'optimisation de nouvelles molécules actives. Cette thèse présente un nouveau logiciel, EADock, mettant en avant une telle précision. Cet algorithme évolutionnaire hybride utilise deux pressions de sélections, combinées à une gestion de la diversité sophistiquée. EADock repose sur CHARMM pour les calculs d'énergie et la gestion des coordonnées atomiques. Sa validation a été effectuée sur 37 complexes protéine-ligand cristallisés, incluant 11 protéines différentes. L'espace de recherche a été étendu à une sphère de 151 de rayon autour du centre de masse du ligand cristallisé, et contrairement aux comparatifs habituels, l'algorithme est parti de solutions optimisées présentant un RMSD jusqu'à 10 R par rapport à la structure cristalline. Cette validation a permis de mettre en évidence l'efficacité de notre heuristique de recherche car des modes d'interactions présentant un RMSD inférieur à 2 R par rapport à la structure cristalline ont été classés premier pour 68% des complexes. Lorsque les cinq meilleures solutions sont prises en compte, le taux de succès grimpe à 78%, et 92% lorsque la totalité de la dernière génération est prise en compte. La plupart des erreurs de prédiction sont imputables à la présence de contacts cristallins. Depuis, EADock a été utilisé pour comprendre les mécanismes moléculaires impliqués dans la régulation de la Na,K ATPase et dans l'activation du peroxisome proliferatoractivated receptor a (PPARa). Il a également permis de décrire l'interaction de polluants couramment rencontrés sur PPARy, ainsi que l'influence de la métabolisation de l'Imatinib (PA anticancéreux) sur la fixation à la kinase Bcr-Abl. Une approche basée sur la prédiction des interactions de fragments moléculaires avec protéine cible est également proposée. Elle a permis la découverte de nouveaux ligands peptidiques de PPARa et de l'intégrine a5ß1. Dans les deux cas, l'activité de ces nouveaux peptides est comparable à celles de ligands bien établis, comme le Wy14,643 pour le premier, et le Cilengitide (PA anticancéreux) pour la seconde.

Effects of verapamil and ryanodine on activity of the embryonic chick heart during anoxia and reoxygenation.

Perturbations of the trans-sarcolemmal and sarcoplasmic Ca2+ transport contribute to the abnormal myocardial activity provoked by anoxia and reoxygenation. Whether Ca2+ pools of the extracellular compartment and sarcoplasmic reticulum (SR) are involved to the same extent in the dysfunction of the anoxic-reoxygenated immature heart has not been investigated. Spontaneously contracting hearts isolated from 4-day-old chick embryos were submitted to repeated anoxia (1 min) followed by reoxygenation (5 min). Heart rate, atrioventricular propagation velocity, ventricular shortening, velocities of contraction and relaxation, and incidence of arrhythmias were studied, recorded continuously. Addition of verapamil (10 nM), which blocks selectively sarcolemmal L-type Ca2+ channels, was expected to protect against excessive entry of extracellular Ca2+, whereas addition of ryanodine (10 nM), which opens the SR Ca2+ release channel, was expected to increase cytosolic Ca2+ concentration. Verapamil (a) had no dromotropic effect by contrast to adult heart, (b) attenuated ventricular contracture induced by repeated anoxia, (c) shortened cardioplegia induced by reoxygenation, and (d) had remarkable antiarrhythmic properties during reoxygenation specially. On the other hand, ryanodine potentiated markedly arrhythmias both during anoxia and at reoxygenation. Thus despite its immaturity, the SR seems to be functional early in the developing chick heart and involved in the reversible dysfunction induced by anoxia-reoxygenation. Moreover, Ca2+ entry through L-type channels appears to worsen arrhythmias especially during reoxygenation. These findings show that the Ca2+-handling systems involved in irregular activity in immature heart, such as the embryonic chick heart, may differ from those in the adult.

Blockade of the renin-angiotensin system and renal tissue oxygenation as measured with BOLD-MRI in patients with type 2 diabetes.

AIM: To assess whether blockade of the renin-angiotensin system (RAS), a recognized strategy to prevent the progression of diabetic nephropathy, affects renal tissue oxygenation in type 2 diabetes mellitus (T2DM) patients. METHODS: Prospective randomized 2-way cross over study; T2DM patients with (micro)albuminuria and/or hypertension underwent blood oxygenation level-dependent magnetic resonance imaging (BOLD-MRI) at baseline, after one month of enalapril (20mgqd), and after one month of candesartan (16mgqd). Each BOLD-MRI was performed before and after the administration of furosemide. The mean R2* (=1/T2*) values in the medulla and cortex were calculated, a low R2* indicating high tissue oxygenation. RESULTS: Twelve patients (mean age: 60±11 years, eGFR: 62±22ml/min/1.73m(2)) completed the study. Neither chronic enalapril nor candesartan intake modified renal cortical or medullary R2* levels. Furosemide significantly decreased cortical and medullary R2* levels suggesting a transient increase in renal oxygenation. Medullary R2* levels correlated positively with urinary sodium excretion and systemic blood pressure, suggesting lower renal oxygenation at higher dietary sodium intake and blood pressure; cortical R2* levels correlated positively with glycemia and HbA1c. CONCLUSION: RAS blockade does not seem to increase renal tissue oxygenation in T2DM hypertensive patients. The response to furosemide and the association with 24h urinary sodium excretion emphasize the crucial role of renal sodium handling as one of the main determinants of renal tissue oxygenation.

Interplay of oxidative, nitrosative/nitrative stress, inflammation, cell death and autophagy in diabetic cardiomyopathy.

Diabetes is a recognized risk factor for cardiovascular diseases and heart failure. Diabetic cardiovascular dysfunction also underscores the development of diabetic retinopathy, nephropathy and neuropathy. Despite the broad availability of antidiabetic therapy, glycemic control still remains a major challenge in the management of diabetic patients. Hyperglycemia triggers formation of advanced glycosylation end products (AGEs), activates protein kinase C, enhances polyol pathway, glucose autoxidation, which coupled with elevated levels of free fatty acids, and leptin have been implicated in increased generation of superoxide anion by mitochondria, NADPH oxidases and xanthine oxidoreductase in diabetic vasculature and myocardium. Superoxide anion interacts with nitric oxide forming the potent toxin peroxynitrite via diffusion limited reaction, which in concert with other oxidants triggers activation of stress kinases, endoplasmic reticulum stress, mitochondrial and poly(ADP-ribose) polymerase 1-dependent cell death, dysregulates autophagy/mitophagy, inactivates key proteins involved in myocardial calcium handling/contractility and antioxidant defense, activates matrix metalloproteinases and redox-dependent pro-inflammatory transcription factors (e.g. nuclear factor kappaB) promoting inflammation, AGEs formation, eventually culminating in myocardial dysfunction, remodeling and heart failure. Understanding the complex interplay of oxidative/nitrosative stress with pro-inflammatory, metabolic and cell death pathways is critical to devise novel targeted therapies for diabetic cardiomyopathy, which will be overviewed in this brief synopsis. This article is part of a Special Issue entitled: Autophagy and protein quality control in cardiometabolic diseases.

CR1Dmod: A Matlab program to model 1D complex resistivity effects in electrical and electromagnetic surveys

We have constructed a forward modelling code in Matlab, capable of handling several commonly used electrical and electromagnetic methods in a 1D environment. We review the implemented electromagnetic field equations for grounded wires, frequency and transient soundings and present new solutions in the case of a non-magnetic first layer. The CR1Dmod code evaluates the Hankel transforms occurring in the field equations using either the Fast Hankel Transform based on digital filter theory, or a numerical integration scheme applied between the zeros of the Bessel function. A graphical user interface allows easy construction of 1D models and control of the parameters. Modelling results are in agreement with other authors, but the time of computation is less efficient than other available codes. Nevertheless, the CR1Dmod routine handles complex resistivities and offers solutions based on the full EM-equations as well as the quasi-static approximation. Thus, modelling of effects based on changes in the magnetic permeability and the permittivity is also possible.

Role of CLOCK and circadian rhythms in calcium homeostasis

Le maintien d'une concentration sanguine constante de calcium est d'une importance cruciale et trois organes participent à la balance calcique normale : les reins, les intestins et les os. La concentration plasmatique de calcium est strictement régulée par l'hormone parathyroïdienne (PTH) et par la vitamine D. Des variations circadiennes de la PTH, de la vitamine D ainsi que du calcium plasmatique ont été décrites précédemment chez l'humain ainsi que chez le rat. Ces rythmes de PTH dans le sérum sont importants pour la régulation du remodelage de l'os. En effet, il a été montré chez les souris C57BL/6J que des injections de PTH une fois par jour mènent à une augmentation de la densité minérale de l'os alors que l'infusion en continu de PTH est associée à une diminution de cette densité. La vitamine D joue également un rôle fondamental dans la physiologie osseuse, car un déficit en vitamine D peut conduire à une ostéomalacie. Cependant la fonction des oscillations de vitamine D au niveau de l'homéostasie osseuse reste inconnue. L'horloge circadienne est un système interne de contrôle biologique du temps générant des rythmes de 24 heures dans l'expression des gènes, ainsi que dans la physiologie et le comportement. Ce contrôle s'opère par des boucles rétroactives positives et négatives de l'expression de gènes circadiens tels que CLOCK, BMAL1, CRY1 et 2 ou PERI et 2. Dans ce travail, nous avons émis l'hypothèse que l'homéostasie calcique est sous le contrôle de l'horloge circadienne. Dans un premier temps, nous avons montré chez les souris C57BL/6J des variations journalières des concentrations de calcium, de PTH et de vitamine D dans le sang, ainsi que de calcium dans les urines. Nous avons également démontré des changements au niveau de l'expression rénale des gènes importants dans l'homéostasie du calcium, tant au niveau de l'ARN messager que des protéines. Ensuite, pour analyser le rôle du système de l'horloge circadienne dans l'homéostasie du calcium, nous avons étudié des souris dans lesquelles a été supprimé le gène CLOCK crucial pour la fonction de l'horloge et nous avons comparé ces souris à des souris de type sauvage de même portée. Les souris CLOCK-I- étaient hypercalciuriques à chaque moment de la journée. Cependant le rythme circadien de l'excrétion de calcium était préservé. Le taux de calcium plasmatique ne différait pas entre les génotypes, mais les souris CLOCK -/- ne montraient pas de variations journalières de ce paramètre. Une perte du rythme journalier était également observée pour les niveaux de vitamine D, perte qui pourrait être une cause de l'altération de la micro-architecture osseuse révélée chez les souris CLOCK-/-. En effet, ces souris montrent une diminution du nombre de trabécules, de leur volume ainsi que de leur surface, ce qui suggère la présence d'ostéoporose. Nous avons également trouvé que le rythme de l'expression de l'ARN messager de CYP27B1 était aboli dans les reins des souris CLOCK -/-, ce qui peut expliquer l'altération du rythme de la vitamine D. Les taux sanguins de PTH étaient comparables entre les souris CLOCK -/- et de type sauvage. Dans les reins, une augmentation de l'expression de l'ARN messager de TRPV5 et NCX1 a été constatée, ce qui suggérerait une augmentation de la réabsorption de calcium dans le tubule convoluté distal et dans le tubule connecteur. Dans les intestins, la réabsorption calcique était diminuée, chez les souris CLOCK-I-, fait confirmé par une diminution des niveaux d'ARN messager de TRPV6 et PMCAL. En résumé, la suppression du gène CLOCK chez les souris a conduit à une hypercalciurie, une altération du rythme des taux plasmatiques de calcium et de vitamine D et à une détérioration de l'architecture osseuse. Pour conclure, ces résultats montrent que l'horloge circadienne est essentielle à l'homéostasie calcique ainsi qu'à la physiologie des os. - L'ostéoporose affecte environ 22 millions de femmes et 5.5 millions d'hommes en Europe, réduisant significativement leur qualité de vie et a causé 3.5 millions de nouvelles fractures en 2010. Les dépenses totales liées à ces fractures ont atteint 37 milliards d'euro et ce coût devrait augmenter de 25% d'ici à 2025. Le nombre de nouvelles fractures dues à l'ostéoporose à travers le monde est estimé à environ 1000 par heure. Parmi les causes de l'ostéoporose, le déficit én calcium et/ou en vitamine D joue un rôle important, mais il existe également des causes génétiques ou liées à des facteurs comme les hormones sexuelles (estrogènes, testostérone), l'âge, le tabac, le poids corporel, certains médicaments,... La vie est rythmique : ceci est dû à l'alternance naturelle du jour et de la nuit et de ses effets sur le corps. La prise alimentaire, par exemple, est un processus qui a lieu pendant la phase active, qui est prévisible (il se produit toujours au même moment) et qui peut être anticipé par le corps. Pour cela, une horloge interne est présente dans chaque cellule du corps et est synchronisée par la lumière du jour, entre autres stimuli. Cette horloge indique la phase du jour et régule l'expression de gènes impliqués dans les différents processus qui nécessitent une anticipation. Pendant mon travail de thèse, je me suis demandé si des îythmes circadiens (c'est-à-dire d'une durée d'environ 24 heures et indépendants des stimuli externes) étaient observables'pour les gènes régulant les flux de calcium dans le corps et si l'interruption de ces rythmes pouvait mener à des altérations de la qualité de l'os. J'ai d'abord travaillé avec des souris normales et j'ai pu montrer la présence de rythmes au niveau du calcium sanguin et urinaire, mais également au niveau des hormones et gènes qui contrôlent le métabolisme du calcium dans le corps, comme la vitamine D et l'hormone parathyroidienne. De manière intéressante, j'ai observé que la plupart de ces gènes ont un rythme synchronisé. J'ai ensuite utilisé un modèle de souris dans lequel l'horloge interne a été génétiquement invalidée et j'ai montré que ces souris présentent une augmentation de leur excrétion urinaire de calcium et un rythme circadien altéré de la vitamine D dans le sang. Ces souris absorbent aussi moins bien le calcium intestinal et présentent une ostéoporose marquée. Ce travail montre donc que l'horloge interne est nécessaire pour établir un rythme circadiens de certains facteurs influant les flux de calcium dans l'organisme, comme la vitamine D, et que la perturbation de ces rythmes mène à une dérégulation du métabolisme osseux. Ainsi, la perturbation de l'horloge interne peut causer une ostéoporose et une hypercalciurie qui pourraient aboutir à la formation de fractures et de calculs rénaux. L'extrapolation de ces observations chez l'homme ou à des changements plus subtiles des rythmes circadiens, comme le décalage horaire, restent à montrer. Cette recherche a démontré que les rythmes circadiens des mécanismes de régulation des flux de calcium dans l'organisme sont essentiels au maintien d'un squelette normal et suggère que les perturbations des rythmes circadiens pourraient être une nouvelle cause de l'ostéoporose. - Maintaining constant calcium concentration in the plasma is of a crucial importance and three organs participate in normal calcium balance - kidney, gut and bone. Plasma calcium concentration is strictly regulated by parathyroid hormone (PTH) and vitamin D. Circadian variations of PTH, vitamin D and plasma calcium were previously described in humans, as well as in rats. Rhythms in serum PTH are important for balanced bone remodelling. Indeed in C57BL/6J mice, PTH injection once per day leads to an increase in bone mineral density (BMD), whilst continuous infusion is associated with decreased BMD. Vitamin D also plays a crucial role in bone physiology, since the deficiency in vitamin D can lead to rickets/osteomalacia. However, the role of vitamin D rhythms in bone homeostasis remains unknown. The circadian clock is an. internal time-keeping system generating rhythms in gene expression with 24h periodicity, in physiology and in behaviour. It is operated by positive- and negative-feedback loops of circadian genes, such as CLOCK, BMAL1, CRY1 and 2 or PERI and 2. In this work, we hypothesized, that calcium homeostasis is under the control of the circadian clock. First, we showed daily variations in urinary calcium and serum calcium, PTH and l,25(OH)2 vitamin D, together with renal mRNA and protein levels of genes involved in calcium homeostasis in C57BL/6J mice. Second, and to investigate the role of the circadian clock system in calcium handling, we studied mice lacking the gene CLOCK crucial for fonction of the clock system and compared them to the WT littermates. CLOCK-/- mice were hypercalciuric at all timepoints of the day. However, the circadian rhythm of calcium excretion was preserved. Serum calcium levels did not differ between the genotypes, but CLOCK-/- mice did not exhibit daily variation for this parameter. Loss of rhythm was observed also for serum l,25(OH)2 vitamin D levels, which may be one of the causes of altered bone microarchitecture that was revealed in CLOCK-/- mice. They displayed increased trabecular separation and decreased trabecular number, trabecular bone volume and trabecular bone surface, suggestive of osteoporosis. We found that the rhythm of the mRNA expression of CYP27B1 was abolished in the kidney of CLOCK-/- mice, which could induce the altered rhythm of l,25(OH)2 vitamin. Serum PTH levels were comparable between CLOCK-/- and WT mice. In the kidney, increased mRNA expression of TRPV5 and NCX1 suggests increased calcium reabsorption in the distal convoluted and connecting tubule. In the gut, intestinal calcium absorption was decreased in CLOCK¬/- mice, confirmed by decreased mRNA levels of TRPV6 and PMCA1. In summary, deletion of the CLOCK gene in mice conducts to hypercalciuria, alteration of the rhythm in serum calcium and l,25(OH)2D levels, and impainnent of their bone microarchitecture. In conclusion, these data show that the circadian clock system is essential in calcium homeostasis and bone physiology.

Detecting, quantifying and adjusting for publication bias in meta-analyses: protocol of a systematic review on methods.

BACKGROUND: Health professionals and policymakers aspire to make healthcare decisions based on the entire relevant research evidence. This, however, can rarely be achieved because a considerable amount of research findings are not published, especially in case of 'negative' results - a phenomenon widely recognized as publication bias. Different methods of detecting, quantifying and adjusting for publication bias in meta-analyses have been described in the literature, such as graphical approaches and formal statistical tests to detect publication bias, and statistical approaches to modify effect sizes to adjust a pooled estimate when the presence of publication bias is suspected. An up-to-date systematic review of the existing methods is lacking. METHODS/DESIGN: The objectives of this systematic review are as follows:âeuro¢ To systematically review methodological articles which focus on non-publication of studies and to describe methods of detecting and/or quantifying and/or adjusting for publication bias in meta-analyses.âeuro¢ To appraise strengths and weaknesses of methods, the resources they require, and the conditions under which the method could be used, based on findings of included studies.We will systematically search Web of Science, Medline, and the Cochrane Library for methodological articles that describe at least one method of detecting and/or quantifying and/or adjusting for publication bias in meta-analyses. A dedicated data extraction form is developed and pilot-tested. Working in teams of two, we will independently extract relevant information from each eligible article. As this will be a qualitative systematic review, data reporting will involve a descriptive summary. DISCUSSION: Results are expected to be publicly available in mid 2013. This systematic review together with the results of other systematic reviews of the OPEN project (To Overcome Failure to Publish Negative Findings) will serve as a basis for the development of future policies and guidelines regarding the assessment and handling of publication bias in meta-analyses.

Eugenol as an anesthetic for juvenile common snook

The objective of this work was to evaluate the efficacy of eugenol as an anesthetic for juvenile common snook, and to determine the minimum effective concentration for use in handling procedures. In the first trial, juvenile common snook were subjected to immersion baths at 25, 50, 75, 100, 125, and 150 mg L-1 eugenol concentrations, after which induction and recovery times were evaluated. In the second experiment, the lethal exposure time (LT50) at 75 mg L-1 was estimated. Minimum effective eugenol concentration was 50 mg L-1, andthe stage of deep anesthesia and recovery were, respectively, reached at 126.3 and 208.8 s. At 75 mg L-1, LT50 was 1,314 s, and induction time and recovery were also satisfactory; however, fish cannot tolerate over 229 s exposure.

EasyStrata: evaluation and visualization of stratified genome-wide association meta-analysis data.

The R package EasyStrata facilitates the evaluation and visualization of stratified genome-wide association meta-analyses (GWAMAs) results. It provides (i) statistical methods to test and account for between-strata difference as a means to tackle gene-strata interaction effects and (ii) extended graphical features tailored for stratified GWAMA results. The software provides further features also suitable for general GWAMAs including functions to annotate, exclude or highlight specific loci in plots or to extract independent subsets of loci from genome-wide datasets. It is freely available and includes a user-friendly scripting interface that simplifies data handling and allows for combining statistical and graphical functions in a flexible fashion. AVAILABILITY: EasyStrata is available for free (under the GNU General Public License v3) from our Web site www.genepi-regensburg.de/easystrata and from the CRAN R package repository cran.r-project.org/web/packages/EasyStrata/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

New directions in Stochastic Multicriteria Acceptability Analysis

Decisions taken in modern organizations are often multi-dimensional, involving multiple decision makers and several criteria measured on different scales. Multiple Criteria Decision Making (MCDM) methods are designed to analyze and to give recommendations in this kind of situations. Among the numerous MCDM methods, two large families of methods are the multi-attribute utility theory based methods and the outranking methods. Traditionally both method families require exact values for technical parameters and criteria measurements, as well as for preferences expressed as weights. Often it is hard, if not impossible, to obtain exact values. Stochastic Multicriteria Acceptability Analysis (SMAA) is a family of methods designed to help in this type of situations where exact values are not available. Different variants of SMAA allow handling all types of MCDM problems. They support defining the model through uncertain, imprecise, or completely missing values. The methods are based on simulation that is applied to obtain descriptive indices characterizing the problem. In this thesis we present new advances in the SMAA methodology. We present and analyze algorithms for the SMAA-2 method and its extension to handle ordinal preferences. We then present an application of SMAA-2 to an area where MCDM models have not been applied before: planning elevator groups for high-rise buildings. Following this, we introduce two new methods to the family: SMAA-TRI that extends ELECTRE TRI for sorting problems with uncertain parameter values, and SMAA-III that extends ELECTRE III in a similar way. An efficient software implementing these two methods has been developed in conjunction with this work, and is briefly presented in this thesis. The thesis is closed with a comprehensive survey of SMAA methodology including a definition of a unified framework.

Materiaalitietokanta prosessilaitetoimituksen laitevalintojen tukena

In this master's thesis, material characteristics which includes material database will be developed, through which will be managed process device deliveries device selections. The analysed bulk materials exact values serve as basic data when choosing glass- and dried product industry raw material equipment solutions. In this study, the material database function structure will be divided into classified bulk materials characteristics and process device guidelines and planning which they are based on. In order to achieve the best possible advantage in use of the material database, the raw material handling with the process equipment has to be controlled with the help of database bulk material information. With the help of material database, the personnel of the company will find the needed material and device information to implement the processes. The structure of the database enables adding of information, and hence the database will be always kept updated. The material database uses free PhpOpenDatabase software programme, which adapts itself flexibly with the software and requirements of the company. The material frame and considered features of the material database were determined by the requirements of the company's device technology. National and international classifications were used as help. The collected material information was analysed and the correspondence between internal study results and found classifications were determined. As presented before, on the grounds of chosen and specified material features, also the controllability of device selecting and solutions was improved with device-specific screening and planning guidelines. As material frame of the material database was chosen ANSI/CEMA Standard 550. Its almost 1000 bulk materials form an extensive database frame. Over 600 of these raw materials have been determined according to classification. However, the most important ones are the materials of glass- and dried material plants, which consist of the last part of material frame. The analysis will be sustained and material information will be completed primarily of material part with test run in one's own material laboratory.

Naudan lietelannan käsittelymenetelmien taloudellinen vertailu

Maataloudessa syntyvä lanta on arvokas lannoite ja maanparannusaine, jonka käsittelystä aiheutuu sekä kustannuksia että ympäristövaikutuksia. Muita haasteita ovat esimerkiksi lannan mikrobit, lannan levitykseen soveltuvan ajankohdan lyhyys, lannan ravinteiden sovittaminen kasvien tarpeisiin ja lannan ravinteiden määrä suhteessa levityskelpoisen peltoalan määrään. Tutkimuksen tavoitteena oli tunnistaa naudan lietelannan käsittelyketjuihin liittyvät kustannukset ja osoittaa eri käsittelyketjujen kustannusten eroavaisuudet. Tavoitteena oli myös tunnistaa ja osoittaa käsittelyketjujen laadulliset erot. Kustannukset selvitettiin kustannuslaskelmin ja laadulliset erot SWOT-menetelmällä. Tutkimuksen kohteeksi valittiin 6 tilakokoluokkaa ja käsittelyketjuiksi lietelanta-, kompostointi- ja mädätysketju. Tutkimuksessa alhaisimmat kustannukset olivat lietelantaketjulla, jonka kustannukset 25 – 250 naudan tilalla olivat 5 200 – 6 600 €/a ja yksikkökustannukset 1 – 9 €/m3. Mädätysketjun kustannukset vastaavissa tilakokoluokissa olivat noin 33 000 – 50 000 €/a ja yksikkökustannukset 8 – 55 €/m3. Kompostointiketjun kustannukset olivat 35 000 – 143 000 €/a ja yksikkökustannukset 24 – 58 €/m3. Lietelantaketjun edullisuus johtui vähäisistä laite- ja rakennusinvestoinneista ja pienistä työmääristä ja kompostointiketjun kalleus suurista tukiaine- ja investointikustannuksista. Käsittelyketjujen asettaminen paremmuusjärjestykseen oli hankalaa. Työn määrä oli pienin lietelantaketjussa ja toiseksi pienin suurilla tiloilla mädätysketjussa. Kompostointiketjulla itse levitykseen kuluva aika oli pienin. Ravinteiden osalta mädätysketju oli parhain ja kompostointiketju huonoin. Ympäristövaikutuksiltaan ja hajuhaitoiltaan kompostointi- ja mädätysketju olivat parhaimmat. Mikrobien osalta parhain oli kompostointiketju.

Analysis of Industrial Granular Flow Applications by Using Advanced Collision Models

Granular flow phenomena are frequently encountered in the design of process and industrial plants in the traditional fields of the chemical, nuclear and oil industries as well as in other activities such as food and materials handling. Multi-phase flow is one important branch of the granular flow. Granular materials have unusual kinds of behavior compared to normal materials, either solids or fluids. Although some of the characteristics are still not well-known yet, one thing is confirmed: the particle-particle interaction plays a key role in the dynamics of granular materials, especially for dense granular materials. At the beginning of this thesis, detailed illustration of developing two models for describing the interaction based on the results of finite-element simulation, dimension analysis and numerical simulation is presented. The first model is used to describing the normal collision of viscoelastic particles. Based on some existent models, more parameters are added to this model, which make the model predict the experimental results more accurately. The second model is used for oblique collision, which include the effects from tangential velocity, angular velocity and surface friction based on Coulomb's law. The theoretical predictions of this model are in agreement with those by finite-element simulation. I n the latter chapters of this thesis, the models are used to predict industrial granular flow and the agreement between the simulations and experiments also shows the validation of the new model. The first case presents the simulation of granular flow passing over a circular obstacle. The simulations successfully predict the existence of a parabolic steady layer and show how the characteristics of the particles, such as coefficients of restitution and surface friction affect the separation results. The second case is a spinning container filled with granular material. Employing the previous models, the simulation could also reproduce experimentally observed phenomena, such as a depression in the center of a high frequency rotation. The third application is about gas-solid mixed flow in a vertically vibrated device. Gas phase motion is added to coherence with the particle motion. The governing equations of the gas phase are solved by using the Large eddy simulation (LES) and particle motion is predicted by using the Lagrangian method. The simulation predicted some pattern formation reported by experiment.

Yhdyskuntalietteen ja biopolttoaineiden käsittelyketjujen ja polton teknistaloudellinen tarkastelu

Työn tavoitteena on kartoittaa yhdyskuntalietteen ja kierrätys- sekä biopolttoaineiden käsittelyä ja polttoa lietteenpolttolaitoksen tarpeita ajatellen. Lietteen käsittelyketjun ja kierrätys- sekä biopolttoaineketjujen tekninen tarkastelu on siis työn keskeinen tavoite. Lisäksi lasketaan polttolaitoksen suurimpia mahdollisia investointikustannuksia eri polttoainevaihtoehdoilla. Työssä tehdään muun ohella case-tarkastelua Kaakkois- Suomen alueeseen liittyen. Tavoitteena on muodostaa tarkoitukseen soveltuva polttoaineratkaisu kullekin tapauk-selle. Työn alkuosassa tutustutaan yleisesti lietteeseen sekä polttoaineen että jätteen roolissa. Tarkastelu sisältää tietoja lietteen ominaisuuksista sekä lietteenkäsittelyssä olennaisista lainsäädännöllisistä seikoista. Samoin katsastetaan hieman lietteen esikäsittelyä, mekaanista vedenerotusta, termistä kuivausta ja polttoa tarkastellaan yleisessä valossa. Lisäksi alkuosassa keskitytään eri bio- jakierrätyspolttoainevaihtoehtoihin tarkastelemalla niiden yleisyyttä polttoaineena sekä esittelemällä niiden käsittelyketjuja. Työn loppupuoliskolla kiinnitetään huomiota case tapausten avulla polttolaitoksesta saataviin tuottoihin sekä millaisen liikkumavaran eri polttoainevaihtoehdot investointien osalta sallivat. Case-tapauksissa pohditaan Kymenlaakson ja Etelä-Karjalan paikallisia lietteen-polttomahdollisuuksia yhdistettynä kierrätys- tai biopolttoaineisiin. Mekaanisesti kuivattua lietettä käsitellään kyseisissä tapauksissa vuosittain 6000 t ja 15 000 t. Lietteen polton tuottama sähkö- ja lämpöteho näyttävät riippuvan voimakkaasti lietteen kuiva-ainepitoisuudesta, eivät niinkään lietteen muista ominaisuuksista. Lisäksi joko bio- tai kierrätyspolttoaineella saadaan sähkön- ja lämmöntuotantoa nostettua huomattavasti.