Molecular Dynamics Simulation of Nitrobenzene Dioxygenase Using AMBER Force Field

Elucidating the Role of O2 Uncoupling for the Adaptation of Bacterial Biodegradation Reactions Catalyzed by Rieske Oxygenases

Oxygenation of aromatic and aliphatic hydrocarbons by Rieske oxygenases is the initial step of various biodegradation pathways for environmental organic contaminants. Microorganisms carrying Rieske oxygenases are able to quickly adapt their substrate spectra to alternative carbon and energy sources that are structurally related to the original target substrate, yet the molecular events responsible for this rapid adaptation are not well understood. Here, we evaluated the hypothesis that reactive oxygen species (ROS) generated by unproductive activation of O2, the so-called O2 uncoupling, in the presence of the alternative substrate exert a selective pressure on the bacterium for increasing the oxygenation efficiency of Rieske oxygenases. To that end, we studied wild-type 2-nitrotoluene dioxygenase from Acidovorax sp. strain JS42 and five enzyme variants that have evolved from adaptive laboratory evolution experiments with 3- and 4-nitrotoluene as alternative growth substrates. The enzyme variants showed a substantially increased oxygenation efficiency toward the new target substrates concomitant with a reduction of ROS production, while mechanisms and kinetics of enzymatic O2 activation remained unchanged. Structural analyses and docking studies suggest that amino acid substitutions in enzyme variants occurred at residues lining both substrate and O2 transport tunnels, enabling tighter binding of the target substrates in the active site. Increased oxygenation efficiencies measured in vitro for the various enzyme (variant)-substrate combinations correlated linearly with in vivo changes in growth rates for evolved Acidovorax strains expressing the variants. Our data suggest that the selective pressure from oxidative stress toward more efficient oxygenation by Rieske oxygenases was most notable when O2 uncoupling exceeded 60%.

Structural and molecular dynamics simulation studies of CBL-interacting protein kinase CIPK and its complexes related to plant salinity stress

CONTEXT

Calcium-dependent signaling in plants is responsible for several major cellular events, including the activation of the salinity-responsive pathways. Calcium binds to calcineurin B-like protein (CBL), and the resulting CBL-Ca2+ complex binds to CBL-interacting protein kinase (CIPK). The CBL-CIPK complex enhances the CIPK interaction with an upstream kinase. The upstream kinase phosphorylates CIPK that, in turn, phosphorylates membrane transporters. Phosphorylation influences transporter activity to kick-start many downstream functions, such as balancing the cytosolic Na+-to-K+ ratio. The CBL-CIPK interaction is pivotal for Ca2+-dependent salinity stress signaling.

METHODS

Computational methods are used to model the entire Arabidopsis thaliana CIPK24 protein structure in its autoinhibited and open-activated states. Arabidopsis thaliana CIPK24-CBL4 complex is predicted based on the protein-protein docking methods. The available structural and functional data support the CIPK24 and the CIPK24-CBL4 complex models. Models are energy-minimized and subjected to molecular dynamics (MD) simulations. MD simulations for 500 ns and 300 ns enabled us to predict the importance of conserved residues of the proteins. Finally, the work is extended to predict the CIPK24-CBL4 complex with the upstream kinase GRIK2. MD simulation for 300 ns on the ternary complex structure enabled us to identify the critical CIPK24-GRIK2 interactions. Together, these data could be used to engineer the CBL-CIPK interaction network for developing salt tolerance in crops.

Semi-rational design of nitroarene dioxygenase for catalytic ability toward 2,4-dichloronitrobenzene

