Simultaneous analysis of plasma and CSF by NMR and hierarchical models fusion.

Because cerebrospinal fluid (CSF) is the biofluid which interacts most closely with the central nervous system, it holds promise as a reporter of neurological disease, for example multiple sclerosis (MScl). To characterize the metabolomics profile of neuroinflammatory aspects of this disease we studied an animal model of MScl-experimental autoimmune/allergic encephalomyelitis (EAE). Because CSF also exchanges metabolites with blood via the blood-brain barrier, malfunctions occurring in the CNS may be reflected in the biochemical composition of blood plasma. The combination of blood plasma and CSF provides more complete information about the disease. Both biofluids can be studied by use of NMR spectroscopy. It is then necessary to perform combined analysis of the two different datasets. Mid-level data fusion was therefore applied to blood plasma and CSF datasets. First, relevant information was extracted from each biofluid dataset by use of linear support vector machine recursive feature elimination. The selected variables from each dataset were concatenated for joint analysis by partial least squares discriminant analysis (PLS-DA). The combined metabolomics information from plasma and CSF enables more efficient and reliable discrimination of the onset of EAE. Second, we introduced hierarchical models fusion, in which previously developed PLS-DA models are hierarchically combined. We show that this approach enables neuroinflamed rats (even on the day of onset) to be distinguished from either healthy or peripherally inflamed rats. Moreover, progression of EAE can be investigated because the model separates the onset and peak of the disease.

Quantitative proteomics and metabolomics analysis of normal human cerebrospinal fluid samples.

The analysis of cerebrospinal fluid (CSF) is used in biomarker discovery studies for various neurodegenerative central nervous system (CNS) disorders. However, little is known about variation of CSF proteins and metabolites between patients without neurological disorders. A baseline for a large number of CSF compounds appears to be lacking. To analyze the variation in CSF protein and metabolite abundances in a number of well-defined individual samples of patients undergoing routine, non-neurological surgical procedures, we determined the variation of various proteins and metabolites by multiple analytical platforms. A total of 126 common proteins were assessed for biological variations between individuals by ESI-Orbitrap. A large spread in inter-individual variation was observed (relative standard deviations [RSDs] ranged from 18 to 148%) for proteins with both high abundance and low abundance. Technical variation was between 15 and 30% for all 126 proteins. Metabolomics analysis was performed by means of GC-MS and nuclear magnetic resonance (NMR) imaging and amino acids were specifically analyzed by LC-MS/MS, resulting in the detection of more than 100 metabolites. The variation in the metabolome appears to be much more limited compared with the proteome: the observed RSDs ranged from 12 to 70%. Technical variation was less than 20% for almost all metabolites. Consequently, an understanding of the biological variation of proteins and metabolites in CSF of neurologically normal individuals appears to be essential for reliable interpretation of biomarker discovery studies for CNS disorders because such results may be influenced by natural inter-individual variations. Therefore, proteins and metabolites with high variation between individuals ought to be assessed with caution as candidate biomarkers because at least part of the difference observed between the diseased individuals and the controls will not be caused by the disease, but rather by the natural biological variation between individuals.

Fusion of metabolomics and proteomics data for biomarkers discovery: case study on the experimental autoimmune encephalomyelitis.

BACKGROUND: Analysis of Cerebrospinal Fluid (CSF) samples holds great promise to diagnose neurological pathologies and gain insight into the molecular background of these pathologies. Proteomics and metabolomics methods provide invaluable information on the biomolecular content of CSF and thereby on the possible status of the central nervous system, including neurological pathologies. The combined information provides a more complete description of CSF content. Extracting the full combined information requires a combined analysis of different datasets i.e. fusion of the data. RESULTS: A novel fusion method is presented and applied to proteomics and metabolomics data from a pre-clinical model of multiple sclerosis: an Experimental Autoimmune Encephalomyelitis (EAE) model in rats. The method follows a mid-level fusion architecture. The relevant information is extracted per platform using extended canonical variates analysis. The results are subsequently merged in order to be analyzed jointly. We find that the combined proteome and metabolome data allow for the efficient and reliable discrimination between healthy, peripherally inflamed rats, and rats at the onset of the EAE. The predicted accuracy reaches 89% on a test set. The important variables (metabolites and proteins) in this model are known to be linked to EAE and/or multiple sclerosis. CONCLUSIONS: Fusion of proteomics and metabolomics data is possible. The main issues of high-dimensionality and missing values are overcome. The outcome leads to higher accuracy in prediction and more exhaustive description of the disease profile. The biological interpretation of the involved variables validates our fusion approach.

