Active Site Protonation State Research Articles

Universality of critical active site glutamate as an acid-base catalyst in serine hydroxymethyltransferase function.

Serine hydroxymethyltransferase (SHMT) is a key enzyme in the one-carbon metabolic pathway, utilizing the vitamin B6 derivative pyridoxal 5'-phosphate (PLP) and vitamin B9 derivative tetrahydrofolate (THF) coenzymes to produce essential biomolecules. Many types of cancer utilize SHMT in metabolic reprogramming, exposing the enzyme as a compelling target for antimetabolite chemotherapies. In pursuit of elucidating the catalytic mechanism of SHMT to aid in the design of SHMT-specific inhibitors, we have used room-temperature neutron crystallography to directly determine the protonation states in a model enzyme Thermus thermophilus SHMT (TthSHMT), which exhibits a conserved active site compared to human mitochondrial SHMT2 (hSHMT2). Here we report the analysis of TthSHMT, with PLP in the internal aldimine form and bound THF-analog, folinic acid (FA), by neutron crystallography to reveal H atom positions in the active site, including PLP and FA. We observed protonated catalytic Glu53 revealing its ability to change protonation state upon FA binding. Furthermore, we obtained X-ray structures of TthSHMT-Gly/FA, TthSHMT-l-Ser/FA, and hSHMT2-Gly/FA ternary complexes with the PLP-Gly or PLP-l-Ser external aldimines to analyze the active site configuration upon PLP reaction with an amino acid substrate and FA binding. Accurate mapping of the active site protonation states together with the structural information gained from the ternary complexes allow us to suggest an essential role of the gating loop conformational changes in the SHMT function and to propose Glu53 as the universal acid-base catalyst in both THF-independent and THF-dependent activities of SHMT.

Insight into the Role of Active Site Protonation States and Water Molecules in the Catalytic Inhibition of DPP4 by Vildagliptin.

Vildagliptin (VIL) is an antidiabetic drug that inhibits dipeptidyl peptidase-4 (DPP4) through a covalent mechanism. The molecular bases for this inhibitory process have been addressed experimentally and computationally. Nevertheless, relevant issues remain unknown such as the roles of active site protonation states and conserved water molecules nearby the catalytic center. In this work, molecular dynamics simulations were applied to examine the structures of 12 noncovalent VIL-DPP4 complexes encompassing all possible protonation states of three noncatalytic residues (His126, Asp663, Asp709) that were inconclusively predicted by different computational tools. A catalytically competent complex structure was only achieved in the system with His126 in its ε-form and nonconventional neutral states for Asp663/Asp709. This complex suggested the involvement of one water molecule in the catalytic process of His740/Ser630 activation, which was confirmed by QM/MM simulations. Our findings support the suitability of a novel water-mediated mechanism in which His740/Ser630 activation occurs concertedly with the nucleophilic attack on VIL and the imidate protonation by Tyr547. Then, the restoration of His740/ Tyr547 protonation states occurs via a two-water hydrogen bonding network in a low-barrier process, thus describing the final step of the catalytic cycle for the first time. Additionally, two hydrolytic mechanisms were proposed based on the hydrogen bonding networks formed by water molecules and the catalytic residues along the inhibitory mechanism. These findings are valuable to unveil the molecular features of the covalent inhibition of DPP4 by VIL and support the future development of novel derivatives with improved structural or mechanistic profiles.

Multiple protonation states in ligand-free SARS-CoV-2 main protease revealed by large-scale quantum molecular dynamics simulations

The main protease (Mpro) in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) catalyzes the cleavage of polyproteins for viral replication. Here, large-scale quantum molecular dynamics and metadynamics simulations for ligand-free Mpro were performed, where all the atoms were treated quantum-mechanically, focusing on elucidation of the controversial active-site protonation state. The simulations clarified that the interconverting multiple protonation states exist in unliganded Mpro, and the catalytically relevant ion-pair state is more stable than the neutral state, which is consistent with neutron crystallography. The results highlight the importance of the ion-pair state for repurposing or discovering antiviral drugs that target Mpro.

Unusual commonality in active site structural features of substrate promiscuous and specialist enzymes

