The HD domain of the Escherichia coli tRNA nucleotidyltransferase has 2',3'-cyclic phosphodiesterase, 2'-nucleotidase, and phosphatase activities

In all mature tRNAs, the 3'-terminal CCA sequence is synthesized or repaired by a template-independent nucleotidyltransferase (ATP(CTP):tRNA nucleotidyltransferase; EC 2.7.7.25). The Escherichia coli enzyme comprises two domains: an N-terminal domain containing the nucleotidyltransferase activity and an uncharacterized C-terminal HD domain. The HD motif defines a superfamily of metal-dependent phosphohydrolases that includes a variety of uncharacterized proteins and domains associated with nucleotidyltransferases and helicases from bacteria, archaea, and eukaryotes. The C-terminal HD domain in E. coli tRNA nucleotidyltransferase demonstrated Ni(2+)-dependent phosphatase activity toward pyrophosphate, canonical 5'-nucleoside tri- and diphosphates, NADP, and 2'-AMP. Assays with phosphodiesterase substrates revealed surprising metal-independent phosphodiesterase activity toward 2',3'-cAMP, -cGMP, and -cCMP. Without metal or in the presence of Mg(2+), the tRNA nucleotidyltransferase hydrolyzed 2',3'-cyclic substrates with the formation of 2'-nucleotides, whereas in the presence of Ni(2+), the protein also produced some 3'-nucleotides. Mutations at the conserved His-255 and Asp-256 residues comprising the C-terminal HD domain of this protein inactivated both phosphodiesterase and phosphatase activities, indicating that these activities are associated with the HD domain. Low concentrations of the E. coli tRNA (10 nm) had a strong inhibiting effect on both phosphatase and phosphodiesterase activities. The competitive character of inhibition by tRNA suggests that it might be a natural substrate for these activities. This inhibition was completely abolished by the addition of Mg(2+), Mn(2+), or Ca(2+), but not Ni(2+). The data suggest that the phosphohydrolase activities of the HD domain of the E. coli tRNA nucleotidyltransferase are involved in the repair of the 3'-CCA end of tRNA.     

Targeted molecular dynamics simulation studies of binding and conformational changes in E. coli MurD

Enzymes involved in the biosynthesis of bacterial peptidoglycan, an essential cell wall polymer unique to prokaryotic cells, represent a highly interesting target for antibacterial drug design. Structural studies of E. coli MurD, a three-domain ATP hydrolysis driven muramyl ligase revealed two inactive open conformations of the enzyme with a distinct C-terminal domain position. It was hypothesized that the rigid body rotation of this domain brings the enzyme to its closed active conformation, a structure, which was also determined experimentally. Targeted molecular dynamics 1 ns-length simulations were performed in order to examine the substrate binding process and gain insight into structural changes in the enzyme that occur during the conformational transitions into the active conformation. The key interactions essential for the conformational transitions and substrate binding were identified. The results of such studies provide an important step toward more powerful exploitation of experimental protein structures in structure-based inhibitor design.     

Structural studies of E. coli topoisomerase III-DNA complexes reveal a novel type IA topoisomerase-DNA conformational intermediate

Escherichia coli DNA topoisomerase III belongs to the type IA family of DNA topoisomerases, which transiently cleave single-stranded DNA (ssDNA) via a 5' phosphotyrosine intermediate. We have solved crystal structures of wild-type E. coli topoisomerase III bound to an eight-base ssDNA molecule in three different pH environments. The structures reveal the enzyme in three distinct conformational states while bound to DNA. One conformation resembles the one observed previously with a DNA-bound, catalytically inactive mutant of topoisomerase III where DNA binding realigns catalytic residues to form a functional active site. Another conformation represents a novel intermediate in which DNA is bound along the ssDNA-binding groove but does not enter the active site, which remains in a catalytically inactive, closed state. A third conformation shows an intermediate state where the enzyme is still in a closed state, but the ssDNA is starting to invade the active site. For the first time, the active site region in the presence of both the catalytic tyrosine and ssDNA substrate is revealed for a type IA DNA topoisomerase, although there is no evidence of ssDNA cleavage. Comparative analysis of the various conformational states suggests a sequence of domain movements undertaken by the enzyme upon substrate binding.     

