Evidence of Negative Cooperativity and Half-Site Reactivity within an F420-Dependent Enzyme: Kinetic Analysis of F420H2:NADP(+) Oxidoreductase

Perspectives on Computational Enzyme Modeling: From Mechanisms to Design and Drug Development

Understanding enzyme mechanisms is essential for unraveling the complex molecular machinery of life. In this review, we survey the field of computational enzymology, highlighting key principles governing enzyme mechanisms and discussing ongoing challenges and promising advances. Over the years, computer simulations have become indispensable in the study of enzyme mechanisms, with the integration of experimental and computational exploration now established as a holistic approach to gain deep insights into enzymatic catalysis. Numerous studies have demonstrated the power of computer simulations in characterizing reaction pathways, transition states, substrate selectivity, product distribution, and dynamic conformational changes for various enzymes. Nevertheless, significant challenges remain in investigating the mechanisms of complex multistep reactions, large-scale conformational changes, and allosteric regulation. Beyond mechanistic studies, computational enzyme modeling has emerged as an essential tool for computer-aided enzyme design and the rational discovery of covalent drugs for targeted therapies. Overall, enzyme design/engineering and covalent drug development can greatly benefit from our understanding of the detailed mechanisms of enzymes, such as protein dynamics, entropy contributions, and allostery, as revealed by computational studies. Such a convergence of different research approaches is expected to continue, creating synergies in enzyme research. This review, by outlining the ever-expanding field of enzyme research, aims to provide guidance for future research directions and facilitate new developments in this important and evolving field.

New Insights into the Cooperativity and Dynamics of Dimeric Enzymes

A survey of protein databases indicates that the majority of enzymes exist in oligomeric forms, with about half of those found in the UniProt database being homodimeric. Understanding why many enzymes are in their dimeric form is imperative. Recent developments in experimental and computational techniques have allowed for a deeper comprehension of the cooperative interactions between the subunits of dimeric enzymes. This review aims to succinctly summarize these recent advancements by providing an overview of experimental and theoretical methods, as well as an understanding of cooperativity in substrate binding and the molecular mechanisms of cooperative catalysis within homodimeric enzymes. Focus is set upon the beneficial effects of dimerization and cooperative catalysis. These advancements not only provide essential case studies and theoretical support for comprehending dimeric enzyme catalysis but also serve as a foundation for designing highly efficient catalysts, such as dimeric organic catalysts. Moreover, these developments have significant implications for drug design, as exemplified by Paxlovid, which was designed for the homodimeric main protease of SARS-CoV-2.

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.

Cofactor F420: an expanded view of its distribution, biosynthesis and roles in bacteria and archaea

Many bacteria and archaea produce the redox cofactor F420. F420 is structurally similar to the cofactors FAD and FMN but is catalytically more similar to NAD and NADP. These properties allow F420 to catalyze challenging redox reactions, including key steps in methanogenesis, antibiotic biosynthesis and xenobiotic biodegradation. In the last 5 years, there has been much progress in understanding its distribution, biosynthesis, role and applications. Whereas F420 was previously thought to be confined to Actinobacteria and Euryarchaeota, new evidence indicates it is synthesized across the bacterial and archaeal domains, as a result of extensive horizontal and vertical biosynthetic gene transfer. F420 was thought to be synthesized through one biosynthetic pathway; however, recent advances have revealed variants of this pathway and have resolved their key biosynthetic steps. In parallel, new F420-dependent biosynthetic and metabolic processes have been discovered. These advances have enabled the heterologous production of F420 and identified enantioselective F420H2-dependent reductases for biocatalysis. New research has also helped resolve how microorganisms use F420 to influence human and environmental health, providing opportunities for tuberculosis treatment and methane mitigation. A total of 50 years since its discovery, multiple paradigms associated with F420 have shifted, and new F420-dependent organisms and processes continue to be discovered.

Chloride promotes refolding of active Vibrio alkaline phosphatase through an inactive dimeric intermediate with an altered interface

Most enzymes are homodimers or higher order multimers. Cold‐active alkaline phosphatase from Vibrio splendidus (VAP) transitions into a dimer with very low activity under mild denaturation conditions. The desire to understand why this dimer fails to efficiently catalyse phosphomonoester hydrolysis led us to investigate interfacial communication between subunits. Here, we studied in detail the unfolding mechanism at two pH values and in the presence or absence of sodium chloride. At pH 8.0, the denaturation model had to include an inactive dimer intermediate and follow the pathway: N₂ → I₂ → 2U. At pH 10.5, the model was of a two‐state nature. Enzyme activity was not recovered under several examined refolding conditions. However, in the presence of 0.5 m NaCl, the enzyme was nearly fully reactivated after urea treatment. Thermal inactivation experiments were biphasic where the inactivation could be detected using CD spectroscopy at 190–200 nm. By incorporating a bimane fluorescence probe at the dimer interface, we could monitor inactivation/denaturation at two distinct sites at the dimer interface. A change in bimane fluorescence at both sites was observed during inactivation, but prior to the global unfolding event. Furthermore, the rate of change in bimane fluorescence correlated with inactivation rates at 40 °C. These results indicate and support the hypothesis that the subunits of VAP are only functional in the dimeric state due to the cooperative nature of the reaction mechanism when proper crosstalk between subunits is facilitated.

