Combined simulation and mutagenesis analyses reveal the involvement of key residues for peroxisome proliferator-activated receptor alpha helix 12 dynamic behavior.

The dynamic properties of helix 12 in the ligand binding domain of nuclear receptors are a major determinant of AF-2 domain activity. We investigated the molecular and structural basis of helix 12 mobility, as well as the involvement of individual residues with regard to peroxisome proliferator-activated receptor alpha (PPARalpha) constitutive and ligand-dependent transcriptional activity. Functional assays of the activity of PPARalpha helix 12 mutants were combined with free energy molecular dynamics simulations. The agreement between the results from these approaches allows us to make robust claims concerning the mechanisms that govern helix 12 functions. Our data support a model in which PPARalpha helix 12 transiently adopts a relatively stable active conformation even in the absence of a ligand. This conformation provides the interface for the recruitment of a coactivator and results in constitutive activity. The receptor agonists stabilize this conformation and increase PPARalpha transcription activation potential. Finally, we disclose important functions of residues in PPARalpha AF-2, which determine the positioning of helix 12 in the active conformation in the absence of a ligand. Substitution of these residues suppresses PPARalpha constitutive activity, without changing PPARalpha ligand-dependent activation potential.

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.

Direct numerical simulations of interface dynamics to link capillary pressure and total surface energy

The flow of two immiscible fluids through a porous medium depends on the complex interplay between gravity, capillarity, and viscous forces. The interaction between these forces and the geometry of the medium gives rise to a variety of complex flow regimes that are difficult to describe using continuum models. Although a number of pore-scale models have been employed, a careful investigation of the macroscopic effects of pore-scale processes requires methods based on conservation principles in order to reduce the number of modeling assumptions. In this work 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 model the transition from stable flow to viscous fingering, we focus on the macroscopic capillary pressure and we compare different definitions of this quantity under quasi-static and dynamic conditions. We show that the difference between the intrinsic phase-average pressures, which is commonly used as definition of Darcy-scale capillary pressure, is subject to several limitations and it is not accurate in presence of viscous effects or trapping. In contrast, a definition based on the variation of the total surface energy provides an accurate estimate of the macroscopic capillary pressure. This definition, which links the capillary pressure to its physical origin, allows a better separation of viscous effects and does not depend on the presence of trapped fluid clusters.

Protein-protein interaction investigated by steered molecular dynamics: the TCR-pMHC complex.

We present a novel steered molecular dynamics scheme to induce the dissociation of large protein-protein complexes. We apply this scheme to study the interaction of a T cell receptor (TCR) with a major histocompatibility complex (MHC) presenting a peptide (p). Two TCR-pMHC complexes are considered, which only differ by the mutation of a single amino acid on the peptide; one is a strong agonist that produces T cell activation in vivo, while the other is an antagonist. We investigate the interaction mechanism from a large number of unbinding trajectories by analyzing van der Waals and electrostatic interactions and by computing energy changes in proteins and solvent. In addition, dissociation potentials of mean force are calculated with the Jarzynski identity, using an averaging method developed for our steering scheme. We analyze the convergence of the Jarzynski exponential average, which is hampered by the large amount of dissipative work involved and the complexity of the system. The resulting dissociation free energies largely underestimate experimental values, but the simulations are able to clearly differentiate between wild-type and mutated TCR-pMHC and give insights into the dissociation mechanism.

Application of a roughness-length representation to parameterize energy loss in 3-D numerical simulations of large rivers