Probing conformational plasticity of the activation domain of trypsin: the role of glycine hinges

Trypsin-like serine proteases play essential roles in diverse physiological processes such as hemostasis, apoptosis, signal transduction, reproduction, immune response, matrix remodeling, development, and differentiation. All of these proteases share an intriguing activation mechanism that involves the transition of an unfolded domain (activation domain) of the zymogen to a folded one in the active enzyme. During this conformational change, activation domain segments move around highly conserved glycine hinges. In the present study, hinge glycines were replaced by alanine residues via site directed mutagenesis. The effects of these mutations on the interconversion of the zymogen-like and active conformations as well as on catalytic activity were studied. Mutant trypsins showed zymogen-like structures to varying extents characterized by increased flexibility of some activation domain segments, a more accessible N-terminus and a deformed substrate binding site. Our results suggest that the trypsinogen to trypsin transition is hindered by the mutations, which results in a shift of the equilibrium between the inactive zymogen-like and active enzyme conformations toward the inactive state. Our data also showed, however, that the inactive conformations of the various mutants differ from each other. Binding of substrate analogues shifted the conformational equilibrium toward the active enzyme since inhibited forms of the trypsin mutants showed similar structural features as the wild-type enzyme. The catalytic activity of the mutants correlated with the proper conformation of the active site, which could be supported by varying conformations of the N-terminus and the autolysis loop. Transient kinetic measurements confirmed the existence of an inactive to active conformational transition occurring prior to substrate binding.     

Triggering loops and enzyme function: identification of loops that trigger and modulate movements

Enzyme function often involves a conformational change. There is a general agreement that loops play a vital role in correctly positioning the catalytically important residues. Nevertheless, predicting the functional loops and most importantly their role in enzyme function remains a difficult task. A major reason for this difficulty is that loops that undergo conformational change are frequently not well conserved in their primary sequence. beta1,4-Galactosyltransferase is one such enzyme. There, the amino acid sequence of a long loop that undergoes a large conformational change upon substrate binding is not well conserved. Our molecular dynamics simulations show that the large conformational change in the long loop is brought about by a second, interacting loop. Interestingly, while the structural change of the second loop is much smaller than that of the long loop, its sequence (particularly glycine residues) is highly conserved. We further examine the generality of the proposition that there are loops that trigger movements but nevertheless show little or no structural changes in crystals. We focus on two other enzymes, enolase and lipase. We chose these enzymes, since they too undergo conformational change upon ligand binding, however, they have different folds and different functions. Through multiple sets of simulations we show that the conformational change of the functional loop(s) is brought about through communication of flexibility by triggering loops that have several glycine residues. We further propose that similar to the conservation of common favorable fold types and structural motifs, evolution has also conserved common "skillful" mechanisms. Mechanisms may be conserved across different folds, sequences and functions, with adaptation to specific enzymatic roles.     

The crystal structure of YloQ, a circularly permuted GTPase essential for Bacillus subtilis viability

