Computational toolbox towards evolutionary domain mapping of membrane proteins

Membrane proteins account for about 20% to 30% of all proteins encoded in a typical genome. They play central roles in multiple cellular processes mediating the interaction of the cell with its surrounding. Over 60% of all drug targets contain a membrane domain. The experimental difficulties of obtaining a crystal structural severely limits our ability or understanding of membrane protein function. Computational evolutionary studies of proteins are crucial for the prediction of 3D structures. In this project, we construct a tool able to quantify the evolutionary positive selective pressure on each residue of membrane proteins through maximum likelihood phylogeny reconstruction. The conservation plot combined with a structural homology model is also a potent tool to predict those residues that have essentials roles in the structure and function of a membrane protein and can be very useful in the design of validation experiments.

MgO/Ag(001) interface structure and STM images from first principles

The interface of MgO/Ag(001) has been studied with density functional theory applied to slabs. We have found that regular MgO films show a small adhesion to the silver substrate, the binding can be increased in off-stoichiometric regimes, either by the presence of O vacancies at the oxide film or by a small excess of O atoms at the interface between the ceramic to the metal. By means of theoretical methods, the scanning tunneling microscopy signatures of these films is also analyzed in some detail. For defect free deposits containing 1 or 2 ML and at low voltages, tunnelling takes place from the surface Ag substrate, and at large positive voltages Mg atoms are imaged. If defects, oxygen vacancies, are present on the surface of the oxide they introduce much easier channels for tunnelling resulting in big protrusions and controlling the shape of the image, the extra O stored at the interface can also be detected for very thin films.

Rockfall susceptibility assessment and remote geological mapping with LiDAR point clouds

