Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Homology modelling”

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

Structural Models and Sequence Alignment Results of the Desulfovibrio vulgaris Proteome

This dataset contains the structural models for the primary transcripts of the Desulfovibrio vulgaris proteome as well as sequence alignment results for a subset of the encoded proteins. For each protein, the five models inferred from AlphaFold 2 are provided. The largest pTM-scoring model for each protein was energy minimized; this minimized structure as well as its AlphaFold pickle output file are also provided. This set of structures represent an alternate source of models for the D. vulgaris proteome to those available in the AlphaFold Protein Structure Database (AFDB). This is a bit more complicated since the proteins reporting in the AFDB originate from an outdated form of the D. vulgaris sequence. The different versions of the D. vulgaris gene annotation are collected in the Chronology subdirectory; further consideration of these changes on the structural space of the proteome are currently underway. For proteins that have been annotated as hypothetical, sequence alignment results from the HHblits and SAdLSA alignment methods are provided. These methods are often more capable to resolve sequence homology than other methods. Therefore, the results from both HHblits and SAdLSA are provided to identify possible homologs for these challenging proteins. Numerous sequence databases are utilized for these alignments. References AlphaFold v2 Multimer: https://doi.org/10.1101/2021.10.04.463034. References HHblits: hhtps://doi.org/10.1186/s12859-019-3019-7. References SAdLSA: hhtps://doi.org/10.3389/fbinf.2021.689960.

59 BASIC BIOLOGICAL SCIENCES↗

cGAS-like receptors sense RNA and control 3'2'-cGAMP signalling in Drosophila

Cyclic GMP–AMP synthase (cGAS) is a cytosolic DNA sensor that produces the second messenger cG[2'–5']pA[3'–5']p (2'3'-cGAMP) and controls activation of innate immunity in mammalian cells. Animal genomes typically encode multiple proteins with predicted homology to cGAS, but the function of these uncharacterized enzymes is unknown. Here we show that cGAS-like receptors (cGLRs) are innate immune sensors that are capable of recognizing divergent molecular patterns and catalysing synthesis of distinct nucleotide second messenger signals. Crystal structures of human and insect cGLRs reveal a nucleotidyltransferase signalling core shared with cGAS and a diversified primary ligand-binding surface modified with notable insertions and deletions. We demonstrate that surface remodelling of cGLRs enables altered ligand specificity and used a forward biochemical screen to identify cGLR1 as a double-stranded RNA sensor in the model organism Drosophila melanogaster. We show that RNA recognition activates Drosophila cGLR1 to synthesize the novel product cG[3'–5']pA[2'–5']p (3'2'-cGAMP). A crystal structure of Drosophila stimulator of interferon genes (dSTING) in complex with 3'2'-cGAMP explains selective isomer recognition, and 3'2'-cGAMP induces an enhanced antiviral state in vivo that protects from viral infection. Similar to radiation of Toll-like receptors in pathogen immunity, our results establish cGLRs as a diverse family of metazoan pattern recognition receptors.

59 BASIC BIOLOGICAL SCIENCES↗

Chromosome-Level Genome Assembly of Mentha longifolia L. Reveals Gene Organization Underlying Disease Resistance and Essential Oil Traits

Mentha longifolia (L.) Huds., a wild, diploid mint species, has been developed as a model for mint genetic and genomic research to aid breeding efforts that target Verticillium wilt disease resistance and essential oil monoterpene composition. Here, we present a near-complete, chromosome-scale mint genome assembly for M. longifolia USDA accession CMEN 585. This new assembly is an update of a previously published genome draft, with dramatic improvements. A total of 42,107 protein-coding genes were annotated and placed on 12 chromosomal scaffolds. One hundred fifty-three genes contained conserved sequence domains consistent with nucleotide binding site-leucine-rich-repeat (NBS-LRR) plant disease resistance genes. Homologs of genes implicated in Verticillium wilt resistance in other plant species were also identified. Multiple paralogs of genes putatively involved in p-menthane monoterpenoid biosynthesis were identified and several cases of gene clustering documented. Heterologous expression of candidate genes, purification of recombinant target proteins, and subsequent enzyme assays allowed us to identify the genes underlying the pathway that leads to the most abundant monoterpenoid volatiles. Here, the bioinformatic and functional analyses presented here are laying the groundwork for using marker-assisted selection in improving disease resistance and essential oil traits in mints.

