PYK-SubstitutionOME: an integrated database containing allosteric coupling, ligand affinity and mutational, structural, pathological, bioinformatic and computational information about pyruvate kinase isozymes.

Interpreting changes in patient genomes, understanding how viruses evolve and engineering novel protein function all depend on accurately predicting the functional outcomes that arise from amino acid substitutions. To that end, the development of first-generation prediction algorithms was guided by historic experimental datasets. However, these datasets were heavily biased toward substitutions at positions that have not changed much throughout evolution (i.e. conserved). Although newer datasets include substitutions at positions that span a range of evolutionary conservation scores, these data are largely derived from assays that agglomerate multiple aspects of function. To facilitate predictions from the foundational chemical properties of proteins, large substitution databases with biochemical characterizations of function are needed. We report here a database derived from mutational, biochemical, bioinformatic, structural, pathological and computational studies of a highly studied protein family-pyruvate kinase (PYK). A centerpiece of this database is the biochemical characterization-including quantitative evaluation of allosteric regulation-of the changes that accompany substitutions at positions that sample the full conservation range observed in the PYK family. We have used these data to facilitate critical advances in the foundational studies of allosteric regulation and protein evolution and as rigorous benchmarks for testing protein predictions. We trust that the collected dataset will be useful for the broader scientific community in the further development of prediction algorithms. Database URL https://github.com/djparente/PYK-DB.

Determinants of cytoplasmic microtubule depolymerization during ciliogenesis in Chlamydomonas reinhardtii.

At the core of cilia are microtubules which are important for establishing length and assisting ciliary assembly and disassembly; however, another role for microtubule regulation on ciliogenesis lies outside of the cilium. The microtubule cytoskeleton is a highly dynamic structure which polymerizes and depolymerizes rapidly to assist in cellular processes. These processes have been studied across various organisms with chemical as well as genetic perturbations. However, these have generated conflicting data in terms of the role of cytoplasmic microtubules (CytoMTs) and free tubulin dynamics during ciliogenesis. Here we look at the relationship between ciliogenesis and cytoplasmic microtubule dynamics in Chlamydomonas reinhardtii using chemical and mechanical perturbations. We find that not only can stabilized CytoMTs allow for normal ciliary assembly, but high calcium concentrations and low pH-induced deciliation cause CytoMTs to depolymerize separately from ciliary shedding. In addition, we find that ciliary shedding through mechanical shearing, cilia regenerate earlier despite intact CytoMTs. Our data suggests that cytoplasmic microtubules are not a sink for a limiting pool of cytoplasmic tubulin in Chlamydomonas, depolymerization that occurs following deciliation is a consequence rather than a requirement for ciliogenesis, and intact CytoMTs in the cytoplasm and the proximal cilium support more efficient ciliary assembly.

The ERK activator, BCI, inhibits ciliogenesis and causes defects in motor behavior, ciliary gating, and cytoskeletal rearrangement.

MAPK pathways are well-known regulators of the cell cycle, but they have also been found to control ciliary length in a wide variety of organisms and cell types from Caenorhabditis elegans neurons to mammalian photoreceptors through unknown mechanisms. ERK1/2 is a MAP kinase in human cells that is predominantly phosphorylated by MEK1/2 and dephosphorylated by the phosphatase DUSP6. We have found that the ERK1/2 activator/DUSP6 inhibitor, (E)-2-benzylidene-3-(cyclohexylamino)-2,3-dihydro-1H-inden-1-one (BCI), inhibits ciliary maintenance in Chlamydomonas and hTERT-RPE1 cells and assembly in Chlamydomonas These effects involve inhibition of total protein synthesis, microtubule organization, membrane trafficking, and KAP-GFP motor dynamics. Our data provide evidence for various avenues for BCI-induced ciliary shortening and impaired ciliogenesis that gives mechanistic insight into how MAP kinases can regulate ciliary length.

Lithium-induced ciliary lengthening sparks Arp2/3 complex-dependent endocytosis.

Ciliary length is highly regulated, but can be disrupted by lithium, which causes ciliary elongation across cell types and organisms. We used the algal system Chlamydomonas reinhardtii to investigate the mechanism behind lithium-induced ciliary elongation. Protein synthesis is not required for lengthening, and the target of lithium, GSK3, has substrates that can influence membrane dynamics. Further, ciliary assembly requires a supply of ciliary membrane as well as protein. Lithium-treated cilia elongate normally with brefeldin treatment, but dynasore treatment produced defective lengthening suggesting a source of membrane from the cell surface rather than the Golgi. Genetic or chemical perturbation of the Arp2/3 complex or dynamin, required for endocytosis, blocks lithium-induced ciliary lengthening. Finally, we found an increase in Arp2/3 complex- and endocytosis-dependent actin filaments near the ciliary base upon lithium treatment. Our results identify a mechanism for lithium-mediated cilium lengthening and demonstrate the endocytic pathway for cilium membrane supply in algae is likely a conserved mechanism given lithium's conserved effects across organisms.

Rheostat functional outcomes occur when substitutions are introduced at nonconserved positions that diverge with speciation.

