Proximal effects in the modulation of nitric oxide synthase reactivity: a QM-MM study

Mechanism and regulation of ferrous heme-nitric oxide (NO) oxidation in NO synthases

Nitric oxide (NO) synthases (NOSs) catalyze the formation of NO from l-arginine. We have shown previously that the NOS enzyme catalytic cycle involves a large number of reactions but can be characterized by a global model with three main rate-limiting steps. These are the rate of heme reduction by the flavin domain (k_r ), of dissociation of NO from the ferric heme-NO complex (k_d ), and of oxidation of the ferrous heme-NO complex (k _ox). The reaction of oxygen with the ferrous heme-NO species is part of a futile cycle that does not directly contribute to NO synthesis but allows a population of inactive enzyme molecules to return to the catalytic cycle, and thus, enables a steady-state NO synthesis rate. Previously, we have reported that this reaction does involve the reaction of oxygen with the NO-bound ferrous heme complex, but the mechanistic details of the reaction, that could proceed via either an inner-sphere or an outer-sphere mechanism, remained unclear. Here, we present additional experiments with neuronal NOS (nNOS) and inducible NOS (iNOS) variants (nNOS W409F and iNOS K82A and V346I) and computational methods to study how changes in heme access and electronics affect the reaction. Our results support an inner-sphere mechanism and indicate that the particular heme-thiolate environment of the NOS enzymes can stabilize an N-bound Fe^III-N(O)OO^- intermediate species and thereby catalyze this reaction, which otherwise is not observed or favorable in proteins like globins that contain a histidine-coordinated heme.

Emerging Roles of Nitric Oxide Synthase in Bacterial Physiology

Nitric oxide (NO) is a potent inhibitor of diverse cellular processes in bacteria. Therefore, it was surprising to discover that several bacterial species, primarily Gram-positive organisms, harboured a gene encoding nitric oxide synthase (NOS). Recent attempts to characterize bacterial NOS (bNOS) have resulted in the discovery of structural features that may allow it to function as a NO dioxygenase and produce nitrate in addition to NO. Consistent with this characterization, investigations into the biological function of bNOS have also emphasized a role for NOS-dependent nitrate and nitrite production in aerobic and microaerobic respiration. In this review, we aim to compare, contrast, and summarize the structure, biochemistry, and biological role of bNOS with mammalian NOS and discuss how recent advances in our understanding of bNOS have enabled efforts at designing inhibitors against it.

QM/MM Models of the O2-Evolving Complex of Photosystem II

This paper introduces structural models of the oxygen-evolving complex of photosystem II (PSII) in the dark-stable S1 state, as well as in the reduced S0 and oxidized S2 states, with complete ligation of the metal-oxo cluster by amino acid residues, water, hydroxide, and chloride. The models are developed according to state-of-the-art quantum mechanics/molecular mechanics (QM/MM) hybrid methods, applied in conjunction with the X-ray crystal structure of PSII from the cyanobacterium Thermosynechococcus elongatus, recently reported at 3.5 Å resolution. Manganese and calcium ions are ligated consistently with standard coordination chemistry assumptions, supported by biochemical and spectroscopic data. Furthermore, the calcium-bound chloride ligand is found to be bound in a position consistent with pulsed electron paramagnetic resonance data obtained from acetate-substituted PSII. The ligation of protein ligands includes monodentate coordination of D1-D342, CP43-E354, and D1-D170 to Mn(1), Mn(3), and Mn(4), respectively; η(2) coordination of D1-E333 to both Mn(3) and Mn(2); and ligation of D1-E189 and D1-H332 to Mn(2). The resulting QM/MM structural models are consistent with available mechanistic data and also are compatible with X-ray diffraction models and extended X-ray absorption fine structure measurements of PSII. It is, therefore, conjectured that the proposed QM/MM models are particularly relevant to the development and validation of catalytic water-oxidation intermediates.

Fourier transform infrared spectroscopy study of ligand photodissociation and migration in inducible nitric oxide synthase

