Exploring unfrequent events with accelerated molecular dynamics simulations: from photochemistry to drug design

La meva incorporació al grup de recerca del Prof. McCammon (University of California San Diego) en qualitat d’investigador post doctoral amb una beca Beatriu de Pinós, va tenir lloc el passat 1 de desembre de 2010; on vaig dur a terme les meves tasques de recerca fins al darrer 1 d’abril de 2012. El Prof. McCammon és un referent mundial en l’aplicació de simulacions de dinàmica molecular (MD) en sistemes biològics d’interès humà. La contribució més important del Prof. McCammon en la simulació de sistemes biològics és el desenvolupament del mètode de dinàmiques moleculars accelerades (AMD). Les simulacions MD convencionals, les quals estan limitades a l’escala de temps del nanosegon (~10-9s), no son adients per l’estudi de sistemes biològics rellevants a escales de temps mes llargues (μs, ms...). AMD permet explorar fenòmens moleculars poc freqüents però que son clau per l’enteniment de molts sistemes biològics; fenòmens que no podrien ser observats d’un altre manera. Durant la meva estada a la “University of California San Diego”, vaig treballar en diferent aplicacions de les simulacions AMD, incloent fotoquímica i disseny de fàrmacs per ordinador. Concretament, primer vaig desenvolupar amb èxit una combinació dels mètodes AMD i simulacions Car-Parrinello per millorar l’exploració de camins de desactivació (interseccions còniques) en reaccions químiques fotoactivades. En segon lloc, vaig aplicar tècniques estadístiques (Replica Exchange) amb AMD en la descripció d’interaccions proteïna-lligand. Finalment, vaig dur a terme un estudi de disseny de fàrmacs per ordinador en la proteïna-G Rho (involucrada en el desenvolupament de càncer humà) combinant anàlisis estructurals i simulacions AMD. Els projectes en els quals he participat han estat publicats (o estan encara en procés de revisió) en diferents revistes científiques, i han estat presentats en diferents congressos internacionals. La memòria inclosa a continuació conté més detalls de cada projecte esmentat.

Orthotropic active strain models for the numerical simulation of cardiac biomechanics

An active strain formulation for orthotropic constitutive laws arising in cardiac mechanics modeling is introduced and studied. The passive mechanical properties of the tissue are described by the Holzapfel-Ogden relation. In the active strain formulation, the Euler-Lagrange equations for minimizing the total energy are written in terms of active and passive deformation factors, where the active part is assumed to depend, at the cell level, on the electrodynamics and on the specific orientation of the cardiac cells. The well-posedness of the linear system derived from a generic Newton iteration of the original problem is analyzed and different mechanical activation functions are considered. In addition, the active strain formulation is compared with the classical active stress formulation from both numerical and modeling perspectives. Taylor-Hood and MINI finite elements are employed to discretize the mechanical problem. The results of several numerical experiments show that the proposed formulation is mathematically consistent and is able to represent the main key features of the phenomenon, while allowing savings in computational costs.

Simulation methods with extended stability for stiff biochemical kinetics

Background: With increasing computer power, simulating the dynamics of complex systems in chemistry and biology is becoming increasingly routine. The modelling of individual reactions in (bio)chemical systems involves a large number of random events that can be simulated by the stochastic simulation algorithm (SSA). The key quantity is the step size, or waiting time, τ, whose value inversely depends on the size of the propensities of the different channel reactions and which needs to be re-evaluated after every firing event. Such a discrete event simulation may be extremely expensive, in particular for stiff systems where τ can be very short due to the fast kinetics of some of the channel reactions. Several alternative methods have been put forward to increase the integration step size. The so-called τ-leap approach takes a larger step size by allowing all the reactions to fire, from a Poisson or Binomial distribution, within that step. Although the expected value for the different species in the reactive system is maintained with respect to more precise methods, the variance at steady state can suffer from large errors as τ grows. Results: In this paper we extend Poisson τ-leap methods to a general class of Runge-Kutta (RK) τ-leap methods. We show that with the proper selection of the coefficients, the variance of the extended τ-leap can be well-behaved, leading to significantly larger step sizes.Conclusions: The benefit of adapting the extended method to the use of RK frameworks is clear in terms of speed of calculation, as the number of evaluations of the Poisson distribution is still one set per time step, as in the original τ-leap method. The approach paves the way to explore new multiscale methods to simulate (bio)chemical systems.

