Uncovering Endolysins against Methicillin-Resistant Staphylococcus aureus Using a Microbial Single-Cell Genome Database.

Endolysins, peptidoglycan hydrolases derived from bacteriophages (phages), are being developed as a promising alternative to conventional antibiotics. To obtain highly active endolysins, a diverse library of these endolysins is vital. We propose here microbial single-cell genome sequencing as an efficient tool to discover dozens of previously unknown endolysins, owing to its culture-independent sequencing method. As a proof of concept, we analyzed and recovered endolysin genes within prophage regions of Staphylococcus single-amplified genomes in human skin microbiome samples. We constructed a library of chimeric endolysins by shuffling domains of the natural endolysins and performed high-throughput screening against Staphylococcus aureus. One of the lead endolysins, bbst1027, exhibited desirable antimicrobial properties, such as rapid bactericidal activity, no detectable resistance development, and in vivo efficacy. We foresee that this endolysin discovery pipeline is in principle applicable to any bacterial target and boost the development of novel antimicrobial agents.

Generative aptamer discovery using RaptGen.

Nucleic acid aptamers are generated by an in vitro molecular evolution method known as systematic evolution of ligands by exponential enrichment (SELEX). Various candidates are limited by actual sequencing data from an experiment. Here we developed RaptGen, which is a variational autoencoder for in silico aptamer generation. RaptGen exploits a profile hidden Markov model decoder to represent motif sequences effectively. We showed that RaptGen embedded simulation sequence data into low-dimensional latent space on the basis of motif information. We also performed sequence embedding using two independent SELEX datasets. RaptGen successfully generated aptamers from the latent space even though they were not included in high-throughput sequencing. RaptGen could also generate a truncated aptamer with a short learning model. We demonstrated that RaptGen could be applied to activity-guided aptamer generation according to Bayesian optimization. We concluded that a generative method by RaptGen and latent representation are useful for aptamer discovery.

RaptRanker: in silico RNA aptamer selection from HT-SELEX experiment based on local sequence and structure information.

Aptamers are short single-stranded RNA/DNA molecules that bind to specific target molecules. Aptamers with high binding-affinity and target specificity are identified using an in vitro procedure called high throughput systematic evolution of ligands by exponential enrichment (HT-SELEX). However, the development of aptamer affinity reagents takes a considerable amount of time and is costly because HT-SELEX produces a large dataset of candidate sequences, some of which have insufficient binding-affinity. Here, we present RNA aptamer Ranker (RaptRanker), a novel in silico method for identifying high binding-affinity aptamers from HT-SELEX data by scoring and ranking. RaptRanker analyzes HT-SELEX data by evaluating the nucleotide sequence and secondary structure simultaneously, and by ranking according to scores reflecting local structure and sequence frequencies. To evaluate the performance of RaptRanker, we performed two new HT-SELEX experiments, and evaluated binding affinities of a part of sequences that include aptamers with low binding-affinity. In both datasets, the performance of RaptRanker was superior to Frequency, Enrichment and MPBind. We also confirmed that the consideration of secondary structures is effective in HT-SELEX data analysis, and that RaptRanker successfully predicted the essential subsequence motifs in each identified sequence.

NMR monitoring of the SELEX process to confirm enrichment of structured RNA.

RNA aptamers are RNA molecules that bind to a target molecule with high affinity and specificity using uniquely-folded tertiary structures. RNA aptamers are selected from an RNA pool typically comprising up to 1015 different sequences generated by iterative steps of selection and amplification known as Systematic Evolution of Ligands by EXponential enrichment (SELEX). Over several rounds of SELEX, the diversity of the RNA pool decreases and the aptamers are enriched. Hence, monitoring of the enrichment of these RNA pools is critical for the successful selection of aptamers, and several methods for monitoring them have been developed. In this study, we measured one-dimensional imino proton NMR spectra of RNA pools during SELEX. The spectrum of the initial RNA pool indicates that the RNAs adopt tertiary structures. The structural diversity of the RNA pools was shown to depend highly on the design of the primer-binding sequence. Furthermore, we demonstrate that enrichment of RNA aptamers can be monitored using NMR. The RNA pools can be recovered from the NMR tube after measurement of NMR spectra. We also can monitor target binding in the NMR tubes. Thus, we propose using NMR to monitor the enrichment of structured aptamers during the SELEX process.

Distinct functions of the two specificity determinants in replication initiation of plasmids ColE2-P9 and ColE3-CA38.

The plasmid ColE2-P9 Rep protein specifically binds to the cognate replication origin to initiate DNA replication. The replicons of the plasmids ColE2-P9 and ColE3-CA38 are closely related, although the actions of the Rep proteins on the origins are specific to the plasmids. The previous chimera analysis identified two regions, regions A and B, in the Rep proteins and two sites, alpha and beta, in the origins as specificity determinants and showed that when each component of the region A-site alpha pair and the region B-site beta pair is derived from the same plasmid, plasmid DNA replication is efficient. It is also indicated that the replication specificity is mainly determined by region A and site alpha. By using an electrophoretic mobility shift assay, we demonstrated that region B and site beta play a critical role for stable Rep protein-origin binding and, furthermore, that 284-Thr in this region of the ColE2 Rep protein and the corresponding 293-Trp of the ColE3 Rep protein mainly determine the Rep-origin binding specificity. On the other hand, region A and site alpha were involved in the efficient unwinding of several nucleotide residues around site alpha, although they were not involved in the stable binding of the Rep protein to the origin. Finally, we discussed how the action of the Rep protein on the origin involving these specificity determinants leads to the plasmid-specific replication initiation.

Characterization of the ColE2-like replicon of plasmid pTT8 from Thermus thermophilus.

We identified the 1.6-kb region of Thermus thermophilus plasmid pTT8 capable of autonomous replication, which shows a significant sequence similarity to the replicon regions of the ColE2-related plasmids. We showed the requirement of DNA polymerase I for pTT8 replication. The putative rep gene coding for the replication initiator protein, Rep, similar to those of the ColE2-related plasmids was cloned into an expression vector. The 6xHis-Rep protein expressed in Escherichia coli was successfully purified by stepwise denaturing with urea and refolding in the presence of glycerol on Ni-resin. We identified the nucleotide sequence recognized by the pTT8 Rep protein by the SELEX experiment using the purified protein, and proposed the existence of the third origin of pTT8 replication different from those predicted previously.

The ColE2-P9 Rep protein binds to the origin DNA as a monomer.

The Rep proteins of some plasmid replicons have two functions. Dimers bind to the operator sequences acting as auto-repressors, whereas monomers bind to the iterons to initiate replication of DNA. The ColE2 Rep proteins are present mostly in a dimeric form with some multimers larger than dimers in solution, while the form of Rep binding to Ori is not known. We used an EMSA-based method to determine the molecular weight of Rep in the Rep-Ori complex. The result suggested that Rep binds to Ori as a monomer. In addition, the result of EMSA using the Rep protein fused with the maltose binding protein and the His6-tag also supported this conclusion. We proposed that dimerization of Rep might probably be involved in keeping the copy number of the ColE2 plasmid at the normal low level by limiting the amount of active monomeric forms of Rep in the host cell.