Structure-based protein design with deep learning

Since the first revelation of proteins functioning as macromolecular machines through their three dimensional structures, researchers have been intrigued by the marvelous ways the biochemical processes are carried out by proteins. The aspiration to understand protein structures has fueled extensive efforts across different scientific disciplines. In recent years, it has been demonstrated that proteins with new functionality or shapes can be designed via structure-based modeling methods, and the design strategies have combined all available information — but largely piece-by-piece — from sequence derived statistics to the detailed atomic-level modeling of chemical interactions. Despite the significant progress, incorporating data-derived approaches through the use of deep learning methods can be a game changer. In this review, we summarize current progress, compare the arc of developing the deep learning approaches with the conventional methods, and describe the motivation and concepts behind current strategies that may lead to potential future opportunities.     

An Eight Amino Acid Segment Controls Oligomerization and Preferred Conformation of the two Non-visual Arrestins

G protein coupled receptors signal through G proteins or arrestins. A long-standing mystery in the field is why vertebrates have two non-visual arrestins, arrestin-2 and arrestin-3. These isoforms are ~75% identical and 85% similar; each binds numerous receptors, and appear to have many redundant functions, as demonstrated by studies of knockout mice. We previously showed that arrestin-3 can be activated by inositol-hexakisphosphate (IP₆). IP₆ interacts with the receptor-binding surface of arrestin-3, induces arrestin-3 oligomerization, and this oligomer stabilizes the active conformation of arrestin-3. Here, we compared the impact of IP₆ on oligomerization and conformational equilibrium of the highly homologous arrestin-2 and arrestin-3 and found that these two isoforms are regulated differently. In the presence of IP₆, arrestin-2 forms “infinite” chains, where each promoter remains in the basal conformation. In contrast, full length and truncated arrestin-3 form trimers and higher-order oligomers in the presence of IP₆; we showed previously that trimeric state induces arrestin-3 activation (Chen et al., 2017). Thus, in response to IP₆, the two non-visual arrestins oligomerize in different ways in distinct conformations. We identified an insertion of eight residues that is conserved across arrestin-2 homologs, but absent in arrestin-3 that likely accounts for the differences in the IP₆ effect. Because IP₆ is ubiquitously present in cells, this suggests physiological consequences, including differences in arrestin-2/3 trafficking and JNK3 activation. The functional differences between two non-visual arrestins are in part determined by distinct modes of their oligomerization. The mode of oligomerization might regulate the function of other signaling proteins.     

What paths do advantageous alleles take during short‐term evolutionary change?

Combining experimental evolution with whole‐genome resequencing is a promising new strategy for investigating the dynamics of evolutionary change. Published studies that have resequenced laboratory‐selected populations of sexual organisms have typically focused on populations sampled at the end of an evolution experiment. These studies have attempted to associate particular alleles with phenotypic change and attempted to distinguish between different theoretical models of adaptation. However, neither the population used to initiate the experiment nor multiple time points sampled during the evolutionary trajectory are generally available for examination. In this issue of Molecular Ecology, Orozco‐terWengel et al. (2012) take a significant step forward by estimating genome‐wide allele frequencies at the start, 15 generations into and at the end of a 37‐generation Drosophila experimental evolution study. The authors identify regions of the genome that have responded to laboratory selection and describe the temporal dynamics of allele frequency change. They identify two common trajectories for putatively adaptive alleles: alleles either gradually increase in frequency throughout the entire 37 generations or alleles plateau at a new frequency by generation 15. The identification of complex trajectories of alleles under selection contributes to a growing body of literature suggesting that simple models of adaptation, whereby beneficial alleles arise and increase in frequency unimpeded until they become fixed, may not adequately describe short‐term response to selection.     

Non-random distribution of individual genetic diversity along an environmental gradient

Improving our knowledge of the links between ecology and evolution is especially critical in the actual context of global rapid environmental changes. A critical step in that direction is to quantify how variation in ecological factors linked to habitat modifications might shape observed levels of genetic variability in wild populations. Still, little is known on the factors affecting levels and distribution of genetic diversity at the individual level, despite its vital underlying role in evolutionary processes. In this study, we assessed the effects of habitat quality on population structure and individual genetic diversity of tree swallows (Tachycineta bicolor) breeding along a gradient of agricultural intensification in southern Québec, Canada. Using a landscape genetics approach, we found that individual genetic diversity was greater in poorer quality habitats. This counter-intuitive result was partly explained by the settlement patterns of tree swallows across the landscape. Individuals of higher genetic diversity arrived earlier on their breeding grounds and settled in the first available habitats, which correspond to intensive cultures. Our results highlight the importance of investigating the effects of environmental variability on individual genetic diversity, and of integrating information on landscape structure when conducting such studies.     

