Modeling Techniques for the High-Resolution Interpretation of Cryo-Electron Microscopy Reconstructions

Essential biological processes are governed by organized, dynamic interactions between multiple biomolecular systems. Complexes are thus formed to enable the biological function and get dissembled as the process is completed. Examples of such processes include the translation of the messenger RNA into protein by the ribosome, the folding of proteins by chaperonins or the entry of viruses in host cells. Understanding these fundamental processes by characterizing the molecular mechanisms that enable then, would allow the (better) design of therapies and drugs. Such molecular mechanisms may be revealed trough the structural elucidation of the biomolecular assemblies at the core of these processes. Various experimental techniques may be applied to investigate the molecular architecture of biomolecular assemblies. High-resolution techniques, such as X-ray crystallography, may solve the atomic structure of the system, but are typically constrained to biomolecules of reduced flexibility and dimensions. In particular, X-ray crystallography requires the sample to form a three dimensional (3D) crystal lattice which is technically di‑cult, if not impossible, to obtain, especially for large, dynamic systems. Often these techniques solve the structure of the different constituent components within the assembly, but encounter difficulties when investigating the entire system. On the other hand, imaging techniques, such as cryo-electron microscopy (cryo-EM), are able to depict large systems in near-native environment, without requiring the formation of crystals. The structures solved by cryo-EM cover a wide range of resolutions, from very low level of detail where only the overall shape of the system is visible, to high-resolution that approach, but not yet reach, atomic level of detail. In this dissertation, several modeling methods are introduced to either integrate cryo-EM datasets with structural data from X-ray crystallography, or to directly interpret the cryo-EM reconstruction. Such computational techniques were developed with the goal of creating an atomic model for the cryo-EM data. The low-resolution reconstructions lack the level of detail to permit a direct atomic interpretation, i.e. one cannot reliably locate the atoms or amino-acid residues within the structure obtained by cryo-EM. Thereby one needs to consider additional information, for example, structural data from other sources such as X-ray crystallography, in order to enable such a high-resolution interpretation. Modeling techniques are thus developed to integrate the structural data from the different biophysical sources, examples including the work described in the manuscript I and II of this dissertation. At intermediate and high-resolution, cryo-EM reconstructions depict consistent 3D folds such as tubular features which in general correspond to alpha-helices. Such features can be annotated and later on used to build the atomic model of the system, see manuscript III as alternative. Three manuscripts are presented as part of the PhD dissertation, each introducing a computational technique that facilitates the interpretation of cryo-EM reconstructions. The first manuscript is an application paper that describes a heuristics to generate the atomic model for the protein envelope of the Rift Valley fever virus. The second manuscript introduces the evolutionary tabu search strategies to enable the integration of multiple component atomic structures with the cryo-EM map of their assembly. Finally, the third manuscript develops further the latter technique and apply it to annotate consistent 3D patterns in intermediate-resolution cryo-EM reconstructions. The first manuscript, titled An assembly model for Rift Valley fever virus, was submitted for publication in the Journal of Molecular Biology. The cryo-EM structure of the Rift Valley fever virus was previously solved at 27Å-resolution by Dr. Freiberg and collaborators. Such reconstruction shows the overall shape of the virus envelope, yet the reduced level of detail prevents the direct atomic interpretation. High-resolution structures are not yet available for the entire virus nor for the two different component glycoproteins that form its envelope. However, homology models may be generated for these glycoproteins based on similar structures that are available at atomic resolutions. The manuscript presents the steps required to identify an atomic model of the entire virus envelope, based on the low-resolution cryo-EM map of the envelope and the homology models of the two glycoproteins. Starting with the results of the exhaustive search to place the two glycoproteins, the model is built iterative by running multiple multi-body refinements to hierarchically generate models for the different regions of the envelope. The generated atomic model is supported by prior knowledge regarding virus biology and contains valuable information about the molecular architecture of the system. It provides the basis for further investigations seeking to reveal different processes in which the virus is involved such as assembly or fusion. The second manuscript was recently published in the of Journal of Structural Biology (doi:10.1016/j.jsb.2009.12.028) under the title Evolutionary tabu search strategies for the simultaneous registration of multiple atomic structures in cryo-EM reconstructions. This manuscript introduces the evolutionary tabu search strategies applied to enable a multi-body registration. This technique is a hybrid approach that combines a genetic algorithm with a tabu search strategy to promote the proper exploration of the high-dimensional search space. Similar to the Rift Valley fever virus, it is common that the structure of a large multi-component assembly is available at low-resolution from cryo-EM, while high-resolution structures are solved for the different components but lack for the entire system. Evolutionary tabu search strategies enable the building of an atomic model for the entire system by considering simultaneously the different components. Such registration indirectly introduces spatial constrains as all components need to be placed within the assembly, enabling the proper docked in the low-resolution map of the entire assembly. Along with the method description, the manuscript covers the validation, presenting the benefit of the technique in both synthetic and experimental test cases. Such approach successfully docked multiple components up to resolutions of 40Å. The third manuscript is entitled Evolutionary Bidirectional Expansion for the Annotation of Alpha Helices in Electron Cryo-Microscopy Reconstructions and was submitted for publication in the Journal of Structural Biology. The modeling approach described in this manuscript applies the evolutionary tabu search strategies in combination with the bidirectional expansion to annotate secondary structure elements in intermediate resolution cryo-EM reconstructions. In particular, secondary structure elements such as alpha helices show consistent patterns in cryo-EM data, and are visible as rod-like patterns of high density. The evolutionary tabu search strategy is applied to identify the placement of the different alpha helices, while the bidirectional expansion characterizes their length and curvature. The manuscript presents the validation of the approach at resolutions ranging between 6 and 14Å, a level of detail where alpha helices are visible. Up to resolution of 12 Å, the method measures sensitivities between 70-100% as estimated in experimental test cases, i.e. 70-100% of the alpha-helices were correctly predicted in an automatic manner in the experimental data. The three manuscripts presented in this PhD dissertation cover different computation methods for the integration and interpretation of cryo-EM reconstructions. The methods were developed in the molecular modeling software Sculptor (http://sculptor.biomachina.org) and are available for the scientific community interested in the multi-resolution modeling of cryo-EM data. The work spans a wide range of resolution covering multi-body refinement and registration at low-resolution along with annotation of consistent patterns at high-resolution. Such methods are essential for the modeling of cryo-EM data, and may be applied in other fields where similar spatial problems are encountered, such as medical imaging.

Multiple GABA receptors in rabbit retina

The task of encoding and processing complex sensory input requires many types of transsynaptic signals. This requirement is served in part by an extensive group of neurotransmitter substances which may include thirty or more different compounds. At the next level of information processing, the existence of multiple receptors for a given neurotransmitter appears to be a widely used mechanism to generate multiple responses to a given first messenger (Snyder and Goodman, 1980). Despite the wealth of published data on GABA receptors, the existence of more than one GABA receptor was in doubt until the mid 1980's. Presently there is still disagreement on the number of types of GABA receptors, estimates for which range from two to four (DeFeudis, 1983; Johnston, 1985). Part of the problem in evaluating data concerning multiple receptor types is the lack of information on the number of gene products and their subsequent supramolecular organization in different neurons. In order to evaluate the question concerning the diversity of GABA receptors in the nervous system, we must rely on indirect information derived from a wide variety of experimental techniques. These include pharmacological binding studies to membrane fractions, electrophysiological studies, localization studies, purification studies, and functional assays. Almost all parts of the central and peripheral nervous system use GABA as a neurotransmitter, and these experimental techniques have therefore been applied to many different parts of the nervous system for the analysis of GABA receptor characteristics. We are left with a large amount of data from a wide variety of techniques derived from many parts of the nervous system. When this project was initiated in 1983, there were only a handful of pharmacological tools to assess the question of multiple GABA receptors. The approach adopted was to focus on a single model system, using a variety of experimental techniques, in order to evaluate the existence of multiple forms of GABA receptors. Using the in vitro rabbit retina, a combination of pharmacological binding studies, functional release studies and partial purification studies were undertaken to examine the GABA receptor composition of this tissue. Three types of GABA receptors were observed: Al receptors coupled to benzodiazepine and barbiturate modulation, and A2 or uncoupled GABA-A receptors, and GABA-B receptors. These results are evaluated and discussed in light of recent findings by others concerning the number and subtypes of GABA receptors in the nervous system. ^

Effect of Acute Administration of Angiopoietin-1 in Experimental Traumatic Spinal Cord Injury: Magnetic Resonance Imaging and Neurobehavioral Studies

Spinal cord injury (SCI) is a devastating condition that affects people in the prime of their lives. A myriad of vascular events occur after SCI, each of which contributes to the evolving pathology. The primary trauma causes mechanical damage to blood vessels, resulting in hemorrhage. The blood-spinal cord barrier (BSCB), a neurovascular unit that limits passage of most agents from systemic circulation to the central nervous system, breaks down, resulting in inflammation, scar formation, and other sequelae. Protracted BSCB disruption may exacerbate cellular injury and hinder neurobehavioral recovery in SCI. In these studies, angiopoietin-1 (Ang1), an agent known to reduce vascular permeability, was hypothesized to attenuate the severity of secondary injuries of SCI. Using longitudinal magnetic resonance imaging (MRI) studies (dynamic contrast-enhanced [DCE]-MRI for quantification of BSCB permeability, highresolution anatomical MRI for calculation of lesion size, and diffusion tensor imaging for assessment of axonal integrity), the acute, subacute, and chronic effects of Ang1 administration after SCI were evaluated. Neurobehavioral assessments were also performed. These non-invasive techniques have applicability to the monitoring of therapies in patients with SCI. In the acute phase of injury, Ang1 was found to reduce BSCB permeability and improve neuromotor recovery. Dynamic contrast-enhanced MRI revealed a persistent compromise of the BSCB up to two months post-injury. In the subacute phase of injury, Ang1’s effect on reducing BSCB permeability was maintained and it was found to transiently reduce axonal integrity. The SCI lesion burden was assessed with an objective method that compared favorably with segmentations from human raters. In the chronic phase of injury, Ang1 resulted in maintained reduction in BSCB permeability, a decrease in lesion size, and improved axonal integrity. Finally, longitudinal correlations among data from the MRI modalities and neurobehavioral assays were evaluated. Locomotor recovery was negatively correlated with lesion size in the Ang1 cohort and positively correlated with diffusion measures in the vehicle cohort. In summary, the results demonstrate a possible role for Ang1 in mitigating the secondary pathologies of SCI during the acute and chronic phases of injury.

{\it In vivo\/} magnetic resonance studies of experimental liver disease: Carbon tetrachloride hepatotoxicity and alcohol-induced fatty liver in rat

Magnetic resonance imaging (MRI) and magnetic resonance spectroscopy (MRS) were used to non-invasively determine if cirrhosis induced by carbon tetrachloride (CCl$\sb4$) and phospholipase-D (PLD) could be distinguished from fatty infiltration in rat. MRS localization and water suppression methods were developed, implemented and evaluated in terms of their application to in vivo proton NMR studies of experimental liver disease. MRS studies were also performed to quantitate fatty infiltration resulting from carbon tetrachloride (CCl$\sb4$) or alcohol (ethanol) administration and the MRS results were confirmed using biochemical total lipid analysis and histology. $\rm T\sb1$ weighted MR images acquired weekly, 48 hours post administration, demonstrated only a slight increase in overall liver intensity with CCl$\sb4$ or alcohol administration, which is consistent with previously reported results. The MR images were able to detect nodules resulting from CCl$\sb4$+PLD induced cirrhosis as hypointense regions, also consistent with previous reports. Localized in vivo water and lipid proton $\rm T\sb1$ relaxation time measurements were performed and demonstrated no statistically significant trends for either agent. In vivo proton spectra were also acquired using stimulated echo techniques to quantitatively follow the changes in liver lipid content. The changes in liver lipid content observed using MRS were verified by total lipid analysis using the Folch technique and histology. The in vivo $\rm T\sb1$ and lipid quantification data str inconsistent with the previous hypothesis that the changes in $\rm T\sb1$ weighted images were the result of increased "free" water content and, therefore, increased water $\rm T\sb1$ relaxation times. These data indicate that the long term changes are more likely the result of changes in lipid content. The data are also shown to agree with the accepted hypothesis that the time course and mechanism of fatty infiltration are different for CCl$\sb4$ and alcohol. The hypothesis that the lipids resulting from either protocol are from the same lipid fraction(s), presumably triglycerides, is also supported. And lastly, on the basis of MR images and quantitative MRS lipid information, it was shown that cirrhosis could be distinguished from fatty infiltration. ^

Experimental dissection and characterization of the functional domains in cardiac troponin C

Contraction of vertebrate cardiac muscle is regulated by the binding of Ca$\sp{2+}$ to the troponin C (cTnC) subunit of the troponin complex. In this study, we have used site-directed mutagenesis and a variety of assay techniques to explore the functional roles of regions in cTnC, including Ca$\sp{2+}$/Mg$\sp{2+}$-binding sites III and IV, the functionally inactive site I, the N-terminal helix, the N-terminal hydrophobic pocket and the two cysteine residues with regard to their ability to form disulfide bonds. Conversion of the first Ca$\sp{2+}$ ligand from Asp to Ala inactivated sites III and IV and decreased the apparent affinity of cTnC for the thin filament. Conversion of the second ligand from Asn to Ala also inactivated these sites in the free protein but Ca$\sp{2+}$-binding was recovered upon association with troponin I and troponin T. The Ca$\sp{2+}$-concentrations required for tight thin filament-binding by proteins containing second-ligand mutations were significantly greater than that required for the wild-type protein. Mutation of site I such that the primary sequence was that of an active site with the first Ca$\sp{2+}$ ligand changed from Asp to Ala resulted in a 70% decrease in maximal Ca$\sp{2\sp+}$ dependent ATPase activity in both cardiac and fast skeletal myofibrils. Thus, the primary sequence of the inactive site I in cTnC is functionally important. Major changes in the sequence of the N-terminus had little effect on the ability of cTnC to recover maximal activity but deletion of the first nine residues resulted in a 60 to 80% decrease in maximal activity with only a minor decrease in the pCa$\sb{50}$ of activation, suggesting that the N-terminal helix must be present but that a specific sequence is not required. The formation of an inter- or intramolecular disulfide bonds caused the exposure of hydrophobic surfaces on cTnC and rendered the protein Ca$\sp{2+}$ independent. Finally, elution patterns from a hydrophobic interactions column suggest that cTnC undergoes a significant change in hydrophobicity upon Ca$\sp{2+}$ binding, the majority of which is caused by site II. These latter data show an interesting correlation between exposure of hydrophobic surfaces on and activation of cTnC. Overall, these results represent significant progress toward the elucidation of the functional roles of a variety of structural regions in cTnC. ^

Assessment of cone beam computed tomography techniques for imaging lung damage in mice in vivo

Lung damage is a common side effect of chemotherapeutic drugs such as bleomycin. This study used a bleomycin mouse model which simulates the lung damage observed in humans. Noninvasive, in vivo cone-beam computed tomography (CBCT) was used to visualize and quantify fibrotic and inflammatory damage over the entire lung volume of mice. Bleomycin was used to induce pulmonary damage in vivo and the results from two CBCT systems, a micro-CT and flat panel CT (fpCT), were compared to histologic measurements, the standard method of murine lung damage quantification. Twenty C57BL/6 mice were given either 3 U/kg of bleomycin or saline intratracheally. The mice were scanned at baseline, before the administration of bleomycin, and then 10, 14, and 21 days afterward. At each time point, a subset of mice was sacrificed for histologic analysis. The resulting CT images were used to assess lung volume. Percent lung damage (PLD) was calculated for each mouse on both the fpCT (PLDfpcT) and the micro-CT (PLDμCT). Histologic PLD (PLDH) was calculated for each histologic section at each time point (day 10, n = 4; day 14, n = 4; day 21, n = 5; control group, n = 5). A linear regression was applied to the PLDfpCT vs. PLDH, PLDμCT vs. PLDH and PLDfpCT vs. PLDμCT distributions. This study did not demonstrate strong correlations between PLDCT and PLDH. The coefficient of determination, R, was 0.68 for PLDμCT vs. PLDH and 0.75 for the PLD fpCT vs. PLDH. The experimental issues identified from this study were: (1) inconsistent inflation of the lungs from scan to scan, (2) variable distribution of damage (one histologic section not representative of overall lung damage), (3) control mice not scanned with each group of bleomycin mice, (4) two CT systems caused long anesthesia time for the mice, and (5) respiratory gating did not hold the volume of lung constant throughout the scan. Addressing these issues might allow for further improvement of the correlation between PLDCT and PLDH. ^

Improved Techniques for Acquisition and Analysis of Dynamic Contrast-Enhanced Magnetic Resonance Imaging for Detecting Vascular Permeability in the Central Nervous System

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) is a noninvasive technique for quantitative assessment of the integrity of blood-brain barrier and blood-spinal cord barrier (BSCB) in the presence of central nervous system pathologies. However, the results of DCE-MRI show substantial variability. The high variability can be caused by a number of factors including inaccurate T1 estimation, insufficient temporal resolution and poor contrast-to-noise ratio. My thesis work is to develop improved methods to reduce the variability of DCE-MRI results. To obtain fast and accurate T1 map, the Look-Locker acquisition technique was implemented with a novel and truly centric k-space segmentation scheme. In addition, an original multi-step curve fitting procedure was developed to increase the accuracy of T1 estimation. A view sharing acquisition method was implemented to increase temporal resolution, and a novel normalization method was introduced to reduce image artifacts. Finally, a new clustering algorithm was developed to reduce apparent noise in the DCE-MRI data. The performance of these proposed methods was verified by simulations and phantom studies. As part of this work, the proposed techniques were applied to an in vivo DCE-MRI study of experimental spinal cord injury (SCI). These methods have shown robust results and allow quantitative assessment of regions with very low vascular permeability. In conclusion, applications of the improved DCE-MRI acquisition and analysis methods developed in this thesis work can improve the accuracy of the DCE-MRI results.