Feedback on hypothalamic TRH transcription is dependent on thyroid hormone receptor N terminus.

The beta thyroid hormone receptor (TRbeta), but not TRalpha1, plays a specific role in mediating T(3)-dependent repression of hypothalamic TRH transcription. To investigate the structural basis of isoform specificity, we compared the transcriptional regulation and DNA binding obtained with chimeric and N-terminally deleted TRs. Using in vivo transfection assays to follow hypothalamic TRH transcription in the mouse brain, we found that TRbeta1 and chimeras with the TRbeta1 N terminus did not affect either transcriptional activation or repression from the rat TRH promoter, whereas N-terminally deleted TRbeta1 impaired T(3)-dependent repression. TRalpha1 or chimeras with the TRalpha1 N terminus reduced T(3)-independent transcriptional activation and blocked T(3)-dependent repression of transcription. Full deletion of the TRalpha1 N terminus restored ligand-independent activation of transcription. No TR isoform specificity was seen after transcription from a positive thyroid hormone response element. Gel mobility assays showed that all TRs tested bound specifically to the main negative thyroid hormone response element in the TRH promoter (site 4). Addition of neither steroid receptor coactivator 1 nor nuclear extracts from the hypothalamic paraventricular nuclei revealed any TR isoform specificity in binding to site 4. Thus N-terminal sequences specify TR T(3)-dependent repression of TRH transcription but not DNA recognition, emphasizing as yet unknown neuron-specific contributions to protein-promoter interactions in vivo.

Revisiting the specificity of the MHC class II transactivator CIITA in vivo.

CIITA is a master regulatory factor for the expression of MHC class II (MHC-II) and accessory genes involved in Ag presentation. It has recently been suggested that CIITA also regulates numerous other genes having diverse functions within and outside the immune system. To determine whether these genes are indeed relevant targets of CIITA in vivo, we studied their expression in CIITA-transgenic and CIITA-deficient mice. In contrast to the decisive control of MHC-II and related genes by CIITA, nine putative non-MHC target genes (Eif3s2, Kpna6, Tap1, Yars, Col1a2, Ctse, Ptprr, Tnfsf6 and Plxna1) were found to be CIITA independent in all cell types examined. Two other target genes, encoding IL-4 and IFN-gamma, were indeed found to be up- and down-regulated, respectively, in CIITA-transgenic CD4(+) T cells. However, there was no correlation between MHC-II expression and this Th2 bias at the level of individual transgenic T cells, indicating an indirect control by CIITA. These results show that MHC-II-restricted Ag presentation, and its indirect influences on T cells, remains the only pathway under direct control by CIITA in vivo. They also imply that precisely regulated MHC-II expression is essential for maintaining a proper Th1-Th2 balance.

New inhibitors of Helicobacter pylori urease holoenzyme selected from phage-displayed peptide libraries.

Urease is an important virulence factor for Helicobacter pylori and is critical for bacterial colonization of the human gastric mucosa. Specific inhibition of urease activity has been proposed as a possible strategy to fight this bacteria which infects billions of individual throughout the world and can lead to severe pathological conditions in a limited number of cases. We have selected peptides which specifically bind and inhibit H. pylori urease from libraries of random peptides displayed on filamentous phage in the context of pIII coat protein. Screening of a highly diverse 25-mer combinatorial library and two newly constructed random 6-mer peptide libraries on solid phase H. pylori urease holoenzyme allowed the identification of two peptides, 24-mer TFLPQPRCSALLRYLSEDGVIVPS and 6-mer YDFYWW that can bind and inhibit the activity of urease purified from H. pylori. These two peptides were chemically synthesized and their inhibition constants (Ki) were found to be 47 microM for the 24-mer and 30 microM for the 6-mer peptide. Both peptides specifically inhibited the activity of H. pylori urease but not that of Bacillus pasteurii.

Role of penicillin-binding proteins in Streptococcus gordonii

