Existence of solutions and asymtptotic limits of the Euler-Poisson and the quantum drift-diffusion systems

My work concerns two different systems of equations used in the mathematical modeling of semiconductors and plasmas: the Euler-Poisson system and the quantum drift-diffusion system. The first is given by the Euler equations for the conservation of mass and momentum, with a Poisson equation for the electrostatic potential. The second one takes into account the physical effects due to the smallness of the devices (quantum effects). It is a simple extension of the classical drift-diffusion model which consists of two continuity equations for the charge densities, with a Poisson equation for the electrostatic potential. Using an asymptotic expansion method, we study (in the steady-state case for a potential flow) the limit to zero of the three physical parameters which arise in the Euler-Poisson system: the electron mass, the relaxation time and the Debye length. For each limit, we prove the existence and uniqueness of profiles to the asymptotic expansion and some error estimates. For a vanishing electron mass or a vanishing relaxation time, this method gives us a new approach in the convergence of the Euler-Poisson system to the incompressible Euler equations. For a vanishing Debye length (also called quasineutral limit), we obtain a new approach in the existence of solutions when boundary layers can appear (i.e. when no compatibility condition is assumed). Moreover, using an iterative method, and a finite volume scheme or a penalized mixed finite volume scheme, we numerically show the smallness condition on the electron mass needed in the existence of solutions to the system, condition which has already been shown in the literature. In the quantum drift-diffusion model for the transient bipolar case in one-space dimension, we show, by using a time discretization and energy estimates, the existence of solutions (for a general doping profile). We also prove rigorously the quasineutral limit (for a vanishing doping profile). Finally, using a new time discretization and an algorithmic construction of entropies, we prove some regularity properties for the solutions of the equation obtained in the quasineutral limit (for a vanishing pressure). This new regularity permits us to prove the positivity of solutions to this equation for at least times large enough.

Contribution of Non-Covalent Interactions and Electronic Effects on the Conformational Landscape and Tautomeric Equilibria of Molecules and Molecular Complexes: Structural and Dynamical Data from Rotational Spectroscopy

This thesis concerns the study of complex conformational surfaces and tautomeric equilibria of molecules and molecular complexes by quantum chemical methods and rotational spectroscopy techniques. In particular, the focus of this research is on the effects of substitution and noncovalent interactions in determining the energies and geometries of different conformers, tautomers or molecular complexes. The Free-Jet Absorption Millimeter Wave spectroscopy and the Pulsed-Jet Fourier Transform Microwave spectroscopy have been applied to perform these studies and the obtained results showcase the suitability of these techniques for the study of conformational surfaces and intermolecular interactions. The series of investigations of selected medium-size molecules and complexes have shown how different instrumental setups can be used to obtain a variety of results on molecular properties. The systems studied, include molecules of biological interest such as anethole and molecules of astrophysical interest such as N-methylaminoethanol. Moreover halogenation effects have been investigated on halogen substituted tautomeric systems (5-chlorohydroxypyridine and 6-chlorohydroxypyridine), where it has shown that the position of the inserted halogen atom affects the prototropic equilibrium. As for fluorination effects, interesting results have been achieved investigating some small complexes where a molecule of water is used as a probe to reveal the changes on the electrostatic potential of different fluorinated compounds: 2-fluoropyridine, 3-fluoropyridine and penta-fluoropyridine. While in the case of the molecular complex between water and 2-fluoropyridine and 3-fluoropyridine the geometry of the complex with one water molecule is analogous to that of pyridine with the water molecule linked to the pyridine nitrogen, the case of pentafluoropyridine reveals the effect of perfluorination and the water oxygen points towards the positive center of the pyridine ring. Additional molecular adducts with a molecule of water have been analyzed (benzylamine-water and acrylic acid-water) in order to reveal the stabilizing driving forces that characterize these complexes.

Modeling and simulation of fibrinogen and its adsorption behavior