59 BASIC BIOLOGICAL SCIENCES↗

The catalytic domains of Streptococcus mutans glucosyltransferases: a structural analysis

Streptococcus mutans, found in the human oral cavity, is a significant contributor to the pathogenesis of dental caries. This bacterium expresses three genetically distinct types of glucosyltransferases named GtfB (GTF-I), GtfC (GTF-SI) and GtfD (GTF-S) that play critical roles in the development of dental plaque. The catalytic domains of GtfB, GtfC and GtfD contain conserved active-site residues for the overall enzymatic activity that relate to hydrolytic glycosidic cleavage of sucrose to glucose and fructose, release of fructose and generation of a glycosyl-enzyme intermediate in the reducing end. In a subsequent transglycosylation step, the glucosyl moiety is transferred to the nonreducing end of an acceptor to form a growing glucan polymer chain made up of glucose molecules. It has been proposed that both sucrose breakdown and glucan synthesis occur in the same active site of the catalytic domain, although the active site does not appear to be large enough to accommodate both functions. These three enzymes belong to glycoside hydrolase family 70 (GH70), which shows homology to glycoside hydrolase family 13 (GH13). GtfC synthesizes both soluble and insoluble glucans (α-1,3 and α-1,6 glycosidic linkages), while GtfB and GtfD synthesize only insoluble or soluble glucans, respectively. Here, crystal structures of the catalytic domains of GtfB and GtfD are reported. These structures are compared with previously determined structures of the catalytic domain of GtfC. With this work, apo structures and inhibitor-complex structures with acarbose are now available for the catalytic domains of GtfC and GtfB. The structure of GtfC with maltose allows further identification and comparison of active-site residues. A model of sucrose binding to GtfB is also included. The new structure of the catalytic domain of GtfD affords a structural comparison of the three S. mutans glycosyltransferases. Unfortunately, the catalytic domain of GtfD is not complete since crystallization resulted in the structure of a truncated protein lacking approximately 200 N-terminal residues of domain IV.

59 BASIC BIOLOGICAL SCIENCES↗

TopFusion: Using Topological Feature Space for Fusion and Imputation in Multi-Modal Data

We present a novel multi-modal data fusion technique using topological features. The method, TopFusion, leverages the flexibility of topological data analysis tools (namely persistent homology and persistence images) to map multi-modal datasets into a common feature space by forming a new multi-channel persistence image. Each channel in the image is representative of a view of the data from a modality-dependent filtration. We demonstrate that the topological perspective we take allows for more effective data reconstruction, i.e. imputation. In particular, by performing imputation in topological feature space we are able to outperform the same imputation techniques applied to raw data or alternatively derived features. We show that TopFusion representations can be used as input to downstream deep learning-based computer vision models and doing so achieves comparable performance to other fusion methods for classification on two multi-modal datasets.

Myers, Audun D.↗

Functional anatomy of zinc finger antiviral protein complexes

Abstract ZAP is an antiviral protein that binds to and depletes viral RNA, which is often distinguished from vertebrate host RNA by its elevated CpG content. Two ZAP cofactors, TRIM25 and KHNYN, have activities that are poorly understood. Here, we show that functional interactions between ZAP, TRIM25 and KHNYN involve multiple domains of each protein, and that the ability of TRIM25 to multimerize via its RING domain augments ZAP activity and specificity. We show that KHNYN is an active nuclease that acts in a partly redundant manner with its homolog N4BP1. The ZAP N-terminal RNA binding domain constitutes a minimal core that is essential for antiviral complex activity, and we present a crystal structure of this domain that reveals contacts with the functionally required KHNYN C-terminal domain. These contacts are remote from the ZAP CpG binding site and would not interfere with RNA binding. Based on our dissection of ZAP, TRIM25 and KHNYN functional anatomy, we could design artificial chimeric antiviral proteins that reconstitute the antiviral function of the intact authentic proteins, but in the absence of protein domains that are otherwise required for activity. Together, these results suggest a model for the RNA recognition and action of ZAP-containing antiviral protein complexes.

