The millions of protein sequences generated by genomics are expected to transform protein engineering and personalized medicine. To achieve these goals, tools for predicting outcomes of amino acid changes must be improved. Currently, advances are hampered by insufficient experimental data about nonconserved amino acid positions. Since the property “nonconserved” is identified using a sequence alignment, we designed experiments to recapitulate that context: Mutagenesis and functional characterization was carried out in 15 LacI/GalR homologs (rows) at 12 nonconserved positions (columns). Multiple substitutions were made at each position, to reveal how various amino acids of a nonconserved column were tolerated in each protein row. Results showed that amino acid preferences of nonconserved positions were highly context-dependent, had few correlations with physico-chemical similarities, and were not predictable from their occurrence in natural LacI/GalR sequences. Further, unlike the “toggle switch” behaviors of conserved positions, substitutions at nonconserved positions could be rank-ordered to show a “rheostatic”, progressive effect on function that spanned several orders of magnitude. Comparisons to various sequence analyses suggested that conserved and strongly co-evolving positions act as functional toggles, whereas other important, nonconserved positions serve as rheostats for modifying protein function. Both the presence of rheostat positions and the sequence analysis strategy appear to be generalizable to other protein families and should be considered when engineering protein modifications or predicting the impact of protein polymorphisms.

Concomitant with the genomic era, many bioinformatics programs have been developed to identify functionally important positions from sequence alignments of protein families. To evaluate these analyses, many have used the LacI/GalR family and determined whether positions predicted to be "important" are validated by published experiments. However, we previously noted that predictions do not identify all of the experimentally important positions present in the linker regions of these homologs. In an attempt to reconcile these differences, we corrected and expanded the LacI/GalR sequence set commonly-used in sequence/function analyses. Next, a variety of analyses were carried out (1) for the entire LacI/GalR sequence set and (2) for a subset of homologs with functionally-important "YxPxxxAxxL" motifs in their linkers. This strategy was devised to determine whether predictions could be improved by knowledge-based sequence sorting and -for some analyses -did increase the number of linker positions identified. However, two functionally important linker positions were not reliably identified by any analysis. Finally, we compared the new predictions to all known experimental data for E. coli LacI and three homologous linkers. From these, we estimate that >50% of positions are important to the functions of the LacI/GalR homologs. In corollary, neutral positions might occur less frequently and might be easier to detect in sequence analyses. Although analyses have successfully guided mutations that partially exchange protein functions, a better experimental understanding of the sequence/ function relationships in protein families would be helpful for uncovering the remaining rules used by nature to evolve new protein functions.

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.

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.

Protein families might evolve paralogous functions on their common tertiary scaffold in two ways. First, the locations of functionally-important sites might be “hard-wired” into the structure, with novel functions evolved by altering the amino acid (e.g. Ala vs Ser) at these positions. Alternatively, the tertiary scaffold might be adaptable, accommodating a unique set of functionally important sites for each paralogous function. To discriminate between these possibilities, we compared the set of functionally important sites in the six largest paralogous subfamilies of the LacI/GalR transcription repressor family. LacI/GalR paralogs share a common tertiary structure, but have low sequence identity (≤30%), and regulate a variety of metabolic processes. Functionally important positions were identified by conservation and co-evolutionary sequence analyses. Results showed that conserved positions use a mixture of the “hard-wired” and “accommodating” scaffold frameworks, but that the co-evolution networks were highly dissimilar between any pair of subfamilies. Therefore, the tertiary structure can accommodate multiple networks of functionally important positions. This possibility should be included when designing and interpreting sequence analyses of other protein families. Software implementing conservation and co-evolution analyses is available at https://sourceforge.net/projects/coevolutils/.

Neuroligins are post-synaptic, cellular adhesion molecules implicated in synaptic formation and function. NLGN2 is strongly linked to inhibitory, GABAergic signaling and is crucial for maintaining the excitation-inhibition balance in the brain. Disruption of the excitation-inhibition balance is associated with neuropsychiatric disease. In animal models, altered NLGN2 expression causes anxiety, developmental delay, motor discoordination, social impairment, aggression, and sensory processing defects. In humans, mutations in NLGN3 and NLGN4 are linked to autism and schizophrenia; NLGN2 missense variants are implicated in schizophrenia. Copy number variants encompassing NLGN2 on 17p13.1 are associated with autism, intellectual disability, metabolic syndrome, diabetes, and dysmorphic features, but an isolated NLGN2 nonsense variant has not yet been described in humans. Here, we describe a 15-year-old male with severe anxiety, obsessive-compulsive behaviors, developmental delay, autism, obesity, macrocephaly, and some dysmorphic features. Exome sequencing identified a heterozygous, de novo, c.441C>A p.(Tyr147Ter) variant in NLGN2 that is predicted to cause loss of normal protein function. This is the first report of an NLGN2 nonsense variant in humans, adding to the accumulating evidence that links synaptic proteins with a spectrum of neurodevelopmental phenotypes. © 2016 Wiley Periodicals, Inc.

As proteins evolve, amino acid positions key to protein structure or function are subject to mutational constraints. These positions can be detected by analyzing sequence families for amino acid conservation or for co-evolution between pairs of positions. Co-evolutionary scores are usually rank-ordered and thresholded to reveal the top pairwise scores, but they also can be treated as weighted networks. Here, we used network analyses to bypass a major complication of co-evolution studies: For a given sequence alignment, alternative algorithms usually identify different, top pairwise scores. We reconciled results from five commonly-used, mathematically divergent algorithms (ELSC, McBASC, OMES, SCA, and ZNMI), using the LacI/GalR and 1,6-bisphosphate aldolase protein families as models. Calculations used unthresholded co-evolution scores from which column-specific properties such as sequence entropy and random noise were subtracted; “central” positions were identified by calculating various network centrality scores. When compared among algorithms, network centrality methods, particularly eigenvector centrality, showed markedly better agreement than comparisons of the top pairwise scores. Positions with large centrality scores occurred at key structural locations and/or were functionally sensitive to mutations. Further, the top central positions often differed from those with top pairwise co-evolution scores: Instead of a few strong scores, central positions often had multiple, moderate scores. We conclude that eigenvector centrality calculations reveal a robust evolutionary pattern of constraints – detectable by divergent algorithms – that occur at key protein locations. Finally, we discuss the fact that multiple patterns co-exist in evolutionary data that, together, give rise to emergent protein functions.

In the goal of interpreting human exomes, when predictive programs assign functional importance to some positions, they implicitly assume the existence of non-important positions -those that accommodate many side chain chemistries without altering function ("neutral"). However, very few (if any) experimental studies have demonstrated the existence of neutral positions. We sought experimental evidence for neutral positions using human liver pyruvate kinase (hLPYK) as a model system. To that end, we used multiple evolutionary criteria to identify 20 possibly neutral positions. Nine positions were further tested with a total of 117 amino acid substitutions. Although "all" potential hLPYK functions can never be explored, we measured effects on 5 parameters associated with substrate/ligand affinities and allosteric coupling. At each position, the aggregate outcomes of multiple variants were used to quantify "neutrality" scores. Three of the nine positions showed perfect neutral scores in all 5 parameters; a fourth position had high neutral scores. Although our strategy for predicting positions had low predictive power for the identification of neutral positions, all positions had neutral scores that were much higher than positions in functional sites. Given this evidence for the existence of neutral positions, similar studies should be carried out for other proteins to generate a database of well characterized neutral positions that can then be available to benchmark and validate predictions about amino acid substitutions.