Characterizing the geological features and structures in three dimensions over inaccessible rock cliffs is needed to assess natural hazards such as rockfalls and rockslides and also to perform investigations aimed at mapping geological contacts and building stratigraphy and fold models. Indeed, the detailed 3D data, such as LiDAR point clouds, allow to study accurately the hazard processes and the structure of geologic features, in particular in vertical and overhanging rock slopes. Thus, 3D geological models have a great potential of being applied to a wide range of geological investigations both in research and applied geology projects, such as mines, tunnels and reservoirs. Recent development of ground-based remote sensing techniques (LiDAR, photogrammetry and multispectral / hyperspectral images) are revolutionizing the acquisition of morphological and geological information. As a consequence, there is a great potential for improving the modeling of geological bodies as well as failure mechanisms and stability conditions by integrating detailed remote data. During the past ten years several large rockfall events occurred along important transportation corridors where millions of people travel every year (Switzerland: Gotthard motorway and railway; Canada: Sea to sky highway between Vancouver and Whistler). These events show that there is still a lack of knowledge concerning the detection of potential rockfalls, making mountain residential settlements and roads highly risky. It is necessary to understand the main factors that destabilize rocky outcrops even if inventories are lacking and if no clear morphological evidences of rockfall activity are observed. In order to increase the possibilities of forecasting potential future landslides, it is crucial to understand the evolution of rock slope stability. Defining the areas theoretically most prone to rockfalls can be particularly useful to simulate trajectory profiles and to generate hazard maps, which are the basis for land use planning in mountainous regions. The most important questions to address in order to assess rockfall hazard are: Where are the most probable sources for future rockfalls located? What are the frequencies of occurrence of these rockfalls? I characterized the fracturing patterns in the field and with LiDAR point clouds. Afterwards, I developed a model to compute the failure mechanisms on terrestrial point clouds in order to assess the susceptibility to rockfalls at the cliff scale. Similar procedures were already available to evaluate the susceptibility to rockfalls based on aerial digital elevation models. This new model gives the possibility to detect the most susceptible rockfall sources with unprecented detail in the vertical and overhanging areas. The results of the computation of the most probable rockfall source areas in granitic cliffs of Yosemite Valley and Mont-Blanc massif were then compared to the inventoried rockfall events to validate the calculation methods. Yosemite Valley was chosen as a test area because it has a particularly strong rockfall activity (about one rockfall every week) which leads to a high rockfall hazard. The west face of the Dru was also chosen for the relevant rockfall activity and especially because it was affected by some of the largest rockfalls that occurred in the Alps during the last 10 years. Moreover, both areas were suitable because of their huge vertical and overhanging cliffs that are difficult to study with classical methods. Limit equilibrium models have been applied to several case studies to evaluate the effects of different parameters on the stability of rockslope areas. The impact of the degradation of rockbridges on the stability of large compartments in the west face of the Dru was assessed using finite element modeling. In particular I conducted a back-analysis of the large rockfall event of 2005 (265'000 m3) by integrating field observations of joint conditions, characteristics of fracturing pattern and results of geomechanical tests on the intact rock. These analyses improved our understanding of the factors that influence the stability of rock compartments and were used to define the most probable future rockfall volumes at the Dru. Terrestrial laser scanning point clouds were also successfully employed to perform geological mapping in 3D, using the intensity of the backscattered signal. Another technique to obtain vertical geological maps is combining triangulated TLS mesh with 2D geological maps. At El Capitan (Yosemite Valley) we built a georeferenced vertical map of the main plutonio rocks that was used to investigate the reasons for preferential rockwall retreat rate. Additional efforts to characterize the erosion rate were made at Monte Generoso (Ticino, southern Switzerland) where I attempted to improve the estimation of long term erosion by taking into account also the volumes of the unstable rock compartments. Eventually, the following points summarize the main out puts of my research: The new model to compute the failure mechanisms and the rockfall susceptibility with 3D point clouds allows to define accurately the most probable rockfall source areas at the cliff scale. The analysis of the rockbridges at the Dru shows the potential of integrating detailed measurements of the fractures in geomechanical models of rockmass stability. The correction of the LiDAR intensity signal gives the possibility to classify a point cloud according to the rock type and then use this information to model complex geologic structures. The integration of these results, on rockmass fracturing and composition, with existing methods can improve rockfall hazard assessments and enhance the interpretation of the evolution of steep rockslopes. -- La caractérisation de la géologie en 3D pour des parois rocheuses inaccessibles est une étape nécessaire pour évaluer les dangers naturels tels que chutes de blocs et glissements rocheux, mais aussi pour réaliser des modèles stratigraphiques ou de structures plissées. Les modèles géologiques 3D ont un grand potentiel pour être appliqués dans une vaste gamme de travaux géologiques dans le domaine de la recherche, mais aussi dans des projets appliqués comme les mines, les tunnels ou les réservoirs. Les développements récents des outils de télédétection terrestre (LiDAR, photogrammétrie et imagerie multispectrale / hyperspectrale) sont en train de révolutionner l'acquisition d'informations géomorphologiques et géologiques. Par conséquence, il y a un grand potentiel d'amélioration pour la modélisation d'objets géologiques, ainsi que des mécanismes de rupture et des conditions de stabilité, en intégrant des données détaillées acquises à distance. Pour augmenter les possibilités de prévoir les éboulements futurs, il est fondamental de comprendre l'évolution actuelle de la stabilité des parois rocheuses. Définir les zones qui sont théoriquement plus propices aux chutes de blocs peut être très utile pour simuler les trajectoires de propagation des blocs et pour réaliser des cartes de danger, qui constituent la base de l'aménagement du territoire dans les régions de montagne. Les questions plus importantes à résoudre pour estimer le danger de chutes de blocs sont : Où se situent les sources plus probables pour les chutes de blocs et éboulement futurs ? Avec quelle fréquence vont se produire ces événements ? Donc, j'ai caractérisé les réseaux de fractures sur le terrain et avec des nuages de points LiDAR. Ensuite, j'ai développé un modèle pour calculer les mécanismes de rupture directement sur les nuages de points pour pouvoir évaluer la susceptibilité au déclenchement de chutes de blocs à l'échelle de la paroi. Les zones sources de chutes de blocs les plus probables dans les parois granitiques de la vallée de Yosemite et du massif du Mont-Blanc ont été calculées et ensuite comparés aux inventaires des événements pour vérifier les méthodes. Des modèles d'équilibre limite ont été appliqués à plusieurs cas d'études pour évaluer les effets de différents paramètres sur la stabilité des parois. L'impact de la dégradation des ponts rocheux sur la stabilité de grands compartiments de roche dans la paroi ouest du Petit Dru a été évalué en utilisant la modélisation par éléments finis. En particulier j'ai analysé le grand éboulement de 2005 (265'000 m3), qui a emporté l'entier du pilier sud-ouest. Dans le modèle j'ai intégré des observations des conditions des joints, les caractéristiques du réseau de fractures et les résultats de tests géoméchaniques sur la roche intacte. Ces analyses ont amélioré l'estimation des paramètres qui influencent la stabilité des compartiments rocheux et ont servi pour définir des volumes probables pour des éboulements futurs. Les nuages de points obtenus avec le scanner laser terrestre ont été utilisés avec succès aussi pour produire des cartes géologiques en 3D, en utilisant l'intensité du signal réfléchi. Une autre technique pour obtenir des cartes géologiques des zones verticales consiste à combiner un maillage LiDAR avec une carte géologique en 2D. A El Capitan (Yosemite Valley) nous avons pu géoréferencer une carte verticale des principales roches plutoniques que j'ai utilisé ensuite pour étudier les raisons d'une érosion préférentielle de certaines zones de la paroi. D'autres efforts pour quantifier le taux d'érosion ont été effectués au Monte Generoso (Ticino, Suisse) où j'ai essayé d'améliorer l'estimation de l'érosion au long terme en prenant en compte les volumes des compartiments rocheux instables. L'intégration de ces résultats, sur la fracturation et la composition de l'amas rocheux, avec les méthodes existantes permet d'améliorer la prise en compte de l'aléa chute de pierres et éboulements et augmente les possibilités d'interprétation de l'évolution des parois rocheuses.

What can ecosystems learn? Expanding evolutionary ecology with learning theory.

BACKGROUND: The structure and organisation of ecological interactions within an ecosystem is modified by the evolution and coevolution of the individual species it contains. Understanding how historical conditions have shaped this architecture is vital for understanding system responses to change at scales from the microbial upwards. However, in the absence of a group selection process, the collective behaviours and ecosystem functions exhibited by the whole community cannot be organised or adapted in a Darwinian sense. A long-standing open question thus persists: Are there alternative organising principles that enable us to understand and predict how the coevolution of the component species creates and maintains complex collective behaviours exhibited by the ecosystem as a whole? RESULTS: Here we answer this question by incorporating principles from connectionist learning, a previously unrelated discipline already using well-developed theories on how emergent behaviours arise in simple networks. Specifically, we show conditions where natural selection on ecological interactions is functionally equivalent to a simple type of connectionist learning, 'unsupervised learning', well-known in neural-network models of cognitive systems to produce many non-trivial collective behaviours. Accordingly, we find that a community can self-organise in a well-defined and non-trivial sense without selection at the community level; its organisation can be conditioned by past experience in the same sense as connectionist learning models habituate to stimuli. This conditioning drives the community to form a distributed ecological memory of multiple past states, causing the community to: a) converge to these states from any random initial composition; b) accurately restore historical compositions from small fragments; c) recover a state composition following disturbance; and d) to correctly classify ambiguous initial compositions according to their similarity to learned compositions. We examine how the formation of alternative stable states alters the community's response to changing environmental forcing, and we identify conditions under which the ecosystem exhibits hysteresis with potential for catastrophic regime shifts. CONCLUSIONS: This work highlights the potential of connectionist theory to expand our understanding of evo-eco dynamics and collective ecological behaviours. Within this framework we find that, despite not being a Darwinian unit, ecological communities can behave like connectionist learning systems, creating internal conditions that habituate to past environmental conditions and actively recalling those conditions. REVIEWERS: This article was reviewed by Prof. Ricard V Solé, Universitat Pompeu Fabra, Barcelona and Prof. Rob Knight, University of Colorado, Boulder.