Science & Technology - Other Topics↗

Anaerobic Fungal Mevalonate Pathway Genomic Biases Lead to Heterologous Toxicity Underpredicted by Codon Adaptation Indices

Anaerobic fungi are emerging biotechnology platforms with genomes rich in biosynthetic potential. Yet, the heterologous expression of their biosynthetic pathways has had limited success in model hosts like E. coli. We find one reason for this is that the genome composition of anaerobic fungi like P. indianae are extremely AT-biased with a particular preference for rare and semi-rare AT-rich tRNAs in E coli, which are not explicitly predicted by standard codon adaptation indices (CAI). Native P. indianae genes with these extreme biases create drastic growth defects in E. coli (up to 69% reduction in growth), which is not seen in genes from other organisms with similar CAIs. However, codon optimization rescues growth, allowing for gene evaluation. In this manner, we demonstrate that anaerobic fungal homologs such as PI.atoB are more active than S. cerevisiae homologs in a hybrid pathway, increasing the production of mevalonate up to 2.5 g/L (more than two-fold) and reducing waste carbon to acetate by ~90% under the conditions tested. This work demonstrates the bioproduction potential of anaerobic fungal enzyme homologs and how the analysis of codon utilization enables the study of otherwise difficult to express genes that have applications in biocatalysis and natural product discovery.

60 APPLIED LIFE SCIENCES↗

gRNA-SeqRET: a universal tool for targeted and genome-scale gRNA design and sequence extraction for prokaryotes and eukaryotes

High-throughput genetic screening is frequently employed to rapidly associate gene with phenotype and establish sequence-function relationships. With the advent of CRISPR technology, and the ability to functionally interrogate previously genetically recalcitrant organisms, non-model organisms can be investigated using pooled guide RNA (gRNA) libraries and sequencing-based assays to quantitatively assess fitness of every targeted locus in parallel. To aid the construction of pooled gRNA assemblies, we have developed an in silico design workflow for gRNA selection using the gRNA Sequence Region Extraction Tool (gRNA-SeqRET). Built upon the previously developed CCTop, gRNA-SeqRET enables automated, scalable design of gRNA libraries that target user-specified regions or whole genomes of any prokaryote or eukaryote. Additionally, gRNA-SeqRET automates the bulk extraction of any regions of sequence relative to genes or other features, aiding in the design of homology arms for insertion or deletion constructs. We also assess in silico the application of a designed gRNA library to other closely related genomes and demonstrate that for very closely related organisms Average Nucleotide Identity (ANI) > 95% a large fraction of the library may be of relevance. The gRNA-SeqRET web application pipeline can be accessed at https://grna.jgi.doe.gov. The source code is comprised of freely available software tools and customized Python scripts, and is available at https://bitbucket.org/berkeleylab/grnadesigner/src/master/ under a modified BSD open-source license (https://bitbucket.org/berkeleylab/grnadesigner).

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↗

Diel Transcriptional Oscillations of a Plastid Antiporter Reflect Increased Resilience of Thalassiosira pseudonana in Elevated CO 2

Acidification of the ocean due to high atmospheric CO 2 levels may increase the resilience of diatoms causing dramatic shifts in abiotic and biotic cycles with lasting implications on marine ecosystems. Here, we report a potential bioindicator of a shift in the resilience of a coastal and centric model diatom Thalassiosira pseudonana under elevated CO 2 . Specifically, we have discovered, through EGFP-tagging, a plastid membrane localized putative Na + (K + )/H + antiporter that is significantly upregulated at >800 ppm CO 2 , with a potentially important role in maintaining pH homeostasis. Notably, transcript abundance of this antiporter gene was relatively low and constant over the diel cycle under contemporary CO 2 conditions. In future acidified oceanic conditions, dramatic oscillation with >10-fold change between nighttime (high) and daytime (low) transcript abundances of the antiporter was associated with increased resilience of T. pseudonana. By analyzing metatranscriptomic data from the Tara Oceans project, we demonstrate that phylogenetically diverse diatoms express homologs of this antiporter across the globe. We propose that the differential between night- and daytime transcript levels of the antiporter could serve as a bioindicator of a shift in the resilience of diatoms in response to high CO 2 conditions in marine environments.