yloQ is one of 11 essential genes in Bacillus subtilis with unknown roles in the physiology of the cell. It encodes a polypeptide of 298 residues with motifs characteristic of GTPases. As a contribution to elucidating its indispensable cellular function, we have solved the crystal structure of YloQ to 1.6 A spacing, revealing a three-domain organisation. At the heart of the molecule is the putative GTPase domain, which exhibits a classical alpha/beta nucleotide-binding fold with a topology very similar to that of Ras and Era. However, as anticipated from the order in which the conserved G protein motifs appear in the sequence, the GTPase domain fold in YloQ is circularly permuted with respect to the classical GTPases. The nucleotide-binding pocket in YloQ is unoccupied, and analysis of the phosphate-binding (P) loop indicates that conformational changes in this region would be needed to accommodate GTP. The GTPase domain is flanked at its N terminus by a beta-barrel domain with an oligonucleotide/oligosaccharide-binding (OB) fold, and at its C terminus by an alpha-helical domain containing a coordinated zinc ion. This combination of protein modules is unique to YloQ and its orthologues. Sequence comparisons reveal a clustering of conserved basic and aromatic residues on one face of the OB domain, perhaps pointing to a role for YloQ in nucleic acid binding. The zinc ion in the alpha-helical domain is coordinated by three cysteine residues and a histidine residue in a novel ligand organisation. The juxtaposition of the switch I and switch II regions of the G domain and the OB and zinc-binding domains suggests that chemical events at the GTPase active site may be transduced into relative movements of these domains. The pattern of conserved residues and electrostatic surface potential calculations suggest that the OB and/or Zn-binding domains participate in nucleic acid binding consistent with a possible role for YloQ at some stage during mRNA translation.     

Structural switching of Staphylococcus aureus Clp protease: a key to understanding protease dynamics

ATP-dependent Clp protease (ClpP) is an attractive new target for the development of anti-infective agents. The ClpP protease consists of two heptameric rings that enclose a large chamber containing 14 proteolytic active sites. Recent studies indicate that ClpP likely undergoes conformational switching between an extended and degraded active state required for substrate proteolysis and a compacted and catalytically inactive state allowing product release. Here, we present the wild-type ClpP structures in two distinct states from Staphylococcus aureus. One structure is very similar to those solved ClpP structures in the extended states. The other is strikingly different from both the extended and the compacted state as observed in ClpP from other species; the handle domain of this structure kinks to take on a compressed conformation. Structural analysis and molecular dynamic simulations show that the handle domain predominantly controls the way in which degradation products exit the chamber through dynamic conformational switching from the extended state to the compressed state. Given the highly conserved sequences among ClpP from different species, this compressed conformation is unexpected and novel, which is potentially valuable for understanding the enzymatic dynamics and the acting mechanisms of ClpP.     

Substrate Binding Specifically Modulates Domain Arrangements in Adenylate Kinase

The enzyme adenylate kinase (ADK) features two substrate binding domains that undergo large-scale motions during catalysis. In the apo state, the enzyme preferentially adopts a globally open state with accessible binding sites. Binding of two substrate molecules (AMP + ATP or ADP + ADP) results in a closed domain conformation, allowing efficient phosphoryl-transfer catalysis. We employed molecular dynamics simulations to systematically investigate how the individual domain motions are modulated by the binding of substrates. Two-dimensional free-energy landscapes were calculated along the opening of the two flexible lid domains for apo and holo ADK as well as for all single natural substrates bound to one of the two binding sites of ADK. The simulations reveal a strong dependence of the conformational ensembles on type and binding position of the bound substrates and a nonsymmetric behavior of the lid domains. Altogether, the ensembles suggest that, upon initial substrate binding to the corresponding lid site, the opposing lid is maintained open and accessible for subsequent substrate binding. In contrast, ATP binding to the AMP-lid induces global domain closing, preventing further substrate binding to the ATP-lid site. This might constitute a mechanism by which the enzyme avoids the formation of a stable but enzymatically unproductive state.     

Molecular dynamics simulation reveals preorganization of the chloroplast FtsY towards complex formation induced by GTP binding

