Orientation-dependent backbone-only residue pair scoring functions for fixed backbone protein design

<p>Abstract</p> <p>Background</p> <p>Empirical scoring functions have proven useful in protein structure modeling. Most such scoring functions depend on protein side chain conformations. However, backbone-only scoring functions do not require computationally intensive structure optimization and so are well suited to protein design, which requires fast score evaluation. Furthermore, scoring functions that account for the distinctive relative position and orientation preferences of residue pairs are expected to be more accurate than those that depend only on the separation distance.</p> <p>Results</p> <p>Residue pair scoring functions for fixed backbone protein design were derived using only backbone geometry. Unlike previous studies that used spherical harmonics to fit 2D angular distributions, Gaussian Mixture Models were used to fit the full 3D (position only) and 6D (position and orientation) distributions of residue pairs. The performance of the 1D (residue separation only), 3D, and 6D scoring functions were compared by their ability to identify correct threading solutions for a non-redundant benchmark set of protein backbone structures. The threading accuracy was found to steadily increase with increasing dimension, with the 6D scoring function achieving the highest accuracy. Furthermore, the 3D and 6D scoring functions were shown to outperform side chain-dependent empirical potentials from three other studies. Next, two computational methods that take advantage of the speed and pairwise form of these new backbone-only scoring functions were investigated. The first is a procedure that exploits available sequence data by averaging scores over threading solutions for homologs. This was evaluated by applying it to the challenging problem of identifying interacting transmembrane alpha-helices and found to further improve prediction accuracy. The second is a protein design method for determining the optimal sequence for a backbone structure by applying Belief Propagation optimization using the 6D scoring functions. The sensitivity of this method to backbone structure perturbations was compared with that of fixed-backbone all-atom modeling by determining the similarities between optimal sequences for two different backbone structures within the same protein family. The results showed that the design method using 6D scoring functions was more robust to small variations in backbone structure than the all-atom design method.</p> <p>Conclusions</p> <p>Backbone-only residue pair scoring functions that account for all six relative degrees of freedom are the most accurate and including the scores of homologs further improves the accuracy in threading applications. The 6D scoring function outperformed several side chain-dependent potentials while avoiding time-consuming and error prone side chain structure prediction. These scoring functions are particularly useful as an initial filter in protein design problems before applying all-atom modeling.</p

Genome Erosion in a Nitrogen-Fixing Vertically Transmitted Endosymbiotic Multicellular Cyanobacterium

Background: An ancient cyanobacterial incorporation into a eukaryotic organism led to the evolution of plastids (chloroplasts) and subsequently to the origin of the plant kingdom. The underlying mechanism and the identities of the partners in this monophyletic event remain elusive. Methodology/Principal Findings: To shed light on this evolutionary process, we sequenced the genome of a cyanobacterium residing extracellularly in an endosymbiosis with a plant, the water-fern Azolla filiculoides Lam. This symbiosis was selected as it has characters which make it unique among extant cyanobacterial plant symbioses: the cyanobacterium lacks autonomous growth and is vertically transmitted between plant generations. Our results reveal features of evolutionary significance. The genome is in an eroding state, evidenced by a large proportion of pseudogenes (31.2%) and a high frequency of transposable elements (,600) scattered throughout the genome. Pseudogenization is found in genes such as the replication initiator dnaA and DNA repair genes, considered essential to free-living cyanobacteria. For some functional categories of genes pseudogenes are more prevalent than functional genes. Loss of function is apparent even within the ‘core’ gene categories of bacteria, such as genes involved in glycolysis and nutrient uptake. In contrast, serving as a critical source of nitrogen for the host, genes related to metabolic processes such as cell differentiation and nitrogen-fixation are well preserved. Conclusions/Significance: This is the first finding of genome degradation in a plant symbiont and phenotypically complex cyanobacterium and one of only a few extracellular endosymbionts described showing signs of reductive genome evolution. Our findings suggest an ongoing selective streamlining of this cyanobacterial genome which has resulted in an organism devoted to nitrogen fixation and devoid of autonomous growth. The cyanobacterial symbiont of Azolla can thus be considered at the initial phase of a transition from free-living organism to a nitrogen-fixing plant entity, a transition process which may mimic what drove the evolution of chloroplasts from a cyanobacterial ancestor

Genome fluctuations in cyanobacteria reflect evolutionary, developmental and adaptive traits

