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.

Natural trans-splicing in carnitine octanoyltransferase pre-mRNAs in rat liver

Carnitine octanoyltransferase (COT) transports medium-chain fatty acids through the peroxisome. During isolation of a COT clone from a rat liver library, a cDNA in which exon 2 was repeated, was characterized. Reverse transcription-PCR amplifications of total RNAs from rat liver showed a three-band pattern. Sequencing of the fragments revealed that, in addition to the canonical exon organization, previously reported [Choi, S. J. et al. (1995) Biochim. Biophys. Acta 1264, 215–222], there were two other forms in which exon 2 or exons 2 and 3 were repeated. The possibility of this exonic repetition in the COT gene was ruled out by genomic Southern blot. To study the gene expression, we analyzed RNA transcripts by Northern blot after RNase H digestion of total RNA. Three different transcripts were observed. Splicing experiments also were carried out in vitro with different constructs that contain exon 2 plus the 5′ or the 3′ adjacent intron sequences. Our results indicate that accurate joining of two exons 2 occurs by a trans-splicing mechanism, confirming the potential of these structures for this process in nature. The trans-splicing can be explained by the presence of three exon-enhancer sequences in exon 2. Analysis by Western blot of the COT proteins by using specific antibodies showed that two proteins corresponding to the expected Mr are present in rat peroxisomes. This is the first time that a natural trans-splicing reaction has been demonstrated in mammalian cells.

A dynamin GTPase mutation causes a rapid and reversible temperature-inducible locomotion defect in C. elegans

Drosophila shibire and its mammalian homologue dynamin regulate an early step in endocytosis. We identified a Caenorhabditis elegans dynamin gene, dyn-1, based upon hybridization to the Drosophila gene. The dyn-1 RNA transcripts are trans-spliced to the spliced leader 1 and undergo alternative splicing to code for either an 830- or 838-amino acid protein. These dyn-1 proteins are highly similar in amino acid sequence, structure, and size to the Drosophila and mammalian dynamins: they contain an N-terminal GTPase, a pleckstrin homology domain, and a C-terminal proline-rich domain. We isolated a recessive temperature-sensitive dyn-1 mutant containing an alteration within the GTPase domain that becomes uncoordinated when shifted to high temperature and that recovers when returned to lower temperatures, similar to D. shibire mutants. When maintained at higher temperatures, dyn-1 mutants become constipated, egg-laying defective, and produce progeny that die during embryogenesis. Using a dyn-1::lacZ gene fusion, a high level of dynamin expression was observed in motor neurons, intestine, and pharyngeal muscle. Our results suggest that dyn-1 function is required during development and for normal locomotion.

Influenza C virus CM2 integral membrane glycoprotein is produced from a polypeptide precursor by cleavage of an internal signal sequence

The influenza C virus CM2 protein is a small glycosylated integral membrane protein (115 residues) that spans the membrane once and contains a cleavable signal sequence at its N terminus. The coding region for CM2 (CM2 ORF) is located at the C terminus of the 342-amino acid (aa) ORF of a colinear mRNA transcript derived from influenza C virus RNA segment 6. Splicing of the colinear transcript introduces a translational stop codon into the ORF and the spliced mRNA encodes the viral matrix protein (CM1) (242 aa). The mechanism of CM2 translation was investigated by using in vitro and in vivo translation of RNA transcripts. It was found that the colinear mRNA derived from influenza C virus RNA segment 6 serves as the mRNA for CM2. Furthermore, CM2 translation does not depend on any of the three in-frame methionine residues located at the beginning of CM2 ORF. Rather, CM2 is a proteolytic cleavage product of the p42 protein product encoded by the colinear mRNA: a cleavage event that involves the recognition and cleavage of an internal signal peptide presumably by signal peptidase resident in the endoplasmic reticulum. Alteration of the predicted signal peptidase cleavage site by mutagenesis blocked generation of CM2. The other polypeptide species resulting from the cleavage of p42, designated p31, contains the CM1 coding region and an additional C-terminal 17 aa (formerly the CM2 signal peptide). Protein p31, in comparison to CM1, displays characteristics of an integral membrane protein.

Capping and methylation of mRNA by purified recombinant VP4 protein of bluetongue virus

