SPIN2: Predicting sequence profiles from protein structures using deep neural networks

Designing protein sequences that can fold into a given structure is a well‐known inverse protein‐folding problem. One important characteristic to attain for a protein design program is the ability to recover wild‐type sequences given their native backbone structures. The highest average sequence identity accuracy achieved by current protein‐design programs in this problem is around 30%, achieved by our previous system, SPIN. SPIN is a program that predicts sequences compatible with a provided structure using a neural network with fragment‐based local and energy‐based nonlocal profiles. Our new model, SPIN2, uses a deep neural network and additional structural features to improve on SPIN. SPIN2 achieves over 34% in sequence recovery in 10‐fold cross‐validation and independent tests, a 4% improvement over the previous version. The sequence profiles generated from SPIN2 are expected to be useful for improving existing fold recognition and protein design techniques. SPIN2 is available at http://sparks-lab.org .     

Marmoset phylogenetics, conservation perspectives, and evolution of the mtDNA control region

Marmosets (genus Callithrix) are a diverse group of platyrrhine primates with 13-15 purported taxa, many of them considered endangered. Morphological analyses constitute most of the basis for recognition of these forms as distinct taxa. The purpose of this study was to provide a molecular view, based on mitochondrial control region sequences, of the evolutionary history of the marmosets, concomitant with a molecular phylogenetic perspective on species diversity within the group. An additional purpose was to provide the first comparative examination of a complete New World monkey control region sequence with those of other mammals. The phylogenetic analyses provide convincing support for a split between the Atlantic forest and Amazonian marmosets, with the inclusion of the pygmy marmoset (Cebuella pygmaea) at the base of the Amazonian clade. The earliest branch of the Atlantic forest group was C. aurita. In the Amazonian group, the analyses do not support the recognition of C. humeralifer and the recently described C mauesi as distinct taxa. They do, however, support a clear distinction between C. argentata and a strongly supported mixed clade of C. humeralifer and C. mauesi. In the Atlantic forest group, the phylogenetic tree suggests mixing between C. penicillata, C. kuhli, and possibly C. jacchus. Most of the sequence features characteristic of other mammal control regions were also evident in marmosets, with the exception that conserved sequence blocks (CSBs) 2 and 3 were not clearly identifiable. Tandem repeat units often associated with heteroplasmy in a variety of other mammals were not evident in the marmoset sequences.      

Persistent nuclear ribosomal DNA sequence polymorphism in the Amelanchier agamic complex (Rosaceae)

Individual plants of several Amelanchier taxa contain many polymorphic nucleotide sites in the internal transcribed spacers (ITS) of nuclear ribosomal DNA (nrDNA). This polymorphism is unusual because it is not recent in origin and thus has resisted homogenization by concerted evolution. Amelanchier ITS sequence polymorphism is hypothesized to be the result of gene flow between two major North American clades resolved by phylogenetic analysis of ITS sequences. Western North American species plus A. humilis and A. sanguinea of eastern North America form one clade (A), and the remaining eastern North American Amelanchier make up clade B. Five eastern North American taxa are polymorphic at many of the nucleotide sites where clades A and B have diverged and are thought to be of hybrid origin, with A. humilis or A. sanguinea as one parent and various members of clade B as the other parent. Morphological evidence suggests that A. humilis is one of the parents of one of the polymorphic taxa, a microspecies that we refer to informally as A. "erecta." Sequences of 21 cloned copies of the ITS1- 5.8S gene-ITS2 region from one A. "erecta" individual are identical to A. humilis sequence or to the clade B consensus sequence, or they are apparent recombinants of A. humilis and clade B ITS repeats. Amelanchier "erecta" and another polymorphic taxon are suspected to be relatively old because both grow several hundred kilometers beyond the range of one of their parents. ITS sequence polymorphisms have apparently persisted in these two taxa perhaps because of polyploidy and/or agamospermy (asexual seed production), which are prevalent in the genus.      

Ground wētā in vines of the Awatere Valley,Marlborough: biology,density and distribution

The endemic New Zealand ground wētā (Hemiandrus sp. ‘promontorius’) has a Naturally Uncommon conservation status. This is because of the paucity of information on its density and distribution. Here, the biology, density and distribution of a population of this wētā found in and around vineyards in the Awatere Valley, Marlborough was studied. Wētā density was assessed in vineyards, paddocks and shrublands in this valley. Soil moisture, penetration resistance, pH and organic matter were recorded at locations with and without wētā. Wētā density in vineyards was significantly higher than in either paddocks or shrub habitats. In vineyards, the density of this insect was significantly higher under-vines than in the inter-rows. Higher numbers of this wētā were found in moist soils that required lower force to burrow. Females laid an average of 55 eggs between March and April, which hatched in September. These findings highlight the intersection between agriculture and conservation.     

