Ab initio modeling and molecular dynamics simulation of the alpha 1b-adrenergic receptor activation.

This work describes the ab initio procedure employed to build an activation model for the alpha 1b-adrenergic receptor (alpha 1b-AR). The first version of the model was progressively modified and complicated by means of a many-step iterative procedure characterized by the employment of experimental validations of the model in each upgrading step. A combined simulated (molecular dynamics) and experimental mutagenesis approach was used to determine the structural and dynamic features characterizing the inactive and active states of alpha 1b-AR. The latest version of the model has been successfully challenged with respect to its ability to interpret and predict the functional properties of a large number of mutants. The iterative approach employed to describe alpha 1b-AR activation in terms of molecular structure and dynamics allows further complications of the model to allow prediction and interpretation of an ever-increasing number of experimental data.

Spatial and temporal dynamics of jasmonate synthesis and accumulation in Arabidopsis in response to wounding.

A new metabolite profiling approach combined with an ultrarapid sample preparation procedure was used to study the temporal and spatial dynamics of the wound-induced accumulation of jasmonic acid (JA) and its oxygenated derivatives in Arabidopsis thaliana. In addition to well known jasmonates, including hydroxyjasmonates (HOJAs), jasmonoyl-isoleucine (JA-Ile), and its 12-hydroxy derivative (12-HOJA-Ile), a new wound-induced dicarboxyjasmonate, 12-carboxyjasmonoyl-l-isoleucine (12-HOOCJA-Ile) was discovered. HOJAs and 12-HOOCJA-Ile were enriched in the midveins of wounded leaves, strongly differentiating them from the other jasmonate metabolites studied. The polarity of these oxylipins at physiological pH correlated with their appearance in midveins. When the time points of accumulation of different jasmonates were determined, JA levels were found to increase within 2-5 min of wounding. Remarkably, these changes occurred throughout the plant and were not restricted to wounded leaves. The speed of the stimulus leading to JA accumulation in leaves distal to a wound is at least 3 cm/min. The data give new insights into the spatial and temporal accumulation of jasmonates and have implications in the understanding of long-distance wound signaling in plants.

Intracellular excision and reintegration dynamics of the ICEclc genomic island of Pseudomonas knackmussii sp. strain B13.

Genomic islands are DNA elements acquired by horizontal gene transfer that are common to a large number of bacterial genomes, which can contribute specific adaptive functions, e.g. virulence, metabolic capacities or antibiotic resistances. Some genomic islands are still self-transferable and display an intricate life-style, reminiscent of both bacteriophages and conjugative plasmids. Here we studied the dynamical process of genomic island excision and intracellular reintegration using the integrative and conjugative element ICEclc from Pseudomonas knackmussii B13 as model. By using self-transfer of ICEclc from strain B13 to Pseudomonas putida and Cupriavidus necator as recipients, we show that ICEclc can target a number of different tRNA(Gly) genes in a bacterial genome, but only those which carry the GCC anticodon. Two conditional traps were designed for ICEclc based on the attR sequence, and we could show that ICEclc will insert with different frequencies in such traps producing brightly fluorescent cells. Starting from clonal primary transconjugants we demonstrate that ICEclc is excising and reintegrating at detectable frequencies, even in the absence of recipient. Recombination site analysis provided evidence to explain the characteristics of a larger number of genomic island insertions observed in a variety of strains, including Bordetella petri, Pseudomonas aeruginosa and Burkholderia.

The effects of physiologically plausible connectivity structure on local and global dynamics in large scale brain models.

Functionally relevant large scale brain dynamics operates within the framework imposed by anatomical connectivity and time delays due to finite transmission speeds. To gain insight on the reliability and comparability of large scale brain network simulations, we investigate the effects of variations in the anatomical connectivity. Two different sets of detailed global connectivity structures are explored, the first extracted from the CoCoMac database and rescaled to the spatial extent of the human brain, the second derived from white-matter tractography applied to diffusion spectrum imaging (DSI) for a human subject. We use the combination of graph theoretical measures of the connection matrices and numerical simulations to explicate the importance of both connectivity strength and delays in shaping dynamic behaviour. Our results demonstrate that the brain dynamics derived from the CoCoMac database are more complex and biologically more realistic than the one based on the DSI database. We propose that the reason for this difference is the absence of directed weights in the DSI connectivity matrix.

