Resolution of ångström-scale protein conformational changes by analyzing fluorescence anisotropy.

Conformational changes within typical protein molecules are rapid and small, making their quantitative resolution challenging. These changes generally involve rotational motions and may thus be resolved by determining changes in the orientation of a fluorescent label that assumes a unique orientation in each conformation. Here, by analyzing fluorescence intensities collected using a polarization microscope at a rate of 50 frames per second, we follow the changes of 10-16° in the orientation of a single bifunctional rhodamine molecule attached to a regulator of conductance to K+ (RCK) domain of the MthK channel, and thus, the transitions between its three conformational states, with effective standard deviation (σ) of 2-5°. Based on available crystal structures, the position of the fluorophore's center differs by 3.4-8.1 Å among the states. Thus, the present approach allows the resolution of protein conformational changes involving ångström-scale displacements.

Energetics of ångström-scale conformational changes in an RCK domain of the MthK K+ channel.

Allosteric proteins transition among different conformational states in a ligand-dependent manner. Upon resolution of a protein's individual states, one can determine the probabilities of these states, thereby dissecting the energetic mechanisms underlying their conformational changes. Here we examine individual regulator of conductance to K+ (RCK) domains that form the regulatory module of the Ca2+-activated MthK channel. Each domain adopts multiple conformational states differing on an ångström scale. The probabilities of these different states of the domain, assessed in different Ca2+ concentrations, allowed us to fully determine a six-state model that is minimally required to account for the energetic characteristics of the Ca2+-dependent conformational changes of an RCK domain. From the energetics of this domain, we deduced, in the framework of statistical mechanics, an analytic model that quantitatively predicts the experimentally observed Ca2+ dependence of the channel's open probability.

Co-cultivation of the anaerobic fungus Anaeromyces robustus with Methanobacterium bryantii enhances transcription of carbohydrate active enzymes.

Anaerobic gut fungi are biomass degraders that form syntrophic associations with other microbes in their native rumen environment. Here, RNA-Seq was used to track and quantify carbohydrate active enzyme (CAZyme) transcription in a synthetic consortium composed of the anaerobic fungus Anaeromyces robustus with methanogen Methanobacterium bryantii. Approximately 5% of total A. robustus genes were differentially regulated in co-culture with M. bryantii relative to cultivation of A. robustus alone. We found that 105 CAZymes (12% of the total predicted CAZymes of A. robustus) were upregulated while 29 were downregulated. Upregulated genes encode putative proteins with a wide array of cellulolytic, xylanolytic, and carbohydrate transport activities; 75% were fused to fungal dockerin domains, associated with a carbohydrate binding module, or both. Collectively, this analysis suggests that co-culture of A. robustus with M. bryantii remodels the transcriptional landscape of CAZymes and associated metabolic pathways in the fungus to aid in lignocellulose breakdown.

[NiFe]-hydrogenases are constitutively expressed in an enriched Methanobacterium sp. population during electromethanogenesis.

Electromethanogenesis is the bioreduction of carbon dioxide (CO2) to methane (CH4) utilizing an electrode as electron donor. Some studies have reported the active participation of Methanobacterium sp. in electron capturing, although no conclusive results are available. In this study, we aimed at determining short-time changes in the expression levels of [NiFe]-hydrogenases (Eha, Ehb and Mvh), heterodisulfide reductase (Hdr), coenzyme F420-reducing [NiFe]-hydrogenase (Frh), and hydrogenase maturation protein (HypD), according to the electron flow in independently connected carbon cloth cathodes poised at- 800 mV vs. standard hydrogen electrode (SHE). Amplicon massive sequencing of cathode biofilm confirmed the presence of an enriched Methanobacterium sp. population (>70% of sequence reads), which remained in an active state (78% of cDNA reads), tagging this archaeon as the main methane producer in the system. Quantitative RT-PCR determinations of ehaB, ehbL, mvhA, hdrA, frhA, and hypD genes resulted in only slight (up to 1.5 fold) changes for four out of six genes analyzed when cells were exposed to open (disconnected) or closed (connected) electric circuit events. The presented results suggested that suspected mechanisms for electron capturing were not regulated at the transcriptional level in Methanobacterium sp. for short time exposures of the cells to connected-disconnected circuits. Additional tests are needed in order to confirm proteins that participate in electron capturing in Methanobacterium sp.

