Uncertainty Quantification Of A Semi-Supervised Support Vector Regression Reservoir Model

Uncertainty quantification of petroleum reservoir models is one of the present challenges, which is usually approached with a wide range of geostatistical tools linked with statistical optimisation or/and inference algorithms. Recent advances in machine learning offer a novel approach to model spatial distribution of petrophysical properties in complex reservoirs alternative to geostatistics. The approach is based of semisupervised learning, which handles both ?labelled? observed data and ?unlabelled? data, which have no measured value but describe prior knowledge and other relevant data in forms of manifolds in the input space where the modelled property is continuous. Proposed semi-supervised Support Vector Regression (SVR) model has demonstrated its capability to represent realistic geological features and describe stochastic variability and non-uniqueness of spatial properties. On the other hand, it is able to capture and preserve key spatial dependencies such as connectivity of high permeability geo-bodies, which is often difficult in contemporary petroleum reservoir studies. Semi-supervised SVR as a data driven algorithm is designed to integrate various kind of conditioning information and learn dependences from it. The semi-supervised SVR model is able to balance signal/noise levels and control the prior belief in available data. In this work, stochastic semi-supervised SVR geomodel is integrated into Bayesian framework to quantify uncertainty of reservoir production with multiple models fitted to past dynamic observations (production history). Multiple history matched models are obtained using stochastic sampling and/or MCMC-based inference algorithms, which evaluate posterior probability distribution. Uncertainty of the model is described by posterior probability of the model parameters that represent key geological properties: spatial correlation size, continuity strength, smoothness/variability of spatial property distribution. The developed approach is illustrated with a fluvial reservoir case. The resulting probabilistic production forecasts are described by uncertainty envelopes. The paper compares the performance of the models with different combinations of unknown parameters and discusses sensitivity issues.

Quantification of abscisic Acid in a single maize root.

Quantitative analyses of abscisic acid in the elongating zone of a single maize root (Zea mays L. cv LG 11) were performed by gas chromatography-mass spectrometry using negative chemical ion ionization. Data showed that the more abscisic acid, the slower the growth, but a large dispersion of individual values was observed. We assume that abscisic acid is perhaps not correlated only to the growth rate.

Analysis and quantification of vitamin D metabolites in serum by ultra-performance liquid chromatography coupled to tandem mass spectrometry and high-resolution mass spectrometry: a method comparison and validation.

RATIONALE: The aim of the work was to develop and validate a method for the quantification of vitamin D metabolites in serum using ultra-high-pressure liquid chromatography coupled to mass spectrometry (LC/MS), and to validate a high-resolution mass spectrometry (LC/HRMS) approach against a tandem mass spectrometry (LC/MS/MS) approach using a large clinical sample set. METHODS: A fast, accurate and reliable method for the quantification of the vitamin D metabolites, 25-hydroxyvitamin D2 (25OH-D2) and 25-hydroxyvitamin D3 (25OH-D3), in human serum was developed and validated. The C3 epimer of 25OH-D3 (3-epi-25OH-D3) was also separated from 25OH-D3. The samples were rapidly prepared via a protein precipitation step followed by solid-phase extraction (SPE) using an HLB μelution plate. Quantification was performed using both LC/MS/MS and LC/HRMS systems. RESULTS: Recovery, matrix effect, inter- and intra-day reproducibility were assessed. Lower limits of quantification (LLOQs) were determined for both 25OH-D2 and 25OH-D3 for the LC/MS/MS approach (6.2 and 3.4 µg/L, respectively) and the LC/HRMS approach (2.1 and 1.7 µg/L, respectively). A Passing & Bablok fit was determined between both approaches for 25OH-D3 on 662 clinical samples (1.11 + 1.06x). It was also shown that results can be affected by the inclusion of the isomer 3-epi-25OH-D3. CONCLUSIONS: Quantification of the relevant vitamin D metabolites was successfully developed and validated here. It was shown that LC/HRMS is an accurate, powerful and easy to use approach for quantification within clinical laboratories. Finally, the results here suggest that it is important to separate 3-epi-25OH-D3 from 25OH-D3. Copyright © 2012 John Wiley & Sons, Ltd.

Handling macromolecule signals in the quantification of the neurochemical profile.

In vivo localized proton magnetic resonance spectroscopy (1H MRS) became a powerful and unique technique to non-invasively investigate brain metabolism of rodents and humans. The main goal of 1H MRS is the reliable quantification of concentrations of metabolites (neurochemical profile) in a well-defined region of the brain. The availability of very high magnetic field strengths combined with the possibility of acquiring spectra at very short echo time have dramatically increased the number of constituents of the neurochemical profile. The quantification of spectra measured at short echo times is complicated by the presence of macromolecule signals of particular importance at high magnetic fields. An error in the macromolecule estimation can lead to substantial errors in the obtained neurochemical profile. The purpose of the present review is to overview methods of high field 1H MRS with a focus on the metabolite quantification, in particular in handling signals of macromolecules. Three main approaches of handling signals of macromolecules are described, namely mathematical estimation of macromolecules, measurement of macromolecules in vivo, and direct acquisition of the in vivo spectrum without the contribution of macromolecules.

