Molecular dynamics simulations of N-terminal peptides from a nucleotide binding protein

Molecular dynamics (MD) simulations of N-terminal peptides from lactate dehydrogenase (LDH) with increasing length and individual secondary structure elements were used to study their stability in relation to folding. Ten simulations of 1–2 ns of different peptides in water starting from the coordinates of the crystal structure were performed. The stability of the peptides was compared qualitatively by analyzing the root mean square deviation (RMSD) from the crystal structure, radius of gyration, secondary and tertiary structure, and solvent accessible surface area. In agreement with earlier MD studies, relatively short (< 15 amino acids) peptides containing individual secondary structure elements were generally found to be unstable; the hydrophobic α₁-helix of the nucleotide binding fold displayed a significantly higher stability, however. Our simulations further showed that the first βαβ supersecondary unit of the characteristic dinucleotide binding fold (Rossmann fold) of LDH is somewhat more stable than other units of similar length and that the α₂-helix, which unfolds by itself, is stabilized by binding to this unit. This finding suggests that the first βαβ unit could function as an N-terminal folding nucleus, upon which the remainder of the polypeptide chain can be assembled. Indeed, simulations with longer units (βαβα and βαβαββ) showed that all structural elements of these units are rather stable. The outcome of our studies is in line with suggestions that folding of the N-terminal portion of LDH in vivo can be a cotranslational process that takes place during the ribosomal peptide synthesis.     

UDP-N-acetylglucosamine 4'-epimerase from the intestinal protozoan Giardia intestinalis lacks UDP-glucose 4'-epimerase activity

The protozoan parasite Giardia intestinalis has a simple life cycle consisting of an intestinal trophozoite stage and an environmentally resistant cyst stage. The cyst is formed when a trophozoite encases itself within an external filamentous covering, the cyst wall, which is crucial to the cyst's survival outside of the host. The filaments in the cyst wall consist mainly of a beta (1-3) polymer of N-acetylgalactosamine. Its precursor, UDP-N-acetylgalactosamine, is synthesized from fructose 6-phosphate by a pathway of five inducible enzymes. The fifth, UDP-N-acetylglucosamine 4'-epimerase, epimerizes UDP-N-acetylglucosamine to UDP-N-acetylgalactosamine reversibly. The epimerase of G. intestinalis lacks UDP-glucose/UDP-galactose 4'-epimerase activity and shows characteristic amino acyl residues to allow binding of only the larger UDP-N-acetylhexosamines. While the Giardia epimerase catalyzes the reversible epimerization of UDP-N-acetylglucosamine to UDP-N-acetylgalactosamine, the reverse reaction apparently is favored. The enzyme has a higher Vmax and a smaller Km in this direction. Therefore, an excess of UDP-N-acetylglucosamine is required to drive the reaction towards the synthesis of UDP-N-acetylgalactosamine, when it is needed for cyst wall formation. This forms the ultimate regulatory step in cyst wall biosynthesis.     

Structure of an antibiotic-synthesizing UDP-glucuronate 4-epimerase MoeE5 in complex with substrate

The epimerase MoeE5 from Streptomyces viridosporus converts UDP-glucuronic acid (UDP-GlcA) to UDP-galacturonic acid (UDP-GalA) to provide the first sugar in synthesizing moenomycin, a potent inhibitor against bacterial peptidoglycan glycosyltransferases. The enzyme belongs to the UDP-hexose 4-epimerase family, and uses NAD⁺ as its cofactor. Here we present the complex crystal structures of MoeE5/NAD⁺/UDP-GlcA and MoeE5/NAD⁺/UDP-glucose, determined at 1.48 Å and 1.66 Å resolution. The cofactor NAD⁺ is bound to the N-terminal Rossmann-fold domain and the substrate is bound to the smaller C-terminal domain. In both crystals the C4 atom of the sugar moiety of the substrate is in close proximity to the C4 atom of the nicotinamide of NAD⁺, and the O4 atom of the sugar is also hydrogen bonded to the side chain of Tyr154, suggesting a productive binding mode. As the first complex structure of this protein family with a bound UDP-GlcA in the active site, it shows an extensive hydrogen-bond network between the enzyme and the substrate. We further built a model with the product UDP-GalA, and found that the unique Arg192 of MoeE5 might play an important role in the catalytic pathway. Consequently, MoeE5 is likely a specific epimerase for UDP-GlcA to UDP-GalA conversion, rather than a promiscuous enzyme as some other family members.     

The crystal structure of the apoenzyme of the iron-sulphur cluster-free hydrogenase

The iron-sulphur cluster-free hydrogenase (Hmd, EC 1.12.98.2) from methanogenic archaea is a novel type of hydrogenase that tightly binds an iron-containing cofactor. The iron is coordinated by two CO molecules, one sulphur and a pyridone derivative, which is linked via a phosphodiester bond to a guanosine base. We report here on the crystal structure of the Hmd apoenzyme from Methanocaldococcus jannaschii at 1.75 A and from Methanopyrus kandleri at 2.4 A resolution. Homodimeric Hmd reveals a unique architecture composed of one central and two identical peripheral globular units. The central unit is composed of the intertwined C-terminal segments of both subunits, forming a novel intersubunit fold. The two peripheral units consist of the N-terminal domain of each subunit. The Rossmann fold-like structure of the N-terminal domain contains a mononucleotide-binding site, which could harbour the GMP moiety of the cofactor. Another binding site for the iron-containing cofactor is most probably Cys176, which is located at the bottom of a deep intersubunit cleft and which has been shown to be essential for enzyme activity. Adjacent to the iron of the cofactor modelled as a ligand to Cys176, an extended U-shaped extra electron density, interpreted as a polyethyleneglycol fragment, suggests a binding site for the substrate methenyltetrahydromethanopterin.     

Computational studies of tryptophanyl-tRNA synthetase: activation of ATP by induced-fit

Catalysis of amino acid activation by Bacillus stearothermophilus tryptophanyl-tRNA synthetase involves three allosteric states: (1) Open; (2) closed pre-transition state (PreTS); and (3) closed products (Product). The interconversions of these states entail significant domain motions driven by ligand binding. We explore the application of molecular dynamics simulations to investigate ligand-linked conformational stability changes associated with this catalytic cycle. Multiple molecular dynamics trajectories (5 ns) for 11 distinct liganded and unliganded monomer configurations show that the PreTS conformation is unstable in the absence of ATP, reverting within approximately 600 ps nearly to the Open conformation. In contrast, Open and Product state trajectories were stable, even without ligands, confirming the previous suggestion that catalysis entails destabilization of the protein conformation, driven by ATP-binding energies developed as the PreTS state assembles during induced-fit. The simulations suggest novel mechanistic details associated with both induced-fit (Open-PreTS) and catalysis (PreTS-Product). Notably, Mg2+ -ATP interactions are coupled to interactions between ATP and active-site lysine side-chains via mechanisms that cannot be captured under the molecular mechanics approximations, and which therefore require restraining potentials for stable simulation. Simulations of Mg2+. ATP-bound PreTS complexes with restraining potentials and with a virtual K111A mutant confirm that these coupling interactions are necessary to sustain the PreTS conformation and, in turn, provide a new model for how the PreTS conformation activates ATP for catalysis. These results emphasize the central role of the PreTS state as a high-energy intermediate structure along the catalytic pathway and suggest that Mg2+ and the KMSKS loop function cooperatively during catalysis.