Enzyme-assisted high throughput sequencing of an expanded genetic alphabet at single base resolution.

With just four building blocks, low sequence information density, few functional groups, poor control over folding, and difficulties in forming compact folds, natural DNA and RNA have been disappointing platforms from which to evolve receptors, ligands, and catalysts. Accordingly, synthetic biology has created "artificially expanded genetic information systems" (AEGIS) to add nucleotides, functionality, and information density. With the expected improvements seen in AegisBodies and AegisZymes, the task for synthetic biologists shifts to developing for expanded DNA the same analytical tools available to natural DNA. Here we report one of these, an enzyme-assisted sequencing of expanded genetic alphabet (ESEGA) method to sequence six-letter AEGIS DNA. We show how ESEGA analyses this DNA at single base resolution, and applies it to optimized conditions for six-nucleotide PCR, assessing the fidelity of various DNA polymerases, and extending this to AEGIS components with functional groups. This supports the renewed exploitation of expanded DNA alphabets in biotechnology.

Prediction and verification of glycosyltransferase activity by bioinformatics analysis and protein engineering.

A significant number of proteins are annotated as functionally uncharacterized proteins. Within this protocol, we describe how to use protein family multiple sequence alignments and structural bioinformatics resources to design loss-of-function mutations of previously uncharacterized proteins within the glycosyltransferase family. We detail approaches to determine target protein active sites using three-dimensional modeling. We generate active site mutants and quantify any changes in enzymatic function by a glycosyltransferase assay. With modifications, this protocol could be applied to other metal-dependent enzymes. For complete details on the use and execution of this protocol, please refer to Ilina et al. (2022).1.

Enzyme-Assisted High Throughput Sequencing of an Expanded Genetic Alphabet at Single Base Resolution.

Many efforts have sought to apply laboratory in vitro evolution (LIVE) to natural nucleic acid (NA) scaffolds to directly evolve functional molecules. However, synthetic biology can move beyond natural NA scaffolds to create molecular systems whose libraries are far richer reservoirs of functionality than natural NAs. For example, "artificially expanded genetic information systems" (AEGIS) add up to eight nucleotides to the four found in standard NA. Even in its simplest 6-letter versions, AEGIS adds functional groups, information density, and folding motifs that natural NA libraries lack. To complete this vision, however, tools are needed to sequence molecules that are created by AEGIS LIVE. Previous sequencing approaches, including approaches from our laboratories, exhibited limited performance and lost many sequences in diverse library mixtures. Here, we present a new approach that enzymatically transforms the target AEGIS DNA. With higher transliteration efficiency and fidelity, this Enzyme-Assisted Sequencing of Expanded Genetic Alphabet (ESEGA) approach produces substantially better sequences of 6-letter (AGCTZP) DNA than previous transliteration approaches. Therefore, ESEGA facilitates precise analysis of libraries, allowing 'next-generation deep sequencing' to accurately quantify the sequences of 6-letter DNA molecules at single base resolution. We then applied ESEGA to three tasks: (a) defining optimal conditions to perform 6-nucleotide PCR (b) evaluating the fidelity of 6-nucleotide PCR with various DNA polymerases, and (c) extending that evaluation to AEGIS components functionalized with alkynyl and aromatic groups. No other approach at present has this scope, allowing this work to be the next step towards exploiting the potential of expanded DNA alphabets in biotechnology.

Enzymatic activity of glycosyltransferase GLT8D1 promotes human glioblastoma cell migration.

Glioblastoma (GBM) is the most aggressive primary brain tumor characterized by infiltrative growth of malignant glioma cells into the surrounding brain parenchyma. In this study, our analysis of GBM patient cohorts revealed a significantly higher expression of Glycosyltransferase 8 domain containing 1 (GLT8D1) compared to normal brain tissue and could be associated with impaired patient survival. Increased in vitro expression of GLT8D1 significantly enhanced migration of two different sphere-forming GBM cell lines. By in silico analysis we predicted the 3D-structure as well as the active site residues of GLT8D1. The introduction of point mutations in the predicted active site reduced its glycosyltransferase activity in vitro and consequently impaired GBM tumor cell migration. Examination of GLT8D1 interaction partners by LC-MS/MS implied proteins associated with cytoskeleton and intracellular transport as potential substrates. In conclusion, we demonstrated that the enzymatic activity of glycosyltransferase GLT8D1 promotes GBM cell migration.