Quantification of 4 antidepressants and a metabolite by LC-MS for therapeutic drug monitoring.

A liquid chromatography method coupled to mass spectrometry was developed for the quantification of bupropion, its metabolite hydroxy-bupropion, moclobemide, reboxetine and trazodone in human plasma. The validation of the analytical procedure was assessed according to Société Française des Sciences et Techniques Pharmaceutiques and the latest Food and Drug Administration guidelines. The sample preparation was performed with 0.5mL of plasma extracted on a cation-exchange solid phase 96-well plate. The separation was achieved in 14min on a C18 XBridge column (2.1mm×100mm, 3.5μm) using a 50mM ammonium acetate pH 9/acetonitrile mobile phase in gradient mode. The compounds of interest were analysed in the single ion monitoring mode on a single quadrupole mass spectrometer working in positive electrospray ionisation mode. Two ions were selected per molecule to increase the number of identification points and to avoid as much as possible any false positives. Since selectivity is always a critical point for routine therapeutic drug monitoring, more than sixty common comedications for the psychiatric population were tested. For each analyte, the analytical procedure was validated to cover the common range of concentrations measured in plasma samples: 1-400ng/mL for reboxetine and bupropion, 2-2000ng/mL for hydroxy-bupropion, moclobemide, and trazodone. For all investigated compounds, reliable performance in terms of accuracy, precision, trueness, recovery, selectivity and stability was obtained. One year after its implementation in a routine process, this method demonstrated a high robustness with accurate values over the wide concentration range commonly observed among a psychiatric population.

Generic approach for the sensitive absolute quantification of large undigested peptides in plasma using a particular liquid chromatography-mass spectrometry setup.

A generic LC-MS approach for the absolute quantification of undigested peptides in plasma at mid-picomolar levels is described. Nine human peptides namely, brain natriuretic peptide (BNP), substance P (SubP), parathyroid hormone 1-34 (PTH), C-peptide, orexines A and B (Orex-A and -B), oxytocin (Oxy), gonadoliberin-1 (gonadothropin releasing-hormone or luteinizing hormone-releasing hormone, LHRH) and α-melanotropin (α-MSH) were targeted. Plasma samples were extracted via a 2-step procedure: protein precipitation using 1vol of acetonitrile followed by ultrafiltration of supernatants on membranes with a MW cut-off of 30 kDa. By applying a specific LC-MS setup, large volumes of filtrates (e.g., 2×750 μL) were injected and the peptides were trapped on a 1mm i.d.×10 mm length C8 column using a 10× on-line dilution. Then, the peptides were back-flushed and a second on-line dilution (2×) was applied during the transfer step. The refocalized peptides were resolved on a 0.3mm i.d. C18 analytical column. Extraction recovery, matrix effect and limits of detection were evaluated. Our comprehensive protocol demonstrates a simple and efficient sample preparation procedure followed by the analysis of peptides with limits of detection in the mid-picomolar range. This generic approach can be applied for the determination of most therapeutic peptides and possibly for endogenous peptides with latest state-of-the-art instruments.

Quantification of typical antipsychotics in human plasma by ultra-high performance liquid chromatography tandem mass spectrometry for therapeutic drug monitoring.

A selective and sensitive method was developed for the simultaneous quantification of seven typical antipsychotic drugs (cis-chlorprothixene, flupentixol, haloperidol, levomepromazine, pipamperone, promazine and zuclopenthixol) in human plasma. Ultra-high performance liquid chromatography (UHPLC) was used for complete separation of the compounds in less than 4.5min on an Acquity UPLC BEH C18 column (2.1mm×50mm; 1.7μm), with a gradient elution of ammonium formate buffer pH 4.0 and acetonitrile at a flow rate of 400μl/min. Detection was performed on a tandem quadrupole mass spectrometer (MS/MS) equipped with an electrospray ionization interface. A simple protein precipitation procedure with acetonitrile was used for sample preparation. Thanks to the use of stable isotope-labeled internal standards for all analytes, internal standard-normalized matrix effects were in the range of 92-108%. The method was fully validated to cover large concentration ranges of 0.2-90ng/ml for haloperidol, 0.5-90ng/ml for flupentixol, 1-450ng/ml for levomepromazine, promazine and zuclopenthixol and 2-900ng/ml for cis-chlorprothixene and pipamperone. Trueness (89.1-114.8%), repeatability (1.8-9.9%), intermediate precision (1.9-16.3%) and accuracy profiles (<30%) were in accordance with the latest international recommendations. The method was successfully used in our laboratory for routine quantification of more than 500 patient plasma samples for therapeutic drug monitoring. To the best of our knowledge, this is the first UHPLC-MS/MS method for the quantification of the studied drugs with a sample preparation based on protein precipitation.