Rieske non-heme dioxygenase family enzymes play an important role in the aerobic biodegradation of nitroaromatic pollutants, but no active dioxygenases are available in nature for initial reactions in the degradation of many refractory pollutants like 2,4-dichloronitrobenzene (24DCNB). Here, we report the engineering of hotspots in 2,3-dichloronitrobenzene dioxygenase from Diaphorobacter sp. strain JS3051, achieved through molecular dynamic simulation analysis and site-directed mutagenesis, with the aim of enhancing its catalytic activity toward 24DCNB. The computationally predicted activity scores were largely consistent with the detected activities in wet experiments. Among them, the two most beneficial mutations (E204M and M248I) were obtained, and the combined mutant reached up to a 62-fold increase in activity toward 24DCNB, generating a single product, 3,5-dichlorocatechol, which is a naturally occurring compound. In silico analysis confirmed that residue 204 affected the substrate preference for meta-substituted nitroarenes, while residue 248 may influence substrate preference by interaction with residue 295. Overall, this study provides a framework for manipulating nitroarene dioxygenases using computational methods to address various nitroarene contamination problems.IMPORTANCEAs a result of human activities, various nitroaromatic pollutants continue to enter the biosphere with poor degradability, and dioxygenation is an important kickoff step to remove toxic nitro-groups and convert them into degradable products. The biodegradation of many nitroarenes has been reported over the decades; however, many others still lack corresponding enzymes to initiate their degradation. Although rieske non-heme dioxygenase family enzymes play extraordinarily important roles in the aerobic biodegradation of various nitroaromatic pollutants, prediction of their substrate specificity is difficult. This work greatly improved the catalytic activity of dioxygenase against 2,4-dichloronitrobenzene by computer-aided semi-rational design, paving a new way for the evolution strategy of nitroarene dioxygenase. This study highlights the potential for using enzyme structure-function information with computational pre-screening methods to rapidly tailor the catalytic functions of enzymes toward poorly biodegradable contaminants.

Unusual catalytic strategy by non-heme Fe(ii)/2-oxoglutarate-dependent aspartyl hydroxylase AspH

Biocatalytic C-H oxidation reactions are of important synthetic utility, provide a sustainable route for selective synthesis of important organic molecules, and are an integral part of fundamental cell processes. The multidomain non-heme Fe(ii)/2-oxoglutarate (2OG) dependent oxygenase AspH catalyzes stereoselective (3R)-hydroxylation of aspartyl- and asparaginyl-residues. Unusually, compared to other 2OG hydroxylases, crystallography has shown that AspH lacks the carboxylate residue of the characteristic two-His-one-Asp/Glu Fe-binding triad. Instead, AspH has a water molecule that coordinates Fe(ii) in the coordination position usually occupied by the Asp/Glu carboxylate. Molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) studies reveal that the iron coordinating water is stabilized by hydrogen bonding with a second coordination sphere (SCS) carboxylate residue Asp721, an arrangement that helps maintain the six coordinated Fe(ii) distorted octahedral coordination geometry and enable catalysis. AspH catalysis follows a dioxygen activation-hydrogen atom transfer (HAT)-rebound hydroxylation mechanism, unusually exhibiting higher activation energy for rebound hydroxylation than for HAT, indicating that the rebound step may be rate-limiting. The HAT step, along with substrate positioning modulated by the non-covalent interactions with SCS residues (Arg688, Arg686, Lys666, Asp721, and Gln664), are essential in determining stereoselectivity, which likely proceeds with retention of configuration. The tetratricopeptide repeat (TPR) domain of AspH influences substrate binding and manifests dynamic motions during catalysis, an observation of interest with respect to other 2OG oxygenases with TPR domains. The results provide unique insights into how non-heme Fe(ii) oxygenases can effectively catalyze stereoselective hydroxylation using only two enzyme-derived Fe-ligating residues, potentially guiding enzyme engineering for stereoselective biocatalysis, thus advancing the development of non-heme Fe(ii) based biomimetic C-H oxidation catalysts, and supporting the proposal that the 2OG oxygenase superfamily may be larger than once perceived.

Effects of Active-Center Reduction of Plant-Type Ferredoxin on Its Structure and Dynamics: Computational Analysis Using Molecular Dynamics Simulations

"Plant-type" ferredoxins (Fds) in the thylakoid membranes of plants, algae, and cyanobacteria possess a single [2Fe-2S] cluster in active sites and mediate light-induced electron transfer from Photosystem I reaction centers to various Fd-dependent enzymes. Structural knowledge of plant-type Fds is relatively limited to static structures, and the detailed behavior of oxidized and reduced Fds has not been fully elucidated. It is important that the investigations of the effects of active-center reduction on the structures and dynamics for elucidating electron-transfer mechanisms. In this study, model systems of oxidized and reduced Fds were constructed from the high-resolution crystal structure of Chlamydomonas reinhardtii Fd1, and three 200 ns molecular dynamics simulations were performed for each system. The force field parameters of the oxidized and reduced active centers were independently obtained using quantum chemical calculations. There were no substantial differences in the global conformations of the oxidized and reduced forms. In contrast, active-center reduction affected the hydrogen-bond network and compactness of the surrounding residues, leading to the increased flexibility of the side chain of Phe61, which is essential for the interaction between Fd and the target protein. These computational results will provide insight into the electron-transfer mechanisms in the Fds.