mcrA sequencing reveals the role of basophilic methanogens in a cathodic methanogenic community.

Cathodic methanogenesis is a promising method for accelerating and stabilising bioenergy recovery in anaerobic processes. The change in composition of microbial (especially methanogenic) communities in response to an applied potential-and especially the associated pH gradient-is critical for achieving this goal, but is not well understood in cathodic biofilms. We found here that the pH-polarised region in the 2 mm surrounding the cathode ranged from 6.9 to 10.1, as determined using a pH microsensor; this substantially affected methane production rate as well as microbial community structure. Miseq sequencing data of a highly conserved region of the mcrA gene revealed a dramatic variation in alpha diversity of methanogens concentrated in electrode biofilms under the applied potential, and confirmed that the dominant microbes at the cathode were hydrogenotrophic methanogens (mostly basophilic Methanobacterium alcaliphilum). These results indicate that regional pH variation in the microenvironment surrounding the electrode is an ecological niche enriched with Methanobacterium.

Structural and mutational analysis of archaeal ATP-dependent RNA ligase identifies amino acids required for RNA binding and catalysis.

An ATP-dependent RNA ligase from Methanobacterium thermoautotrophicum (MthRnl) catalyzes intramolecular ligation of single-stranded RNA to form a closed circular RNA via covalent ligase-AMP and RNA-adenylylate intermediate. Here, we report the X-ray crystal structures of an MthRnl•ATP complex as well as the covalent MthRnl-AMP intermediate. We also performed structure-guided mutational analysis to survey the functions of 36 residues in three component steps of the ligation pathway including ligase-adenylylation (step 1), RNA adenylylation (step 2) and phosphodiester bond synthesis (step 3). Kinetic analysis underscored the importance of motif 1a loop structure in promoting phosphodiester bond synthesis. Alanine substitutions of Thr117 or Arg118 favor the reverse step 2 reaction to deadenylate the 5'-AMP from the RNA-adenylate, thereby inhibiting step 3 reaction. Tyr159, Phe281 and Glu285, which are conserved among archaeal ATP-dependent RNA ligases and are situated on the surface of the enzyme, are required for RNA binding. We propose an RNA binding interface of the MthRnl based on the mutational studies and two sulfate ions that co-crystallized at the active site cleft in the MthRnl-AMP complex.

Dynamical study, hydrogen bond analysis, and constant pH simulations of the beta carbonic anhydrase of Methanobacterium thermoautotrophicum.

Within the five classes (α, ß, γ, δ, and ζ) of carbonic anhydrases (CAs) the first two, containing mammal and plant representatives, are the most studied among all CAs. In this study, we have focused our investigation on the beta-class carbonic anhydrase of Methanobacterium thermoautotrophicum. We investigated both the importance of the Asp-Arg dyad near the catalytic zinc-bound water and the possible roles that water molecules within the active site and residues near the entrance of the catalytic cleft have on the first step of the enzyme's reaction mechanism. Hydrogen-bonding analysis of selected residues within the active site and constant pH replica exchange molecular dynamics constant pH replica exchange simulations were performed. The latter was done in order to evaluate the pKa values of possible proton acceptors. We found an intricate hydrogen-bonding network involving two acidic residues within the active site, Asp16 and Asp34, and the catalytic water molecule. We also observed a very strong interaction between the zinc-bound water and residues Asp34 and Arg36. This interaction was not significantly affected by the change in the protonation state of both the catalytic water and aspartate residue 34. The pKa analysis show that the effect of the R36A mutation affects not only the possible proton acceptors, but also the catalytic water itself.

Identification and characterization of functional homologs of nitrogenase cofactor biosynthesis protein NifB from methanogens.