A Social Network Analysis of Twitter: Mapping the Digital Humanities Community

Defining digital humanities might be an endless debate if we stick to the discussion about the boundaries of this concept as an academic "discipline". In an attempt to concretely identify this field and its actors, this paper shows that it is possible to analyse them through Twitter, a social media widely used by this "community of practice". Based on a network analysis of 2,500 users identified as members of this movement, the visualisation of the "who's following who?" graph allows us to highlight the structure of the network's relationships, and identify users whose position is particular. Specifically, we show that linguistic groups are key factors to explain clustering within a network whose characteristics look similar to a small world.

Statistical Inversion Theory in X-ray Tomography

In this thesis the X-ray tomography is discussed from the Bayesian statistical viewpoint. The unknown parameters are assumed random variables and as opposite to traditional methods the solution is obtained as a large sample of the distribution of all possible solutions. As an introduction to tomography an inversion formula for Radon transform is presented on a plane. The vastly used filtered backprojection algorithm is derived. The traditional regularization methods are presented sufficiently to ground the Bayesian approach. The measurements are foton counts at the detector pixels. Thus the assumption of a Poisson distributed measurement error is justified. Often the error is assumed Gaussian, altough the electronic noise caused by the measurement device can change the error structure. The assumption of Gaussian measurement error is discussed. In the thesis the use of different prior distributions in X-ray tomography is discussed. Especially in severely ill-posed problems the use of a suitable prior is the main part of the whole solution process. In the empirical part the presented prior distributions are tested using simulated measurements. The effect of different prior distributions produce are shown in the empirical part of the thesis. The use of prior is shown obligatory in case of severely ill-posed problem.

