Molecular structure of antibodies against small ligands

Antibodies which bind bioactive ligands can serve as a template for the generation of a second antibody which may react with the physiological receptor. This phenomenon of molecular mimicry by antibodies has been described in a variety of systems. In order to understand the chemical and molecular mechanisms involved in these interactions, monoclonal antibodies directed against two pharmacologically active alkaloids, morphine and nicotine, were carefully studied using experimental and theoretical molecular modeling techniques. The molecular characterization of these antibodies involved binding studies with ligand analogs and determination of the variable region amino acid sequence. A three-dimensional model of the anti-morphine binding site was constructed using computational and graphics display techniques. The antibody response in BALB/c mice to morphine appears relatively restricted, in that all of the antibodies examined in this study contained a $\lambda$ light chain, which is normally found in only 5% of mouse immunoglobulins. This study represents the first use of theoretical and experimental modeling techniques to describe the antigen binding site of a mouse Fv region containing a $\lambda$ light chain. The binding site model indicates that a charged glutamic acid residue and aromatic side chains are key features in ionic and hydrophobic interactions with the ligand morphine. A glutamic acid residue is found in the identical position in the anti-nicotine antibody and may play a role in binding nicotine. ^

Stress-induced targeting of molecular chaperones in the yeast Saccharomyces cerevisiae

The eukaryotic stress response is an essential mechanism that helps protect cells from a variety of environmental stresses. Cell death can result if cells are not able to properly adapt and protect themselves against adverse stress conditions. Failure to properly deal with stress has implications in human diseases including neurodegenerative disorders and distinct cancers, emphasizing the importance of understanding the eukaryotic stress response in detail. As part of this response, expression of a battery of heat shock proteins (HSP) is induced, which act as molecular chaperones to assist in the repair or triage of unfolded proteins. The 90-kDa HSP (Hsp90) operates in the context of a multi-chaperone complex to promote the maturation of nuclear and cytoplasmic clients. I have discovered that Hsp90 and the co-chaperone Sba1 accumulate in the nucleus of quiescent Saccharomyces cerevisiae cells in a karyopherin-dependent manner. I isolated nuclear accumulation- defective HSP82 mutant alleles to probe the nature of this targeting event and identified a mutant with a single amino acid substitution (I578F) sufficient to prevent nuclear accumulation of Hsp90 in quiescent cells. Diploid hsp82-I578F cells exhibited pronounced defects in spore wall construction and maturation, resulting in catastrophic sporulation. The mislocalization and sporulation phenotypes were shared by another previously identified HSP82 mutant allele, further linking localization to Hsp90 functional status. Pharmacological inhibition of Hsp90 with macbecin in sporulating diploid cells also blocked spore formation, underscoring the importance of this chaperone in this developmental program. The yeast molecular chaperone Hsp104 is a member of the Hsp100 superfamily of AAA+ ATPases. Unlike the Hsp90 family of chaperones, Hsp104 is not restricted to a specific set of client proteins, but rather assists in reactivating stress-denatured proteins by solubilizing protein aggregates. I have discovered that Hsp104, along with the Hsp70 chaperone, Ssa1, and the sHSP Hsp26 accumulate into RNA processing bodies (P- bodies) and stress granules, sites of mRNA metabolism. I found that Hsp104 recruits both Ssa1 and Hsp26 to P-bodies and that these three chaperones are required for stress granule formation. These findings suggest a possible role for chaperones in mRNA metabolism by aiding in the assembly, disassembly or conversion of these enigmatic mRNP complexes. Taken together, the work presented in this dissertation serves to better understand the eukaryotic stress response by illustrating the importance of subcellular-chaperone localization in key biological processes.

Molecular epidemiology of vancomycin-resistant Enterococcus faecium: a prospective, multicenter study in South American hospitals.