The ELIXIR Human Copy Number Variations Community: building bioinformatics infrastructure for research.

Copy number variations (CNVs) are major causative contributors both in the genesis of genetic diseases and human neoplasias. While "High-Throughput" sequencing technologies are increasingly becoming the primary choice for genomic screening analysis, their ability to efficiently detect CNVs is still heterogeneous and remains to be developed. The aim of this white paper is to provide a guiding framework for the future contributions of ELIXIR's recently established human CNV Community, with implications beyond human disease diagnostics and population genomics. This white paper is the direct result of a strategy meeting that took place in September 2018 in Hinxton (UK) and involved representatives of 11 ELIXIR Nodes. The meeting led to the definition of priority objectives and tasks, to address a wide range of CNV-related challenges ranging from detection and interpretation to sharing and training. Here, we provide suggestions on how to align these tasks within the ELIXIR Platforms strategy, and on how to frame the activities of this new ELIXIR Community in the international context.

Comparative functional potency of DNA vaccines encoding Plasmodium falciparum transmission blocking target antigens Pfs48/45 and Pfs25 administered alone or in combination by in vivo electroporation in rhesus macaques.

Antibodies recognizing conformational epitopes in Pfs48/45, an antigen expressed on the surface of Plasmodium falciparum gametes and zygotes, have firmly established Pfs48/45 as a promising transmission blocking vaccine (TBV) candidate. However, it has been difficult to reproducibly express Pfs48/45 in a variety of recombinant expression systems. The goal of our studies was to evaluate functional immunogenicity of Pfs48/45 using DNA vaccine format in rhesus macaques. An additional goal was to ensure that when used in combination with another malarial antigen, specific immunity to both antigens was not compromised. For testing combination vaccines, we employed Pfs25 DNA plasmids that have previously undergone evaluations in rodents and nonhuman primates. Pfs25 is expressed on the surface of parasites after fertilization and is also a lead TBV candidate. DNA plasmids based on codon-optimized sequences of Pfs48/45 and Pfs25 were administered by in vivo electroporation, followed by a final recombinant protein boost. Our studies demonstrate that Pfs48/45 encoded by DNA plasmids is capable of inducing potent transmission blocking antibody responses, and such transmission blocking immune potency of Pfs48/45 was not compromised when tested in combination with Pfs25, These findings provide the evidence in favor of further studies on Pfs48/45 and Pfs25, either alone or in combination with other known malaria vaccine candidates for developing effective vaccines capable of interrupting malaria transmission.

Immunogenicity and malaria transmission reducing potency of Pfs48/45 and Pfs25 encoded by DNA vaccines administered by intramuscular electroporation.

Pfs48/45 and Pfs25 are leading candidates for the development of Plasmodium falciparum transmission blocking vaccines (TBV). Expression of Pfs48/45 in the erythrocytic sexual stages and presentation to the immune system during infection in the human host also makes it ideal for natural boosting. However, it has been challenging to produce a fully folded, functionally active Pfs48/45, using various protein expression platforms. In this study, we demonstrate that full-length Pfs48/45 encoded by DNA plasmids is able to induce significant transmission reducing immune responses. DNA plasmids encoding Pfs48/45 based on native (WT), codon optimized (SYN), or codon optimized and mutated (MUT1 and MUT2), to prevent any asparagine (N)-linked glycosylation were compared with or without intramuscular electroporation (EP). EP significantly enhanced antibody titers and transmission blocking activity elicited by immunization with SYN Pfs48/45 DNA vaccine. Mosquito membrane feeding assays also revealed improved functional immunogenicity of SYN Pfs48/45 (N-glycosylation sites intact) as compared to MUT1 or MUT2 Pfs48/45 DNA plasmids (all N-glycosylation sites mutated). Boosting with recombinant Pfs48/45 protein after immunization with each of the different DNA vaccines resulted in significant boosting of antibody response and improved transmission reducing capabilities of all four DNA vaccines. Finally, immunization with a combination of DNA plasmids (SYN Pfs48/45 and SYN Pfs25) also provides support for the possibility of combining antigens targeting different life cycle stages in the parasite during transmission through mosquitoes.