Enzyme promiscuity is the ability of (some) enzymes to perform alternate reactions or catalyze non-cognate substrate(s). The latter is referred to as substrate promiscuity, widely studied for its biotechnological applications and understanding enzyme evolution. Insights into the structural basis of substrate promiscuity would greatly benefit the design and engineering of enzymes. Previous studies on some enzymes have suggested that flexibility, hydrophobicity, and active site protonation state could play an important role in enzyme promiscuity. However, it is not known yet whether substrate promiscuous enzymes have distinctive structural characteristics compared to specialist enzymes, which are specific for a substrate. In pursuit to address this, we have systematically compared substrate/catalytic binding site structural features of substrate promiscuous with those of specialist enzymes. For this, we have carefully constructed dataset of substrate promiscuous and specialist enzymes. On careful analysis, surprisingly, we found that substrate promiscuous and specialist enzymes are similar in various binding/catalytic site structural features such as flexibility, surface area, hydrophobicity, depth, and secondary structures. Recent studies have also alluded that promiscuity is widespread among enzymes. Based on these observations, we propose that substrate promiscuity could be defined as a continuum feature that varies from narrow (specialist) to broad range of substrate preferences. Moreover, diversity of conformational states of an enzyme accessible for ligand binding may possibly regulate its substrate preferences.

Imaging active site chemistry and protonation states: NMR crystallography of the tryptophan synthase α-aminoacrylate intermediate

NMR-assisted crystallography-the integrated application of solid-state NMR, X-ray crystallography, and first-principles computational chemistry-holds significant promise for mechanistic enzymology: by providing atomic-resolution characterization of stable intermediates in enzyme active sites, including hydrogen atom locations and tautomeric equilibria, NMR crystallography offers insight into both structure and chemical dynamics. Here, this integrated approach is used to characterize the tryptophan synthase α-aminoacrylate intermediate, a defining species for pyridoxal-5'-phosphate-dependent enzymes that catalyze β-elimination and replacement reactions. For this intermediate, NMR-assisted crystallography is able to identify the protonation states of the ionizable sites on the cofactor, substrate, and catalytic side chains as well as the location and orientation of crystallographic waters within the active site. Most notable is the water molecule immediately adjacent to the substrate β-carbon, which serves as a hydrogen bond donor to the ε-amino group of the acid-base catalytic residue βLys87. From this analysis, a detailed three-dimensional picture of structure and reactivity emerges, highlighting the fate of the L-serine hydroxyl leaving group and the reaction pathway back to the preceding transition state. Reaction of the α-aminoacrylate intermediate with benzimidazole, an isostere of the natural substrate indole, shows benzimidazole bound in the active site and poised for, but unable to initiate, the subsequent bond formation step. When modeled into the benzimidazole position, indole is positioned with C3 in contact with the α-aminoacrylate Cβ and aligned for nucleophilic attack. Here, the chemically detailed, three-dimensional structure from NMR-assisted crystallography is key to understanding why benzimidazole does not react, while indole does.

Neutron structures of Leishmania mexicana triosephosphate isomerase in complex with reaction-intermediate mimics shed light on the proton-shuttling steps.

Triosephosphate isomerase (TIM) is a key enzyme in glycolysis that catalyses the interconversion of glyceraldehyde 3-phosphate and dihydroxy-acetone phosphate. This simple reaction involves the shuttling of protons mediated by protolysable side chains. The catalytic power of TIM is thought to stem from its ability to facilitate the deprotonation of a carbon next to a carbonyl group to generate an enediolate intermediate. The enediolate intermediate is believed to be mimicked by the inhibitor 2-phosphoglycolate (PGA) and the subsequent enediol intermediate by phosphoglycolohydroxamate (PGH). Here, neutron structures of Leishmania mexicana TIM have been determined with both inhibitors, and joint neutron/X-ray refinement followed by quantum refinement has been performed. The structures show that in the PGA complex the postulated general base Glu167 is protonated, while in the PGH complex it remains deprotonated. The deuteron is clearly localized on Glu167 in the PGA-TIM structure, suggesting an asymmetric hydrogen bond instead of a low-barrier hydrogen bond. The full picture of the active-site protonation states allowed an investigation of the reaction mechanism using density-functional theory calculations.

Direct Observation of Protonation State Modulation in SARS-CoV-2 Main Protease upon Inhibitor Binding with Neutron Crystallography.

