Multimethodological study of molecular recognition phenomena

The study of the bio-recognition phenomena behind a biological process is nowadays considered a useful tool to deeply understand physiological mechanisms allowing the discovery of novel biological target and the development of new lead candidates. Moreover, understanding this kind of phenomena can be helpful in characterizing absorption, distribution, metabolism, elimination and toxicity properties of a new drug (ADMET parameters). Recent estimations show that about half of all drugs in development fail to make it to the market because of ADMET deficiencies; thus a rapid determination of ADMET parameters in early stages of drug discovery would save money and time, allowing to choose the better compound and to eliminate any losers. The monitoring of drug binding to plasma proteins is becoming essential in the field of drug discovery to characterize the drug distribution in human body. Human serum albumin (HSA) is the most abundant protein in plasma playing a fundamental role in the transport of drugs, metabolites and endogenous factors; so the study of the binding mechanism to HSA has become crucial to the early characterization of the pharmacokinetic profile of new potential leads. Furthermore, most of the distribution experiments carried out in vivo are performed on animals. Hence it is interesting to determine the binding of new compounds to albumins from different species to evaluate the reliability of extrapolating the distribution data obtained in animals to humans. It is clear how the characterization of interactions between proteins and drugs determines a growing need of methodologies to study any specific molecular event. A wide variety of biochemical techniques have been applied to this purpose. High-performance liquid affinity chromatography, circular dichroism and optical biosensor represent three techniques that can be able to elucidate the interaction of a new drug with its target and with others proteins that could interfere with ADMET parameters.

Optimization of molecular and crystalline forms of drugs, agrochemicals, pesticides in relation to activity, bioavailability, patentability and to the fabrication of polymorphs, solvates, co-crystals with green chemistry methods

This doctorate was funded by the Regione Emilia Romagna, within a Spinner PhD project coordinated by the University of Parma, and involving the universities of Bologna, Ferrara and Modena. The aim of the project was: - Production of polymorphs, solvates, hydrates and co-crystals of active pharmaceutical ingredients (APIs) and agrochemicals with green chemistry methods; - Optimization of molecular and crystalline forms of APIs and pesticides in relation to activity, bioavailability and patentability. In the last decades, a growing interest in the solid-state properties of drugs in addition to their solution chemistry has blossomed. The achievement of the desired and/or the more stable polymorph during the production process can be a challenge for the industry. The study of crystalline forms could be a valuable step to produce new polymorphs and/or co-crystals with better physical-chemical properties such as solubility, permeability, thermal stability, habit, bulk density, compressibility, friability, hygroscopicity and dissolution rate in order to have potential industrial applications. Selected APIs (active pharmaceutical ingredients) were studied and their relationship between crystal structure and properties investigated, both in the solid state and in solution. Polymorph screening and synthesis of solvates and molecular/ionic co-crystals were performed according to green chemistry principles. Part of this project was developed in collaboration with chemical/pharmaceutical companies such as BASF (Germany) and UCB (Belgium). We focused on on the optimization of conditions and parameters of crystallization processes (additives, concentration, temperature), and on the synthesis and characterization of ionic co-crystals. Moreover, during a four-months research period in the laboratories of Professor Nair Rodriguez-Hormedo (University of Michigan), the stability in aqueous solution at the equilibrium of ionic co-crystals (ICCs) of the API piracetam was investigated, to understand the relationship between their solid-state and solution properties, in view of future design of new crystalline drugs with predefined solid and solution properties.

Systematic functionalization of molecules for molecular self-assembly