Elucidating the Role of O2 Uncoupling in the Oxidative Biodegradation of Organic Contaminants by Rieske Non-heme Iron Dioxygenases

Oxygenations of aromatic soil and water contaminants with molecular O₂ catalyzed by Rieske dioxygenases are frequent initial steps of biodegradation in natural and engineered environments. Many of these non-heme ferrous iron enzymes are known to be involved in contaminant metabolism, but the understanding of enzyme–substrate interactions that lead to successful biodegradation is still elusive. Here, we studied the mechanisms of O₂ activation and substrate hydroxylation of two nitroarene dioxygenases to evaluate enzyme- and substrate-specific factors that determine the efficiency of oxygenated product formation. Experiments in enzyme assays of 2-nitrotoluene dioxygenase (2NTDO) and nitrobenzene dioxygenase (NBDO) with methyl-, fluoro-, chloro-, and hydroxy-substituted nitroaromatic substrates reveal that typically 20–100% of the enzyme’s activity involves unproductive paths of O₂ activation with generation of reactive oxygen species through so-called O₂ uncoupling. The ¹⁸O and ¹³C kinetic isotope effects of O₂ activation and nitroaromatic substrate hydroxylation, respectively, suggest that O₂ uncoupling occurs after generation of Fe^III-(hydro)peroxo species in the catalytic cycle. While 2NTDO hydroxylates ortho-substituted nitroaromatic substrates more efficiently, NBDO favors meta-substituted, presumably due to distinct active site residues of the two enzymes. Our data implies, however, that the O₂ uncoupling and hydroxylation activity cannot be assessed from simple structure–reactivity relationships. By quantifying O₂ uncoupling by Rieske dioxygenases, our work provides a mechanistic link between contaminant biodegradation, the generation of reactive oxygen species, and possible adaptation strategies of microorganisms to the exposure of new contaminants.

Substrate-Specific Coupling of O2 Activation to Hydroxylations of Aromatic Compounds by Rieske Non-heme Iron Dioxygenases

Rieske dioxygenases catalyze the initial steps in the hydroxylation of aromatic compounds and are critical for the metabolism of xenobiotic substances. Because substrates do not bind to the mononuclear non-heme Fe^II center, elementary steps leading to O₂ activation and substrate hydroxylation are difficult to delineate, thus making it challenging to rationalize divergent observations on enzyme mechanisms, reactivity, and substrate specificity. Here, we show for nitrobenzene dioxygenase, a Rieske dioxygenase capable of transforming nitroarenes to nitrite and substituted catechols, that unproductive O₂ activation with the release of the unreacted substrate and reactive oxygen species represents an important path in the catalytic cycle. Through correlation of O₂ uncoupling for a series of substituted nitroaromatic compounds with ¹⁸O and ¹³C kinetic isotope effects of dissolved O₂ and aromatic substrates, respectively, we show that O₂ uncoupling occurs after the rate-limiting formation of Fe^III-(hydro)peroxo species from which substrates are hydroxylated. Substituent effects on the extent of O₂ uncoupling suggest that the positioning of the substrate in the active site rather than the susceptibility of the substrate for attack by electrophilic oxygen species is responsible for unproductive O₂ uncoupling. The proposed catalytic cycle provides a mechanistic basis for assessing the very different efficiencies of substrate hydroxylation vs unproductive O₂ activation and generation of reactive oxygen species in reactions catalyzed by Rieske dioxygenases.

Rational Design, Synthesis, and Biological Investigations of N-Methylcarbamoylguanidinyl Azamacrolides as a Novel Chitinase Inhibitor