In this thesis different approaches for the modeling and simulation of the blood protein fibrinogen are presented. The approaches are meant to systematically connect the multiple time and length scales involved in the dynamics of fibrinogen in solution and at inorganic surfaces. The first part of the thesis will cover simulations of fibrinogen on an all atom level. Simulations of the fibrinogen protomer and dimer are performed in explicit solvent to characterize the dynamics of fibrinogen in solution. These simulations reveal an unexpectedly large and fast bending motion that is facilitated by molecular hinges located in the coiled-coil region of fibrinogen. This behavior is characterized by a bending and a dihedral angle and the distribution of these angles is measured. As a consequence of the atomistic detail of the simulations it is possible to illuminate small scale behavior in the binding pockets of fibrinogen that hints at a previously unknown allosteric effect. In a second step atomistic simulations of the fibrinogen protomer are performed at graphite and mica surfaces to investigate initial adsorption stages. These simulations highlight the different adsorption mechanisms at the hydrophobic graphite surface and the charged, hydrophilic mica surface. It is found that the initial adsorption happens in a preferred orientation on mica. Many effects of practical interest involve aggregates of many fibrinogen molecules. To investigate such systems, time and length scales need to be simulated that are not attainable in atomistic simulations. It is therefore necessary to develop lower resolution models of fibrinogen. This is done in the second part of the thesis. First a systematically coarse grained model is derived and parametrized based on the atomistic simulations of the first part. In this model the fibrinogen molecule is represented by 45 beads instead of nearly 31,000 atoms. The intra-molecular interactions of the beads are modeled as a heterogeneous elastic network while inter-molecular interactions are assumed to be a combination of electrostatic and van der Waals interaction. A method is presented that determines the charges assigned to beads by matching the electrostatic potential in the atomistic simulation. Lastly a phenomenological model is developed that represents fibrinogen by five beads connected by rigid rods with two hinges. This model only captures the large scale dynamics in the atomistic simulations but can shed light on experimental observations of fibrinogen conformations at inorganic surfaces.

Correlation between Accurate Electron Density and Linear Optical Properties in Amino Acid Derivatives: L-Histidinium Hydrogen Oxalate

The accurate electron density and linear optical properties of L-histidinium hydrogen oxalate are discussed. Two high-resolution single crystal X-ray diffraction experiments were performed and compared with density functional calculations in the solid state as well as in the gas phase. The crystal packing and the hydrogen bond network are accurately investigated using topological analysis based on quantum theory of atoms in molecules, Hirshfeld surface analysis, and electrostatic potential mapping. The refractive indices are computed from couple perturbed Kohn-Sham calculations and measured experimentally. Moreover, distributed atomic polarizabilities are used to analyze the origin of the linear susceptibility in the crystal, in order to separate molecular and intermolecular causes. The optical properties are also correlated with the electron density distribution. This compound also offers the possibility to test the electron density building block approach for material science and different refinement schemes for accurate positions and displacement parameters of hydrogen atoms, in the absence of neutron diffraction data.

Calculation of crystal optical properties from molecular electron density

We have recently developed a method to obtain distributed atomic polarizabilities adopting a partitioning of the molecular electron density (for example, the Quantum Theory of Atoms in Molecules, [1]), calculated with or without an applied electric field. The procedure [2] allows to obtained atomic polarizability tensors, which are perfectly exportable, because quite representative of an atom in a given functional group. Among the many applications of this idea, the calculation of crystal susceptibility is easily available, either from a rough estimation (the polarizability of the isolated molecule is used) or from a more precise estimation (the polarizability of a molecule embedded in a cluster representing the first coordination sphere is used). Lorentz factor is applied to include the long range effect of packing, which is enhancing the molecular polarizability. Simple properties like linear refractive index or the gyration tensor can be calculated at relatively low costs and with good precision. This approach is particularly useful within the field of crystal engineering of organic/organometallic materials, because it would allow a relatively easy prediction of a property as a function of the packing, thus allowing "reverse crystal engineering". Examples of some amino acid crystals and salts of amino acids [3] will be illustrated, together with other crystallographic or non-crystallographic applications. For example, the induction and dispersion energies of intermolecular interactions could be calculated with superior precision (allowing anisotropic van der Waals interactions). This could allow revision of some commonly misunderstood intermolecular interactions, like the halogen bonding (see for example the recent remarks by Stone or Gilli [4]). Moreover, the chemical reactivity of coordination complexes could be reinvestigated, by coupling the conventional analysis of the electrostatic potential (useful only in the circumstances of hard nucleophilic/electrophilic interaction) with the distributed atomic polarizability. The enhanced reactivity of coordinated organic ligands would be better appreciated. [1] R. F. W. Bader, Atoms in Molecules: A Quantum Theory. Oxford Univ. Press, 1990. [2] A. Krawczuk-Pantula, D. Pérez, K. Stadnicka, P. Macchi, Trans. Amer. Cryst. Ass. 2011, 1-25 [3] A. S. Chimpri1, M. Gryl, L. H.R. Dos Santos1, A. Krawczuk, P. Macchi Crystal Growth & Design, in the press. [4] a) A. J. Stone, J. Am. Chem. Soc. 2013, 135, 7005−7009; b) V. Bertolasi, P. Gilli, G. Gilli Crystal Growth & Design, 2013, 12, 4758-4770.