Effects of isoleucine 135 side chain length on the cofactor donor-acceptor distance within F420H2:NADP+ oxidoreductase: A kinetic analysis

F₄₂₀H₂:NADP⁺ Oxidoreductase (Fno) catalyzes the reversible reduction of NADP⁺ to NADPH by transferring a hydride from the reduced F₄₂₀ cofactor. Here, we have employed binding studies, steady-state and pre steady-state kinetic methods upon wtFno and isoleucine 135 (I135) Fno variants in order to study the effects of side chain length on the donor-acceptor distance between NADP⁺ and the F₄₂₀ precursor, FO. The conserved I135 residue of Fno was converted to a valine, alanine and glycine, thereby shortening the side chain length. The steady-state kinetic analysis of wtFno and the variants showed classic Michaelis-Menten kinetics with varying FO concentrations. The data revealed a decreased k_cat as side chain length decreased, with varying FO concentrations. The steady-state plots revealed non-Michaelis-Menten kinetic behavior when NADPH was varied. The double reciprocal plot of the varying NADPH concentrations displays a downward concave shape, while the NADPH binding curves gave Hill coefficients of less than 1. These data suggest that negative cooperativity occurs between the two identical monomers. The pre steady-state Abs₄₂₀ versus time trace revealed biphasic kinetics, with a fast phase (hydride transfer) and a slow phase. The fast phase displayed an increased rate constant as side chain length decreased. The rate constant for the second phase, remained ~2 s^-1 for each variant. Our data suggest that I135 plays a key role in sustaining the donor-acceptor distance between the two cofactors, thereby regulating the rate at which the hydride is transferred from FOH₂ to NADP⁺. Therefore, Fno is a dynamic enzyme that regulates NADPH production.

Investigating the Reaction Mechanism of F420-Dependent Glucose-6-phosphate Dehydrogenase from Mycobacterium tuberculosis: Kinetic Analysis of the Wild-Type and Mutant Enzymes

F₄₂₀-dependent glucose-6-phosphate dehydrogenase (FGD) catalyzes the conversion of glucose-6-phosphate (G6P) to 6-phosphogluconolactone, using F₄₂₀ cofactor as the hydride transfer acceptor, within mycobacteria. A previous crystal structure of wild-type FGD led to a proposed mechanism suggesting that the active site residues His40, Trp44, and Glu109 could be involved in catalysis. We have characterized the wild-type FGD and five FGD variants (H40A, W44F, W44Y, W44A, and E109Q) by fluorescence binding assays and steady-state and pre-steady-state kinetic experiments. Compared to wild-type FGD, all the variants had lower binding affinities for F₄₂₀, thus suggesting that Trp44, His40, and Glu109 aid in F₄₂₀ binding. While all the variants had decreased catalytic efficiencies, FGD H40A and W44A were the least efficient, having lost ∼1000- and ∼2000-fold activity, respectively. This confirms a crucial catalytic role for His40 in the FGD reaction and suggests that aromaticity at residue 44 aids catalysis. To investigate the proposed roles of Glu109 and His40 in acid-base catalysis, the pH dependence of kinetic parameters has been determined for the E109Q and H40A mutants and compared to those of the wild-type enzyme. The log k_cat-pH profile of wild-type FGD and E109Q revealed two ionizable residues in the enzyme-substrate complex, while H40A displayed only one ionization event. The FGD E109Q variant displayed pH-dependent kinetic cooperativity with respect to the F₄₂₀ cofactor. The multiple-turnover pre-steady-state kinetics were biphasic for wild-type FGD, W44F, W44Y, and E109Q, while the H40A and W44A variants displayed only a single phase because of their reduced catalytic efficiency.

Optimization of Expression and Purification of Recombinant Archeoglobus fulgidus F420H2:NADP+ Oxidoreductase, an F420 Cofactor Dependent Enzyme

Methanogens play a critical role in carbon cycling and contain a number of intriguing biosynthetic pathways. One unusual cofactor found in methanogenic and sulfate reducing archaea is Factor 420 (F420), which can be interconverted between its reduced and oxidized forms by the F420H2:NADP(+) oxidoreductase (Fno) through hydride transfer mechanisms. Here, we report an optimized expression and purification method for recombinant Fno derived from the extreme thermophile Archeoglobus fulgidus. An expression vector that is codon-optimized for heterologous expression in Escherichia coli, modified growth conditions, and a modified purification protocol involving a key polyethyleneimine precipitation step results in a highly purified, homogeneous preparation of Fno that displays high catalytic activity with a truncated F420 analog. This method should accelerate studies on how Fno uses the unusual F420 cofactor during catalysis.