Chitinase is one of the most important glycoside hydrolyases, widely existing in bacteria, fungi, insects, and plants. It is involved in fungal cell wall remodeling and insect molting. Chitinase inhibitors are an effective means of controlling pathogens and pests. Natural product argifin is a 17-membered pentapeptide that exhibits efficient chitinase inhibitory activity. However, the complexity of the synthetic process results in a lot of restrictions for wide range of applications. In this work, we designed a series of azamacrolide chitinase inhibitors based on the structural features of argifin that have high inhibitory activities against bacterial and insectile chitinase. The most potent chitinase inhibitor compound 19c exhibited IC₅₀ values of 56 nM and 110 nM against OfChi-h and SmChiB, respectively. The molecular docking and molecular dynamics simulations revealed that all inhibitors were bound to the -1 subsite of chitinases via N-methylcarbamoylguanidinyl as well as argifin. Finally, a bioactivity assay against pests was carried out. Compound 18a showed 80% mortality for Mythimna separata at a concentration of 50 mg/L. Besides, insecticides 19b and 19c exhibited high mortality against Plutella xylostella (76 and 73% mortalities at 50 mg/L, respectively).

Promising inhibitors of nsp2 of CHIKV using molecular docking and temperature-dependent molecular dynamics simulations

Infection due to the Chikungunya virus (CHIKV) has taken the life of lots of people; and researchers are working to find the vaccine or promisng drug candidates against this viral infection. In this work, the authors have designed one component reaction based on the thia-/oxa-azolidineone and created a library of 2000 molecules based on the product obtained. Further, the compounds were screened through the docking using iGemdock against the non-structural protein 2 (nsp2) of CHIKV. Molecular docking gives the binding energy (BE) or energy for the formation of the complex between the designed compound and nsp2 of CHIKV; and CMPD222 gave the lowest energy. This is based on the energy obtained from van der Waal's interaction, hydrogen bonding and electrostatic instructions. Further, molecular dynamics simulations (MDS) of nsp2 of CHIKV with and without screened compound (222) were performed to validate the docking results and the change in free energy for the formation of the complex is -10.8327 kcal/mol. To explore the potential of CMPD222, the MDS of the CMPD222-nsp2 of CHIKV were performed at different temperatures (325, 350, 375 and 400 K) to understand the inhibition of the protease. MM-GBSA calculations were performed to determined change in entropy, change in enthalpy and change in free energy to understand the inhibition. Maximum inhibition of nsp2 of CHIKV with CMPD222 is observed at 375 K with a change in free energy of -19.3754 kcal/mol.Communicated by Ramaswamy H. Sarma.

Role of Structural Dynamics in Selectivity and Mechanism of Non-heme Fe(II) and 2-Oxoglutarate-Dependent Oxygenases Involved in DNA Repair

AlkB and its human homologue AlkBH2 are Fe(II)- and 2-oxoglutarate (2OG)-dependent oxygenases that repair alkylated DNA bases occurring as a consequence of reactions with mutagenic agents. We used molecular dynamics (MD) and combined quantum mechanics/molecular mechanics (QM/MM) methods to investigate how structural dynamics influences the selectivity and mechanisms of the AlkB- and AlkBH2-catalyzed demethylation of 3-methylcytosine (m3C) in single (ssDNA) and double (dsDNA) stranded DNA. Dynamics studies reveal the importance of the flexibility in both the protein and DNA components in determining the preferences of AlkB for ssDNA and of AlkBH2 for dsDNA. Correlated motions, including of a hydrophobic β-hairpin, are involved in substrate binding in AlkBH2-dsDNA. The calculations reveal that 2OG rearrangement prior to binding of dioxygen to the active site Fe is preferred over a ferryl rearrangement to form a catalytically productive Fe(IV)=O intermediate. Hydrogen atom transfer proceeds via a σ-channel in AlkBH2-dsDNA and AlkB-dsDNA; in AlkB-ssDNA, there is a competition between σ- and π-channels, implying that the nature of the complexed DNA has potential to alter molecular orbital interactions during the substrate oxidation. Our results reveal the importance of the overall protein-DNA complex in determining selectivity and how the nature of the substrate impacts the mechanism.

Discovery of novel acylhydrazone neuraminidase inhibitors