Polymerase Interactions with Wobble Mismatches in Synthetic Genetic Systems and Their Evolutionary Implications.

In addition to completing the Watson-Crick nucleobase matching "concept" (big pairs with small, hydrogen bond donors pair with hydrogen bond acceptors), artificially expanded genetic information systems (AEGIS) also challenge DNA polymerases with a complete set of mismatches, including wobble mismatches. Here, we explore wobble mismatches with AEGIS with DNA polymerase 1 from Escherichia coli. Remarkably, we find that the polymerase tolerates an AEGIS:standard wobble that has the same geometry as the G:T wobble that polymerases have evolved to exclude but excludes a wobble geometry that polymerases have never encountered in natural history. These results suggest certain limits to "structural analogy" and "evolutionary guidance" as tools to help synthetic biologists expand DNA alphabets.

Evolution of functional six-nucleotide DNA.

Axiomatically, the density of information stored in DNA, with just four nucleotides (GACT), is higher than in a binary code, but less than it might be if synthetic biologists succeed in adding independently replicating nucleotides to genetic systems. Such addition could also add functional groups not found in natural DNA, but useful for molecular performance. Here, we consider two new nucleotides (Z and P, 6-amino-5-nitro-3-(1'-ß-D-2'-deoxyribo-furanosyl)-2(1H)-pyridone and 2-amino-8-(1'-ß-D-2'-deoxyribofuranosyl)-imidazo[1,2-a]-1,3,5-triazin-4(8H)-one). These are designed to pair via complete Watson-Crick geometry. These were added to a library of oligonucleotides used in a laboratory in vitro evolution (LIVE) experiment; the GACTZP library was challenged to deliver molecules that bind selectively to liver cancer cells, but not to untransformed liver cells. Unlike in classical in vitro selection, low levels of mutation allow this system to evolve to create binding molecules not necessarily present in the original library. Over a dozen binding species were recovered. The best had Z and/or P in their sequences. Several had multiple, nearby, and adjacent Zs and Ps. Only the weaker binders contained no Z or P at all. This suggests that this system explored much of the sequence space available to this genetic system and that GACTZP libraries are richer reservoirs of functionality than standard libraries.

Three-dimensional topology of the SMC2/SMC4 subcomplex from chicken condensin I revealed by cross-linking and molecular modelling.

SMC proteins are essential components of three protein complexes that are important for chromosome structure and function. The cohesin complex holds replicated sister chromatids together, whereas the condensin complex has an essential role in mitotic chromosome architecture. Both are involved in interphase genome organization. SMC-containing complexes are large (more than 650 kDa for condensin) and contain long anti-parallel coiled-coils. They are thus difficult subjects for conventional crystallographic and electron cryomicroscopic studies. Here, we have used amino acid-selective cross-linking and mass spectrometry combined with structure prediction to develop a full-length molecular draft three-dimensional structure of the SMC2/SMC4 dimeric backbone of chicken condensin. We assembled homology-based molecular models of the globular heads and hinges with the lengthy coiled-coils modelled in fragments, using numerous high-confidence cross-links and accounting for potential irregularities. Our experiments reveal that isolated condensin complexes can exist with their coiled-coil segments closely apposed to one another along their lengths and define the relative spatial alignment of the two anti-parallel coils. The centres of the coiled-coils can also approach one another closely in situ in mitotic chromosomes. In addition to revealing structural information, our cross-linking data suggest that both H2A and H4 may have roles in condensin interactions with chromatin.