This thesis describes the investigation of systematically varied organic molecules for use in molecular self-assembly processes. All experiments were performed using high-resolution non-contact atomic force microscopy under UHV conditions and at room temperature. Using this technique, three different approaches for influencing intermolecular and molecule-surface interaction on the insulating calcite(10.4) surface were investigated by imaging the structure formation at the molecular scale. I first demonstrated the functionalization of shape-persistent oligo(p-benzamide)s that was engineered by introducing different functional groups and investigating their effect on the structural formation on the sample surface. The molecular core was designed to provide significant electrostatic anchoring towards the surface, while at the same time maintaining the flexibility to fine-tune the resulting structure by adjusting the intermolecular cohesion energy. The success of this strategy is based on a clear separation of the molecule-substrate interaction from the molecule-molecule interaction. My results show that sufficient molecule-surface anchoring can be achieved without restricting the structural flexibility that is needed for the design of complex molecular systems. Three derivatives of terephthalic acid (TPA) were investigated in chapter 7. Here, the focus was on changing the adhesion to the calcite surface by introducing different anchor functionalities to the TPA backbone. For all observed molecules, the strong substrate templating effect results in molecular structures that are strictly oriented along the calcite main crystal directions. This templating is especially pronounced in the case of 2-ATPA where chain formation on the calcite surface is observed in contrast to the formation of molecular layers in the bulk. At the same time, the amino group of 2-ATPA proved an efficient anchor functionality, successfully stabilizing the molecular chains on the sample surface. These findings emphasizes, once again, the importance of balancing and fine-tuning molecule-molecule and molecule-surface interactions in order to achieve stable, yet structurally flexible molecular arrangements on the sample surface. In the last chapter, I showed how the intrinsic property of molecular chirality decisively influences the structure formation in molecular self-assembly. This effect is especially pronounced in the case of the chiral heptahelicene-2-carboxylic acid. Deposition of the enantiopure molecules results in the formation of homochiral islands on the sample surface which is in sharp contrast to the formation of uni-directional double rows upon deposition of the racemate onto the same surface. While it remained uncertain from these previous experiments whether the double rows are composed of hetero- or homochiral molecules, I could clearly answer that question here and demonstrate that the rows are of heterochiral origin. Chirality, thus, proves to be another important parameter to steer the intermolecular interaction on surfaces. Altogether, the results of this thesis demonstrate that, in order to successfully control the structure formation in molecular self-assembly, the correct combination of molecule and surface properties is crucial. This is of special importance when working on substrates that exhibit a strong influence on the structure formation, such as the calcite(10.4) surface. Through the systematic variation of functional groups several important parameters that influence the balance between molecule-surface and molecule-molecule interaction were identified here, and the results of this thesis can, thus, act as a guideline for the rational design of molecules for use in molecular self-assembly.

Molecular plating of thin lanthanide layers with improved material properties for nuclear applications