Neuraminidase (NA) plays a crucial role in the replication and transmission of influenza virus. NA inhibitors have been developed as effective treatments for influenza A and B infections. In this paper, a new lead neuraminidase inhibitor 6a (IC₅₀ = 7.10 ± 0.2 μM) was discovered by ligand-based virtual screening, receptor-based virtual screening, molecular dynamics simulation (MD), and bioassay validation. MD simulation indicates that the morpholinyl group of 6a could be embedded in 430-loop of NA. To exploit the 430-loop in the active site, a series of novel acylhydrazone NA inhibitors 6b-6g were designed and synthesized based on the lead compound 6a. Compound 6e exerts the most potency, with IC₅₀ value of 2.37 ± 0.5 μM against NA, which is lower than that of oseltamivir carboxylate (OC) (IC₅₀ = 3.84 μM). Overall, this work provided unique insights in the discovery of potent inhibitors against NA.

Resolving Discrepancy between Theory and Experiment in 4-Nitrotoluene Oxidation

We have performed calculations of possible oxidation pathways of 4-nitrotoluene (4NT) by permanganate anion and evaluated relative contributions of oxidation of the methyl group and aromatic ring. Although a few theory levels matched the experimental results obtained by compound specific isotope analysis (CSIA) for 4NT, they failed in reproducing results for other nitrotoluene derivatives studied previously [Wijker, R.S.; Adamczyk, P.; Bolotin, J.; Paneth, P.; Hofstetter, T.B. Environ. Sci. Technol., 2013, 47, 13459-13468]. This discrepancy prompted us to reevaluate the experimental isotopic fractionation of carbon and hydrogen for 4NT on which the relative contributions of the oxidation channels has been based. Using position specific isotope analysis (PSIA) for hydrogen isotopic fractionation we have found that the previously determined value was incorrect. Reexamination of theory levels that are in agreement with these new findings indicated that while a better agreement for this particular case can be reached, overall, the previously used B3LYP functional expressed in the 6-31+G(d,p) basis set with inclusion of the polarized continuum model (PCM) of aqueous solution remains the theoretical level of choice in modeling oxidation of nitrotoluene derivatives by permanganate.

Simulation of the Bottleneck Controlling Access into a Rieske Active Site: Predicting Substrates of Naphthalene 1,2-Dioxygenase

Naphthalene 1,2-dioxygenase (NDO) has been computationally understudied despite the extensive experimental knowledge obtained for this enzyme, including numerous crystal structures and over 100 demonstrated substrates. In this study, we have developed a substrate prediction model that moves away from the traditional active-site-centric approach to include the energetics of substrate entry into the active site. By comparison with experimental data, the accuracy of the model for predicting substrate oxidation is 92%, with a positive predictive value of 93% and a negative predictive value of 98%. Also, the present analysis has revealed that the amino acid residues that provided the largest energetic barrier for compounds entering the active site are residues F224, L227, P234, and L235. In addition, F224 is proposed to play a role in controlling ligand entrance via π-π stacking stabilization as well as providing stabilization via T-shaped π-π interactions once the ligand has reached the active-site cavity. Overall, we present a method capable of being scaled to computationally discover thousands of substrates of NDO, and we present parameters to be used for expanding the prediction method to other members of the Rieske non-heme iron oxygenase family.

Substrate and Enzyme Specificity of the Kinetic Isotope Effects Associated with the Dioxygenation of Nitroaromatic Contaminants

Compound-specific isotope analysis (CSIA) is a promising approach for tracking biotransformation of organic pollutants, but isotope fractionation associated with aromatic oxygenations is only poorly understood. We investigated the dioxygenation of a series of nitroaromatic compounds to the corresponding catechols by two enzymes, namely, nitrobenzene and 2-nitrotoluene dioxygenase (NBDO and 2NTDO) to elucidate the enzyme- and substrate-specificity of C and H isotope fractionation. While the apparent (13)C- and (2)H-kinetic isotope effects of nitrobenzene, nitrotoluene isomers, 2,6-dinitrotoluene, and naphthalene dioxygenation by NBDO varied considerably, the correlation of C and H isotope fractionation revealed a common mechanism for nitrobenzene and nitrotoluenes. Similar observations were made for the dioxygenation of these substrates by 2NTDO. Evaluation of reaction kinetics, isotope effects, and commitment-to-catalysis based on experiment and theory showed that rates of dioxygenation are determined by the enzymatic O2 activation and aromatic C oxygenation. The contribution of enzymatic O2 activation to the reaction rate varies for different nitroaromatic substrates of NBDO and 2NTDO. Because aromatic dioxygenation by nonheme iron dioxygenases is frequently the initial step of biodegradation, O2 activation kinetics may also have been responsible for the minor isotope fractionation reported for the oxygenation of other aromatic contaminants.