Recent technological advances in remote sensing have enabled investigation of the morphodynamics and hydrodynamics of large rivers. However, measuring topography and flow in these very large rivers is time consuming and thus often constrains the spatial resolution and reach-length scales that can be monitored. Similar constraints exist for computational fluid dynamics (CFD) studies of large rivers, requiring maximization of mesh-or grid-cell dimensions and implying a reduction in the representation of bedform-roughness elements that are of the order of a model grid cell or less, even if they are represented in available topographic data. These ``subgrid'' elements must be parameterized, and this paper applies and considers the impact of roughness-length treatments that include the effect of bed roughness due to ``unmeasured'' topography. CFD predictions were found to be sensitive to the roughness-length specification. Model optimization was based on acoustic Doppler current profiler measurements and estimates of the water surface slope for a variety of roughness lengths. This proved difficult as the metrics used to assess optimal model performance diverged due to the effects of large bedforms that are not well parameterized in roughness-length treatments. However, the general spatial flow patterns are effectively predicted by the model. Changes in roughness length were shown to have a major impact upon flow routing at the channel scale. The results also indicate an absence of secondary flow circulation cells in the reached studied, and suggest simpler two-dimensional models may have great utility in the investigation of flow within large rivers. Citation: Sandbach, S. D. et al. (2012), Application of a roughness-length representation to parameterize energy loss in 3-D numerical simulations of large rivers, Water Resour. Res., 48, W12501, doi: 10.1029/2011WR011284.

Glutamatergic regulation of mitochondrial ion dynamics in astrocytes