In der vorliegenden Arbeit werden Experimente beschrieben, die zu einem vertieften Verständnis fundamentaler Prozesse bei der elektrochemischen Herstellung von Dünnschichten, sog. Targets, für kernphysikalische und -chemische Studien führten. Targets wurden mittels 'Molecular Plating' (MP) hergestellt, indem eine Elektrodeposition aus organischem Medium in der Regel bei konstantem Strom in Zwei-Elektroden-Zellen. Die Resultate erlaubten, optimierte Herstellungs-bedingungen zu ermitteln, welche die Produktion deutlich verbesserter Targets erlaubten. MP bei konstantem Strom ist ein massentransportkontrollierter Prozess. Der angelegte Strom wird durch einen konstanten Fluss elektroaktiver Spezies zur Kathode – auf der die Schicht wächst – und Anode aufrechterhalten. Die Untersuchungen zeigten, dass das Zellenpotential des Elektrodepositionsystems immer durch den Ohm'schen Spannungsabfall auf Grund des Widerstandes der verwendeten Lösung dominiert wurde. Dies erlaubte die Herleitung einer Beziehung zwischen dem Zellenpotential und der Konzentration der elektroaktiven Spezies. Die Beziehung erlaubt die Erklärung des gemessenen zeitlichen Verlaufs des Zellenpotentials während der Abscheidung als Funktion der Elektrolytkonzentration. Dies dient als Basis, auf der nun ein umfassenderes Bild der Prozesse, die für die charakteristischen Minima im Potentialverlauf einer Abscheidung verantwortlich sind, gewonnen werden kann. Es konnte gezeigt werden, dass die Minima mit der fast vollständigen Entfernung (durch Abscheidung) der aus einem gelösten Salz erzeugten Nd-Ionen korrespondieren. Die abgeschiedene Spezies wurde als Nd3+ identifiziert, vermutlich als Carboxylat, Oxid oder Hydroxid, was auf Grund der hohen negative Werte des Standardredoxpotentials der Lanthanide verständlich erscheint. Von den vorliegenden elektroaktiven Spezies tragen die Nd3+ Ionen nur zu knapp 20% zum Gesamtstrom bei. Durch Elektrolyse tragen auch die Lösungsmittelkomponenten zu diese Strom bei. Die Gegenwart von elektrolysiertem Lösungsmittel wurde in Analysen der Dünnschichten bestätigt. Diese waren immer mit chemi- und physisorbierten Lösungsmittelmolekülen bedeckt. Die Analyse der Dünnschichten zeigte, dass die Oberflächen von einem furchenartiges Netz durchzogen waren, und dass diese während des Trocknen der Schichten nach dem MP entstanden. Ob die Schichten an Luft oder in inerter Atmosphäre trockneten, hatte keinen Einfluss. Es wurden Experimente mit mehreren Lösungsmitteln durchgeführt, die sich deutlich in ihren physikalischen Eigenschaften, v.a. dem Siedepunkt, unterschieden. Furchenfreie Dünnschichten konnten insbesondere bei MP in N,N-dimethylformamide (DMF) erzeugt werden. Die Verwendung von DMF in Kombination mit einer Abscheidung auf sehr glatten Substraten erlaubte die Produktion von sehr homogenen, glatten und defektfreien Schichten. Diese waren vermutlich geringeren inneren Spannungen während des Trocknens ausgesetzt, als Schichten auf raueren Substraten oder solche, die aus flüchtigeren Lösungsmitteln hergestellt wurden. Die Oberflächenrauigkeit des Substrats und das gewählte Lösungsmittel wurden so als Schlüsselfaktoren für die Produktion hochqualitativer Schichten identifiziert. Es konnte gezeigt werden, dass mit MP eine sehr effiziente Methode zur Herstellung homogener Schichten mit exzellenter Ausbeute ist. In weiteren Experimenten mit dem primordialen Alpha-Emitter 147Sm als Modellisotop wurde die Eignung solcher Schichten als Alpha-Quelle untersucht. Sowohl die Energieauflösung als auch der Anteil der Alpha-Teilchen, die den Detektor erreichten, waren von den Quelleneigenschaften abhängig. Die Effekte wurden verschiedenen Variablen der Dünnschicht zugeordnet, welche die Alpha-Spektren beeinflussten. Dominant war die Wahl des Lösungsmittels und die Rauigkeit des Substrats. Dies beeinflusste Schichtdicke und -morphologie sowie die Art des Schichtwachstums und veränderte die Detektionseffizienz in Alpha-Messungen bis zu 15%. Nur homogene, ebene Schichten, die aus DMF auf glatten Substraten abgeschieden wurden, eignen sich optimal als Alpha-Quelle. Die gewonnenen Ergebnisse erlauben die optimierte Herstellung nuklearer Targets durch MP. Künftige Anwendungen beinhalten insbesondere die Herstellung von Targets für neutroneninduzierte Spaltexperimente und untergrundarmeAlpha-Messungen sehr kleiner Aktivitäten.

Novel algorithms for electronic structure based molecular dynamics with linear system-size scaling

