Search NASA⌕ Search

SEARCH · Search NASA

Results for “Large deviations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Canonical Structure and Orthogonality of Forces and Currents in Irreversible Markov Chains

We discuss a canonical structure that provides a unifying description of dynamical large deviations for irreversible finite state Markov chains (continuous time), Onsager theory, and Macroscopic Fluctuation Theory (MFT). For Markov chains, this theory involves a non-linear relation between probability currents and their conjugate forces. Within this framework, we show how the forces can be split into two components, which are orthogonal to each other, in a generalised sense. This splitting allows a decomposition of the pathwise rate function into three terms, which have physical interpretations in terms of dissipation and convergence to equilibrium. Similar decompositions hold for rate functions at level 2 and level 2.5. These results clarify how bounds on entropy production and fluctuation theorems emerge from the underlying dynamical rules. We discuss how these results for Markov chains are related to similar structures within MFT, which describes hydrodynamic limits of such microscopic models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Front Propagation for Reaction–Diffusion Equations in Composite Structures

We consider asymptotic problems concerning the motion of interface separating the regions of large and small values of the solution of a reaction–diffusion equation in the media consisting of domains with different characteristics (composites). Under certain conditions, the motion can be described by the Huygens principle in the appropriate Finsler (e.g., Riemannian) metric. In general, the motion of the interface has, in a sense, non-local nature. In particular, the interface may move by jumps. We are mostly concerned with the nonlinear term that is of KPP type. The results are based on limit theorems for large deviations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Constitutive Modeling of the High-Temperature Flow Behavior of α-Ti Alloy Tube

In the hot metal gas forming process, the deformation conditions, such as temperature, strain rate and deformation degree, are often prominently changed. The understanding of the flow behavior of α-Ti seamless tubes over a relatively wide range of temperatures and strain rates is important. In this study, the stress–strain curves in the temperature range of 973-1123 K and the initial strain rate range of 0.0004-0.4 s{sup −1} were measured by isothermal tensile tests to conduct a constitutive analysis and a deformation behavior analysis. The results show that the flow stress decreases with the decrease in the strain rate and the increase of the deformation temperature. The Fields–Backofen model and Fields–Backofen–Zhang model were used to describe the stress–strain curves. The Fields–Backofen–Zhang model shows better predictability on the flow stress than the Fields–Backofen model, but there exists a large deviation in the deformation condition of 0.4 s{sup −1}. A modified Fields–Backofen–Zhang model is proposed, in which a strain rate term is introduced. This modified Fields–Backofen–Zhang model gives a more accurate description of the flow stress variation under hot forming conditions with a higher strain rate up to 0.4 s{sup −1}. Accordingly, it is reasonable to adopt the modified Fields–Backofen–Zhang model for the hot forming process which is likely to reach a higher strain rate, such as 0.4 s{sup −1}.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Triple junction drag effects during topological changes in the evolution of polycrystalline microstructures

Experiments, theory and atomistic simulations show that finite triple junction mobility results in non-equilibrium triple junction angles in evolving polycrystalline systems. These angles have been predicted and verified for cases where grain boundary migration is steady-state. Yet, steady-state never occurs during the evolution of polycrystalline microstructures as a result of changing grain size and topological events (e.g., grain face/edge switching - “T 1 ” process, or grain disappearance “T 2 ” or “T 3 ” processes). We examine the non-steady evolution of the triple junction angle in the vicinity of topological events and show that large deviations from equilibrium and/or steady-state angles occur. We analyze $\tau$ the characteristic relaxation time of triple junction angles τ by consideration of a pair of topological events, beginning from steady-state migration. Using numerical results and theoretical analysis we predict how the triple junction angle varies with time and how τ varies with triple junction mobility. We argue that it is precisely those cases where grain boundaries are moving quickly (e.g., topological process in nanocrystalline materials), that the classical steady-state prediction of the triple junction angle about finite triple junction mobility is inapplicable and may only be applied qualitatively.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Exponential number of equilibria and depinning threshold for a directed polymer in a random potential

