Operations research: from computational biology to sensor network

In this dissertation we discuss the deployment of combinatorial optimization methods for modeling and solve real life problemS, with a particular emphasis to two biological problems arising from a common scenario: the reconstruction of the three-dimensional shape of a biological molecule from Nuclear Magnetic Resonance (NMR) data. The fi rst topic is the 3D assignment pathway problem (APP) for a RNA molecule. We prove that APP is NP-hard, and show a formulation of it based on edge-colored graphs. Taking into account that interactions between consecutive nuclei in the NMR spectrum are diff erent according to the type of residue along the RNA chain, each color in the graph represents a type of interaction. Thus, we can represent the sequence of interactions as the problem of fi nding a longest (hamiltonian) path whose edges follow a given order of colors (i.e., the orderly colored longest path). We introduce three alternative IP formulations of APP obtained with a max flow problem on a directed graph with packing constraints over the partitions, which have been compared among themselves. Since the last two models work on cyclic graphs, for them we proposed an algorithm based on the solution of their relaxation combined with the separation of cycle inequalities in a Branch & Cut scheme. The second topic is the discretizable distance geometry problem (DDGP), which is a formulation on discrete search space of the well-known distance geometry problem (DGP). The DGP consists in seeking the embedding in the space of a undirected graph, given a set of Euclidean distances between certain pairs of vertices. DGP has two important applications: (i) fi nding the three dimensional conformation of a molecule from a subset of interatomic distances, called Molecular Distance Geometry Problem, and (ii) the Sensor Network Localization Problem. We describe a Branch & Prune (BP) algorithm tailored for this problem, and two versions of it solving the DDGP both in protein modeling and in sensor networks localization frameworks. BP is an exact and exhaustive combinatorial algorithm that examines all the valid embeddings of a given weighted graph G=(V,E,d), under the hypothesis of existence of a given order on V. By comparing the two version of BP to well-known algorithms we are able to prove the e fficiency of BP in both contexts, provided that the order imposed on V is maintained

Compressed Sensing Accelerated Magnetic Resonance Spectroscopic Imaging

abstract: Magnetic resonance spectroscopic imaging (MRSI) is a valuable technique for assessing the in vivo spatial profiles of metabolites like N-acetylaspartate (NAA), creatine, choline, and lactate. Changes in metabolite concentrations can help identify tissue heterogeneity, providing prognostic and diagnostic information to the clinician. The increased uptake of glucose by solid tumors as compared to normal tissues and its conversion to lactate can be exploited for tumor diagnostics, anti-cancer therapy, and in the detection of metastasis. Lactate levels in cancer cells are suggestive of altered metabolism, tumor recurrence, and poor outcome. A dedicated technique like MRSI could contribute to an improved assessment of metabolic abnormalities in the clinical setting, and introduce the possibility of employing non-invasive lactate imaging as a powerful prognostic marker. However, the long acquisition time in MRSI is a deterrent to its inclusion in clinical protocols due to associated costs, patient discomfort (especially in pediatric patients under anesthesia), and higher susceptibility to motion artifacts. Acceleration strategies like compressed sensing (CS) permit faithful reconstructions even when the k-space is undersampled well below the Nyquist limit. CS is apt for MRSI as spectroscopic data are inherently sparse in multiple dimensions of space and frequency in an appropriate transform domain, for e.g. the wavelet domain. The objective of this research was three-fold: firstly on the preclinical front, to prospectively speed-up spectrally-edited MRSI using CS for rapid mapping of lactate and capture associated changes in response to therapy. Secondly, to retrospectively evaluate CS-MRSI in pediatric patients scanned for various brain-related concerns. Thirdly, to implement prospective CS-MRSI acquisitions on a clinical magnetic resonance imaging (MRI) scanner for fast spectroscopic imaging studies. Both phantom and in vivo results demonstrated a reduction in the scan time by up to 80%, with the accelerated CS-MRSI reconstructions maintaining high spectral fidelity and statistically insignificant errors as compared to the fully sampled reference dataset. Optimization of CS parameters involved identifying an optimal sampling mask for CS-MRSI at each acceleration factor. It is envisioned that time-efficient MRSI realized with optimized CS acceleration would facilitate the clinical acceptance of routine MRSI exams for a quantitative mapping of important biomarkers.Dissertation/ThesisDoctoral Dissertation Bioengineering 201