Nitrogenase biosynthesis protein NifB catalyzes the radical S-adenosyl-L-methionine (SAM)-dependent insertion of carbide into the M cluster, the cofactor of the molybdenum nitrogenase from Azotobacter vinelandii. Here, we report the identification and characterization of two naturally "truncated" homologs of NifB from Methanosarcina acetivorans (NifB(Ma)) and Methanobacterium thermoautotrophicum (NifB(Mt)), which contain a SAM-binding domain at the N terminus but lack a domain toward the C terminus that shares homology with NifX, an accessory protein in M cluster biosynthesis. NifB(Ma) and NifB(Mt) are monomeric proteins containing a SAM-binding [Fe4S4] cluster (designated the SAM cluster) and a [Fe4S4]-like cluster pair (designated the K cluster) that can be processed into an [Fe8S9] precursor to the M cluster (designated the L cluster). Further, the K clusters in NifB(Ma) and NifB(Mt) can be converted to L clusters upon addition of SAM, which corresponds to their ability to heterologously donate L clusters to the biosynthetic machinery of A. vinelandii for further maturation into the M clusters. Perhaps even more excitingly, NifB(Ma) and NifB(Mt) can catalyze the removal of methyl group from SAM and the abstraction of hydrogen from this methyl group by 5'-deoxyadenosyl radical that initiates the radical-based incorporation of methyl-derived carbide into the M cluster. The successful identification of NifB(Ma) and NifB(Mt) as functional homologs of NifB not only enabled classification of a new subset of radical SAM methyltransferases that specialize in complex metallocluster assembly, but also provided a new tool for further characterization of the distinctive, NifB-catalyzed methyl transfer and conversion to an iron-bound carbide.

Nicotinamide mononucleotide adenylyltransferase displays alternate binding modes for nicotinamide nucleotides.

Nicotinamide mononucleotide adenylyltransferase (NMNAT) catalyzes the biosynthesis of NAD(+) and NaAD(+). The crystal structure of NMNAT from Methanobacterium thermoautotrophicum complexed with NAD(+) and SO4(2-) revealed the active-site residues involved in binding and catalysis. Site-directed mutagenesis was used to further characterize the roles played by several of these residues. Arg11 and Arg136 were implicated in binding the phosphate groups of the ATP substrate. Both of these residues were mutated to lysine individually. Arg47 does not interact with either NMN or ATP substrates directly, but was deemed to play a role in binding as it is proximal to Arg11 and Arg136. Arg47 was mutated to lysine and glutamic acid. Surprisingly, when expressed in Escherichia coli all of these NMNAT mutants trapped a molecule of NADP(+) in their active sites. This NADP(+) was bound in a conformation that was quite different from that displayed by NAD(+) in the native enzyme complex. When NADP(+) was co-crystallized with wild-type NMNAT, the same structural arrangement was observed. These studies revealed a different conformation of NADP(+) in the active site of NMNAT, indicating plasticity of the active site.

Polynucleotide 3'-terminal phosphate modifications by RNA and DNA ligases.