<p>Abstract</p> <p>Background</p> <p>Cyanobacteria belong to an ancient group of photosynthetic prokaryotes with pronounced variations in their cellular differentiation strategies, physiological capacities and choice of habitat. Sequencing efforts have shown that genomes within this phylum are equally diverse in terms of size and protein-coding capacity. To increase our understanding of genomic changes in the lineage, the genomes of 58 contemporary cyanobacteria were analysed for shared and unique orthologs.</p> <p>Results</p> <p>A total of 404 protein families, present in all cyanobacterial genomes, were identified. Two of these are unique to the phylum, corresponding to an AbrB family transcriptional regulator and a gene that escapes functional annotation although its genomic neighbourhood is conserved among the organisms examined. The evolution of cyanobacterial genome sizes involves a mix of gains and losses in the clade encompassing complex cyanobacteria, while a single event of reduction is evident in a clade dominated by unicellular cyanobacteria. Genome sizes and gene family copy numbers evolve at a higher rate in the former clade, and multi-copy genes were predominant in large genomes. Orthologs unique to cyanobacteria exhibiting specific characteristics, such as filament formation, heterocyst differentiation, diazotrophy and symbiotic competence, were also identified. An ancestral character reconstruction suggests that the most recent common ancestor of cyanobacteria had a genome size of approx. 4.5 Mbp and 1678 to 3291 protein-coding genes, 4%-6% of which are unique to cyanobacteria today.</p> <p>Conclusions</p> <p>The different rates of genome-size evolution and multi-copy gene abundance suggest two routes of genome development in the history of cyanobacteria. The expansion strategy is driven by gene-family enlargment and generates a broad adaptive potential; while the genome streamlining strategy imposes adaptations to highly specific niches, also reflected in their different functional capacities. A few genomes display extreme proliferation of non-coding nucleotides which is likely to be the result of initial expansion of genomes/gene copy number to gain adaptive potential, followed by a shift to a life-style in a highly specific niche (e.g. symbiosis). This transition results in redundancy of genes and gene families, leading to an increase in junk DNA and eventually to gene loss. A few orthologs can be correlated with specific phenotypes in cyanobacteria, such as filament formation and symbiotic competence; these constitute exciting exploratory targets.</p

Genomic epidemiology reveals multiple introductions of Zika virus into the United States

Zika virus (ZIKV) is causing an unprecedented epidemic linked to severe congenital abnormalities. In July 2016, mosquito-borne ZIKV transmission was reported in the continental United States; since then, hundreds of locally acquired infections have been reported in Florida. To gain insights into the timing, source, and likely route(s) of ZIKV introduction, we tracked the virus from its first detection in Florida by sequencing ZIKV genomes from infected patients and Aedes aegypti mosquitoes. We show that at least 4 introductions, but potentially as many as 40, contributed to the outbreak in Florida and that local transmission is likely to have started in the spring of 2016-several months before its initial detection. By analysing surveillance and genetic data, we show that ZIKV moved among transmission zones in Miami. Our analyses show that most introductions were linked to the Caribbean, a finding corroborated by the high incidence rates and traffic volumes from the region into the Miami area. Our study provides an understanding of how ZIKV initiates transmission in new regions

Technologies of sleep research

Sleep is investigated in many different ways, many different species and under many different circumstances. Modern sleep research is a multidisciplinary venture. Therefore, this review cannot give a complete overview of all techniques used in sleep research and sleep medicine. What it will try to do is to give an overview of widely applied techniques and exciting new developments. Electroencephalography has been the backbone of sleep research and sleep medicine since its first application in the 1930s. The electroencephalogram is still used but now combined with many different techniques monitoring body and brain temperature, changes in brain and blood chemistry, or changes in brain functioning. Animal research has been very important for progress in sleep research and sleep medicine. It provides opportunities to investigate the sleeping brain in ways not possible in healthy volunteers. Progress in genomics has brought new insights in sleep regulation, the best example being the discovery of hypocretin/orexin deficiency as the cause of narcolepsy. Gene manipulation holds great promise for the future since it is possible not only to investigate the functions of different genes under normal conditions, but also to mimic human pathology in much greater detail

Search for direct top squark pair production in final states with two leptons in √s=13 TeV pp collisions with the ATLAS detector

The results of a search for direct pair production of top squarks in events with two opposite-charge leptons (electrons or muons) are reported, using 36.1 fb−1 of integrated luminosity from proton–proton collisions at √s=13 TeV collected by the ATLAS detector at the Large Hadron Collider. To cover a range of mass differences between the top squark t~ and lighter supersymmetric particles, four possible decay modes of the top squark are targeted with dedicated selections: the decay t~→bχ~1± into a b-quark and the lightest chargino with χ~1±→Wχ~10 , t~→tχ~10 into an on-shell top quark and the lightest neutralino, the three-body decay t~→bWχ~10 and the four-body decay t~→bℓνχ~10. No significant excess of events is observed above the Standard Model background for any selection, and limits on top squarks are set as a function of the t~ and χ~01 masses. The results exclude at 95% confidence level t~ masses up to about 720 GeV, extending the exclusion region of supersymmetric parameter space covered by previous searches