Mechanisms of apoptosis in chronic lymphocytic leukemia (CLL)

Chronic lymphocytic leukemia (CLL) is an incurable disease characterized by the accumulation of terminally differentiated, mature B cells that do not progress beyond the G1 stage of cell cycle, suggesting that these cells possess intrinsic defects in apoptosis. Treatment relies heavily on chemotherapy (primarily nucleoside analogs and glucocorticoids) that may initially be effective in patients, but ultimately give rise to refractory, untreatable disease. The purpose of this study was to determine whether key components of the apoptotic machinery were intact in CLL lymphocytes, especially in patients refractory to therapy. ^ Activation of proteases has been shown to be at the core of the apoptotic pathway and this work demonstrates that protease activation is required for glucocorticoid and nucleoside analog-induced apoptosis in CLL cells. Inhibitors of serine proteases as well as caspase inhibitors blocked induced DNA fragmentation, and a peptide inhibitor of the nuclear scaffold (NS) protease completely suppressed both induced and spontaneous apoptosis. However, the NS protease inhibitor actually promoted several pro-apoptotic events, such as caspase activation, exposure of surface phosphatidylserine, and loss of mitochondrial membrane potential. These results suggested that the NS protease may interact with the apoptotic program in CLL cells at two separate points. ^ In order to further investigate the role of the NS protease in CLL, patient isolates were treated with proteasome inhibitors because of previous results suggesting that the ISIS protease might be a β subunit of the proteasome. Proteasome inhibitors induced massive DNA fragmentation in every patient tested, even in those resistant to the effects of glucocorticoid and nucleoside analogs in vitro. Several other features of apoptosis were also promoted by the proteasome inhibitor, including mitochondrial alterations such as release of cytochrome c and drops in mitochondrial membrane potential. Proteasome inhibitor-induced apoptosis was associated with inhibition of NFκB, a proteasome-regulated transcription factor that has been implicated in the suppression of apoptosis in a number of systems. The NS protease inhibitor also caused a decrease in active NFκB, suggesting that the proapoptotic effects of this agent might be due to depletion of NFκB. ^ Given these findings, the role of NFκB, in conferring survival in CLL was investigated. Glucocorticoid hormone treatment was shown to cause decreases in the activity of the transcription factor, while phorbol dibutyrate, which blocks glucocorticoid-induced DNA fragmentation, was capable of upregulating NFκB. Compellingly, introduction of an undegradable form of the constitutive NFκB inhibitor, IκB, caused DNA fragmentation in several patient isolates, some of which were resistant to glucocorticoid in vitro. Transcription of anti-apoptotic proteins by NFκB was postulated to be responsible for its effects on survival, but Bcl-2 levels did not fluctuate with glucocorticoid or proteasome inhibitor treatment. ^ The in vitro values generated from these studies were organized into a database containing numbers for over 250 patients. Correlation of relevant clinical parameters revealed that levels of spontaneous apoptosis in vitro differ significantly between Rai stages. Importantly, in vitro resistance to nucleoside analogs or glucocorticoids predicted resistance to chemotherapy in vivo, and inability to achieve remission. ^

Cosmogenic activation of xenon and copper