Enterococcus faecium has emerged as an important nosocomial pathogen worldwide, and this trend has been associated with the dissemination of a genetic lineage designated clonal cluster 17 (CC17). Enterococcal isolates were collected prospectively (2006 to 2008) from 32 hospitals in Colombia, Ecuador, Perú, and Venezuela and subjected to antimicrobial susceptibility testing. Genotyping was performed with all vancomycin-resistant E. faecium (VREfm) isolates by pulsed-field gel electrophoresis (PFGE) and multilocus sequence typing. All VREfm isolates were evaluated for the presence of 16 putative virulence genes (14 fms genes, the esp gene of E. faecium [espEfm], and the hyl gene of E. faecium [hylEfm]) and plasmids carrying the fms20-fms21 (pilA), hylEfm, and vanA genes. Of 723 enterococcal isolates recovered, E. faecalis was the most common (78%). Vancomycin resistance was detected in 6% of the isolates (74% of which were E. faecium). Eleven distinct PFGE types were found among the VREfm isolates, with most belonging to sequence types 412 and 18. The ebpAEfm-ebpBEfm-ebpCEfm (pilB) and fms11-fms19-fms16 clusters were detected in all VREfm isolates from the region, whereas espEfm and hylEfm were detected in 69% and 23% of the isolates, respectively. The fms20-fms21 (pilA) cluster, which encodes a putative pilus-like protein, was found on plasmids from almost all VREfm isolates and was sometimes found to coexist with hylEfm and the vanA gene cluster. The population genetics of VREfm in South America appear to resemble those of such strains in the United States in the early years of the CC17 epidemic. The overwhelming presence of plasmids encoding putative virulence factors and vanA genes suggests that E. faecium from the CC17 genogroup may disseminate in the region in the coming years.

Molecular mechanisms underlying a cellular analog of operant reward learning.

Operant conditioning is a ubiquitous but mechanistically poorly understood form of associative learning in which an animal learns the consequences of its behavior. Using a single-cell analog of operant conditioning in neuron B51 of Aplysia, we examined second-messenger pathways engaged by activity and reward and how they may provide a biochemical association underlying operant learning. Conditioning was blocked by Rp-cAMP, a peptide inhibitor of PKA, a PKC inhibitor, and by expressing a dominant-negative isoform of Ca2+-dependent PKC (apl-I). Thus, both PKA and PKC were necessary for operant conditioning. Injection of cAMP into B51 mimicked the effects of operant conditioning. Activation of PKC also mimicked conditioning but was dependent on both cAMP and PKA, suggesting that PKC acted at some point upstream of PKA activation. Our results demonstrate how these molecules can interact to mediate operant conditioning in an individual neuron important for the expression of the conditioned behavior.

Molecular basis of maturation promoting factor

Maturation promoting factor (MPF), which is functionally defined by its ability to induce Xenopus oocyte maturation, is an M phase (meiosis and mitosis) specific activity that is present in all species tested. It was hypothesized that MPF is a universal trigger of the interphase to M phase transition during the cell cycle. The current model for the molecular basis of MPF is that MPF is a protein kinase having the cdc2 protein as its catalytic subunit and is identical to the M phase-specific histone H1 kinase. In the present study, I have shown that more than just cdc2 kinase contributes to MPF activity, and M phase-specific H1 kinase is composed of at least two entities, instead of just cdc2 kinase. Therefore, the simple model of MPF = cdc2 kinase = M phase-specific H1 kinase should be ruled out.^ My study began with the characterization of the mitosis-specific monoclonal antibody MPM-2. MPM-2 reacts specifically with M phase cells from different species by recognizing a discrete set of proteins once they are phosphorylated at the G$\sb2$/M transition. I found that phosphorylation of MPM-2 antigens coincided with the appearance of MPF activity during oocyte maturation stimulated by progesterone. If MPM-2 was injected into oocytes before the stimulation, MPF activity failed to appear, and the oocytes could not mature. Furthermore, MPM-2 was able to deplete MPF activity from M phase extracts. These results identified MPM-2 as a probe that recognizes either MPF itself or a regulator of MPF.^ Since M phase-specific H1 kinase was believed to be identical to cdc2 kinase and MPF, I proceeded to determine whether MPM-2 recognized the M phase-specific H1 kinase. I found that MPM-2 did recognize an M phase-specific H1 kinase. However, this kinase was not cdc2 kinase. This kinase (MPM-2 kinase) is present in a latent form in immature oocytes and is activated in tandem with the activation of MPF during oocyte maturation. It appears to accelerate progesterone-induced oocyte maturation. Therefore, MPM-2 kinase may be a novel positive regulator of MPF activation.^ MPM-2 depletes MPF activity, but not cdc2 kinase activity. This discrepancy caused me to question the equivalency of MPF with cdc2 kinase. I found that when a high percentage of MPF activity was recovered from gel filtration of mature oocyte extract, the recovered MPF activity was due to two factors, cdc2 kinase and a factor recognized by MPM-2. This factor might activate and stabilize cdc2 kinase. Identification of this factor in the present study may contribute to the understanding of the autoactivation of MPF. ^