Die vorliegende Arbeit behandelt die Entwicklung und Verbesserung von linear skalierenden Algorithmen für Elektronenstruktur basierte Molekulardynamik. Molekulardynamik ist eine Methode zur Computersimulation des komplexen Zusammenspiels zwischen Atomen und Molekülen bei endlicher Temperatur. Ein entscheidender Vorteil dieser Methode ist ihre hohe Genauigkeit und Vorhersagekraft. Allerdings verhindert der Rechenaufwand, welcher grundsätzlich kubisch mit der Anzahl der Atome skaliert, die Anwendung auf große Systeme und lange Zeitskalen. Ausgehend von einem neuen Formalismus, basierend auf dem großkanonischen Potential und einer Faktorisierung der Dichtematrix, wird die Diagonalisierung der entsprechenden Hamiltonmatrix vermieden. Dieser nutzt aus, dass die Hamilton- und die Dichtematrix aufgrund von Lokalisierung dünn besetzt sind. Das reduziert den Rechenaufwand so, dass er linear mit der Systemgröße skaliert. Um seine Effizienz zu demonstrieren, wird der daraus entstehende Algorithmus auf ein System mit flüssigem Methan angewandt, das extremem Druck (etwa 100 GPa) und extremer Temperatur (2000 - 8000 K) ausgesetzt ist. In der Simulation dissoziiert Methan bei Temperaturen oberhalb von 4000 K. Die Bildung von sp²-gebundenem polymerischen Kohlenstoff wird beobachtet. Die Simulationen liefern keinen Hinweis auf die Entstehung von Diamant und wirken sich daher auf die bisherigen Planetenmodelle von Neptun und Uranus aus. Da das Umgehen der Diagonalisierung der Hamiltonmatrix die Inversion von Matrizen mit sich bringt, wird zusätzlich das Problem behandelt, eine (inverse) p-te Wurzel einer gegebenen Matrix zu berechnen. Dies resultiert in einer neuen Formel für symmetrisch positiv definite Matrizen. Sie verallgemeinert die Newton-Schulz Iteration, Altmans Formel für beschränkte und nicht singuläre Operatoren und Newtons Methode zur Berechnung von Nullstellen von Funktionen. Der Nachweis wird erbracht, dass die Konvergenzordnung immer mindestens quadratisch ist und adaptives Anpassen eines Parameters q in allen Fällen zu besseren Ergebnissen führt.

Numerical analysis of the relation between interactions and structure in a molecular fluid

Coarse graining is a popular technique used in physics to speed up the computer simulation of molecular fluids. An essential part of this technique is a method that solves the inverse problem of determining the interaction potential or its parameters from the given structural data. Due to discrepancies between model and reality, the potential is not unique, such that stability of such method and its convergence to a meaningful solution are issues.rnrnIn this work, we investigate empirically whether coarse graining can be improved by applying the theory of inverse problems from applied mathematics. In particular, we use the singular value analysis to reveal the weak interaction parameters, that have a negligible influence on the structure of the fluid and which cause non-uniqueness of the solution. Further, we apply a regularizing Levenberg-Marquardt method, which is stable against the mentioned discrepancies. Then, we compare it to the existing physical methods - the Iterative Boltzmann Inversion and the Inverse Monte Carlo method, which are fast and well adapted to the problem, but sometimes have convergence problems.rnrnFrom analysis of the Iterative Boltzmann Inversion, we elaborate a meaningful approximation of the structure and use it to derive a modification of the Levenberg-Marquardt method. We engage the latter for reconstruction of the interaction parameters from experimental data for liquid argon and nitrogen. We show that the modified method is stable, convergent and fast. Further, the singular value analysis of the structure and its approximation allows to determine the crucial interaction parameters, that is, to simplify the modeling of interactions. Therefore, our results build a rigorous bridge between the inverse problem from physics and the powerful solution tools from mathematics. rn

Different molecular and structural adaptations with eccentric and conventional strength training in elderly men and women

Reprogramming of gene expression contributes to structural and functional adaptation of muscle tissue in response to altered use. The aim of this study was to investigate mechanisms for observed improvements in leg extension strength, gain in relative thigh muscle mass and loss of body and thigh fat content in response to eccentric and conventional strength training in elderly men (n = 14) and women (n = 14; average age of the men and women: 80.1 ± 3.7 years) by means of structural and molecular analyses. Biopsies were collected from m. vastus lateralis in the resting state before and after 12 weeks of training with two weekly resistance exercise sessions (RET) or eccentric ergometer sessions (EET). Gene expression was analyzed using custom-designed low-density PCR arrays. Muscle ultrastructure was evaluated using EM morphometry. Gain in thigh muscle mass was paralleled by an increase in muscle fiber cross-sectional area (hypertrophy) with RET but not with EET, where muscle growth is likely occurring by the addition of sarcomeres in series or by hyperplasia. The expression of transcripts encoding factors involved in muscle growth, repair and remodeling (e.g., IGF-1, HGF, MYOG, MYH3) was increased to a larger extent after EET than RET. MicroRNA 1 expression was decreased independent of the training modality, and was paralleled by an increased expression of IGF-1 representing a potential target. IGF-1 is a potent promoter of muscle growth, and its regulation by microRNA 1 may have contributed to the gain of muscle mass observed in our subjects. EET depressed genes encoding mitochondrial and metabolic transcripts. The changes of several metabolic and mitochondrial transcripts correlated significantly with changes in mitochondrial volume density. Intramyocellular lipid content was decreased after EET concomitantly with total body fat. Changes in intramyocellular lipid content correlated with changes in body fat content with both RET and EET. In the elderly, RET and EET lead to distinct molecular and structural adaptations which might contribute to the observed small quantitative differences in functional tests and body composition parameters. EET seems to be particularly convenient for the elderly with regard to improvements in body composition and strength but at the expense of reducing muscular oxidative capacity.

