Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sequence Alignment”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

ModBase, a database of annotated comparative protein structure models and associated resources

ModBase (http://salilab.org/modbase) is a database of annotated comparative protein structure models. The models are calculated by ModPipe, an automated modeling pipeline that relies primarily on Modeller for fold assignment, sequence-structure alignment, model building and model assessment (http://salilab.org/modeller/). ModBase currently contains almost 30 million reliable models for domains in 4.7 million unique protein sequences. ModBase allows users to compute or update comparative models on demand, through an interface to the ModWeb modeling server (http://salilab.org/ modweb). ModBase models are also available through the Protein Model Portal (http://www.prote inmodelportal.org/). Recently developed associated resources include the AllosMod server for modeling ligand-induced protein dynamics (http://salilab.org/ allosmod), the AllosMod-FoXS server for predicting a structural ensemble that fits an SAXS profile (http://salilab.org/allosmod-foxs), the FoXSDock server for protein–protein docking filtered by an SAXS profile (http://salilab.org/foxsdock), the SAXS Merge server for automatic merging of SAXS profiles (http://salilab.org/saxsmerge) and the Pose & Rank server for scoring protein–ligand complexes (http://salilab.org/poseandrank). In this update, we also highlight two applications of ModBase: a PSI:Biology initiative to maximize the structural coverage of the human alpha-helical transmembrane proteome and a determination of structural determinants of human immunodeficiency virus-1 protease specificity.

59 BASIC BIOLOGICAL SCIENCES↗

Multiplex detection and identification of viral, bacterial, and protozoan pathogens in human blood and plasma using an expanded high-density resequencing microarray platform

Introduction: Nucleic acid tests for blood donor screening have improved the safety of the blood supply; however, increasing numbers of emerging pathogen tests are burdensome. Multiplex testing platforms are a potential solution. Methods: The Blood Borne Pathogen Resequencing Microarray Expanded (BBP-RMAv.2) can perform multiplex detection and identification of 80 viruses, bacteria and parasites. This study evaluated pathogen detection in human blood or plasma. Samples spiked with selected pathogens, each with one of 6 viruses, 2 bacteria and 5 protozoans were tested on this platform. The nucleic acids were extracted, amplified using multiplexed sets of primers, and hybridized to a microarray. The reported sequences were aligned to a database to identify the pathogen. To directly compare the microarray to an emerging molecular approach, the amplified nucleic acids were also submitted to nanopore next generation sequencing (NGS). Results: The BBP-RMAv.2 detected viral pathogens at a concentration as low as 100 copies/ml and a range of concentrations from 1,000 to 100,000 copies/ml for all the spiked pathogens. Coded specimens were identified correctly demonstrating the effectiveness of the platform. The nanopore sequencing correctly identified most samples and the results of the two platforms were compared. Discussion: These results indicated that the BBP-RMAv.2 could be employed for multiplex detection with potential for use in blood safety or disease diagnosis. The NGS was nearly as effective at identifying pathogens in blood and performed better than BBP-RMAv.2 at identifying pathogen-negative samples.

59 BASIC BIOLOGICAL SCIENCES↗

Structure, function, and evolution of bacterial ATP-binding cassette systems

The ATP-binding cassette (ABC) systems constitute one of the largest superfamilies of paralogous sequences. All ABC systems share a highly conserved ATP-hydrolyzing domain or protein (the ABC; also referred to as a nucleotide-binding domain [NBD]) that is unequivocally characterized by three short sequence motifs (Fig. 1): these are the Walker A and Walker B motifs, indicative of the presence of a nucleotide-binding site, and the signature motif, unique to ABC proteins, located upstream of the Walker B motif (426). Other motifs diagnostic of ABC proteins are also indicated in Fig. 1. The biological significance of these motifs is discussed in Structure, Function, and Dynamics of the ABC. ABC systems are widespread among living organisms and have been detected in all genera of the three kingdoms of life, with remarkable conservation in the primary sequence of the cassette and in the organization of the constitutive domains or subunits (203, 420). ABC systems couple the energy of ATP hydrolysis to an impressively large variety of essential biological phenomena, comprising not only transmembrane (TM) transport, for which they are best known, but also several non-transport-related processes, such as translation elongation (62) and DNA repair (174). Although ABC systems deserve much attention because they are involved in severe human inherited diseases (107), they were first discovered and characterized in detail in prokaryotes, as early as the 1970s (13, 148, 238, 468). The most extensively analyzed systems were the high-affinity histidine and maltose uptake systems of Salmonella enterica serovar Typhimurium and Escherichia coli. Over 2 decades ago, after the completion of the nucleotide sequences encoding these transporters in the respective laboratories of Giovanna Ames and Maurice Hofnung, Hiroshi Nikaido and colleagues noticed that the two systems displayed a global similarity in the nature of their components and, moreover, that the primary sequences of MalK and HisP, the proteins suspected to energize these transporters, shared as much as 32% identity in amino acid residues when their sequences were aligned (171). Later, it was found that several bacterial proteins involved in uptake of nutrients, export of toxins, cell division, bacterial nodulation of plants, and DNA repair displayed the same similarity in their sequences (127, 196). This led to the notion that the conserved protein, which had been shown to bind ATP (198, 201), would probably energize the systems mentioned above by coupling the energy of ATP hydrolysis to transport. The latter was demonstrated with the maltose and histidine transporters by use of isolated membrane vesicles (105, 379) and purified transporters reconstituted into proteoliposomes (30, 98). The determination of the sequence of the first eukaryotic protein strongly similar to these bacterial transporters (the P-glycoprotein, involved in resistance of cancer cells to multiple drugs) (169, 179) demonstrated that these proteins were not restricted to prokaryotes. Two names, 'traffic ATPases' (15) and the more accepted name 'ABC transporters' (193, 218), were proposed for members of this new superfamily. ABC systems can be divided into three main functional categories, as follows. Importers mediate the uptake of nutrients in prokaryotes. The nature of the substrates that are transported is very wide, including mono- and oligosaccharides, organic and inorganic ions, amino acids, peptides, ironsiderophores, metals, polyamine cations, opines, and vitamins. Exporters are involved in the secretion of various molecules, such as peptides, lipids, hydrophobic drugs, polysaccharides, and proteins, including toxins such as hemolysin. The third category of systems is apparently not involved in transport, with some members being involved in translation of mRNA and in DNA repair. Despite the large, diverse population of substrates handled and the difference in the polarity of transport, importers and exporters share a common organization made of two hydrophobic membrane-spanning or integral membrane (IM) domains and two hydrophilic domains carrying the ABC peripherally associated with the IM domains on the cytosolic side of the membrane (26). In importers, these four domains are almost always independent polypeptide chains that come together to form a multimeric complex. In most exporters, including the E. coli hemolysin exporter HlyB, the N-terminal IM and the C-terminal ABC domains are fused as a single polypeptide chain (IM-ABC). An inverted organization in which the IM domain is C-terminal with respect to the ABC domain (ABC-IM) exists, such as in the MacB protein, involved in macrolide resistance in E. coli. No IM domain partners have been identified for ABC proteins falling into the third category, and these proteins consist of two ABCs fused together (ABC2).

59 BASIC BIOLOGICAL SCIENCES↗

Scalable Parallel Methods for Analyzing Metagenomics Data at Extreme Scale

The field of bioinformatics and computational biology is currently experiencing a data revolution. The exciting prospect of making fundamental biological discoveries is fueling the rapid development and deployment of numerous cost-effective, high-throughput next-generation sequencing technologies. The result is that the DNA and protein sequence repositories are being bombarded with new sequence information. Databases are continuing to report a Moore’s law-like growth trajectory in their database sizes, roughly doubling every 18 months. In what seems to be a paradigm-shift, individual projects are now capable of generating billions of raw sequence data that need to be analyzed in the presence of already annotated sequence information. While it is clear that data-driven methods, such as sequencing homology detection, are becoming the mainstay in the field of computational life sciences, the algorithmic advancements essential for implementing complex data analytics at scale have mostly lagged behind. Sequence homology detection is central to a number of bioinformatics applications including genome sequencing and protein family characterization. Given millions of sequences, the goal is to identify all pairs of sequences that are highly similar (or “homologous”) on the basis of alignment criteria. While there are optimal alignment algorithms to compute pairwise homology, their deployment for large-scale is currently not feasible; instead, heuristic methods are used at the expense of quality. In this dissertation, we present the design and evaluation of a parallel implementation for conducting optimal homology detection on distributed memory supercomputers. Our approach uses a combination of techniques from asynchronous load balancing (viz. work stealing, dynamic task counters), data replication, and exact-matching filters to achieve homology detection at scale. Results for a collection of 2.56M sequences show parallel efficiencies of ~75-100% on up to 8K cores, representing a time-to-solution of 33 seconds. We extend this work with a detailed analysis of single-node sequence alignment performance using the latest CPU vector instruction set extensions. Preliminary results reveal that current sequence alignment algorithms are unable to fully utilize widening vector registers.

59 BASIC BIOLOGICAL SCIENCES↗

Classification of bacterial plasmid and chromosome derived sequences using machine learning

Plasmids are important genetic elements that facilitate horizonal gene transfer between bacteria and contribute to the spread of virulence and antimicrobial resistance. Most bacterial genome sequences in the public archives exist in draft form with many contigs, making it difficult to determine if a contig is of chromosomal or plasmid origin. Using a training set of contigs comprising 10,584 chromosomes and 10,654 plasmids from the PATRIC database, we evaluated several machine learning models including random forest, logistic regression, XGBoost, and a neural network for their ability to classify chromosomal and plasmid sequences using nucleotide k-mers as features. Based on the methods tested, a neural network model that used nucleotide 6-mers as features that was trained on randomly selected chromosomal and plasmid subsequences 5kb in length achieved the best performance, outperforming existing out-of-the-box methods, with an average accuracy of 89.38% ± 2.16% over a 10-fold cross validation. The model accuracy can be improved to 92.08% by using a voting strategy when classifying holdout sequences. In both plasmids and chromosomes, subsequences encoding functions involved in horizontal gene transfer—including hypothetical proteins, transporters, phage, mobile elements, and CRISPR elements—were most likely to be misclassified by the model. This study provides a straightforward approach for identifying plasmid-encoding sequences in short read assemblies without the need for sequence alignment-based tools.

59 BASIC BIOLOGICAL SCIENCES↗

SCOPe: improvements to the structural classification of proteins – extended database to facilitate variant interpretation and machine learning

Abstract The Structural Classification of Proteins—extended (SCOPe, https://scop.berkeley.edu) knowledgebase aims to provide an accurate, detailed, and comprehensive description of the structural and evolutionary relationships amongst the majority of proteins of known structure, along with resources for analyzing the protein structures and their sequences. Structures from the PDB are divided into domains and classified using a combination of manual curation and highly precise automated methods. In the current release of SCOPe, 2.08, we have developed search and display tools for analysis of genetic variants we mapped to structures classified in SCOPe. In order to improve the utility of SCOPe to automated methods such as deep learning classifiers that rely on multiple alignment of sequences of homologous proteins, we have introduced new machine-parseable annotations that indicate aberrant structures as well as domains that are distinguished by a smaller repeat unit. We also classified structures from 74 of the largest Pfam families not previously classified in SCOPe, and we improved our algorithm to remove N- and C-terminal cloning, expression and purification sequences from SCOPe domains. SCOPe 2.08-stable classifies 106 976 PDB entries (about 60% of PDB entries).

59 BASIC BIOLOGICAL SCIENCES↗

Data from: Generation of transcript assemblies and identification of single nucleotide polymorphisms from seven lowland and upland cultivars of switchgrass

Switchgrass (Panicum virgatum) is a North American perennial prairie species that has been used as a rangeland and forage crop and has recently been targeted as a potential biofuel feedstock species. Switchgrass, which occurs as tetraploid and octoploid forms, is classified into lowland or upland ecotypes that differ in growth phenotypes and adaptation to distinct habitats. Using RNA-sequencing reads derived from crown, young shoot and leaf tissues, we generated sequence data from seven switchgrass cultivars, three lowland and four upland, to enable comparative analyses between switchgrass cultivars and to identify single nucleotide polymorphisms (SNPs) for use in breeding and genetic analysis. We also generated individual transcript assemblies for each of the cultivars. Transcript data indicate that subgenomes of octoploid switchgrass are not substantially different from subgenomes of tetraploids as expected for an autopolyploid origin of switchgrass octoploids. Using RNA-sequencing reads aligned to the switchgrass Release 0 AP13 reference genome, we identified 1,305,976 high-confidence SNPs. Of these SNPs, 438,464 were unique to lowland cultivars, but only 12,002 were found in all lowlands. Conversely, 723,678 SNPs were unique to upland cultivars with only 34,665 observed in all uplands. Comparison of our high-confidence transcriptome-derived SNPs with SNPs previously identified in a genotyping-by-sequencing study of an association panel revealed limited overlap between the two methods, highlighting the utility of transcriptome-based SNP discovery in augmenting genome diversity polymorphism datasets. The transcript and SNP data described here provide a useful resource for switchgrass gene annotation and marker-based analyses of the switchgrass genome.

09 BIOMASS FUELS↗

Analysis of genomic signatures associated with Variovorax endosphere colonization

This repository contains the analysis code and supporting datasets associated with the study “Genomic signatures in Variovorax enabling colonization of the Populus endosphere.” Beals DG, Carper DL, Hochanadel LH, Jawdy SS, Klingeman DM, Piatkowski BT, Weston DJ, Doktycz MJ, Pelletier DA. 2026. Genomic signatures in Variovorax enabling colonization of the Populus endosphere. mSystems 11:e01605-25. https://doi.org/10.1128/msystems.01605-25 The scripts are organized sequentially (01–07) and document the workflows used for: Sequence-read alignment and feature counting Orthogroup and KEGG Ortholog annotation Count normalization Statistical analysis and aggregation Generation of manuscript figures and tables Repository contents The uncompressed files are the finalized, formatted datasets used to generate the figures and tables reported in the study, including the supplemental CSV files referenced in the manuscript. The accompanying ZIP archive contains the complete codebase and example data_input/ and data_output/ directories illustrating the organization and execution of the analytical workflow. Individual scripts identify the corresponding manuscript analyses and figure panels. Raw sequencing data Raw sequencing reads are available through the NCBI Sequence Read Archive under BioProject accession PRJNA1322484.

Beals, Delaney [ORNL] (ORCID:0000000306274574)↗

Structures of TOG1 and TOG2 from the human microtubule dynamics regulator CLASP1

Tubulin-binding TOG domains are found arrayed in a number of proteins that regulate microtubule dynamics. While much is known about the structure and function of TOG domains from the XMAP215 microtubule polymerase family, less in known about the TOG domain array found in animal CLASP family members. The animal CLASP TOG array promotes microtubule pause, potentiates rescue, and limits catastrophe. How structurally distinct the TOG domains of animal CLASP are from one another, from XMAP215 family TOG domains, and whether a specific order of structurally distinct TOG domains in the TOG array is conserved across animal CLASP family members is poorly understood. We present the x-ray crystal structures of Homo sapiens (H.s.) CLASP1 TOG1 and TOG2. The structures of H.s. CLASP1 TOG1 and TOG2 are distinct from each other and from the previously determined structure of Mus musculus (M.m.) CLASP2 TOG3. Comparative analyses of CLASP family TOG domain structures determined to date across species and paralogs supports a conserved CLASP TOG array paradigm in which structurally distinct TOG domains are arrayed in a specific order. H.s. CLASP1 TOG1 bears structural similarity to the free-tubulin binding TOG domains of the XMAP215 family but lacks many of the key tubulin-binding determinants found in XMAP215 family TOG domains. This aligns with studies that report that animal CLASP family TOG1 domains cannot bind free tubulin or microtubules. In contrast, animal CLASP family TOG2 and TOG3 domains have reported microtubule-binding activity but are structurally distinct from the free-tubulin binding TOG domains of the XMAP215 family. H.s. CLASP1 TOG2 has a convex architecture, predicted to engage a hyper-curved tubulin state that may underlie its ability to limit microtubule catastrophe and promote rescue. M.m. CLASP2 TOG3 has unique structural elements in the C-terminal half of its α-solenoid domain that our modeling studies implicate in binding to laterally-associated tubulin subunits in the microtubule lattice in a mode similar to, yet distinct from those predicted for the XMAP215 family TOG4 domain. The potential ability of the animal CLASP family TOG3 domain to engage lateral tubulin subunits may underlie the microtubule rescue activity ascribed to the domain. These findings highlight the structural diversity of TOG domains within the CLASP family TOG array and provide a molecular foundation for understanding CLASP-dependent effects on microtubule dynamics.

59 BASIC BIOLOGICAL SCIENCES↗

Standardized Residue Numbering and Secondary Structure Nomenclature in the Class D β-Lactamases

Over 1370 class D β-lactamases are currently known, and they pose a serious threat to the effective treatment of many infectious diseases, particularly in some pathogenic bacteria where evolving carbapenemase activity has been reported. Detailed understanding of their molecular biology, enzymology, and structural biology are critically important, but the lack of a standardized residue numbering scheme and inconsistent secondary structure annotation has made comparative analyses sometimes difficult and cumbersome. Compounding this, in the post-AlphaFold world where we currently find ourselves, an extraordinary wealth of detailed structural information on these enzymes is literally at our fingertips; therefore it is vitally important that a standard numbering system is in place to facilitate the accurate and straightforward analysis of their structures. In conclusion, here we present a residue numbering and secondary structure scheme for the class D enzymes based on the sequence and structure of OXA-48 and apply it to test targets to demonstrate the ease with which it can be used.

59 BASIC BIOLOGICAL SCIENCES↗

Proteome-wide analysis of HIV-specific naive and memory CD4 + T cells in unexposed blood donors

The preexisting HIV-1–specific T cell repertoire must influence both the immunodominance of T cells after infection and immunogenicity of vaccines. We directly compared two methods for measuring the preexisting CD4 + T cell repertoire in healthy HIV-1–negative volunteers, the HLA-peptide tetramer enrichment and T cell library technique, and show high concordance (r = 0.989). Using the library technique, we examined whether naive, central memory, and/or effector memory CD4 + T cells specific for overlapping peptides spanning the entire HIV-1 proteome were detectable in 10 HLA diverse, HIV-1–unexposed, seronegative donors. HIV-1–specific cells were detected in all donors at a mean of 55 cells/million naive cells and 38.9 and 34.1 cells/million in central and effector memory subsets. Remarkably, peptide mapping showed most epitopes recognized by naive (88%) and memory (56%) CD4 + T cells had been previously reported in natural HIV-1 infection. Furthermore, 83% of epitopes identified in preexisting memory subsets shared epitope length matches (8–12 amino acids) with human microbiome proteins, suggestive of a possible cross-reactive mechanism. These results underline the power of a proteome-wide analysis of peptide recognition by human T cells for the identification of dominant antigens and provide a baseline for optimizing HIV-1–specific helper cell responses by vaccination.

59 BASIC BIOLOGICAL SCIENCES↗

The repABC plasmids with quorum-regulated transfer systems in members of the Rhizobiales divide into two structurally and separately evolving groups

The large repABC plasmids of the order Rhizobiales with Class I quorum-regulated conjugative transfer systems often define the nature of the bacterium that harbors them. These otherwise diverse plasmids contain a core of highly conserved genes for replication and conjugation raising the question of their evolutionary relationships. In an analysis of 18 such plasmids these elements fall into two organizational classes, Group I and Group II, based on the sites at which cargo DNA is located. Cladograms constructed from proteins of the transfer and quorum-sensing components indicated that those of the Group I plasmids, while coevolving, have diverged from those coevolving proteins of the Group II plasmids. Moreover, within these groups the phylogenies of the proteins usually occupy similar, if not identical, tree topologies. Remarkably, such relationships were not seen among proteins of the replication system; although RepA and RepB coevolve, RepC does not. Nor do the replication proteins coevolve with the proteins of the transfer and quorum-sensing systems. Functional analysis was mostly consistent with phylogenies. TraR activated promoters from plasmids within its group, but not between groups and dimerized with TraR proteins from within but not between groups. However, oriT sequences, which are highly conserved, were processed by the transfer system of plasmids regardless of group. Here, we conclude that these plasmids diverged into two classes based on the locations at which cargo DNA is inserted, that the quorum-sensing and transfer functions are coevolving within but not between the two groups, and that this divergent evolution extends to function.

59 BASIC BIOLOGICAL SCIENCES↗

Tandem repeats in giant archaeal Borg elements undergo rapid evolution and create new intrinsically disordered regions in proteins

Borgs are huge, linear extrachromosomal elements associated with anaerobic methane-oxidizing archaea. Striking features of Borg genomes are pervasive tandem direct repeat (TR) regions. Here, we present six new Borg genomes and investigate the characteristics of TRs in all ten complete Borg genomes. We find that TR regions are rapidly evolving, recently formed, arise independently, and are virtually absent in host Methanoperedens genomes. Flanking partial repeats and A-enriched character constrain the TR formation mechanism. TRs can be in intergenic regions, where they might serve as regulatory RNAs, or in open reading frames (ORFs). TRs in ORFs are under very strong selective pressure, leading to perfect amino acid TRs (aaTRs) that are commonly intrinsically disordered regions. Proteins with aaTRs are often extracellular or membrane proteins, and functionally similar or homologous proteins often have aaTRs composed of the same amino acids. We propose that Borg aaTR-proteins functionally diversify Methanoperedens and all TRs are crucial for specific Borg–host associations and possibly cospeciation.

59 BASIC BIOLOGICAL SCIENCES↗

A taxonomy of bacterial microcompartment loci constructed by a novel scoring method

Bacterial microcompartments (BMCs) are proteinaceous organelles involved in both autotrophic and heterotrophic metabolism. All BMCs share homologous shell proteins but differ in their complement of enzymes; these are typically encoded adjacent to shell protein genes in genetic loci, or operons. To enable the identification and prediction of functional (sub)types of BMCs, we developed LoClass, an algorithm that finds putative BMC loci and inventories, weights, and compares their constituent pfam domains to construct a locus similarity network and predict locus (sub)types. In addition to using LoClass to analyze sequences in the Non-redundant Protein Database, we compared predicted BMC loci found in seven candidate bacterial phyla (six from single-cell genomic studies) to the LoClass taxonomy. Together, these analyses resulted in the identification of 23 different types of BMCs encoded in 30 distinct locus (sub)types found in 23 bacterial phyla. These include the two carboxysome types and a divergent set of metabolosomes, BMCs that share a common catalytic core and process distinct substrates via specific signature enzymes. Furthermore, many Candidate BMCs were found that lack one or more core metabolosome components, including one that is predicted to represent an entirely new paradigm for BMC-associated metabolism, joining the carboxysome and metabolosome. By placing these results in a phylogenetic context, we provide a framework for understanding the horizontal transfer of these loci, a starting point for studies aimed at understanding the evolution of BMCs. This comprehensive taxonomy of BMC loci, based on their constituent protein domains, foregrounds the functional diversity of BMCs and provides a reference for interpreting the role of BMC gene clusters encoded in isolate, single cell, and metagenomic data. Many loci encode ancillary functions such as transporters or genes for cofactor assembly; this expanded vocabulary of BMC-related functions should be useful for design of genetic modules for introducing BMCs in bioengineering applications.

59 BASIC BIOLOGICAL SCIENCES↗

Contact-dependent growth inhibition (CDI) systems deploy a large family of polymorphic ionophoric toxins for inter-bacterial competition

Contact-dependent growth inhibition (CDI) is a widespread form of inter-bacterial competition mediated by CdiA effector proteins. CdiA is presented on the inhibitor cell surface and delivers its toxic C-terminal region (CdiA-CT) into neighboring bacteria upon contact. Inhibitor cells also produce CdiI immunity proteins, which neutralize CdiA-CT toxins to prevent auto-inhibition. Here, we describe a diverse group of CDI ionophore toxins that dissipate the transmembrane potential in target bacteria. These CdiA-CT toxins are composed of two distinct domains based on AlphaFold2 modeling. The C-terminal ionophore domains are all predicted to form five-helix bundles capable of spanning the cell membrane. The N-terminal "entry" domains are variable in structure and appear to hijack different integral membrane proteins to promote toxin assembly into the lipid bilayer. The CDI ionophores deployed by E. coli isolates partition into six major groups based on their entry domain structures. Comparative sequence analyses led to the identification of receptor proteins for ionophore toxins from groups 1 & 3 (AcrB), group 2 (SecY) and groups 4 (YciB). Using forward genetic approaches, we identify novel receptors for the group 5 and 6 ionophores. Group 5 exploits homologous putrescine import proteins encoded by puuP and plaP, and group 6 toxins recognize di/tripeptide transporters encoded by paralogous dtpA and dtpB genes. Finally, we find that the ionophore domains exhibit significant intra-group sequence variation, particularly at positions that are predicted to interact with CdiI. Accordingly, the corresponding immunity proteins are also highly polymorphic, typically sharing only ~30% sequence identity with members of the same group. Competition experiments confirm that the immunity proteins are specific for their cognate ionophores and provide no protection against other toxins from the same group. The specificity of this protein interaction network provides a mechanism for self/nonself discrimination between E. coli isolates.

59 BASIC BIOLOGICAL SCIENCES↗

Genetic Variability of West Nile Virus in U.S. Blood Donors from the 2012 Epidemic Season

West Nile virus (WNV) is an arbovirus maintained in nature in a bird-mosquito enzootic cycle which can also infect other vertebrates including humans. WNV is now endemic in the United States (U.S.), causing yearly outbreaks that have resulted in an estimated total of 4–5 million human infections. Over 41,700 cases of West Nile disease, including 18,810 neuroinvasive cases and 1,765 deaths, were reported to the CDC between 1999 and 2014. In 2012, the second largest West Nile outbreak in the U.S. was reported, which caused 5,674 cases and 286 deaths. WNV continues to evolve, and three major WNV lineage I genotypes (NY99, WN02, and SW/WN03) have been described in the U.S. since introduction of the virus in 1999. We report here the WNV sequences obtained from 19 human samples acquired during the 2012 U.S. outbreak and our examination of the evolutionary dynamics in WNV isolates sequenced from 1999–2012. Maximum-likelihood and Bayesian methods were used to perform the phylogenetic analyses. Selection pressure analyses were performed with the HyPhy package using the Datamonkey web-server. Using different codon-based and branch-site selection models, we detected a number of codons subjected to positive pressure in WNV genes. Thirteen of the 19 completely sequenced isolates from 10 U.S. states were genetically similar, sharing up to 55 nucleotide mutations and 4 amino acid substitutions when compared with the prototype isolate WN-NY99. Altogether, these analyses showed that following a brief contraction in 2008–2009, WNV genetic divergence in the U.S. continued to increase in 2012, and that closely related variants were found across a broad geographic range of the U.S., coincident with the second-largest WNV outbreak in U.S. history.

59 BASIC BIOLOGICAL SCIENCES↗

Crystal Structure of the 6-Hydroxymethyl-7,8-Dihydropterin Pyrophosphokinase•Dihydropteroate Synthase Bifunctional Enzyme from Francisella tularensis

The 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) and dihydropteroate synthase (DHPS) enzymes catalyze sequential metabolic reactions in the folate biosynthetic pathway of bacteria and lower eukaryotes. Both enzymes represent validated targets for the development of novel anti-microbial therapies. We report herein that the genes which encode Ft HPPK and Ft DHPS from the biowarfare agent Francisella tularensis are fused into a single polypeptide. The potential of simultaneously targeting both modules with pterin binding inhibitors prompted us to characterize the molecular details of the multifunctional complex. Our high resolution crystallographic analyses reveal the structural organization between Ft HPPK and Ft DHPS which are tethered together by a short linker. Additional structural analyses of substrate complexes reveal that the active sites of each module are virtually indistinguishable from those of the monofunctional enzymes. The fused bifunctional enzyme therefore represents an excellent vehicle for finding inhibitors that engage the pterin binding pockets of both modules that have entirely different architectures. To demonstrate that this approach has the potential of producing novel two-hit inhibitors of the folate pathway, we identify and structurally characterize a fragment-like molecule that simultaneously engages both active sites. Our study provides a molecular framework to study the enzyme mechanisms of HPPK and DHPS, and to design novel and much needed therapeutic compounds to treat infectious diseases.

59 BASIC BIOLOGICAL SCIENCES↗

Viral Linkage in HIV-1 Seroconverters and Their Partners in an HIV-1 Prevention Clinical Trial

Characterization of viruses in HIV-1 transmission pairs will help identify biological determinants of infectiousness and evaluate candidate interventions to reduce transmission. Although HIV-1 sequencing is frequently used to substantiate linkage between newly HIV-1 infected individuals and their sexual partners in epidemiologic and forensic studies, viral sequencing is seldom applied in HIV-1 prevention trials. The Partners in Prevention HSV/HIV Transmission Study (ClinicalTrials.gov #NCT00194519) was a prospective randomized placebo-controlled trial that enrolled serodiscordant heterosexual couples to determine the efficacy of genital herpes suppression in reducing HIV-1 transmission; as part of the study analysis, HIV-1 sequences were examined for genetic linkage between seroconverters and their enrolled partners. We obtained partial consensus HIV-1 env and gag sequences from blood plasma for 151 transmission pairs and performed deep sequencing of env in some cases. We analyzed sequences with phylogenetic techniques and developed a Bayesian algorithm to evaluate the probability of linkage. For linkage, we required monophyletic clustering between enrolled partners’ sequences and a Bayesian posterior probability of $50%. Adjudicators classified each seroconversion, finding 108 (71.5%) linked, 40 (26.5%) unlinked, and 3 (2.0%) indeterminate transmissions, with linkage determined by consensus env sequencing in 91 (84%). Male seroconverters had a higher frequency of unlinked transmissions than female seroconverters. The likelihood of transmission from the enrolled partner was related to time on study, with increasing numbers of unlinked transmissions occurring after longer observation periods. Finally, baseline viral load was found to be significantly higher among linked transmitters. In this first use of HIV-1 sequencing to establish endpoints in a large clinical trial, more than one-fourth of transmissions were unlinked to the enrolled partner, illustrating the relevance of these methods in the design of future HIV-1 prevention trials in serodiscordant couples. A hierarchy of sequencing techniques, analysis methods, and expert adjudication contributed to the linkage determination process.

59 BASIC BIOLOGICAL SCIENCES↗