RNA and DNA ligases catalyze the formation of a phosphodiester bond between the 5'-phosphate and 3'-hydroxyl ends of nucleic acids. In this work, we describe the ability of the thermophilic RNA ligase MthRnl from Methanobacterium thermoautotrophicum to recognize and modify the 3'-terminal phosphate of RNA and single-stranded DNA (ssDNA). This ligase can use an RNA 3'p substrate to generate an RNA 2',3'-cyclic phosphate or convert DNA3'p to ssDNA(3')pp(5')A. An RNA ligase from the Thermus scotoductus bacteriophage TS2126 and a predicted T4 Rnl1-like protein from Thermovibrio ammonificans, TVa, were also able to adenylate ssDNA 3'p. These modifications of RNA and DNA 3'-phosphates are similar to the activities of RtcA, an RNA 3'-phosphate cyclase. The initial step involves adenylation of the enzyme by ATP, which is then transferred to either RNA 3'p or DNA 3'p to generate the adenylated intermediate. For RNA (3')pp(5')A, the third step involves attack of the adjacent 2' hydroxyl to generate the RNA 2',3'-cyclic phosphate. These steps are analogous to those in classical 5' phosphate ligation. MthRnl and TS2126 RNA ligases were not able to modify a 3'p in nicked double-stranded DNA. However, T4 DNA ligase and RtcA can use 3'-phosphorylated nicks in double-stranded DNA to produce a 3'-adenylated product. These 3'-terminal phosphate-adenylated intermediates are substrates for deadenylation by yeast 5'Deadenylase. Our findings that classic ligases can duplicate the adenylation and phosphate cyclization activity of RtcA suggests that they have an essential role in metabolism of nucleic acids with 3'-terminal phosphates.

Quantum and classical simulations of orotidine monophosphate decarboxylase: support for a direct decarboxylation mechanism.

Orotidine 5'-monophosphate (OMP) decarboxylase (ODCase) catalyzes the decarboxylation of OMP to uridine 5'-monophosphate (UMP). Numerous studies of this reaction have suggested a plethora of mechanisms including covalent addition, ylide or carbene formation, and concerted or stepwise protonation. Recent experiments and simulations present strong evidence for a direct decarboxylation mechanism, although direct comparison between experiment and theory is still lacking. In the current work we present hybrid quantum mechanics-molecular mechanics simulations that address the detailed decarboxylation mechanisms for OMP and 5-fluoro-OMP by ODCase. Multidimensional potentials of mean force are computed as functions of structural progress coordinates for the Methanobacterium thermoautotrophicum ODCase reaction: the decarboxylation reaction coordinate, an orbital rehybridization coordinate, and the proton transfer coordinate between Lys72 and the substrate. The computed free energy profiles are in accord with the available experimental data. To facilitate further direct comparison with experiment, we compute the kinetic isotope effects (KIEs) for the enzyme-catalyzed reactions using a mass-perturbation-based path-integral method. The computed KIE provide further support for a direct decarboxylation mechanism. In agreement with experiment, the data suggest a role for Lys72 in stabilizing the transition state in the catalysis of OMP and, to a somewhat lesser extent, in 5-fluoro-OMP.

Structure-function analysis of Methanobacterium thermoautotrophicum RNA ligase - engineering a thermostable ATP independent enzyme.

BACKGROUND: RNA ligases are essential reagents for many methods in molecular biology including NextGen RNA sequencing. To prevent ligation of RNA to itself, ATP independent mutant ligases, defective in self-adenylation, are often used in combination with activated pre-adenylated linkers. It is important that these ligases not have de-adenylation activity, which can result in activation of RNA and formation of background ligation products. An additional useful feature is for the ligase to be active at elevated temperatures. This has the advantage or reducing preferences caused by structures of single-stranded substrates and linkers. RESULTS: To create an RNA ligase with these desirable properties we performed mutational analysis of the archaeal thermophilic RNA ligase from Methanobacterium thermoautotrophicum. We identified amino acids essential for ATP binding and reactivity but dispensable for phosphodiester bond formation with 5' pre-adenylated donor substrate. The motif V lysine mutant (K246A) showed reduced activity in the first two steps of ligation reaction. The mutant has full ligation activity with pre-adenylated substrates but retained the undesirable activity of deadenylation, which is the reverse of step 2 adenylation. A second mutant, an alanine substitution for the catalytic lysine in motif I (K97A) abolished activity in the first two steps of the ligation reaction, but preserved wild type ligation activity in step 3. The activity of the K97A mutant is similar with either pre-adenylated RNA or single-stranded DNA (ssDNA) as donor substrates but we observed two-fold preference for RNA as an acceptor substrate compared to ssDNA with an identical sequence. In contrast, truncated T4 RNA ligase 2, the commercial enzyme used in these applications, is significantly more active using pre-adenylated RNA as a donor compared to pre-adenylated ssDNA. However, the T4 RNA ligases are ineffective in ligating ssDNA acceptors. CONCLUSIONS: Mutational analysis of the heat stable RNA ligase from Methanobacterium thermoautotrophicum resulted in the creation of an ATP independent ligase. The K97A mutant is defective in the first two steps of ligation but retains full activity in ligation of either RNA or ssDNA to a pre-adenylated linker. The ability of the ligase to function at 65°C should reduce the constraints of RNA secondary structure in RNA ligation experiments.

DNA repair and replication fork helicases are differentially affected by alkyl phosphotriester lesion.

DNA helicases are directly responsible for catalytically unwinding duplex DNA in an ATP-dependent and directionally specific manner and play essential roles in cellular nucleic acid metabolism. It has been conventionally thought that DNA helicases are inhibited by bulky covalent DNA adducts in a strand-specific manner. However, the effects of highly stable alkyl phosphotriester (PTE) lesions that are induced by chemical mutagens and refractory to DNA repair have not been previously studied for their effects on helicases. In this study, DNA repair and replication helicases were examined for unwinding a forked duplex DNA substrate harboring a single isopropyl PTE specifically positioned in the helicase-translocating or -nontranslocating strand within the double-stranded region. A comparison of SF2 helicases (RecQ, RECQ1, WRN, BLM, FANCJ, and ChlR1) with a SF1 DNA repair helicase (UvrD) and two replicative helicases (MCM and DnaB) demonstrates unique differences in the effect of the PTE on the DNA unwinding reactions catalyzed by these enzymes. All of the SF2 helicases tested were inhibited by the PTE lesion, whereas UvrD and the replication fork helicases were fully tolerant of the isopropyl backbone modification, irrespective of strand. Sequestration studies demonstrated that RECQ1 helicase was trapped by the PTE lesion only when it resided in the helicase-translocating strand. Our results are discussed in light of the current models for DNA unwinding by helicases that are likely to encounter sugar phosphate backbone damage during biological DNA transactions.

Simple and efficient synthesis of 5' pre-adenylated DNA using thermostable RNA ligase.

We report a simple method of enzymatic synthesis of pre-adenylated DNA linkers/adapters for next-generation sequencing using thermostable RNA ligase from Methanobacterium thermoautotrophicum (MthRnl). Using RNA ligase for the reaction instead of the existing chemical or T4 DNA ligase-based methods allows quantitative conversion of 5'-phosphorylated single-stranded DNA (ssDNA) to the adenylated form. The MthRnl adenylation reaction is specific for ATP and either ssDNA or RNA. In the presence of Mg(+2), the reaction has a pH optimum of 6.0-6.5. Unlike reactions that use T4 DNA ligase, this protocol does not require synthesis of a template strand for adenylation. The high yield of the reaction simplifies isolation and purification of the adenylated product. Conducting the adenylation reaction at the elevated temperature (65°C) reduces structural constraints, while increased ATP concentrations allow quantitative adenylation of DNA with a 3'-unprotected end.

[Cloning of genes, purification and properties investigation of recombinant DNA ligases from the thermophilic archaeon Pyrococcus abyssi and Methanobacterium thermoautotrophicum].

The genes encoding of DNA ligases from the thermophilic archaeon Pyrococcus abyssi (PabDNA ligase) and Methanobacterium thermoautotrophicum (MthDNA ligase) were cloned and expressed in Escherichia coli. The activity of purified enzymes was studied by ligation of two oligonucleotides, one of which had preformed hairpin structure. In the used system the maximal output of reaction products for both DNA ligases was observed near 70 degrees C that is explained by substrate thermostability. At stoichiometric ratio of enzymes and substrate the output of a product reaches of plateau at 70-75% of theoretical ones. Investigated DNA ligases showed different thermostability. The half-time life of PabDNA ligase was about 60 min at 90 degrees C. MthDNA ligase was completely inactivated at this temperature during 10 min. Recombinant DNA ligases from P. abyssi and M. thermoautotrophicum possessed high stability during a storage at 4 degrees C.

Lys314 is a nucleophile in non-classical reactions of orotidine-5'-monophosphate decarboxylase.

Orotidine-5'-monophosphate decarboxylase (OMPD) catalyzes the decarboxylation of orotidine-5'-monophosphate (OMP) to uridine-5'-monophosphate (UMP) in an extremely proficient manner. The reaction does not require any cofactors and proceeds by an unknown mechanism. In addition to decarboxylation, OMPD is able to catalyze other reactions. We show that several C6-substituted UMP derivatives undergo hydrolysis or substitution reactions that depend on a lysine residue (Lys314) in the OMPD active site. 6-Cyano-UMP is converted to UMP, and UMP derivatives with good leaving groups inhibit OMPD by a suicide mechanism in which Lys314 covalently binds to the substrate. These non-classical reactivities of human OMPD were characterized by cocrystallization and freeze-trapping experiments with wild-type OMPD and two active-site mutants by using substrate and inhibitor nucleotides. The structures show that the C6-substituents are not coplanar with the pyrimidine ring. The extent of this substrate distortion is a function of the substituent geometry. Structure-based mechanisms for the reaction of 6-substituted UMP derivatives are extracted in accordance with results from mutagenesis, mass spectrometry, and OMPD enzyme activity. The Lys314-based mechanisms explain the chemodiversity of OMPD, and offer a strategy to design mechanism-based inhibitors that could be used for antineoplastic purposes for example.

Structural characterization of the molecular events during a slow substrate-product transition in orotidine 5'-monophosphate decarboxylase.

Crystal structures of substrate-product complexes of Methanobacterium thermoautotrophicum orotidine 5'-monophosphate decarboxylase, obtained at various steps in its catalysis of the unusual transformation of 6-cyano-uridine 5'-monophosphate (UMP) into barbituric acid ribosyl monophosphate, show that the cyano substituent of the substrate, when bound to the active site, is first bent significantly from the plane of the pyrimidine ring and then replaced by an oxygen atom. Although the K72A and D70A/K72A mutants are either catalytically impaired or even completely inactive, they still display bending of the C6 substituent. Interestingly, high-resolution structures of the D70A and D75N mutants revealed a covalent bond between C6 of UMP and the Lys72 side chain after the -CN moiety's release. The same covalent bond was observed when the native enzyme was incubated with 6-azido-UMP and 6-iodo-UMP; in contrast, the K72A mutant transformed 6-iodo-UMP to barbituric acid ribosyl 5'-monophosphate. These results demonstrate that, given a suitable environment, native orotidine 5'-monophosphate decarboxylase and several of its mutants are not restricted to the physiologically relevant decarboxylation; they are able to catalyze even nucleophilic substitution reactions but consistently maintain distortion on the C6 substituent as an important feature of catalysis.

Archaeal RNA ligase is a homodimeric protein that catalyzes intramolecular ligation of single-stranded RNA and DNA.

RNA ligases participate in repair, splicing and editing pathways that either reseal broken RNAs or alter their primary structure. Here, we report the characterization of an RNA ligase from the thermophilic archaeon, Methanobacterium thermoautotrophicum. The 381-amino acid Methanobacterium RNA ligase (MthRnl) catalyzes intramolecular ligation of 5'-PO(4) single-strand RNA to form a covalently closed circular RNA molecule through ligase-adenylylate and RNA-adenylylate (AppRNA) intermediates. At the optimal temperature of 65 degrees C, AppRNA was predominantly ligated to a circular product. In contrast, at 35 degrees C, phosphodiester bond formation was suppressed and the majority of the AppRNA was deadenylylated. Sedimentation analysis indicates that MthRnl is a homodimer in solution. The C-terminal 127-amino acid segment is required for dimerization, is itself capable of oligomeization and acts in trans to inhibit the ligation activity of native MthRnl. MthRnl can also join single-stranded DNA to form a circular molecule. The lack of specificity for RNA and DNA by MthRnl may exemplify an undifferentiated ancestral stage in the evolution of ATP-dependent ligases.

Mechanism of OMP decarboxylation in orotidine 5'-monophosphate decarboxylase.

Despite extensive experimental and theoretical studies, the detailed catalytic mechanism of orotidine 5'-monophosphate decarboxylase (ODCase) remains controversial. In particular simulation studies using high level quantum mechanics have failed to reproduce experimental activation free energy. One common feature of many previous simulations is that there is a water molecule in the vicinity of the leaving CO2 group whose presence was only observed in the inhibitor bound complex of ODCase/BMP. Various roles have even been proposed for this water molecule from the perspective of stabilizing the transition state and/or intermediate state. We hypothesize that this water molecule is not present in the active ODCase/OMP complex. Based on QM/MM minimum free energy path simulations with accurate density functional methods, we show here that in the absence of this water molecule the enzyme functions through a simple direct decarboxylation mechanism. Analysis of the interactions in the active site indicates multiple factors contributing to the catalysis, including the fine-tuned electrostatic environment of the active site and multiple hydrogen-bonding interactions. To understand better the interactions between the enzyme and the inhibitor BMP molecule, simulations were also carried out to determine the binding free energy of this special water molecule in the ODCase/BMP complex. The results indicate that the water molecule in the active site plays a significant role in the binding of BMP by contributing approximately -3 kcal/mol to the binding free energy of the complex. Therefore, the complex of BMP plus a water molecule, instead of the BMP molecule alone, better represents the tight binding transition state analogue of ODCase. Our simulation results support the direct decarboxylation mechanism and highlight the importance of proper recognition of protein bound water molecules in the protein-ligand binding and the enzyme catalysis.