When amino acids vary during evolution, the outcome can be functionally neutral or biologically-important. We previously found that substituting a subset of nonconserved positions, "rheostat" positions, can have surprising effects on protein function. Since changes at rheostat positions can facilitate functional evolution or cause disease, more examples are needed to understand their unique biophysical characteristics. Here, we explored whether "phylogenetic" patterns of change in multiple sequence alignments (such as positions with subfamily specific conservation) predict the locations of functional rheostat positions. To that end, we experimentally tested eight phylogenetic positions in human liver pyruvate kinase (hLPYK), using 10-15 substitutions per position and biochemical assays that yielded five functional parameters. Five positions were strongly rheostatic and three were non-neutral. To test the corollary that positions with low phylogenetic scores were not rheostat positions, we combined these phylogenetic positions with previously-identified hLPYK rheostat, "toggle" (most substitution abolished function), and "neutral" (all substitutions were like wild-type) positions. Despite representing 428 variants, this set of 33 positions was poorly statistically powered. Thus, we turned to the in vivo phenotypic dataset for E. coli lactose repressor protein (LacI), which comprised 12-13 substitutions at 329 positions and could be used to identify rheostat, toggle, and neutral positions. Combined hLPYK and LacI results show that positions with strong phylogenetic patterns of change are more likely to exhibit rheostat substitution outcomes than neutral or toggle outcomes. Furthermore, phylogenetic patterns were more successful at identifying rheostat positions than were co-evolutionary or eigenvector centrality measures of evolutionary change.

Separable roles for RanGTP in nuclear and ciliary trafficking of a kinesin-2 subunit.

Kinesin is part of the microtubule-binding motor protein superfamily, which serves important roles in cell division and intraorganellar transport. The heterotrimeric kinesin-2, consisting of the heterodimeric motor subunits, kinesin family member 3A/3B (KIF3A/3B), and kinesin-associated protein 3 (KAP3), is highly conserved across species from the unicellular eukaryote Chlamydomonas to humans. It plays diverse roles in cargo transport including anterograde (base to tip) trafficking in cilia. However, the molecular determinants mediating trafficking of heterotrimeric kinesin-2 itself are poorly understood. It has been previously suggested that ciliary transport is analogous to nuclear transport mechanisms. Using Chlamydomonas and human telomerase reverse transcriptase-retinal pigment epithelial cell line, we show that RanGTP, a small GTPase that dictates nuclear transport, regulates ciliary trafficking of KAP3, a key component for functional kinesin-2. We found that the armadillo-repeat region 6 to 9 (ARM6-9) of KAP3, required for its nuclear translocation, is also necessary and sufficient for its targeting to the ciliary base. Given that KAP3 is essential for cilium formation and there are the emerging roles for RanGTP/importin ß in ciliary protein targeting, we further investigated the effect of RanGTP in cilium formation and maintenance. We found that precise control of RanGTP levels, revealed by different Ran mutants, is crucial for cilium formation and maintenance. Most importantly, we were able to provide orthogonal support in an algal model system that segregates RanGTP regulation of ciliary protein trafficking from its nuclear roles. Our work provides important support for the model that nuclear import mechanisms have been co-opted for independent roles in ciliary import.

Identification of biochemically neutral positions in liver pyruvate kinase.

Understanding how each residue position contributes to protein function has been a long-standing goal in protein science. Substitution studies have historically focused on conserved protein positions. However, substitutions of nonconserved positions can also modify function. Indeed, we recently identified nonconserved positions that have large substitution effects in human liver pyruvate kinase (hLPYK), including altered allosteric coupling. To facilitate a comparison of which characteristics determine when a nonconserved position does vs does not contribute to function, the goal of the current work was to identify neutral positions in hLPYK. However, existing hLPYK data showed that three features commonly associated with neutral positions-high sequence entropy, high surface exposure, and alanine scanning-lacked the sensitivity needed to guide experimental studies. We used multiple evolutionary patterns identified in a sequence alignment of the PYK family to identify which positions were least patterned, reasoning that these were most likely to be neutral. Nine positions were tested with a total of 117 amino acid substitutions. Although exploring all potential functions is not feasible for any protein, five parameters associated with substrate/effector affinities and allosteric coupling were measured for hLPYK variants. For each position, the aggregate functional outcomes of all variants were used to quantify a "neutrality" score. Three positions showed perfect neutral scores for all five parameters. Furthermore, the nine positions showed larger neutral scores than 17 positions located near allosteric binding sites. Thus, our strategy successfully enriched the dataset for positions with neutral and modest substitutions.

RheoScale: A tool to aggregate and quantify experimentally determined substitution outcomes for multiple variants at individual protein positions.

Human mutations often cause amino acid changes (variants) that can alter protein function or stability. Some variants fall at protein positions that experimentally exhibit "rheostatic" mutation outcomes (different amino acid substitutions lead to a range of functional outcomes). In ongoing studies of rheostat positions, we encountered the need to aggregate experimental results from multiple variants, to describe the overall roles of individual positions. Here, we present "RheoScale" which generates quantitative scores to discriminate rheostat positions from those with "toggle" (most substitutions abolish function) or "neutral" (most substitutions have wild-type function) outcomes. RheoScale scores facilitate correlations of experimental data (such as binding affinity or stability) with structural and bioinformatic analyses. The RheoScale calculator is encoded into a Microsoft Excel workbook and an R script. Example analyses are shown for three model protein systems, including one assessed via deep mutational scanning. The RheoScale calculator quickly and efficiently provided quantitative descriptions that were in good agreement with prior qualitative observations. As an example application, scores were compared to the example proteins' structures; strong rheostat positions tended to occur in dynamic locations. In the future, RheoScale scores can be easily integrated into computational studies to facilitate improved algorithms for predicting outcomes of human variants.