Résumé Streptococcus gordonii est une bactérie colonisatrice naturelle de la cavité buccale de l'homme. Bien que normalement commensale, elle peut causer des infections graves, telles que des bactériémies ou des endocardites infectieuses. La pénicilline étant un des traitements privilégiés dans de tels cas, l'augmentation rapide et globale des résistances à cet antibiotique devient inquiétante. L'étude de la physiologie et des bases génétiques de ces résistances chez S. gordonii s'avère donc importante. Les cibles moléculaires privilégiées de la pénicilline G et des β-lactames sont les penicilllin-binding proteins (PBPs). Ces enzymes associées à la membrane ont pour rôle de catalyser les réactions de transpeptidation et de transglycosylation, qui constituent les dernières étapes de la biosynthèse du peptidoglycan (PG). Elles sont définies comme classe A ou B selon leur capacité d'assurer soit les deux réactions, soit uniquement la transpeptidation. Les β-lactames inhibent le domaine transpeptidase de toutes les PBPs, entraînant l'inhibition de la synthèse du PG, l'inhibition de la croissance, et finalement la mort cellulaire. Chez les streptocoques, les PBPs sont aussi les premiers déterminants de la résistance à la pénicilline. De plus, elles sont impliquées dans la morphologie bactérienne, en raison de leur rôle crucial dans la formation du PG. Le but de ce travail était de caractériser les PBPs de S. gordonii et d'étudier leurs fonctions dans la vie végétative de la bactérie ainsi que durant le développement de la résistance à la pénicilline. Premièrement, des mutants auxquels il manque une ou deux PBP(s) ont été construits. Leur étude - au niveau physiologique, biochimique et morphologique - a montré le caractère essentiel ou dispensable de chaque protéine, ainsi que certaines de leurs fonctions potentielles. Deuxièmement, des mutants résistants à la pénicilline ont été générés. Leur caractérisation a montré l'importance des mutations dans les PBPs ainsi que dans d'autres gènes encore inconnus, de même que le rôle crucial des PBPs de classe A dans le développement de la résistance à la pénicilline. Des expériences supplémentaires sur des isolats résistants ont aussi prouvé que la résistance a un coût en terme de fitness, coût que S. gordonii parvient à compenser par des mécanismes d'adaptation. Finalement, les promoteurs des gènes des PBPs ont été déterminés et leur expression a été étudiée grâce au gène de luciférase. Il a ainsi été montré que la résistance à la pénicilline entraîne non seulement des altérations au niveau des protéines, mais aussi au niveau de la régulation des gènes. De plus, la pénicilline génère directement des modifications dans l'expression de PBPs spécifiques. Summary Streptococcus gordonii is a normal inhabitant of the human oral cavity and a pioneer colonizer of teeth. Although usually considered as a commensal, this organism can cause life-threatening infections such as bacteraemia or endocarditis. Since penicillin is one of the preferential treatments for such pathologies, the rapid and general increase of antibiotic resistance in the overall population becomes an issue. Thus, studying the physiologic and genetic bases of such a resistance in S. gordonii is of interest. The primary molecular targets of penicillin G and other β-lactams are the so called penicillin-binding proteins (PBPs). These are membrane-associated proteins that catalyze the last steps in peptidoglycan (PG) biosynthesis, namely transpeptidation and transglycosylation. Depending on their capacity to catalyze either reactions or only transpeptidation, they are considered as class A or class B PBPs, respectively. β-lactam antibiotics inhibit the transpeptidase domain of both of these classes of enzymes, resulting in inhibition of PG assembly, inhibition of bacterial growth, and ultimately leading to cell death. In streptococci, PBPs are also the primary determinants of penicillin-resistance. Moreover, because of their crucial role in PG formation, they are implicated in fundamental aspects of cell morphology. The goal of this work was thus to characterize S. gordonii PBPs and to explore their functions in terms of vegetative life and penicillin-resistance development. First, single and double PBP-inactivated mutants were generated and their effect on the bacterial physiology, cell wall biochemistry and ultrastructural morphology was assessed. This demonstrated the essentiality or dispensability of each protein for bacterial life. Second, penicillin-resistant mutants were generated by cyclic exposure to increasing concentrations of the drug. Characterization of these mutants pointed out the importance of both PBP and non-PBP mutations, as well as the crucial role of the class A PBPs in the development of penicillin-resistance. Further experiments on resistant isolates demonstrated the fitness cost of this resistance, but also the capacity of S. gordonii to adapt and regain the fitness of the wild-type. Finally, the promoters of PBP genes were determined and their expression was monitored using luciferase fusions. This showed that penicillin-resistance, in addition to modifications at the level of the protein, also triggered genetic alterations. Moreover, penicillin itself generated modifications in the expression of specific PBPs.