Rare event search experiments using liquid xenon as target and detection medium require ultra-low background levels to fully exploit their physics potential. Cosmic ray induced activation of the detector components and, even more importantly, of the xenon itself during production, transportation and storage at the Earth's surface, might result in the production of radioactive isotopes with long half-lives, with a possible impact on the expected background. We present the first dedicated study on the cosmogenic activation of xenon after 345 days of exposure to cosmic rays at the Jungfraujoch research station at 3470m above sea level, complemented by a study of copper which has been activated simultaneously. We have directly observed the production of 7Be, 101Rh, 125Sb, 126I and 127Xe in xenon, out of which only 125Sb could potentially lead to background for a multi-ton scale dark matter search. The production rates for five out of eight studied radioactive isotopes in copper are in agreement with the only existing dedicated activation measurement, while we observe lower rates for the remaining ones. The specific saturation activities for both samples are also compared to predictions obtained with commonly used software packages, where we observe some underpredictions, especially for xenon activation.

Low Work-Function Thermionic Emission and Orbital-Motion-Limited Ion Collection at Bare-Tether Cathodic Contact

Una amarra electrodinámica (electrodynamic tether) opera sobre principios electromagnéticos intercambiando momento con la magnetosfera planetaria e interactuando con su ionosfera. Es un subsistema pasivo fiable para desorbitar etapas de cohetes agotadas y satélites al final de su misión, mitigando el crecimiento de la basura espacial. Una amarra sin aislamiento captura electrones del plasma ambiente a lo largo de su segmento polarizado positivamente, el cual puede alcanzar varios kilómetros de longitud, mientras que emite electrones de vuelta al plasma mediante un contactor de plasma activo de baja impedancia en su extremo catódico, tal como un cátodo hueco (hollow cathode). En ausencia de un contactor catódico activo, la corriente que circula por una amarra desnuda en órbita es nula en ambos extremos de la amarra y se dice que ésta está flotando eléctricamente. Para emisión termoiónica despreciable y captura de corriente en condiciones limitadas por movimiento orbital (orbital-motion-limited, OML), el cociente entre las longitudes de los segmentos anódico y catódico es muy pequeño debido a la disparidad de masas entre iones y electrones. Tal modo de operación resulta en una corriente media y fuerza de Lorentz bajas en la amarra, la cual es poco eficiente como dispositivo para desorbitar. El electride C12A7 : e−, que podría presentar una función de trabajo (work function) tan baja como W = 0.6 eV y un comportamiento estable a temperaturas relativamente altas, ha sido propuesto como recubrimiento para amarras desnudas. La emisión termoiónica a lo largo de un segmento así recubierto y bajo el calentamiento de la operación espacial, puede ser más eficiente que la captura iónica. En el modo más simple de fuerza de frenado, podría eliminar la necesidad de un contactor catódico activo y su correspondientes requisitos de alimentación de gas y subsistema de potencia, lo que resultaría en un sistema real de amarra “sin combustible”. Con este recubrimiento de bajo W, cada segmento elemental del segmento catódico de una amarra desnuda de kilómetros de longitud emitiría corriente como si fuese parte de una sonda cilíndrica, caliente y uniformemente polarizada al potencial local de la amarra. La operación es similar a la de una sonda de Langmuir 2D tanto en los segmentos catódico como anódico. Sin embargo, en presencia de emisión, los electrones emitidos resultan en carga espacial (space charge) negativa, la cual reduce el campo eléctrico que los acelera hacia fuera, o incluso puede desacelerarlos y hacerlos volver a la sonda. Se forma una doble vainas (double sheath) estable con electrones emitidos desde la sonda e iones provenientes del plasma ambiente. La densidad de corriente termoiónica, variando a lo largo del segmento catódico, podría seguir dos leyes distintas bajo diferentes condiciones: (i) la ley de corriente limitada por la carga espacial (space-charge-limited, SCL) o (ii) la ley de Richardson-Dushman (RDS). Se presenta un estudio preliminar sobre la corriente SCL frente a una sonda emisora usando la teoría de vainas (sheath) formada por la captura iónica en condiciones OML, y la corriente electrónica SCL entre los electrodos cilíndricos según Langmuir. El modelo, que incluye efectos óhmicos y el efecto de transición de emisión SCL a emisión RDS, proporciona los perfiles de corriente y potencial a lo largo de la longitud completa de la amarra. El análisis muestra que en el modo más simple de fuerza de frenado, bajo condiciones orbitales y de amarras típicas, la emisión termoiónica proporciona un contacto catódico eficiente y resulta en una sección catódica pequeña. En el análisis anterior, tanto la transición de emisión SCL a RD como la propia ley de emisión SCL consiste en un modelo muy simplificado. Por ello, a continuación se ha estudiado con detalle la solución de vaina estacionaria de una sonda con emisión termoiónica polarizada negativamente respecto a un plasma isotrópico, no colisional y sin campo magnético. La existencia de posibles partículas atrapadas ha sido ignorada y el estudio incluye tanto un estudio semi-analítico mediante técnica asintóticas como soluciones numéricas completas del problema. Bajo las tres condiciones (i) alto potencial, (ii) R = Rmax para la validez de la captura iónica OML, y (iii) potencial monotónico, se desarrolla un análisis asintótico auto-consistente para la estructura de plasma compleja que contiene las tres especies de cargas (electrones e iones del plasma, electrones emitidos), y cuatro regiones espaciales distintas, utilizando teorías de movimiento orbital y modelos cinéticos de las especies. Aunque los electrones emitidos presentan carga espacial despreciable muy lejos de la sonda, su efecto no se puede despreciar en el análisis global de la estructura de la vaina y de dos capas finas entre la vaina y la región cuasi-neutra. El análisis proporciona las condiciones paramétricas para que la corriente sea SCL. También muestra que la emisión termoiónica aumenta el radio máximo de la sonda para operar dentro del régimen OML y que la emisión de electrones es mucho más eficiente que la captura iónica para el segmento catódico de la amarra. En el código numérico, los movimientos orbitales de las tres especies son modelados para potenciales tanto monotónico como no-monotónico, y sonda de radio R arbitrario (dentro o más allá del régimen de OML para la captura iónica). Aprovechando la existencia de dos invariante, el sistema de ecuaciones Poisson-Vlasov se escribe como una ecuación integro-diferencial, la cual se discretiza mediante un método de diferencias finitas. El sistema de ecuaciones algebraicas no lineal resultante se ha resuelto de con un método Newton-Raphson paralelizado. Los resultados, comparados satisfactoriamente con el análisis analítico, proporcionan la emisión de corriente y la estructura del plasma y del potencial electrostático. ABSTRACT An electrodynamic tether operates on electromagnetic principles and exchanges momentum through the planetary magnetosphere, by continuously interacting with the ionosphere. It is a reliable passive subsystem to deorbit spent rocket stages and satellites at its end of mission, mitigating the growth of orbital debris. A tether left bare of insulation collects electrons by its own uninsulated and positively biased segment with kilometer range, while electrons are emitted by a low-impedance active device at the cathodic end, such as a hollow cathode, to emit the full electron current. In the absence of an active cathodic device, the current flowing along an orbiting bare tether vanishes at both ends and the tether is said to be electrically floating. For negligible thermionic emission and orbital-motion-limited (OML) collection throughout the entire tether (electron/ion collection at anodic/cathodic segment, respectively), the anodic-to-cathodic length ratio is very small due to ions being much heavier, which results in low average current and Lorentz drag. The electride C12A7 : e−, which might present a possible work function as low as W = 0.6 eV and moderately high temperature stability, has been proposed as coating for floating bare tethers. Thermionic emission along a thus coated cathodic segment, under heating in space operation, can be more efficient than ion collection and, in the simplest drag mode, may eliminate the need for an active cathodic device and its corresponding gas-feed requirements and power subsystem, which would result in a truly “propellant-less” tether system. With this low-W coating, each elemental segment on the cathodic segment of a kilometers-long floating bare-tether would emit current as if it were part of a hot cylindrical probe uniformly polarized at the local tether bias, under 2D probe conditions that are also applied to the anodic-segment analysis. In the presence of emission, emitted electrons result in negative space charge, which decreases the electric field that accelerates them outwards, or even reverses it, decelerating electrons near the emitting probe. A double sheath would be established with electrons being emitted from the probe and ions coming from the ambient plasma. The thermionic current density, varying along the cathodic segment, might follow two distinct laws under different con ditions: i) space-charge-limited (SCL) emission or ii) full Richardson-Dushman (RDS) emission. A preliminary study on the SCL current in front of an emissive probe is presented using the orbital-motion-limited (OML) ion-collection sheath and Langmuir’s SCL electron current between cylindrical electrodes. A detailed calculation of current and bias profiles along the entire tether length is carried out with ohmic effects considered and the transition from SCL to full RDS emission is included. Analysis shows that in the simplest drag mode, under typical orbital and tether conditions, thermionic emission provides efficient cathodic contact and leads to a short cathodic section. In the previous analysis, both the transition between SCL and RDS emission and the current law for SCL condition have used a very simple model. To continue, considering an isotropic, unmagnetized, colissionless plasma and a stationary sheath, the probe-plasma contact is studied in detail for a negatively biased probe with thermionic emission. The possible trapped particles are ignored and this study includes both semianalytical solutions using asymptotic analysis and complete numerical solutions. Under conditions of i) high bias, ii) R = Rmax for ion OML collection validity, and iii) monotonic potential, a self-consistent asymptotic analysis is carried out for the complex plasma structure involving all three charge species (plasma electrons and ions, and emitted electrons) and four distinct spatial regions using orbital motion theories and kinetic modeling of the species. Although emitted electrons present negligible space charge far away from the probe, their effect cannot be neglected in the global analysis for the sheath structure and two thin layers in between the sheath and the quasineutral region. The parametric conditions for the current to be space-chargelimited are obtained. It is found that thermionic emission increases the range of probe radius for OML validity and is greatly more effective than ion collection for cathodic contact of tethers. In the numerical code, the orbital motions of all three species are modeled for both monotonic and non-monotonic potential, and for any probe radius R (within or beyond OML regime for ion collection). Taking advantage of two constants of motion (energy and angular momentum), the Poisson-Vlasov equation is described by an integro differential equation, which is discretized using finite difference method. The non-linear algebraic equations are solved using a parallel implementation of the Newton-Raphson method. The results, which show good agreement with the analytical results, provide the results for thermionic current, the sheath structure, and the electrostatic potential.