Dynamical crises,multistability and the influence of the duration of immunity in a seasonally-forced model of disease transmission

Background

Highly successful strategies to make populations more resilient to infectious diseases, such as childhood vaccinations programs, may nonetheless lead to unpredictable outcomes due to the interplay between seasonal variations in transmission and a population’s immune status.

Methods

Motivated by the study of diseases such as pertussis we introduce a seasonally-forced susceptible-infectious-recovered model of disease transmission with waning and boosting of immunity. We study the system’s dynamical properties using a combination of numerical simulations and bifurcation techniques, paying particular attention to the properties of the initial condition space.

Results

We find that highly unpredictable behaviour can be triggered by changes in biologically relevant model parameters such as the duration of immunity. In the particular system we analyse — previously used in the literature to study pertussis dynamics — we identify the presence of an initial-condition landscape containing three coexisting attractors. The system’s response to interventions which perturb population immunity (e.g. vaccination "catch-up" campaigns) is therefore difficult to predict.

Conclusion

Given the increasing use of models to inform policy decisions regarding vaccine introduction and scheduling and infectious diseases intervention policy more generally, our findings highlight the importance of thoroughly investigating the dynamical properties of those models to identify key areas of uncertainty. Our findings suggest that the often stated tension between capturing biological complexity and utilising mathematically simple models is perhaps more nuanced than generally suggested. Simple dynamical models, particularly those which include forcing terms, can give rise to incredibly complex behaviour.

     

Diminished transmission of drug resistant HIV-1 variants with reduced replication capacity in a human transmission model

Background

Different patterns of drug resistance are observed in treated and therapy naïve HIV-1 infected populations. Especially the NRTI-related M184I/V variants, which are among the most frequently encountered mutations in treated patients, are underrepresented in the antiretroviral naïve population. M184I/V mutations are known to have a profound effect on viral replication and tend to revert over time in the new host. However it is debated whether a diminished transmission efficacy of HIV variants with a reduced replication capacity can also contribute to the observed discrepancy in genotypic patterns.As dendritic cells (DCs) play a pivotal role in HIV-1 transmission, we used a model containing primary human Langerhans cells (LCs) and DCs to compare the transmission efficacy M184 variants (HIV-M184V/I/T) to HIV wild type (HIV-WT). As control, we used HIV harboring the NNRTI mutation K103N (HIV-K103N) which has a minor effect on replication and is found at a similar prevalence in treated and untreated individuals.

Results

In comparison to HIV-WT, the HIV-M184 variants were less efficiently transmitted to CCR5⁺ Jurkat T cells by both LCs and DCs. The transmission rate of HIV-K103N was slightly reduced to HIV-WT in LCs and even higher than HIV-WT in DCs. Replication experiments in CCR5⁺ Jurkat T cells revealed no apparent differences in replication capacity between the mutant viruses and HIV-WT. However, viral replication in LCs and DCs was in concordance with the transmission results; replication by the HIV-M184 variants was lower than replication by HIV-WT, and the level of replication of HIV-K103N was intermediate for LCs and higher than HIV-WT for DCs.

Conclusions

Our data demonstrate that drug resistant M184-variants display a reduced replication capacity in LCs and DCs which directly impairs their transmission efficacy. As such, diminished transmission efficacy may contribute to the lower prevalence of drug resistant variants in therapy naive individuals.

     

The Salmonella typhimurium virulence plasmid complement resistance gene rck is homologous to a family of virulence-related outer membrane protein genes, including pagC and ail.

A fragment of the Salmonella typhimurium virulence plasmid containing the rck locus, when cloned in the recombinant cosmid pADE016, was shown previously to confer high-level complement resistance on both rough and smooth Escherichia coli, Salmonella minnesota, and S. typhimurium and was associated with the production of an outer membrane protein. We determined the nucleotide sequence of the fragment containing the rck locus. Mutations in the two major open reading frames confirmed that the complement resistance mediated by pADE016 was due to a single 555-bp rck gene encoding a 17-kDa outer membrane protein. Analysis of the rck gene revealed that the Rck outer membrane protein consisted of 185 amino acid residues, with a calculated postcleavage molecular mass of 17.4 kDa. Rck is homologous to a family of outer membrane proteins expressed in gram-negative bacteria, two of which have been associated with virulence-related phenotypes: PagC, required by S. typhimurium for survival in macrophages and for virulence in mice; and Ail, a product of the Yersinia enterocolitica chromosome capable of mediating bacterial adherence to and invasion of epithelial cell lines. Rck, most closely related to PagC, represents the third outer membrane protein in this five-member family with a distinct virulence-associated phenotype.