Molecular basis of retinoid action: The role of retinoic acid receptors in retinoic acid-induced tissue transglutaminase expression

Retinoic acid has profound effects on the cellular growth and differentiation of a variety of cells. However, the molecular basis of retinoic acid action has, until recently, not been well understood. The identification of retinoic acid receptors which bear a high degree of homology to members of the steroid receptor super-family has dramatically altered our understanding of the biology of retinoids. The focus of this dissertation has been toward identification of retinoic acid binding proteins responsible for the effects of this molecule on gene expression.^ We have characterized in detail the retinoic acid-dependent induction of tissue transglutaminase gene expression in a myeloid cell line, human promyelocytic leukemia cells (HL-60 cells). Using cDNA probes specific for tissue transglutaminase, we have determined that the retinoic acid induced increase in enzyme level is due to an increase in the level of tissue transglutaminase mRNA. We have used this model as a probe to investigate the molecular basis of retinoid regulated gene expression.^ This thesis demonstrates that retinoic acid receptors are expressed in cells which induce tissue transglutaminase expression in response to retinoic acid. In Hl-60 cells retinoic acid-induced transglutaminase expression is associated with saturable nuclear retonic acid binding. Transcripts for both the alpha and beta forms of the retinoic acid receptors can be detected in these cells. Pretreatment of HL-60 cells with agents that potentiate retinoic acid-induced transglutaminase expression also modestly induced the alpha form of the retinoic acid receptor. Studies in macrophages and umbilical vein endothelial cells have also associated expression of the beta form of the retinoic acid with retinoic acid induced tissue transglutaminase expression.^ To investigate directly if retinoic acid receptors regulate retinoic acid-induced tissue transglutaminase expression we developed a series of stably transfected Balb-c 3T3 cells expressing different levels of the beta or gamma form of the retinoic acid receptor. These studies indicated that either the beta or gamma receptor can stimulate endogenous tissue transglutaminase expression in response to retinoic acid. These are among the first studies in the steroid field to describe regulation of an endogenous gene by a transfected receptor. ^

Mapping and cloning of genes from portions of the human genome using somatic cell hybrids

Human x rodent somatic cell hybrids have played an important role in human genetics research. They have been especially useful for assigning genes to chromosomes and isolating DNA markers from specific regions of the human genome.^ By employing a combination of somatic cell genetic, recombinant DNA, and cytogenetic techniques, human DNA excision repair gene ERCC4 was mapped regionally to human 16p13.13-13.2, even though the gene has not been cloned. Human x Chinese hamster ovary (CHO) cell hybrids selected for human ERCC4 activity and containing 16p13.1-p13.3 as the only human genetic material were identified. These hybrids were used to order DNA markers located in 16p13.1-p13.3. New DNA markers physically close to ERCC4 were isolated from such hybrids. Using amplified human DNA from the hybrids as probe in fluorescent in situ hybridization, the short arm breakpoint in the chromosome 16 inversion associated with acute myelomonocytic leukemia (AMML) was found to be physically close to the ERCC4 gene. The physical mapping and eventually, the cloning of the ERCC4 gene, will benefit the understanding of the DNA repair system and the study of other important biomedical problems such as tumorigenesis.^ To facilitate the cloning of ERCC4 gene and, in general, the cloning of genes from any defined regions of the human genome, a method was developed for the direct isolation of human transcribed genes ffom somatic cell hybrids. cDNA was prepared from human x rodent hybrid by using consensus 5$\sp\prime$ splice site sequences as primers. These primers were designed to select immature, unspliced messenger RNA (still retaining species specific repeat sequences) as templates. Screening of a derived cDNA library for human repeat sequences resulted in the isolation of human clones at the anticipated frequency with characteristics expected of exons of transcribed human genes. The usefulness of the splice site specific primers was analyzed and the cDNA synthesis conditions with these primers were optimized. The procedure was shown to be sensitive enough to clone weakly expressed genes. Studying the expression of the represented genes with the isolated clones was shown to be feasible. Such regional specific human gene fragments will be very valuable for many human genetic studies such as the search of inherited disease genes and the construction of a cDNA map of the human genome. ^