Résumé grand public :Le cerveau se compose de cellules nerveuses appelées neurones et de cellules gliales dont font partie les astrocytes. Les neurones communiquent entre eux par signaux électriques et en libérant des molécules de signalisation comme le glutamate. Les astrocytes ont eux pour charge de capter le glucose depuis le sang circulant dans les vaisseaux sanguins, de le transformer et de le transmettre aux neurones pour qu'ils puissent l'utiliser comme source d'énergie. L'astrocyte peut ensuite utiliser ce glucose de deux façons différentes pour produire de l'énergie : la première s'opère dans des structures appelées mitochondries qui sont capables de produire plus de trente molécules riches en énergie (ATP) à partir d'une seule molécule de glucose ; la seconde possibilité appelée glycolyse peut produire deux molécules d'ATP et un dérivé du glucose appelé lactate. Une théorie couramment débattue propose que lorsque les astrocytes capturent le glutamate libéré par les neurones, ils libèrent en réponse du lactate qui servirait de base énergétique aux neurones. Cependant, ce mécanisme n'envisage pas une augmentation de l'activité des mitochondries des astrocytes, ce qui serait pourtant bien plus efficace pour produire de l'énergie.En utilisant la microscopie par fluorescence, nous avons pu mesurer les changements de concentrations ioniques dans les mitochondries d'astrocytes soumis à une stimulation glutamatergique. Nous avons démontré que les mitochondries des astrocytes manifestent des augmentations spontanées et transitoires de leur concentrations ioniques, dont la fréquence était diminuée au cours d'une stimulation avec du glutamate. Nous avons ensuite montré que la capture de glutamate augmentait la concentration en sodium et acidifiait les mitochondries des astrocytes. En approfondissant ces mécanismes, plusieurs éléments ont suggéré que l'acidification induite diminuerait le potentiel de synthèse d'énergie d'origine mitochondriale et la consommation d'oxygène dans les astrocytes. En résumé, l'ensemble de ces travaux suggère que la signalisation neuronale impliquant le glutamate dicte aux astrocytes de sacrifier temporairement l'efficacité de leur métabolisme énergétique, en diminuant l'activité de leurs mitochondries, afin d'augmenter la disponibilité des ressources énergétiques utiles aux neurones.Résumé :La remarquable efficacité du cerveau à compiler et propager des informations coûte au corps humain 20% de son budget énergétique total. Par conséquent, les mécanismes cellulaires responsables du métabolisme énergétique cérébral se sont adéquatement développés pour répondre aux besoins énergétiques du cerveau. Les dernières découvertes en neuroénergétique tendent à démontrer que le site principal de consommation d'énergie dans le cerveau est situé dans les processus astrocytaires qui entourent les synapses excitatrices. Un nombre croissant de preuves scientifiques a maintenant montré que le transport astrocytaire de glutamate est responsable d'un coût métabolique important qui est majoritairement pris en charge par une augmentation de l'activité glycolytique. Cependant, les astrocytes possèdent également un important métabolisme énergétique de type mitochondrial. Par conséquent, la localisation spatiale des mitochondries à proximité des transporteurs de glutamate suggère l'existence d'un mécanisme régulant le métabolisme énergétique astrocytaire, en particulier le métabolisme mitochondrial.Afin de fournir une explication à ce paradoxe énergétique, nous avons utilisé des techniques d'imagerie par fluorescence pour mesurer les modifications de concentrations ioniques spontanées et évoquées par une stimulation glutamatergique dans des astrocytes corticaux de souris. Nous avons montré que les mitochondries d'astrocytes au repos manifestaient des changements individuels, spontanés et sélectifs de leur potentiel électrique, de leur pH et de leur concentration en sodium. Nous avons trouvé que le glutamate diminuait la fréquence des augmentations spontanées de sodium en diminuant le niveau cellulaire d'ATP. Nous avons ensuite étudié la possibilité d'une régulation du métabolisme mitochondrial astrocytaire par le glutamate. Nous avons montré que le glutamate initie dans la population mitochondriale une augmentation rapide de la concentration en sodium due à l'augmentation cytosolique de sodium. Nous avons également montré que le relâchement neuronal de glutamate induit une acidification mitochondriale dans les astrocytes. Nos résultats ont indiqué que l'acidification induite par le glutamate induit une diminution de la production de radicaux libres et de la consommation d'oxygène par les astrocytes. Ces études ont montré que les mitochondries des astrocytes sont régulées individuellement et adaptent leur activité selon l'environnement intracellulaire. L'adaptation dynamique du métabolisme énergétique mitochondrial opéré par le glutamate permet d'augmenter la quantité d'oxygène disponible et amène au relâchement de lactate, tous deux bénéfiques pour les neurones.Abstract :The remarkable efficiency of the brain to compute and communicate information costs the body 20% of its total energy budget. Therefore, the cellular mechanisms responsible for brain energy metabolism developed adequately to face the energy needs. Recent advances in neuroenergetics tend to indicate that the main site of energy consumption in the brain is the astroglial process ensheating activated excitatory synapses. A large body of evidence has now shown that glutamate uptake by astrocytes surrounding synapses is responsible for a significant metabolic cost, whose metabolic response is apparently mainly glycolytic. However, astrocytes have also a significant mitochondrial oxidative metabolism. Therefore, the location of mitochondria close to glutamate transporters raises the question of the existence of mechanisms for tuning their energy metabolism, in particular their mitochondrial metabolism.To tackle these issues, we used real time imaging techniques to study mitochondrial ionic alterations occurring at resting state and during glutamatergic stimulation of mouse cortical astrocytes. We showed that mitochondria of intact resting astrocytes exhibited individual spontaneous and selective alterations of their electrical potential, pH and Na+ concentration. We found that glutamate decreased the frequency of mitochondrial Na+ transient activity by decreasing the cellular level of ATP. We then investigated a possible link between glutamatergic transmission and mitochondrial metabolism in astrocytes. We showed that glutamate triggered a rapid Na+ concentration increase in the mitochondrial population as a result of plasma-membrane Na+-dependent uptake. We then demonstrated that neuronally released glutamate also induced a mitochondrial acidification in astrocytes. Glutamate induced a pH-mediated and cytoprotective decrease of mitochondrial metabolism that diminished oxygen consumption. Taken together, these studies showed that astrocytes contain mitochondria that are individually regulated and sense the intracellular environment to modulate their own activity. The dynamic regulation of astrocyte mitochondrial energy output operated by glutamate allows increasing oxygen availability and lactate production both being beneficial for neurons.

NO rebinding to myoglobin: a reactive molecular dynamics study.

