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

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.

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

The rise of the computer as a powerful tool for model building and refinement has revolutionized the field of structure determination for large biomolecular systems. Despite the wide availability of robust experimental methods capable of resolving structural details across a range of spatiotemporal resolutions, computational hybrid methods have the unique ability to integrate the diverse data from multimodal techniques such as X-ray crystallography and electron microscopy into consistent, fully atomistic structures. Here, commonly employed strategies for computational real-space structural refinement are reviewed, and their specific applications are illustrated for several large macromolecular complexes: ribosome, virus capsids, chemosensory array, and photosynthetic chromatophore. The increasingly important role of computational methods in large-scale structural refinement, along with current and future challenges, is discussed.

Virus capsids are protein shells that package the viral genome. Although their morphology and biological functions can vary markedly, capsids often play critical roles in regulating viral infection pathways. A detailed knowledge of virus capsids, including their dynamic structure, interactions with cellular factors, and the specific roles that they play in the replication cycle, is imperative for the development of antiviral therapeutics. The following Perspective introduces an emerging area of computational biology that focuses on the dynamics of virus capsids and capsid–protein assemblies, with particular emphasis on the effects of small-molecule drug binding on capsid structure, stability, and allosteric pathways. When performed at chemical detail, molecular dynamics simulations can reveal subtle changes in virus capsids induced by drug molecules a fraction of their size. Here, the current challenges of performing all-atom capsid–drug simulations are discussed, along with an outlook on the applicability of virus capsid simulations to reveal novel drug targets.

The cellular process responsible for providing energy for most life on Earth, namely photosynthetic light-harvesting, requires the cooperation of hundreds of proteins across an organelle, involving length and time scales spanning several orders of magnitude over quantum and classical regimes. Simulation and visualization of this fundamental energy conversion process pose many unique methodological and computational challenges. We present, in two accompanying movies, light-harvesting in the photosynthetic apparatus found in purple bacteria, the so-called chromatophore. The movies are the culmination of three decades of modeling efforts, featuring the collaboration of theoretical, experimental, and computational scientists. We describe the techniques that were used to build, simulate, analyze, and visualize the structures shown in the movies, and we highlight cases where scientific needs spurred the development of new parallel algorithms that efficiently harness GPU accelerators and petascale computers.

Defining the molecular interaction between Gag proteins in an assembled hexagonal lattice of immature retrovirus particles is crucial for elucidating the mechanisms of virus assembly and maturation. Recent advances in cryo-electron microscopy have yielded subnanometer structural information on the morphology of immature Gag lattices, making computational modeling and simulations feasible for investigating the Gag-Gag interactions at the atomic level. We have examined the structure of Rous sarcoma virus (RSV) using all-atom molecular dynamics simulations and in vitro assembly, to create the first all-atom model of an immature retroviral lattice. Microseconds-long replica exchange molecular dynamics simulation of the spacer peptide (SP)-nucleocapsid (NC) subdomains results in a 6-helix bundle with amphipathic properties. The resulting model of the RSV Gag lattice shows features and dynamics of the capsid protein (CA) with implications for the maturation process, and confirms the stabilizing role of the upstream and downstream regions of Gag, namely p10 and SP-NC.