Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sequence comparison”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Shaping Nanobodies and Intrabodies against Proteoforms

Proteoforms expand genomic diversity and direct developmental processes. While high-resolution mass spectrometry has accelerated characterization of proteoforms, molecular techniques working to bind, and disrupt the function of specific proteoforms have lagged behind. In this study, we worked to develop intrabodies capable of binding specific proteoforms. Herein we employed a synthetic camelid nanobody library expressed in yeast to identify nanobody binders of different SARS-CoV-2 receptor binding domain (RBD) proteoforms. Importantly, employment of the positive and negative selection mechanism inherent to the synthetic system allowed for amplification of nanobody-expressing yeast that bind to the original (Wuhan strain RBD) but not the E484K (Beta variant) mutation. Nanobodies raised against specific RBD proteoforms were validated by yeast-2-hybrid analysis and sequence comparisons. These results provide a clear framework for development of proteoform-specific nanobodies and intrabodies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Structural basis for expanded substrate specificities of human long chain acyl-CoA dehydrogenase and related acyl-CoA dehydrogenases

Crystal structures of human long-chain acyl-CoA dehydrogenase (LCAD) and the catalytically inactive Glu291Gln mutant, have been determined. These structures suggest that LCAD harbors functions beyond its historically defined role in mitochondrial β-oxidation of long and medium-chain fatty acids. LCAD is a homotetramer containing one FAD per 43 kDa subunit with Glu291 as the catalytic base. The substrate binding cavity of LCAD reveals key differences which makes it specific for longer and branched chain substrates. The presence of Pro132 near the start of the E helix leads to helix unwinding that, together with adjacent smaller residues, permits binding of bulky substrates such as 3α, 7α, l2α-trihydroxy-5β-cholestan-26-oyl-CoA. This structural element is also utilized by ACAD11, a eucaryotic ACAD of unknown function, as well as bacterial ACADs known to metabolize sterol substrates. Sequence comparison suggests that ACAD10, another ACAD of unknown function, may also share this substrate specificity. These results suggest that LCAD, ACAD10, ACAD11 constitute a distinct class of eucaryotic acyl CoA dehydrogenases.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Transmission of Multiple HIV-1 Subtype C Transmitted/founder Viruses into the Same Recipients Was not Determined by Modest Phenotypic Differences

A severe bottleneck exists during HIV-1 mucosal transmission. However, viral properties that determine HIV-1 transmissibility are not fully elucidated. We identified multiple transmitted/founder (T/F) viruses in six HIV-1-infected subjects by analyzing whole genome sequences. Comparison of biological phenotypes of different T/F viruses from the same individual allowed us to more precisely identify critical determinants for viral transmissibility since they were transmitted under similar conditions. All T/F viruses used coreceptor CCR5, while no T/F viruses used CXCR4 or GPR15. However, the efficiency for different T/F viruses from the same individual to use CCR5 was significantly variable, and the differences were even more significant for usage of coreceptors FPRL1, CCR3 and APJ. Resistance to IFN-α was also different between T/F viruses in 2 of 3 individuals. The relative fitness between T/F viruses from the same subject was highly variable (2–6%). Importantly, the levels of coreceptor usage efficiency, resistance to IFN-α and viral fitness were not associated with proportions of T/F viruses in each individual during acute infection. Our results show that the modest but significant differences in coreceptor usage efficiency, IFN-α sensitivity and viral fitness each alone may not play a critical role in HIV-1 transmission.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Unusual respiratory capacity and nitrogen metabolism in a Parcubacterium (OD1) of the Candidate Phyla Radiation

The Candidate Phyla Radiation (CPR) is a large group of bacteria, the scale of which approaches that of all other bacteria. CPR organisms are inferred to depend on other community members for many basic cellular building blocks and all appear to be obligate anaerobes. To date, there has been no evidence for any significant respiratory capacity in an organism from this radiation. Here we report a curated draft genome for Candidatus Parcunitrobacter nitroensis' a member of the Parcubacteria (OD1) superphylum of the CPR. The genome encodes versatile energy pathways, including fermentative and respiratory capacities, nitrogen and fatty acid metabolism, as well as the first complete electron transport chain described for a member of the CPR. The sequences of all of these enzymes are highly divergent from sequences found in other organisms, suggesting that these capacities were not recently acquired from non-CPR organisms. Although the wide respiration-based repertoire points to a different lifestyle compared to other CPR bacteria, we predict similar obligate dependence on other organisms or the microbial community. The results substantially expand the known metabolic potential of CPR bacteria, although sequence comparisons indicate that these capacities are very rare in members of this radiation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