59 BASIC BIOLOGICAL SCIENCES↗

Acarbose impairs gut Bacteroides growth by targeting intracellular glucosidases

ABSTRACT Acarbose is a type 2 diabetes medicine that prevents dietary starch breakdown into glucose by inhibiting host amylase and glucosidase enzymes. Numerous gut species in theBacteroidesgenus enzymatically break down starch and change in relative abundance within the gut microbiome in acarbose-treated individuals. To mechanistically explain this observation, we used two model starch-degradingBacteroides,Bacteroides ovatus(Bo), andBacteroides thetaiotaomicron(Bt). Bt growth on starch polysaccharides is severely impaired by acarbose, whereas Bo growth is much less affected by the drug. TheBacteroidesuse a starch utilization system (Sus) to grow on starch. We hypothesized that Bo and Bt Sus enzymes are differentially inhibited by acarbose. Instead, we discovered that although acarbose primarily targets the Sus periplasmic GH97 enzymes in both organisms, the drug affects starch processing at multiple other points. Acarbose competes for transport through the TonB-dependent SusC proteins and binds to the Sus transcriptional regulators. Furthermore, Bo expresses a non-Sus GH97 (BoGH97D) when grown in starch with acarbose. The Bt homolog, BtGH97H, is not expressed in the same conditions, nor can overexpression of BoGH97D complement the Bt growth inhibition in the presence of acarbose. This work informs us about unexpected complexities of Sus function and regulation inBacteroides, including variation between related species. Furthermore, this indicates that the gut microbiome may be a source of variable response to acarbose treatment for diabetes. IMPORTANCE Acarbose is a type 2 diabetes medication that works primarily by stopping starch breakdown into glucose in the small intestine. This is accomplished by the inhibition of host enzymes, leading to better blood sugar control via reduced ability to derive glucose from dietary starches. The drug and undigested starch travel to the large intestine where acarbose interferes with the ability of some bacteria to grow on starch. However, little is known about how gut bacteria interact with acarbose, including microbes that can use starch as a carbon source. Here, we show that two gut species,Bacteroides ovatus(Bo) andBacteroides thetaiotaomicron(Bt), respond differently to acarbose: Bt growth is inhibited by acarbose, while Bo growth is less affected. We reveal a complex set of mechanisms involving differences in starch import and sensing behind the different Bo and Bt responses. This indicates the gut microbiome may be a source of variable response to acarbose treatment for diabetes via complex mechanisms in common gut microbes.

Microbiology↗

Impact of uranium oxide (UO 2 ) fuel with molybdenum (Mo) inserts on pressurized water reactor performance and safety

This work investigates nuclear reactor performance and safety characteristics of UO 2 with high thermal conductivity Mo insert structures by using multiphysics modeling techniques. Additionally, the purpose of this study is to use scoping analyses to quantify the impact of using Mo inserts from neutronic and heat transfer standpoints. Attention is given to reactor performance parameters, such as cycle length, maximum fuel temperature, temperature gradients in the fuel, and stored energy in the fuel. The finite-element code BISON and the Monte Carlo particle transport code Serpent were used to perform sensitivity analyses on the Mo insert geometry to optimize the insert design and inform larger scale modeling that required the homogenization of the UO 2 and Mo. Although BISON is often used as a fuel performance analysis tool, it is used in this context for heat transfer analysis only. Fuel performance optimization is outside the scope of the current study, but would be important for future work focused on this concept. The results showed that the insert had little impact on neutronic performance and that homogenizing the UO 2 and Mo was acceptable for reactor physics calculations. Reactivity temperature coefficients calculated using homogeneous UO 2 -Mo were shown to be relatively similar to UO 2 , but higher Mo content and 235 U enrichment can reduce the worth of soluble boron and control rods. The effect of insert geometry on heat transfer was much greater, and an approximately 15–20% difference in maximum fuel temperature was predicted between the best and worst performing heat transfer geometries. Furthermore, thermal conductivity calibration based on the finite element analysis results was performed to improve the accuracy of temperature predictions in reactor analysis models that homogenized the UO 2 -Mo fuel. Compared with UO 2 in a pressurized water reactor (PWR), the optimized UO 2 -Mo design increased the margin to fuel melt by 13–32% across the fuel cycle, but it requires the 235U enrichment to exceed 5% to match the cycle length of conventional UO 2 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Systematic characterization of gene function in the photosynthetic alga Chlamydomonas reinhardtii