By extending the Kac–Rice approach to manifolds of finite internal dimension, we show that the mean number ⟨N{sub tot}⟩ of all possible equilibria (i.e. force-free configurations, a.k.a. equilibrium points) of an elastic line (directed polymer), confined in a harmonic well and submitted to a quenched random Gaussian potential in dimension d=1+1, grows exponentially ⟨N{sub tot}⟩ ~exp(r L) with its length L. The growth rate r is found to be directly related to the generalized Lyapunov exponent (GLE) which is a moment-generating function characterizing the large-deviation type fluctuations of the solution to the initial value problem associated with the random Schrödinger operator of the 1D Anderson localization problem. For strong confinement, the rate r is small and given by a non-perturbative (instanton, Lifshitz tail-like) contribution to GLE. For weak confinement, the rate r is found to be proportional to the inverse Larkin length of the pinning theory. As an application, identifying the depinning with a landscape “topology trivialization” phenomenon, we obtain an upper bound for the depinning threshold f{sub c}, in the presence of an applied force, for elastic lines and d-dimensional manifolds, expressed through the mean modulus of the spectral determinant of the Laplace operators with a random potential. We also discuss the question of counting of stable equilibria. Finally, we extend the method to calculate the asymptotic number of equilibria at fixed energy (elastic, potential and total), and obtain the (annealed) distribution of the energy density over these equilibria (i.e. force-free configurations). Some connections with the Larkin model are also established.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

On the Mechanism of Bilayer Separation by Extrusion, or Why Your LUVs Are Not Really Unilamellar

Extrusion through porous filters is a widely used method for preparing biomimetic model membranes. Of primary importance in this approach is the efficient production of single bilayer (unilamellar) vesicles that eliminate the influence of interlamellar interactions and strictly define the bilayer surface area available to external reagents such as proteins. Sub-microscopic vesicles produced using extrusion are widely assumed to be unilamellar, and large deviations from this assumption would dramatically impact interpretations from many model membrane experiments. In this work, using three probe-free methods—small-angle X-ray and neutron scattering (SAXS and SANS) and cryogenic electron microscopy (cryoEM)—we report unambiguous evidence of extensive multilamellarity in extruded vesicles composed of neutral phosphatidylcholine lipids, including for the common case of neutral lipids dispersed in physiological buffer and extruded through 100 nm diameter pores. In such preparations, only ~35% of lipids are externally accessible, and this fraction is highly dependent on preparation conditions. Charged lipids promote unilamellarity, as does decreasing solvent ionic strength, indicating the importance of electrostatic interactions in determining the lamellarity of extruded vesicles. Smaller extrusion pore sizes also robustly increase the fraction of unilamellar vesicles, suggesting a role for membrane bending. Taken together, these observations suggest a mechanistic model for extrusion, wherein formation of unilamellar vesicles involves competition between bilayer bending and adhesion energies. The findings presented here have wide-ranging implications for the design and interpretation of model membrane studies, especially ensemble-averaged observations relying on the assumption of unilamellarity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Counterflow ignition of n-butanol at atmospheric and elevated pressures

Vital to the development of predictive combustion models is a robust understanding of the coupled effects of chemical kinetics and convective–diffusive transport at both atmospheric and elevated pressures. The current study describes a new variable-pressure non-premixed counterflow ignition experiment designed to address the need for well-characterized reference data to validate such models under conditions sensitive to both chemical and transport processes. A comprehensive characterization of system behavior is provided to demonstrate boundary condition and ignition quality as well as adherence to the assumption of quasi-one-dimensionality, and suggest limitations and best practices for counterflow ignition experiments. This effort reveals that the counterflow ignition experiment requires special attention to ignition location in order to ensure that the assumption of quasi-one-dimensionality is valid, particularly at elevated pressures. This experimental tool is then applied to the investigation of n-butanol for pressures of 1–4 atm, pressure-weighted strain rates of 200–400 s -1 , and fuel mole fractions of 0.05–0.25. Results are simulated using two n-butanol models available in the literature and used to validate and assess model performance. Comparison of experimental and numerical ignition results for n-butanol demonstrates that while existing models largely capture the trends observed with varying pressure, strain rate, and fuel loading, the models universally over-predict experimental ignition temperatures. While several transport coefficients are found to exhibit order-of-magnitude or greater sensitivities relative to reaction rates, variation of transport parameters is not able to account for the large deviations observed between experimental and numerical results. Additional comparison of ignition kernel structure and fuel breakdown pathways between two literature models suggests that an under-prediction in the radical pool growth with respect to temperature variation may be responsible for both the deviation from the experimental results and the discrepancy in ignition temperature results observed between models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Robust control of the current profile and plasma energy in EAST