Inducible nitric oxide synthase (iNOS) is a homodimeric heme enzyme that catalyzes the formation of nitric oxide (NO) from dioxygen and L-arginine (L-Arg) in a two-step process. The produced NO can either diffuse out of the heme pocket into the surroundings or it can rebind to the heme iron and inhibit enzyme action. Here we have employed Fourier transform infrared (FTIR) photolysis difference spectroscopy at cryogenic temperatures, using the carbon monoxide (CO) and NO stretching bands as local probes of the active site of iNOS. Characteristic changes were observed in the spectra of the heme-bound ligands upon binding of the cofactors. Unlike photolyzed CO, which becomes trapped in well-defined orientations, as indicated by sharp photoproduct bands, photoproduct bands of NO photodissociated from the ferric heme iron were not visible, indicating that NO does not reside in the protein interior in a well-defined location or orientation. This may be favorable for NO release from the enzyme during catalysis because it reduces self-inhibition. Moreover, we used temperature derivative spectroscopy (TDS) with FTIR monitoring to explore the dynamics of NO and carbon monoxide (CO) inside iNOS after photodissociation at cryogenic temperatures. Only a single kinetic photoproduct state was revealed, but no secondary docking sites as in hemoglobins. Interestingly, we observed that intense illumination of six-coordinate ferrous iNOS oxy-NO ruptures the bond between the heme iron and the proximal thiolate to yield five-coordinate ferric iNOS oxy-NO, demonstrating the strong trans effect of the heme-bound NO.

Fourier transform infrared spectroscopy study of ligand photodissociation and migration in inducible nitric oxide synthase

Inducible nitric oxide synthase (iNOS) is a homodimeric heme enzyme that catalyzes the formation of nitric oxide (NO) from dioxygen and L-arginine (L-Arg) in a two-step process. The produced NO can either diffuse out of the heme pocket into the surroundings or it can rebind to the heme iron and inhibit enzyme action. Here we have employed Fourier transform infrared (FTIR) photolysis difference spectroscopy at cryogenic temperatures, using the carbon monoxide (CO) and NO stretching bands as local probes of the active site of iNOS. Characteristic changes were observed in the spectra of the heme-bound ligands upon binding of the cofactors. Unlike photolyzed CO, which becomes trapped in well-defined orientations, as indicated by sharp photoproduct bands, photoproduct bands of NO photodissociated from the ferric heme iron were not visible, indicating that NO does not reside in the protein interior in a well-defined location or orientation. This may be favorable for NO release from the enzyme during catalysis because it reduces self-inhibition. Moreover, we used temperature derivative spectroscopy (TDS) with FTIR monitoring to explore the dynamics of NO and carbon monoxide (CO) inside iNOS after photodissociation at cryogenic temperatures. Only a single kinetic photoproduct state was revealed, but no secondary docking sites as in hemoglobins. Interestingly, we observed that intense illumination of six-coordinate ferrous iNOS oxy-NO ruptures the bond between the heme iron and the proximal thiolate to yield five-coordinate ferric iNOS oxy-NO, demonstrating the strong trans effect of the heme-bound NO.

Multireference and multiconfiguration ab initio methods in heme-related systems: what have we learned so far?

This work reviews the recent applications of ab initio multireference/multiconfiguration (MR/MC) electronic structure methods to heme-related systems, involving tetra-, penta-, and hexa-coordinate species, as well as the high-valent iron-oxo species. The current accuracy of these methods in the various systems is discussed, with special attention to potential sources of systematic errors. Thus, the review summarizes and tries to rationalize the key elements of MR/MC calculations, namely, the choice of the employed active space, especially the so-called double-shell effect that has already been recognized to be important in transition-metal-containing systems, and the impact of these elements on the spin-state energetics of heme species, as well as on the bonding mechanism of small molecules to the heme. It is shown that expansion of the MC wave function into one based on localized orbitals provides a compact and insightful view on some otherwise complex electronic structures. The effects of protein environment on the MR/MC results are summarized for the few available quantum mechanical/molecular mechanical (QM/MM) studies. Comparisons with corresponding DFT results are also made wherever available. Potential future directions are proposed.

Is the bound substrate in nitric oxide synthase protonated or neutral and what is the active oxidant that performs substrate hydroxylation?

