Modeling of biocatalytic reactions: A workflow for model calibration, selection and validation using Bayesian statistics

We present a workflow for kinetic modeling of biocatalytic reactions which combines methods from Bayesian learning and uncertainty quantification for model calibration, model selection, evaluation, and model reduction in a consistent statistical framework. Our workflow is particularly tailored to sparse data settings in which a considerable variability of the parameters remains after the models have been adapted to available data, a ubiquitous problem in many real‐world applications. Our workflow is exemplified on an enzyme‐catalyzed two‐substrate reaction mechanism describing the symmetric carboligation of 3,5‐dimethoxy‐benzaldehyde to (R )‐3,3′,5,5′‐tetramethoxybenzoin catalyzed by benzaldehyde lyase from Pseudomonas fluorescens . Results indicate a substrate‐dependent inactivation of enzyme, which is in accordance with other recent studies

Alteration of the Route to Menaquinone towards Isochorismate-Derived Metabolites

Chorismate and isochorismate constitute branch-point intermediates in the biosynthesis of many aromatic metabolites in microorganisms and plants. To obtain unnatural compounds, we modified the route to menaquinone in Escherichia coli. We propose a model for the binding of isochorismate to the active site of MenD ((1R,2S, 5S,6S)-2-succinyl-5-enolpyruvyl-6-hydroxycyclohex-3-ene-1-carboxylate (SEPHCHC) synthase) that explains the outcome of the native reaction with α-ketoglutarate. We have rationally designed variants of MenD for the conversion of several isochorismate analogues. The double-variant Asn117Arg–Leu478Thr preferentially converts (5S,6S)-5,6-dihydroxycyclohexa-1,3-diene-1-carboxylate (2,3-trans-CHD), the hydrolysis product of isochorismate, with a >70-fold higher ratio than that for the wild type. The single-variant Arg107Ile uses (5S,6S)-6-amino-5-hydroxycyclohexa-1,3-diene-1-carboxylate (2,3-trans-CHA) as substrate with >6-fold conversion compared to wild-type MenD. The novel compounds have been made accessible in vivo (up to 5.3 g L−1). Unexpectedly, as the identified residues such as Arg107 are highly conserved (>94 %), some of the designed variations can be found in wild-type SEPHCHC synthases from other bacteria (Arg107Lys, 0.3 %). This raises the question for the possible natural occurrence of as yet unexplored branches of the shikimate pathway.Fil: Fries, Alexander Erich. Consejo Nacional de Investigaciones Científicas y Técnicas. Instituto de Ciencias de la Tierra y Ambientales de La Pampa. Universidad Nacional de La Pampa. Facultad de Ciencias Exactas y Naturales. Instituto de Ciencias de la Tierra y Ambientales de La Pampa; Argentina. Albert Ludwigs University of Freiburg; AlemaniaFil: Mazzaferro, Laura. Consejo Nacional de Investigaciones Científicas y Técnicas. Instituto de Ciencias de la Tierra y Ambientales de La Pampa. Universidad Nacional de La Pampa. Facultad de Ciencias Exactas y Naturales. Instituto de Ciencias de la Tierra y Ambientales de La Pampa; Argentina. Albert Ludwigs University of Freiburg; AlemaniaFil: Grüning, Björn. Albert Ludwigs University of Freiburg; AlemaniaFil: Bisel, Philippe. Albert Ludwigs University of Freiburg; AlemaniaFil: Stibal, Karin. Albert Ludwigs University of Freiburg; AlemaniaFil: Buchholz, Patrick C. F.. University of Stuttgart; AlemaniaFil: Pleiss, Jürgen. Universität Stuttgart;Fil: Sprenger, Georg A.. Universität Stuttgart;Fil: Müller, Michael. Albert Ludwigs University of Freiburg; Alemani

Search for the standard model Higgs boson decaying to a bb pair in events with one charged lepton and large missing transverse energy using the full CDF data set