Integrated control of the toroidal current density profile, or alternatively the q-profile, and plasma stored energy is essential to achieve advanced plasma scenarios characterized by high plasma confinement, magnetohydrodynamics stability, and noninductively driven plasma current. The q-profile evolution is closely related to the evolution of the poloidal magnetic flux profile, whose dynamics is modeled by a nonlinear partial differential equation (PDE) referred to as the magnetic-flux diffusion equation (MDE). The MDE prediction depends heavily on the chosen models for the electron temperature, plasma resistivity, and non-inductive current drives. To aid control synthesis, control-oriented models for these plasma quantities are necessary to make the problem tractable. However, a relatively large deviation between the predictions by these control-oriented models and experimental data is not uncommon. For this reason, the electron temperature, plasma resistivity, and non-inductive current drives are modeled for control synthesis in this work as the product of an “uncertain” reference profile and a nonlinear function of the different auxiliary heating and current-drive (H&CD) source powers and the total plasma current. The uncertainties are quantified in such a way that the family of models arising from the modeling process is able to capture the q-profile and plasma stored energy dynamics from a typical EAST shot. A control-oriented nonlinear PDE model is developed by combining the MDE with the “uncertain” models for the electron temperature, plasma resistivity, and non-inductive current drives. This model is then rewritten into a control framework to design a controller that is robust against the modeled uncertainties. The resulting controller utilizes EAST's H&CD powers and total plasma current to regulate the q profile and plasma stored energy even when mismatches between modeled and actual dynamics are present. The effectiveness of the controller is demonstrated through nonlinear simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Modeling and Experimental Validation of a Direct-Contact Counter-Flow Fluidized Bed Heat Exchanger for Thermal Energy Storage (TES) Applications

