Hybrid molecular dynamics-quantum mechanics simulations of solute spectral properties in the condensed phase: evaluation of simulation parameters

We have explored the impact of a number of basic simulation parameters on the results of a recently developed hybrid molecular dynamics-quantum mechanics (MD-QM) method (Mercer et al., J Phys Chem B 1999, 103, 7720). The method utilizes MD simulations to explore the ground-state configuration space of the system and QM evaluation of those structures to yield the time-dependent electronic transition energy, which is transformed into the optical line-broadening function using the second-order cumulant expansion. Both linear and nonlinear optical spectra can then be generated for comparison to experiment. The dependence of the resulting spectra on the length of the MD trajectory, the QM sampling rate, and the QM model chemistry have all been examined. In particular, for the system of oxazine-4 in methanol studied here, at least 20 ps of MD trajectory are needed for qualitative convergence of linear spectral properties, and >100 ps is needed for quantitative convergence. Surprisingly, little difference is found between the 3-21G and 6-31G(d) basis sets, and the CIS and TD-B3LYP methods yield remarkably similar spectra. The semiempirical INDO/s method yields the most accurate results, reproducing the experimental Stokes shift to within 5% and the FWHM to within 20%. Nonlinear 3-pulse photon echo peak shift (3PEPS) decays have also been simulated. Decays are generally poorly reproduced, though the initial peak shift which depends on the overall coupling of motions to the solute transition energy is within 15% of experiment for all model chemistries other than those using the STO-3G basis.     

Trapping of a dopaquinone intermediate in the TPQ cofactor biogenesis in a copper-containing amine oxidase from Arthrobacter globiformis

The biogenesis of the topaquinone (TPQ) cofactor of copper amine oxidase (CAO) is self-catalyzed and requires copper and molecular oxygen. A dopaquinone intermediate has been proposed to undergo 1,4-addition of a copper-associated water molecule to form the reduced form of TPQ (TPQ(red)), followed by facile oxidation by O(2) to yield the mature TPQ (TPQ(ox)). In this study, we have incorporated a lysine residue in the active site of Arthrobacter globiformis CAO (AGAO) by site-directed mutagenesis to produce D298K-AGAO. The X-ray crystal structure of D298K-AGAO at 1.7-A resolution revealed that a covalent linkage formed between the epsilon-amino side chain of Lys298 and the C2 position of a dopaquinone derived from Tyr382, a precursor to TPQ(ox). We assigned the species as an iminoquinone tautomer (LTI) of lysine tyrosylquinone (LTQ), the organic cofactor of lysyl oxidase (LOX). The time course of the formation of LTI at pH 6.8 was followed by UV/vis and resonance Raman spectroscopies. In the early phase of the reaction, an LTQ-like intermediate was observed. This intermediate then slowly converted to LTI in an isosbestic manner. Not only is the presence of a dopaquinone intermediate in the TPQ biogenesis confirmed, but it also provides strong support for the proposed intermediacy of a dopaquinone in the biogenesis of LTQ in LOX. Further, this study indicates that the dopaquinone intermediate in AGAO is mobile and can swing from the copper site into the active-site wedge to react with Lys298.     

On the mechanism of reaction of radicals with tirapazamine

