Grappling with anisotropic data, pseudo-merohedral twinning and pseudo-translational noncrystallographic symmetry: a case study involving pyruvate kinase

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.

Characterizing the Effect of the Lysine Deacetylation Modification on Enzyme Activity of Pyruvate Kinase I and Pathogenicity of Vibrio alginolyticus

Pyruvate kinase I (PykF) is one of the key enzymes of glycolysis and plays a crucial role in bacterial metabolism. Several acetylation sites of Vibrio alginolyticus PykF were reported in previous studies and then 11 sites were first verified in this study, however, the specific roles of PykF acetylation remains unclear. Overlap-PCR and homologous recombination were implied to delete V. alginolyticus pykF gene and constructed complementary strains of site-directed mutagenesis for the further research focus on the deacetylation regulation on PykF. The results showed that the pyruvate kinase activity was sharply suppressed in the deacetylation status of K52, K68, and K317 of PykF, as well as the extracellular protease activity was significantly decreased in the deacetylation status of K52 and K68, but not induced with K317. Moreover, the growth rates of V. alginolyticus were not influenced with these three deacetylation sites. The ΔpykF mutant exhibited a 6-fold reduction in virulence to zebrafish. Site-directed mutations of K52R and K68R also showed reduced virulence while mutations of K317R didn't. The in vitro experiments showed that PykF was acetylated by acetyl phosphate (AcP), with the increase of incubation time by AcP, the acetylation level of PykF increased while the enzyme activity of PykF decreased correspondingly. Besides, PykF was deacetylated by CobB deacetylase and in result that the deacetylation was significantly down-regulated while the pyruvate kinase activity of PykF increased. Moreover, deletion of cobB gene had no significant difference in pyruvate kinase activity. These results confirm that CobB can regulate the acetylation level and pyruvate kinase activity of PykF. In summary, the results of this study provide a theoretical basis for further understanding of the deacetylation modification of PykF. It provides a new idea for the prevention and cure of vibriosis.

High-throughput characterization of mutations in genes that drive clonal evolution using multiplex adaptome capture sequencing

Understanding how cells are likely to evolve can guide medical interventions and bioengineering efforts that must contend with unwanted mutations. The adaptome of a cell-the neighborhood of genetic changes that are most likely to drive adaptation in a given environment-can be mapped by tracking rare beneficial variants during the early stages of clonal evolution. We used multiplex adaptome capture sequencing (mAdCap-seq), a procedure that combines unique molecular identifiers and hybridization-based enrichment, to characterize mutations in eight Escherichia coli genes known to be under selection in a laboratory environment. We tracked 301 mutations at frequencies as low as 0.01% and inferred the fitness effects of 240 of these mutations. There were distinct molecular signatures of selection on protein structure and function for the three genes with the most beneficial mutations. Our results demonstrate how mAdCap-seq can be used to deeply profile a targeted portion of a cell's adaptome.

The lid domain is important, but not essential, for catalysis of Escherichia coli pyruvate kinase

Pyruvate kinase catalyses the final step of the glycolytic pathway in central energy metabolism. The monomeric structure comprises three domains: a catalytic TIM-barrel, a regulatory domain involved in allosteric activation, and a lid domain that encloses the substrates. The lid domain is thought to close over the TIM-barrel domain forming contacts with the substrates to promote catalysis and may be involved in stabilising the activated state when the allosteric activator is bound. However, it remains unknown whether the lid domain is essential for pyruvate kinase catalytic or regulatory function. To address this, we removed the lid domain of Escherichia coli pyruvate kinase type 1 (PK^TIM+Reg) using protein engineering. Biochemical analyses demonstrate that, despite the absence of key catalytic residues in the lid domain, PK^TIM+Reg retains a low level of catalytic activity and has a reduced binding affinity for the substrate phosphoenolpyruvate. The enzyme retains allosteric activation, but the regulatory profile of the enzyme is changed relative to the wild-type enzyme. Analytical ultracentrifugation and small-angle X-ray scattering data show that, beyond the loss of the lid domain, the PK^TIM+Reg structure is not significantly altered and is consistent with the wild-type tetramer that is assembled through interactions at the TIM and regulatory domains. Our results highlight the contribution of the lid domain for facilitating pyruvate kinase catalysis and regulation, which could aid in the development of small molecule inhibitors for pyruvate kinase and related lid-regulated enzymes.

Comparing the Conformational Stability of Pyruvate Kinase in the Gas Phase and in Solution