Pooled assembly of marine metagenomic datasets: enriching annotation through chimerism.

MOTIVATION: Despite advances in high-throughput sequencing, marine metagenomic samples remain largely opaque. A typical sample contains billions of microbial organisms from thousands of genomes and quadrillions of DNA base pairs. Its derived metagenomic dataset underrepresents this complexity by orders of magnitude because of the sparseness and shortness of sequencing reads. Read shortness and sequencing errors pose a major challenge to accurate species and functional annotation. This includes distinguishing known from novel species. Often the majority of reads cannot be annotated and thus cannot help our interpretation of the sample. RESULTS: Here, we demonstrate quantitatively how careful assembly of marine metagenomic reads within, but also across, datasets can alleviate this problem. For 10 simulated datasets, each with species complexity modeled on a real counterpart, chimerism remained within the same species for most contigs (97%). For 42 real pyrosequencing ('454') datasets, assembly increased the proportion of annotated reads, and even more so when datasets were pooled, by on average 1.6% (max 6.6%) for species, 9.0% (max 28.7%) for Pfam protein domains and 9.4% (max 22.9%) for PANTHER gene families. Our results outline exciting prospects for data sharing in the metagenomics community. While chimeric sequences should be avoided in other areas of metagenomics (e.g. biodiversity analyses), conservative pooled assembly is advantageous for annotation specificity and sensitivity. Intriguingly, our experiment also found potential prospects for (low-cost) discovery of new species in 'old' data. CONTACT: dgerloff@ffame.org SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

MaGnET: Malaria Genome Exploration Tool.

SUMMARY: The Malaria Genome Exploration Tool (MaGnET) is a software tool enabling intuitive 'exploration-style' visualization of functional genomics data relating to the malaria parasite, Plasmodium falciparum. MaGnET provides innovative integrated graphic displays for different datasets, including genomic location of genes, mRNA expression data, protein-protein interactions and more. Any selection of genes to explore made by the user is easily carried over between the different viewers for different datasets, and can be changed interactively at any point (without returning to a search). AVAILABILITY AND IMPLEMENTATION: Free online use (Java Web Start) or download (Java application archive and MySQL database; requires local MySQL installation) at http://malariagenomeexplorer.org CONTACT: joanna.sharman@ed.ac.uk or dgerloff@ffame.org SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Protein domains of unknown function are essential in bacteria.

UNLABELLED: More than 20% of all protein domains are currently annotated as "domains of unknown function" (DUFs). About 2,700 DUFs are found in bacteria compared with just over 1,500 in eukaryotes. Over 800 DUFs are shared between bacteria and eukaryotes, and about 300 of these are also present in archaea. A total of 2,786 bacterial Pfam domains even occur in animals, including 320 DUFs. Evolutionary conservation suggests that many of these DUFs are important. Here we show that 355 essential proteins in 16 model bacterial species contain 238 DUFs, most of which represent single-domain proteins, clearly establishing the biological essentiality of DUFs. We suggest that experimental research should focus on conserved and essential DUFs (eDUFs) for functional analysis given their important function and wide taxonomic distribution, including bacterial pathogens. IMPORTANCE: The functional units of proteins are domains. Typically, each domain has a distinct structure and function. Genomes encode thousands of domains, and many of the domains have no known function (domains of unknown function [DUFs]). They are often ignored as of little relevance, given that many of them are found in only a few genomes. Here we show that many DUFs are essential DUFs (eDUFs) based on their presence in essential proteins. We also show that eDUFs are often essential even if they are found in relatively few genomes. However, in general, more common DUFs are more often essential than rare DUFs.

BRCT domains: A little more than kin, and less than kind.