Molecular epidemiology of pseudomonas aeruginosae in intensive care units (ICUs) over a 10-year period (1998-2007)

Pseudomonas aeruginosa, une bactérie environnementale ubiquitaire, est un des pathogènes nosocomiaux les plus fréquents aux soins intensifs. La source de ce microorganisme peut être soit endogène, 2,6 à 24 % des patients hospitalisés étant colonisés au niveau digestif, soit exogène. La proportion des cas d'infections à P. aeruginosa d'origine exogène, donc secondaires à une transmission par manuportage ou par l'eau du réseau utilisée pour la toilette ou d'autres soins, reste débattue. Or une meilleure évaluation du taux d'infections exogènes est importante pour la mise en place de mesures de contrôle appropriées. Le but de cette étude était de déterminer sur une période de 10 ans les rôles respectifs des sources exogènes (robinets, autres patients) et endogène dans la colonisation et/ou l'infection par P.aeruginosa chez les patients des Soins Intensifs, ainsi que de documenter les variations épidémiologiques au cours du temps. L'étude a été menée dans les unités de Soins Intensifs du Centre Hospitalier Universitaire Vaudois (CHUV). Les patients colonisés et/ou infectés par P. aeruginosa entre 1998 et 2007ont été identifiés via la base de données du laboratoire de microbiologie. Ils ont été inclus dans l'étude s'ils étaient hospitalisés dans une des unités de Soins Intensifs, Durant cette période, des prélèvements pour recherche de P. aeruginosa ont été effectués sur des robinets des soins intensifs. Un typage moléculaire a été effectué sur toutes les souches cliniques et environnementales isolées en 1998, 2000, 2003, 2004 et 2007. Les patients inclus dans l'étude ont été répartis en quatre catégories (A-D) selon le résultat du typage moléculaire leur souche de P. aeruginosa. La catégorie A inclut les cas pour lesquels le génotype de P. aeruginosa est identique à un des génotypes retrouvé dans l'environnement. La catégorie B comprend les cas pour lesquels le génotype est identique à celui d'au moins un autre patient. La catégorie C comprend les cas avec un génotype unique et la catégorie D comprend les cas pour lesquels la souche était non disponible pour le typage. Les cas des catégories A et B sont considérés comme ayant une origine exogène. Au cours des années de l'étude, le nombre d'admissions aux soins intensifs est resté stable. En moyenne, 86 patients par année ont été identifiés colonisés ou infectés par P. aeruginosa aux Soins Intensifs. Durant la première année d'investigation, un grand nombre de patients colonisés par une souche de P. aeruginosa identique à une de celles retrouvées dans l'environnement a été mis en évidence. Par la suite, possiblement suite à l'augmentation de la température du réseau d'eau chaude, le nombre de cas dans la catégorie A a diminué. Dans la catégorie B, le nombre de cas varie de 1,9 à 20 cas/1000 admissions selon les années. Ce nombre est supérieur à 10 cas/1000 admissions en 1998, 2003 et 2007 et correspond à des situations épidémiques transitoires. Tout au long des 10 ans de l'étude, le nombre de cas dans la catégorie C (source endogène) est demeuré stable et indépendant des variations du nombre de cas dans les catégories A et B. En conclusion, la contribution relative des réservoirs endogène et exogène dans la colonisation et/ou l'infection des patients de soins Intensifs varie au cours du temps. Les facteurs principaux qui contribuent à de telles variations sont probablement le degré de contamination de l'environnement, la compliance des soignants aux mesures de contrôle des infections et la génétique du pathogène lui-même. Etant donné que ce germe est ubiquitaire dans l'environnement aqueux et colonise jusqu'à 15% des patients hospitalisés, la disparition de son réservoir endogène semble difficile. Cependant, cette étude démontre que son contrôle est possible dans l'environnement, notamment dans les robinets en augmentant la température de l'eau. De plus, si une souche multi-résistante est retrouvée de manière répétée dans l'environnement, des efforts doivent être mis en place pour éliminer cette souche. Des efforts doivent être également entrepris afin de limiter la transmission entre les patients, qui est une cause importante et récurrente de contamination exogène. - Pseudomonas aeruginosa is one of the leading nosocomial pathogens in intensive care units (ICUs). The source of this microorganism can be either endogenous or exogenous. The proportion of cases as a result of transmission is still debated, and its elucidation is important for implementing appropriate control measures. To understand the relative importance of exogenous vs. endogenous sources of P. aeru¬ginosa, molecular typing was performed on all available P. aeruginosa isolated from ICU clinical and environmental specimens in 1998, 2000, 2003, 2004 and 2007. Patient samples were classified according to their P. aeruginosa genotypes into three categories: (A) identical to isolate from faucet; (B) identical to at least one other patient sample and not found in faucet; and (C) unique genotype. Cases in cat¬egories A and Β were considered as possibly exogenous, and cases in category C as possibly endogenous. A mean of 34 cases per 1000 admissions per year were found to be colonized or infected by P. aeruginosa. Higher levels of faucet contamination were correlated with a higher number of cases in category A. The number of cases in category Β varied from 1.9 to 20 cases per 1000 admissions. This num¬ber exceeded 10/1000 admissions on three occasions and was correlated with an outbreak on one occasion. The number of cases con¬sidered as endogenous (category C) was stable and independent of the number of cases in categories A and B. The present study shows that repeated molecular typing can help identify variations in the epidemiology of P. aeruginosa in ICU patients and guide infection control measures.