Most genes in photosynthetic organisms remain functionally uncharacterized. Here, using a barcoded mutant library of the model eukaryotic alga Chlamydomonas reinhardtii, we determined the phenotypes of more than 58,000 mutants under more than 121 different environmental growth conditions and chemical treatments. A total of 59% of genes are represented by at least one mutant that showed a phenotype, providing clues to the functions of thousands of genes. Mutant phenotypic profiles place uncharacterized genes into functional pathways such as DNA repair, photosynthesis, the CO 2 -concentrating mechanism and ciliogenesis. We illustrate the value of this resource by validating phenotypes and gene functions, including three new components of an actin cytoskeleton defense pathway. The data also inform phenotype discovery in land plants; mutants in Arabidopsis thaliana genes exhibit phenotypes similar to those we observed in their Chlamydomonas homologs. We anticipate that this resource will guide the functional characterization of genes across the tree of life.

59 BASIC BIOLOGICAL SCIENCES↗

Functional role of myosin-binding protein H in thick filaments of developing vertebrate fast-twitch skeletal muscle

Myosin-binding protein H (MyBP-H) is a component of the vertebrate skeletal muscle sarcomere with sequence and domain homology to myosin-binding protein C (MyBP-C). Whereas skeletal muscle isoforms of MyBP-C (fMyBP-C, sMyBP-C) modulate muscle contractility via interactions with actin thin filaments and myosin motors within the muscle sarcomere “C-zone,” MyBP-H has no known function. This is in part due to MyBP-H having limited expression in adult fast-twitch muscle and no known involvement in muscle disease. Quantitative proteomics reported here reveal that MyBP-H is highly expressed in prenatal rat fast-twitch muscles and larval zebrafish, suggesting a conserved role in muscle development and prompting studies to define its function. We take advantage of the genetic control of the zebrafish model and a combination of structural, functional, and biophysical techniques to interrogate the role of MyBP-H. Transgenic, FLAG-tagged MyBP-H or fMyBP-C both localize to the C-zones in larval myofibers, whereas genetic depletion of endogenous MyBP-H or fMyBP-C leads to increased accumulation of the other, suggesting competition for C-zone binding sites. Does MyBP-H modulate contractility in the C-zone? Globular domains critical to MyBP-C’s modulatory functions are absent from MyBP-H, suggesting that MyBP-H may be functionally silent. However, our results suggest an active role. In vitro motility experiments indicate MyBP-H shares MyBP-C’s capacity as a molecular “brake.” These results provide new insights and raise questions about the role of the C-zone during muscle development.

59 BASIC BIOLOGICAL SCIENCES↗

Co-locality to co-functionality: Eukaryotic gene neighborhoods as a resource for function

Diverging from the classic paradigm of random gene order in eukaryotes, gene proximity can be leveraged to systematically identify functionally related gene neighborhoods in eukaryotes, utilizing techniques pioneered in bacteria. Current methods of identifying gene neighborhoods typically rely on sequence similarity to characterized gene products. However, this approach is not robust for non model organisms like algae, which are evolutionarily distant from well-characterized model organisms. Here, we utilize a comparative genomic approach to identify evolutionarily conserved Proximal Orthologous Gene (POG) pairs conserved across at least two taxonomic classes of green algae. A total of 317 gene neighborhoods were identified. In some cases, gene proximity appears to have been conserved since before the streptophyte-chlorophyte split, 1,000 million years ago. Using functional inferences derived from reconstructed evolutionary relationships, we identified several novel functional clusters. A putative mycosporine-like amino acid (MAA), "sunscreen", neighborhood contains genes similar to either vertebrate or cyanobacterial pathways, suggesting a novel mosaic biosynthetic pathway in green algae. One of two putative arsenic-detoxification neighborhoods includes an organoarsenical transporter (ArsJ), a glyceraldehyde 3-phosphate dehydrogenase-like gene, homologs of which are involved in arsenic detoxification in bacteria, and a novel algal-specific phosphoglycerate kinase-like gene (PGK). Mutants of the ArsJ-like transporter and PGK-like genes in Chlamydomonas reinhardtii were found to be sensitive to arsenate, providing experimental support for the role of these identified neighbors in resistance to arsenate. Potential evolutionary origins of neighborhoods are discussed, and updated annotations for formerly poorly annotated genes are presented, highlighting the potential of this strategy for functional annotation.