Small but powerful: Top predator local extinction affects ecosystem structure and function in an intermittent stream.

Top predator loss is a major global problem, with a current trend in biodiversity loss towards high trophic levels that modifies most ecosystems worldwide. Most research in this area is focused on large-bodied predators, despite the high extinction risk of small-bodied freshwater fish that often act as apex consumers. Consequently, it remains unknown if intermittent streams are affected by the consequences of top-predators' extirpations. The aim of our research was to determine how this global problem affects intermittent streams and, in particular, if the loss of a small-bodied top predator (1) leads to a 'mesopredator release', affects primary consumers and changes whole community structures, and (2) triggers a cascade effect modifying the ecosystem function. To address these questions, we studied the topdown effects of a small endangered fish species, Barbus meridionalis (the Mediterranean barbel), conducting an enclosure/exclosure mesocosm experiment in an intermittent stream where B. meridionalis became locally extinct following a wildfire.We found that top predator absence led to 'mesopredator release', and also to 'prey release' despite intraguild predation, which contrasts with traditional food web theory. In addition, B. meridionalis extirpation changed whole macroinvertebrate community composition and increased total macroinvertebrate density. Regarding ecosystem function, periphyton primary production decreased in apex consumer absence. In this study, the apex consumer was functionally irreplaceable; its local extinction led to the loss of an important functional role that resulted in major changes to the ecosystem's structure and function. This study evidences that intermittent streams can be affected by the consequences of apex consumers' extinctions, and that the loss of small-bodied top predators can lead to large ecosystem changes. We recommend the reintroduction of small-bodied apex consumers to systems where they have been extirpated, to restore ecosystem structure and function.

Parallel characterizations of a generalized shapley value and a generalized banzhaf value for cooperative games with level structure of cooperation