NMR and pattern recognition can distinguish neuroinflammation and peripheral inflammation.

Multiple Sclerosis (MScl) is a neurodegenerative disease of the CNS, associated with chronic neuroinflammation. Cerebrospinal fluid (CSF), being in closest interaction with CNS, was used to profile neuroinflammation to discover disease-specific markers. We used the commonly accepted animal model for the neuroinflammatory aspect of MScl: the experimental autoimmune/allergic encephalomyelitis (EAE). A combination of advanced (1)H NMR spectroscopy and pattern recognition methods was used to establish the metabolic profile of CSF of EAE-affected rats (representing neuroinflammation) and of two control groups (healthy and peripherally inflamed) to detect specific markers for early neuroinflammation. We found that the CSF metabolic profile for neuroinflammation is distinct from healthy and peripheral inflammation and characterized by changes in concentrations of metabolites such as creatine, arginine, and lysine. Using these disease-specific markers, we were able to detect early stage neuroinflammation, with high accuracy in a second independent set of animals. This confirms the predictive value of these markers. These findings from the EAE model may help to develop a molecular diagnosis for the early stage MScl in humans.

Application of fluorine NMR for structure identification of steroids.

Fluorinated steroids were examined using 1D and 2D homo- and heteronuclear (19)F NMR, such as (19)F-(1) H and (19)F-(13)C. The utilization of fluorine NMR accounted for spectral simplification and resulted in a straightforward pathway for the determination of structures including the configuration of these compounds; these steroids present an illustrative example for other types of fluorinated compounds, which are increasingly encountered in drug discovery. The potential of (19)F NMR is elaborated on in detail for two compounds containing diastereotopic fluorines with different coupling patterns. The analysis of the coupling patterns and the through-space interactions resulted in the determination of the structure and configuration. Heteronuclear correlation experiments, i.e. (19)F-(1)H HETCOR, (19)F-(13)C HMQC and HMBC, and (19)F-(1)H HOESY, were applied to determine first the relative stereochemistry and then the molecular configuration at C4 and C5 of a steroidal compound bearing a fused three-membered ring with two fluorine substituents. These examples proved (19)F NMR to be a useful addition to the extensively used (1)H and (13)C NMR within structure elucidation and configuration determination of small molecules.

The unstable part of the apical stem of duck hepatitis B virus epsilon shows enhanced base pair opening but not pico- to nanosecond dynamics and is essential for reverse transcriptase binding.

Hepatitis B virus (HBV) replication starts with binding of reverse transcriptase (RT) to the apical stem-loop region of epsilon, a conserved element of the RNA pregenome. For duck HBV, an in vitro replication system has provided molecular details of this interaction. Further insights can be obtained from the structure and dynamics of the duck and human apical stem-loops. Previously, we reported these for the human apical stem-loop. Here, we present the same for the duck counterpart. Unlike its human counterpart, the duck apical stem is unstable in its middle/upper part and contains noncanonical base pairs. This dynamics study is the first of an unstable RNA-DNA stem. Similar to the human stem, the duck apical stem comprises two helical segments with a bend angle of ca. 10 degrees , separated by a nonpaired mobile U residue. It is capped by a well-structured conserved UGUU loop with two residues mobile on the pico- to nanosecond time scale, one of which is involved in RT binding. Remarkably, the unstable middle/upper part of the stem does not show enhanced pico- to nanosecond time scale dynamics. Instead, adenine dispersion relaxation studies indicate enhanced millisecond time scale dynamics involving base pair opening. It can then be concluded that base pair opening is essential for epsilon-RT binding, because stabilization of the stem abolishes binding. We hypothesize that binding occurs by conformational capture of bases in the base pair open state. The unstable secondary structure of the apical stem-loop makes duck epsilon-RT binding unusual in light of recent classifications of RNA target interactions that assume stable secondary structures.