Structures of fibrils formed by α-synuclein hereditary disease mutant H50Q reveal new polymorphs.

Deposits of amyloid fibrils of α-synuclein are the histological hallmarks of Parkinson's disease, dementia with Lewy bodies and multiple system atrophy, with hereditary mutations in α-synuclein linked to the first two of these conditions. Seeing the changes to the structures of amyloid fibrils bearing these mutations may help to understand these diseases. To this end, we determined the cryo-EM structures of α-synuclein fibrils containing the H50Q hereditary mutation. We find that the H50Q mutation results in two previously unobserved polymorphs of α-synuclein: narrow and wide fibrils, formed from either one or two protofilaments, respectively. These structures recapitulate conserved features of the wild-type fold but reveal new structural elements, including a previously unobserved hydrogen-bond network and surprising new protofilament arrangements. The structures of the H50Q polymorphs help to rationalize the faster aggregation kinetics, higher seeding capacity in biosensor cells and greater cytotoxicity that we observe for H50Q compared to wild-type α-synuclein

Bacteriophages and their structural organisation

Viruses are extremely small infectious particles that are not visible in a light microscope, and are able to pass through fine porcelain filters. They exist in a huge variety of forms and infect practically all living systems: animals, plants, insects and bacteria. All viruses have a genome, typically only one type of nucleic acid, but it could be one or several molecules of DNA or RNA, which is surrounded by a protective stable coat (capsid) and sometimes by additional layers which may be very complex and contain carbohydrates, lipids, and additional proteins. The viruses that have only a protein coat are named “naked”, or non- enveloped viruses. Many viruses have an envelope (enveloped viruses) that wraps around the protein capsid. This envelope is formed from a lipid membrane of the host cell during the release of a virus out of the cell

Updates in metabolomics tools and resources: 2014-2015

Data processing and interpretation represent the most challenging and time-consuming steps in high-throughput metabolomic experiments, regardless of the analytical platforms (MS or NMR spectroscopy based) used for data acquisition. Improved machinery in metabolomics generates increasingly complex datasets that create the need for more and better processing and analysis software and in silico approaches to understand the resulting data. However, a comprehensive source of information describing the utility of the most recently developed and released metabolomics resources—in the form of tools, software, and databases—is currently lacking. Thus, here we provide an overview of freely-available, and open-source, tools, algorithms, and frameworks to make both upcoming and established metabolomics researchers aware of the recent developments in an attempt to advance and facilitate data processing workflows in their metabolomics research. The major topics include tools and researches for data processing, data annotation, and data visualization in MS and NMR-based metabolomics. Most in this review described tools are dedicated to untargeted metabolomics workflows; however, some more specialist tools are described as well. All tools and resources described including their analytical and computational platform dependencies are summarized in an overview Table

Altered Neurocircuitry in the Dopamine Transporter Knockout Mouse Brain