Design and analysis of ChIP-Seq experiments for nuclear factor I DNA-binding proteins

SUMMARY : Eukaryotic DNA interacts with the nuclear proteins using non-covalent ionic interactions. Proteins can recognize specific nucleotide sequences based on the sterical interactions with the DNA and these specific protein-DNA interactions are the basis for many nuclear processes, e.g. gene transcription, chromosomal replication, and recombination. New technology termed ChIP-Seq has been recently developed for the analysis of protein-DNA interactions on a whole genome scale and it is based on immunoprecipitation of chromatin and high-throughput DNA sequencing procedure. ChIP-Seq is a novel technique with a great potential to replace older techniques for mapping of protein-DNA interactions. In this thesis, we bring some new insights into the ChIP-Seq data analysis. First, we point out to some common and so far unknown artifacts of the method. Sequence tag distribution in the genome does not follow uniform distribution and we have found extreme hot-spots of tag accumulation over specific loci in the human and mouse genomes. These artifactual sequence tags accumulations will create false peaks in every ChIP-Seq dataset and we propose different filtering methods to reduce the number of false positives. Next, we propose random sampling as a powerful analytical tool in the ChIP-Seq data analysis that could be used to infer biological knowledge from the massive ChIP-Seq datasets. We created unbiased random sampling algorithm and we used this methodology to reveal some of the important biological properties of Nuclear Factor I DNA binding proteins. Finally, by analyzing the ChIP-Seq data in detail, we revealed that Nuclear Factor I transcription factors mainly act as activators of transcription, and that they are associated with specific chromatin modifications that are markers of open chromatin. We speculate that NFI factors only interact with the DNA wrapped around the nucleosome. We also found multiple loci that indicate possible chromatin barrier activity of NFI proteins, which could suggest the use of NFI binding sequences as chromatin insulators in biotechnology applications. RESUME : L'ADN des eucaryotes interagit avec les protéines nucléaires par des interactions noncovalentes ioniques. Les protéines peuvent reconnaître les séquences nucléotidiques spécifiques basées sur l'interaction stérique avec l'ADN, et des interactions spécifiques contrôlent de nombreux processus nucléaire, p.ex. transcription du gène, la réplication chromosomique, et la recombinaison. Une nouvelle technologie appelée ChIP-Seq a été récemment développée pour l'analyse des interactions protéine-ADN à l'échelle du génome entier et cette approche est basée sur l'immuno-précipitation de la chromatine et sur la procédure de séquençage de l'ADN à haut débit. La nouvelle approche ChIP-Seq a donc un fort potentiel pour remplacer les anciennes techniques de cartographie des interactions protéine-ADN. Dans cette thèse, nous apportons de nouvelles perspectives dans l'analyse des données ChIP-Seq. Tout d'abord, nous avons identifié des artefacts très communs associés à cette méthode qui étaient jusqu'à présent insoupçonnés. La distribution des séquences dans le génome ne suit pas une distribution uniforme et nous avons constaté des positions extrêmes d'accumulation de séquence à des régions spécifiques, des génomes humains et de la souris. Ces accumulations des séquences artéfactuelles créera de faux pics dans toutes les données ChIP-Seq, et nous proposons différentes méthodes de filtrage pour réduire le nombre de faux positifs. Ensuite, nous proposons un nouvel échantillonnage aléatoire comme un outil puissant d'analyse des données ChIP-Seq, ce qui pourraient augmenter l'acquisition de connaissances biologiques à partir des données ChIP-Seq. Nous avons créé un algorithme d'échantillonnage aléatoire et nous avons utilisé cette méthode pour révéler certaines des propriétés biologiques importantes de protéines liant à l'ADN nommés Facteur Nucléaire I (NFI). Enfin, en analysant en détail les données de ChIP-Seq pour la famille de facteurs de transcription nommés Facteur Nucléaire I, nous avons révélé que ces protéines agissent principalement comme des activateurs de transcription, et qu'elles sont associées à des modifications de la chromatine spécifiques qui sont des marqueurs de la chromatine ouverte. Nous pensons que lés facteurs NFI interagir uniquement avec l'ADN enroulé autour du nucléosome. Nous avons également constaté plusieurs régions génomiques qui indiquent une éventuelle activité de barrière chromatinienne des protéines NFI, ce qui pourrait suggérer l'utilisation de séquences de liaison NFI comme séquences isolatrices dans des applications de la biotechnologie.