Molecular mechanisms of Ca$\sp{2+}$ signal transduction from cardiac troponin C to troponin I

$\rm Ca\sp{2+}$-dependent exposure of an N-terminal hydrophobic region in troponin C (TnC) is thought to be important for the regulation of contraction in striated muscle. To study these conformational changes in cardiac troponin (cTnC), the $\varepsilon$C and $\varepsilon$H chemical shifts for all 10 Met residues in cTnC were sequence-specific assigned on NMR spectra using a combination of two dimensional NMR techniques and site-directed mutagenesis. The assigned methyl-Met chemical shifts were used as structural markers to monitor conformational changes induced by $\rm Ca\sp{2+}.$ The results showed that binding of $\rm Ca\sp{2+}$ to the regulatory site in the N-domain induced large changes in the $\varepsilon$H and $\varepsilon$C chemical shifts of Met 45, Met 80, Met 81 in the predicted N-terminal hydrophobic region, but had no effect on the chemical shifts of Met residues located in the C-domain. These results suggest that the $\rm Ca\sp{2+}$-dependent functions of cTnC are mainly through N-terminal domain of cTnC.^ To further define the molecular mechanism by which TnC regulates muscle contraction, single Cys residues were engineered at positions 45, 81, 84 or 85 in the N-terminal hydrophobic region of cTnC to provide sites for attachment of specific blocking groups. Blocking groups were coupled to these Cys residues in cTnC mutants and the covalent adducts were tested for activity in TnC-extracted myofibrils. Covalent modification of cTnC(C45) had no effect on maximal myofibril ATPase activity. Greatly decreased myofibril ATPase activity resulted when the peptide or biotin was conjugated to residue 81 in cTnC(C81), while less inhibition resulted from covalent modification of cTnC(C84) or cTnC(C85). The results suggest that limited sites of the N-terminal hydrophobic region in cTnC are important for transducing the $\rm Ca\sp{2+}$ signal to troponin I (TnI) and are sensitive to modification, while other regions are less important or can adapt to steric hindrances introduced by bulky blocking groups.^ Although the exposed TnI interaction site in the N-terminal hydrophobic region of TnC is crucial for function of TnC, other regions in the N-domain of TnC may also participate in transducing the $\rm Ca\sp{2+}$ signal and conferring the maximal activation of actomyosin ATPase. The interactions between the B-/C-helices of cTnC and cTnI were characterized using a combination of site-directed mutagenesis, fluorescence and covalent modification. The results suggest that the $\rm Ca\sp{2+}$-dependent interactions of the B-/C-helices of cTnC with TnI may be required for the maximal activation of muscle contraction. ^

STUDIES INTO THE MOLECULAR MECHANISMS REGULATING MUSCULAR ATROPHY RESULTING FROM HINDLIMB IMMOBILIZATION IN THE RAT (PROTEIN, SYNTHESIS, RNA)

The purpose of the work performed in this dissertation was to examine some of the possible regulatory mechanisms involved in the initiation of muscular atrophy during periods of decreased muscle utilization resulting from hindlimb immobilization in the rat. A 37% decrease in the rate of total muscle protein synthesis which has been observed to occur in the first 6 h of immobilization contributes significantly to the observed loss of protein during immobilization.^ The rates of cytochrome c and actin synthesis were determined in adult rat red vastus lateralis and gastrocnemius muscles, respectively, by the constant infusion and incorporation of ('3)H-tyrosine into protein. The fractional synthesis rates of both actin and cytochrome c were significantly decreased (P < 0.05) in the 6th h of hindlimb immobilization.^ RHA was extracted from adult rat gastrocnemius muscle by modification of the phenol: chloroform: SDS extraction procedures commonly used for preparation of RNA for hybridization analysis from other mammalian tissues. RNA content of rat gastrocnemius muscle, as determined by this method of extraction and its subsequent quantification by UV absorbance and orcinol assay, was significantly greater than the RNA content previously determined for adult rat gastrocnemius by other commonly employed methods.^ RNA extracted by this method from gastrocnemius muscles of control and 6h immobilized rats was subjected to "dot blot" hybridization to ('32)P-labelled probe from plasmid p749, containing a cDNA sequence complementary to (alpha)-actin mRNA and from rat skeletal muscle. (alpha)-Actin specific mRNA content as estimated by this procedure is not significantly decreased in rat gastrocnemius following 6h or hindlimb immobilization. However, (alpha)-actin specific mRNA content is significantly decreased (P < 0.05) in adult rat gastrocnemius (alpha)-actin specific mRNA is not decreased in adult rat gastrocnemius muscle following 6h of immobilization, a time when actin synthesis is significantly decreased, it is concluded that a change in (alpha)-actin specific mRNA content is not the initiating event responsible for the early decrease in actin synthesis observed in the 6th h of immobilization. ^