A microfluidic high-resolution NMR flow probe.

A microfluidic high-resolution NMR flow probe based on a novel stripline detector chip is demonstrated. This tool is invaluable for the in situ monitoring of reactions performed in microreactors. As an example, the acetylation of benzyl alcohol with acetyl chloride was monitored. Because of the uncompromised (sub-Hz) resolution, this probe holds great promise for metabolomics studies, as shown by an analysis of 600 nL of human cerebrospinal fluid.

Thermodynamics and NMR studies on Duck, Heron and Human HBV encapsidation signals.

Hepatitis B virus (HBV) replication is initiated by binding of its reverse transcriptase (P) to the apical stem-loop (AL) and primer loop (PL) of epsilon, a highly conserved RNA element at the 5'-end of the RNA pregenome. Mutation studies on duck/heron and human in vitro systems have shown similarities but also differences between their P-epsilon interaction. Here, NMR and UV thermodynamic data on AL (and PL) from these three species are presented. The stabilities of the duck and heron ALs were found to be similar, and much lower than that of human. NMR data show that this low stability stems from an 11-nt internal bulge destabilizing the stem of heron AL. In duck, although structured at low temperature, this region also forms a weak point as its imino resonances broaden to disappearance between 30 and 35 degrees C well below the overall AL melting temperature. Surprisingly, the duck- and heron ALs were both found to be capped by a stable well-structured UGUU tetraloop. All avian ALs are expected to adhere to this because of their conserved sequence. Duck PL is stable and structured and, in view of sequence similarities, the same is expected for heron - and human PL.

Fluorine detected 2D NMR experiments for the practical determination of size and sign of homonuclear F-F and heteronuclear C-F multiple bond J-coupling constants in multiple fluorinated compounds.

The use of fluorine in molecules obtained from chemical synthesis has become increasingly important within the pharmaceutical and agricultural industry. NMR characterization of these compounds is of great value with respect to their structure elucidation, their screening in metabolomics investigations and binding studies. The favorable NMR properties of the fluorine nucleus make NMR with fluorine detection of great value in this respect. A suite of NMR 2D F-F- and F-C-correlation experiments with fluorine detection was applied to the assignment of resonances, (n)J(CF)- and (n)J(FF)-couplings as well as the determination of their size and sign. The utilization of this experiment suite was exemplarily demonstrated for a highly fluorinated vinyl alkyl ether. Especially F-C HSQC and J-scaled F-C HMBC experiments allowed determining the size of the J-couplings of this compound. The relative sign of its homo- and heteronuclear couplings was achieved by different combinations of 2D NMR experiments, including non-selective and F2-selective F-C XLOC, F2-selective F-C HMQC, and F-F COSY. The F2-one/two-site selective F-C XLOC versions were found highly useful, as they led to simplifications of the common E.COSY patterns and resulted in a higher confidence level of the assignment by using selective excitation. The combination of F2-one/two-site selective F-C XLOC experiments with a F2-one-site selective F-C HMQC experiment provided the signs of all (n)J(CF)- and (n)J(FF)-couplings in the vinyl moiety of the test compound. Other combinations of experiments were found useful as well for special purposes when focusing for example on homonuclear couplings a combination of F-F COSY-10 with a F2-one-site selective F-C HMQC could be used. The E.COSY patterns in the spectra demonstrated were analyzed by use of the spin-selective displacement vectors, and in case of the XLOC also by use of the DQ- and ZQ-displacement vectors. The variety of experiments presented shall contribute to facilitate the interpretation of F-C correlations as well as to open alternative pathways for the determination of size and signs of homo- and heteronuclear couplings of multiply fluorinated small molecules.

Interpretation and visualization of non-linear data fusion in kernel space: study on metabolomic characterization of progression of multiple sclerosis.

