Comparative genomics of Burkholderia multivorans, a ubiquitous pathogen with a highly conserved genomic structure

The natural environment serves as a reservoir of opportunistic pathogens. A well-established method for studying the epidemiology of such opportunists is multilocus sequence typing, which in many cases has defined strains predisposed to causing infection. Burkholderia multivorans is an important pathogen in people with cystic fibrosis (CF) and its epidemiology suggests that strains are acquired from non-human sources such as the natural environment. This raises the central question of whether the isolation source (CF or environment) or the multilocus sequence type (ST) of B. multivorans better predicts their genomic content and functionality. We identified four pairs of B. multivorans isolates, representing distinct STs and consisting of one CF and one environmental isolate each. All genomes were sequenced using the PacBio SMRT sequencing technology, which resulted in eight high-quality B. multivorans genome assemblies. The present study demonstrated that the genomic structure of the examined B. multivorans STs is highly conserved and that the B. multivorans genomic lineages are defined by their ST. Orthologous protein families were not uniformly distributed among chromosomes, with core orthologs being enriched on the primary chromosome and ST-specific orthologs being enriched on the second and third chromosome. The ST-specific orthologs were enriched in genes involved in defense mechanisms and secondary metabolism, corroborating the strain-specificity of these virulence characteristics. Finally, the same B. multivorans genomic lineages occur in both CF and environmental samples and on different continents, demonstrating their ubiquity and evolutionary persistence

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

MODBASE (http://salilab.org/modbase) is a database of annotated comparative protein structure models. The models are calculated by MODPIPE, an automated modeling pipeline that relies primarily on MODELLER for fold assignment, sequence-structure alignment, model building and model assessment (http:/salilab.org/modeller). MODBASE currently contains 5,152,695 reliable models for domains in 1,593,209 unique protein sequences; only models based on statistically significant alignments and/or models assessed to have the correct fold are included. MODBASE also allows users to calculate comparative models on demand, through an interface to the MODWEB modeling server (http://salilab.org/modweb). Other resources integrated with MODBASE include databases of multiple protein structure alignments (DBAli), structurally defined ligand binding sites (LIGBASE), predicted ligand binding sites (AnnoLyze), structurally defined binary domain interfaces (PIBASE) and annotated single nucleotide polymorphisms and somatic mutations found in human proteins (LS-SNP, LS-Mut). MODBASE models are also available through the Protein Model Portal (http://www.proteinmodelportal.org/)

FFAS server: novel features and applications.

The Fold and Function Assignment System (FFAS) server [Jaroszewski et al. (2005) FFAS03: a server for profile-profile sequence alignments. Nucleic Acids Research, 33, W284-W288] implements the algorithm for protein profile-profile alignment introduced originally in [Rychlewski et al. (2000) Comparison of sequence profiles. Strategies for structural predictions using sequence information. Protein Science: a Publication of the Protein Society, 9, 232-241]. Here, we present updates, changes and novel functionality added to the server since 2005 and discuss its new applications. The sequence database used to calculate sequence profiles was enriched by adding sets of publicly available metagenomic sequences. The profile of a user's protein can now be compared with ∼20 additional profile databases, including several complete proteomes, human proteins involved in genetic diseases and a database of microbial virulence factors. A newly developed interface uses a system of tabs, allowing the user to navigate multiple results pages, and also includes novel functionality, such as a dotplot graph viewer, modeling tools, an improved 3D alignment viewer and links to the database of structural similarities. The FFAS server was also optimized for speed: running times were reduced by an order of magnitude. The FFAS server, http://ffas.godziklab.org, has no log-in requirement, albeit there is an option to register and store results in individual, password-protected directories. Source code and Linux executables for the FFAS program are available for download from the FFAS server

Atomic Resolution Structure of the Oncolytic Parvovirus LuIII by Electron Microscopy and 3D Image Reconstruction.

LuIII, a protoparvovirus pathogenic to rodents, replicates in human mitotic cells, making it applicable for use to kill cancer cells. This virus group includes H-1 parvovirus (H-1PV) and minute virus of mice (MVM). However, LuIII displays enhanced oncolysis compared to H-1PV and MVM, a phenotype mapped to the major capsid viral protein 2 (VP2). This suggests that within LuIII VP2 are determinants for improved tumor lysis. To investigate this, the structure of the LuIII virus-like-particle was determined using single particle cryo-electron microscopy and image reconstruction to 3.17 Å resolution, and compared to the H-1PV and MVM structures. The LuIII VP2 structure, ordered from residue 37 to 587 (C-terminal), had the conserved VP topology and capsid morphology previously reported for other protoparvoviruses. This includes a core β-barrel and α-helix A, a depression at the icosahedral 2-fold and surrounding the 5-fold axes, and a single protrusion at the 3-fold axes. Comparative analysis identified surface loop differences among LuIII, H-1PV, and MVM at or close to the capsid 2- and 5-fold symmetry axes, and the shoulder of the 3-fold protrusions. The 2-fold differences cluster near the previously identified MVM sialic acid receptor binding pocket, and revealed potential determinants of protoparvovirus tumor tropism

BOOL-AN: A method for comparative sequence analysis and phylogenetic reconstruction

A novel discrete mathematical approach is proposed as an additional tool for molecular systematics which does not require prior statistical assumptions concerning the evolutionary process. The method is based on algorithms generating mathematical representations directly from DNA/RNA or protein sequences, followed by the output of numerical (scalar or vector) and visual characteristics (graphs). The binary encoded sequence information is transformed into a compact analytical form, called the Iterative Canonical Form (or ICF) of Boolean functions, which can then be used as a generalized molecular descriptor. The method provides raw vector data for calculating different distance matrices, which in turn can be analyzed by neighbor-joining or UPGMA to derive a phylogenetic tree, or by principal coordinates analysis to get an ordination scattergram. The new method and the associated software for inferring phylogenetic trees are called the Boolean analysis or BOOL-AN

Atomic structures and deletion mutant reveal different capsid-binding patterns and functional significance of tegument protein pp150 in murine and human cytomegaloviruses with implications for therapeutic development.

Cytomegalovirus (CMV) infection causes birth defects and life-threatening complications in immunosuppressed patients. Lack of vaccine and need for more effective drugs have driven widespread ongoing therapeutic development efforts against human CMV (HCMV), mostly using murine CMV (MCMV) as the model system for preclinical animal tests. The recent publication (Yu et al., 2017, DOI: 10.1126/science.aam6892) of an atomic model for HCMV capsid with associated tegument protein pp150 has infused impetus for rational design of novel vaccines and drugs, but the absence of high-resolution structural data on MCMV remains a significant knowledge gap in such development efforts. Here, by cryoEM with sub-particle reconstruction method, we have obtained the first atomic structure of MCMV capsid with associated pp150. Surprisingly, the capsid-binding patterns of pp150 differ between HCMV and MCMV despite their highly similar capsid structures. In MCMV, pp150 is absent on triplex Tc and exists as a "Λ"-shaped dimer on other triplexes, leading to only 260 groups of two pp150 subunits per capsid in contrast to 320 groups of three pp150 subunits each in a "Δ"-shaped fortifying configuration. Many more amino acids contribute to pp150-pp150 interactions in MCMV than in HCMV, making MCMV pp150 dimer inflexible thus incompatible to instigate triplex Tc-binding as observed in HCMV. While pp150 is essential in HCMV, our pp150-deletion mutant of MCMV remained viable though with attenuated infectivity and exhibiting defects in retaining viral genome. These results thus invalidate targeting pp150, but lend support to targeting capsid proteins, when using MCMV as a model for HCMV pathogenesis and therapeutic studies