STUDIES OF THE MOLECULAR MECHANISMS OF CHROMOSOME ABERRATION FORMATION (DNA REPAIR, CROSSLINKS, MITOMYCIN C)

The objective of this research has been to study the molecular basis for chromosome aberration formation. Predicated on a variety of data, Mitomycin C (MMC)-induced DNA damage has been postulated to cause the formation of chromatid breaks (and gaps) by preventing the replication of regions of the genome prior to mitosis. The basic protocol for these experiments involved treating synchronized Hela cells in G(,1)-phase with a 1 (mu)g/ml dose of MMC for one hour. After removing the drug, cells were then allowed to progress to mitosis and were harvested for analysis by selective detachment. Utilizing the alkaline elution assay for DNA damage, evidence was obtained to support the conclusion that Hela cells can progress through S-phase into mitosis with intact DNA-DNA interstrand crosslinks. A higher level of crosslinking was observed in those cells remaining in interphase compared to those able to reach mitosis at the time of analysis. Dual radioisotope labeling experiments revealed that, at this dose, these crosslinks were associated to the same extent with both parental and newly replicated DNA. This finding was shown not to be the result of a two-step crosslink formation mechanism in which crosslink levels increase with time after drug treatment. It was also shown not to be an artefact of the double-labeling protocol. Using neutral CsCl density gradient ultracentrifugation of mitotic cells containing BrdU-labeled newly replicated DNA, control cells exhibited one major peak at a heavy/light density. However, MMC-treated cells had this same major peak at the heavy/light density, in addition to another minor peak at a density characteristic for light/light DNA. This was interpreted as indicating either: (1) that some parental DNA had not been replicated in the MMC treated sample or; (2) that a recombination repair mechanism was operational. To distinguish between these two possibilities, flow cytometric DNA fluorescence (i.e., DNA content) measurements of MMC-treated and control cells were made. These studies revealed that the mitotic cells that had been treated with MMC while in G(,1)-phase displayed a 10-20% lower DNA content than untreated control cells when measured under conditions that neutralize chromosome condensation effects (i.e., hypotonic treatment). These measurements were made under conditions in which the binding of the drug, MMC, was shown not to interfere with the stoichiometry of the ethidium bromide-mithramycin stain. At the chromosome level, differential staining techniques were used in an attempt to visualize unreplicated regions of the genome, but staining indicative of large unreplicated regions was not observed. These results are best explained by a recombinogenic mechanism. A model consistent with these results has been proposed.^

Signal relay between two integral membrane proteins: The sensory rhodopsin I/HtrI transducer molecular complex

The molecular complex of sensory rhodopsin I (SRI) and its transducer HtrI mediate color-sensitive phototaxis in the archaeon Halobacterium salinarum. Orange light causes an attractant response by a one-photon reaction and white light causes a repellent response by a two-photon reaction. Three aspects of this molecular complex were explored: (i) We determined the stoichiometry of SRI and HtrI to be 2:2 by gene fusion analysis. A SRI-HtrI fusion protein was expressed in H. salinarum and shown to mediate 1-photon and 2-photon phototaxis responses comparable to wild-type complex. Disulfide crosslinking demonstrated that the fusion protein is a homodimer in the membrane. Measurement of photochemical reaction kinetics and pH titration of absorption spectra established that both SRI domains are complexed to HtrI in the fusion protein, and therefore the stoichiometry is 2:2. (ii) Cytoplasmic channel closure of SRI by HtrI, an important aspect of their interaction, was investigated by incremental HtrI truncation. We found that binding of the membrane-embedded portion of HtrI is insufficient for channel closure, whereas cytoplasmic extension of the second HtrI transmembrane helix by 13 residues blocks proton conduction through the channel as well as full-length HtrI. The closure activity is localized to 5 specific residues, each of which incrementally contributes to reduction of proton conductivity. Moreover, these same residues in the dark incrementally and proportionally increase the pKa of the Asp76 counterion to the protonated Schiff base chromophore. We conclude that this critical region of HtrI alters the dark conformation of SRI as well as light-induced channel opening. (iii) We developed a procedure for reconstituting HtrI-free SRI and the SRI/HtrI complex into liposomes, which exhibit photocycles with opened and closed cytoplasmic channels, respectively, as in the membrane. This opens the way for study of the light-induced conformational change and the interaction in vitro by fluorescence and spin-labeling. Single-cysteine mutations were introduced into helix F of SRI, labeled with a nitroxide spin probe and a fluorescence probe, reconstituted into proteoliposomes, and light-induced conformational changes detected in the complex. The probe signals can now be used as the readout of signaling to analyze mutants and the kinetics of signal relay. ^