Particle-based thermal energy storage (TES) systems are an emerging energy storage technology. The technological advances have reduced costs, making TES more competitive and reliable in the marketplace but an efficient and reliable operation is heavily dependent on coherent heat transfer between air to particles or vice versa. The particle-based TES technologies provide an intermediate system that can store energy for short (0-10 h), long (10-200 h) and seasonal (> 200 h) timescales. The TES systems store energy by converting electricity to thermal energy; electricity can be directly sourced intermittent generation technologies and/or the grid, helping manage peak loads and other mismatches in supply and demand. The overall efficiency of the TES system depends on the performance of system components (particle storage silos and particle transfer mechanism etc.). The particle heat exchanger is one of the key system components that affects the system efficiency. The pressurized fluidized bed heat exchanger (PFB HX) performance is challenging to predict due to the chaotic behavior of particle and fluid interaction. This research presents a computational study of a novel direct-contact, counter-flow and air-to-particles PFB HX, that contributes in advancing the particle-based long-duration TES technologies. For the current analysis an unsteady Eulerian-Eulerian CFD model was developed and validated against experiments performed at the National Laboratory of the Rockies for two particle sizes (600 ..mu..m and 825 ..mu..m ). Following validation, parametric simulations were conducted to evaluate the effects of interphase drag models (Syamlal-O'Brien and Gidaspow), particle size, bed height and the influence of a frictional-viscosity term on hydrodynamics and heat transfer between the air & particles. The key findings from the analysis are: (1) for the studied operating window Syamlal-O'Brien provides superior agreement with measured gas temperatures (errors generally < 10%) while Gidaspow shows large deviations for the coarse particle case; (2) model predictions are most sensitive in the lower 0.2 m above the air distributor where bubble initiation and local mixing dominate interphase heat transfer; (3) representation of the distributor (number of inlet ports) materially affects predicted local mixing and temperature stratification; and (4) the Eulerian-Eulerian framework reproduces bulk thermal trends but shows regime dependent limitations for coarse particles, motivating mesoscale informed closures for scale-up analysis for future studies. These results provide validated guidance for drag selection and distributor design in particle-based thermal energy storage applications. Collectively, the validated model and parametric results quantify key drivers of PHB-HX performance and provide practical guidance for design and optimization. The results provide confidence in the model predictability and provide a step forward to improve on heat exchange performance. The demonstrated performance and modeling approach support the deployment and further development of this novel PHB-HX concept for robust, particle-based long-duration thermal energy storage systems.

25 ENERGY STORAGE↗
DOE OSTI2021

Isotopic dependence of the frequency of optical vibrations in molybdenum monohydride

We show that it is currently known that three hydrides – PdH x , MoH x , and TiH x – exhibit an inverse isotope effect in superconductivity. Namely, the phase with a heavier hydrogen isotope, deuterium, has a higher critical temperature. Hydrides and deuterides of palladium have intensively been studied both experimentally and theoretically, but the origin of the isotope effect has not been established with certainty. The commonly accepted explanation is that the effect is likely to be due to the strong anharmonicity of the optical hydrogen vibrations, which was considered to be responsible for the large deviation of the ratio of the fundamental optical frequencies ω H /ω D = 1.51 from the harmonic value $\sqrt{2}$ ≈ 1.41. In the present paper, powder samples of MoH 1.1(1) and MoD 1.07(3) were synthesized under a hydrogen / deuterium pressure of several gigapascals and studied by inelastic neutron scattering (INS) at ambient pressure and T = 10 K. The INS study demonstrated that optical vibrations of H atoms in MoH 1.1 and D atoms in MoD 1.07 are harmonic and the ratio of fundamental optical frequencies ω H /ω D = 1.44 is close to the harmonic value $\sqrt{2}$ ≈ 1.41. This shows that anharmonicity is not a necessary condition for the presence of the inverse isotope effect. The MoD 1.07 sample was additionally studied by neutron diffraction (ND) at ambient pressure and T = 100 K. In agreement with previous ND results for MoH 1.2 , the ND study of MoD 1.07 showed that deuterium atoms occupied almost all octahedral interstitial sites in its hexagonal close-packed metal lattice and formed a NiAs-type crystal structure with the composition close to MoD. The overstoichiometric composition MoD 1.07 of the deuteride is likely to result from a small fraction D/Mo ~ 0.07 of deuterium atoms partially occupying the tetrahedral interstices.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Drag of a heated sphere at low Reynolds numbers in the absence of buoyancy

Fully resolved simulations are used to quantify the effects of heat transfer in the absence of buoyancy on the drag of a spatially fixed heated spherical particle at low Reynolds numbers ( $Re$ ) in the range $10^{-3}\leqslant Re\leqslant 10$ in a variable-property fluid. The case where buoyancy is present is analysed in a subsequent paper. This analysis is carried out without making any assumptions on the amount of heat addition from the sphere and thus encompasses both the heating regime where the Boussinesq approximation holds and the regime where it breaks down. The particle is assumed to have a low Biot number, which means that the particle is uniformly at the same temperature and has no internal temperature gradients. Large deviations in the value of the drag coefficient as the temperature of the sphere increases are observed. When $Re , these deviations are explained using a low-Mach-number perturbation analysis as irrotational corrections to a Stokes–Oseen base flow. Correlations for the drag and Nusselt number of a heated sphere are proposed for the range of Reynolds numbers $10^{-3}\leqslant Re\leqslant 10$ which fit the computationally obtained values with less than 1 % and 3 % errors, respectively. These correlations can be used in simulations of gas–solid flows where the accuracy of the drag law affects the prediction of the overall flow behaviour. Finally, an analogy to incompressible flow over a modified sphere is demonstrated.

Mechanics↗
DOE OSTI2019

Linear theory of electron-plasma waves at arbitrary collisionality

We report the dynamics of electron-plasma waves is described at arbitrary collisionality by considering the full Coulomb collision operator. The description is based on a Hermite–Laguerre decomposition of the velocity dependence of the electron distribution function. The damping rate, frequency and eigenmode spectrum of electron-plasma waves are found as functions of the collision frequency and wavelength. A comparison is made between the collisionless Landau damping limit, the Lenard–Bernstein and Dougherty collision operators and the electron–ion collision operator, finding large deviations in the damping rates and eigenmode spectra. A purely damped entropy mode, characteristic of a plasma where pitch-angle scattering effects are dominant with respect to collisionless effects, is shown to emerge numerically, and its dispersion relation is analytically derived. It is shown that such a mode is absent when simplified collision operators are used, and that like-particle collisions strongly influence the damping rate of the entropy mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Biased Spin-State Energetics of Fe(II) Molecular Complexes within Density-Functional Theory and the Linear-Response Hubbard U Correction

The spin-state energetics of six Fe(II) molecular complexes are computed using the linear-response Hubbard U approach within DFT. The adiabatic energy differences, ΔE H-L , between the high-spin (S = 2) and the low-spin (S = 0) states are computed and compared with accurate-coupled cluster-corrected CASPT2 results. We show that DFT+U fails in correctly capturing the ground state for strong-field ligands yielding ΔE H-L that are almost constant throughout the molecular series. This bias toward high spin together with the metal/ligand charge transfer upon U correction are here quantified and explained using molecular orbital diagrams involving both σ- and π-bonding interactions. Furthermore, with increasing ligand-field strengths this bias also increases owing to the stronger molecular character of the metal/ligand Kohn–Sham orbitals thus resulting in large deviations from the reference larger than 4 eV. Smaller values of U can be employed to mitigate this effect and recover the right energetics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Correlating Thermodynamic and Kinetic Hydricities of Rhenium Hydrides

The kinetics of hydride transfer from Re( R bpy)(CO) 3 H (bpy = 4,4'-R-2,2'-bipyridine; R = OMe, t Bu, Me, H, Br, COOMe, CF 3 ) to CO 2 and seven different cationic N-heterocycles were determined. Additionally, the thermodynamic hydricities of complexes of the type Re( R bpy)(CO) 3 H were established primarily using computational methods. Linear free-energy relationships (LFERs) derived by correlating thermodynamic and kinetic hydricities indicate that, in general, the rate of hydride transfer increases as the thermodynamic driving force for the reaction increases. Kinetic isotope effects range from inverse for hydride transfer reactions with a small driving force to normal for reactions with a large driving force. Hammett analysis indicates that hydride transfer reactions with greater thermodynamic driving force are less sensitive to changes in the electronic properties of the metal hydride, presumably because there is less buildup of charge in the increasingly early transition state. Bronsted α values were obtained for a range of hydride transfer reactions and along with DFT calculations suggest the reactions are concerted, which enables the use of Marcus theory to analyze hydride transfer reactions involving transition metal hydrides. It is notable, however, that even slight perturbations in the steric properties of the Re hydride or the hydride acceptor result in large deviations in the predicted rate of hydride transfer based on thermodynamic driving forces. This indicates that thermodynamic considerations alone cannot be used to predict the rate of hydride transfer, which has implications for catalyst design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Supersonic plasma turbulence in the laboratory

The properties of supersonic, compressible plasma turbulence determine the behavior of many terrestrial and astrophysical systems. In the interstellar medium and molecular clouds, compressible turbulence plays a vital role in star formation and the evolution of our galaxy. Observations of the density and velocity power spectra in the Orion B and Perseus molecular clouds show large deviations from those predicted for incompressible turbulence. Hydrodynamic simulations attribute this to the high Mach number in the interstellar medium (ISM), although the exact details of this dependence are not well understood. Here we investigate experimentally the statistical behavior of boundary-free supersonic turbulence created by the collision of two laser-driven high-velocity turbulent plasma jets. The Mach number dependence of the slopes of the density and velocity power spectra agree with astrophysical observations, and supports the notion that the turbulence transitions from being Kolmogorov-like at low Mach number to being more Burgers-like at higher Mach numbers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Non–linear bonding trends in maleonitrile-1,2–dithiolate complexes of the transuranium actinides

The trivalent actinides are produced in the nuclear fuel cycle during power production and provide the largest long-term radiation dose in used nuclear fuel. It is ideal for these elements to be removed from used nuclear fuel for disposal and a necessity for fuel recycling. A key challenge to this is the similarity of chemical behavior of the trivalent actinides to the lanthanides that are also present as fission products in used fuel. Thus far, some of the most effective separations of actinides from lanthanides utilise chelating agents containing sulfur moieties such as dithiophosphinates that selectively bind to actinide ions because of a greater bond covalency relative to lanthanide ions. Typically, greater differences between actinide and lanthanide ions are observable the more ligands and chelators bonds have a covalent character. Here, a series of complexes of the trivalent actinides Np(III) through Cf(III) (excluding Bk(III)) with maleonitrile-1,2-dithiolate (mnt 2– ) are synthesized along with their lanthanide counterparts (La(III) – Nd(III), Sm(III) – Gd(III), Dy(III)), in order to characterize the nature of chemical bonds with these metal ions and a polarizable, non-innocent, sulfur-donor ligand. The metal-sulfur bonds in these complexes trend shorter than measured for lanthanides with equivalent ionic radii. However, particularly large deviations are observed in the neptunium and plutonium complexes in both structure and bonding, resulting in a nonlinear bond length trendline for the actinide series. Density Functional Theory (DFT) calculations with Quantum Theory of Atoms in Molecules (QTAIM) and Natural Bond Order (NBO) analyses indicate that for the neptunium and plutonium complexes, the presence of increased 5f-orbital participation, energy degeneracy of the metal and ligand orbitals, and the structure packing result in shortened M–S bonds. The stabilization of the energy of the 5f-orbitals and the decrease in f-contribution to bonding orbitals in the later actinides results in structural properties more similar to the lanthanide complexes.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2026

Stability Constants of Lanthanide-nitrate Complexes in Aqueous Solutions: A Theoretical Study

Calculating stability constants for lanthanide--nitrate complexes in aqueous solution is challenging due to the complex free-energy landscapes of the participating species. In this work, we compare cluster-continuum solvation and condensed-phase approaches using universal machine learning interatomic potentials (MLIPs) for determining the stability constants of lanthanide--nitrate complexes in aqueous solutions. Within the cluster--continuum solvation framework at the B3LYP level of theory, reactions involving lanthanide coordination numbers of both 8 and 9 are found to be relevant. After an empirical linear free-energy correction, the cluster--continuum results fall on the same order-of-magnitude scale as the experimental stability constants. By contrast, MACE-MP0 MLIP underestimates lanthanide hydration numbers and gives PMF-derived stability constants with large deviations from experiment, whereas MACE-MATPES-R2SCAN improves both hydration structure and the stability-constant scale but still does not quantitatively reproduce the detailed lanthanide trend. Across both approaches, nitrate binding is best viewed as a labile coordination motif rather than a fixed mono- or bidentate structure, with hydration-shell structure influencing which configurations are favored. Overall, the cluster--continuum calculations provide a practical semi-quantitative baseline for the experimental stability-constant scale, while the explicit-solvent MLIP benchmarks show clear progress from MACE-MP0 to MACE-MATPES-r2SCAN but also highlight the need for lanthanide-targeted training, fine-tuning, or improved long-range and polarization treatments to obtain predictive thermodynamics for complex aqueous lanthanide chemistry.

Dinpajooh, Mohammadhasan↗
DOE OSTI2021

Structure and function of aerotolerant, multiple-turnover THI4 thiazole synthases

Plant and fungal THI4 thiazole synthases produce the thiamin thiazole moiety in aerobic conditions via a single-turnover suicide reaction that uses an active-site Cys residue as sulfur donor. Multiple-turnover (i.e. catalytic) THI4s lacking an active-site Cys (non-Cys THI4s) that use sulfide as sulfur donor have been biochemically characterized -- but only from archaeal methanogens that are anaerobic, O 2 -sensitive hyperthermophiles from sulfide-rich habitats. These THI4s prefer iron as cofactor. A survey of prokaryote genomes uncovered non-Cys THI4s in aerobic mesophiles from sulfide-poor habitats, suggesting that multiple-turnover THI4 operation is possible in aerobic, mild, low-sulfide conditions. This was confirmed by testing 23 representative non-Cys THI4s for complementation of an Escherichia coli ΔthiG thiazole auxotroph in aerobic conditions. Sixteen were clearly active, and more so when intracellular sulfide level was raised by supplying Cys, demonstrating catalytic function in the presence of O 2 at mild temperatures and indicating use of sulfide or a sulfide metabolite as sulfur donor. Comparative genomic evidence linked non-Cys THI4s with proteins from families that bind, transport, or metabolize cobalt or other heavy metals. The crystal structure of the aerotolerant bacterial Thermovibrio ammonificans THI4 was determined to probe the molecular basis of aerotolerance. The structure suggested no large deviations compared with the structures of THI4s from O 2 -sensitive methanogens, but is consistent with an alternative catalytic metal. Together with complementation data, use of cobalt rather than iron was supported. We conclude that catalytic THI4s can indeed operate aerobically and that the metal cofactor inserted is a likely natural determinant of aerotolerance.

59 BASIC BIOLOGICAL SCIENCES↗