The main protease (3CL Mpro) from severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes COVID-19, is an essential enzyme for viral replication with no human counterpart, making it an attractive drug target. To date, no small-molecule clinical drugs are available that specifically inhibit SARS-CoV-2 Mpro. To aid rational drug design, we determined a neutron structure of Mpro in complex with the α-ketoamide inhibitor telaprevir at near-physiological (22 °C) temperature. We directly observed protonation states in the inhibitor complex and compared them with those in the ligand-free Mpro, revealing modulation of the active-site protonation states upon telaprevir binding. We suggest that binding of other α-ketoamide covalent inhibitors can lead to the same protonation state changes in the Mpro active site. Thus, by studying the protonation state changes induced by inhibitors, we provide crucial insights to help guide rational drug design, allowing precise tailoring of inhibitors to manipulate the electrostatic environment of SARS-CoV-2 Mpro.

Active-site protonation states in an acyl-enzyme intermediate of a class A β-lactamase

Active-site protonation states in an acyl-enzyme intermediate of a class A β-lactamase

Active-Site Protonation States in an Acyl-Enzyme Intermediate of a Class A β-Lactamase with a Monobactam Substrate.

ABSTRACTThe monobactam antibiotic aztreonam is used to treat cystic fibrosis patients with chronic pulmonary infections colonized by Pseudomonas aeruginosa strains expressing CTX-M extended-spectrum β-lactamases. The protonation states of active-site residues that are responsible for hydrolysis have been determined previously for the apo form of a CTX-M β-lactamase but not for a monobactam acyl-enzyme intermediate. Here we used neutron and high-resolution X-ray crystallography to probe the mechanism by which CTX-M extended-spectrum β-lactamases hydrolyze monobactam antibiotics. In these first reported structures of a class A β-lactamase in an acyl-enzyme complex with aztreonam, we directly observed most of the hydrogen atoms (as deuterium) within the active site. Although Lys 234 is fully protonated in the acyl intermediate, we found that Lys 73 is neutral. These findings are consistent with Lys 73 being able to serve as a general base during the acylation part of the catalytic mechanism, as previously proposed.

A crystallographic study of the Toho-1 β-lactamase acylation mechanism

β-lactam antibiotics have been used effectively over several decades against many types of highly virulent bacteria. The predominant cause of resistance to these antibiotics in Gram-negative bacterial pathogens is the production of serine β-lactamase enzymes. A key aspect of the class A serine β-lactamase mechanism that remains unresolved and controversial is the identity of the residue acting as the catalytic base during the acylation reaction. Multiple mechanisms have been proposed for the formation of the acyl-enzyme intermediate that are predicated on understanding the protonation states and hydrogen-bonding interactions among the important residues involved in substrate binding and catalysis of these enzymes. For resolving a controversy of this nature surrounding the catalytic mechanism, neutron crystallography is a powerful complement to X-ray crystallography that can explicitly determine the location of deuterium atoms in proteins, thereby directly revealing the hydrogen-bonding interactions of important amino acid residues. Neutron crystallography was used to unambiguously reveal the ground-state active site protonation states and the resulting hydrogen-bonding network in two ligand-free Toho-1 β-lactamase mutants which provided remarkably clear pictures of the active site region prior to substrate binding and subsequent acylation [1,2] and an acylation transition-state analog, benzothiophene-2-boronic acid (BZB), which was also isotopically enriched with 11B. The neutron structure revealed the locations of all deuterium atoms in the active site region and clearly indicated that Glu166 is protonated in the BZB transition-state analog complex. As a result, the complete hydrogen-bonding pathway throughout the active site region could then deduced for this protein-ligand complex that mimics the acylation tetrahedral intermediate [3].

Can Molecular Dynamics and QM/MM Solve the Penicillin Binding Protein Protonation Puzzle?

Benzylpenicillin, a member of the β-lactam antibiotic class, has been widely used to combat bacterial infections since 1947. The general mechanism is well-known: a serine protease enzyme (i.e., DD-peptidase) forms a long lasting intermediate with the lactam ring of the antibiotic known as acylation, effectively preventing biosynthesis of the bacterial cell wall. Despite this overall mechanistic understanding, many details of binding and catalysis are unclear. Specifically, there is ongoing debate about active site protonation states and the role of general acids/bases in the reaction. Herein, a unique combination of MD simulations, QM/MM minimizations, and QM/MM orbital analyses is combined with systematic variation of active site residue protonation states. Critical interactions that maximize the stability of the bound inhibitor are examined and used as metrics. This approach was validated by examining cefoxitin interactions in the CTX-M β-lactamase from E. coli and compared to an ultra high-resolution (0.88 Å) crystal structure. Upon confirming the approach used, an investigation of the preacylated Streptomyces R61 active site with bound benzylpenicillin was performed, varying the protonation states of His298 and Lys65. We concluded that protonated His298 and deprotonated Lys65 are most likely to exist in the R61 active site.