Two GTPases in the signal recognition particle (SRP) and SRP receptor (SR) interact with one another to mediate the cotranslational protein targeting pathway. Previous studies have shown that a universally conserved SRP RNA facilitates an efficient SRP–SR interaction in the presence of a signal sequence bound to SRP. However, a remarkable exception has been found in chloroplast SRP (cpSRP) pathway, in which the SRP RNA is missing. Based on biochemical and structural analyses, it is proposed that free cpSRP receptor (cpFtsY) has already been preorganized into a closed state for efficient cpSRP–cpFtsY association. However, no direct evidence has been reported to support this postulation thus far. In this study, we characterized the structural dynamics of cpFtsY and its conformational rearrangements induced by GTP binding using molecular dynamics (MD) simulations. Our results showed that the GTP-binding event triggered substantial conformational changes in free cpFtsY, including the relative orientation of N–G domain and several conserved motifs that are critical in complex formation. These rearrangements enabled the cpFtsY to relax into a preorganized ‘closed’ state that favored the formation of a stable complex with cpSRP54. Interestingly, the intrinsic flexibility of αN1 helix facilitated these rearrangements. In addition, GTP binding in cpFtsY was mediated by conserved residues that have been shown in other SRP GTPases. These findings suggested that GTP-bound cpFtsY could fluctuate into conformations that are favorable to form the stable complex, providing explanation of why SRP–SR interaction bypasses the requirement of the SRP RNA at a molecular level.     

Molecular dynamics simulation of the interaction between the complex iron-sulfur flavoprotein glutamate synthase and its substrates

Glutamate synthase (GltS) is a complex iron-sulfur flavoprotein that catalyzes the reductive transfer of L-glutamine amide group to the C2 carbon of 2-oxoglutarate yielding two molecules of L-glutamate. Molecular dynamics calculations in explicit solvent were carried out to gain insight into the conformational flexibility of GltS and into the role played by the enzyme substrates in regulating the catalytic cycle. We have modelled the free (unliganded) form of Azospirillum brasilense GltS alpha subunit and the structure of the reduced enzyme in complex with the L-glutamine and 2-oxoglutarate substrates starting from the crystallographically determined coordinates of the GltS alpha subunit in complex with L-methionine sulphone and 2-oxoglutarate. The present 4-ns molecular dynamics calculations reveal that the GltS glutaminase site may exist in a catalytically inactive conformation unable to bind glutamine, and in a catalytically competent conformation, which is stabilized by the glutamine substrate. Substrates binding also induce (1) closure of the loop formed by residues 263-271 with partial shielding of the glutaminase site from solvent, and (2) widening of the ammonia tunnel entrance at the glutaminase end to allow for ammonia diffusion toward the synthase site. The Q-loop of glutamate synthase, which acts as an active site lid in other amidotransferases, seems to maintain an open conformation. Finally, binding of L-methionine sulfone, a glutamine analog that mimics the tetrahedral transient species occurring during its hydrolysis, causes a coordinated rigid-body motion of segments of the glutaminase domain that results in the inactive conformation observed in the crystal structure of GltS alpha subunit.     

Combining protein sequence,structure, and dynamics: A novel approach for functional evolution analysis of PAS domain superfamily

PAS domains are widespread in archaea, bacteria, and eukaryota, and play important roles in various functions. In this study, we aim to explore functional evolutionary relationship among proteins in the PAS domain superfamily in view of the sequence‐structure‐dynamics‐function relationship. We collected protein sequences and crystal structure data from RCSB Protein Data Bank of the PAS domain superfamily belonging to three biological functions (nucleotide binding, photoreceptor activity, and transferase activity). Protein sequences were aligned and then used to select sequence‐conserved residues and build phylogenetic tree. Three‐dimensional structure alignment was also applied to obtain structure‐conserved residues. The protein dynamics were analyzed using elastic network model (ENM) and validated by molecular dynamics (MD) simulation. The result showed that the proteins with same function could be grouped by sequence similarity, and proteins in different functional groups displayed statistically significant difference in their vibrational patterns. Interestingly, in all three functional groups, conserved amino acid residues identified by sequence and structure conservation analysis generally have a lower fluctuation than other residues. In addition, the fluctuation of conserved residues in each biological function group was strongly correlated with the corresponding biological function. This research suggested a direct connection in which the protein sequences were related to various functions through structural dynamics. This is a new attempt to delineate functional evolution of proteins using the integrated information of sequence, structure, and dynamics.     

Structure and dynamics of pin1 during catalysis by NMR