We present a search for the standard model Higgs boson produced in association with a W boson in sqrt(s) = 1.96 TeV p-pbar collision data collected with the CDF II detector at the Tevatron corresponding to an integrated luminosity of 9.45 fb-1. In events consistent with the decay of the Higgs boson to a bottom-quark pair and the W boson to an electron or muon and a neutrino, we set 95% credibility level upper limits on the WH production cross section times the H->bb branching ratio as a function of Higgs boson mass. At a Higgs boson mass of 125 GeV/c2 we observe (expect) a limit of 4.9 (2.8) times the standard model value.Comment: Submitted to Phys. Rev. Lett (v2 contains clarifications suggested by PRL

Search for the standard model Higgs boson decaying to a bbˉb\bar{b} pair in events with no charged leptons and large missing transverse energy using the full CDF data set

We report on a search for the standard model Higgs boson produced in association with a vector boson in the full data set of proton-antiproton collisions at $\sqrt{s} = 1.96$ TeV recorded by the CDF II detector at the Tevatron, corresponding to an integrated luminosity of 9.45 fb$^{-1}$. We consider events having no identified charged lepton, a transverse energy imbalance, and two or three jets, of which at least one is consistent with originating from the decay of a $b$ quark. We place 95% credibility level upper limits on the production cross section times standard model branching fraction for several mass hypotheses between 90 and $150 \mathrm{GeV}/c^2$. For a Higgs boson mass of $125 \mathrm{GeV}/c^2$, the observed (expected) limit is 6.7 (3.6) times the standard model prediction.Comment: Accepted by Phys. Rev. Let

Search for the standard model Higgs boson decaying to a bb pair in events with two oppositely-charged leptons using the full CDF data set

We present a search for the standard model Higgs boson produced in association with a Z boson in data collected with the CDF II detector at the Tevatron, corresponding to an integrated luminosity of 9.45/fb. In events consistent with the decay of the Higgs boson to a bottom-quark pair and the Z boson to electron or muon pairs, we set 95% credibility level upper limits on the ZH production cross section times the H -> bb branching ratio as a function of Higgs boson mass. At a Higgs boson mass of 125 GeV/c^2 we observe (expect) a limit of 7.1 (3.9) times the standard model value.Comment: To be submitted to Phys. Rev. Let

Measurement of the cosmic ray spectrum above 4×10184{\times}10^{18} eV using inclined events detected with the Pierre Auger Observatory

A measurement of the cosmic-ray spectrum for energies exceeding $4{\times}10^{18}$ eV is presented, which is based on the analysis of showers with zenith angles greater than $60^{\circ}$ detected with the Pierre Auger Observatory between 1 January 2004 and 31 December 2013. The measured spectrum confirms a flux suppression at the highest energies. Above $5.3{\times}10^{18}$ eV, the "ankle", the flux can be described by a power law $E^{-\gamma}$ with index $\gamma=2.70 \pm 0.02 \,\text{(stat)} \pm 0.1\,\text{(sys)}$ followed by a smooth suppression region. For the energy ($E_\text{s}$) at which the spectral flux has fallen to one-half of its extrapolated value in the absence of suppression, we find $E_\text{s}=(5.12\pm0.25\,\text{(stat)}^{+1.0}_{-1.2}\,\text{(sys)}){\times}10^{19}$ eV.Comment: Replaced with published version. Added journal reference and DO

Energy Estimation of Cosmic Rays with the Engineering Radio Array of the Pierre Auger Observatory

The Auger Engineering Radio Array (AERA) is part of the Pierre Auger Observatory and is used to detect the radio emission of cosmic-ray air showers. These observations are compared to the data of the surface detector stations of the Observatory, which provide well-calibrated information on the cosmic-ray energies and arrival directions. The response of the radio stations in the 30 to 80 MHz regime has been thoroughly calibrated to enable the reconstruction of the incoming electric field. For the latter, the energy deposit per area is determined from the radio pulses at each observer position and is interpolated using a two-dimensional function that takes into account signal asymmetries due to interference between the geomagnetic and charge-excess emission components. The spatial integral over the signal distribution gives a direct measurement of the energy transferred from the primary cosmic ray into radio emission in the AERA frequency range. We measure 15.8 MeV of radiation energy for a 1 EeV air shower arriving perpendicularly to the geomagnetic field. This radiation energy -- corrected for geometrical effects -- is used as a cosmic-ray energy estimator. Performing an absolute energy calibration against the surface-detector information, we observe that this radio-energy estimator scales quadratically with the cosmic-ray energy as expected for coherent emission. We find an energy resolution of the radio reconstruction of 22% for the data set and 17% for a high-quality subset containing only events with at least five radio stations with signal.Comment: Replaced with published version. Added journal reference and DO

Measurement of the Radiation Energy in the Radio Signal of Extensive Air Showers as a Universal Estimator of Cosmic-Ray Energy

We measure the energy emitted by extensive air showers in the form of radio emission in the frequency range from 30 to 80 MHz. Exploiting the accurate energy scale of the Pierre Auger Observatory, we obtain a radiation energy of 15.8 \pm 0.7 (stat) \pm 6.7 (sys) MeV for cosmic rays with an energy of 1 EeV arriving perpendicularly to a geomagnetic field of 0.24 G, scaling quadratically with the cosmic-ray energy. A comparison with predictions from state-of-the-art first-principle calculations shows agreement with our measurement. The radiation energy provides direct access to the calorimetric energy in the electromagnetic cascade of extensive air showers. Comparison with our result thus allows the direct calibration of any cosmic-ray radio detector against the well-established energy scale of the Pierre Auger Observatory.Comment: Replaced with published version. Added journal reference and DOI. Supplemental material in the ancillary file