152 research outputs found

    Development of 1H-NMR Serum Profiling Methods for High-Throughput Metabolomics

    Get PDF
    El perfilat de sèrum per ressonància magnètica nuclear de protó (1H-RMN) està especialment indicat per a anàlisi a gran escala en estudis epidemiològics, nutricionals o farmacològics. L’espectroscòpia 1H-RMN requereix mínima manipulació de mostra i gràcies a la seva resposta quantitativa permet la comparació directa entre laboratoris. Un perfilat complet de sèrum per 1H-RMN requereix de tres mesures que es corresponen amb tres espècies moleculars diferents: lipoproteïnes, metabòlits de baix pes molecular i lípids. El perfilat de sèrum per 1H-RMN permet obtenir informació de grandària, nombre de partícules i contingut lipídic de les subfraccions lipoproteiques, així com l'abundància d'aminoàcids, productes de la glicòlisi, cossos cetònics, àcids grassos i fosfolípids, entre d'altres. No obstant això, la complexitat espectral afavoreix la inclusió d'errors en l'anàlisi manual de les dades, mentre que les múltiples interaccions moleculars en el sèrum comprometen la seva precisió quantitativa. És per tant necessari desenvolupar mètodes robustos de perfilat metabòlic per consolidar la 1H-RMN en la pràctica clínica. Per a això, aquesta tesi presenta diverses estratègies metodològiques i computacionals. En el primer treball, es van desenvolupar mètodes de regressió dels lípids del perfil lipídic clàssic, generalitzables a mostres de població sana i amb valors de lípids i lipoproteïnes anormals. Aquests lípids representen els principals indicadors de risc cardiovascular i els objectius terapèutics primaris. En el segon estudi caracteritzem els errors de quantificació en el perfilat 1H-RMN de metabòlits clínicament rellevants, que són deguts a la seva agregació a la proteïna sanguínia. També proposem un mètode que fomenta la competició per l'agregació i que ens permet obtenir quantificacions dels nostres metabòlits properes a les absolutes. Finalment, el tercer treball presenta LipSpin: una eina bioinformàtica de codi obert específicament dissenyada per al perfilat de lípids per 1H-RMN. A més, aquest estudi exposa alguns aspectes metodològics per millorar l'anàlisi de lípids per RMN.El perfilado de suero por resonancia magnética nuclear de protón (1H-RMN) está especialmente indicado para el análisis a gran escala en estudios epidemiológicos, nutricionales o farmacológicos. La espectroscopía 1H-RMN requiere mínima manipulación de muestra y gracias a su respuesta cuantitativa permite la comparación directa entre laboratorios. Un perfilado completo de suero por 1H-RMN requiere de tres mediciones que se corresponden con tres especies moleculares distintas: lipoproteínas, metabolitos de bajo peso molecular y lípidos. El perfilado de suero por 1H-RMN permite obtener información de tamaño, número de partículas y contenido lipídico de las subfracciones lipoproteicas, así como la abundancia de aminoácidos, productos de la glicólisis, cuerpos cetónicos, ácidos grasos y fosfolípidos, entre otros. Sin embargo, la complejidad espectral favorece la inclusión de errores en el análisis manual de los datos, mientras que las múltiples interacciones moleculares en el suero comprometen su precisión cuantitativa. Es por tanto necesario desarrollar métodos robustos de perfilado metabólico para consolidar la 1H-RMN en la práctica clínica. Para ello, esta tesis presenta varias estrategias metodológicas y computacionales. En el primer trabajo, se desarrollaron métodos de regresión de los lípidos del perfil lipídico clásico, generalizables a muestras de población sana y con valores de lípidos y lipoproteínas anormales. Estos lípidos representan los principales indicadores de riesgo cardiovascular y los objetivos terapéuticos primarios. En el segundo estudio caracterizamos los errores de cuantificación en el perfilado 1H-RMN de metabolitos clínicamente relevantes, que son debidos a su agregación a la proteína sanguínea. También proponemos un método que fomenta la competición por la agregación y que nos permite obtener cuantificaciones de nuestros metabolitos cercanas a las absolutas. Por último, el tercer trabajo presenta LipSpin: una herramienta bioinformática de código abierto específicamente diseñada para el perfilado de lípidos por 1H-RMN. Además, este estudio expone algunos aspectos metodológicos para mejorar el análisis de lípidos por RMN.1H-NMR serum profiling is especially suitable for high-throughput epidemiological studies and large-scale nutritional studies and drug monitoring. It requires minimal sample manipulation and its quantitative response allows inter-laboratory comparison. A comprehensive 1H-NMR serum profiling consists of three measurements encoding different molecular species: lipoproteins, low-molecular-weight metabolites and lipids. 1H-NMR serum profiling provides information of size, particle number and lipid content of lipoprotein subclasses, as well as abundance of amino acids, glycolysis-related metabolites, ketone bodies, fatty acids and phospholipids, among others. However, the spectral complexity promotes errors in manual data analysis and the multiple molecular interactions within the sample compromise reliable quantifications. Developing robust methods of metabolite serum profiling is therefore desirable to consolidate high-throughput 1H-NMR in the clinical practice. This thesis presents several methodological and computational strategies to that end. In the first study, we developed generalizable regression methods for lipids in routine clinical practice (known as “lipid panel”), to be applied in healthy population and in a wide spectrum of lipid and lipoprotein abnormalities. These standard lipids are still the main measurements of cardiovascular risk and therapy targets. In the second study, we characterised the quantitative errors introduced by protein binding in 1H-NMR profiling of clinically-relevant LMWM in native serum. Then, we proposed a competitive binding strategy to achieve quantifications closer to absolute concentrations, being fully compatible with high-throughput NMR. Finally, the third study presents LipSpin: an open source bioinformatics tool specifically designed for 1H-NMR profiling of serum lipids. Moreover, some methodological aspects to improve NMR-based serum lipid analysis are discussed

    Dolphin and whale: development, evaluation and application of novel bioinformatics tools for metabolite profiling in high throughput 1H-NMR analysis

    Get PDF
    El perfilat de metabòlits es la tasca més difícil dins l'anàlisi espectral de RMN. El seu objectiu es comprendre els processos biològics que tenen lloc en un moment concret mitjançant la identificació i quantificació dels metabòlits presents en mescles d' RMN complexes. Un espectre de RMN està compost per ressonàncies d'un gran nombre de metabòlits, i aquestes se solen solapar entre elles, canviar de posició depenent del pH de la mostra i poden quedar emmascarades per senyals de macromolècules. Tots aquests problemes compliquen la identificació i quantificació de metabòlits, pel que obtenir un perfil de metabòlits curat en una mostra pot ser un gran repte inclús per usuaris experts. En aquest context, la motivació d'aquesta tesi va néixer amb l'objectiu de donar automatismes i funcions fàcils de fer servir per al perfilat de metabòlits en RMN, millorant la qualitat dels resultats i reduint el temps d'anàlisi. Per fer-ho, es van implementar un conjunt d'algoritmes que van acabar empaquetats en dos programes, Dolphin i Whale.El perfilado de metabolitos es la tarea más difícil dentro del análisis espectral de RMN. Su objetivo es comprender los procesos biológicos que tienen lugar en un momento concreto a través de la identificación y cuantificación de los metabolitos presentes en mezclas de RMN complejas. Un espectro de RMN está compuesto por resonancias de un gran numero de metabolitos, y éstas a menudo se solapan entre ellas, cambian de posición dependiendo del pH de la muestra y pueden quedar enmascaradas por señales de macromoléculas. Todos estos problemas complican la identificación y cuantificación de metabolitos, por lo que obtener un perfilado de metabolitos curado en una muestra puede ser un gran reto incluso para usuarios expertos. En este contexto, la motivación de esta tesis nació con el objetivo de dar automatismos y funciones fáciles de usar para el perfilado de metabolitos en RMN, mejorando la calidad de los resultados y reduciendo el tiempo de análisis. Para hacerlo, se implementaron un conjunto de algoritmos que acabaron empaquetados en dos programas, Dolphin y Whale.Metabolite profiling is the most challenging approach in NMR spectral analysis. It aims to comprehend biological processes occurring in a certain moment through identifying and quantifying metabolites present in complex NMR mixtures. An NMR spectrum is composed by resonances of a huge number of metabolites, and these resonances often overlap between them, shift position depending on the sample pH and can be masked by macromolecules signals. All these drawbacks hinder metabolite identification and quantification, so obtaining a cured metabolite profile of a sample can be a very big issue even for expert users. In this context, the motivation of this thesis was born with the aim to provide automatisms and user-friendly interactive functions for NMR metabolite profiling, improving the quality of the results and reducing the time span of the analysis. To do so, several algorisms were implemented and embedded into two software packages, Dolphin and Whale

    On Nature of the Gradient Echo MR Signal and Its Application to Monitoring Multiple Sclerosis

    Get PDF
    Multiple Sclerosis is a common disease, affecting 2.5 million people world-wide. The clinical course is heterogeneous, ranging from benign disease in which patients live an almost normal life to severe and devastating disease that may shorten life. Despite much research, a fully effective treatment for MS is still unavailable and diagnostic techniques for monitoring MS disease evolution are much needed. As a non-invasive tool, Magnetic resonance imaging: MRI) plays a key role in MS diagnosis. Numerous MRI techniques have been proposed over the years. Among most widely used are conventional T1-weighted: T1W), T2-weighted: T2W) and FLuid Attenuated Inversion Recovery: FLAIR) imaging techniques. However their results do not correlate well with neurological findings. Several advanced MRI techniques are also used as research tools to study MS. Among them are magnetization transfer contrast imaging: MT), MR spectroscopy: MRS), and Diffusion Tensor Imaging: DTI) but they have not penetrated to clinical arena yet. Gradient Echo Plural Contrast Imaging: GEPCI) developed in our laboratory is a post processing technique based on multi-echo gradient echo sequence. It offers basic contrasts such as T1W images and T2* maps obtained from magnitude of GEPCI signal, and frequency maps obtained from GEPCI signal phase. Phase information of Gradient Echo MR signal has recently attracted much attention of the MR community since it manifests superior gray matter/ white matter contrast and sub-cortical contrast, especially at high field: 7 T) MRI. However the nature of this contrast is under intense debates. Our group proposed a theoretical framework - Generalized Lorentzian Approach - which emphasizes that, contrary to a common-sense intuition, phase contrast in brain tissue is not directly proportional to the tissue bulk magnetic susceptibility but is rather determined by the geometrical arrangement of brain tissue components: lipids, proteins, iron, etc.) at the cellular and sub-cellular levels - brain tissue magnetic architecture . In this thesis we have provide first direct prove of this hypothesis by measurement of phase contrast in isolated optic nerve. We have also provided first quantitative measurements of the contribution to phase contrast from the water-macromolecule exchange effect. Based on our measurement in protein solutions, we demonstrated that the magnitude of exchange effect is 1/2 of susceptibility effect and to the opposite sign. GEPCI technique also offers a scoring method for monitoring Multiple Sclerosis based on the quantitative T2* maps generated from magnitude information of gradient echo signal. Herein we demonstrated a strong agreement between GEPCI quantitative scores and traditional lesion load assessment. We also established a correlation between GEPCI scores and clinical tests for MS patients. We showed that this correlation is stronger than that found between traditional lesion load and clinical tests. Such studies will be carried out for longer period and on MS subjects with broader range of disease severity in the future. We have also demonstrated that the magnitude and phase information available from GEPCI experiment can be combined in multiple ways to generate novel contrasts that can help with visualization of neurological brain abnormalities beyond Multiple Sclerosis. In summary, in this study, we 1) propose novel contrasts for GEPCI from its basic images; 2) investigate the biophysical mechanisms behind phase contrast; 3) evaluate the benefits of quantitative T2* map offered by GEPCI in monitoring disease of Multiple Sclerosis by comparing GEPCI results to clinical standard techniques; 4) apply our theoretical framework - Generalized Lorentzian Approach - to better understand phase contrast in MS lesions

    Image Analysis for the Life Sciences - Computer-assisted Tumor Diagnostics and Digital Embryomics

    Get PDF
    Current research in the life sciences involves the analysis of such a huge amount of image data that automatization is required. This thesis presents several ways how pattern recognition techniques may contribute to improved tumor diagnostics and to the elucidation of vertebrate embryonic development. Chapter 1 studies an approach for exploiting spatial context for the improved estimation of metabolite concentrations from magnetic resonance spectroscopy imaging (MRSI) data with the aim of more robust tumor detection, and compares against a novel alternative. Chapter 2 describes a software library for training, testing and validating classification algorithms that estimate tumor probability based on MRSI. It allows flexible adaptation towards changed experimental conditions, classifier comparison and quality control without need for expertise in pattern recognition. Chapter 3 studies several models for learning tumor classifiers that allow for the common unreliability of human segmentations. For the first time, models are used for this task that additionally employ the objective image information. Chapter 4 encompasses two contributions to an image analysis pipeline for automatically reconstructing zebrafish embryonic development based on time-resolved microscopy: Two approaches for nucleus segmentation are experimentally compared, and a procedure for tracking nuclei over time is presented and evaluated

    Nowe metody przetwarzania losowo próbkowanych wielowymiarowych eksperymentów NMR

    Get PDF
    The topic of this dissertation is a new algorithm for processing of sparsely sampled data sets from multidimensional nuclear magnetic resonance (NMR) experiments. NMR remains one of the major experimental technique for studying biological macromolecules. However, increasing size of investigated objects poses a challenge for NMR due to rapidly decreasing sensitivity and increasing signal crowding. The first chapter focuses on recent advances in sensitivity enhancements and summarises a few solutions for resolution of spectral overlap. Subsequently, one describes the crucial and limiting problem of signal sampling in multidimensional NMR, which, up to recently, has impeded the widespread use of high-dimensional NMR methods. Major fast acquisition and non-uniform sampling (NUS) approaches are presented. The particular emphasis was put on detailed discussion of competetive approaches to processing of data from NUS experiments. In chapter 3 the new iterative algorithm is proposed for artefact suppression in high-resolution NMR spectra. The detailed description of its design and implementation is given, and followed by comparison with selected processing methods. The efficacy of the algorithm is demonstrated on model synthetic and experimental data. The last chapter of the thesis shows various applications of the proposed method to existing and new four- and five-dimensional NMR experiments. The algorithm is proven most beneficial in challenging applications including spectra for assignment of sidechain resonances in protein and nucleic acids, NOESY spectra for structural analysis, and cross-correlated relaxation measurements for proteins. // Niniejsza praca jest poświecona nowej metodzie przetwarzania danych pochodzących z oszczędnie próbkowanych wielowymiarowych eksperymentów jądrowego rezonansu magnetycznego (ang. Nuclear Magnetic Resonance, NMR). Technika ta jest, obok krystalografii rentgenowskiej, główną eksperymentalną metodą badawczą pozwalającą na określenie struktury i dynamiki makromolekuł o znaczeniu biologicznym. Jednakże NMR napotyka dwie istotne przeszkody w odniesieniu do dużych biomolekuł, a mianowicie gwałtownie pogarszającą się czułość oraz krytyczne zatłoczenie sygnałów w widmach. W rozdziale pierwszym przedstawiono ostatnie osiagnięcia w poprawie czułości technik NMR oraz rozwiązania służące podniesieniu rozdzielczości widm. Następnie opisano kluczowy problem próbkowania wielowymiarowych sygnałów NMR, który do niedawna uniemożliwiał wykorzystanie pełnego potencjału tych technik do rozdzielenia sygnałów. Omówiono pokrótce współczesne podejścia do szybkiej akwizycji i oszczędnego próbkowania sygnałów NMR (ang. non-uniform sampling, NUS). Szczególny nacisk położono na porównanie i dyskusje wad i zalet stosowanych obecnie metod przetwarzania sygnałów niejednorodnie próbkowanych. W rozdziale 3-cim opisano nowy iteracyjny algorytm oparty o transformacje Fouriera, usuwający artefakty oszczędnego próbkowania w wysokorozdzielczych widmach NMR. Szczegółowo omówiono schemat algorytmu oraz jego programową implementację. Rozdział uzupełnia porównanie wyników algorytmu oraz wybranych metod przetwarzania na wysymulowanych oraz modelowych danych eksperymentalnych. W ostatnim rozdziale pracy zademonstrowano użyteczność nowej metody do literaturowych oraz nowych cztero- i pieciowymiarowych eksperymentów NMR. Wśród proponowanych zastosowań wymienić można widma do przypisania sygnałów w łańcuchach bocznych aminokwasów (w białkach) i pierścieniach rybozy (w kwasach rybonukleinowych), widma NOESY służące określeniu struktury trójwymiarowej biomolekuł, oraz pomiary szybkości relaksacji skorelowanej w łańcuchach głównych białek

    A review of machine learning applications for the proton MR spectroscopy workflow

    Get PDF
    This literature review presents a comprehensive overview of machine learning (ML) applications in proton MR spectroscopy (MRS). As the use of ML techniques in MRS continues to grow, this review aims to provide the MRS community with a structured overview of the state-of-the-art methods. Specifically, we examine and summarize studies published between 2017 and 2023 from major journals in the MR field. We categorize these studies based on a typical MRS workflow, including data acquisition, processing, analysis, and artificial data generation. Our review reveals that ML in MRS is still in its early stages, with a primary focus on processing and analysis techniques, and less attention given to data acquisition. We also found that many studies use similar model architectures, with little comparison to alternative architectures. Additionally, the generation of artificial data is a crucial topic, with no consistent method for its generation. Furthermore, many studies demonstrate that artificial data suffers from generalization issues when tested on in vivo data. We also conclude that risks related to ML models should be addressed, particularly for clinical applications. Therefore, output uncertainty measures and model biases are critical to investigate. Nonetheless, the rapid development of ML in MRS and the promising results from the reviewed studies justify further research in this field.</p

    Quantitative Susceptibility Imaging of Tissue Microstructure Using Ultra-High Field MRI

    Get PDF
    This thesis has used ultra-high field (UHF) magnetic resonance imaging (MRI) to investigate the fundamental relationships between tissue microstructure and such susceptibility-based contrast parameters as the apparent transverse relaxation rate (R2*), the local Larmor frequency shift (LFS) and quantitative volume magnetic susceptibility (QS). The interaction of magnetic fields with biological tissues results in shifts in the LFS which can be used to distinguish underlying cellular architecture. The LFS is also linked to the relaxation properties of tissues in a gradient echo MRI sequence. Equally relevant, histological analysis has identified iron and myelin as two major sources of the LFS. As a result, computation of LFS and the associated volume magnetic susceptibility from MRI phase data may serve as a significant method for in vivo monitoring of changes in iron and myelin associated with normal, healthy aging, as well as neurological disease processes. In this research, the cellular level underpinnings of the R2* and LFS signals were examined in a model rat brain system using 9.4 T MRI. The study was carried out using biophysical modeling and correlation with quantitative histology. For the first time, multiple biophysical modeling schemes were compared in both gray and white matter of excised rat brain tissue. Suprisingly, R2* dependence on tissue orientation has not been fully understood. Accordingly, scaling relations were derived for calculating the reversible, mesoscopic magnetic field component, R2\u27, of the apparent transverse relaxation rate from the orientation dependence in gray and white matter. Our results demonstrate that the orientation dependence of R2* and LFS in both white and cortical gray matter has a sinusoidal dependence on tissue orientation and a linear dependence on the volume fraction of myelin in the tissue. A susceptibility processing pipeline was also developed and applied to the calculation of phase-combined LFS and QS maps. The processing pipeline was subsequently used to monitor myelin and iron changes in multiple sclerosis (MS) patients compared to healthy, age and gender-matched controls. With the use of QS and R2* mapping, evidence of statistically significant increases in iron deposition in sub-cortical gray matter, as well as myelin degeneration along the white matter skeleton, were identified in MS patients. The magnetic susceptibility-based MRI methods were then employed as potential clinical biomarkers for disease severity monitoring of MS. It was demonstrated that the combined use of R2* and QS, obtained from multi-echo gradient echo MRI, could serve as an improved metric for monitoring both gray and white matter changes in early MS

    Development and evaluation of a novel advanced lipoprotein test based on 2d diffusion orderen 1h nmr spectroscopy

    Get PDF
    La determinació de la mida i el nombre de lipoproteïnes utilitzant tests avançats de lipoproteïnes és d'un gran interès clínic ja que el nombre de partícules LDL s'ha posicionat com a millor predictor de risc cardiovascular que el colesterol LDL. Tanmateix, els tests avançats de lipoproteïnes actuals encara no s'han introduït en l'àmbit clínic en gran part per la falta d'una estandarització. En aquesta tesi presentem el test LipoScale, un nou test avançat de lipoproteïnes basat en espectroscopia de RMN de difusió 2D. Amb aquest test es pretén obtenir una millor caracterització de les lipoproteïnes plasmàtiques, tant el seu contingut lipídic com la seva mida i nombre de partícules, de manera que amb ell s'aconsegueixi una millor predicció del risc cardiovascular. Durant el desenvolupament del test s’han estudiat diferents patologies i cohorts dins del marc de les malalties metabòliques (les quals són un factor de risc de les malalties cardiovasculars). Entre les malalties estudiades destaquem la diabetis, la dislipèmia aterògena i la síndrome de l’ovari poliquístic (PCOS). A més, també s’han monitoritzat canvis en el perfil de les lipoproteïnes deguts a intervencions nutricionals i a l’exercici. La principal diferència entre la nostra aproximació i la dels mètodes actuals és que aquests últims utilitzen mètodes de RMN 1D estàndards, mentre que el nostre test està basat en l'ús de gradients de camps magnètic, els quals generen espectres 2D amb els que es pot obtenir informació directa i objectiva de la mida de les partícules lipoproteiques. Aquesta tesi ha generat diferents publicacions científiques així com també s'ha fet la sol•licitud d'una patent europea i s'ha creat una spin-off per comercialitzar el test.La determinación del tamaño y el número de lipoproteínas utilizando tests avanzados de lipoproteínas es de un gran interés clínico ya que el número de partículas LDL se ha posicionado como mejor predictor de riesgo cardiovascular que el colesterol LDL. Sin embargo, los tests avanzados de lipoproteínas actuales aún no se han introducido en el ámbito clínico en gran parte por la falta de una estandarización. En esta tesis presentamos el test LipoScale, un nuevo test avanzado de lipoproteínas basado en espectroscopía de RMN de difusión 2D. Con este test se pretende obtener una mejor caracterización de las lipoproteínas plasmáticas, tanto su contenido lipídico como su tamaño y número de partículas, por lo que con él se consiga una mejor predicción del riesgo cardiovascular. Durante el desarrollo del test se han estudiado diferentes patologías y cohortes dentro del marco de las enfermedades metabólicas (las cuales son un factor de riesgo de las enfermedades cardiovasculares). Entre las enfermedades estudiadas destacamos la diabetes, la dislipemia aterògena y el síndrome del ovario poliquístico (PCOS). Además, también se han monitorizado cambios en el perfil de las lipoproteínas debidos a intervenciones nutricionales y el ejercicio. La principal diferencia entre nuestra aproximación y la de los métodos actuales es que estos últimos utilizan métodos de RMN 1D estándar, mientras que nuestro test está basado en el uso de gradientes de campo magnético, los cuales generan espectros 2D con los que se puede obtener información directa y objetiva del tamaño de las partículas lipoproteicas. Esta tesis a generado diferentes publicaciones científicas así como también se ha hecho la solicitud de una patente europea y se ha creado una spin-off para comercializar el test.Determination of lipoprotein particle size and particle number using advanced lipoprotein analyses is of particular interest since the LDL particle number has been shown to improve cardiovascular disease risk prediction. Advanced lipoprotein tests (ALT), however, are not yet routinely introduced in clinical practice partly due to the lack of standardization. This thesis presents the LipoScale test, a novel advanced lipoprotein test based on 2D diffusion-ordered 1H NMR spectroscopy. This test is to obtain a better characterization of plasma lipoproteins in terms of their lipid content, particle size and particle number that will allow a better assessment of cardiovascular risk. During the development of the test various diseases and cohorts were studied in the context of metabolic diseases (which are a risk factor for cardiovascular disease). Among the diseases studied we highlight diabetes, atherogenic dyslipidemia and polycystic ovary syndrome (PCOS). In addition, changes were also monitored in the lipoprotein profile due to nutritional interventions and exercise. The main difference between our approach and the current NMR methods is that the latter use standard 1D methods, whereas our test is based on the use of magnetic field gradients, which generate the 2D spectra that can be used to get direct and objective information on lipoprotein particle sizes. This thesis generated various scientific publications, includes an application for a European patent and a spin-off has been created to commercialize the test

    Multi-parameter optimisation of quantum optical systems

    Get PDF
    Quantum optical systems are poised to become integral components of technologies of the future. While there is growing commercial interest in these systems---for applications in information processing, secure communication and precision metrology---there remain significant technical challenges to overcome before widespread adoption is possible. In this thesis we consider the general problem of optimising quantum optical systems, with a focus on sensing and information processing applications. We investigate four different classes of system with varying degrees of generality and complexity, and demonstrate four corresponding optimisation techniques. At the most specific end of the spectrum---where behaviour is best understood---we consider the problem of interferometric sensitivity enhancement, specifically in the context of long-baseline gravitational wave detectors. We investigate the use of an auxiliary optomechanical system to generate squeezed light exhibiting frequency-dependent quadrature rotation. Such rotation is necessary to evade the effect of quantum back action and achieve broadband sensitivity beyond the standard quantum limit. We find that a cavity optomechanical system is generally unsuitable for this purpose, since the quadrature rotation occurs in the opposite direction to that required for broadband sensitivity improvement. Next we introduce a general technique to engineer arbitrary optical spring potentials in cavity optomechanical systems. This technique has the potential to optimise many types of sensors relying on the optical spring effect. As an example, we show that this technique could yield an enhancement in sensitivity by a factor of 5 when applied to a certain gravitational sensor based on a levitated cavity mirror. We then consider a particular nanowire-based optomechanical system with potential applications in force sensing. We demonstrate a variety of ways to improve its sensitivity to transient forces. We first apply a non-stationary feedback cooling protocol to the system, and achieve an improvement in peak signal-to-noise ratio by a factor of 3, corresponding to a force resolution of 0.2fN. We then implement two non-stationary estimation schemes, which involve post-processing data taken in the absence of physical feedback cooling, to achieve a comparable enhancement in performance without the need for additional experimental complexity. Finally, to address the most complex of systems, we present a general-purpose machine learning algorithm capable of automatically modelling and optimising arbitrary physical systems without human input. To demonstrate the potential of the algorithm we apply it to a magneto-optical trap used for a quantum memory, and achieve an improvement in optical depth from 138 to 448. The four techniques presented differ significantly in their style and the types of systems to which they are applicable. Successfully harnessing the full range of such optimisation procedures will be vital in unlocking the potential of quantum optical systems in the technologies of the futur

    Application of Singular Spectrum Analysis (SSA), Independent Component Analysis (ICA) and Empirical Mode Decomposition (EMD) for automated solvent suppression and automated baseline and phase correction from multi-dimensional NMR spectra

    Get PDF
    A common problem on protein structure determination by NMR spectroscopy is due to the solvent artifact. Typically, a deuterated solvent is used instead of normal water. However, several experimental methods have been developed to suppress the solvent signal in the case that one has to use a protonated solvent or if the signals of the remaining protons even in a highly deuterated sample are still too strong. For a protein dissolved in 90% H2O / 10% D2O, the concentration of solvent protons is about five orders of magnitude greater than the concentration of the protons of interest in the solute. Therefore, the evaluation of multi-dimensional NMR spectra may be incomplete since certain resonances of interest (e.g. Hα proton resonances) are hidden by the solvent signal and since signal parts of the solvent may be misinterpreted as cross peaks originating from the protein. The experimental solvent suppression procedures typically are not able to recover these significant protein signals. Many post-processing methods have been designed in order to overcome this problem. In this work, several algorithms for the suppression of the water signal have been developed and compared. In particular, it has been shown that the Singular Spectrum Analysis (SSA) can be applied advantageously to remove the solvent artifact from NMR spectra of any dimensionality both digitally and analogically acquired. In particular, the investigated time domain signals (FIDs) are decomposed into water and protein related components by means of an initial embedding of the data in the space of time-delayed coordinates. Eigenvalue decomposition is applied on these data and the component with the highest variance (typically represented by the dominant solvent signal) is neglected before reverting the embedding. Pre-processing (group delay management and signal normalization) and post-processing (inverse normalization, Fourier transformation and phase and baseline corrections) of the NMR data is mandatory in order to obtain a better performance of the suppression. The optimal embedding dimension has been empirically determined in accordance to a specific qualitative and quantitative analysis of the extracted components applied on a back-calculated two-dimensional spectrum of HPr protein from Staphylococcus aureus. Moreover, the investigation of experimental data (three-dimensional 1H13C HCCH-TOCSY spectrum of Trx protein from Plasmodium falciparum and two-dimensional NOESY and TOCSY spectra of HPr protein from Staphylococcus aureus) has revealed the ability of the algorithm to recover resonances hidden underneath the water signal. Pathological diseases and the effects of drugs and lifestyle can be detected from NMR spectroscopy applied on samples containing biofluids (e.g. urine, blood, saliva). The detection of signals of interest in such spectra can be hampered by the solvent as well. The SSA has also been successfully applied to one-dimensional urine, blood and cell spectra. The algorithm for automated solvent suppression has been introduced in the AUREMOL software package (AUREMOL_SSA). It is optionally followed by an automated baseline correction in the frequency domain (AUREMOL_ALS) that can be also used out the former algorithm. The automated recognition of baseline points is differently performed in dependence on the dimensionality of the data. In order to investigate the limitations of the SSA, it has been applied to spectra whose dominant signal is not the solvent (as in case of watergate solvent suppression and in case of back-calculated data not including any experimental water signal) determining the optimal solvent-to-solute ratio. The Independent Component Analysis (ICA) represents a valid alternative for water suppression when the solvent signal is not the dominant one in the spectra (when it is smaller than the half of the strongest solute resonance). In particular, two components are obtained: the solvent and the solute. The ICA needs as input at least as many different spectra (mixtures) as the number of components (source signals), thus the definition of a suitable protocol for generating a dataset of one-dimensional ICA-tailored inputs is straightforward. The ICA has revealed to overcome the SSA limitations and to be able to recover resonances of interest that cannot be detected applying the SSA. The ICA avoids all the pre- and post-processing steps, since it is directly applied in the frequency domain. On the other hand, the selection of the component to be removed is automatically detected in the SSA case (having the highest variance). In the ICA, a visual inspection of the extracted components is still required considering that the output is permutable and scale and sign ambiguities may occur. The Empirical Mode Decomposition (EMD) has revealed to be more suitable for automated phase correction than for solvent suppression purposes. It decomposes the FID into several intrinsic mode functions (IMFs) whose frequency of oscillation decreases from the first to the last ones (that identifies the solvent signal). The automatically identified non-baseline regions in the Fourier transform of the sum of the first IMFs are separately evaluated and genetic algorithms are applied in order to determine the zero- and first-order terms suitable for an optimal phase correction. The SSA and the ALS algorithms have been applied before assigning the two-dimensional NOESY spectrum (with the program KNOWNOE) of the PSCD4-domain of the pleuralin protein in order to increase the number of already existing distance restraints. A new routine to derive 3JHNHα couplings from torsion angles (Karplus relation) and vice versa, has been introduced in the AUREMOL software. Using the newly developed tools a refined three-dimensional structure of the PSCD4-domain could be obtained
    corecore