Temporal and spatial dynamics of benthic invertebrate communities in an alpine glacier-fed alluvial system

Summary The present thesis work focused on the ecology of benthic invertebrates in the proglacial floodplain of the Rhone in the Swiss Alps. The main glacial Rhone River and a smaller glacial tributary, the Mutt River, joined and entered a braiding multi-thread area. A first part concentrated on the disruption of the longitudinal patterns of environmental conditions and benthic invertebrate fauna in the Rhone by its tributary the Mutt. The Mutt had less harsh environmental conditions, higher taxonomic richness and more abundant zoobenthos compared to the Rhone upstream of the confluence. Although the habitat conditions in the main stream were little modified by the tributary, the fauna was richer and more diverse below the confluence. Colonisation from the Mutt induced the occurrence of faunal elements uncommon of glacial streams in the upper Rhone, where water temperature remains below 4°C. Although the glacial Rhone dominated the system with regard to hydrology and certain environmental conditions, the Mutt tributary has to be seen as the faunal driver of the system. The second part of the study concerned the spatio-temporal differentiation of the habitats and the benthic communities along and across the flood plain. No longitudinal differentiation was found. The spatial transversal differentiation of three habitat types with different environmental characteristics was successfully reflected in the spatial variability of benthic assemblages. This typology separated marginal sites of the flood plain, left bank sites under the influence of the Mutt, and the right bank sites under the influence of the Rh6ne. Faunistic spatial differences were emphasized by the quantitative structure of the fauna, richness, abundances and Simpson index of diversity. Seasonal environmental variability was positively related with Simpson index of diversity and the total richness per site. Low flow conditions were the most favourable season for the fauna and November was characterized by low spatial environmental heterogeneity, high spatial heterogeneity of faunal assemblage, maximum taxonomic richness, a particular taxonomic composition, highest abundances, as well as the highest primary food resources. The third part studied the egg development of three species of Ephemeroptera in the laboratory at 1.5 to 7°C and the ecological implications in the field. Species revealed very contrasting development strategies. Baetis alpinus has a synchronous and efficient egg development, which is faster in warmer habitats, enabling it to exploit short periods of favourable conditions in the floodplain. Ecdyonurus picteti has a very long development time slightly decreasing in warmer conditions. The high degree of individual variation suggests a genetic determination of the degree-days demand. Combined with the glacial local conditions, this strategy leads to an extreme delay of hatching and allows it to develop in very unpredictable habitats. Rhithrogena nivata is the second cold adapted species in Ephemeroptera. The incubation duration is long and success largely depends on the timing of hatching and the discharge conditions. This species is able to exploit extremely unstable and cold habitats where other species are limited by low water temperatures. The fourth part dealt with larval development in different habitats of the floodplain. Addition of data on egg development allowed the description of the life histories of the species from oviposition until emergence. Rhithrogena nivata and loyolaea generally have a two-year development, with the first winter passed as eggs and the second one as larvae. Development of Ecdyonurus picteti is difficult to document but appears to be efficient in a harsh and unpredictable environment. Baetis alpinus was studied separately in four habitats of the floodplain system with contrasting thermal regimes. Differences in success and duration of larval development and in growth rates are emphasised. Subvention mechanisms between habitats by migration of young or grown larvae were demonstrated. Development success and persistence of the populations in the system were thus increased. Emergence was synchronised to the detriment of the optimisation of the adult's size and fecundity. These very different development strategies induce a spatial and temporal distribution in the use of food resources and ecological niches. The last part of this work aimed at the synthesis of the characteristics and the ecological features of three distinct compartments of the system that are the upper Rhone, the Mutt and the floodplain. Their particular role as well as their inter-dependence concerning the structure and the dynamics of the benthic communities was emphasised. Résumé Ce travail de thèse est consacré à l'écologie des invertébrés benthiques dans la zone alluviale proglaciaire du Rhône dans les Alpes suisses. Le Rhône, torrent glaciaire principal, reçoit les eaux de la Mutt, affluent glaciaire secondaire, puis pénètre dans une zone de tressage formée de plusieurs bras. La première partie de l'étude se concentre sur la disruption par la Mutt des processus longitudinaux, tant environnementaux que faunistiques, existants dans le Rhône. Les conditions environnementales régnant dans la Mutt sont moins rudes, la richesse taxonomique plus élevée et le zoobenthos plus abondant que dans le Rhône en amont de la confluence. Bien que les conditions environnementales dans le torrent principal soient peu modifiées par l'affluent, la faune s'avère être plus riche et plus diversifiée en aval de la confluence. La colonisation depuis la Mutt permet l'occurrence de taxons inhabituels dans le Rhône en amont de la confluence, où la température de l'eau se maintient en dessous de 4°C. Bien que le Rhône, torrent glaciaire principal, domine le système du point de vu de l'hydrologie et de certains paramètres environnementaux, l'affluent Mutt doit être considéré comme l'élément structurant la faune dans le système. La deuxième partie concerne la différentiation spatiale et temporelle des habitats et des communautés benthiques à travers la plaine alluviale. Aucune différentiation longitudinale n'a été mise en évidence. La différentiation transversale de trois types d'habitats sur la base des caractéristiques environnementales a été confirmée par la variabilité spatiale de la faune. Cette typologie sépare les sites marginaux de la plaine alluviale, ceux sous l'influence de la Mutt (en rive gauche) et ceux sous l'influence du Rhône amont (en rive droite). Les différences spatiales de la faune sont mises en évidence par la structure quantitative de la faune, la richesse, les abondances et l'indice de diversité de Simpson. La variabilité saisonnière du milieu est positivement liée avec l'indice de diversité de Simpson et la richesse totale par site. L'étiage correspond à la période la plus favorable pour la faune et novembre réunit des conditions de faible hétérogénéité spatiale du milieu, de forte hétérogénéité spatiale de la faune, une richesse taxonomique maximale, une composition faunistique particulière, les abondances ainsi que les ressources primaires les plus élevées. La troisième partie est consacrée à l'étude du développement des oeufs de trois espèces d'Ephémères au laboratoire à des températures de 1.5 à 7°C, ainsi qu'aux implications écologiques sur le terrain. Ces espèces présentent des stratégies de développement très contrastées. Baetis alpinus a un développement synchrone et efficace, plus rapide en milieu plus chaud et lui permettant d'exploiter les courtes périodes de conditions favorables. Ecdyonurus picteti présente une durée de développement très longue, diminuant légèrement dans des conditions plus chaudes. L'importante variation interindividuelle suggère un déterminisme génétique de la durée de développement. Cette stratégie, associée aux conditions locales, conduit à un décalage extrême des éclosions et permet à l'espèce de se développer dans des habitats imprévisibles. Rhithrogena nivata est la seconde espèce d'Ephémères présentant une adaptation au froid. L'incubation des oeufs est longue et son succès dépend de la période des éclosions et des conditions hydrologiques. Cette espèce est capable d'exploiter des habitats extrêmement instables et froids, où la température est facteur limitant pour d'autres espèces. La quatrième partie traite du développement larvaire dans différents habitats de la plaine alluviale. Le développement complet est décrit pour les espèces étudiées de la ponte jusqu'à l'émergence. Rhithrogena nivata et loyolaea atteignent généralement le stade adulte en deux ans, le premier hiver étant passé sous forme d'oeuf et le second sous forme de larve. Le développement de Ecdyonurus picteti est difficile à documenter, mais s'avère cependant efficace dans un environnement rude et imprévisible. Baetis alpinus a été étudié séparément dans quatre habitats de la plaine ayant des régimes thermiques contrastés. La réussite et la durée du développement embryonnaire ainsi que les taux de croissance y sont variables. Des mécanismes de subvention entre habitats sont possibles par la migration de larves juvéniles ou plus développées, augmentant ainsi la réussite du développement et le maintien des populations dans le système. L'émergence devient synchrone, au détriment de l'optimisation de la taille et de la fécondité des adultes. Ces stratégies très différentes induisent une distribution spatiale et temporelle dans l'usage des ressources et des niches écologiques. La dernière partie synthétise les caractéristiques écologiques des trois compartiments du système que sont le Rhône amont, la Mutt et la zone alluviale. Leurs rôles particuliers et leurs interdépendances du point de vue de la structure et de la dynamique des communautés benthiques sont mis en avant.