The rebinding of NO to myoglobin after photolysis is studied using the 'reactive molecular dynamics' method. In this approach the energy of the system is evaluated on two potential energy surfaces that include the heme-ligand interactions which change between liganded and unliganded myoglobin. This makes it possible to take into account in a simple way, the high dimensionality of the transition seam connecting the reactant and product states. The dynamics of the dissociated NO molecules are examined, and the geometrical and energetic properties of the transition seam are studied. Analysis of the frequency of recrossing shows that the height of the effective rebinding barrier is dependent on the time after photodissociation. This effect is due mainly to protein relaxation and may contribute to the experimentally observed non-exponential rebinding rate of NO, as has been suggested previously.

Gender and weight shape brain dynamics during food viewing.

Hemodynamic imaging results have associated both gender and body weight to variation in brain responses to food-related information. However, the spatio-temporal brain dynamics of gender-related and weight-wise modulations in food discrimination still remain to be elucidated. We analyzed visual evoked potentials (VEPs) while normal-weighted men (n = 12) and women (n = 12) categorized photographs of energy-dense foods and non-food kitchen utensils. VEP analyses showed that food categorization is influenced by gender as early as 170 ms after image onset. Moreover, the female VEP pattern to food categorization co-varied with participants' body weight. Estimations of the neural generator activity over the time interval of VEP modulations (i.e. by means of a distributed linear inverse solution [LAURA]) revealed alterations in prefrontal and temporo-parietal source activity as a function of image category and participants' gender. However, only neural source activity for female responses during food viewing was negatively correlated with body-mass index (BMI) over the respective time interval. Women showed decreased neural source activity particularly in ventral prefrontal brain regions when viewing food, but not non-food objects, while no such associations were apparent in male responses to food and non-food viewing. Our study thus indicates that gender influences are already apparent during initial stages of food-related object categorization, with small variations in body weight modulating electrophysiological responses especially in women and in brain areas implicated in food reward valuation and intake control. These findings extend recent reports on prefrontal reward and control circuit responsiveness to food cues and the potential role of this reactivity pattern in the susceptibility to weight gain.

Integrating receptor interactions and dynamics and Interference with endocrine disruptors in the mechanisms of action of PPAR nuclear receptors