Binding under Conflict Conditions: State-Space Analysis of Multivariate EEG Synchronization.

Real-world objects are often endowed with features that violate Gestalt principles. In our experiment, we examined the neural correlates of binding under conflict conditions in terms of the binding-by-synchronization hypothesis. We presented an ambiguous stimulus ("diamond illusion") to 12 observers. The display consisted of four oblique gratings drifting within circular apertures. Its interpretation fluctuates between bound ("diamond") and unbound (component gratings) percepts. To model a situation in which Gestalt-driven analysis contradicts the perceptually explicit bound interpretation, we modified the original diamond (OD) stimulus by speeding up one grating. Using OD and modified diamond (MD) stimuli, we managed to dissociate the neural correlates of Gestalt-related (OD vs. MD) and perception-related (bound vs. unbound) factors. Their interaction was expected to reveal the neural networks synchronized specifically in the conflict situation. The synchronization topography of EEG was analyzed with the multivariate S-estimator technique. We found that good Gestalt (OD vs. MD) was associated with a higher posterior synchronization in the beta-gamma band. The effect of perception manifested itself as reciprocal modulations over the posterior and anterior regions (theta/beta-gamma bands). Specifically, higher posterior and lower anterior synchronization supported the bound percept, and the opposite was true for the unbound percept. The interaction showed that binding under challenging perceptual conditions is sustained by enhanced parietal synchronization. We argue that this distributed pattern of synchronization relates to the processes of multistage integration ranging from early grouping operations in the visual areas to maintaining representations in the frontal networks of sensory memory.

Mini-review: Specificity and expression of CIITA, the master regulator of MHC class II genes.

The class II transactivator (CIITA) has been referred to as the "master control factor" for the expression of MHC class II (MHCII) genes. As our knowledge on the specificity and function of CIITA grows, it is becoming increasingly evident that this sobriquet is entirely justified. First, despite extensive investigations, the major target genes of CIITA remain those implicated in the presentation of antigenic peptides by MHCII molecules. Although other putative target genes have been reported, the contribution of CIITA to their expression remains indirect, controversial or comparatively minor relative to its decisive role as a regulator of MHCII and related genes. Second, the most important parameter dictating MHCII expression is by far the expression pattern of the gene encoding CIITA (MHC2TA). The vast majority of signals that activate or repress MHCII expression under physiological and pathological situations converge on one or more of the three alternative promoters that drive transcription of the MHC2TA gene. In short, with respect to its specificity and its exquisitely controlled pattern of expression, CIITA is by a long stretch the single most important transcription factor for the regulation of genes required for MHCII-restricted antigen-presentation.

Fatty acids and eicosanoids regulate gene expression through direct interactions with peroxisome proliferator-activated receptors alpha and gamma.

Peroxisome proliferator-activated receptors (PPARs) alpha and gamma are key regulators of lipid homeostasis and are activated by a structurally diverse group of compounds including fatty acids, eicosanoids, and hypolipidemic drugs such as fibrates and thiazolidinediones. While thiazolidinediones and 15-deoxy-Delta12, 14-prostaglandin J2 have been shown to bind to PPARgamma, it has remained unclear whether other activators mediate their effects through direct interactions with the PPARs or via indirect mechanisms. Here, we describe a novel fibrate, designated GW2331, that is a high-affinity ligand for both PPARalpha and PPARgamma. Using GW2331 as a radioligand in competition binding assays, we show that certain mono- and polyunsaturated fatty acids bind directly to PPARalpha and PPARgamma at physiological concentrations, and that the eicosanoids 8(S)-hydroxyeicosatetraenoic acid and 15-deoxy-Delta12,14-prostaglandin J2 can function as subtype-selective ligands for PPARalpha and PPARgamma, respectively. These data provide evidence that PPARs serve as physiological sensors of lipid levels and suggest a molecular mechanism whereby dietary fatty acids can modulate lipid homeostasis.