The link between internal enzyme motions and catalysis is poorly understood. Correlated motions in the microsecond-to-millisecond timescale may be critical for enzyme function. We have characterized the backbone dynamics of the peptidylprolyl isomerase (Pin1) catalytic domain in the free state and during catalysis. Pin1 is a prolyl isomerase of the parvulin family and specifically catalyzes the isomerization of phosphorylated Ser/Thr-Pro peptide bonds. Pin1 has been shown to be essential for cell-cycle progression and to interact with the neuronal tau protein inhibiting its aggregation into fibrillar tangles as found in Alzheimer's disease. (15)N relaxation dispersion measurements performed on Pin1 during catalysis reveal conformational exchange processes in the microsecond timescale. A subset of active site residues undergo kinetically similar exchange processes even in the absence of a substrate, suggesting that this area is already "primed" for catalysis. Furthermore, structural data of the turning-over enzyme were obtained through inter- and intramolecular nuclear Overhauser enhancements. This analysis together with a characterization of the substrate concentration dependence of the conformational exchange allowed the distinguishing of regions of the enzyme active site that are affected primarily by substrate binding versus substrate isomerization. Together these data suggest a model for the reaction trajectory of Pin1 catalysis.     

Conformational plasticity and dynamics in the generic protein folding catalyst SlyD unraveled by single-molecule FRET

The relation between conformational dynamics and chemistry in enzyme catalysis recently has received increasing attention. While, in the past, the mechanochemical coupling was mainly attributed to molecular motors, nowadays, it seems that this linkage is far more general. Single-molecule fluorescence methods are perfectly suited to directly evidence conformational flexibility and dynamics. By labeling the enzyme SlyD, a member of peptidyl-prolyl cis-trans isomerases of the FK506 binding protein type with an inserted chaperone domain, with donor and acceptor fluorophores for single-molecule fluorescence resonance energy transfer, we directly monitor conformational flexibility and conformational dynamics between the chaperone domain and the FK506 binding protein domain. We find a broad distribution of distances between the labels with two main maxima, which we attribute to an open conformation and to a closed conformation of the enzyme. Correlation analysis demonstrates that the conformations exchange on a rate in the 100 Hz range. With the aid from Monte Carlo simulations, we show that there must be conformational flexibility beyond the two main conformational states. Interestingly, neither the conformational distribution nor the dynamics is significantly altered upon binding of substrates or other known binding partners. Based on these experimental findings, we propose a model where the conformational dynamics is used to search the conformation enabling the chemical step, which also explains the remarkable substrate promiscuity connected with a high efficiency of this class of peptidyl-prolyl cis-trans isomerases.     

Redesigning the PheA domain of gramicidin synthetase leads to a new understanding of the enzyme's mechanism and selectivity

The PheA domain of gramicidin synthetase A, a non-ribosomal peptide synthetase, selectively binds phenylalanine along with ATP and Mg2+ and catalyzes the formation of an aminoacyl adenylate. In this study, we have used a novel protein redesign algorithm, K*, to predict mutations in PheA that should exhibit improved binding for tyrosine. Interestingly, the introduction of two predicted mutations to PheA did not significantly improve KD, as measured by equilibrium fluorescence quenching. However, the mutations improved the specificity of the enzyme for tyrosine (as measured by kcat/KM), primarily driven by a 56-fold improvement in KM, although the improvement did not make tyrosine the preferred substrate over phenylalanine. Using stopped-flow fluorometry, we examined binding of different amino acid substrates to the wild-type and mutant enzymes in the pre-steady state in order to understand the improvement in KM. Through these investigations, it became evident that substrate binding to the wild-type enzyme is more complex than previously described. These experiments show that the wild-type enzyme binds phenylalanine in a kinetically selective manner; no other amino acids tested appeared to bind the enzyme in the early time frame examined (500 ms). Furthermore, experiments with PheA, phenylalanine, and ATP reveal a two-step binding process, suggesting that the PheA-ATP-phenylalanine complex may undergo a conformational change toward a catalytically relevant intermediate on the pathway to adenylation; experiments with PheA, phenylalanine, and other nucleotides exhibit only a one-step binding process. The improvement in KM for the mutant enzyme toward tyrosine, as predicted by K*, may indicate that redesigning the side-chain binding pocket allows the substrate backbone to adopt productive conformations for catalysis but that further improvements may be afforded by modeling an enzyme:ATP:substrate complex, which is capable of undergoing conformational change.     