BRCT domains are versatile protein modular domains found as single units or as multiple copies in more than 20 different proteins in the human genome. Interestingly, most BRCT-containing proteins function in the same biological process, the DNA damage response network, but show specificity in their molecular interactions. BRCT domains have been found to bind a wide array of ligands from proteins, phosphorylated linear motifs, and DNA. Here we discuss the biology of BRCT domains and how a domain-centric analysis can aid in the understanding of signal transduction events in the DNA damage response network.

Cell contact-dependent outer membrane exchange in myxobacteria: genetic determinants and mechanism.

Biofilms are dense microbial communities. Although widely distributed and medically important, how biofilm cells interact with one another is poorly understood. Recently, we described a novel process whereby myxobacterial biofilm cells exchange their outer membrane (OM) lipoproteins. For the first time we report here the identification of two host proteins, TraAB, required for transfer. These proteins are predicted to localize in the cell envelope; and TraA encodes a distant PA14 lectin-like domain, a cysteine-rich tandem repeat region, and a putative C-terminal protein sorting tag named MYXO-CTERM, while TraB encodes an OmpA-like domain. Importantly, TraAB are required in donors and recipients, suggesting bidirectional transfer. By use of a lipophilic fluorescent dye, we also discovered that OM lipids are exchanged. Similar to lipoproteins, dye transfer requires TraAB function, gliding motility and a structured biofilm. Importantly, OM exchange was found to regulate swarming and development behaviors, suggesting a new role in cell-cell communication. A working model proposes TraA is a cell surface receptor that mediates cell-cell adhesion for OM fusion, in which lipoproteins/lipids are transferred by lateral diffusion. We further hypothesize that cell contact-dependent exchange helps myxobacteria to coordinate their social behaviors.

The interface of protein structure, protein biophysics, and molecular evolution.

Abstract The interface of protein structural biology, protein biophysics, molecular evolution, and molecular population genetics forms the foundations for a mechanistic understanding of many aspects of protein biochemistry. Current efforts in interdisciplinary protein modeling are in their infancy and the state-of-the art of such models is described. Beyond the relationship between amino acid substitution and static protein structure, protein function, and corresponding organismal fitness, other considerations are also discussed. More complex mutational processes such as insertion and deletion and domain rearrangements and even circular permutations should be evaluated. The role of intrinsically disordered proteins is still controversial, but may be increasingly important to consider. Protein geometry and protein dynamics as a deviation from static considerations of protein structure are also important. Protein expression level is known to be a major determinant of evolutionary rate and several considerations including selection at the mRNA level and the role of interaction specificity are discussed. Lastly, the relationship between modeling and needed high-throughput experimental data as well as experimental examination of protein evolution using ancestral sequence resurrection and in vitro biochemistry are presented, towards an aim of ultimately generating better models for biological inference and prediction.

OpenKnowledge for peer-to-peer experimentation in protein identification by MS/MS.

BACKGROUND: Traditional scientific workflow platforms usually run individual experiments with little evaluation and analysis of performance as required by automated experimentation in which scientists are being allowed to access numerous applicable workflows rather than being committed to a single one. Experimental protocols and data under a peer-to-peer environment could potentially be shared freely without any single point of authority to dictate how experiments should be run. In such environment it is necessary to have mechanisms by which each individual scientist (peer) can assess, locally, how he or she wants to be involved with others in experiments. This study aims to implement and demonstrate simple peer ranking under the OpenKnowledge peer-to-peer infrastructure by both simulated and real-world bioinformatics experiments involving multi-agent interactions. METHODS: A simulated experiment environment with a peer ranking capability was specified by the Lightweight Coordination Calculus (LCC) and automatically executed under the OpenKnowledge infrastructure. The peers such as MS/MS protein identification services (including web-enabled and independent programs) were made accessible as OpenKnowledge Components (OKCs) for automated execution as peers in the experiments. The performance of the peers in these automated experiments was monitored and evaluated by simple peer ranking algorithms. RESULTS: Peer ranking experiments with simulated peers exhibited characteristic behaviours, e.g., power law effect (a few dominant peers dominate), similar to that observed in the traditional Web. Real-world experiments were run using an interaction model in LCC involving two different types of MS/MS protein identification peers, viz., peptide fragment fingerprinting (PFF) and de novo sequencing with another peer ranking algorithm simply based on counting the successful and failed runs. This study demonstrated a novel integration and useful evaluation of specific proteomic peers and found MASCOT to be a dominant peer as judged by peer ranking. CONCLUSION: The simulated and real-world experiments in the present study demonstrated that the OpenKnowledge infrastructure with peer ranking capability can serve as an evaluative environment for automated experimentation.