Neutron diffraction studies towards deciphering the protonation state of catalytic residues in the bacterial KDN9P phosphatase.

The enzyme 2-keto-3-deoxy-9-O-phosphonononic acid phosphatase (KDN9P phosphatase) functions in the pathway for the production of 2-keto-3-deoxy-D-glycero-D-galacto-nononic acid, a sialic acid that is important for the survival of commensal bacteria in the human intestine. The enzyme is a member of the haloalkanoate dehalogenase superfamily and represents a good model for the active-site protonation state of family members. Crystals of approximate dimensions 1.5 × 1.0 × 1.0 mm were obtained in space group P2(1)2(1)2, with unit-cell parameters a = 83.1, b = 108.9, c = 75.7 Å. A complete neutron data set was collected from a medium-sized H/D-exchanged crystal at BIODIFF at the Heinz Maier-Leibnitz Zentrum (MLZ), Garching, Germany in 18 d. Initial refinement to 2.3 Å resolution using only neutron data showed significant density for catalytically important residues.

On the Active Site Protonation State in Aspartic Proteases: Implications for Drug Design

Aspartic proteases (AP) are a family of important hydrolytic enzymes in medicinal chemistry, since many of its members have become therapeutical targets for a wide variety of diseases from AIDS to Alzheimer. The enzymatic activity of these proteins is driven by the Asp dyad, a pair of active site Asp residues that participate in the hydrolysis of peptides. Hence, the protonation state of these and other acidic residues present in these enzymes determines the catalytic rate and the affinity for an inhibitor at a given pH. In the present work we have reviewed the effect of the protonation states of the titratable residues in AP's both on catalysis and inhibition in this family of enzymes. The first section focuses on the details of the catalytic reaction mechanism picture brought about by a large number of kinetic, crystallographic and computational chemistry analyses. The results indicate that although the mechanism is similar in both retroviral and eukaryotic enzymes, there are some clear differences. For instance, while in the former family branch the binding of the substrate induces a mono-ionic charge state for the Asp dyad, this charge state seems to be already present in the unbound state of the eukaryotic enzymes. In this section we have explored as well the possible existence of low barrier hydrogen bonds (LBHB's) in the enzymatic path. Catalytic rate enhancement in AP's could in part be explained by the lowering of the barrier for proton transfer in a hydrogen bond from donor to acceptor, which is a typical feature of LBHB's. Review of the published work indicates that the experimental support for this type of bonds is rather scarce and it may be more probable in the first stages of the hydrolytic mechanism in retroviral proteases. The second section deals with the effect of active site protonation state on inhibitor binding. The design of highly potent AP inhibitors, that could be the basis for drug leads require a deep knowledge of the protonation state of the active site residues induced by their presence. This vital issue has been tackled by experimental techniques like NMR, X-ray crystallography, calorimetric and binding kinetic techniques. Recently, we have developed a protocol that combines monitoring the pH effect on binding affinities by SPR methods and rationalization of the results by molecular mechanics based calculations. We have used this combined method on BACE-1 and HIV-1 PR, two important therapeutic targets. Our calculations are able to reproduce the inhibitor binding trends to either enzyme upon a pH increase. The results indicate that inhibitors that differ in the Asp dyad binding fragments will present different binding affinity trends upon a pH increase. Our calculations have enabled us to predict the protonation states at different pH values that underlie the above mentioned trends. We have found out that these results have many implications not only for in silico hit screening campaigns aimed at finding high affinity binders, but also (in the case of BACE-1) for the discovery of cell active compounds.

A QM/MM study on the catalytic mechanism of pyruvate decarboxylase

Pyruvate decarboxylase (PDC) is a typical thiamin diphosphate (ThDP)-dependent enzyme with widespread applications in industry. Though studies regarding the reaction mechanism of PDC have been reported, they are mainly focused on the formation of ThDP ylide and some elementary steps in the catalytic cycle, studies about the whole catalytic cycle of PDC are still not completed. In these previous studies, a major controversy is whether the key active residues (Glu473, Glu50′, Asp27′, His113′, His114′) are protonated or ionized during the reaction. To explore the catalytic mechanism and the role of key residues in the active site, three whole-enzyme models were considered, and the combined QM/MM calculations on the nonoxidative decarboxylation of pyruvate to acetaldehyde catalyzed by PDC were performed. According to our computational results, the fundamental reaction pathways, the complete energy profiles of the whole catalytic cycle, and the specific role of key residues in the common steps were obtained. It is also found that the same residue with different protonation states will lead to different reaction pathways and energy profiles. The mechanism derived from the model in which the residues (Glu473, Glu50′, Asp27′, His113′, His114′) are in their protonated states is most consistent with experimental observations. Therefore, extreme care must be taken when assigning the protonation states in the mechanism study. Because the experimental determination of protonation state is currently difficult, the combined QM/MM method provides an indirect means for determining the active-site protonation state.