The plasma membrane transporters for the monoamine neurotransmitters dopamine, serotonin, and norepinephrine modulate the dynamics of these monoamine neurotransmitters. Thus, activity of these transporters has significant consequences for monoamine activity throughout the brain and for a number of neurological and psychiatric disorders. Gene knockout (KO) mice that reduce or eliminate expression of each of these monoamine transporters have provided a wealth of new information about the function of these proteins at molecular, physiological and behavioral levels. In the present work we use the unique properties of magnetic resonance imaging (MRI) to probe the effects of altered dopaminergic dynamics on meso-scale neuronal circuitry and overall brain morphology, since changes at these levels of organization might help to account for some of the extensive pharmacological and behavioral differences observed in dopamine transporter (DAT) KO mice. Despite the smaller size of these animals, voxel-wise statistical comparison of high resolution structural MR images indicated little morphological change as a consequence of DAT KO. Likewise, proton magnetic resonance spectra recorded in the striatum indicated no significant changes in detectable metabolite concentrations between DAT KO and wild-type (WT) mice. In contrast, alterations in the circuitry from the prefrontal cortex to the mesocortical limbic system, an important brain component intimately tied to function of mesolimbic/mesocortical dopamine reward pathways, were revealed by manganese-enhanced MRI (MEMRI). Analysis of co-registered MEMRI images taken over the 26 hours after introduction of Mn^(2+) into the prefrontal cortex indicated that DAT KO mice have a truncated Mn^(2+) distribution within this circuitry with little accumulation beyond the thalamus or contralateral to the injection site. By contrast, WT littermates exhibit Mn^(2+) transport into more posterior midbrain nuclei and contralateral mesolimbic structures at 26 hr post-injection. Thus, DAT KO mice appear, at this level of anatomic resolution, to have preserved cortico-striatal-thalamic connectivity but diminished robustness of reward-modulating circuitry distal to the thalamus. This is in contradistinction to the state of this circuitry in serotonin transporter KO mice where we observed more robust connectivity in more posterior brain regions using methods identical to those employed here

Protease Cleavage Leads to Formation of Mature Trimer Interface in HIV-1 Capsid

During retrovirus particle maturation, the assembled Gag polyprotein is cleaved by the viral protease into matrix (MA), capsid (CA), and nucleocapsid (NC) proteins. To form the mature viral capsid, CA rearranges, resulting in a lattice composed of hexameric and pentameric CA units. Recent structural studies of assembled HIV-1 CA revealed several inter-subunit interfaces in the capsid lattice, including a three-fold interhexamer interface that is critical for proper capsid stability. Although a general architecture of immature particles has been provided by cryo-electron tomographic studies, the structural details of the immature particle and the maturation pathway remain unknown. Here, we used cryo-electron microscopy (cryoEM) to determine the structure of tubular assemblies of the HIV-1 CA-SP1-NC protein. Relative to the mature assembled CA structure, we observed a marked conformational difference in the position of the CA-CTD relative to the NTD in the CA-SP1-NC assembly, involving the flexible hinge connecting the two domains. This difference was verified via engineered disulfide crosslinking, revealing that inter-hexamer contacts, in particular those at the pseudo three-fold axis, are altered in the CA-SP1-NC assemblies compared to the CA assemblies. Results from crosslinking analyses of mature and immature HIV-1 particles containing the same Cys substitutions in the Gag protein are consistent with these findings. We further show that cleavage of preassembled CA-SP1-NC by HIV-1 protease in vitro leads to release of SP1 and NC without disassembly of the lattice. Collectively, our results indicate that the proteolytic cleavage of Gag leads to a structural reorganization of the polypeptide and creates the three-fold interhexamer interface, important for the formation of infectious HIV-1 particles. © 2012 Meng et al

23Na multi-quantum coherences: from cellular spectroscopy to clinical imaging development