BISC: binary subcomplexes in proteins database.

Binary subcomplexes in proteins database (BISC) is a new protein-protein interaction (PPI) database linking up the two communities most active in their characterization: structural biology and functional genomics researchers. The BISC resource offers users (i) a structural perspective and related information about binary subcomplexes (i.e. physical direct interactions between proteins) that are either structurally characterized or modellable entries in the main functional genomics PPI databases BioGRID, IntAct and HPRD; (ii) selected web services to further investigate the validity of postulated PPI by inspection of their hypothetical modelled interfaces. Among other uses we envision that this resource can help identify possible false positive PPI in current database records. BISC is freely available at http://bisc.cse.ucsc.edu.

Highly sensitive detection of individual HEAT and ARM repeats with HHpred and COACH.

BACKGROUND: HEAT and ARM repeats occur in a large number of eukaryotic proteins. As these repeats are often highly diverged, the prediction of HEAT or ARM domains can be challenging. Except for the most clear-cut cases, identification at the individual repeat level is indispensable, in particular for determining domain boundaries. However, methods using single sequence queries do not have the sensitivity required to deal with more divergent repeats and, when applied to proteins with known structures, in some cases failed to detect a single repeat. METHODOLOGY AND PRINCIPAL FINDINGS: Testing algorithms which use multiple sequence alignments as queries, we found two of them, HHpred and COACH, to detect HEAT and ARM repeats with greatly enhanced sensitivity. Calibration against experimentally determined structures suggests the use of three score classes with increasing confidence in the prediction, and prediction thresholds for each method. When we applied a new protocol using both HHpred and COACH to these structures, it detected 82% of HEAT repeats and 90% of ARM repeats, with the minimum for a given protein of 57% for HEAT repeats and 60% for ARM repeats. Application to bona fide HEAT and ARM proteins or domains indicated that similar numbers can be expected for the full complement of HEAT/ARM proteins. A systematic screen of the Protein Data Bank for false positive hits revealed their number to be low, in particular for ARM repeats. Double false positive hits for a given protein were rare for HEAT and not at all observed for ARM repeats. In combination with fold prediction and consistency checking (multiple sequence alignments, secondary structure prediction, and position analysis), repeat prediction with the new HHpred/COACH protocol dramatically improves prediction in the twilight zone of fold prediction methods, as well as the delineation of HEAT/ARM domain boundaries. SIGNIFICANCE: A protocol is presented for the identification of individual HEAT or ARM repeats which is straightforward to implement. It provides high sensitivity at a low false positive rate and will therefore greatly enhance the accuracy of predictions of HEAT and ARM domains.

Welcome to Automated Experimentation: a new open access journal.

Modern experimental science provides more opportunities for yet larger series of experiments. Demand for experimental results also has become more diverse, requiring results that have direct connections to systems outside the laboratory. With this has come an ability to automate many areas of experimental science, not only the experiments themselves but also the larger processes that contribute to experimentation and analysis more broadly. As automated experimentation becomes more widely used and understood, we launch this journal to provide a proper publication channel for this new breed of interdisciplinary research as well as a bridge to all significant groundwork research that would facilitate possible automated experimentation. With this in mind, we are interested in publishing all kinds of research into scientific experimentation, including research where the potential for automation is at proof or concept or early deployment stage.