Abstract Peroxisome Proliferator-Activated Receptors (PPARs) form a family of three nuclear receptors regulating important cellular and metabolic functions. PPARs control gene expression by directly binding to target promoters as heterodimers with the Retinoid X Receptor (RXR), and their transcriptional activity is enhanced upon activation by natural or pharmacological ligands. The binding of PPAR/RXR heterodimers on target promoters allows the anchoring of a series of coactivators and corepressors involved in promoter remodeling and the recruitment of the transcription machinery. The transcriptional output finally depends on a complex interplay between (i) the respective expression levels of PPARs, RXRs and of other nuclear receptors competing for DNA binding and RXR recruitment, (ii) the availability and the nature of PPAR and RXR ligands, (iii) the expression levels and the nature of the different coactivators and corepressors and (iv) the sequence and the epigenetic status of the promoter. Understanding how all these factors and signals integrate and fine-tune transcription remains a challenge but is necessary to understand the specificity of the physiological functions regulated by PPARs. The work presented herein focuses on the molecular mechanisms of PPAR action and aims at understanding how the interactions and mobility of the receptor modulate transcription in the physiological context of a living cell: Such observations in vivo rely on the use of engineered fluorescent protein chimeras and require the development and the application of complementary imaging techniques such as Fluorescence Recovery After Photobleaching (FRAP), Fluorescence Resonance Energy Transfer (FRET) and Fluorescence Correlation Spectroscopy (FCS). Using such techniques, PPARs are shown to reside solely in the nucleus where they are constitutively associated with RXR but transcriptional activation by ligand binding -does not promote the formation of sub-nuclear structures as observed with other nuclear receptors. In addition, the engagement of unliganded PPARs in large complexes of cofactors in living cells provides a molecular basis for their ligand-independent activity. Ligand binding reduces receptor diffusion by promoting the recruitment of coactivators which further enlarge the size of PPAR complexes to acquire full transcriptional competence. Using these molecular approaches, we deciphered the molecular mechanisms through which phthalates, a class of pollutants from the plastic industry, interfere with PPARγ signaling. Mono-ethyl-hexyl-phthalate (MEHP) binding induces the recruitment of a specific subset of cofactors and translates into the expression of a specific subset of target genes, the transcriptional output being strongly conditioned by the differentiation status of the cell. This selective PPARγ modulation induces limited adipogenic effects in cellular models while exposure to phthalates in animal models leads to protective effects on glucose tolerance and diet-induced obesity. These results demonstrate that phthalates influence lipid and carbohydrate metabolism through complex mechanisms which most likely involve PPARγ but also probably PPARα and PPARß, Altogether, the molecular and physiological demonstration of the interference of pollutants with PPAR action outlines an important role of chemical exposure in metabolic regulations. Résumé Les PPARs (Peroxisome Proliferator-Activated Receptors) forment une famille de récepteurs nucléaires qui régulent des fonctions cellulaires et métaboliques importantes. Les PPARs contrôlent l'expression des gènes en se liant directement à leurs promoteurs sous forme d'hétérodimères avec les récepteurs RXR (Retinoid X Receptor), et leur activité transcriptionnelle est stimulée par la liaison de ligands naturels ou pharmacologiques. L'association des hétérodimères PPAR/RXR avec les promoteurs des gènes cibles permet le recrutement de coactivateurs et de corépresseurs qui vont permettre le remodelage de la chromatine et le recrutement de la machinerie transcriptionnelle. Les actions transcriptionnelles du récepteur dépendent toutefois d'interactions complexes qui sont régulées par (i) le niveau d'expression des PPARs, des RXRs et d'autres récepteurs nucléaires entrant en compétition pour la liaison à l'ADN et l'association avec RXR, (ii) la disponibilité et la nature de ligands de PPAR et de RXR, (iii) les niveaux d'expression et la nature des différents coactivateurs et corépresseurs et (iv) la séquence et le marquage épigénétique des promoteurs. La compréhension des mécanismes qui permettent d'intégrer ces aspects pour assurer une régulation fine de l'activité transcriptionnelle est un défi qu'il est nécessaire de relever pour comprendre la spécificité des fonctions physiologiques régulées par les PPARs. Ce travail concerne l'étude des mécanismes d'action moléculaire des PPARs et vise à mieux comprendre comment les interactions du récepteur avec d'autres protéines ainsi que la mobilité de ce dernier régulent son activité transcriptionnelle dans le contexte physiologique des cellules vivantes. De telles observations reposent sur l'emploi de protéines fusionnées à des protéines fluorescentes ainsi que sur le développement et l'utilisation de techniques d'imagerie complémentaires telles que le FRAP (Fluorescence Recovery After Photobleaching), le FRET (Fluorescence Resonance Energy Transfer) ou la FCS (Fluorescence Corrélation Spectroscopy). En appliquant ces méthodes, nous avons pu montrer que les PPARs résident toujours dans le noyau où ils sont associés de manière constitutive à RXR, mais que l'ajout de ligand n'induit pas la formation de structures sub-nucléaires comme cela a pu être décrit pour d'autres récepteurs nucléaires. De plus, les PPARs sont engagés dans de larges complexes protéiques de cofacteurs en absence de ligand, ce qui procure une explication moléculaire à leur activité ligand-indépendante. La liaison du ligand réduit la vitesse de diffusion du récepteur en induisant le recrutement de coactivateurs qui augmente encore plus la taille des complexes afin d'acquérir un potentiel d'activation maximal. En utilisant ces approches moléculaires, nous avons pu caractériser les mécanismes permettant aux phtalates, une classe de polluants provenant de l'industrie plastique, d'interférer avec PPARγ. La liaison du mono-ethyl-hexyl-phtalate (NERF) à PPARγ induit un recrutement sélectif de cofacteurs, se traduisant par l'induction spécifique d'un sous-ensemble de gènes qui varie en fonction du niveau de différentiation cellulaire. La modulation sélective de PPARγ par le MEHP provoque une adipogenèse modérée dans des modèles cellulaires alors que l'exposition de modèles animaux aux phtalates induit des effets bénéfiques sur la tolérance au glucose et sur le développement de l'obésité. Toutefois, les phtalates ont une action complexe sur le métabolisme glucido-lipidique en faisant intervenir PPARγ mais aussi probablement PPARα et PPARß. Cette démonstration moléculaire et physiologique de l'interférence des polluants avec les récepteurs nucléaires PPAR souligne un rôle important de l'exposition à de tels composés dans les régulations métaboliques.