Non-invasive biological tissue information is vital for medical diagnostics and treat- ment monitoring. While standard 1H magnetic resonance imaging (MRI) methods show detailed morphological information, 23Na MRI provides additional biochemical information about the tissue. Sodium nuclei have spin 3/2 and, therefore, can exhibit higher quantum coherence signals. Multi-quantum (MQ) imaging offers additional information compared to standard SQ sodium, which focuses on tissue sodium concentration (TSC), e.g., it is hard to discern edema and tumor regions, both exhibiting higher TSC. 23Na triple-quantum (TQ) signals are of high interest to probe the molecular environment in tissues and alleviate this problem. The focus of this thesis was on the 23Na TQ signal from preclinical investigations on cells via a simulation study and, finally, a transfer of preclinical and simulation find- ings into the development of an optimal clinical MQ imaging sequence, CRISTINA. First, the physiological importance was studied at ultra-high field, 9.4 T, to gain in- sight into TQ signal changes under different cellular conditions. An MR-compatible bioreactor setup allowed for finely tunable TQ signal monitoring of cell lines human liver cells (HEP G2) and neonatal cardiomyocytes of mice. Both cell lines showed a TQ signal in vital state, under standard perfusion [0.26%,15σ], normalized to the SQ signal. Hypertrophy was simulated with oxygen and nutrient stop and resulted in a TQ signal to 56% of the initial value [0.15%,24σ]. Re-perfusion resulted in a come back of the TQ signal to 92% of the initial value. Reference measurements without cells as well as dead cells showed a TQ signal of [0.06%,1σ and 0.016%,1σ]. Further, the long-standing debate of TQ signal connection to intracellular space was investigated based on liposomal cell-phantoms and it was shown that the TQ signal in liposomes was related to the interaction of the sodium ions with the double lipid membrane, which is constituted of negatively charged fatty acids. A single-voxel localization technique was developed on the preclinical system as the first step in the direction of a clinical sequence and tested on phantoms and in-vivo rat. Second, simulation of different phase-cycle schemes of the standard three-pulses coherence transfer technique was performed. A unified framework was developed to compare and find an optimum. Destructive effects of B0 inhomogeneity were investigated and verified for Ω(Hz) = (kπ + ξ)/(τ1), k ∈ Z. Stimulated echo signal stood as further potential biases but resulted in a continuous offset after Fourier Transformation. Third, knowledge from part I and II was transferred to develop an efficient clinical vii MQ imaging method: CRISTINA, a 2D Cartesian MQ, multi-echo imaging sequence for clinical use. A Multi-parameter fit routine provided T2 relaxations maps and ratio of TQ to SQ signals which could be of interest to monitor pathologies in future. CRISTINA was tested and optimized on phantoms and in vivo on 5 healthy brain volunteers. A linear relationship was found for the ratio TQ over SQ signal against agar concentrations (R2 =0.87, p-value = 0.0007) as well as for the SQ signal against TSC (R2 =0.75, p-value = 0.006)

Imaging and 3D reconstruction of membrane protein complexes by cryo-electron microscopy and single particle analysis

Cryo-electron microscopy (cryo-EM) in combination with single particle image processing and volume reconstruction is a powerful technology to obtain medium-resolution structures of large protein complexes, which are extremely difficult to crystallize and not amenable to NMR studies due to size limitation. Depending on the stability and stiffness as well as on the symmetry of the complex, three-dimensional reconstructions at a resolution of 10-30 ˚ can be achieved. In this range of resolution, we may not be able to answer A chemical questions at the level of atomic interactions, but we can gain detailed insight into the macromolecular architecture of large multi-subunit complexes and their mechanisms of action. In this thesis, several prevalently large membrane protein complexes of great physiological importance were examined by various electron microscopy techniques and single particle image analysis. The core part of my work consists in the imaging of a mammalian V-ATPase, frozen-hydrated in amorphous ice and of the completion of the first volume reconstruction of this type of enzyme, derived from cryo-EM images. This ubiquitous rotary motor is essential in every eukaryotic cell and is of high medical importance due to its implication in various diseases such as osteoporosis, skeletal cancer and kidney disorders. My contribution to the second and third paper concerns the volume reconstruction of two bacterial outer membrane pore complexes from cryo-EM images recorded by my colleague Mohamed Chami. PulD from Klebsiella oxytoca constitutes a massive translocating pore capable of transporting a fully folded cell surface protein PulA through the membrane. It is part of the Type II secretion system, which is common for Gram-negative bacteria. The second volume regards ClyA, a pore-forming heamolytic toxin of virulent Escherichia coli and Salmonella enterica strains that kill target cells by inserting pores into their membranes. To the last two papers, I contributed with cryo-negative stain imaging of the cell division protein DivIVA from Bacillus subtilis and with image processing of the micrographs displaying the siderophore receptor FrpB from Neisseria meningitidis