We present parallel characterizations of two different values in the framework of restricted cooperation games. The restrictions are introduced as a finite sequence of partitions defined on the player set, each of them being coarser than the previous one, hence forming a structure of different levels of a priori unions. On the one hand, we consider a value first introduced in Ref. [18], which extends the Shapley value to games with different levels of a priori unions. On the other hand, we introduce another solution for the same type of games, which extends the Banzhaf value in the same manner. We characterize these two values using logically comparable properties.

Self-consistent theory of current and voltage noise in multimode ballistic conductors

Electron transport in a self-consistent potential along a ballistic two-terminal conductor has been investigated. We have derived general formulas which describe the nonlinear current-voltage characteristics, differential conductance, and low-frequency current and voltage noise assuming an arbitrary distribution function and correlation properties of injected electrons. The analytical results have been obtained for a wide range of biases: from equilibrium to high values beyond the linear-response regime. The particular case of a three-dimensional Fermi-Dirac injection has been analyzed. We show that the Coulomb correlations are manifested in the negative excess voltage noise, i.e., the voltage fluctuations under high-field transport conditions can be less than in equilibrium.

Capital Structure of Finnish SMEs and Financial Constraints

This thesis studies capital structure of Finnish small and medium sized enterprises. The specific object of the study is to test whether financial constraints have an effect on capital structure. In addition influences of several other factors were studied. Capital structure determinants are formulated based on three capital structure theories. The tradeoff theory and the agency theory concentrate on the search of optimal capital structure. The pecking order theory concerns favouring on financing source over another. The data of this study consists of financial statement data and results of corporate questionnaire. Regression analysis was used to find out the effects of several determinants. Regression models were formed based on the presented theories. Short and long term debt ratios were considered separately. The metrics of financially constrained firms was included in all models. It was found that financial constrains have a negative and significant effect to short term debt ratios. The effect was negative also to long term debt ratio but not statistically significant. Other considerable factors that influenced debt ratios were fixed assets, age, profitability, single owner and sufficiency of internal financing.

Hybrid LC Filter for Power Electronic Drives: Theory and Implementation

Power electronic converter drives use, for the sake of high efficiency, pulse-width modulation that results in sequences of high-voltage high-frequency steep-edged pulses. Such a signal contains a set of high harmonics not required for control purposes. Harmonics cause reflections in the cable between the motor and the inverter leading to faster winding insulation ageing. Bearing failures and problems with electromagnetic compatibility may also result. Electrical du/dt filters provide an effective solution to problems caused by pulse-width modulation, thereby increasing the performance and service life of the electrical machines. It is shown that RLC filters effectively decrease the reflection phenomena in the cable. Improved (simple, but effective) solutions are found for both differential- and common-mode signals; these solutions use a galvanic connection between the RLC filter star point and the converter DC link. Foil chokes and film capacitors are among the most widely used components in high-power applications. In actual applications they can be placed in different parts of the cabinet. This fact complicates the arrangement of the cabinet and decreases the reliability of the system. In addition, the inductances of connection wires may prevent filtration at high frequencies. This thesis introduces a new hybrid LC filter that uses a natural capacitance between the turns of the foil choke based on integration of an auxiliary layer into it. The main idea of the hybrid LC filter results from the fact that both the foil choke and the film capacitors have the same roll structure. Moreover, the capacitance between the turns (“intra capacitance”) of the foil inductors is the reason for the deterioration of their properties at high frequencies. It is shown that the proposed filter has a natural cancellation of the intra capacitance. A hybrid LC filter may contain two or more foil layers isolated from each other and coiled on a core. The core material can be iron or even air as in the filter considered in this work. One of the foils, called the main foil, can be placed between the inverter and the motor cable. Other ones, called auxiliary foils, may be connected in star to create differential-mode noise paths, and then coupled to the DC link midpoint to guarantee a traveling path, especially for the common-mode currents. This way, there is a remarkable capacitance between the main foil and the auxiliary foil. Investigations showed that such a system can be described by a simple equivalent LC filter in a wide range of frequencies. Because of its simple hybrid construction, the proposed LC filter can be a cost-effective and competitive solution for modern power drives. In the thesis, the application field of the proposed filter is considered and determined. The basics of hybrid LC filter design are developed further. High-frequency behaviour of the proposed filter is analysed by simulations. Finally, the thesis presents experimental data proving that the hybrid LC filter can be used for du/dt of PWM pulses and reduction of common-mode currents.