A plasmid borne, functionally novel glycoside hydrolase family 30 subfamily 8 endoxylanase from solventogenic Clostridium

Glycoside hydrolase family 30 subfamily 8 (GH30-8) β-1,4-endoxylanases are known for their appendage-dependent function requiring recognition of an α-1,2-linked glucuronic acid (GlcA) common to glucuronoxylans for hydrolysis. Structural studies have indicated that the GlcA moiety of glucuronoxylans is coordinated through six hydrogen bonds and a salt bridge. These GlcA-dependent endoxylanases do not have significant activity on xylans that do not bear GlcA substitutions such as unsubstituted linear xylooligosaccharides or cereal bran arabinoxylans. In the present study, we present the structural and biochemical characteristics of xylanase 30A from Clostridium acetobutylicum (CaXyn30A) which was originally selected for study due to predicted structural differences within the GlcA coordination loops. Amino acid sequence comparisons indicated that this Gram-positive-derived GH30-8 more closely resembles Gram-negative derived forms of these endoxylanases: a hypothesis borne out in the developed crystallographic structure model of the CaXyn30A catalytic domain (CaXyn30A-CD). CaXyn30A-CD hydrolyzes xylans to linear and substituted oligoxylosides showing the greatest rate with the highly arabinofuranose (Araf)-substituted cereal arabinoxylans. CaXyn30A-CD hydrolyzes xylooligosaccharides larger than xylotriose and shows an increased relative rate of hydrolysis for xylooligosaccharides containing α-1,2-linked arabinofuranose substitutions. Biochemical analysis confirms that CaXyn30A benefits from five xylose-binding subsites which extend from the -3 subsite to the +2 subsite of the binding cleft. These studies indicate that CaXyn30A is a GlcA-independent endoxylanase that may have evolved for the preferential recognition of α-1,2-Araf substitutions on xylan chains.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Functional analysis of (4 S )-limonene synthase mutants reveals determinants of catalytic outcome in a model monoterpene synthase

We used crystal structural data for (4S)-limonene synthase [(4S)-LS] of spearmint (Mentha spicata L.) to infer which amino acid residues are in close proximity to the substrate and carbocation intermediates of the enzymatic reaction. Alanine-scanning mutagenesis of 48 amino acids combined with enzyme fidelity analysis [percentage of (-)-limonene produced] indicated which residues are most likely to constitute the active site. Furthermore, the mutation of residues W324 and H579 caused a significant drop in enzyme activity and formation of products (myrcene, linalool, and terpineol) characteristic of a premature termination of the reaction. A double mutant (W324A/H579A) had no detectable enzyme activity, indicating that either substrate binding or the terminating reaction was impaired. Exchanges to other aromatic residues (W324H, W324F, W324Y, H579F, H579Y, and H579W) resulted in enzyme catalysts with significantly reduced activity. Sequence comparisons across the angiosperm lineage provided evidence that W324 is a conserved residue, whereas the position equivalent to H579 is occupied by aromatic residues (H, F, or Y). Our results are consistent with a critical role of W324 and H579 in the stabilization of carbocation intermediates. Finally, the potential of these residues to serve as the catalytic base facilitating the terminal deprotonation reaction is discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2012

Structural and functional evolution of isopropylmalate dehydrogenases in the leucine and glucosinolate pathways of Arabidopsis thaliana