Repeated, but not acute, stress suppresses inflammatory plasma extravasation

Clinical findings suggest that inflammatory disease symptoms are aggravated by ongoing, repeated stress, but not by acute stress. We hypothesized that, compared with single acute stressors, chronic repeated stress may engage different physiological mechanisms that exert qualitatively different effects on the inflammatory response. Because inhibition of plasma extravasation, a critical component of the inflammatory response, has been associated with increased disease severity in experimental arthritis, we tested for a potential repeated stress-induced inhibition of plasma extravasation. Repeated, but not single, exposures to restraint stress produced a profound inhibition of bradykinin-induced synovial plasma extravasation in the rat. Experiments examining the mechanism of inhibition showed that the effect of repeated stress was blocked by adrenalectomy, but not by adrenal medullae denervation, suggesting that the adrenal cortex mediates this effect. Consistent with known effects of stress and with mediation by the adrenal cortex, restraint stress evoked repeated transient elevations of plasma corticosterone levels. This elevated corticosterone was necessary and sufficient to produce inhibition of plasma extravasation because the stress-induced inhibition was blocked by preventing corticosterone synthesis and, conversely, induction of repeated transient elevations in plasma corticosterone levels mimicked the effects of repeated stress. These data suggest that repetition of a mild stressor can induce changes in the physiological state of the animal that enable a previously innocuous stressor to inhibit the inflammatory response. These findings provide a potential explanation for the clinical association between repeated stress and aggravation of inflammatory disease symptoms and provide a model for study of the biological mechanisms underlying the stress-induced aggravation of chronic inflammatory diseases.