Cell-specific localization of monocarboxylate transporters, MCT1 and MCT2, in the adult mouse brain revealed by double immunohistochemical labeling and confocal microscopy.

Recent evidence suggests that lactate could be a preferential energy substrate transferred from astrocytes to neurons. This would imply the presence of specific transporters for lactate on both cell types. We have investigated the immunohistochemical localization of two monocarboxylate transporters, MCT1 and MCT2, in the adult mouse brain. Using specific antibodies raised against MCT1 and MCT2, we found strong immunoreactivity for each transporter in glia limitans, ependymocytes and several microvessel-like elements. In addition, small processes distributed throughout the cerebral parenchyma were immunolabeled for monocarboxylate transporters. Double immunofluorescent labeling and confocal microscopy examination of these small processes revealed no co-localization between glial fibrillary acidic protein and monocarboxylate transporters, although many glial fibrillary acidic protein-positive processes were often in close apposition to elements labeled for monocarboxylate transporters. In contrast, several elements expressing the S100beta protein, another astrocytic marker found to be located in distinct parts of the same cell when compared with glial fibrillary acidic protein, were also strongly immunoreactive for MCT1, suggesting expression of this transporter by astrocytes. In contrast, MCT2 was expressed in a small subset of microtubule-associated protein-2-positive elements, indicating a neuronal localization. In conclusion, these observations are consistent with the possibility that lactate, produced and released by astrocytes (via MCT1), could be taken up (via MCT2) and used by neurons as an energy substrate.

Angiotensin-converting enzyme inhibition by hydroxamic zinc-binding idrapril in humans.

The new angiotensin-converting enzyme (ACE) inhibitor idrapril acts by binding the catalytically important zinc ion to a hydroxamic group. We investigated its pharmacodynamic and pharmacokinetic properties in 8 healthy men: Increasing doses of 1, 5, and 25 mg idrapril as well as placebo or 5 mg captopril were administered intravenously (i.v.) at 1-week intervals. Six of the subjects received 100 mg idrapril orally (p.o.) last, and two ingested oral placebo as a double-blind control. Blood pressure (BP) and heart rate (HR) remained unchanged. No serious side effects were observed. ACE inhibition in vivo was evaluated by changes in the ratio of specifically measured plasma angiotensin II (AngII) and AngI concentrations determined by high-performance liquid chromatography/radioimmunoassay (HPLC/RIA) techniques. Plasma ACE activity in vitro was estimated by radioenzymatic assay; it was suppressed by > or = 93% at 15 min after injection of 25 mg idrapril or 5 mg captopril and by 96% 2 h after idrapril intake. Mean AngII levels were decreased dose dependently at 15 min after idrapril injections. At the same time, plasma renin activity (PRA) and AngI increased according to the doses. The AngII/AngI ratio was clearly related to plasma idrapril levels (r = -0.88, n = 60). Oral idrapril inhibited ACE maximally at 1-4 h after dosing, when < 7% of initial ACE activity was observed in vitro and in vivo. Idrapril is a safe and efficient ACE inhibitor in human subjects. It is well absorbed orally. Besides having a slightly slower onset of action, idrapril has pharmacodynamic effects comparable to those of captopril.

Expression of human estrogen receptor mutants in Xenopus oocytes: correlation between transcriptional activity and ability to form protein-DNA complexes.