BACKGROUND: In the last decade data fusion has become widespread in the field of metabolomics. Linear data fusion is performed most commonly. However, many data display non-linear parameter dependences. The linear methods are bound to fail in such situations. We used proton Nuclear Magnetic Resonance and Gas Chromatography-Mass Spectrometry, two well established techniques, to generate metabolic profiles of Cerebrospinal fluid of Multiple Sclerosis (MScl) individuals. These datasets represent non-linearly separable groups. Thus, to extract relevant information and to combine them a special framework for data fusion is required. METHODOLOGY: The main aim is to demonstrate a novel approach for data fusion for classification; the approach is applied to metabolomics datasets coming from patients suffering from MScl at a different stage of the disease. The approach involves data fusion in kernel space and consists of four main steps. The first one is to extract the significant information per data source using Support Vector Machine Recursive Feature Elimination. This method allows one to select a set of relevant variables. In the next step the optimized kernel matrices are merged by linear combination. In step 3 the merged datasets are analyzed with a classification technique, namely Kernel Partial Least Square Discriminant Analysis. In the final step, the variables in kernel space are visualized and their significance established. CONCLUSIONS: We find that fusion in kernel space allows for efficient and reliable discrimination of classes (MScl and early stage). This data fusion approach achieves better class prediction accuracy than analysis of individual datasets and the commonly used mid-level fusion. The prediction accuracy on an independent test set (8 samples) reaches 100%. Additionally, the classification model obtained on fused kernels is simpler in terms of complexity, i.e. just one latent variable was sufficient. Finally, visualization of variables importance in kernel space was achieved.

Determination of size and sign of hetero-nuclear coupling constants from 2D 19F-13C correlation spectra.

Fluorinated organic compounds have become increasingly important within the polymer and the pharmaceutical industry as well as for clinical applications. For the structural elucidation of such compounds, NMR experiments with fluorine detection are of great value due to the favorable NMR properties of the fluorine nucleus. For the investigation of three fluorinated compounds, triple resonance 2D HSQC and HMBC experiments were adopted to fluorine detection with carbon and/or proton decoupling to yield F-C, F-C{H}, F-C{C(acq)} and F-C{H,C(acq)} variants. Analysis of E.COSY type cross-peak patterns in the F-C correlation spectra led, apart from the chemical shift assignments, to determination of size and signs of the J(CH), J(CF), and J(HF) coupling constants. In addition, the fully coupled F-C HMQC spectrum of steroid 1 was interpreted in terms of E.COSY type patterns. This example shows how coupling constants due to different nuclei can be determined together with their relative signs from a single spectrum. The analysis of cross-peak patterns, as presented here, not only provides relatively straightforward routes to the determination of size and sign of hetero-nuclear J-couplings in fluorinated compounds, it also provides new and easy ways for the determination of residual dipolar couplings and thus for structure elucidation. The examples and results presented in this study may contribute to a better interpretation and understanding of various F-C correlation experiments and thereby stimulate their utilization.

Determination and identification of estrogenic compounds generated with biosynthetic enzymes using hyphenated screening assays, high resolution mass spectrometry and off-line NMR.

This paper describes the determination and identification of active and inactive estrogenic compounds produced by biosynthetic methods. A hyphenated screening assay towards the human estrogen receptor ligand binding domain (hER)alpha and hERbeta integrating target-ligand interactions and liquid chromatography-high resolution mass spectrometry was used. With this approach, information on both biologic activity and structure identity of compounds produced by bacterial mutants of cytochrome P450s was obtained in parallel. Initial structure identification was achieved by high resolution MS/MS, while for full structure determination, P450 incubations were scaled up and the produced entities were purified using preparative liquid chromatography with automated fraction collection. NMR spectroscopy was performed on all fractions for 3D structure analysis; this included 1D-(1)H, 2D-COSY, 2D-NOESY, and (1)H-(13)C-HSQC experiments. This multidimensional screening approach enabled the detection of low abundant biotransformation products which were not suitable for detection in either one of its single components. In total, the analytical scale biosynthesis produced over 85 compounds from 6 different starting templates. Inter- and intra-day variation of the biochemical signals in the dual receptor affinity detection system was less than 5%. The multi-target screening approach combined with full structure characterization based on high resolution MS(/MS) and NMR spectroscopy demonstrated in this paper can generally be applied to e.g. metabolism studies and compound-library screening.