Molecular interactions and signal transfer between sensory rhodopsin-I and its cognate transducer

SRI is unique among known photoreceptors in that it produces opposite signals depending on the color of light stimuli. Absorption of orange light (587 nm) triggers an attractant response by the cell, whereas absorption of orange light followed by near-UV light (373 run) triggers a repellent response. Using behavioral mutants that exhibit aberrant color-sensing ability, we tested a two-conformation equilibrium model, using FRET and EPR spectroscopy. The essence of the model applied to SRI-HtrI is that the complex exists in a metastable two-conformer equilibrium which is shifted in one direction by orange light absorption (producing an attractant signal) and in the opposite direction by a second UV-violet photon (producing a repellent signal). First, by FRET we found that the E-F cytoplasmic loop of SRI moves toward the RAMP domain of the HtrI transducer during the formation of the orange-light activated signaling state of the complex. This is the first localization of a change in the physical relationship between the receptor and transducer subunits of the complex and provides a structural property of the two proposed conformers that we can monitor. Second, EPR spectra of a spin label probe at this cytoplasmic position showed shifts in the dark in the mutants toward shorter or longer EF loop-RAMP distances, explaining their behavior in terms of their mutations causing pre-stimulus shifts into one or the other conformer. ^ Next, we applied a novel electrophysiological method for monitoring the directionality of proton movement during photoactivation of SRI, to investigate the process of proton transfer in the photoactive site from the chromophore to proton acceptors on both the wildtype and aberrant color-response mutants. We observed an unexpected and critical difference in the two signaling conformations of the SRI-HtrI complex. The finding is that the vectoriality (i.e. movement away or toward the cytoplasm) of the light-induced proton transfer from the chromophore to the protein is opposite in formation of the two conformations. Retinylidene proton transfer is a common critical process in rhodopsins and these results are the first to show differences in vectoriality in a rhodopsin receptor, and to demonstrate functional importance of the direction of proton transfer. ^

Characterization of fecal contamination targeting 16S rRNA bacteroides molecular markers in the Santa Cruz River, Arizona

Understanding the origins, transport and fate of contamination is essential to effective management of water resources and public health. Individuals and organizations with management responsibilities need to understand the risks to ecosystems and to humans from contact with contamination. Managers also need to understand how key contaminants vary over time and space in order to design and prioritize mitigation strategies. Tumacacori National Historic Park (NHP) is responsible for management of its water resources for the benefit of the park and for the health of its visitors. The existence of microbial contaminants in the park poses risks that must be considered in park planning and operations. The water quality laboratory at the Maricopa Agricultural Center (in collaboration with stakeholder groups and individuals located in the ADEQ-targeted watersheds) identified biological changes in surface water quality in impaired reaches rivers to determine the sources of Escherichia coli (E. coli); bacteria utilizing innovative water quality microbial/bacterial source tracking methods. The end goal was to support targeted watershed groups and ADEQ towards E. coli reductions. In the field monitoring was conducted by the selected targeted watershed groups in conjunction with The University of Arizona Maricopa Agricultural Center Water Quality Laboratory. This consisted of collecting samples for Bacteroides testing from multiple locations on select impaired reaches, to determine contamination resulting from cattle, human recreation, and other contributions. Such testing was performed in conjunction with high flow and base flow conditions in order to accurately portray water quality conditions and variations. Microbial monitoring was conducted by The University of Arizona Water Quality Laboratory at the Maricopa Agricultural Center using genetic typing to differentiate among two categories of Bacteroides: human and all (total). Testing used microbial detection methodologies and molecular source tracking techniques.^

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.