The human estrogen receptor (hER) is a trans-acting regulatory protein composed of a series of discrete functional domains. We have microinjected an hER expression vector (HEO) into Xenopus oocyte nuclei and demonstrate, using Western blot assay, that the hER is synthesized. When nuclear extracts from oocytes were prepared and incubated in the presence of a 2.7 kb DNA fragment comprising the 5' end of the vitellogenin gene B2, formation of estrogen-dependent complexes could be visualized by electron microscopy over the estrogen responsive element (ERE). Of crucial importance is the observation that the complex formation is inhibited by the estrogen antagonist tamoxifen, is restored by the addition of the hormone and does not take place with extracts from control oocytes injected with the expression vector lacking the sequences encoding the receptor. The presence of the biologically active hER is confirmed in co-injection experiments, in which HEO is co-introduced with a CAT reporter gene under the control of a vitellogenin promoter containing or lacking the ERE. CAT assays and primer extensions analyses reveal that both the receptor and the ERE are essential for estrogen induced stimulation of transcription. The same approach was used to analyze selective hER mutants. We find that the DNA binding domain (region C) is essential for protein--DNA complex formation at the ERE but is not sufficient by itself to activate transcription from the reporter gene. In addition to region C, both the hormone binding (region E) and amino terminal (region A/B) domains are needed for an efficient transcription activation.(ABSTRACT TRUNCATED AT 250 WORDS)

Ligand binding transmits conformational changes across the membrane-spanning region to the intracellular side of the 5-HT3 serotonin receptor.

Conformational changes of channel activation: Five enhanced green fluorescent protein (EGFP) molecules (green cylinders) were integrated into the intracellular part of the homopentameric ionotropic 5-HT3 receptor. This allowed the detection of extracellular binding of fluorescent ligands (?) to EGFP by FRET, and also enabled the quantification of agonist-induced conformational changes in the intracellular region of the receptor by homo-FRET between EGFPs. The approach opens novel ways for probing receptor activation and functional screening of therapeutic compounds.

Nucleic acid-binding properties of the RRM-containing protein RDM1.

RDM1 (RAD52 Motif 1) is a vertebrate protein involved in the cellular response to the anti-cancer drug cisplatin. In addition to an RNA recognition motif, RDM1 contains a small amino acid motif, named RD motif, which it shares with the recombination and repair protein, RAD52. RDM1 binds to single- and double-stranded DNA, and recognizes DNA distortions induced by cisplatin adducts in vitro. Here, we have performed an in-depth analysis of the nucleic acid-binding properties of RDM1 using gel-shift assays and electron microscopy. We show that RDM1 possesses acidic pH-dependent DNA-binding activity and that it binds RNA as well as DNA, and we present evidence from competition gel-shift experiments that RDM1 may be capable of discrimination between the two nucleic acids. Based on reported studies of RAD52, we have generated an RDM1 variant mutated in its RD motif. We find that the L119GF --> AAA mutation affects the mode of RDM1 binding to single-stranded DNA.

Differential involvement of peroxisome-proliferator-activated receptors alpha and delta in fibrate and fatty-acid-mediated inductions of the gene encoding liver fatty-acid-binding protein in the liver and the small intestine.

Liver fatty-acid-binding protein (L-FABP) is a cytoplasmic polypeptide that binds with strong affinity especially to long-chain fatty acids (LCFAs). It is highly expressed in both the liver and small intestine, where it is thought to have an essential role in the control of the cellular fatty acid (FA) flux. Because expression of the gene encoding L-FABP is increased by both fibrate hypolipidaemic drugs and LCFAs, it seems to be under the control of transcription factors, termed peroxisome-proliferator-activated receptors (PPARs), activated by fibrate or FAs. However, the precise molecular mechanism by which these regulations take place remain to be fully substantiated. Using transfection assays, we found that the different PPAR subtypes (alpha, gamma and delta) are able to mediate the up-regulation by FAs of the gene encoding L-FABP in vitro. Through analysis of LCFA- and fibrate-mediated effects on L-FABP mRNA levels in wild-type and PPARalpha-null mice, we have found that PPARalpha in the intestine does not constitute a dominant regulator of L-FABP gene expression, in contrast with what is known in the liver. Only the PPARdelta/alpha agonist GW2433 is able to up-regulate the gene encoding L-FABP in the intestine of PPARalpha-null mice. These findings demonstrate that PPARdelta can act as a fibrate/FA-activated receptor in tissues in which it is highly expressed and that L-FABP is a PPARdelta target gene in the small intestine. We propose that PPARdelta contributes to metabolic adaptation of the small intestine to changes in the lipid content of the diet.

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.