Conformational Entropy in PEPCK Catalysis: Dynamic Motions Critical to Function

Enzymes increase the rate at which chemical reactions occur. How they achieve this rate increase has been an area of intense research for many years. Multiple studies have shown that dynamic motions of individual protein segments can play a role in enzyme function. Our enzyme of interest, phosphoenolpyruvate carboxykinase (PEPCK), has multiple dynamic motions that work in concert during the catalytic cycle. One of these mobile elements, termed the omega-loop due to its architecture, acts like a lid that closes over the active site and has been hypothesized to serve multiple functions during the catalytic cycle. Furthermore, the act of opening and closing the lid domain is thought to be a delicate free energy balance, meaning the lid pays an energetic penalty to close, which is offset by the energetic contribution from substrate binding. To investigate the role of the loop, specifically the open/closed transition, and to probe for unappreciated roles, multiple mutant enzymes were created. In the first set of mutations the lid region was removed to investigate other potential roles for the omega-loop. In the second set two different amino acid residues (Ala-467 and Glu-89) were mutated to either increase the energetic penalty of lid closure or decrease the energetic contribution from substrate binding. The resultant mutants were characterized via structure/function experiments. The data revealed two new roles for the omega-loop lid, confirmed the delicate free energy balance, and shed light on the energetic pathway between the active site and the lid. Taking into account the data from the WT and all seven mutant PEPCK enzymes, a revised model for the role of the omega-loop lid during catalysis is proposed: The role of the lid during the catalytic cycle is to correctly position the substrates in the active site, stabilize the R-loop and N/C-terminal lobe closure, and sequester and protect the reaction intermediate

Molecular principles underlying dual RNA specificity in the Drosophila SNF protein

The first RNA recognition motif of the Drosophila SNF protein is an example of an RNA binding protein with multi-specificity. It binds different RNA hairpin loops in spliceosomal U1 or U2 small nuclear RNAs, and only in the latter case requires the auxiliary U2A′ protein. Here we investigate its functions by crystal structures of SNF alone and bound to U1 stem-loop II, U2A′ or U2 stem-loop IV and U2A′, SNF dynamics from NMR spectroscopy, and structure-guided mutagenesis in binding studies. We find that different loop-closing base pairs and a nucleotide exchange at the tips of the loops contribute to differential SNF affinity for the RNAs. U2A′ immobilizes SNF and RNA residues to restore U2 stem-loop IV binding affinity, while U1 stem-loop II binding does not require such adjustments. Our findings show how U2A′ can modulate RNA specificity of SNF without changing SNF conformation or relying on direct RNA contacts

Computational Investigations on Polymerase Actions in Gene Transcription and Replication Combining Physical Modeling and Atomistic Simulations

Polymerases are protein enzymes that move along nucleic acid chains and catalyze template-based polymerization reactions during gene transcription and replication. The polymerases also substantially improve transcription or replication fidelity through the non-equilibrium enzymatic cycles. We briefly review computational efforts that have been made toward understanding mechano-chemical coupling and fidelity control mechanisms of the polymerase elongation. The polymerases are regarded as molecular information motors during the elongation process. It requires a full spectrum of computational approaches from multiple time and length scales to understand the full polymerase functional cycle. We keep away from quantum mechanics based approaches to the polymerase catalysis due to abundant former surveys, while address only statistical physics modeling approach and all-atom molecular dynamics simulation approach. We organize this review around our own modeling and simulation practices on a single-subunit T7 RNA polymerase, and summarize commensurate studies on structurally similar DNA polymerases. For multi-subunit RNA polymerases that have been intensively studied in recent years, we leave detailed discussions on the simulation achievements to other computational chemical surveys, while only introduce very recently published representative studies, including our own preliminary work on structure-based modeling on yeast RNA polymerase II. In the end, we quickly go through kinetic modeling on elongation pauses and backtracking activities. We emphasize the fluctuation and control mechanisms of the polymerase actions, highlight the non-equilibrium physical nature of the system, and try to bring some perspectives toward understanding replication and transcription regulation from single molecular details to a genome-wide scale

Substrate Induced Structural and Dynamics Changes in Human Phosphomevalonate Iinase and Implications for Mechanism

Phosphomevalonate kinase (PMK) catalyzes an essential step in the mevalonate pathway, which is the only pathway for synthesis of isoprenoids and steroids in humans. PMK catalyzes transfer of the γ-phosphate of ATP to mevalonate 5-phosphate (M5P) to form mevalonate 5-diphosphate. Bringing these phosphate groups in proximity to react is especially challenging, given the high negative charge density on the four phosphate groups in the active site. As such, conformational and dynamics changes needed to form the Michaelis complex are of mechanistic interest. Herein, we report the characterization of substrate induced changes (Mg-ADP, M5P, and the ternary complex) in PMK using NMR-based dynamics and chemical shift perturbation measurements. Mg-ADP and M5P Kd\u27s were 6–60 μM in all complexes, consistent with there being little binding synergy. Binding of M5P causes the PMK structure to compress (τc = 13.5 nsec), whereas subsequent binding of Mg-ADP opens the structure up (τc = 15.6 nsec). The overall complex seems to stay very rigid on the psec-nsec timescale with an average NMR order parameter of S2 ∼0.88. Data are consistent with addition of M5P causing movement around a hinge region to permit domain closure, which would bring the M5P domain close to ATP to permit catalysis. Dynamics data identify potential hinge residues as H55 and R93, based on their low order parameters and their location in extended regions that connect the M5P and ATP domains in the PMK homology model. Likewise, D163 may be a hinge residue for the lid region that is homologous to the adenylate kinase lid, covering the “Walker-A” catalytic loop. Binding of ATP or ADP appears to cause similar conformational changes; however, these observations do not indicate an obvious role for γ-phosphate binding interactions. Indeed, the role of γ-phosphate interactions may be more subtle than suggested by ATP/ADP comparisons, because the conservative O to NH substitution in the β-γ bridge of ATP causes a dramatic decrease in affinity and induces few chemical shift perturbations. In terms of positioning of catalytic residues, binding of M5P induces a rigidification of Gly21 (adjacent to the catalytically important Lys22), although exchange broadening in the ternary complex suggests some motion on a slower timescale does still occur. Finally, the first nine residues of the N-terminus are highly disordered, suggesting that they may be part of a cleavable signal or regulatory peptide sequence. Proteins 2009. © 2008 Wiley-Liss, Inc

Mining electron density for functionally relevant protein polysterism in crystal structures.

This review focuses on conceptual and methodological advances in our understanding and characterization of the conformational heterogeneity of proteins. Focusing on X-ray crystallography, we describe how polysterism, the interconversion of pre-existing conformational substates, has traditionally been analyzed by comparing independent crystal structures or multiple chains within a single crystal asymmetric unit. In contrast, recent studies have focused on mining electron density maps to reveal previously 'hidden' minor conformational substates. Functional tests of the importance of minor states suggest that evolutionary selection shapes the entire conformational landscape, including uniquely configured conformational substates, the relative distribution of these substates, and the speed at which the protein can interconvert between them. An increased focus on polysterism may shape the way protein structure and function is studied in the coming years