VR-SCOSMO: A smooth conductor-like screening model with charge-dependent radii for modeling chemical reactions

To better represent the solvation effects observed along reaction pathways, and of ionic species in general, a charge-dependent variable-radii smooth conductor-like screening model (VR-SCOSMO) is developed. This model is implemented and parameterized with a third order density-functional tight binding quantum model, DFTB3/3OB-OPhyd, a quantum method which was developed for organic and biological compounds, utilizing a specific parameterization for phosphate hydrolysis reactions. Unlike most other applications with the DFTB3/3OB model, an auxiliary set of atomic multipoles is constructed from the underlying DFTB3 density matrix which is used to interact the solute with the solvent response surface. The resulting method is variational, produces smooth energies, and has analytic gradients. As a baseline, a conventional SCOSMO model with fixed radii is also parameterized. The SCOSMO and VR-SCOSMO models shown have comparable accuracy in reproducing neutral-molecule absolute solvation free energies; however, the VR-SCOSMO model is shown to reduce the mean unsigned errors (MUEs) of ionic compounds by half (about 2-3 kcal/mol). The VR-SCOSMO model presents similar accuracy as a charge-dependent Poisson-Boltzmann model introduced by Hou et al. [J. Chem. Theory Comput. 6, 2303 (2010)]. VR-SCOSMO is then used to examine the hydrolysis of trimethylphosphate and seven other phosphoryl transesterification reactions with different leaving groups. Two-dimensional energy landscapes are constructed for these reactions and calculated barriers are compared to those obtained from ab initio polarizable continuum calculations and experiment. Results of the VR-SCOSMO model are in good agreement in both cases, capturing the rate-limiting reaction barrier and the nature of the transition state.

Isotope effects of enzymatic dioxygenation of nitrobenzene and 2-nitrotoluene by nitrobenzene dioxygenase

Oxygenation of aromatic rings is a frequent initial step in the biodegradation of persistent contaminants, and the accompanying isotope fractionation is increasingly used to assess the extent of transformation in the environment. Here, we systematically investigated the dioxygenation of two nitroaromatic compounds (nitrobenzene and 2-nitrotoluene) by nitrobenzene dioxygenase (NBDO) to obtain insights into the factors governing its C, H, and N isotope fractionation. Experiments were carried out at different levels of biological complexity from whole bacterial cells to pure enzyme. C, H, and N isotope enrichment factors and kinetic isotope effects (KIEs) were derived from the compound-specific isotope analysis of nitroarenes, whereas C isotope fractionation was also quantified in the oxygenated reaction products. Dioxygenation of nitrobenzene to catechol and 2-nitrotoluene to 3-methylcatechol showed large C isotope enrichment factors, ϵC, of -4.1 ± 0.2‰ and -2.5 ± 0.2‰, respectively, and was observed consistently in the substrates and dioxygenation products. ϵH- and ϵN-values were smaller, that is -5.7 ± 1.3‰ and -1.0 ± 0.3‰, respectively. C isotope fractionation was also identical in experiments with whole bacterial cells and pure enzymes. The corresponding (13)C-KIEs for the dioxygenation of nitrobenzene and 2-nitrotoluene were 1.025 ± 0.001 and 1.018 ± 0.001 and suggest a moderate substrate specificity. Our study illustrates that dioxygenation of nitroaromatic contaminants exhibits a large C isotope fractionation, which is not masked by substrate transport and uptake processes and larger than dioxygenation of other aromatic hydrocarbons.

A DFT and ONIOM study of C–H hydroxylation catalyzed by nitrobenzene 1,2-dioxygenase

The DFT and ONIOM calculations show that C–H hydroxylation by nitrobenzene 1,2-dioxygenase proceeds through a HO–Fe^VO intermediate.