The core of bluetongue virus (BTV) is a multienzyme complex composed of two major proteins (VP7 and VP3) and three minor proteins (VP1, VP4, and VP6) in addition to the viral genome. The core is transcriptionally active and produces capped mRNA from which all BTV proteins are translated, but the relative role of each core component in the overall reaction process remains unclear. Previously we showed that the 76-kDa VP4 protein possesses guanylyltransferase activity, a necessary part of the RNA capping reaction. Here, through the use of highly purified (>95%) VP4 and synthetic core-like particles containing VP4, we have investigated the extent to which this protein is also responsible for other activities associated with cap formation. We show that VP4 catalyzes the conversion of unmethylated GpppG or in vitro-produced uncapped BTV RNA transcripts to m7GpppGm in the presence of S-adenosyl-l-methionine. Analysis of the methylated products of the reaction by HPLC identified both methyltransferase type 1 and type 2 activities associated with VP4, demonstrating that the complete BTV capping reaction is associated with this one protein.

Ribosomes can slide over and beyond “hungry” codons, resuming protein chain elongation many nucleotides downstream

In cells subjected to moderate aminoacyl-tRNA limitation, the peptidyl-tRNA–ribosome complex stalled at the “hungry” codon can slide well beyond it on the messenger RNA and resume translation further downstream. This behavior is proved by unequivocal amino acid sequence data, showing a protein that lacks the bypassed sequence encoded between the hungry codon and specific landing sites. The landing sites are codons cognate to the anticodon of the peptidyl-tRNA. The efficiency of this behavior can be as high as 10–20% but declines with the length of the slide. Interposition of “trap” sites (nonproductive landing sites) in the bypassed region reduces the frequency of successful slides, confirming that the ribosome–peptidyl-tRNA complex passes through the untranslated region of the message. This behavior appears to be quite general: it can occur at the two kinds of hungry codons tested, AUA and AAG; the sliding peptidyl-tRNA can be any of three species tested, phenylalanine, tyrosine, or leucine tRNA; the peptidyl component can be either of two very different peptide sequences; and translation can resume at any of the three codons tested.

Inhibition of mRNA export in vertebrate cells by nuclear export signal conjugates

Leucine-rich nuclear export signals (NESs) are recognized by the NES receptor exportin 1 and are central to the export of multiple shuttling proteins and RNAs. The export of messenger RNA in vertebrates was, however, thought to occur by a different pathway, because inhibition by injection of a synthetic Rev NES conjugate could not be demonstrated. Here we find that peptide conjugates composed of the NES of either protein kinase A inhibitor protein (PKI) or the HIV-1 Rev protein, when coupled to human serum albumin, are potent inhibitors of mRNA and small nuclear RNA export. These results provide direct evidence that mRNA export in vertebrates depends on interactions between an NES and its cognate NES receptors. PKI NES conjugates are significantly more efficient at inhibiting RNA export than are REV NES conjugates, indicating that different NESs may have different abilities to promote protein and RNA export. Surprisingly, an expected control conjugate containing the mutant Rev NES sequence M10 strongly inhibited the export of intronless dihydrofolate reductase mRNA. Nuclear injection of NES peptide conjugates led to mislocalization to the nucleus of 10–20% of the cytoplasmic Ran GTPase-binding protein (RanBP1) indicating that RanBP1 shuttles between the nucleus and the cytoplasm via an NES pathway. These results demonstrate that in vertebrates the export of mRNA, like that of small nuclear RNA, 5S rRNA, and transport factors such as RanBP1, employs NES-mediated molecular machinery.

In vivo analysis of the 3′ untranslated region of the hepatitis C virus after in vitro mutagenesis of an infectious cDNA clone

Large sections of the 3′ untranslated region (UTR) of hepatitis C virus (HCV) were deleted from an infectious cDNA clone, and the RNA transcripts from seven deletion mutants were tested sequentially for infectivity in a chimpanzee. Mutants lacking all or part of the 3′ terminal conserved region or the poly(U–UC) region were unable to infect the chimpanzee, indicating that both regions are critical for infectivity in vivo. However, the third region, the variable region, was able to tolerate a deletion that destroyed the two putative stem–loop structures within this region. Mutant VR-24 containing a deletion of the proximal 24 nt of the variable region of the 3′ UTR was viable in the chimpanzee and seemed to replicate as well as the undeleted parent virus. The chimpanzee became viremic 1 week after inoculation with mutant VR-24, and the HCV genome titer increased over time during the early acute infection. Therefore, the poly(U–UC) region and the conserved region, but not the variable region, of the 3′ UTR seem to be critical for in vivo infectivity of HCV.

