Asymmetric ligand binding in homodimeric Enterobacter cloacae nitroreductase yields the Michaelis complex for nitroaromatic substrates

Structure and Dynamics of Three Escherichia coli NfsB Nitro-Reductase Mutants Selected for Enhanced Activity with the Cancer Prodrug CB1954

Escherichia coli NfsB has been studied extensively for its potential for cancer gene therapy by reducing the prodrug CB1954 to a cytotoxic derivative. We have previously made several mutants with enhanced activity for the prodrug and characterised their activity in vitro and in vivo. Here, we determine the X-ray structure of our most active triple and double mutants to date, T41Q/N71S/F124T and T41L/N71S. The two mutant proteins have lower redox potentials than wild-type NfsB, and the mutations have lowered activity with NADH so that, in contrast to the wild-type enzyme, the reduction of the enzyme by NADH, rather than the reaction with CB1954, has a slower maximum rate. The structure of the triple mutant shows the interaction between Q41 and T124, explaining the synergy between these two mutations. Based on these structures, we selected mutants with even higher activity. The most active one contains T41Q/N71S/F124T/M127V, in which the additional M127V mutation enlarges a small channel to the active site. Molecular dynamics simulations show that the mutations or reduction of the FMN cofactors of the protein has little effect on its dynamics and that the largest backbone fluctuations occur at residues that flank the active site, contributing towards its broad substrate range.

The 3D-structure, kinetics and dynamics of the E. coli nitroreductase NfsA with NADP+ provide glimpses of its catalytic mechanism

Nitroreductases activate nitroaromatic antibiotics and cancer prodrugs to cytotoxic hydroxylamines and reduce quinones to quinols. Using steady-state and stopped-flow kinetics, we show that the Escherichia coli nitroreductase NfsA is 20-50 fold more active with NADPH than with NADH and that product release may be rate-limiting. The crystal structure of NfsA with NADP⁺ shows that a mobile loop forms a phosphate-binding pocket. The nicotinamide ring and nicotinamide ribose are mobile, as confirmed in molecular dynamics (MD) simulations. We present a model of NADPH bound to NfsA. Only one NADP⁺ is seen bound to the NfsA dimers, and MD simulations show that binding of a second NADP(H) cofactor is unfavourable, suggesting that NfsA and other members of this protein superfamily may have a half-of-sites mechanism.

NTO Degradation by Nitroreductase: A DFT Study

NTO (5-nitro-1,2,4-triazol-3-one), an energetic material used in military applications, may be released to the environment during manufacturing, transportation, storage, training, and disposal. A detailed investigation of the possible mechanism for all steps of reduction of NTO by oxygen-insensitive nitroreductase, as one of the pathways for NTO environmental degradation, was performed by computational study at the PCM(Pauling)/M06-2X/6-311++G(d,p) level. Obtained results reveal an overall sequence for NTO transformation into ATO (5-amino-1,2,4-triazol-3-one) with the flavin mononucleotide (FMN) cofactor of nitroreductase. Reduction of the nitro group to the nitroso group and the nitroso group to the hydroxylamino group follow a similar mechanism that consists of the sequential electron and proton transfer from the flavin cofactor. The hydride transfer mechanism may contribute to reduction of the nitroso group by the anionic form of the reduced flavin cofactor. Reduction of 5-(hydroxylamino)-1,2,4-triazol-3-one by the neutral form of the reduced flavin is impossible, whereas reduction of the hydroxylamino group to the amino group occurs with the anionic form of the reduced cofactor by a mechanism involving an initial proton transfer from the hydroxonium ion followed by two electrons and one proton transfers from the flavin cofactor. Small activation energies and high exothermicity support the significant contribution of oxygen-insensitive nitroreductase and other enzymes, containing FMN as a cofactor, to NTO degradation in the environment.

Flavin oxidation state impacts on nitrofuran antibiotic binding orientation in nitroreductases