We present here results of a series of density functional theory (DFT) studies on enzyme active site models of nitric oxide synthase (NOS) and address the key steps in the catalytic cycle whereby the substrate (L-arginine) is hydroxylated to N(omega)-hydroxo-arginine. It has been proposed that the mechanism follows a cytochrome P450-type catalytic cycle; however, our calculations find an alternative low energy pathway whereby the bound L-arginine substrate has two important functions in the catalytic cycle, namely first as a proton donor and later as the substrate in the reaction mechanism. Thus, the DFT studies show that the oxo-iron active species (compound I) cannot abstract a proton and neither a hydrogen atom from protonated L-arginine due to the strength of the N-H bonds of the substrate. However, the hydroxylation of neutral arginine by compound I and its one electron reduced form (compound II) requires much lower barriers and is highly exothermic. Detailed analysis of proton transfer mechanisms shows that the basicity of the dioxo dianion and the hydroperoxo-iron (compound 0) intermediates in the catalytic cycle are larger than that of arginine, which makes it likely that protonated arginine donates one of the two protons needed during the first catalytic cycle of NOS. Therefore, DFT predicts that in NOS enzymes arginine binds to the active site in its protonated form, but is deprotonated during the oxygen activation process in the catalytic cycle by either the dioxo dianion species or compound 0. As a result of the low ionization potential of neutral arginine, the actual hydroxylation reaction starts with an initial electron transfer from the substrate to compound I to create compound II followed by a concerted hydrogen abstraction/radical rebound from the substrate. These studies indicate that compound II is the actual oxidant in NOS enzymes that performs the hydroxylation reaction of arginine, which is in sharp contrast with the cytochromes P450 where compound II was shown to be a sluggish oxidant. This is the first example of an enzyme where compound II is able to participate in the reaction mechanism. Moreover, arginine hydroxylation by NOS enzymes is catalyzed in a significantly different way from the cytochromes P450 although the active sites of the two enzyme classes are very similar in structure. Detailed studies of environmental effects on the reaction mechanism show that environmental perturbations as appear in the protein have little effect and do not change the energies of the reaction. Finally, a valence bond curve crossing model has been set up to explain the obtained reaction mechanisms for the hydrogen abstraction processes in P450 and NOS enzymes.

QM/MM calculations with DFT for taking into account protein effects on the EPR and optical spectra of metalloproteins. Plastocyanin as a case study

A detailed study of the influence of the surrounding protein on magnetic and optical spectra of metalloproteins is presented using the quantum-mechanical/molecular mechanical (QM/MM) approach. The well-studied type I copper site in plastocyanin in the cupric oxidation state is taken as a test case because its spectroscopic properties have been extensively studied and are well understood. The calculations have been performed using nonrelativistic and scalar relativistic (at the level of the zeroth order regular approximation, ZORA) calculations (B3LYP functional). Linear response theory has been used to calculate first- and second-order properties, namely the EPR g-tensor, the central metal hyperfine couplings (HFCs), the HFCs of the directly coordinating ligands, as well as superhyperfine couplings (1H, 14N) from remote nuclei, transition energies, and oscillator strengths. Two different model systems have been defined that do not and do include important amino acids from the second coordination sphere, respectively. For comparison, calculations have been carried out in the gas phase and in a dielectric continuum (conductor like screening model, COSMO) with a dielectric constant of four. The best results were obtained at the scalar relativistic ZORA level for the largest model in conjunction with explicit modeling of the protein environment through the QM/MM procedure, which is also considered to be the highest level of theory used in this work. The protein effects beyond the second coordination sphere were found to be quite substantial (up to 30% changes on some properties), and were found to require an explicit treatment of the protein beyond the second coordination sphere. In addition, the embedding water cage was found to have a nonnegligible influence on the calculated spectroscopic data, which is of the same order as the influence of the protein backbone charges. However, while qualitatively satisfactory, the errors in the calculated spectroscopic parameters are still substantial, and can all be traced back to the fact that the linear-response of the presently available functionals is "too stiff" with respect to the external perturbations at least for the model systems studied here. Ligand field-based approaches are used to correct for systematic errors in the DFT procedures. As a consequence, we propose a new breakdown of the copper hyperfine interaction into Fermi-contact, spin-dipolar and spin-orbit contributions.

Investigation of the Catalytic Mechanism of Farnesyl Pyrophosphate Synthase by Computer Simulation

Farnesyl pyrophosphate synthase (FPPS) catalyses the formation of a key cellular intermediate in isoprenoid metabolic pathways, farnesyl pyrophosphate, by the sequential head-to-tail condensation of two molecules of isopentenyl diphosphate (IPP) with dimethylallyl diphosphate (DMAPP). Recently, FPPS has been shown to represent an important target for the treatment of parasitic diseases such as Chagas disease and African trypanosomiasis. Bisphosphonates, pyrophosphate analogues in which the oxygen bridge between the two phosphorus atoms has been replaced by a carbon substituted with different side chains, are able to inhibit the FPPS enzyme. Moreover, nitrogen-containing bisphosphonates have been proposed as carbocation transition state analogues of FPPS. On the basis of structural and kinetic data, different catalytic mechanisms have been proposed for FPPS. By analyzing different reaction coordinates we propose that the reaction occurs in one step through a carbocationic transition state and the subsequent transfer of a hydrogen atom from IPP to the pyrophosphate moiety of DMAPP. Moreover, we have analyzed the role of the active site amino acids on the activation barrier and the reaction mechanism. The structure of the active site is well conserved in the isoprenyl diphosphate synthase family; thus, our results are relevant for the understanding of this important class of enzymes and for the design of more potent and specific inhibitors for the treatment of parasitic diseases.