Manual reaction times and brain dynamics after 'awake surgery' of slow-growing tumors invading the parietal area. A case report

Primary objectives: Awake surgeries of slow-growing tumours invading the brain and guided by direct electrical stimulation induce major brain reorganizations accompanied with slight impairments post-operatively. In most cases, these deficits are so slight after a few days that they are often not detectable on classical neuropsychological evaluations. Consequently, this study investigated whether simple visuo-manual reaction time paradigms would sign some level of functional asymmetries between both hemispheres. Importantly, the visual stimulus was located in the saggital plane in order to limit attentional biases and to focus mainly on the inter-hemispheric asymmetry. Methods and procedures: Three patients (aged 41, 59 and 59 years) after resections in parietal regions and a control group (age¼44, SD¼6.9) were compared during simple uni- and bimanual reaction times (RTs). Main outcomes and results: Longer RTs were observed for the contralesional compared to the ipsilesional hand in the unimanual condition. This asymmetry was reversed for the bimanual condition despite longer RTs. Conclusion and clinical implications: Reaction time paradigms are useful in these patients to monitor more precisely their functional deficits, especially their level of functional asymmetry, and to understand brain (re)organization following slowgrowing lesions.

The Caribbean-North America-Cocos Triple Junction and the dynamics of the Polochic-Motagua fault systems: Pull-up and zipper models