Nitroreductases catalyse the NAD(P)H-dependent nitro reduction in nitrofuran antibiotics, which activates them into cytotoxic molecules leading to cell death. The design of new effective nitrofuran antibiotics relies on knowledge of the kinetic mechanism and nitrofuran binding mode of microbial nitroreductases NfsA and NfsB. This has been hampered by multiple co-crystallisation studies revealing ligand binding in non-electron transfer competent states. In a recent study by Day et al. (2021) the authors investigated the likely reaction mechanism and mode of nitrofurantoin binding to NfsA using potentiometry, global kinetics analysis, crystallography and molecular dynamics simulations. Their findings suggest nitrofurantoin reduction proceeds via a direct hydride transfer from reduced FMN, while the crystallographic binding orientation is an inhibitory complex. Molecular dynamics simulations suggest ligand binding orientations is dependent on the oxidation state of the FMN. This study highlights the importance of utilising computational studies alongside traditional crystallographic approaches, when multiple stable ligand binding orientations can occur.

Systematic Comparison of the Structural and Dynamic Properties of Commonly Used Water Models for Molecular Dynamics Simulations

Water is a unique solvent that is ubiquitous in biology and present in a variety of solutions, mixtures, and materials settings. It therefore forms the basis for all molecular dynamics simulations of biological phenomena, as well as for many chemical, industrial, and materials investigations. Over the years, many water models have been developed, and it remains a challenge to find a single water model that accurately reproduces all experimental properties of water simultaneously. Here, we report a comprehensive comparison of structural and dynamic properties of 30 commonly used 3-point, 4-point, 5-point, and polarizable water models simulated using consistent settings and analysis methods. For the properties of density, coordination number, surface tension, dielectric constant, self-diffusion coefficient, and solvation free energy of methane, models published within the past two decades consistently show better agreement with experimental values compared to models published earlier, albeit with some notable exceptions. However, no single model reproduced all experimental values exactly, highlighting the need to carefully choose a water model for a particular study, depending on the phenomena of interest. Finally, machine learning algorithms quantified the relationship between the water model force field parameters and the resulting bulk properties, providing insight into the parameter-property relationship and illustrating the challenges of developing a water model that can accurately reproduce all properties of water simultaneously.

The structures of E. coli NfsA bound to the antibiotic nitrofurantoin; to 1,4-benzoquinone and to FMN

NfsA is a dimeric flavoprotein that catalyses the reduction in nitroaromatics and quinones by NADPH. This reduction is required for the activity of nitrofuran antibiotics. The crystal structure of free Escherichia coli NfsA and several homologues have been determined previously, but there is no structure of the enzyme with ligands. We present here crystal structures of oxidised E. coli NfsA in the presence of several ligands, including the antibiotic nitrofurantoin. Nitrofurantoin binds with the furan ring, rather than the nitro group that is reduced, near the N5 of the FMN. Molecular dynamics simulations show that this orientation is only favourable in the oxidised enzyme, while potentiometry suggests that little semiquinone is formed in the free protein. This suggests that the reduction occurs by direct hydride transfer from FMNH⁻ to nitrofurantoin bound in the reverse orientation to that in the crystal structure. We present a model of nitrofurantoin bound to reduced NfsA in a viable hydride transfer orientation. The substrate 1,4-benzoquinone and the product hydroquinone are positioned close to the FMN N5 in the respective crystal structures with NfsA, suitable for reaction, but are mobile within the active site. The structure with a second FMN, bound as a ligand, shows that a mobile loop in the free protein forms a phosphate-binding pocket. NfsA is specific for NADPH and a similar conformational change, forming a phosphate-binding pocket, is likely to also occur with the natural cofactor.

Backbone assignment of E. coli NfsB and the effects of addition of the cofactor analogue nicotinic acid

E. coli nitroreductase NfsB (also called NfnB) has been studied extensively, largely due to its potential for cancer gene therapy. A homodimeric flavoprotein of 216 residues, it catalyses the reduction of nitroaromatics to cytotoxic hydroxylamines by NADH and NADPH and also the reduction of quinones to hydroxyquinones. Its role in vivo is not known but it is postulated to be involved in reducing oxidative stress. The crystal structures of the wild type protein and several homologues have been determined in the absence and presence of ligands, including nicotinate as a mimic of the headpiece of the nicotinamide cofactors. There is little effect on the overall structure of the protein on binding ligands, but, from the B factors, there appears to be a decrease in mobility of 2 helices near the active site. As a first step towards examining the dynamics of the protein in solution with and without ligand, we have assigned the backbone ¹³C, ¹⁵N, and ¹H_N resonances of NfsB and examined the effect of the binding of nicotinate on the amide ¹⁵N, and ¹H_N shifts.