Quantification of myocardial blood flow with 82Rb positron emission tomography: clinical validation with 15O-water.

PURPOSE: Quantification of myocardial blood flow (MBF) with generator-produced (82)Rb is an attractive alternative for centres without an on-site cyclotron. Our aim was to validate (82)Rb-measured MBF in relation to that measured using (15)O-water, as a tracer 100% of which can be extracted from the circulation even at high flow rates, in healthy control subject and patients with mild coronary artery disease (CAD). METHODS: MBF was measured at rest and during adenosine-induced hyperaemia with (82)Rb and (15)O-water PET in 33 participants (22 control subjects, aged 30 ± 13 years; 11 CAD patients without transmural infarction, aged 60 ± 13 years). A one-tissue compartment (82)Rb model with ventricular spillover correction was used. The (82)Rb flow-dependent extraction rate was derived from (15)O-water measurements in a subset of 11 control subjects. Myocardial flow reserve (MFR) was defined as the hyperaemic/rest MBF. Pearson's correlation r, Bland-Altman 95% limits of agreement (LoA), and Lin's concordance correlation ρ (c) (measuring both precision and accuracy) were used. RESULTS: Over the entire MBF range (0.66-4.7 ml/min/g), concordance was excellent for MBF (r = 0.90, [(82)Rb-(15)O-water] mean difference ± SD = 0.04 ± 0.66 ml/min/g, LoA = -1.26 to 1.33 ml/min/g, ρ(c) = 0.88) and MFR (range 1.79-5.81, r = 0.83, mean difference = 0.14 ± 0.58, LoA = -0.99 to 1.28, ρ(c) = 0.82). Hyperaemic MBF was reduced in CAD patients compared with the subset of 11 control subjects (2.53 ± 0.74 vs. 3.62 ± 0.68 ml/min/g, p = 0.002, for (15)O-water; 2.53 ± 1.01 vs. 3.82 ± 1.21 ml/min/g, p = 0.013, for (82)Rb) and this was paralleled by a lower MFR (2.65 ± 0.62 vs. 3.79 ± 0.98, p = 0.004, for (15)O-water; 2.85 ± 0.91 vs. 3.88 ± 0.91, p = 0.012, for (82)Rb). Myocardial perfusion was homogeneous in 1,114 of 1,122 segments (99.3%) and there were no differences in MBF among the coronary artery territories (p > 0.31). CONCLUSION: Quantification of MBF with (82)Rb with a newly derived correction for the nonlinear extraction function was validated against MBF measured using (15)O-water in control subjects and patients with mild CAD, where it was found to be accurate at high flow rates. (82)Rb-derived MBF estimates seem robust for clinical research, advancing a step further towards its implementation in clinical routine.