The Polochic-Motagua fault systems (PMFS) are part of the sinistral transform boundary between the North American and Caribbean plates. To the west, these systems interact with the subduction zone of the Cocos plate, forming a subduction-subduction-transform triple junction. The North American plate moves westward relative to the Caribbean plate. This movement does not affect the geometry of the subducted Cocos plate, which implies that deformation is accommodated entirely in the two overriding plates. Structural data, fault kinematic analysis, and geomorphic observations provide new elements that help to understand the late Cenozoic evolution of this triple junction. In the Miocene, extension and shortening occurred south and north of the Motagua fault, respectively. This strain regime migrated northward to the Polochic fault after the late Miocene. This shift is interpreted as a ``pull-up'' of North American blocks into the Caribbean realm. To the west, the PMFS interact with a trench-parallel fault zone that links the Tonala fault to the Jalpatagua fault. These faults bound a fore-arc sliver that is shared by the two overriding plates. We propose that the dextral Jalpatagua fault merges with the sinistral PMFS, leaving behind a suturing structure, the Tonala fault. This tectonic ``zipper'' allows the migration of the triple junction. As a result, the fore-arc sliver comes into contact with the North American plate and helps to maintain a linear subduction zone along the trailing edge of the Caribbean plate. All these processes currently make the triple junction increasingly diffuse as it propagates eastward and inland within both overriding plates.

Three essays on short-term macroeconomics : business fluctuations, large devaluations and inflation dynamics

Résumé: At least since the Great Depression, explaining why there are business fluctuations has been one of the biggest challenges that the science of economics has had to face. The hope is that if we could better understand recessions, then we could also be more successful in overcoming them. This dissertation consists of three papers that are part of the general endeavor of economists to understand these fluctuations. The first paper discusses, for a particular model, whether a result related to fluctuations would still hold if time were modeled as continuous rather than discrete. The two other papers focus on price stickiness. The second paper discusses why, after a large devaluation, prices of non-tradables may change by only a small amount in comparison to the magnitude of the devaluation. The third paper examines price adjustment in a model in which information is imperfect and it is costly to change prices.

A thematic analysis of delusion with religious content in schizophrenia: open, closed, and mixed dynamics

The aim of the present study was to elicit how patients with delusions with religious contents conceptualized or experienced their spirituality and religiousness. Sixty-two patients with present or past religious delusions went through semistructured interviews, which were analyzed using the three coding steps described in the grounded theory. Three major themes were found in religious delusions: ''spiritual identity,'' ''meaning of illness,'' and ''spiritual figures.'' One higher-order concept was found: ''structure of beliefs.'' We identified dynamics that put these personal beliefs into a constant reconstruction through interaction with the world and others (i.e., open dynamics) and conversely structural dynamics that created a complete rupture with the surrounding world and others (i.e., closed structural dynamics); those dynamics may coexist. These analyses may help to identify psychological functions of delusions with religious content and, therefore, to better conceptualize interventions when dealing with it in psychotherapy.

Effects of anxiety on visuospatial pattern discrimination and decision-making processes: what can we learn from brain dynamics ?