Role of the Iron Axial Ligands of Heme Carrier HasA in Heme Uptake and Release

The hemophore protein HasA from Serratia marcescens cycles between two states as follows: the heme-bound holoprotein, which functions as a carrier of the metal cofactor toward the membrane receptor HasR, and the heme-free apoprotein fishing for new porphyrin to be taken up after the heme has been delivered to HasR. Holo- and apo-forms differ for the conformation of the two loops L1 and L2, which provide the axial ligands of the iron through His(32) and Tyr(75), respectively. In the apo-form, loop L1 protrudes toward the solvent far away from loop L2; in the holoprotein, closing of the loops on the heme occurs upon establishment of the two axial coordination bonds. We have established that the two variants obtained via single point mutations of either axial ligand (namely H32A and Y75A) are both in the closed conformation. The presence of the heme and one out of two axial ligands is sufficient to establish a link between L1 and L2, thanks to the presence of coordinating solvent molecules. The latter are stabilized in the iron coordination environment by H-bond interactions with surrounding protein residues. The presence of such a water molecule in both variants is revealed here through a set of different spectroscopic techniques. Previous studies had shown that heme release and uptake processes occur via intermediate states characterized by a Tyr(75)-iron-bound form with open conformation of loop L1. Here, we demonstrate that these states do not naturally occur in the free protein but can only be driven by the interaction with the partner proteins.

Thermodynamic and Kinetic Stabilities of Active Site Protonation States of Class C β-Lactamase

By employing computationally intensive molecular dynamics simulations using hybrid quantum-mechanical/molecular-mechanical approach, we analyze here the kinetic and thermodynamic stabilities of various active site protonation states of a fully solvated class C β-lactamase. We report the detailed mechanism of proton transfer between catalytically important active site residues and the associated free energy barriers. In the apoenzyme, significant structural changes are associated with the proton transfer, and the orientations of active site residues are distinctly different for various protonation states. Among several propositions on the protonation state of the apoprotein, we find that the one with Tyr150 deprotonated and both Lys67 and Lys315 residues being protonated is the most stable one, both thermodynamically and kinetically. However, the equilibrium structure at room temperature is a dynamic one, with Lys315Hζ delocalized between Tyr150Oη and Lys315Nζ. Of great importance, the kinetic and thermodynamic stability of protonation states are significantly affected on noncovalently complexing with cephalothin, an antibiotic molecule. The equilibrium structure of the enzyme-substrate (precovalent) complex has a dynamic protonation state where a proton shuttles frequently between the Tyr150Oη and Lys67Nζ. We examine here the genesis of the manifold change in stability at the molecular level. The importance of our observations toward understanding the reactivity of the enzyme is discussed and experimental observations are rationalized.

The active site protonation states of perdeuterated Toho-1 β-lactamase determined by neutron diffraction support a role for Glu166 as the general base in acylation

Room temperature neutron diffraction data of the fully perdeuterated Toho-1 R274N/R276N double mutant β-lactamase in the apo form were used to visualize deuterium atoms within the active site of the enzyme. This perdeuterated neutron structure of the Toho-1 R274N/R276N reveals the clearest picture yet of the ground-state active site protonation states and the complete hydrogen-bonding network in a β-lactamase enzyme. The ground-state active site protonation states detailed in this neutron diffraction study are consistent with previous high-resolution X-ray studies that support the role of Glu166 as the general base during the acylation reaction in the class A β-lactamase reaction pathway.

How active-site protonation state influences the reactivity and ligation of the heme in chlorite dismutase.

