Hypothesis-free phenotype prediction within a genetics-first framework.

Cohort-wide sequencing studies have revealed that the largest category of variants is those deemed 'rare', even for the subset located in coding regions (99% of known coding variants are seen in less than 1% of the population. Associative methods give some understanding how rare genetic variants influence disease and organism-level phenotypes. But here we show that additional discoveries can be made through a knowledge-based approach using protein domains and ontologies (function and phenotype) that considers all coding variants regardless of allele frequency. We describe an ab initio, genetics-first method making molecular knowledge-based interpretations for exome-wide non-synonymous variants for phenotypes at the organism and cellular level. By using this reverse approach, we identify plausible genetic causes for developmental disorders that have eluded other established methods and present molecular hypotheses for the causal genetics of 40 phenotypes generated from a direct-to-consumer genotype cohort. This system offers a chance to extract further discovery from genetic data after standard tools have been applied.

Recognition of sites of functional specialisation in all known eukaryotic protein kinase families.

The conserved function of protein phosphorylation, catalysed by members of protein kinase superfamily, is regulated in different ways in different kinase families. Further, differences in activating triggers, cellular localisation, domain architecture and substrate specificity between kinase families are also well known. While the transfer of γ-phosphate from ATP to the hydroxyl group of Ser/Thr/Tyr is mediated by a conserved Asp, the characteristic functional and regulatory sites are specialized at the level of families or sub-families. Such family-specific sites of functional specialization are unknown for most families of kinases. In this work, we systematically identify the family-specific residue features by comparing the extent of conservation of physicochemical properties, Shannon entropy and statistical probability of residue distributions between families of kinases. An integrated discriminatory score, which combines these three features, is developed to demarcate the functionally specialized sites in a kinase family from other sites. We achieved an area under ROC curve of 0.992 for the discrimination of kinase families. Our approach was extensively tested on well-studied families CDK and MAPK, wherein specific protein interaction sites and substrate recognition sites were successfully detected (p-value < 0.05). We also find that the known family-specific oncogenic driver mutation sites were scored high by our method. The method was applied to all known kinases encompassing 107 families from diverse eukaryotic organisms leading to a comprehensive list of family-specific functional sites. Apart from other uses, our method facilitates identification of specific protein interaction sites and drug target sites in a kinase family.

Interaction of the Saccharomyces cerevisiae RING-domain protein Nse1 with Nse3 and the Smc5/6 complex is required for chromosome replication and stability.

Genomic stability is maintained by the concerted actions of numerous protein complexes that participate in chromosomal duplication, repair, and segregation. The Smc5/6 complex is an essential multi-subunit complex crucial for repair of DNA double-strand breaks. Two of its subunits, Nse1 and Nse3, are homologous to the RING-MAGE complexes recently described in human cells. We investigated the contribution of the budding yeast Nse1 RING-domain by isolating a mutant nse1-103 bearing substitutions in conserved Zinc-coordinating residues of the RING-domain that is hypersensitive to genotoxic stress and temperature. The nse1-103 mutant protein was defective in interaction with Nse3 and other Smc5/6 complex subunits, Nse4 and Smc5. Chromosome loss was enhanced, accompanied by a delay in the completion of replication and a modest defect in sister chromatid cohesion, in nse1-103. The nse1-103 mutant was synthetic sick with rrm3∆ (defective in fork passage through pause sites), this defect was rescued by inactivation of Tof1, a subunit of the fork protection complex that enforces pausing. The temperature sensitivity of nse1-103 was partially suppressed by deletion of MPH1, encoding a DNA-helicase. Homology modeling of the structure of the budding yeast Nse1-Nse3 heterodimer based on the human Nse1-MAGEG1 structure suggests a similar organization and indicates that perturbation of the Zn-coordinating cluster has the potential to allosterically alter structural elements at the Nse1/Nse3 interaction interface that may abrogate their association. Our findings demonstrate that the budding yeast Nse1 RING-domain organization is important for interaction with Nse3, which is crucial for completion of chromosomal replication, cohesion, and maintenance of chromosome stability.

Sendai virus recruits cellular villin to remodel actin cytoskeleton during fusion with hepatocytes.

Reconstituted Sendai viral envelopes (virosomes) are well recognized for their promising potential in membrane fusion-mediated delivery of bioactive molecules to liver cells. Despite the known function of viral envelope glycoproteins in catalyzing fusion with cellular membrane, the role of host cell proteins remains elusive. Here, we used two-dimensional differential in-gel electrophoresis to analyze hepatic cells in early response to virosome-induced membrane fusion. Quantitative mass spectrometry together with biochemical analysis revealed that villin, an actin-modifying protein, is differentially up-regulated and phosphorylated at threonine 206-an early molecular event during membrane fusion. We found that villin influences actin dynamics and that this influence, in turn, promotes membrane mixing through active participation of Sendai viral envelope glycoproteins. Modulation of villin in host cells also resulted in a discernible effect on the entry and egress of progeny Sendai virus. Taken together, these results suggest a novel mechanism of regulated viral entry in animal cells mediated by host factor villin.

Conservation of structural fluctuations in homologous protein kinases and its implications on functional sites.

Our aim is to explore the similarities in structural fluctuations of homologous kinases. Gaussian Network Model based Normal Mode Analysis was performed on 73 active conformation structures in Ser/Thr/Tyr kinase superfamily. Categories of kinases with progressive evolutionary divergence, viz. (i) Same kinase with many crystal structures, (ii) Within-Subfamily, (iii) Within-Family, (iv) Within-Group, and (v) Across-Group, were analyzed. We identified a flexibility signature conserved in all kinases involving residues in and around the catalytic loop with consistent low-magnitude fluctuations. However, the overall structural fluctuation profiles are conserved better in closely related kinases (Within-Subfamily and Within-family) than in distant ones (Within-Group and Across-Group). A substantial 65.4% of variation in flexibility was not accounted by variation in sequences or structures. Interestingly, we identified substructural residue-wise fluctuation patterns characteristic of kinases of different categories. Specifically, we recognized statistically significant fluctuations unique to families of protein kinase A, cyclin-dependent kinases, and nonreceptor tyrosine kinases. These fluctuation signatures localized to sites known to participate in protein-protein interactions typical of these kinase families. We report for the first time that residues characterized by fluctuations unique to the group/family are involved in interactions specific to the group/family. As highlighted for Src family, local regions with differential fluctuations are proposed as attractive targets for drug design. Overall, our study underscores the importance of consideration of fluctuations, over and above sequence and structural features, in understanding the roles of sites characteristic of kinases. Proteins 2016; 84:957-978. © 2016 Wiley Periodicals, Inc.

A Gaussian network model study suggests that structural fluctuations are higher for inactive states than active states of protein kinases.

We performed Gaussian network model based normal mode analysis of 3-dimensional structures of multiple active and inactive forms of protein kinases. In 14 different kinases, a more number of residues (1095) show higher structural fluctuations in inactive states than those in active states (525), suggesting that, in general, mobility of inactive states is higher than active states. This statistically significant difference is consistent with higher crystallographic B-factors and conformational energies for inactive than active states, suggesting lower stability of inactive forms. Only a small number of inactive conformations with the DFG motif in the "in" state were found to have fluctuation magnitudes comparable to the active conformation. Therefore our study reports for the first time, intrinsic higher structural fluctuation for almost all inactive conformations compared to the active forms. Regions with higher fluctuations in the inactive states are often localized to the αC-helix, αG-helix and activation loop which are involved in the regulation and/or in structural transitions between active and inactive states. Further analysis of 476 kinase structures involved in interactions with another domain/protein showed that many of the regions with higher inactive-state fluctuation correspond to contact interfaces. We also performed extensive GNM analysis of (i) insulin receptor kinase bound to another protein and (ii) holo and apo forms of active and inactive conformations followed by multi-factor analysis of variance. We conclude that binding of small molecules or other domains/proteins reduce the extent of fluctuation irrespective of active or inactive forms. Finally, we show that the perceived fluctuations serve as a useful input to predict the functional state of a kinase.