Hepatitis C virus lacking the hypervariable region 1 of the second envelope protein is infectious and causes acute resolving or persistent infection in chimpanzees

Persistent infection with hepatitis C virus (HCV) is among the leading causes of chronic liver disease. Previous studies suggested that genetic variation in hypervariable region 1 (HVR1) of the second envelope protein, possibly in response to host immune pressure, influences the outcome of HCV infection. In the present study, a chimpanzee transfected intrahepatically with RNA transcripts of an infectious HCV clone (pCV-H77C) from which HVR1 was deleted became infected; the ΔHVR1 virus was subsequently transmitted to a second chimpanzee. Infection with ΔHVR1 virus resulted in persistent infection in the former chimpanzee and in acute resolving infection in the latter chimpanzee. Both chimpanzees developed hepatitis. The ΔHVR1 virus initially replicated to low titers, but virus titer increased significantly after mutations appeared in the viral genome. Thus, wild-type HCV without HVR1 was apparently attenuated, suggesting a functional role of HVR1. However, our data indicate that HVR1 is not essential for the viability of HCV, the resolution of infection, or the progression to chronicity.

Production of zebrafish germ-line chimeras from embryo cell cultures

Although the zebrafish possesses many characteristics that make it a valuable model for genetic studies of vertebrate development, one deficiency of this model system is the absence of methods for cell-mediated gene transfer and targeted gene inactivation. In mice, embryonic stem cell cultures are routinely used for gene transfer and provide the advantage of in vitro selection for rare events such as homologous recombination and targeted mutation. Transgenic animals possessing a mutated copy of the targeted gene are generated when the selected cells contribute to the germ line of a chimeric embryo. Although zebrafish embryo cell cultures that exhibit characteristics of embryonic stem cells have been described, successful contribution of the cells to the germ-cell lineage of a host embryo has not been reported. In this study, we demonstrate that short-term zebrafish embryo cell cultures maintained in the presence of cells from a rainbow trout spleen cell line (RTS34st) are able to produce germ-line chimeras when introduced into a host embryo. Messenger RNA encoding the primordial germ-cell marker, vasa, was present for more than 30 days in embryo cells cocultured with RTS34st cells or their conditioned medium and disappeared by 5 days in the absence of the spleen cells. The RTS34st cells also inhibited melanocyte and neuronal cell differentiation in the embryo cell cultures. These results suggest that the RTS34st splenic–stromal cell line will be a valuable tool in the development of a cell-based gene transfer approach to targeted gene inactivation in zebrafish.

Molecular Cloning and Expression Analysis of the Mitochondrial Pyruvate Dehydrogenase from Maize1

Four cDNAs, one encoding an α-subunit and three encoding β-subunits of the mitochondrial pyruvate dehydrogenase, were isolated from maize (Zea mays L.) libraries. The deduced amino acid sequences of both α- and β-subunits are approximately 80% identical with Arabidopsis and pea (Pisum sativum L.) homologs. The mature N terminus was determined for the β-subunit by microsequencing the protein purified from etiolated maize shoot mitochondria and was resolved by two-dimensional gel electrophoresis. This single isoelectric species comprised multiple isoforms. Both α- and β-subunits are encoded by multigene families in maize, as determined by Southern-blot analyses. RNA transcripts for both α- and β-subunits were more abundant in roots than in young leaves or etiolated shoots. Pyruvate dehydrogenase activity was also higher in roots (5-fold) compared with etiolated shoots and leaves. Both subunits were present at similar levels in all tissues examined, indicating coordinated gene regulation. The protein levels were highest in heterotrophic organs and in pollen, which contained about 2-fold more protein than any other organ examined. The relative abundance of these proteins in nonphotosynthetic tissues may reflect a high cellular content of mitochondria, a high level of respiratory activity, or an extra plastidial requirement for acetate.