Ketyl radicals produced by photolysis of ketones or di-tert-butyl peroxide (DTBP) in alcohol solvents react rapidly with tirapazamine (TPZ). The acetone ketyl radical (ACOH) reacts with TPZ with an absolute second-order rate constant of (9.7 +/- 0.4) x 108 M-1 s-1. The reaction kinetics can be followed by monitoring the bleaching of TPZ absorption at 475 nm or the formation of a reaction product which absorbs at 320 and 410 nm. The ACOD radical reacts with TPZ in 2-propanol-OD with an absolute rate constant of (6.7 +/- 0.5) x 108 M-1 s-1, corresponding to a kinetic isotope effect (KIE) of 1.4. Deuteration of the radical on carbon (ACOH-d6) retards the reaction of the radical with TPZ even further (absolute rate constant = (4.8 +/- 0.04) x 108 M-1 s-1). This result corresponds to a KIE of 2.0. Radicals derived from dioxane and diisopropyl ether by flash photolysis of DTBP in ethereal solvent react with TPZ more slowly than do ketyl radicals. It is concluded that ketyl radicals react, in part, with TPZ in organic solvents by transfer of a hydrogen atom from the OH and CH3 groups of the ketyl radical to the oxygen atom at the N4 position of TPZ to form acetone or acetone enol and a radical derivative of TPZ (TPZH). The latter species absorbs at 320 and 405 nm, has a lifetime of hundreds of microseconds in alcohol solvents, and decays by disproportionation to form TPZ and a reduced heterocycle. The reduced heterocycle eventually forms a desoxytirapazamine by a polar mechanism. The results are supported by density functional theory calculations. It is proposed that dioxanyl radical will also react, in part, with TPZ by transfer of a hydrogen atom from the carbon adjacent to the radical center to the oxygen atom at the N4 position of TPZ. This produces the enol ether and the previously mentioned TPZH radical. It is further posited that ether radicals react a bit more slowly than ketyl radicals because they lack the second mode of hydrogen transfer (from the OH group) that is present in the ACOH radical. Our data are permissive of the possibility that ether radicals add to TPZ at a rate that is competitive with beta-hydrogen atom transfer.     

Redox-state-dependent complex formation between pseudoazurin and nitrite reductase

Bacterial copper-containing nitrite reductase catalyzes the reduction of nitrite to nitric oxide as part of the denitrification process. Pseudoazurin interacts with nitrite reductase in a transient fashion to supply the necessary electrons. The redox-state dependence of complex formation between pseudoazurin and nitrite reductase was studied by nuclear magnetic resonance spectroscopy and isothermal titration calorimetry. Binding of pseudoazurin in the reduced state is characterized by the presence of two binding modes, a slow and a fast exchange mode, with a K(d)(app) of 100 microM. In the oxidized state of pseudoazurin, binding occurs in a single fast exchange mode with a similar affinity. Metal-substituted proteins have been used to show that the mode of binding of pseudoazurin is independent of the metal charge of nitrite reductase. Contrary to what was found for other cupredoxins, protonation of the exposed His ligand to the copper of pseudoazurin, His81, does not appear to be involved directly in the dual binding mode of the reduced form. A model assuming the presence of a minor form of pseudoazurin is proposed to explain the behavior of the complex in the reduced state.     

A stabilized beta-oxaphosphoniumbetaine: an elusive zwitterion

A beta-oxaphosphoniumbetaine stabilized by two tert-butyl groups at the phosphonium part of the betaine and two trifluoromethyl groups at the carbon adjoined to the oxa part of the betaine was isolated and structurally characterized. Additional stabilization results from the solvation of the betaine by a parent phosphonium salt, as in 4, or by protonation with methanol, as in 10. According to X-ray analysis, the betaine exhibits a sterically strained "gauche" conformeric form, with torsion angles in the P-C-C-O moiety of 32.2 degrees for 4 and 28.1 degrees for 10. The P...O separations of 3.121 A for 4 and 3.086 A for 10 are just under the sum of the O and P van der Waals radii (3.32 A). The fast H/D exchange of alpha hydrogen atoms to the phosphonium center was observed in the solution of 10 in deuterium oxide and MeOH-d4. The beta-oxaphosphoniumbetaine has long been suspected as an intermediate in the Wittig reaction, but this is the first time a stabilized derivative has been isolated.     

(+)-saxitoxin: a first and second generation stereoselective synthesis

A stereoselective synthesis of the bis-guanidinium toxin (+)-saxitoxin (STX), the agent infamously associated with red tides and paralytic shellfish poisoning, is described. Our approach to this unique natural product advances through an unusual nine-membered ring guanidine intermediate 39 en route to the tricyclic skeleton that defines STX. The effectiveness of this strategy is notable, as only four steps are needed to transform 39 into the target molecule, including a four-electron alkene oxidation catalyzed by OsCl3. Construction of the critical monocyclic guanidine has been achieved through two channels, the first of which makes use of Rh-catalyzed C-H amination and highlights a novel class of heterocyclic N,O-acetals as iminium ion equivalents for crafting functionalized amines. A second route to 39 relies on a stereoselective acetylide dianion addition to a serine-based nitrone, thereby facilitating the preparation of STX in just 14 linear steps from commercial material.