How T cell receptors interact with peptide-MHCs: A multiple steered molecular dynamics study.

The T-cell receptor (TCR) interaction with antigenic peptides (p) presented by the major histocompatibility complex (MHC) molecule is a key determinant of immune response. In addition, TCR-pMHC interactions offer examples of features more generally pertaining to protein-protein recognition: subtle specificity and cross-reactivity. Despite their importance, molecular details determining the TCR-pMHC binding remain unsolved. However, molecular simulation provides the opportunity to investigate some of these aspects. In this study, we perform extensive equilibrium and steered molecular dynamics simulations to study the unbinding of three TCR-pMHC complexes. As a function of the dissociation reaction coordinate, we are able to obtain converged H-bond counts and energy decompositions at different levels of detail, ranging from the full proteins, to separate residues and water molecules, down to single atoms at the interface. Many observed features do not support a previously proposed two-step model for TCR recognition. Our results also provide keys to interpret experimental point-mutation results. We highlight the role of water both in terms of interface resolvation and of water molecules trapped in the bound complex. Importantly, we illustrate how two TCRs with similar reactivity and structures can have essentially different binding strategies. Proteins 2011; © 2011 Wiley-Liss, Inc.

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.

On the use of spatially discrete data to compute energy and mass balance

In many practical applications the state of field soils is monitored by recording the evolution of temperature and soil moisture at discrete depths. We theoretically investigate the systematic errors that arise when mass and energy balances are computed directly from these measurements. We show that, even with no measurement or model errors, large residuals might result when finite difference approximations are used to compute fluxes and storage term. To calculate the limits set by the use of spatially discrete measurements on the accuracy of balance closure, we derive an analytical solution to estimate the residual on the basis of the two key parameters: the penetration depth and the distance between the measurements. When the thickness of the control layer for which the balance is computed is comparable to the penetration depth of the forcing (which depends on the thermal diffusivity and on the forcing period) large residuals arise. The residual is also very sensitive to the distance between the measurements, which requires accurately controlling the position of the sensors in field experiments. We also demonstrate that, for the same experimental setup, mass residuals are sensitively larger than the energy residuals due to the nonlinearity of the moisture transport equation. Our analysis suggests that a careful assessment of the systematic mass error introduced by the use of spatially discrete data is required before using fluxes and residuals computed directly from field measurements.

Simulation of interface dynamics across pore discontinuities

We study the dynamics of a water-oil meniscus moving from a smaller to a larger pore. The process is characterised by an abrupt change in the configuration, yielding a sudden energy release. A theoretic study for static conditions provides analytical solutions of the surface energy content of the system. Although the configuration after the sudden energy release is energetically more convenient, an energy barrier must be overcome before the process can happen spontaneously. The energy barrier depends on the system geometry and on the flow parameters. The analytical results are compared to numerical simulations that solve the full Navier-Stokes equation in the pore space and employ the Volume Of Fluid (VOF) method to track the evolution of the interface. First, the numerical simulations of a quasi-static process are validated by comparison with the analytical solutions for a static meniscus, then numerical simulations with varying injection velocity are used to investigate dynamic effects on the configuration change. During the sudden energy jump the system exhibits an oscillatory behaviour. Extension to more complex geometries might elucidate the mechanisms leading to a dynamic capillary pressure and to bifurcations in final distributions of fluid phases in porous