Chemistry, isotope values (δD, δ18O, δ34S(SO4)) and temperatures of the water inflows in two Gotthard tunnels, Swiss Alps

Water inflows in the Gotthard Highway Tunnel and in the Gotthard Exploration Tunnel are meteoric waters infiltrating at different elevations, on both sides of an important orographic divide. Limited interaction of meteoric waters with gneissic rocks produces Ca-HCO3 and Na-Ca-HCO3 waters, whereas prolonged interaction of meteoric waters with the same rocks generates Na-HCO3 to Na-SO4 waters. Waters circulating in Triassic carbonate-evaporite rocks have a Ca-SO4 composition. Calcium-Na-SO4 waters are also present. They can be produced through interaction of either Na-HCO3 waters with anhydrite or Ca-SO4 waters with a local gneissic rock, as suggested by reaction path modeling. An analogous simulation indicates that Na-HCO3 waters are generated through interaction of Ca-HCO3 waters with a local gneissic rock. The two main SO4-sources present in the Alps are leaching of upper Triassic sulfate minerals and oxidative dissolution of sulfide minerals of crystalline rocks. Values of delta S-34(SO4) < <similar to>+ 9 parts per thousand, are due to oxidative dissolution of sulfide minerals, whereas delta S-34(SO4) > similar to+ 9 parts per thousand are controlled either by bacterial SO4 reduction or leaching of upper Triassic sulfate minerals. Most waters have temperatures similar to the expected values for a geothermal gradient of 22 degreesC/km and are close to thermal equilibrium with rocks. However relatively large, descending flows of cold waters and ascending flows of warm waters are present in both tunnels and determine substantial cooling and heating, respectively, of the interacting rocks. The most import upflow zone of warm, Na-rich waters is below Guspisbach, in the Gotthard Highway Tunnel, at 6.2-9.0 km from the southern portal. These warm waters have equilibrium temperatures of 65-75 degreesC and therefore constitute an important low-enthalpy geothermal resource. (C) 2001 Elsevier Science Ltd. All rights reserved.

Molecular and cellular analysis of genodermatoses