Interaction between rheumatoid arthritis and pregnancy: correlation of molecular data with clinical disease activity measures

OBJECTIVE: The factors that induce remission of RA during pregnancy and the relapse occurring after delivery remain an enigma. In a previous study, we investigated gene-expression profiles of peripheral blood mononuclear cells (PBMC) in patients with RA and healthy women in late pregnancy and postpartum. Profiles of samples from both groups were similar in late pregnancy with elevated monocyte and decreased lymphocyte signatures. Postpartum, in RA PBMC the high level of monocyte transcripts persisted. Further increase was observed in adhesion, migration and signalling processes related to monocytes but also in lymphocytes despite similar clinical activity due to intensified drug treatment. This prompted us to investigate correlations between clinical parameters of disease activity and gene profiles. METHODS: Transcriptome data were correlated with RADAI, CRP, monocyte and lymphocyte counts. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway annotations, monocytes and lymphocytes signatures were used as reference information. RESULTS: Comparative analysis of PBMC expression profiles from RA patients during and after pregnancy with RADAI and CRP revealed a correlation of these disease activity parameters predominantly with monocyte transcripts. Genes related to cellular programs of adhesion, migration and response to infections were upregulated. Comparing clinically active and not-active RA patients postpartum revealed a cluster of 19 genes that could also identify active disease during pregnancy. CONCLUSION: The data suggest that an increase of the RADAI and an elevation of CRP is a consequence of molecular activation of monocytes. Furthermore, they indicate that molecular activation of T lymphocytes may remain clinically unrecognized postpartum. It is conceivable that a set of 19 genes may qualify as molecular disease activity marker.

Predicting the Mechanical Properties of Carbon-Based Materials Using Molecular Dynamics

The development of innovative carbon-based materials can be greatly facilitated by molecular modeling techniques. Although molecular modeling has been used extensively to predict elastic properties of materials, modeling of more complex phenomenon such as fracture has only recently been possible with the development of new force fields such as ReaxFF, which is used in this work. It is not fully understood what molecular modeling parameters such as thermostat type, thermostat coupling, time step, system size, and strain rate are required for accurate modeling of fracture. Selection of modeling parameters to model fracture can be difficult and non-intuitive compared to modeling elastic properties using traditional force fields, and the errors generated by incorrect parameters may be non-obvious. These molecular modeling parameters are systematically investigated and their effects on the fracture of well-known carbon materials are analyzed. It is determined that for coupling coefficients of 250 fs and greater do not result in substantial differences in the stress-strain response of the materials using any thermostat type. A time step of 0.5 fs of smaller is required for accurate results. Strain rates greater than 2.2 ns-1 are sufficient to obtain repeatable results with slower strain rates for the materials studied. The results of this study indicate that further refinement of the Chenoweth parameter set is required to accurately predict the mechanical response of carbon-based systems. The ReaxFF has been used extensively to model systems in which bond breaking and formation occur. In particular ReaxFF has been used to model reactions of small molecules. Some elastic and fracture properties have been successfully modeled using ReaxFF in materials such as silicon and some metals. However, it is not clear if current parameterizations for ReaxFF are able to accurately reproduce the elastic and fracture properties of carbon materials. The stress-strain response of a new ReaxFF parameterization is compared to the previous parameterization and density functional theory results for well-known carbon materials. The new ReaxFF parameterization makes xv substantial improvements to the predicted mechanical response of carbon materials, and is found to be suitable for modeling the mechanical response of carbon materials. Finally, a new material composed of carbon nanotubes within an amorphous carbon (AC) matrix is modeled using the ReaxFF. Various parameters that may be experimentally controlled are investigated such as nanotube bundling, comparing multi-walled nanotube with single-walled nanotubes, and degree of functionalization of the nanotubes. Elastic and fracture properties are investigated for the composite systems and compared to results of pure-nanotube and pure-AC models. It is found that the arrangement of the nanotubes and degree of crosslinking may substantially affect the properties of the systems, particularly in the transverse directions.