Characterization of the ligand-binding site of the serotonin 5-HT3 receptor: the role of glutamate residues 97, 224, AND 235.

Ligand-gated ion channels of the Cys loop family are receptors for small amine-containing neurotransmitters. Charged amino acids are strongly conserved in the ligand-binding domain of these receptor proteins. To investigate the role of particular residues in ligand binding of the serotonin 5-HT3AS receptor (5-HT3R), glutamate amino acid residues at three different positions, Glu97, Glu224, and Glu235, in the extracellular N-terminal domain were substituted with aspartate and glutamine using site-directed mutagenesis. Wild type and mutant receptor proteins were expressed in HEK293 cells and analyzed by electrophysiology, radioligand binding, fluorescence measurements, and immunochemistry. A structural model of the ligand-binding domain of the 5-HT3R based on the acetylcholine binding protein revealed the position of the mutated amino acids. Our results demonstrate that mutations of Glu97, distant from the ligand-binding site, had little effect on the receptor, whereas mutations Glu224 and Glu235, close to the predicted binding site, are indeed important for ligand binding. Mutations E224Q, E224D, and E235Q decreased EC50 and Kd values 5-20-fold, whereas E235D was functionally expressed at a low level and had a more than 100-fold increased EC50 value. Comparison of the fluorescence properties of a fluorescein-labeled antagonist upon binding to wild type 5-HT3R and E235Q, allowed us to localize Glu235 within a distance of 1 nm around the ligand-binding site, as proposed by our model.

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

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

Differential effects of glutamate transporter inhibitors on the global electrophysiological response of astrocytes to neuronal stimulation

Astrocytes are responsible for regulating extracellular levels of glutamate and potassium during neuronal activity. Glutamate clearance is handled by glutamate transporter subtypes glutamate transporter 1 and glutamate-aspartate transporter in astrocytes. DL-threo-beta-benzyloxyaspartate (TBOA) and dihydrokainate (DHK) are extensively used as inhibitors of glial glutamate transport activity. Using whole-cell recordings, we characterized the effects of both transporter inhibitors on afferent-evoked astrocyte currents in acute cortical slices of 3-week-old rats. When neuronal afferents were stimulated, passive astrocytes responded by a rapid inward current followed by a persistent tail current. The first current corresponded to a glutamate transporter current. This current was inhibited by both inhibitors and by tetrodotoxin. The tail current is an inward potassium current as it was blocked by barium. Besides inhibiting transporter currents, TBOA strongly enhanced the tail current. This effect was barium-sensitive and might be due to a rise in extracellular potassium level and increased glial potassium uptake. Unlike TBOA, DHK did not enhance the tail current but rather inhibited it. This result suggests that, in addition to inhibiting glutamate transport, DHK prevents astrocyte potassium uptake, possibly by blockade of inward-rectifier channels. This study revealed that, in brain slices, glutamate transporter inhibitors exert complex effects that cannot be attributed solely to glutamate transport inhibition.

Non-invasive quantification of brain glycogen absolute concentration.

The only currently available method to measure brain glycogen in vivo is 13C NMR spectroscopy. Incorporation of 13C-labeled glucose (Glc) is necessary to allow glycogen measurement, but might be affected by turnover changes. Our aim was to measure glycogen absolute concentration in the rat brain by eliminating label turnover as variable. The approach is based on establishing an increased, constant 13C isotopic enrichment (IE). 13C-Glc infusion is then performed at the IE of brain glycogen. As glycogen IE cannot be assessed in vivo, we validated that it can be inferred from that of N-acetyl-aspartate IE in vivo: After [1-13C]-Glc ingestion, glycogen IE was 2.2 +/- 0.1 fold that of N-acetyl-aspartate (n = 11, R(2) = 0.77). After subsequent Glc infusion, glycogen IE equaled brain Glc IE (n = 6, paired t-test, p = 0.37), implying isotopic steady-state achievement and complete turnover of the glycogen molecule. Glycogen concentration measured in vivo by 13C NMR (mean +/- SD: 5.8 +/- 0.7 micromol/g) was in excellent agreement with that in vitro (6.4 +/- 0.6 micromol/g, n = 5). When insulin was administered, the stability of glycogen concentration was analogous to previous biochemical measurements implying that glycogen turnover is activated by insulin. We conclude that the entire glycogen molecule is turned over and that insulin activates glycogen turnover.