Highly stable expression of a foreign gene from rabies virus vectors.

A reverse genetics approach was applied to generate a chimeric nonsegmented negative strand RNA virus, rabies virus (RV) of the Rhabdoviridae family, that expresses a foreign protein. DNA constructs containing the entire open reading frame of the bacterial chloramphenicol acetyltransferase (CAT) gene and an upstream RV cistron border sequence were inserted either into the nontranslated pseudogene region of a full-length cDNA copy of the RV genome or exchanged with the pseudogene region. After intracellular T7 RNA polymerase-driven expression of full-length antigenome RNA transcripts and RV nucleoprotein, phosphoprotein and polymerase from transfected plasmids, RVs transcribing novel monocistronic mRNAs and expressing CAT at high levels, were recovered. The chimeric viruses possessed the growth characteristics of standard RV and were genetically stable upon serial cell culture passages. CAT activity was still observed in cell cultures infected with viruses passaged for more than 25 times. Based on the unprecedented stability of the chimeric RNA genomes, which is most likely due to the structure of the rhabdoviral ribonucleoprotein complex, we predict the successful future use of recombinant rhabdovirus vectors for displaying foreign antigens or delivering therapeutic genes.

Long-term culture of lymphohematopoietic stem cells.

Pluripotent hematopoietic stem cells (PHSCs) show self-renewal and give rise to all blood cell types. The extremely low number of these cells in primary hematopoietic organs and the lack of culture systems that support proliferation of undifferentiated PHSCs have precluded the study of both the biology of these cells and their clinical application. We describe here cell lines and clones derived from PHSCs that were established from hematopoietic cells from the fetal liver or bone marrow of normal and p53-deficient mice with a combination of four growth factors. Most cell lines were Sca-1+, c-Kit+, PgP-1+, HSA+, and Lin- (B-220-, Joro 75-, 8C5-, F4/80-, CD4-, CD8-, CD3-, IgM-, and TER 119-negative) and expressed three new surface markers: Joro 177, Joro 184, and Joro 96. They did not synthesize RNA transcripts for several genes expressed at early stages of lymphocyte and myeloid/erythroid cell development. The clones were able to generate lymphoid, myeloid, and erythroid hematopoietic cells and to reconstitute the hematopoietic system of irradiated mice for a long time. The availability of lymphohematopoietic stem cell lines should facilitate the analysis of the molecular mechanisms that control self-renewal and differentiation and the development of efficient protocols for somatic gene therapy.

Challenging the dogma: the hidden layer of non-protein-coding RNAs in complex organisms

The central dogma of biology holds that genetic information normally flows from DNA to RNA to protein. As a consequence it has been generally assumed that genes generally code for proteins, and that proteins fulfil not only most structural and catalytic but also most regulatory functions, in all cells, from microbes to mammals. However, the latter may not be the case in complex organisms. A number of startling observations about the extent of non-protein-coding RNA (ncRNA) transcription in the higher eukaryotes and the range of genetic and epigenetic phenomena that are RNA-directed suggests that the traditional view of the structure of genetic regulatory systems in animals and plants may be incorrect. ncRNA dominates the genomic output of the higher organisms and has been shown to control chromosome architecture, mRNA turnover and the developmental timing of protein expression, and may also regulate transcription and alternative splicing. This paper re-examines the available evidence and suggests a new framework for considering and understanding the genomic programming of biological complexity, autopoletic development and phenotypic variation. BioEssays 25:930-939,2003. (C) 2003 Wiley Periodicals, Inc.

Vitamin D3 and brain development

Evidence for the presence of the vitamin D receptor in brain implies this vitamin may have some function in this organ. This study investigates whether vitamin D-3 acts during brain development. We demonstrate that rats born to vitamin D-3-deficient mothers had profound alterations in the brain at birth. The cortex was longer but not wider, the lateral ventricles were enlarged, the cortex was proportionally thinner and there was more cell proliferation throughout the brain. There were reductions in brain content of nerve growth factor and glial cell line-derived neurotrophic factor and reduced expression of p75(NTR), the low-affinity neurotrophin receptor. Our findings would suggest that low maternal vitamin D3 has important ramifications for the developing brain. (C) 2003 IBRO. Published by Elsevier Science Ltd. All rights reserved.