Structure of the adenylation domain of an NAD+-dependent DNA ligase.

BACKGROUND: DNA ligases catalyse phosphodiester bond formation between adjacent bases in nicked DNA, thereby sealing the nick. A key step in the catalytic mechanism is the formation of an adenylated DNA intermediate. The adenyl group is derived from either ATP (in eucaryotes and archaea) or NAD+4 (in bacteria). This difference in cofactor specificity suggests that DNA ligase may be a useful antibiotic target. RESULTS: The crystal structure of the adenylation domain of the NAD+-dependent DNA ligase from Bacillus stearothermophilus has been determined at 2.8 A resolution. Despite a complete lack of detectable sequence similarity, the fold of the central core of this domain shares homology with the equivalent region of ATP-dependent DNA ligases, providing strong evidence for the location of the NAD+-binding site. CONCLUSIONS: Comparison of the structure of the NAD+4-dependent DNA ligase with that of ATP-dependent ligases and mRNA-capping enzymes demonstrates the manifold utilisation of a conserved nucleotidyltransferase domain within this family of enzymes. Whilst this conserved core domain retains a common mode of nucleotide binding and activation, it is the additional domains at the N terminus and/or the C terminus that provide the alternative specificities and functionalities in the different members of this enzyme superfamily.     

X-Ray structures of the universal translation initiation factor IF2/eIF5B: conformational changes on GDP and GTP binding

X-ray structures of the universal translation initiation factor IF2/eIF5B have been determined in three states: free enzyme, inactive IF2/eIF5B.GDP, and active IF2/eIF5B.GTP. The "chalice-shaped" enzyme is a GTPase that facilitates ribosomal subunit joining and Met-tRNA(i) binding to ribosomes in all three kingdoms of life. The conserved core of IF2/eIF5B consists of an N-terminal G domain (I) plus an EF-Tu-type beta barrel (II), followed by a novel alpha/beta/alpha-sandwich (III) connected via an alpha helix to a second EF-Tu-type beta barrel (IV). Structural comparisons reveal a molecular lever, which amplifies a modest conformational change in the Switch 2 region of the G domain induced by Mg(2+)/GTP binding over a distance of 90 A from the G domain active center to domain IV. Mechanisms of GTPase function and ribosome binding are discussed.     

Mobility in the structure of E.coli recQ helicase upon substrate binding as seen from molecular dynamics simulations

RecQ helicases feature multiple domains in their structure, of which the helicase domain, the RecQ-Ct domain and the HRDC domains are well conserved among the SF2 helicases. The helicase domain and the RecQ-Ct domain constitute the catalytic core of the enzyme. The domain interfaces are the DNA binding sites which display significant conformational changes in our molecular dynamics simulation studies. The preferred conformational states of the DNA bound and unbound forms of RecQ appear to be quite different from each other. DNA binding induces inter-domain flexibility leading to hinge mobility between the domains. The divergence in the dynamics of the two structures is caused by changes in the interactions at the domain interface, which seems to propagate along the whole protein structure. This could be essential in ssDNA binding after strand separation, as well as aiding translocation of the RecQ protein like an inch-worm.     

Structure of the N-terminal domain of Escherichia coli glutamine synthetase adenylyltransferase