Low temperature in situ UV irradiation of [(eta 5-C5H5)Co(C2H4)2] in the NMR probe: formation of Co(III) silyl hydride complexes.

Low temperature in situ UV irradiation of [(eta(5)-C(5)H(5))Co(C(2)H(4))(2)] in the presence of silanes enables the characterisation of unstable fluxional Co(III) silyl hydride complexes [(eta(5)-C(5)H(5))Co(SiR(3))(H)(C(2)H(4))] (SiR(3) = SiEt(3), SiMe(3) or SiHEt(2)) by NMR spectroscopy; the reaction of [Co(eta(5)-C(5)H(5))(C(2)H(4))(2)] with HSiR(3) proceeds thermally to reach an equilibrium when SiR(3) = Si(OMe)(3) or SiClMePh.

Photochemical isomerization of N-heterocyclic carbene ruthenium hydride complexes: in situ photolysis, parahydrogen, and computational studies.

Low-temperature UV irradiation of the N-heterocyclic carbene complex Ru(IEt2Me2)(PPh3)2(CO)H2 (IEt2Me2 = 1,3-bis(ethyl)-4,5-dimethylimidazol-2-ylidene) leads to a remarkable photoisomerization reaction. By combining in situ photolysis and parahydrogen experiments to characterize the ultimate photoproducts and DFT calculations to interrogate the structures of the key 16-electron intermediates, the importance of both PPh3 and H2 loss pathways has been established.

Photochemical reactions of (eta 5-cyclopentadienyl)bis(t-butylacrylate) rhodium with silanes: dynamics of isomer interconversion via Rh(eta 2-silane) species.

The compound CpRh(C(2)H(3)CO(2)(t)Bu)(2) 1 has been synthesised as a mixture of two pairs of interconverting isomers which differ in the relative orientations of the alkene substituents. The four isomers have been fully characterised by NMR spectroscopy. When complex 1 is photolysed in the presence of a silane, HSiR(2)R'R(2)R'= Et(3), Me(3), HEt(2), (OMe)(3) and Me(2)Cl] the corresponding Si-H oxidative addition products CpRh(SiR(2)R')(H)(C(2)H(3)CO(2)(t)Bu) and CpRh(H)(2)(SiR(2)R')(2) are formed. The Rh(III) complexes CpRh(SiR(2)R')(H)(C(2)H(3)CO(2)(t)Bu) exist in two isomeric forms of comparable energy which interconvert in an intramolecular process that does not involve a reversible [1,3] hydride or [1,3] silyl migration. The hydride (1)H NMR resonances for these species consequently broaden before coalescing into a single peak. For R(2)R'= Et(3), the activation parameters for interchange from the major to minor isomer were Delta H++= 60.2 +/- 2 kJ mol(-1) and Delta S++= 8 +/- 9 J mol(-1) K(-1), while for R(2)R'= Me(3) and Et(2)H, Delta H++= 61.5 +/- 1 kJ mol(-1), Delta S++= 6 +/- 5 J mol(-1) K(-1), and Delta H++= 61.8 +/- 3 kJ mol(-1), Delta S++= 12 +/- 9 J mol(-1) K(-1) respectively for conversion from the major isomer to the minor. For these complexes an eta(2)-Rh-H-Si transition state or intermediate is consistent with the evidence. When R(2)R'=(OMe)(3) and Me(2)Cl the change in appearance of the hydride resonances is more complex, with the activation parameters for interchange from the major to minor isomer for the former species being Delta H++= 78.3 +/- 2 kJ mol(-1) and Delta S++= 30 +/- 7 J mol(-1) K(-1) while for Me(2)Cl the barrier proved too high to measure before decomposition occurred. The complex spectral changes could be simulated when a discrete eta(2)-Rh-H-Si intermediate was involved in the isomer interconversion process and hence silane rotation in all these systems is proposed to involve two isomers of CpRh(eta(2)-HSiR(2)R')(C(2)H(3)CO(2)(t)Bu).