The methionine chain-elongation pathway is required for aliphatic glucosinolate biosynthesis in plants and evolved from leucine biosynthesis. In Arabidopsis thaliana, three 3-isopropylmalate dehydrogenases (AtIPMDHs) play key roles in methionine chain-elongation for the synthesis of aliphatic glucosinolates (e.g. AtIPMDH1) and leucine (e.g. AtIPMDH2 and AtIPMDH3). Here we elucidate the molecular basis underlying the metabolic specialization of these enzymes. The 2.25 {angstrom} resolution crystal structure of AtIPMDH2 was solved to provide the first detailed molecular architecture of a plant IPMDH. Modeling of 3-isopropylmalate binding in the AtIPMDH2 active site and sequence comparisons of prokaryotic and eukaryotic IPMDH suggest that substitution of one active site residue may lead to altered substrate specificity and metabolic function. Site-directed mutagenesis of Phe-137 to a leucine in AtIPMDH1 (AtIPMDH1-F137L) reduced activity toward 3-(2'-methylthio)ethylmalate by 200-fold, but enhanced catalytic efficiency with 3-isopropylmalate to levels observed with AtIPMDH2 and AtIPMDH3. Conversely, the AtIPMDH2-L134F and AtIPMDH3-L133F mutants enhanced catalytic efficiency with 3-(2'-methylthio)ethylmalate {approx}100-fold and reduced activity for 3-isopropylmalate. Furthermore, the altered in vivo glucosinolate profile of an Arabidopsis ipmdh1 T-DNA knock-out mutant could be restored to wild-type levels by constructs expressing AtIPMDH1, AtIPMDH2-L134F, or AtIPMDH3-L133F, but not by AtIPMDH1-F137L. These results indicate that a single amino acid substitution results in functional divergence of IPMDH in planta to affect substrate specificity and contributes to the evolution of specialized glucosinolate biosynthesis from the ancestral leucine pathway.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2013

The Catalytic Site Atlas 2.0: cataloging catalytic sites and residues identified in enzymes