We report the crystal structure of the N-terminal domain of Escherichia coli adenylyltransferase that catalyzes the reversible nucleotidylation of glutamine synthetase (GS), a key enzyme in nitrogen assimilation. This domain (AT-N440) catalyzes the deadenylylation and subsequent activation of GS. The structure has been divided into three subdomains, two of which bear some similarity to kanamycin nucleotidyltransferase (KNT). However, the orientation of the two domains in AT-N440 differs from that in KNT. The active site of AT-N440 has been identified on the basis of structural comparisons with KNT, DNA polymerase beta, and polyadenylate polymerase. AT-N440 has a cluster of metal binding residues that are conserved in polbeta-like nucleotidyl transferases. The location of residues conserved in all ATase sequences was found to cluster around the active site. Many of these residues are very likely to play a role in catalysis, substrate binding, or effector binding.     

Multiple substrate binding states and chiral recognition in cofactor-independent glutamate racemase: a molecular dynamics study

Glutamate racemase (MurI) catalyzes the racemization of glutamate; two cysteine residues serve as catalytic acid and base. On the basis of the crystal structure of MurI from the hyperthermophilic bacterium Aquifex pyrophilus, we performed molecular dynamics (MD) simulations of six different systems to investigate stereochemistry, substrate ligation, and active site protonation state. The catalytic competence of individual systems was assessed by the abundance of reactive conformers. Only systems in which Cys70 is poised to deprotonate d-Glu were found to be catalytically competent (idem Cys178/l-Glu), in agreement with the experimentally observed stereochemistry of Lactobacillus fermentii MurI [Tanner, M. E. et al. (1993) Biochemistry 32, 3998-4006]. Only systems in which the alpha-amino group of l/d-Glu and the imidazole moiety of His are deprotonated are catalytically competent. The active site of MurI displays an unusual flexibility in substrate ligation, and several transitions between stable binding patterns were observed. In catalytically competent binding states, the conserved threonine residues 72, 114, and 117 ligate the alpha-carboxylate of Glu and the Asn71 amides ligate the alpha-amino group of Glu, whereas the delta-carboxylate of Glu is steered by electrostatic repulsion from the Asp7 and Glu147 side chain carboxylates. A network of hydrogen bonds controls the positioning of each thiol/thiolate. In what we term substrate flipping, Glu suddenly rotates into a binding pattern that resembles the post-racemization state of the other enantiomer, i.e., each enantiomer can be bound in two distinct states. Substrate flipping and unfavorable substrate binding successively trigger dissociation of the substrate, accompanied by an opening of the active site channel. We explain how the weak binding of Glu contributes to catalysis and suggest a mechanism by which binding mismatches are propagated into an opening of the active site.     

Molecular dynamics simulations of viral RNA polymerases link conserved and correlated motions of functional elements to fidelity

The viral RNA-dependent RNA polymerase (RdRp) is essential for multiplication of all RNA viruses. The sequence diversity of an RNA virus population contributes to its ability to infect the host. This diversity emanates from errors made by the RdRp during RNA synthesis. The physical basis for RdRp fidelity is unclear but is linked to conformational changes occurring during the nucleotide-addition cycle. To understand RdRp dynamics that might influence RdRp function, we have analyzed all-atom molecular dynamics simulations on the nanosecond timescale of four RdRps from the picornavirus family that exhibit 30-74% sequence identity. Principal component analysis showed that the major motions observed during the simulations derived from conserved structural motifs and regions of known function. The dynamics of residues participating in the same biochemical property, for example, RNA binding, nucleotide binding or catalysis, were correlated even when spatially distant on the RdRp structure. The conserved and correlated dynamics of functional structural elements suggest coevolution of dynamics with structure and function of the RdRp. Crystal structures of all picornavirus RdRps exhibit a template-nascent RNA duplex channel too small to fully accommodate duplex RNA. Simulations revealed opening and closing motions of the RNA and nucleoside triphosphate channels, which might be relevant to nucleoside triphosphate entry, inorganic pyrophosphate exit and translocation. A role for nanosecond timescale dynamics in RdRp fidelity is supported by the altered dynamics of the high-fidelity G64S derivative of PV RdRp relative to wild-type enzyme.