Theoretical investigation of the [1,2]-sigmatropic hydrogen migration in the mechanism of oxidation of 2-aminobenzoyl-CoA by 2-aminobenzoyl-CoA monooxygenase/reductase

The flavin hydroperoxide at the active site of the mixed-function oxidase 2-aminobenzoyl-CoA monooxygenase/reductase (Azoarcus evansii) transfers an oxygen to the 5-position of the 2-aminobenzoyl-CoA substrate to provide the alkoxide intermediate II−. Hydrogen migration from C5 to C6 follows this monooxygenation. The nature of the monooxygenation intermediate and plausible competing reactions leading to hydrogen migration have been considered. Ab initio molecular orbital theory has been used to calculate structures and electron distributions in intermediate and transition state structures. Electrostatic potential surface calculations establish that the transition state and product, associated with the C5 to C6 hydrogen transfer, are stabilized by electron distribution to the benzoyl-CoA thioester carbonyl oxygen. This is not so for the transition state and product associated with hydrogen transfer from C5 to C4. The activation energy for the 5,6-shift is 2.5 kcal/mol lower than that for the 5,4-shift. In addition, the product of the hydrogen 5,6-shift is more stable than is the product of the hydrogen 5,4-shift, by ≈6 kcal/mol. These results explain why only the shift of hydrogen from C5 to C6 is observed experimentally. Oxygen transfer and hydrogen migration almost coincide in the gas phase (activation energy of ≈0.6 kcal/mol, equivalent to a single bond vibration). Enzymatic formation of alkoxide II− requires its stabilization; thus, the rate constant for its breakdown would be slower than in the gas phase.