Summary Skin is the essential interface between our body and its environment; not only does it prevent water loss and protect us from external insults it also plays an essential role in the central nervous system acting as a major sense organ primarily for touch and pain. The main cell type present in skin, keratinocyte, undergoes a differentiation process leading to the formation of this protecting barrier. This work is intended to contribute to the understanding of how keratinocyte differentiates and skin functions. To do this, we studied two genetic skin diseases: Erythrokeratodermia variabilis and Mal de Meleda. Our approach was to examine the expression and localization of proteins implicated in these two pathologies in normal and diseased tissues and to determine the influence of mutant proteins at the molecular and cellular levels. Connexins are major components of gap junctions, channels allowing direct communication between cells. Our laboratory has identified mutations in both connexin 30.3 (Cx30.3) and 31 (Cx31) to be causally involved in erythrokeratodermia variabilis (EKV), an autosomal dominant disorder of keratinization. In the first chapter, we show a new mutation of Cx31, L209P-Cx31, in 3 EKV patients, extending the field of EKV-causing mutations although the mechanism by which connexin mutations lead to the disease is unclear. In the second chapter, we studied the effect of F137L-Cx30.3 on expression, trafficking and localization of cotransfected Cx31 and Cx30.3 in connexin-deficient HeLa cells. The F137 amino acid, highly conserved in connexin family, is oriented towards the channel pore and F137L mutation in either Cx30.3 or Cx31 lead to EKV. As two genes can lead to EKV when mutated, our hypothesis was that Cx31 and Cx30.3 might cooperate at a molecular level. We were able to demonstrate a physical interaction between Cx31 and Cx30.3. The presence of F137L-Cx30.3 disturbed the trafficking of both connexins, less connexins were integrated into gap junctions and thus, the coupling between cell was diminished. Connexins formed in the presence of F137L-Cx30.3 are degraded at their exit from the endoplasmic reticulum. In conclusion, our results indicate that the genetic heterogeneity of EKV is due to mutations in two interacting proteins. F137L-Cx30.3 has a dominant negative effect and affects Cx31, disturbing cellular communication in epidermal cells. Mal de Meleda is an autosomal recessive inflammatory and a keratotic palmoplantar skin disorder due to mutations in SLURP1 (secreted LY6/PLAUR-related protein 1). SLURP1 belongs to the LY6/PLAUR family of proteins and has the particularity of being secreted instead of being GPI-anchored. The high degree of structural similarity between SLURP1 and the three fingers motif of snake neurotoxins and LYNX 1-C suggests that this protein could interact with the neuronal acetylcholine receptors. In the third chapter, we show that SLURP1 potentiates responses of the a7 nicotinic acetylcholine receptor (nAchR) to acetylcholine. These results identify SLURP1 as a secreted epidermal neuromodulator that is likely to be essential for palmoplantar skin. In the fourth chapter, we show that SLURP1 is expressed in the granular layer of the epidermis but is absent from skin biopsies of Mal de Meleda patients. SLURP1 is also present in secretions such as sweat, tears or saliva. An in vitro analysis on two mutant of SLURP-I demonstrates that W15R-SLURP1 is absent in cells while G86R-SLURP1 is expressed and secreted, suggesting that SLURP1 can lead to the disease by either an absent or an abnormal protein. Finally, in the fifth chapter, we analyse the expression and biological properties of other LY6/PLAUR members, clustered around SLURP] on chromosome 8. Their GPI-anchored or secreted status were analysed in vitro. SLURP1, LYNX1-A and -B are secreted while LYPDC2 and LYNX 1-C are GPI anchored. Three of these proteins are expressed in the epidermis and in cultured keratinocytes. These results suggest that these LY6/PLAUR members may have an important role in skin homeostasis. Résumé Résumé La peau est la barrière essentielle entre notre corps et l'environnement, nous protégeant des agressions extérieures, de la déshydratation et assurant aussi un rôle dans le système nerveux central en tant qu'organe du toucher et de la douleur. Le principal type de cellules présent dans la peau est le kératinocyte qui suit un processus de différenciation aboutissant à la formation de cette barrière protectrice. Ce travail est destiné à comprendre la différenciation des kératinocytes et le fonctionnement de la peau. Pour cela, nous avons étudié deux maladies génodermatoses : l'Erthrokeratodermia Variabilis (EKV) et le Mal de Meleda. Nous avons examiné l'expression et la localisation des protéines impliquées dans ces deux pathologies dans des tissus normaux et malades puis déterminé l'influence des protéines mutantes aux niveaux moléculaires et cellulaires. Les connexines (Cx) sont les composants majeurs des jonctions communicantes, canaux permettant la communication directe entre les cellules. Notre laboratoire a identifié des mutations dans les Cx30.3 et Cx31 comme responsables de l'EKV, génodermatose de transmission autosomique dominante. Dans le ler chapitre, nous décrivons une nouvelle mutation de Cx31, L209-Cx31, et contribuons à l'établissement du catalogue des mutations de Cx31 entraînant cette maladie. Cependant, le mécanisme par lequel les mutations de Cx31 et C3x0.3 provoquent l'EKV est inconnu. Dans le 2ème chapitre, nous étudions les effets de la mutation F137L-Cx30.3 sur l'expression, le trafic et la localisation des Cx31 et Cx30.3 transfectées dans des cellules HeLa, déficientes en connexines. Comme deux gènes peuvent causer une EKV quand ils sont mutés, notre hypothèse était que Cx31 et Cx30.3 pourraient coopérer au niveau moléculaire. Nous avons montré l'existence d'une interaction physique entre ces deux connexines. La présence de la mutation F137L-Cx30.3 perturbe le trafic des deux connexines, moins de connexines sont intégrées dans les jonctions communicantes et donc le couplage entre les cellules est diminué. Les connexons formés en présence de cette mutation sont dégradés à leur sortie du réticulum endoplasmique. En conclusion, nos résultats indiquent que l'hétérogénéité génétique de EKV est due à des mutations dans deux protéines qui interagissent. F137L-Cx30.3 a un effet dominant négatif et affecte Cx31, perturbant la communication entre les cellules épidermiques. Le Mal de Meleda est une maladie récessive de la peau palmoplantaire due à des mutations dans SLURP1. SLURP1 appartient à la famille des protéines contenant un domaine LY6/PLAUR et a la particularité d'être sécrétée. La grande homologie de structure existant entre SLURP1, les neurotoxines de serpent et LYNX1-C suggère que la protéine pourrait interagir avec des récepteurs à acétylcholine (Ach). Dans le 3ème chapitre, nous montrons que SLURP1 module la réponse à l'Ach du récepteur nicotinique α7. Ces résultats identifient SLURP1 comme un neuromodulateur épidermique sécrété, probablement essentiel pour la peau palmoplantaire. Dans le 4ème chapitre, nous montrons que SLURP1 est exprimé dans la couche granuleuse de l'épiderme et qu'il est absent des biopsies des patients. SLURP1 a aussi été détecté dans des sécrétions telles que la sueur, les lamies et la salive. Une analyse in vitro de deux mutants de SLURP1 a montré que W15R-SLURP1 est absent des cellules tandis que G86R-SLURP1 est exprimé et sécrété, suggérant qu'une absence ou une anomalie de SLURP1 peuvent causer la maladie. Finalement, dans le 5ème chapitre, nous analysons l'expression et les propriétés biologiques d'autres membres de la famille LY6/PLAUR localisés autour de SLURP1 sur le chromosome 8. Leur statut de protéines sécrétées ou liées à la membrane par une ancre GPI est analysé in vitro. SLURP1, LYNXI-A et -B sont sécrétées alors que LYPDC2 et LYNX1-C sont liés à la membrane. Trois de ces protéines sont exprimées dans l'épiderme et dans des kératinocytes cultivés. Ces résultats suggèrent que la famille LY6/PLAUR pourrait avoir un rôle important dans l'homéostasie de la peau.