Male-biased sex ratios in mature damselfly populations: real or artefact?

1. There is ongoing controversy about whether biased sex ratios in diploid insect populations are real or an artefact caused by different behaviours and/or different catchability of the sexes. This was tested by monitoring two field and three semi-natural populations of the damselfly Lestes sponsa. 2. Capture–mark–recapture data showed that population size estimates were about twice as large for males as for females at both field sites. Independent estimates of the sex ratios based on total numbers of males and females captured supported the male bias. 3. Males had higher recapture probabilities than females due to longer times between successive visits in females. Because the same pattern was found in the semi-natural populations, the longer intervals in females are no artefact due to their lower detectability. 4. Theoretical models show that the strong temporary emigration of females tends, if anything, to overestimate female population sizes and that the heterogeneity of recapture probabilities observed in males tends to underestimate male population sizes. Hence, behavioural differences between the sexes do not cause an artificially male-biased sex ratio. 5. Spatial data show that operational sex ratios are male biased at the pond but become female biased in the plots further away from the shoreline; however because of the decrease in densities away from the shoreline, this does not result in a global even sex ratio. 6. Spatial data, temporary emigration patterns, and independent estimates suggest strongly that the male-biased sex ratios in mature damselfly populations are real.     

Random walks and the gas model: spacing behaviour of Grey-Cheeked Mangabeys

1. Simple mechanical models (random walks and the 'gas model') were used to investigate the movement patterns and intergroup encounter rates of a high- and low-density subpopulation of Grey-Cheeked Mangabeys ( Cercocebus albigena johnstoni, Lydeker) at two sites in the Kibale Forest, Western Uganda.
2. Random walk simulations showed that the presence of conspecific groups could act as 'social barriers' which constrained group movements, and promoted high levels of site attachment to a specific home range area.
3. Encounter rate models showed that in the low-density subpopulation (Kanyawara), intergroup encounters occurred no more frequently than expected if groups were moving randomly and independently of each other. This was in contrast to previous analyses which suggested that Grey-Cheeked Mangabeys employed a social spacing mechanism (mediated by male loud calls) in order to reduce the frequency of encounter to below chance levels.
4. Encounters in the high-density subpopulation (Ngogo) were found to occur less frequently than expected at short range (within 100 m). This was suggested to be due to the operation of a territorial mate defence strategy by males at this site.     

Habitat preference in the critically endangered yellow‐tailed woolly monkey (Lagothrix flavicauda) at La Esperanza,Peru

Habitat loss is one of the main threats to wildlife. Therefore, knowledge of habitat use and preference is essential for the design of conservation strategies and identification of priority sites for the protection of endangered species. The yellow‐tailed woolly monkey (Lagothrix flavicauda Humboldt, 1812), categorized as Critically Endangered on the IUCN Red List, is endemic to montane forests in northern Peru where its habitat is greatly threatened. We assessed how habitat use and preference in L. flavicauda are linked to forest structure and composition. The study took place near La Esperanza, in the Amazonas region, Peru. Our objective was to identify characteristics of habitat most utilized by L. flavicauda to provide information that will be useful for the selection of priority sites for conservation measures. Using presence records collected from May 2013 to February 2014 for one group of L. flavicauda, we classified the study site into three different use zones: low‐use, medium‐use, and high‐use. We assessed forest structure and composition for all use zones using 0.1 ha Gentry vegetation transects. Results show high levels of variation in plant species composition across the three use zones. Plants used as food resources had considerably greater density, dominance, and ecological importance in high‐use zones. High‐use zones presented similar structure to medium‐ and low‐use zones; thus it remains difficult to assess the influence of forest structure on habitat preference. We recommend focusing conservation efforts on areas with a similar floristic composition to the high‐use zones recorded in this study and suggest utilizing key alimentation species for reforestation efforts.     

Solution structure of the N‐terminal domain of DC‐UbP/UBTD2 and its interaction with ubiquitin

DC‐UbP/UBTD2 is a ubiquitin (Ub) domain‐containing protein first identified from dendritic cells, and is implicated in ubiquitination pathway. The solution structure and backbone dynamics of the C‐terminal Ub‐like (UbL) domain were elucidated in our previous work. To further understand the biological function of DC‐UbP, we then solved the solution structure of the N‐terminal domain of DC‐UbP (DC‐UbP_N) and studied its Ub binding properties by NMR techniques. The results show that DC‐UbP_N holds a novel structural fold and acts as a Ub‐binding domain (UBD) but with low affinity. This implies that the DC‐UbP protein, composing of a combination of both UbL and UBD domains, might play an important role in regulating protein ubiquitination and delivery of ubiquitinated substrates in eukaryotic cells.