Bifurcation and singularity analysis of a molecular network for the induction of long-term memory.

Withdrawal reflexes of the mollusk Aplysia exhibit sensitization, a simple form of long-term memory (LTM). Sensitization is due, in part, to long-term facilitation (LTF) of sensorimotor neuron synapses. LTF is induced by the modulatory actions of serotonin (5-HT). Pettigrew et al. developed a computational model of the nonlinear intracellular signaling and gene network that underlies the induction of 5-HT-induced LTF. The model simulated empirical observations that repeated applications of 5-HT induce persistent activation of protein kinase A (PKA) and that this persistent activation requires a suprathreshold exposure of 5-HT. This study extends the analysis of the Pettigrew model by applying bifurcation analysis, singularity theory, and numerical simulation. Using singularity theory, classification diagrams of parameter space were constructed, identifying regions with qualitatively different steady-state behaviors. The graphical representation of these regions illustrates the robustness of these regions to changes in model parameters. Because persistent protein kinase A (PKA) activity correlates with Aplysia LTM, the analysis focuses on a positive feedback loop in the model that tends to maintain PKA activity. In this loop, PKA phosphorylates a transcription factor (TF-1), thereby increasing the expression of an ubiquitin hydrolase (Ap-Uch). Ap-Uch then acts to increase PKA activity, closing the loop. This positive feedback loop manifests multiple, coexisting steady states, or multiplicity, which provides a mechanism for a bistable switch in PKA activity. After the removal of 5-HT, the PKA activity either returns to its basal level (reversible switch) or remains at a high level (irreversible switch). Such an irreversible switch might be a mechanism that contributes to the persistence of LTM. The classification diagrams also identify parameters and processes that might be manipulated, perhaps pharmacologically, to enhance the induction of memory. Rational drug design, to affect complex processes such as memory formation, can benefit from this type of analysis.

Microscopic theory of cooperative spin crossover: Interaction of molecular modes with phonons

In this article, we present a new microscopic theoretical approach to the description of spin crossover in molecular crystals. The spin crossover crystals under consideration are composed of molecular fragments formed by the spin-crossover metal ion and its nearest ligand surrounding and exhibiting well defined localized (molecular) vibrations. As distinguished from the previous models of this phenomenon, the developed approach takes into account the interaction of spin-crossover ions not only with the phonons but also a strong coupling of the electronic shells with molecular modes. This leads to an effective coupling of the local modes with phonons which is shown to be responsible for the cooperative spin transition accompanied by the structural reorganization. The transition is characterized by the two order parameters representing the mean values of the products of electronic diagonal matrices and the coordinates of the local modes for the high- and low-spin states of the spin crossover complex. Finally, we demonstrate that the approach provides a reasonable explanation of the observed spin transition in the [Fe(ptz)6](BF4)2 crystal. The theory well reproduces the observed abrupt low-spin → high-spin transition and the temperature dependence of the high-spin fraction in a wide temperature range as well as the pronounced hysteresis loop. At the same time within the limiting approximations adopted in the developed model, the evaluated high-spin fraction vs. T shows that the cooperative spin-lattice transition proves to be incomplete in the sense that the high-spin fraction does not reach its maximum value at high temperature.

Seasonal variation of high-molecular-weight compounds in the water-soluble fraction of organic urban aerosols