59 BASIC BIOLOGICAL SCIENCES↗

Discovery, characterization, and metabolic engineering of Rieske non-heme iron monooxygenases for guaiacol O-demethylation

Aryl-O-demethylation is a common rate-limiting step in the catabolism of lignin-related compounds, including guaiacol. In this work, we used randomly barcoded transposon insertion sequencing (RB-TnSeq) in the bacterium Novosphingobium aromaticivorans to identify a Rieske-type guaiacol O-demethylase, GdmA. Similarity searches identified GdmA homologs in other bacteria, along with candidate reductase partners, denoted GdmB. GdmAB combinations were biochemically characterized for activity with several lignin-related substrates. Structural and sequence comparisons of vanillate- and guaiacol-specific O-demethylase active sites revealed conserved hallmarks of substrate specificity. GdmAB combinations were also evaluated in Pseudomonas putida KT2440, which does not natively utilize guaiacol. GdmAB from Cupriavidus necator N-1 demonstrated the highest rate of guaiacol turnover in vitro and in engineered P. putida strains and notably higher catalytic efficiency than a cytochrome P450 system (GcoAB) and the vanillate Rieske-type O-demethylase from P. putida (VanAB). The GdmAB O-demethylases described here expand the suite of options for microbial conversion of a model lignin-derived substrate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An atypical BRCT–BRCT interaction with the XRCC1 scaffold protein compacts human DNA Ligase IIIα within a flexible DNA repair complex

The XRCC1–DNA ligase IIIα complex (XL) is critical for DNA single-strand break repair, a key target for PARP inhibitors in cancer cells deficient in homologous recombination. Here, we combined biophysical approaches to gain insights into the shape and conformational flexibility of the XL as well as XRCC1 and DNA ligase IIIα (LigIIIα) alone. Structurally-guided mutational analyses based on the crystal structure of the human BRCT–BRCT heterodimer identified the network of salt bridges that together with the N-terminal extension of the XRCC1 C-terminal BRCT domain constitute the XL molecular interface. Coupling size exclusion chromatography with small angle X-ray scattering and multiangle light scattering (SEC-SAXS–MALS), we determined that the XL is more compact than either XRCC1 or LigIIIα, both of which form transient homodimers and are highly disordered. The reduced disorder and flexibility allowed us to build models of XL particles visualized by negative stain electron microscopy that predict close spatial organization between the LigIIIα catalytic core and both BRCT domains of XRCC1. Together our results identify an atypical BRCT–BRCT interaction as the stable nucleating core of the XL that links the flexible nick sensing and catalytic domains of LigIIIα to other protein partners of the flexible XRCC1 scaffold.

59 BASIC BIOLOGICAL SCIENCES↗

Distributed Many-to-Many Protein Sequence Alignment using Sparse Matrices

Identifying similar protein sequences is a core step in many computational biology pipelines such as detection of homologous protein sequences, generation of similarity protein graphs for downstream analysis, functional annotation, and gene location. Performance and scalability of protein similarity search have proven to be a bottleneck in many bioinformatics pipelines due to increase in cheap and abundant sequencing data. This work presents a new distributed-memory software PASTIS. PASTIS relies on sparse matrix computations for efficient identification of possibly similar proteins. We use distributed sparse matrices for scalability and show that the sparse matrix infrastructure is a great fit for protein similarity search when coupled with a fully-distributed dictionary of sequences that allow remote sequence requests to be fulfilled. Our algorithm incorporates the unique bias in amino acid sequence substitution in search without altering basic sparse matrix model, and in turn, achieves ideal scaling up to millions of protein sequences.

97 MATHEMATICS AND COMPUTING↗