Detection and discrimination of visuospatial input involve at least extracting, selecting and encoding relevant information and decision-making processes allowing selecting a response. These two operations are altered, respectively, by attentional mechanisms that change discrimination capacities, and by beliefs concerning the likelihood of uncertain events. Information processing is tuned by the attentional level that acts like a filter on perception, while decision-making processes are weighed by subjective probability of risk. In addition, it has been shown that anxiety could affect the detection of unexpected events through the modification of the level of arousal. Consequently, purpose of this study concerns whether and how decision-making and brain dynamics are affected by anxiety. To investigate these questions, the performance of women with either a high (12) or a low (12) STAI-T (State-Trait Anxiety Inventory, Spielberger, 1983) was examined in a decision-making visuospatial task where subjects have to recognize a target visual pattern from non-target patterns. The target pattern was a schematic image of furniture arranged in such a way as to give the impression of a living room. Non-target patterns were created by either the compression or the dilatation of the distances between objects. Target and non-target patterns were always presented in the same configuration. Preliminary behavioral results show no group difference in reaction time. In addition, visuo-spatial abilities were analyzed trough the signal detection theory for quantifying perceptual decisions in the presence of uncertainty (Green and Swets, 1966). This theory treats detection of a stimulus as a decision-making process determined by the nature of the stimulus and cognitive factors. Astonishingly, no difference in d' (corresponding to the distance between means of the distributions) and c (corresponds to the likelihood ratio) indexes was observed. Comparison of Event-related potentials (ERP) reveals that brain dynamics differ according to anxiety. It shows differences in component latencies, particularly a delay in anxious subjects over posterior electrode sites. However, these differences are compensated during later components by shorter latencies in anxious subjects compared to non-anxious one. These inverted effects seem indicate that the absence of difference in reaction time rely on a compensation of attentional level that tunes cortical activation in anxious subjects, but they have to hammer away to maintain performance.

Inference of chromosomal inversion dynamics from Pool-Seq data in natural and laboratory populations of Drosophila melanogaster.

Sequencing of pools of individuals (Pool-Seq) represents a reliable and cost-effective approach for estimating genome-wide SNP and transposable element insertion frequencies. However, Pool-Seq does not provide direct information on haplotypes so that, for example, obtaining inversion frequencies has not been possible until now. Here, we have developed a new set of diagnostic marker SNPs for seven cosmopolitan inversions in Drosophila melanogaster that can be used to infer inversion frequencies from Pool-Seq data. We applied our novel marker set to Pool-Seq data from an experimental evolution study and from North American and Australian latitudinal clines. In the experimental evolution data, we find evidence that positive selection has driven the frequencies of In(3R)C and In(3R)Mo to increase over time. In the clinal data, we confirm the existence of frequency clines for In(2L)t, In(3L)P and In(3R)Payne in both North America and Australia and detect a previously unknown latitudinal cline for In(3R)Mo in North America. The inversion markers developed here provide a versatile and robust tool for characterizing inversion frequencies and their dynamics in Pool-Seq data from diverse D. melanogaster populations.

Stable hydrogen and oxygen isotope composition of waters from mine tailings in different climatic environments

The stable isotope composition of waters (delta H-2, delta O-18) can be used as a natural tracer of hydrologic processes in systems affected by acid mine drainage. We investigated the delta H-2 and delta O-18 values of pore waters from four oxidizing sulfidic mine tailings impoundments in different climatic regions of Chile (Piuquenes at La Andina with Alpine climate, Cauquenes and Caren at El Teniente with Mediterranean climate, and Talabre at the Chuquicamata deposit with hyperarid climate). No clear relationship was found between altitude and isotopic composition. The observed displacement of the tailings pore waters from the local meteoric water line toward higher delta O-18 values (by similar to +2% delta O-18 relative to delta H-2) is partly due to water-rock interaction processes, including hydration and O-isotope exchange with sulfates and Fe(III) oxyhydroxides produced by pyrite oxidation. In most tailings, from the saturated zone toward the surface, isotopically different zones can be distinguished. Zone I is characterized by an upward depletion of H-2 and O-18 in the pore waters from the saturated zone and the lowermost vadose zone, due to ascending diffused isotopically light water triggered by the constant loss of water vapor by evaporation at the surface. In zone II, the capillary flow of a mix of vapor and liquid water causes an evaporative isotopic enrichment in H-2 and O-18. At the top of the tailings in dry climate a zone III between the capillary zone and the surface contains isotopically light diffused and atmospheric water vapor. In temperate climates, the upper part of the profile is affected by recent rainfall and zone III may not differ isotopically from zone II.