Solution structure of the antiapoptotic protein bcl-2

The structures of two isoforms of Bcl-2 that differ by two amino acids have been determined by NMR spectroscopy. Because wild-type Bcl-2 behaved poorly in solution, the structures were determined by using Bcl-2/Bcl-xL chimeras in which part of the putative unstructured loop of Bcl-2 was replaced with a shortened loop from Bcl-xL. These chimeric proteins have a low pI compared with the wild-type protein and are soluble. The structures of the two Bcl-2 isoforms consist of 6 α-helices with a hydrophobic groove on the surface similar to that observed for the homologous protein, Bcl-xL. Comparison of the Bcl-2 structures to that of Bcl-xL shows that although the overall fold is the same, there are differences in the structural topology and electrostatic potential of the binding groove. Although the structures of the two isoforms of Bcl-2 are virtually identical, differences were observed in the ability of the proteins to bind to a 25-residue peptide from the proapoptotic Bad protein and a 16-residue peptide from the proapoptotic Bak protein. These results suggest that there are subtle differences in the hydrophobic binding groove in Bcl-2 that may translate into differences in antiapoptotic activity for the two isoforms.

Crystal structure at 2.6-A resolution of human macrophage migration inhibitory factor.

Macrophage migration inhibitory factor (MIF) was the first cytokine to be described, but for 30 years its role in the immune response remained enigmatic. In recent studies, MIF has been found to be a novel pituitary hormone and the first protein identified to be released from immune cells on glucocorticoid stimulation. Once secreted, MIF counterregulates the immunosuppressive effects of steroids and thus acts as a critical component of the immune system to control both local and systemic immune responses. We report herein the x-ray crystal structure of human MIF to 2.6 angstrom resolution. The protein is a trimer of identical subunits. Each monomer contains two antiparallel alpha-helices that pack against a four-stranded beta-sheet. The monomer has an additional two beta-strands that interact with the beta-sheets of adjacent subunits to form the interface between monomers. The three beta-sheets are arranged to form a barrel containing a solvent-accessible channel that runs through the center of the protein along a molecular 3-fold axis. Electrostatic potential maps reveal that the channel has a positive potential, suggesting that it binds negatively charged molecules. The elucidated structure for MIF is unique among cytokines or hormonal mediators, and suggests that this counterregulator of glucocorticoid action participates in novel ligand-receptor interactions.

Transdifferentiation of myoblasts by the adipogenic transcription factors PPAR gamma and C/EBP alpha.

Skeletal muscle and adipose tissue development often has a reciprocal relationship in vivo, particularly in myodystrophic states. We have investigated whether determined myoblasts with no inherent adipogenic potential can be induced to transdifferentiate into mature adipocytes by the ectopic expression of two adipogenic transcription factors, PPAR gamma and C/EBP alpha. When cultured under optimal conditions for muscle differentiation, murine G8 myoblasts expressing PPAR gamma and C/EBP alpha show markedly reduced levels of the myogenic basic helix-loop-helix proteins MyoD, myogenin, MRF4, and myf5 and are completely unable to differentiate into myotubes. Under conditions permissive for adipogenesis including a PPAR activator, these cells differentiate into mature adipocytes that express molecular markers characteristic of this lineage. Our results demonstrate that a developmental switch between these two related but highly specialized cell types can be controlled by the expression of key adipogenic transcription factors. These factors have an ability to inhibit myogenesis that is temporally and functionally separate from their ability to stimulate adipogenesis.