Modulating the electronic structure of amino acids: interaction of model lewis acids with anthranilic acid

On the basis of theoretical B3LYP calculations, Yáñez and co-workers (J. Chem. Theory Comput. 2012, 8, 2293) illustrated that beryllium ions are capable of significantly modulating (changing) the electronic structures of imidazole. In this computational organic chemistry study, the interaction of this β-amino acid and five model Lewis acids (BeF1+, Be2+, AlF2(1+), AlF2+, and Al3+) were investigated. Several aspects were addressed: natural bond orbitals, including second order perturbation analysis of intra-molecular charge delocalization and the natural population analysis atomic charges; molecular geometries; selected infrared stretching frequencies (C-N, C-O, and N-H), and selected ¹H-NMR chemical shifts. The data illustrate that this interaction can weaken the H-O bond and goes beyond strengthening the intra-molecular hydrogen bond (N...H-O) to cause a spontaneous transfer of the proton to the nitrogen atom in five cases generating zwitterion structures. Many new features are observed. Most importantly, the zwitterion structures include a stabilizing hydrogen bond (N-H...O) that varies in relative strength according to the Lewis acid. These findings explain the experimental observations of α-amino acids (for example: J. Am. Chem. Soc. 2001, 123, 3577) and are the first reported fundamental electronic structure characterization of β-amino acids in zwitterion form.

Analysis of the structure and vibrational spectra of glucose and fructose

Molecular modelling using semiempirical methods AM1, PM3, PM5 and, MINDO as well as the Density Functional Theory method BLYP/DZVP respectively were used to calculate the structure and vibrational spectra of d-glucose and d-fructose in their open chain, alpha-anomer and beta-anomer monohydrate forms. The calculated data show that both molecules are not linear; ground state and the number for the point-group C is equal to 1. Generally, the results indicate that there are similarities in bond lengths and vibrational modes of both molecules. It is concluded that DFT could be used to study both the structural and vibrational spectra of glucose and fructose.

Evaluation of group electronegativities and hardness (softness) of group 14 elements and containing functional groups through density functional theory and correlation with NMR spectra data

Quantum Chemical calculations for group 14 elements of Periodic Table (C, Si, Ge, Sn, Pb) and their functional groups have been carried out using Density Functional Theory (DFT) based reactivity descriptors such as group electronegativities, hardness and softness. DFT calculations were performed for a large series of tetracoordinated Sn compounds of the CH3SnRR'X type, where X is a halogen and R and R' are alkyl, halogenated alkyl, alkoxy, or alkyl thio groups. The results were interpreted in terms of calculated electronegativity and hardness of the SnRR'X groups, applying a methodology previously developed by Geerlings and coworkers (J. Phys. Chem. 1993, 97, 1826). These calculations allowed to see the regularities concerning the influence of the nature of organic groups RR' and inorganic group X on electronegativities and hardness of the SnRR'X groups; in this case, it was found a very good correlation between the electronegativity of the fragment and experimental 119Sn chemical shifts, a property that sensitively reflects the change in the valence electronic structure of molecules. This work was complemented with the study of some compounds of the EX and ER types, where E= C, Si, Ge, Sn and R= CH3, H, which was performed to study the influence that the central atom has on the electronegativity and hardness of molecules, or whether these properties are mainly affected for the type of ligand bound to the central atom. All these calculations were performed using the B3PW91 functional together with the 6-311++G** basis set level for H, C, Si, Ge, F, Cl and Br atoms and the 3-21G for Sn and I atoms.