Relations of the job demands-control-support model of job strain with personality attributes : a cross-cultural study in Switzerland and South Africa

Among the various work stress models, one of the most popular to date is the job demands-­‐control (JDC) model developed by Karasek (1979), which postulates that work-­‐related strain will be the highest under work conditions characterized by high demands and low autonomy. The absence of social support at work will further increase negative outcomes. However, this model does not apply equally to all individuals and to all cultures. In the following studies, we assessed work characteristics, personality traits, culture-­‐driven individual attributes, and work-­‐related health outcomes, through the administration of questionnaires. The samples consist of Swiss (n = 622) and South African (n = 879) service-­‐oriented employees (from health, finance, education and commerce sectors) and aged from 18 to 65 years old. Results generally confirm the universal contribution of high psychological demands, low decision latitude and low supervisor support at work, as well as high neuroticism predict the worse health outcomes among employees in both countries. Furthermore, low neuroticism plays a moderating role between psychological demands and burnout, while high openness and high conscientiousness each play a moderating role between decision latitude and burnout in South Africa. Results also reveal that culture-­‐driven individual attributes play a role in both countries, but in a unique manner and according to the ethnic group of belonging. Given that organizations are increasingly characterized with multicultural employees as well as increasingly adverse and complex job conditions, our results help in identifying more updated and refined dynamics that are key between the employee and the work environment in today's context. -- L'un des modèles sur le stress au travail des plus répandus est celui développé par Karasek (1979), qui postule qu'une mauvaise santé chez les employés résulte d'une combinaison de demandes psychologiques élevées, d'une latitude décisionnelle faible et de l'absence de soutien social au travail. Néanmoins, ce modèle ne s'applique pas de façon équivalente chez tous les individus et dans toutes les cultures. Dans les études présentées, nous avons mesuré les caractéristiques de travail, les traits de personnalité, les traits culturels et les effets lies à la santé à l'aide de questionnaires. L'échantillon provient de la Suisse (n = 622) et de l'Afrique du Sud (n = 879) et comprend des employés de domaines divers en lien avec le service (notamment des secteurs de la santé, finance, éducation et commerce) tous âgés entre 18 et 65 ans. Les résultats confirment l'universalité des effets directs des demandes au travail, la latitude décisionnelle faible, le soutien social faible provenant du supérieur hiérarchique, ainsi que le névrosisme élevé qui contribuent à un niveau de santé faible au travail, et ce, dans les deux pays. De plus, un niveau faible de névrosisme a un effet de modération entre les demandes au travail et l'épuisement professionnel, alors que l'ouverture élevée et le caractère consciencieux élevé modèrent la relation entre la latitude décisionnelle et l'épuisement professionnel en Afrique du Sud. Nous avons aussi trouvé que les traits culturels jouent un rôle dans les deux pays, mais de façon unique et en fonction du groupe ethnique d'appartenance. Sachant que les organisations sont de plus en plus caractérisées par des employés d'origine ethnique variées, et que les conditions de travail se complexifient, nos résultats contribuent à mieux comprendre les dynamiques entre l'employé et l'environnement de travail contemporain. personnalité, différences individuelles, comparaisons culturelles, culture, stress au travail, épuisement professionnel, santé des employés.

In silico prediction of human carboxylesterase-1 (hCES1) metabolism combining docking analyses and MD simulations.

Metabolic problems lead to numerous failures during clinical trials, and much effort is now devoted in developing in silico models predicting metabolic stability and metabolites. Such models are well known for cytochromes P450 and some transferases, whereas little has been done to predict the hydrolytic activity of human hydrolases. The present study was undertaken to develop a computational approach able to predict the hydrolysis of novel esters by human carboxylesterase hCES1. The study involves both docking analyses of known substrates to develop predictive models, and molecular dynamics (MD) simulations to reveal the in situ behavior of substrates and products, with particular attention being paid to the influence of their ionization state. The results emphasize some crucial properties of the hCES1 catalytic cavity, confirming that as a trend with several exceptions, hCES1 prefers substrates with relatively smaller and somewhat polar alkyl/aryl groups and larger hydrophobic acyl moieties. The docking results underline the usefulness of the hydrophobic interaction score proposed here, which allows a robust prediction of hCES1 catalysis, while the MD simulations show the different behavior of substrates and products in the enzyme cavity, suggesting in particular that basic substrates interact with the enzyme in their unprotonated form.