iLOGP: A Simple, Robust, and Efficient Description of n-Octanol/Water Partition Coefficient for Drug Design Using the GB/SA Approach.

The n-octanol/water partition coefficient (log Po/w) is a key physicochemical parameter for drug discovery, design, and development. Here, we present a physics-based approach that shows a strong linear correlation between the computed solvation free energy in implicit solvents and the experimental log Po/w on a cleansed data set of more than 17,500 molecules. After internal validation by five-fold cross-validation and data randomization, the predictive power of the most interesting multiple linear model, based on two GB/SA parameters solely, was tested on two different external sets of molecules. On the Martel druglike test set, the predictive power of the best model (N = 706, r = 0.64, MAE = 1.18, and RMSE = 1.40) is similar to six well-established empirical methods. On the 17-drug test set, our model outperformed all compared empirical methodologies (N = 17, r = 0.94, MAE = 0.38, and RMSE = 0.52). The physical basis of our original GB/SA approach together with its predictive capacity, computational efficiency (1 to 2 s per molecule), and tridimensional molecular graphics capability lay the foundations for a promising predictor, the implicit log P method (iLOGP), to complement the portfolio of drug design tools developed and provided by the SIB Swiss Institute of Bioinformatics.

Amoeba-resisting Chlamydiales present in domestic drinking water and their potential role in pneumonia

Chlamydia-related bacteria classified in the Chlamydiales order, are strictly intracellular bacteria and are able for the most to replicate in free-living amoebae. Amoebae, ubiquitous in the environment and especially in water, are very resistant to disinfection used in drinking water production. Thus, amoebae may reach easily the distribution and domestic water system, potentially sheltering amoeba-resisting bacteria including Legionella, mycobacteria and Chlamydiales. Indeed, some of these amoeba-resisting bacteria have been shown to cause respiratory infections in people inhaling contaminated water. Therefore, an environmental and clinical study was conducted to determine if Chlamydiales bacteria are also involved in respiratory infections and if a transmission through domestic drinking water could occur. First, large scale molecular and serological tools specific of Chlamydia-related bacteria were developed and then were applied on clinical samples from patients with and without pneumonia. Simultaneously, water and biofilm samples from households of the same patients were investigated using molecular and culture methods for the presence of Chlamydiales bacteria. Chlamydiales were detected in the nasopharyngeal flora from patients with and without pneumonia. However, no significant difference was observed between both groups. Conversely, serological investigations showed that antibody reactivity against members of the Criblamydiaceae was associated with pneumonia. The thesis provided very efficient tools that showed the presence of Chlamydiales in human nasopharyngeal flora as well as in the majority of the domestic drinking water. However, no transmission from domestic drinking water to human could be demonstrated. These tools will help in the future specifying the ecology and pathogenicity of the Chlamydia-re\ated bacteria and especially of the species belonging to the Criblamydiaceae family.