Measurements of long-range azimuthal anisotropies and associated Fourier coefficients for pp collisions at √s=5.02 and 13 TeV and p+Pb collisions at √sNN=5.02 TeV with the ATLAS detector

ATLAS measurements of two-particle correlations are presented for √s=5.02 and 13 TeV ppcollisions and for √sNN=5.02 TeV p+Pb collisions at the LHC. The correlation functions are measured as a function of relative azimuthal angle Δϕ, and pseudorapidity separation Δη, using charged particles detected within the pseudorapidity interval |η|2, is studied using a template fitting procedure to remove a “back-to-back” contribution to the correlation function that primarily arises from hard-scattering processes. In addition to the elliptic, cos (2Δϕ), modulation observed in a previous measurement, the pp correlation functions exhibit significant cos (3Δϕ) and cos (4Δϕ) modulation. The Fourier coefficients vn, n associated with the cos (nΔϕ) modulation of the correlation functions for n=2–4 are measured as a function of charged-particle multiplicity and charged-particle transverse momentum. The Fourier coefficients are observed to be compatible with cos (nϕ) modulation of per-event single-particle azimuthal angle distributions. The single-particle Fourier coefficients vn are measured as a function of charged-particle multiplicity, and charged-particle transverse momentum for n=2–4. The integrated luminosities used in this analysis are, 64nb−1 for the √s=13 TeV pp data, 170 nb−1 for the √ s = 5.02 TeV pp data, and 28 nb−1 for the √sNN = 5.02 TeV p+Pb data

Search for supersymmetry at √s=13 TeV in final states with jets and two same-sign leptons or three leptons with the ATLAS detector

A search for strongly produced supersymmetric particles is conducted using signatures involving multiple energetic jets and either two isolated leptons (e or μ μ) with the same electric charge or at least three isolated leptons. The search also utilises b-tagged jets, missing transverse momentum and other observables to extend its sensitivity. The analysis uses a data sample of proton–proton collisions at √s=13 TeV recorded with the ATLAS detector at the Large Hadron Collider in 2015 corresponding to a total integrated luminosity of 3.2 fb −1. No significant excess over the Standard Model expectation is observed. The results are interpreted in several simplified supersymmetric models and extend the exclusion limits from previous searches. In the context of exclusive production and simplified decay modes, gluino masses are excluded at 95% 95% confidence level up to 1.1–1.3 TeV for light neutralinos (depending on the decay channel), and bottom squark masses are also excluded up to 540 GeV. In the former scenarios, neutralino masses are also excluded up to 550–850 GeV for gluino masses around 1 TeV

Search for long-lived charginos based on a disappearing-track signature in pp collisions at s √ =13 s=13 TeV with the ATLAS detector

This paper presents a search for direct electroweak gaugino or gluino pair production with a chargino nearly mass-degenerate with a stable neutralino. It is based on an integrated luminosity of 36.1 fb−1 of pp collisions at s √ =13 s=13 TeV collected by the ATLAS experiment at the LHC. The final state of interest is a disappearing track accompanied by at least one jet with high transverse momentum from initial-state radiation or by four jets from the gluino decay chain. The use of short track segments reconstructed from the innermost tracking layers significantly improves the sensitivity to short chargino lifetimes. The results are found to be consistent with Standard Model predictions. Exclusion limits are set at 95% confidence level on the mass of charginos and gluinos for different chargino lifetimes. For a pure wino with a lifetime of about 0.2 ns, chargino masses up to 460 GeV are excluded. For the strong production channel, gluino masses up to 1.65 TeV are excluded assuming a chargino mass of 460 GeV and lifetime of 0.2 ns

Searches for the Zγ decay mode of the Higgs boson and for new high-mass resonances in pp collisions at √s=13 TeV with the ATLAS detector

This article presents searches for the Zγ decay of the Higgs boson and for narrow high-mass resonances decaying to Zγ, exploiting Z boson decays to pairs of electrons or muons. The data analysis uses 36.1 fb−1 of pp collisions at √s=13 recorded by the ATLAS detector at the CERN Large Hadron Collider. The data are found to be consistent with the expected Standard Model background. The observed (expected — assuming Standard Model pp → H → Zγ production and decay) upper limit on the production cross section times the branching ratio for pp → H → Zγ is 6.6. (5.2) times the Standard Model prediction at the 95% confidence level for a Higgs boson mass of 125.09 GeV. In addition, upper limits are set on the production cross section times the branching ratio as a function of the mass of a narrow resonance between 250 GeV and 2.4 TeV, assuming spin-0 resonances produced via gluon-gluon fusion, and spin-2 resonances produced via gluon-gluon or quark-antiquark initial states. For high-mass spin-0 resonances, the observed (expected) limits vary between 88 fb (61 fb) and 2.8 fb (2.7 fb) for the mass range from 250 GeV to 2.4 TeV at the 95% confidence level