Chlorite dismutase catalyzes O(2) release from chlorite with exquisite efficiency and specificity. The spectroscopic properties, ligand binding affinities, and steady-state kinetics of chlorite dismutase from Dechloromonas aromatica were examined over pH 3-11.5 to gain insight into how the protonation state of the heme environment influences dioxygen formation. An acid-base transition was observed by UV/visible and resonance Raman (rR) spectroscopy with a pK(a) of 8.7, 2-3 pH units below analogous transitions observed in typical His-ligated peroxidases. This transition marks the conversion of a five-coordinate high-spin Fe(III) to a mixed high/low-spin ferric hydroxide, as confirmed by rR spectroscopy. The two Fe-OH stretching frequencies are quite low, consistent with a weak Fe-OH bond, despite the nearly neutral imidazole side chain of the proximal histidine ligand. The hydroxide is proposed to interact strongly with a distal H-bond donor, thereby weakening the Fe-OH bond. The rR spectra of Cld-CO as a function of pH reveal two forms of the complex, one in which there is minimal interaction of distal residues with the carbonyl oxygen and another, acidic form in which the oxygen is under the influence of positive charge. Recent crystallographic data reveal arginine 183 as the lone H-bond-donating residue in the distal pocket. It is likely that this Arg is the strong, positively charged H-bond donor implicated by vibrational data to interact with exogenous axial heme ligands. The same Arg in its neutral (pK(a) approximately 6.5) form also appears to act as the active-site base in binding reactions of protonated ligands, such as HCN, to ferric Cld. The steady-state profile for the rate of chlorite decomposition is characterized by these same pK(a) values. The five-coordinate high-spin acidic Cld is more active than the alkaline hydroxide-bound form. The acid form decomposes chlorite most efficiently when the distal Arg is protonated/cationic (maximum k(cat) = 2.0(+/-0.6) x 10(5) s(-1), k(cat)/K(M) = 3.2(+/-0.4) x 10(7) M(-1) s(-1), pH 5.2, 4 degrees C) and to a somewhat lesser extent when it acts as a H-bond donor to the axial hydroxide ligand under alkaline conditions.

Neutron Diffraction Studies of a Class A β-Lactamase Toho-1 E166A/R274N/R276N Triple Mutant

β-Lactam antibiotics have been used effectively over several decades against many types of bacterial infectious diseases. However, the most common cause of resistance to the β-lactam antibiotics is the production of β-lactamase enzymes that inactivate β-lactams by rapidly hydrolyzing the amide group of the β-lactam ring. Specifically, the class A extended-spectrum β-lactamases (ESBLs) and inhibitor-resistant enzymes arose that were capable of hydrolyzing penicillins and the expanded-spectrum cephalosporins and monobactams in resistant bacteria, which lead to treatment problems in many clinical settings. A more complete understanding of the mechanism of catalysis of these ESBL enzymes will impact current antibiotic drug discovery efforts. Here, we describe the neutron structure of the class A, CTX-M-type ESBL Toho-1 E166A/R274N/R276N triple mutant in its apo form, which is the first reported neutron structure of a β-lactamase enzyme. This neutron structure clearly reveals the active-site protonation states and hydrogen-bonding network of the apo Toho-1 ESBL prior to substrate binding and subsequent acylation. The protonation states of the active-site residues Ser70, Lys73, Ser130, and Lys234 in this neutron structure are consistent with the prediction of a proton transfer pathway from Lys73 to Ser130 that is likely dependent on the conformation of Lys73, which has been hypothesized to be coupled to the protonation state of Glu166 during the acylation reaction. Thus, this neutron structure is in agreement with a proposed mechanism for acylation that identifies Glu166 as the general base for catalysis.

Compound I of Nitric Oxide Synthase: The Active Site Protonation State

A quantum mechanical/molecular mechanical (QM/MM) study of the formation of the elusive active species Compound I (Cpd I) of nitric oxide synthase (NOS) from the oxyferrous intermediate shows that two protons have to be provided to produce a reaction that is reasonably exothermic and that leads to the appearance of a radical on the tetrahydrobiopterin cofactor. Molecular dynamics and energy considerations show that a possible source of proton is the water H-bond chain formed from the surface to the active site, but that a water molecule by itself cannot be the source of the proton; an H3O+ species that is propagated along the chain is more likely. The QM/MM calculations demonstrate that Cpd I and H2O are formed from the ferric-hydrogen peroxide complex in a unique heterolytic O-O cleavage mechanism. The properties of the so-formed Cpd I are compared with those of the known species of chloroperoxidase, and the geometry and spin densities are found to be compatible. The Mössbauer parameters are calculated and may serve as experimental probes in attempts to characterize NOS Cpd I.