Collision induced unfolding (CIU) is increasingly used to characterize protein complexes in the gas phase and is often employed to detect ligand binding-induced conformational stabilization. However, the extent to which gas-phase conformational stabilities measured by CIU reflect analogous parameters in solution is not yet clear, particularly for systems where conformational and protein complex stability are modulated by point mutation. Here, we compare CIU-derived relative stabilities of four point mutants of the homotetramer pyruvate kinase to solution stabilities measured by differential scanning fluorimetry (DSF) and solution conformational dynamics measured by time-resolved electrospray ionization hydrogen-deuterium exchange (TRESI-HDX). Our results demonstrate that both destabilization of the tetrameric state and generally reduced conformational stability of the monomer in solution are well correlated to lower onset energies for specific unfolding transitions observed in CIU. However, this correlation not fully retained when comparing CIU to HDX data, where the latter measurement is strongly impacted by conformational dynamics within the tetramer.

Optimizing the refinement of merohedrally twinned P61 HIV-1 protease-inhibitor cocrystal structures

Twinning is a crystal-growth anomaly in which protein monomers exist in different orientations but are related in a specific way, causing diffraction reflections to overlap. Twinning imposes additional symmetry on the data, often leading to the assignment of a higher symmetry space group. Specifically, in merohedral twinning, reflections from each monomer overlap and require a twin law to model unique structural data from overlapping reflections. Neglecting twinning in the crystallographic analysis of quasi-rotationally symmetric homo-oligomeric protein structures can mask the degree of structural non-identity between monomers. In particular, any deviations from perfect symmetry will be lost if higher than appropriate symmetry is applied during crystallographic analysis. Such cases warrant choosing between the highest symmetry space group possible or determining whether the monomers have distinguishable structural asymmetries and thus require a lower symmetry space group and a twin law. Using hexagonal cocrystals of HIV-1 protease, a C2-symmetric homodimer whose symmetry is broken by bound ligand, it is shown that both assigning a lower symmetry space group and applying a twin law during refinement are critical to achieving a structural model that more accurately fits the electron density. By re-analyzing three recently published HIV-1 protease structures, improvements in nearly every crystallographic metric are demonstrated. Most importantly, a procedure is demonstrated where the inhibitor can be reliably modeled in a single orientation. This protocol may be applicable to many other homo-oligomers in the PDB.

Effects of Beneficial Mutations in pykF Gene Vary over Time and across Replicate Populations in a Long-Term Experiment with Bacteria

The fitness effects of mutations can depend on the genetic backgrounds in which they occur and thereby influence future opportunities for evolving populations. In particular, mutations that fix in a population might change the selective benefit of subsequent mutations, giving rise to historical contingency. We examine these effects by focusing on mutations in a key metabolic gene, pykF, that arose independently early in the history of 12 Escherichia coli populations during a long-term evolution experiment. Eight different evolved nonsynonymous mutations conferred similar fitness benefits of ∼10% when transferred into the ancestor, and these benefits were greater than the one conferred by a deletion mutation. In contrast, the same mutations had highly variable fitness effects, ranging from ∼0% to 25%, in evolved clones isolated from the populations at 20,000 generations. Two mutations that were moved into these evolved clones conferred similar fitness effects in a given clone, but different effects between the clones, indicating epistatic interactions between the evolved pykF alleles and the other mutations that had accumulated in each evolved clone. We also measured the fitness effects of six evolved pykF alleles in the same populations in which they had fixed, but at seven time points between 0 and 50,000 generations. Variation in fitness effects was high at intermediate time points, and declined to a low level at 50,000 generations, when the mean fitness effect was lowest. Our results demonstrate the importance of genetic context in determining the fitness effects of different beneficial mutations even within the same gene.

Conformational Dynamics and Allostery in Pyruvate Kinase

Pyruvate kinase catalyzes the final step in glycolysis and is allosterically regulated to control flux through the pathway. Two models are proposed to explain how Escherichia coli pyruvate kinase type 1 is allosterically regulated: the "domain rotation model" suggests that both the domains within the monomer and the monomers within the tetramer reorient with respect to one another; the "rigid body reorientation model" proposes only a reorientation of the monomers within the tetramer causing rigidification of the active site. To test these hypotheses and elucidate the conformational and dynamic changes that drive allostery, we performed time-resolved electrospray ionization mass spectrometry coupled to hydrogen-deuterium exchange studies followed by mutagenic analysis to test the activation mechanism. Global exchange experiments, supported by thermostability studies, demonstrate that fructose 1,6-bisphosphate binding to the allosteric domain causes a shift toward a globally more dynamic ensemble of conformations. Mapping deuterium exchange to peptides within the enzyme highlight site-specific regions with altered conformational dynamics, many of which increase in conformational flexibility. Based upon these and mutagenic studies, we propose an allosteric mechanism whereby the binding of fructose 1,6-bisphosphate destabilizes an α-helix that bridges the allosteric and active site domains within the monomeric unit. This destabilizes the β-strands within the (β/α)8-barrel domain and the linked active site loops that are responsible for substrate binding. Our data are consistent with the domain rotation model but inconsistent with the rigid body reorientation model given the increased flexibility at the interdomain interface, and we can for the first time explain how fructose 1,6-bisphosphate affects the active site.