Understanding which are the catalytic residues in an enzyme and what function they perform is crucial to many biology studies, particularly those leading to new therapeutics and enzyme design. The original version of the Catalytic Site Atlas (CSA) (http://www.ebi.ac.uk/thornton-srv/databases/CSA) published in 2004, which catalogs the residues involved in enzyme catalysis in experimentally determined protein structures, had only 177 curated entries and employed a simplistic approach to expanding these annotations to homologous enzyme structures. Here we present a new version of the CSA (CSA 2.0), which greatly expands the number of both curated (968) and automatically annotated catalytic sites in enzyme structures, utilizing a new method for annotation transfer. The curated entries are used, along with the variation in residue type from the sequence comparison, to generate 3D templates of the catalytic sites, which in turn can be used to find catalytic sites in new structures. To ease the transfer of CSA annotations to other resources a new ontology has been developed: the Enzyme Mechanism Ontology, which has permitted the transfer of annotations to Mechanism, Annotation and Classification in Enzymes (MACiE) and UniProt Knowledge Base (UniProtKB) resources. The CSA database schema has been re-designed and both the CSA data and search capabilities are presented in a new modern web interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

A temperature-sensitive FERONIA mutant allele that alters root hair growth

Abstract In plants, root hairs undergo a highly polarized form of cell expansion called tip-growth, in which cell wall deposition is restricted to the root hair apex. In order to identify essential cellular components that might have been missed in earlier genetic screens, we identified conditional temperature-sensitive (ts) root hair mutants by ethyl methanesulfonate mutagenesis in Arabidopsis thaliana. Here, we describe one of these mutants, feronia-temperature sensitive (fer-ts). Mutant fer-ts seedlings were unaffected at normal temperatures (20°C), but failed to form root hairs at elevated temperatures (30°C). Map based-cloning and whole-genome sequencing revealed that fer-ts resulted from a G41S substitution in the extracellular domain of FERONIA (FER). A functional fluorescent fusion of FER containing the fer-ts mutation localized to plasma membranes, but was subject to enhanced protein turnover at elevated temperatures. While tip-growth was rapidly inhibited by addition of rapid alkalinization factor 1 (RALF1) peptides in both wild-type and fer-ts mutants at normal temperatures, root elongation of fer-ts seedlings was resistant to added RALF1 peptide at elevated temperatures. Additionally, at elevated temperatures fer-ts seedlings displayed altered reactive oxygen species (ROS) accumulation upon auxin treatment and phenocopied constitutive fer mutant responses to a variety of plant hormone treatments. Molecular modeling and sequence comparison with other Catharanthus roseus receptor-like kinase 1L (CrRLK1L) receptor family members revealed that the mutated glycine in fer-ts is highly conserved, but is not located within the recently characterized RALF23 and LORELI-LIKE-GLYCOPROTEIN 2 binding domains, perhaps suggesting that fer-ts phenotypes may not be directly due to loss of binding to RALF1 peptides.

Plant Sciences↗
DOE OSTI2021

High-Performance Deep Learning Toolbox for Genome-Scale Prediction of Protein Structure and Function

Computational biology is one of many scientific disciplines ripe for innovation and acceleration with the advent of high-performance computing (HPC). In recent years, the field of machine learning has also seen significant benefits from adopting HPC practices. In this work, we present a novel HPC pipeline that incorporates various machine-learning approaches for structure-based functional annotation of proteins on the scale of whole genomes. Our pipeline makes extensive use of deep learning and provides computational insights into best practices for training advanced deep-learning models for high-throughput data such as proteomics data. We showcase methodologies our pipeline currently supports and detail future tasks for our pipeline to envelop, including large-scale sequence comparison using SAdLSA and prediction of protein tertiary structures using AlphaFold2.

Gao, Mu↗
DOE OSTI2025

Amino acid substrate specificities and tissue expression profiles of the nine CYP79A encoding genes in Sorghum bicolor

Cytochrome P450s of the CYP79 family catalyze two N-hydroxylation reactions, converting a selected number of amino acids into the corresponding oximes. The sorghum genome (Sorghum bicolor) harbours nine CYP79A encoding genes, and here sequence comparisons of the CYP79As along with their substrate recognition sites (SRSs) are provided. The substrate specificity of previously uncharacterized CYP79As was investigated by transient expression in Nicotiana benthamiana and subsequent transformation of the oximes formed into the corresponding stable oxime glucosides catalyzed by endogenous UDPG-glucosyltransferases (UGTs). CYP79A61 uses phenylalanine as a substrate, whereas CYP79A91, CYP79A93, and CYP79A95 use valine and isoleucine as substrates, with CYP79A93 showing the ability also to use phenylalanine. CYP79A94 uses isoleucine as a substrate. Analysis of 249 sorghum transcriptomes from two different sorghum cultivars showed the expression levels and tissue-specific expression of the CYP79As. CYP79A1 is the committed gene in dhurrin formation and was the highest expressed gene in most tissues/organs. CYP79A61 was primarily expressed in fully developed leaf blades and leaf sheaths. CYP79A91 and CYP79A92 were expressed mainly in roots >200 cm below ground, while CYP79A93 and CYP79A94 were most highly expressed in the leaf collar and leaf sheath, respectively. Here, the possible signalling effects of the oximes and their metabolites produced in different sorghum tissues are discussed.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Disruption of the tagF Orthologue in the epa Locus Variable Region of Enterococcus faecalis Causes Cell Surface Changes and Suppresses an eep -Dependent Lysozyme Resistance Phenotype

The disease-producing capacity of the opportunistic pathogen Enterococcus faecalis is enhanced by the ability of the bacterium to evade killing by antimicrobial agents. Survival of E. faecalis in the presence of the human antimicrobial enzyme lysozyme is mediated in part by the site 2 metalloprotease Eep; however, a complete model of enterococcal lysozyme resistance has not been elucidated. To better understand the molecular basis for lysozyme resistance in E. faecalis, we analyzed Δeep suppressor mutants that acquire resistance to lysozyme through mutation of the gene OG1RF_11713, a predicted teichoic acid biosynthesis-encoding gene located within the variable region of the enterococcal polysaccharide antigen (epa) locus. Sequence comparisons revealed that OG1RF_11713 is most similar to the cytidine-5'-diphosphate (CDP)-glycerol:poly-(glycerolphosphate)glycerophosphotransferase TagF from Staphylococcus epidermidis. Inactivation of OG1RF_11713 in both the wild-type and Δeep genetic backgrounds was sufficient to increase the resistance of E. faecalis OG1RF to lysozyme. Minimal amounts of N-acetylgalactosamine were detectable in cell wall carbohydrate extracts of OG1RF_11713 deletion mutants, and this was associated with a reduction in negative cell surface charge. Targeted disruption of OG1RF_11713 was also associated with increased susceptibility to the antibiotic polymyxin B and membrane-targeting detergents and decreased susceptibility to the lantibiotic nisin. Furthermore, this work implicates OG1RF_11713 as a major determinant of cell envelope integrity and provides further validation that lysozyme resistance is intrinsically linked to the modification of enterococcal cell wall polysaccharides.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Lateral Gene Transfer in a Heavy Metal-Contaminated-Groundwater Microbial Community

ABSTRACT Unraveling the drivers controlling the response and adaptation of biological communities to environmental change, especially anthropogenic activities, is a central but poorly understood issue in ecology and evolution. Comparative genomics studies suggest that lateral gene transfer (LGT) is a major force driving microbial genome evolution, but its role in the evolution of microbial communities remains elusive. To delineate the importance of LGT in mediating the response of a groundwater microbial community to heavy metal contamination, representative Rhodanobacter reference genomes were sequenced and compared to shotgun metagenome sequences. 16S rRNA gene-based amplicon sequence analysis indicated that Rhodanobacter populations were highly abundant in contaminated wells with low pHs and high levels of nitrate and heavy metals but remained rare in the uncontaminated wells. Sequence comparisons revealed that multiple geochemically important genes, including genes encoding Fe 2+ /Pb 2+ permeases, most denitrification enzymes, and cytochrome c 553 , were native to Rhodanobacter and not subjected to LGT. In contrast, the Rhodanobacter pangenome contained a recombinational hot spot in which numerous metal resistance genes were subjected to LGT and/or duplication. In particular, Co 2+ /Zn 2+ /Cd 2+ efflux and mercuric resistance operon genes appeared to be highly mobile within Rhodanobacter populations. Evidence of multiple duplications of a mercuric resistance operon common to most Rhodanobacter strains was also observed. Collectively, our analyses indicated the importance of LGT during the evolution of groundwater microbial communities in response to heavy metal contamination, and a conceptual model was developed to display such adaptive evolutionary processes for explaining the extreme dominance of Rhodanobacter populations in the contaminated groundwater microbiome. IMPORTANCE Lateral gene transfer (LGT), along with positive selection and gene duplication, are the three main mechanisms that drive adaptive evolution of microbial genomes and communities, but their relative importance is unclear. Some recent studies suggested that LGT is a major adaptive mechanism for microbial populations in response to changing environments, and hence, it could also be critical in shaping microbial community structure. However, direct evidence of LGT and its rates in extant natural microbial communities in response to changing environments is still lacking. Our results presented in this study provide explicit evidence that LGT played a crucial role in driving the evolution of a groundwater microbial community in response to extreme heavy metal contamination. It appears that acquisition of genes critical for survival, growth, and reproduction via LGT is the most rapid and effective way to enable microorganisms and associated microbial communities to quickly adapt to abrupt harsh environmental stresses.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Classification of Complete Proteomes of Different Organisms and Protein Sets Based on Their Protein Distributions in Terms of Some Key Attributes of Proteins

The existence of complete genome sequences makes it important to develop different approaches for classification of large-scale data sets and to make extraction of biological insights easier. Here, we propose an approach for classification of complete proteomes/protein sets based on protein distributions on some basic attributes. We demonstrate the usefulness of this approach by determining protein distributions in terms of two attributes: protein lengths and protein intrinsic disorder contents (ID). The protein distributions based on L and ID are surveyed for representative proteome organisms and protein sets from the three domains of life. The two-dimensional maps (designated as fingerprints here) from the protein distribution densities in the LD space defined by ln( L ) and ID are then constructed. The fingerprints for different organisms and protein sets are found to be distinct with each other, and they can therefore be used for comparative studies. As a test case, phylogenetic trees have been constructed based on the protein distribution densities in the fingerprints of proteomes of organisms without performing any protein sequence comparison and alignments. The phylogenetic trees generated are biologically meaningful, demonstrating that the protein distributions in the LD space may serve as unique phylogenetic signals of the organisms at the proteome level.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Structural insights into the mechanism defining substrate affinity in Arabidopsis thaliana dUTPase: the role of tryptophan 93 in ligand orientation

Deoxyuridine triphosphate nucleotidohydrolase (dUTPase) hydrolyzes dUTP to dUMP and pyrophosphate to maintain the cellular thymine-uracil ratio. dUTPase is also a target for cancer chemotherapy. However, the mechanism defining its substrate affinity remains unclear. Sequence comparisons of various dUTPases revealed that Arabidopsis thaliana dUTPase has a unique tryptophan at position 93, which potentially contributes to its degree of substrate affinity. To better understand the roles of tryptophan 93, A. thaliana dUTPase was studied. Enzyme assays showed that A. thaliana dUTPase belongs to a high-affinity group of isozymes, which also includes the enzymes from Escherichia coli and Mycobacterium tuberculosis. Enzymes from Homo sapiens and Saccharomyces cerevisiae are grouped as low-affinity dUTPases. The structure of the homo-trimeric A. thaliana dUTPase showed three active sites, each with a different set of ligand interactions between the amino acids and water molecules. On an α-helix, tryptophan 93 appears to keep serine 89 in place via a water molecule and to specifically direct the ligand. Upon being oriented in the active site, the C-terminal residues close the active site to promote the reaction. In the high-affinity group, the prefixed direction of the serine residues was oriented by a positively charged residue located four amino acids away, while low-affinity enzymes possess small hydrophobic residues at the corresponding sites.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Mouse IgG2c Fc loop residues promote greater receptor-binding affinity than mouse IgG2b or human IgG1

The structures of non-human antibodies are largely unstudied despite the potential for the identification of alternative structural motifs and physical properties that will benefit a basic understanding of protein and immune system evolution as well as highlight unexplored motifs to improve therapeutic monoclonal antibody. Here we probe the structure and receptor-binding properties of the mouse IgG2c crystallizable fragment (Fc) to compare to mouse IgG2b and human IgG1 Fcs. Models of mIgG2c Fc determined by x-ray crystallography with a complex-type biantennary (to 2.05 Å) or a truncated (1)GlcNAc asparagine-linked (N)-glycan attached (to 2.04 Å) show differences in key regions related to mouse Fc γ receptor IV (mFcγRIV) binding. Mouse IgG2c forms different non-bonded interactions between the BC, DE and FG loops than the highly-conserved mIgG2b and binds to FcγRIV with 4.7-fold greater affinity in the complex-type glycoform. Secondary structural elements surrounding the Asn297 site of glycosylation form longer beta strands in the truncated mIgG2c Fc glycoform when compared to mIgG2c with the complex-type N-glycan. Solution NMR spectroscopy of the N-linked (1)GlcNAc residues show differences between mIgG2b, 2c and hIgG1 Fc that correlate to receptor binding affinity. Mutations targeting differences in mIgG2 DE and FG loops decreased affinity of mIgG2c for FcγRIV and increased affinity of mIgG2b. Changes in NMR spectra of the mutant Fc proteins mirrored these changes in affinity. Our studies identified structural and functional differences in highly conserved molecules that were not predicted from primary sequence comparison.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Investigation of protein quaternary structure via stoichiometry and symmetry information

The Protein Data Bank (PDB) is the single worldwide archive of experimentally-determined three-dimensional (3D) structures of proteins and nucleic acids. As of January 2017, the PDB housed more than 125,000 structures and was growing by more than 11,000 structures annually. Since the 3D structure of a protein is vital to understand the mechanisms of biological processes, diseases, and drug design, correct oligomeric assembly information is of critical importance. Unfortunately, the biologically relevant oligomeric form of a 3D structure is not directly obtainable by X-ray crystallography, whilst in solution methods (NMR or single particle EM) it is known from the experiment. Instead, this information may be provided by the PDB Depositor as metadata coming from additional experiments, be inferred by sequence-sequence comparisons with similar proteins of known oligomeric state, or predicted using software, such as PISA (Proteins, Interfaces, Structures and Assemblies) or EPPIC (Evolutionary Protein Protein Interface Classifier). Despite significant efforts by professional PDB Biocurators during data deposition, there remain a number of structures in the archive with incorrect quaternary structure descriptions (or annotations). Further investigation is, therefore, needed to evaluate the correctness of quaternary structure annotations. In this study, we aim to identify the most probable oligomeric states for proteins represented in the PDB. Our approach evaluated the performance of four independent prediction methods, including text mining of primary publications, inference from homologous protein structures, and two computational methods (PISA and EPPIC). Aggregating predictions to give consensus results outperformed all four of the independent prediction methods, yielding 83% correct, 9% wrong, and 8% inconclusive predictions, when tested with a well-curated benchmark dataset. We have developed a freely-available web-based tool to make this approach accessible to researchers and PDB Biocurators.

59 BASIC BIOLOGICAL SCIENCES↗