Protein-ligand binding free energy estimation using molecular mechanics and continuum electrostatics. Application to HIV-1 protease inhibitors.

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. Conformational sampling of the protein-ligand complex is performed by molecular dynamics (MD) in vacuo and the solvent effect is calculated a posteriori by solving the Poisson or the Poisson-Boltzmann equation for selected frames of the trajectory. The binding free energy is written as a linear combination of the buried surface upon complexation, SASbur, the electrostatic interaction energy between the ligand and the protein, Eelec, and the difference of the solvation free energies of the complex and the isolated ligand and protein, deltaGsolv. The method uses the buried surface upon complexation to account for the non-polar contribution to the binding free energy because it is less sensitive to the details of the structure than the van der Waals interaction energy. The parameters of the method are developed for a training set of 16 HIV-1 protease-inhibitor complexes of known 3D structure. A correlation coefficient of 0.91 was obtained with an unsigned mean error of 0.8 kcal/mol. When applied to a set of 25 HIV-1 protease-inhibitor complexes of unknown 3D structures, the method provides a satisfactory correlation between the calculated binding free energy and the experimental pIC5o without reparametrization.

Roles of astrocytic sodium dynamics in the neurometabolic coupling : a study by cellular imaging

RESUME LARGE PUBLIC Le système nerveux central est principalement composé de deux types de cellules :les neurones et les cellules gliales. Ces dernières, bien que l'emportant en nombre sur les neurones, ont longtemps été considérées comme des cellules sans intérêts par les neuroscientifiques. Hors, les connaissances modernes à leurs sujets indiquent qu'elles participent à la plupart des tâches physiologiques du cerveau. Plus particulièrement, elles prennent part aux processus énergétiques cérébraux. Ceux-ci, en plus d'être vitaux, sont particulièrement intrigants puisque le cerveau représente seulement 2 % de la masse corporelle mais consomme environ 25 % du glucose (substrat énergétique) corporel. Les astrocytes, un type de cellules gliales, jouent un rôle primordial dans cette formidable utilisation de glucose par le cerveau. En effet, l'activité neuronale (transmission de l'influx nerveux) est accompagnée d'une augmentation de la capture de glucose, issu de la circulation sanguine, par les astrocytes. Ce phénomène est appelé le «couplage neurométabolique » entre neurones et astrocytes. L'ion sodium fait partie des mécanismes cellulaires entrant en fonction lors de ces processus. Ainsi, dans le cadre de cette thèse, les aspects dynamiques de la régulation du sodium astrocytaire et leurs implications dans le couplage neurométabolique ont été étudiés par des techniques d'imagerie cellulaires. Ces études ont démontré que les mitochondries, machineries cellulaires convertissant l'énergie contenue dans le glucose, participent à la régulation du sodium astrocytaire. De plus, ce travail de thèse a permis de découvrir que les astrocytes sont capables de se transmettre, sous forme de vagues de sodium se propageant de cellules en cellules, un message donnant l'ordre d'accroître leur consommation d'énergie. Cette voie de signalisation leur permettrait de fournir de l'énergie aux neurones suite à leur activation. RESUME Le glutamate libéré dans la fente synaptique pendant l'activité neuronale, est éliminé par les astrocytes environnants. Le glutamate est co-transporté avec des ions sodiques, induisant une augmentation intracellulaire de sodium (Na+i) dans les astrocytes. Cette élévation de Na+i déclenche une cascade de mécanismes moléculaires qui aboutissent à la production de substrats énergétiques pouvant être utilisés par les neurones. Durant cette thèse, la mesure simultanée du sodium mitochondrial (Na+mit) et cytosolique par des techniques d'imagerie utilisant des sondes fluorescentes spécifiques, a indiqué que les variations de Na+i induites par le transport du glutamate sont transmises aux mitochondries. De plus, les voies d'entrée et de sortie du sodium mitochondrial ont été identifiées. L'échangeur de Na+ et de Ca2+ mitochondrial semble jouer un rôle primordial dans l'influx de Na+mit, alors que l'efflux de Na+mit est pris en charge par l'échangeur de Na+ et de H+ mitochondrial. L'étude du Na+mit a nécessité l'utilisation d'un système de photoactivation. Les sources de lumière ultraviolette (UV) classiques utilisées à cet effet (lasers, lampes à flash) ayant plusieurs désavantages, une alternative efficace et peu coûteuse a été développée. Il s'agit d'un système compact utilisant une diode électroluminescente (LED) à haute puissance et de longueur d'onde de 365nm. En plus de leurs rôles dans le couplage neurométabolique, les astrocytes participent à la signalisation multicellulaire en transmettant des vagues intercellulaires de calcium. Ce travail de thèse démontre également que des vagues intercellulaires de sodium peuvent être évoquées en parallèle à ces vagues calciques. Le glutamate, suite à sa libération par un mécanisme dépendent du calcium, est réabsorbé par les transporteurs au glutamate. Ce mécanisme a pour conséquence la génération de vagues sodiques se propageant de cellules en cellules. De plus, ces vagues sodiques sont corrélées spatialement avec une consommation accrue de glucose par les astrocytes. En conclusion, ce travail de thèse a permis de montrer que le signal sodique astrocytaire, déclenché en réponse au glutamate, se propage à la fois de façon intracellulaire aux mitochondries et de façon intercellulaire. Ces résultats suggèrent que les astrocytes fonctionnent comme un réseau de cellules nécessaire au couplage énergétique concerté entre neurones et astrocytes et que le sodium est un élément clé dans les mécanismes de signalisations cellulaires sous-jacents. SUMMARY Glutamate, released in the synaptic cleft during neuronal activity, is removed by surrounding astrocytes. Glutamate is taken-up with Na+ ions by specific transporters, inducing an intracellular Na+ (Na+i) elevation in astrocytes which triggers a cascade of molecular mechanisms that provides metabolic substrates to neurons. Thus, astrocytic Na+i homeostasis represents a key component of the so-called neurometabolic coupling. In this context, the first part of this thesis work was aimed at investigating whether cytosolic Na+ changes are transmitted to mitochondria, which could therefore influence their function and contribute to the overall intracellular Na+ regulation. Simultaneous monitoring of both mitochondrial Na+ (Na+mit) and cytosolic Na+ changes with fluorescent dyes revealed that glutamate-evoked cytosolic Na+ elevations are indeed transmitted to mitochondria. The mitochondrial Na+/Ca2+ exchangers have a prominent role in the regulation of Na+mit influx pathway, and Na+mit extrusion appears to be mediated by Na+/H+ exchangers. To demonstrate the implication of Na+/Ca2+ exchangers, this study has required the technical development of an UV-flash photolysis system. Because light sources for flash photolysis have to be powerful and in the near UV range, the use of UV lasers or flash lamps is usually required. As an alternative to these UV sources that have several drawbaks, we developped a compact, efficient and lowcost flash photolysis system which employs a high power 365nm light emitting diode. In addition to their role in neurometabolic coupling, astrocytes participate in multicellular signaling by transmitting intercellular Ca2+ waves. The third part of this thesis show that intercellular Na+ waves can be evoked in parallel to Ca2+ waves. Glutamate released by a Ca2+ wave-dependent mechanism is taken up by glutamate transporters, resulting in a regenerative propagation of cytosolic Na+ increases. Na+ waves in turn lead to a spatially correlated increase in glucose uptake. In conclusion, the present thesis demonstrates that glutamate-induced Na+ changes occurring in the cytosol of astrocytes propagate to both the mitochondrial matrix and the astrocytic network. These results furthermore support the view that astrocytic Na+ is a signal coupled to the brain energy metabolism.