Connecting dynamics to structural data from diverse experimental sources, molecular dynamics simulations permit the exploration of biological phenomena in unparalleled detail. Advances in simulations are moving the atomic resolution descriptions of biological systems into the million-to-billion atom regime, in which numerous cell functions reside. In this opinion, we review the progress, driven by large-scale molecular dynamics simulations, in the study of viruses, ribosomes, bioenergetic systems, and other diverse applications. These examples highlight the utility of molecular dynamics simulations in the critical task of relating atomic detail to the function of supramolecular complexes, a task that cannot be achieved by smaller-scale simulations or existing experimental approaches alone.

Asparaginyl endopeptidases (AEPs) are cysteine proteases which break Asx (Asn/Asp)-Xaa bonds in acidic conditions. Despite sharing a conserved overall structure with AEPs, certain plant enzymes such as butelase 1 act as a peptide asparaginyl ligase (PAL) and catalyze Asx-Xaa bond formation in near-neutral conditions. PALs also serve as macrocyclases in the biosynthesis of cyclic peptides. Here, we address the question of how a PAL can function as a ligase rather than a protease. Based on sequence homology of butelase 1, we identified AEPs and PALs from the cyclic peptideproducing plants Viola yedoensis (Vy) and Viola canadensis (Vc) of the Violaceae family. Using a crystal structure of a PAL obtained at 2.4-Å resolution coupled to mutagenesis studies, we discovered ligase-activity determinants flanking the S1 site, namely LAD1 and LAD2 located around the S2 and S1′ sites, respectively, which modulate ligase activity by controlling the accessibility of water or amine nucleophile to the S-ester intermediate. Recombinantly expressed VyPAL1-3, predicted to be PALs, were confirmed to be ligases by functional studies. In addition, mutagenesis studies on VyPAL1-3, VyAEP1, and VcAEP supported our prediction that LAD1 and LAD2 are important for ligase activity. In particular, mutagenesis targeting LAD2 selectively enhanced the ligase activity of VyPAL3 and converted the protease VcAEP into a ligase. The definition of structural determinants required for ligation activity of the asparaginyl ligases presented here will facilitate genomic identification of PALs and engineering of AEPs into PALs. peptide ligase | data mining | ligase-activity determinant

Butelase-1 is a peptide asparaginyl ligase (PAL) that efficiently catalyzes peptide bond formation after Asn/Asp (Asx), whereas its homolog butelase-2 is an asparaginyl endopeptidase (AEP) that catalyzes the reverse reaction, hydrolyzing Asx-peptide bonds. Since PALs and AEPs share essentially similar overall structures, we surmised that the S2 and S1' substrate-binding pockets immediately flanking the catalytic S1 site constitute the major ligase activity determinants (LADs) that control the catalytic directionality of these enzymes. Here, we report the successful conversion of butelase-2 into butelase-1-like ligases based on the LAD hypothesis. We prepared 23 LAD-mutants by mutating residues of the S2 pocket (LAD1) and/or the S1' pocket (LAD2) of butelase-2 into homologous residues in PALs. These LAD-mutants markedly diminished protease activity and increased ligase activity. In contrast, substituting 12 non-LAD residues to the corresponding residues in butelase-1 did not change their protease profiles. At physiological pH, mutations targeting both LADs resulted in shifting the catalytic directionality from 95% hydrolysis to >95% peptide ligation. This results in the engineering of efficient recombinant peptide ligases that were demonstrated to be useful for macrocyclization

Pyruvate kinase (PYK) is an essential glycolytic enzyme that controls glycolytic flux and is critical for ATP production in all organisms, with tight regulation by multiple metabolites. Yet the allosteric mechanisms governing PYK activity in bacterial pathogens are poorly understood. Here we report biochemical, structural and metabolomic evidence that Mycobacterium tuberculosis (Mtb) PYK uses AMP and glucose-6-phosphate (G6P) as synergistic allosteric activators that function as a molecular “OR logic gate” to tightly regulate energy and glucose metabolism. G6P was found to bind to a previously unknown site adjacent to the canonical site for AMP. Kinetic data and structural network analysis further show that AMP and G6P work synergistically as allosteric activators. Importantly, metabolome profiling in the Mtb surrogate, Mycobacterium bovis BCG, reveals significant changes in AMP and G6P levels during nutrient deprivation, which provides insights into how a PYK OR gate would function during the stress of Mtb infection.