Biochemical characterization and DNA binding properties of the MocR family member GabR, a PLP-dependent transcription factor

GabR è un fattore di trascrizione chimerico appartenente alla famiglia dei MocR/GabR, costituito da un dominio N-terminale elica-giro-elica di legame al DNA e un dominio effettore e/o di oligomerizzazione al C-terminale. I due domini sono connessi da un linker flessibile di 29 aminoacidi. Il dominio C-terminale è strutturalmente omologo agli enzimi aminotransferasici fold-type I, i quali, utilizzando il piridossal-5’-fosfato (PLP) come cofattore, sono direttamente coinvolti nel metabolismo degli aminoacidi. L’interazione contemporanea di PLP e acido γ-aminobutirrico (GABA) a GabR fa sì che questa promuova la trascrizione di due geni, gabT e gabD, implicati nel metabolismo del GABA. GabR cristallizza come un omodimero con una configurazione testa-coda. Il legame con la regione promotrice gabTD avviene attraverso il riconoscimento specifico di due sequenze dirette e ripetute (ATACCA), separate da uno spacer di 34 bp. In questo studio sono state indagate le proprietà biochimiche, strutturali e di legame al DNA della proteina GabR di Bacillus subtilis. L’analisi spettroscopica dimostra che GabR interagisce con il PLP formando l’aldimina interna, mentre in presenza di GABA si ottiene l’aldimina esterna. L’interazione fra il promotore gabTD e le forme holo e apo di GabR è stata monitorata mediante Microscopia a Forza atomica (AFM). In queste due condizioni di legame è stata stimata una Kd di circa 40 ηM. La presenza di GABA invece, determinava un incremento di circa due volte della Kd, variazioni strutturali nei complessi GabR-DNA e una riduzione del compattamento del DNA alla proteina, indipendentemente dalla sequenza del promotore in esame. Al fine di valutare il ruolo delle caratteristiche topologiche del promotore, sono state inserite cinque e dieci bp all’interno della regione spacer che separa le due sequenze ripetute dirette riconosciute da GabR. I significativi cambiamenti topologici riscontrati nel frammento aggiunto di cinque bp si riflettono anche sulla forte riduzione dell’affinità di legame verso la proteina. Al contrario, l’inserzione di 10 bp provoca solamente l’allontanamento delle sequenze ripetute dirette. L’assenza quindi di cambiamenti significativi nella topologia di questo promotore fa sì che l’affinità di legame per GabR rimanga pressoché inalterata rispetto al promotore non mutato. L’analisi del potenziale elettrostatico superficiale di GabR mostra la presenza di una fascia carica positivamente che si estende lungo un’intera faccia della proteina. Per verificare l’importanza di questa caratteristica di GabR nel meccanismo di interazione al DNA, sono stati preparati ed indagati i mutanti R129Q e K362-366Q, in cui la carica positiva superficiale risultava indebolita. L’affinità di legame dei mutanti di GabR per il DNA era inferiore rispetto alla proteina non mutata, in particolar modo nel mutante K362-366Q. Le evidenze acquisite suggeriscono che la curvatura intrinseca del promotore ed il corretto orientamento delle sequenze sulla doppia elica, più della distanza che le separa, siano critici per sostenere l’interazione con GabR. Oltre a questo, la superficie positiva di GabR è richiesta per accomodare la curvatura del DNA sul corpo della proteina. Alla luce di questo, l’interazione GabR-gabTD è un esempio di come il riconoscimento specifico di sequenze, la topologia del DNA e le caratteristiche strutturali della proteina siano contemporaneamente necessarie per sostenere un’interazione proteina-DNA stabile.

Why H Atom Prefers the On-Top Site and Alkali Metals Favor the Middle Hollow Site on the Basal Plane of Graphite

In this work, the different adsorption properties of H and alkali metal atoms on the basal plane of graphite are studied and compared using a density functional method on the same model chemistry level. The results show that H prefers the on-top site while alkali metals favor the middle hollow site of graphite basal plane due to the unique electronic structures of H, alkali metals, and graphite. H has a higher electronegativity than carbon, preferring to form a covalent bond with C atoms, whereas alkaline metals have lower electronegativity, tending to adsorb on the highest electrostatic potential sites. During adsorption, there are more charges transferred from alkali metal to graphite than from H to graphite.