Aerosol samples were collected in Zurich, Switzerland, at an urban background site and were analyzed with size exclusion chromatography (SEC) and laser/desorption ionization mass spectrometry (LDI-MS) for water-soluble organic compounds with high molecular weight. Daily samples were collected during two campaigns in winter and summer, for 1 month each. The concentration of high-molecular-weight compounds (humic-like substances (HULIS)) was between 0.4 and 4 μg/m3 in winter and summer. The most intense signals in the LDI-MS mass spectra were measured between m/z150 and 500, comparing well with the mode of the two main high mass peaks determined with SEC corresponding to masses between 200 and 600 Da. For the maximum molecular weight, however, different results were obtained by the two techniques: whereas a maximum molecular weight between 1300 and 3300 Da was found with SEC, hardly any peaks above m/z700 were measured with LDI-MS. During summer the maximum molecular weight of HULIS (determined with SEC) correlates positively with several parameters such as ozone and increased temperature indicative of enhanced atmospheric photo-oxidation. The HULIS concentration also correlates positively with the oxalic acid concentration in the particles. This suggests that HULIS are generated by secondary processes in summer. The lack of such correlations during winter suggests that other sources and processes might be important during colder seasons.

Thyroid disruption in zebrafish (Danio rerio) larvae: Different molecular response patterns lead to impaired eye development and visual functions

The vertebrate thyroid system is important for multiple developmental processes, including eye development. Thus, its environmentally induced disruption may impact important fitness-related parameters like visual capacities and behaviour. The present study investigated the relation between molecular effects of thyroid disruption and morphological and physiological changes of eye development in zebrafish (Danio rerio). Two test compounds representing different molecular modes of thyroid disruption were used: propylthiouracil (PTU), which is an enzyme-inhibitor of thyroid hormone synthesis, and tetrabromobisphenol A (TBBPA), which interacts with the thyroid hormone receptors. Both chemicals significantly altered transcript levels of thyroid system-related genes (TRα, TRβ, TPO, TSH, DIO1, DIO2 and DIO3) in a compound-specific way. Despite these different molecular response patterns, both treatments resulted in similar pathological alterations of the eyes such as reduced size, RPE cell diameter and pigmentation, which were concentration-dependent. The morphological changes translated into impaired visual performance of the larvae: the optokinetic response was significantly and concentration-dependently decreased in both treatments, together with a significant increase of light preference of PTU-treated larvae. In addition, swimming activity was impacted. This study provides first evidence that different modes of molecular action of the thyroid disruptors can be associated with uniform apical responses. Furthermore, this study is the first to show that pathological eye development, as it can be induced by exposure to thyroid disruptors, indeed translates into impaired visual capacities of zebrafish early life stages.

Parameterizing a Simple Potential Energy Function Appropriate for Vibrational Frequencies Using Molecular Dynamics

A Maple scheme for quickly parameterizing vibrational potential energy functions is presented. As an example, the potential energy function's parameters for the vibrational motions in H_2O_2 are obtained assuming the simplest potential energy function. This paper was originally written as a research paper, but rejected by the referees. It is therefore being edited into an ``educational'' paper for student usage.

Response of the calcifying coccolithophore Emiliania huxleyi to low pH/high pCO2: from physiology to molecular level, 2011

The emergence of ocean acidification as a significant threat to calcifying organisms in marine ecosystems creates a pressing need to understand the physiological and molecular mechanisms by which calcification is affected by environmental parameters. We report here, for the first time, changes in gene expression induced by variations in pH/pCO2 in the widespread and abundant coccolithophore Emiliania huxleyi. Batch cultures were subjected to increased partial pressure of CO2 (pCO2; i.e. decreased pH), and the changes in expression of four functional gene classes directly or indirectly related to calcification were investigated. Increased pCO2 did not affect the calcification rate and only carbonic anhydrase transcripts exhibited a significant down-regulation. Our observation that elevated pCO2 induces only limited changes in the transcription of several transporters of calcium and bicarbonate gives new significant elements to understand cellular mechanisms underlying the early response of E. huxleyi to CO2-driven ocean acidification.