Molecular determinants of desensitization in an ENaC/degenerin channel.

Epithelial Na(+) channel (ENaC)/degenerin family members are involved in mechanosensation, blood pressure control, pain sensation, and the expression of fear. Several of these channel types display a form of desensitization that allows the channel to limit Na(+) influx during prolonged stimulation. We used site-directed mutagenesis and chemical modification, functional analysis, and molecular dynamics simulations to investigate the role of the lower palm domain of the acid-sensing ion channel 1, a member of the ENaC/degenerin family. The lower palm domains of this trimeric channel are arranged around a central vestibule, at ∼20 Å above the plasma membrane and are covalently linked to the transmembrane channel parts. We show that the lower palm domains approach one another during desensitization. Residues in the palm co-determine the pH dependence of desensitization, its kinetics, and the stability of the desensitized state. Mutations of palm residues impair desensitization by preventing the closing movement of the palm. Overexpression of desensitization-impaired channel mutants in central neurons allowed--in contrast to overexpression of wild type--a sustained signaling response to rapid pH fluctuations. We identify and describe here the function of an important regulatory domain that most likely has a conserved role in ENaC/degenerin channels.

Molecular and isotopic characterization of organic samples from the wreck of the Saint-Etienne merchant ship (XVIIIth century): Identification of pitch, fat, hair and sulfur

The wreck U Pezzo, excavated within the Saint Florent Gulf in northern Corsica was identified as the pink, Saint Etienne, a merchant ship which sank on January 31, 1769. In order to determine the composition of organic materials used to coat the hull or to waterproof different parts of the pink, a study of several samples, using molecular biomarker and carbon isotopic analysis, was initiated. The results revealed that the remarkable yellow coat, covering the outside planks of the ship's bottom under the water line, is composed of sulfur, tallow (of ox and not of cetacean origin) and black pitch which corresponds to a mixture called ``couroi'' or ``stuff'. Onboard ropes had been submitted to a tarring treatment with pitch. Hairs mixed with pitch were identified in samples collected between the two layers of the hull or under the sheathing planking. The study also provides a key model for weathering of pitch, as different degrees of degradation were found between the surface and the heart of several samples. Accordingly, molecular parameters for alteration were proposed. Furthermore novel mixed esters between terpenic and diterpenic alcohols and the free major fatty acids (C(14:0), C(16:0), C(18:0)) were detected in the yellow coat. (C) 2009 Elsevier Ltd. All rights reserved.

Direct simulation of interface dynamics: linking capillary pressure, interfacial area and surface energy

We perform direct numerical simulations of drainage by solving Navier- Stokes equations in the pore space and employing the Volume Of Fluid (VOF) method to track the evolution of the fluid-fluid interface. After demonstrating that the method is able to deal with large viscosity contrasts and to model the transition from stable flow to viscous fingering, we focus on the definition of macroscopic capillary pressure. When the fluids are at rest, the difference between inlet and outlet pressures and the difference between the intrinsic phase average pressure coincide with the capillary pressure. However, when the fluids are in motion these quantities are dominated by viscous forces. In this case, only a definition based on the variation of the interfacial energy provides an accurate measure of the macroscopic capillary pressure and allows separating the viscous from the capillary pressure components.

Homeostatic and circadian contribution to EEG and molecular state variables of sleep regulation.

STUDY OBJECTIVES: Besides their well-established role in circadian rhythms, our findings that the forebrain expression of the clock-genes Per2 and Dbp increases and decreases, respectively, in relation to time spent awake suggest they also play a role in the homeostatic aspect of sleep regulation. Here, we determined whether time of day modulates the effects of elevated sleep pressure on clock-gene expression. Time of day effects were assessed also for recognized electrophysiological (EEG delta power) and molecular (Homer1a) markers of sleep homeostasis. DESIGN: EEG and qPCR data were obtained for baseline and recovery from 6-h sleep deprivation starting at ZT0, -6, -12, or -18. SETTING: Mouse sleep laboratory. PARTICIPANTS: Male mice. INTERVENTIONS: Sleep deprivation. RESULTS: The sleep-deprivation induced changes in Per2 and Dbp expression importantly varied with time of day, such that Per2 could even decrease during sleep deprivations occurring at the decreasing phase in baseline. Dbp showed similar, albeit opposite dynamics. These unexpected results could be reliably predicted assuming that these transcripts behave according to a driven damped harmonic oscillator. As expected, the sleep-wake distribution accounted for a large degree of the changes in EEG delta power and Homer1a. Nevertheless, the sleep deprivation-induced increase in delta power varied also with time of day with higher than expected levels when recovery sleep started at dark onset. CONCLUSIONS: Per2 and delta power are widely used as exclusive state variables of the circadian and homeostatic process, respectively. Our findings demonstrate a considerable cross-talk between these two processes. As Per2 in the brain responds to both sleep loss and time of day, this molecule is well positioned to keep track of and to anticipate homeostatic sleep need. CITATION: Curie T; Mongrain V; Dorsaz S; Mang GM; Emmenegger Y; Franken P. Homeostatic and circadian contribution to EEG and molecular state variables of sleep regulation. SLEEP 2013;36(3):311-323.

Brownian dynamics simulation of DNA condensation.

DNA condensation observed in vitro with the addition of polyvalent counterions is due to intermolecular attractive forces. We introduce a quantitative model of these forces in a Brownian dynamics simulation in addition to a standard mean-field Poisson-Boltzmann repulsion. The comparison of a theoretical value of the effective diameter calculated from the second virial coefficient in cylindrical geometry with some experimental results allows a quantitative evaluation of the one-parameter attractive potential. We show afterward that with a sufficient concentration of divalent salt (typically approximately 20 mM MgCl(2)), supercoiled DNA adopts a collapsed form where opposing segments of interwound regions present zones of lateral contact. However, under the same conditions the same plasmid without torsional stress does not collapse. The condensed molecules present coexisting open and collapsed plectonemic regions. Furthermore, simulations show that circular DNA in 50% methanol solutions with 20 mM MgCl(2) aggregates without the requirement of torsional energy. This confirms known experimental results. Finally, a simulated DNA molecule confined in a box of variable size also presents some local collapsed zones in 20 mM MgCl(2) above a critical concentration of the DNA. Conformational entropy reduction obtained either by supercoiling or by confinement seems thus to play a crucial role in all forms of condensation of DNA.

A novel human aquaporin-4 splice variant exhibits a dominant-negative activity: a new mechanism to regulate water permeability.

Two major isoforms of aquaporin-4 (AQP4) have been described in human tissue. Here we report the identification and functional analysis of an alternatively spliced transcript of human AQP4, AQP4-Δ4, that lacks exon 4. In transfected cells AQP4-Δ4 is mainly retained in the endoplasmic reticulum and shows no water transport properties. When AQP4-Δ4 is transfected into cells stably expressing functional AQP4, the surface expression of the full-length protein is reduced. Furthermore, the water transport activity of the cotransfectants is diminished in comparison to transfectants expressing only AQP4. The observed down-regulation of both the expression and water channel activity of AQP4 is likely to originate from a dominant-negative effect caused by heterodimerization between AQP4 and AQP4-Δ4, which was detected in coimmunoprecipitation studies. In skeletal muscles, AQP4-Δ4 mRNA expression inversely correlates with the level of AQP4 protein and is physiologically associated with different types of skeletal muscles. The expression of AQP4-Δ4 may represent a new regulatory mechanism through which the cell-surface expression and therefore the activity of AQP4 can be physiologically modulated.