The Determination of Free Energy of Hydration of Water Ions from First Principles.

We model the autoionization of water by determining the free energy of hydration of the major intermediate species of water ions. We represent the smallest ions─the hydroxide ion OH-, the hydronium ion H3O+, and the Zundel ion H5O2+─by bonded models and the more extended ionic structures by strong nonbonded interactions (e.g., the Eigen H9O4+ = H3O+ + 3(H2O) and the Stoyanov H13O6+ = H5O2+ + 4(H2O)). Our models are faithful to the precise QM energies and their components to within 1% or less. Using the calculated free energies and atomization energies, we compute the pKa of pure water from first principles as a consistency check and arrive at a value within 1.3 log units of the experimental one. From these calculations, we conclude that the hydronium ion, and its hydrated state, the Eigen cation, are the dominant species in the water autoionization process.

Neural Network Corrections to Intermolecular Interaction Terms of a Molecular Force Field Capture Nuclear Quantum Effects in Calculations of Liquid Thermodynamic Properties.

We incorporate nuclear quantum effects (NQE) in condensed matter simulations by introducing short-range neural network (NN) corrections to the ab initio fitted molecular force field ARROW. Force field NN corrections are fitted to average interaction energies and forces of molecular dimers, which are simulated using the Path Integral Molecular Dynamics (PIMD) technique with restrained centroid positions. The NN-corrected force field allows reproduction of the NQE for computed liquid water and methane properties such as density, radial distribution function (RDF), heat of evaporation (HVAP), and solvation free energy. Accounting for NQE through molecular force field corrections circumvents the need for explicit computationally expensive PIMD simulations in accurate calculations of the properties of chemical and biological systems. The accuracy and locality of pairwise NN NQE corrections indicate that this approach could be applicable to complex heterogeneous systems, such as proteins.

Combining Force Fields and Neural Networks for an Accurate Representation of Bonded Interactions.

We present a formalism of a neural network encoding bonded interactions in molecules. This intramolecular encoding is consistent with the models of intermolecular interactions previously designed by this group. Variants of the encoding fed into a corresponding neural network may be used to economically improve the representation of torsional degrees of freedom in any force field. We test the accuracy of the reproduction of the ab initio potential energy surface on a set of conformations of two dipeptides, methyl-capped ALA and ASP, in several scenarios. The encoding, either alone or in conjunction with an analytical potential, improves agreement with ab initio energies that are on par with those of other neural network-based potentials. Using the encoding and neural nets in tandem with an analytical model places the agreements firmly within "chemical accuracy" of ±0.5 kcal/mol.

Combining Force Fields and Neural Networks for an Accurate Representation of Chemically Diverse Molecular Interactions.

A key goal of molecular modeling is the accurate reproduction of the true quantum mechanical potential energy of arbitrary molecular ensembles with a tractable classical approximation. The challenges are that analytical expressions found in general purpose force fields struggle to faithfully represent the intermolecular quantum potential energy surface at close distances and in strong interaction regimes; that the more accurate neural network approximations do not capture crucial physics concepts, e.g., nonadditive inductive contributions and application of electric fields; and that the ultra-accurate narrowly targeted models have difficulty generalizing to the entire chemical space. We therefore designed a hybrid wide-coverage intermolecular interaction model consisting of an analytically polarizable force field combined with a short-range neural network correction for the total intermolecular interaction energy. Here, we describe the methodology and apply the model to accurately determine the properties of water, the free energy of solvation of neutral and charged molecules, and the binding free energy of ligands to proteins. The correction is subtyped for distinct chemical species to match the underlying force field, to segment and reduce the amount of quantum training data, and to increase accuracy and computational speed. For the systems considered, the hybrid ab initio parametrized Hamiltonian reproduces the two-body dimer quantum mechanics (QM) energies to within 0.03 kcal/mol and the nonadditive many-molecule contributions to within 2%. Simulations of molecular systems using this interaction model run at speeds of several nanoseconds per day.

Harnessing the Substrate Promiscuity of Dioxygenase AsqJ and Developing Efficient Chemoenzymatic Synthesis for Quinolones.

Nature has developed complexity-generating reactions within natural product biosynthetic pathways. However, direct utilization of these pathways to prepare compound libraries remains challenging due to limited substrate scopes, involvement of multiple-step reactions, and moderate robustness of these sophisticated enzymatic transformations. Synthetic chemistry, on the other hand, offers an alternative approach to prepare natural product analogs. However, owing to complex and diverse functional groups appended on the targeted molecules, dedicated design and development of synthetic strategies are typically required. Herein, by leveraging the power of chemo-enzymatic synthesis, we report an approach to bridge the gap between biological and synthetic strategies in the preparation of quinolone alkaloid analogs. Leading by in silico analysis, the predicted substrate analogs were chemically synthesized. The AsqJ-catalyzed asymmetric epoxidation of these substrate analogues was followed by an Lewis Acid-triggered ring contraction to complete the viridicatin formation. We evaluated the robustness of this method in gram-scale reactions. Lastly, through chemoenzymatic cascades, a library of quinolone alkaloids is effectively prepared.

Epoxidation Catalyzed by the Nonheme Iron(II)- and 2-Oxoglutarate-Dependent Oxygenase, AsqJ: Mechanistic Elucidation of Oxygen Atom Transfer by a Ferryl Intermediate.

Mechanisms of enzymatic epoxidation via oxygen atom transfer (OAT) to an olefin moiety is mainly derived from the studies on thiolate-heme containing epoxidases, such as cytochrome P450 epoxidases. The molecular basis of epoxidation catalyzed by nonheme-iron enzymes is much less explored. Herein, we present a detailed study on epoxidation catalyzed by the nonheme iron(II)- and 2-oxoglutarate-dependent (Fe/2OG) oxygenase, AsqJ. The native substrate and analogues with different para substituents ranging from electron-donating groups (e.g., methoxy) to electron-withdrawing groups (e.g., trifluoromethyl) were used to probe the mechanism. The results derived from transient-state enzyme kinetics, Mössbauer spectroscopy, reaction product analysis, X-ray crystallography, density functional theory calculations, and molecular dynamic simulations collectively revealed the following mechanistic insights: (1) The rapid O2 addition to the AsqJ Fe(II) center occurs with the iron-bound 2OG adopting an online-binding mode in which the C1 carboxylate group of 2OG is trans to the proximal histidine (His134) of the 2-His-1-carboxylate facial triad, instead of assuming the offline-binding mode with the C1 carboxylate group trans to the distal histidine (His211); (2) The decay rate constant of the ferryl intermediate is not strongly affected by the nature of the para substituents of the substrate during the OAT step, a reactivity behavior that is drastically different from nonheme Fe(IV)-oxo synthetic model complexes; (3) The OAT step most likely proceeds through a stepwise process with the initial formation of a C(benzylic)-O bond to generate an Fe-alkoxide species, which is observed in the AsqJ crystal structure. The subsequent C3-O bond formation completes the epoxide installation.

Modeling Electronic Polarizability Changes in the Course of a Magnesium Ion Water Ligand Exchange Process.

This paper introduces explicit dependence of atomic polarizabilities on intermolecular interactions within the framework of a polarizable force field AMOEBA. Polarizable models used in biomolecular simulations often poorly describe molecular electrostatic induction in condensed phase, in part, due to neglect of a strong dependency of molecular electronic polarizability on intermolecular interactions at short distances. Our variable polarizability model parameters are derived from quantum chemical calculations of small clusters of atoms and molecules, and can be applied in simulations in condensed phase without additional scaling factors. The variable polarizability model is applied to simulate a ligand exchange reaction for a Mg(2+) ion solvated in water. Explicit dependence of water polarizability on a distance between a water oxygen and Mg(2+) is derived from in vacuum MP2 calculations of Mg(2+)-water dimer. The simulations yield a consistent description of the energetics of the Mg(2+)-water clusters of different size. Simulations also reproduce thermodynamics of ion solvation as well as kinetics of a water ligand exchange reaction. In contrast, simulations that used the additive force field or that used the constant polarizability models were not able to consistently and quantitatively describe the properties of the solvated Mg(2+) ion.

Donor/acceptor coupling shortcuts in electron transfer within ruthenium-modified derivatives of cytochrome b(562).

Quantitative theoretical studies of long-range electron transfer are still rare, and reliable computational methods to analyze these reactions are still being developed. We re-examined electron transfer reactions in ruthenium-modified cytochrome b562 derivatives focusing on accurate calculation of statistical average of electron transfer rates that are dominated by a small fraction of accessible protein conformations. We performed a series of ab initio calculations of donor/acceptor interactions over protein fragments sampled from long molecular dynamic trajectories and compared computed electron transfer rates to available experimental data. Our approach takes into account cofactor electronic structure and effects of solvation on the donor-acceptor interactions. It allows predicting absolute values of electron transfer rates in contrast to other computational methodologies that give only qualitative results. Our calculations reproduced with a good accuracy experimental electron transfer rates. We also found that electron transfer in some of the cytochrome b562 derivatives is dominated by "shortcut" conformations, where donor/acceptor interactions are mediated by nonbonded interactions of Ru ligands with protein surface groups. Several derivatives adopt long-lived conformations with the Ru complex interacting with negatively charged protein residues that are characterized by shorter Ru-Fe distances and higher ET rates. We argue that quantitative theoretical analysis is essential for detailed understanding of protein electron transfer and mechanisms of biological redox reactions.

A glance at the reactivity of osma(II)cycles [Os(C-N)x(bpy)3-x](m+) (x=0-3) Covering a 1.8 V potential range toward peroxidase through Monte Carlo Simulations ((-)C-N=o-2-phenylpyridinato, bpy=2,2'-bipyridine).

Three cyclometalated and one coordination compounds [Os(C-N)x(bpy)3-x](m) (x/m=0/2+ (4); 1/1+ (3); 2/1+ (2); 3/0 (1); (-)C-N=2-phenylpyridinato, bpy=2,2'-bipyridine) with drastically different reduction potentials have been used for analyzing the second-order rate constants for one-electron, metal-based osmium(II) to osmium(III) oxidation of the complexes by compound I (k2) and compound II (k3) of horseradish peroxidase. Previously unknown k2 and k3 have been determined by digital simulation of cyclic voltammograms measured in phosphate buffer of pH7.6 and 21 ± 1°C. Osmium(II) species derived from osmium(III) complexes 1 and 2 were generated electrochemically in situ. Under the conditions used the reduction potentials for the Os(III/II) feature equal -0.90, -0.095, 0.23 and 0.85V versus NHE (normal hydrogen electrode) for 1-4, respectively. The rate constants k2 equal ~5 × 10(7), 6 × 10(8), 2 × 10(6) and 1 × 10(5)M(-1)s(-1) and the rate constants k3 equal ~9 × 10(6), 4× 10(7), 1 ×10(6) and 1 × 10(5)M(-1)s(-1) for complexes 1-4, respectively. Both rate constants k2 and k3 first increase with increasing the reaction driving force on going from 4 to 2 but then both decrease on going to complex 1 though the reaction driving force is the highest in this case. The system described has been explored theoretically using docking Monte Carlo simulations.

Control of cytochrome c redox reactivity through off-pathway modifications in the protein hydrogen-bonding network.

Measurements of photoinduced Fe(2+)-to-Ru(3+) electron transfer (ET), supported by theoretical analysis, demonstrate that mutations off the dominant ET pathways can strongly influence the redox reactivity of cytochrome c. The effects arise from the change in the protein dynamics mediated by the intraprotein hydrogen-bonding network.

pH-triggered conformational switching of the diphtheria toxin T-domain: the roles of N-terminal histidines.

pH-induced conformational switching is essential for functioning of diphtheria toxin, which undergoes a membrane insertion/translocation transition triggered by endosomal acidification as a key step of cellular entry. In order to establish the sequence of molecular rearrangements and side-chain protonation accompanying the formation of the membrane-competent state of the toxin's translocation (T) domain, we have developed and applied an integrated approach that combines multiple techniques of computational chemistry [e.g., long-microsecond-range, all-atom molecular dynamics (MD) simulations; continuum electrostatics calculations; and thermodynamic integration (TI)] with several experimental techniques of fluorescence spectroscopy. TI calculations indicate that protonation of H257 causes the greatest destabilization of the native structure (6.9 kcal/mol), which is consistent with our early mutagenesis results. Extensive equilibrium MD simulations with a combined length of over 8 µs demonstrate that histidine protonation, while not accompanied by the loss of structural compactness of the T-domain, nevertheless results in substantial molecular rearrangements characterized by the partial loss of secondary structure due to unfolding of helices TH1 and TH2 and the loss of close contact between the C- and N-terminal segments. The structural changes accompanying the formation of the membrane-competent state ensure an easier exposure of the internal hydrophobic hairpin formed by helices TH8 and TH9, in preparation for its subsequent transmembrane insertion.

Electron tunneling through sensitizer wires bound to proteins.

We report a quantitative theoretical analysis of long-range electron transfer through sensitizer wires bound in the active-site channel of cytochrome P450cam. Each sensitizer wire consists of a substrate group with high binding affinity for the enzyme active site connected to a ruthenium-diimine through a bridging aliphatic or aromatic chain. Experiments have revealed a dramatic dependence of electron transfer rates on the chemical composition of both the bridging group and the substrate. Using combined molecular dynamics simulations and electronic coupling calculations, we show that electron tunneling through perfluorinated aromatic bridges is promoted by enhanced superexchange coupling through virtual reduced states. In contrast, electron flow through aliphatic bridges occurs by hole-mediated superexchange. We have found that a small number of wire conformations with strong donor-acceptor couplings can account for the observed electron tunneling rates for sensitizer wires terminated with either ethylbenzene or adamantane. In these instances, the rate is dependent not only on electronic coupling of the donor and acceptor but also on the nuclear motion of the sensitizer wire, necessitating the calculation of average rates over the course of a molecular dynamics simulation. These calculations along with related recent findings have made it possible to analyze the results of many other sensitizer-wire experiments that in turn point to new directions in our attempts to observe reactive intermediates in the catalytic cycles of P450 and other heme enzymes.

Kinetic and theoretical comprehension of diverse rate laws and reactivity gaps in Coriolus hirsutus laccase-catalyzed oxidation of acido and cyclometalated Ru(II) complexes.

The reactivity of the acido Ru(II) complexes cis-[RuCl(2)(LL)(2)], [RuCO(3)(LL)(2)], cis-[RuCO(3)-(bquin)(2)] (LL = 2,2'-bipyridine (bpy) and 1,10-phenanthroline (phen); bquin = 2,2'-biquinoline) and cyclometalated Ru(II) derivatives of 2-phenylpyridine and 4-(2-tolyl)pyridine [Ru(o-C(6)H(4)-2-py)(phen)(2)]PF(6) (1), [Ru(o-C(6)H(3)-p-R-2-py)(bpy)(MeCN)(2)]PF(6) (2), and [Ru(o-C(6)H(3)-p-R-2-py)(phen)(MeCN)(2)]PF(6) (3) (R = H (a), Me (b)) toward laccase from Coriolus hirsutus has been investigated by conventional UV-vis spectroscopy at pH 3-7 and 25 degrees C. The acido and cyclometalated complexes are readily oxidized into the corresponding Ru(III) species, but the two types of complexes differ substantially in reactivity and obey different rate laws. The acido complexes are oxidized more slowly and the second-order kinetics, first-order in laccase and Ru(II), holds with the rate constants around 5 x 10(4) M(-1) s(-1) at pH 4.5 and 25 degrees C. The cyclometalated complexes 1-3 react much faster and the hyperbolic Michaelis-Menten kinetics holds. However, it is not due to formation of an enzyme-substrate complex but rather because of the ping-pong mechanism of catalysis, viz. E(ox) + Ru(II) --> E(red) + Ru(III) (k(1)); E(red) + 1/4O(2) --> E(ox) (k(2)), with the rate constants k(1) in the range (2-9) x 10(7) M(-1) s(-1) under the same conditions. The huge values of k(1) move the enzymatic oxidation toward a kinetic regime when the dioxygen half-reaction becomes the rate-limiting step. Cyclometalated compounds 1-3 can therefore be used for routine estimation of k(2), that is, the rate constant for reoxidation for laccases by dioxygen. The mechanism proposed was confirmed by the direct stopped-flow measurements of the k(2) rate constant (8.1 x 10(5) M(-1) s(-1) at 26 degrees C) and supported by the theoretical modeling of interaction between the bpy analogue of 1 and Coriolus hirsutes laccase using Monte Carlo simulations.

Cardiolipin deficiency leads to decreased cardiolipin peroxidation and increased resistance of cells to apoptosis.

Cardiolipin (CL), a unique mitochondrial phospholipid synthesized by CL synthase (CLS), plays important, yet not fully understood, roles in mitochondria-dependent apoptosis. We manipulated CL levels in HeLa cells by knocking down CLS using RNA interference and selected a clone of CL-deficient cells with approximately 45% of its normal content. ESI-MS analysis showed that the CL molecular species were the same in CL-deficient and CL-sufficient cells. CL deficiency did not change mitochondrial functions (membrane potential, reactive oxygen species generation, cellular ATP levels) but conferred resistance to apoptosis induced by actinomycin D (ActD), rotenone, or gamma-irradiation. During ActD-induced apoptosis, decreased CL peroxidation along with suppressed cytochrome (cyt) c release was observed in CL-deficient cells, whereas Bax translocation to mitochondria remained similar to that in CL-sufficient HeLa cells. The amounts of loosely bound cyt c (releasable under high ionic strength conditions) were the same in CL-deficient and CL-sufficient cells. Given that CL peroxidation during apoptosis is catalyzed by CL/cyt c complexes and CL oxidation products are essential for cyt c release from mitochondria, our results suggest that CL deficiency prevents adequate assembly of productive CL/cyt c complexes and CL peroxidation, resulting in increased resistance to apoptosis.

The hierarchy of structural transitions induced in cytochrome c by anionic phospholipids determines its peroxidase activation and selective peroxidation during apoptosis in cells.

Activation of peroxidase catalytic function of cytochrome c (cyt c) by anionic lipids is associated with destabilization of its tertiary structure. We studied effects of several anionic phospholipids on the protein structure by monitoring (1) Trp59 fluorescence, (2) Fe-S(Met80) absorbance at 695 nm, and (3) EPR of heme nitrosylation. Peroxidase activity was probed using several substrates and protein-derived radicals. Peroxidase activation of cyt c did not require complete protein unfolding or breakage of the Fe-S(Met80) bond. The activation energy of cyt c peroxidase changed in parallel with stability energies of structural regions of the protein probed spectroscopically. Cardiolipin (CL) and phosphatidic acid (PA) were most effective in inducing cyt c peroxidase activity. Phosphatidylserine (PS) and phosphatidylinositol bisphosphate (PIP2) displayed a significant but much weaker capacity to destabilize the protein and induce peroxidase activity. Phosphatidylinositol trisphosphate (PIP3) appeared to be a stronger inducer of cyt c structural changes than PIP2, indicating a role for the negatively charged extra phosphate group. Comparison of cyt c-deficient HeLa cells and mouse embryonic cells with those expressing a full complement of cyt c demonstrated the involvement of cyt c peroxidase activity in selective catalysis of peroxidation of CL, PS, and PI, which corresponded to the potency of these lipids in inducing cyt c's structural destabilization.

Cardiolipin switch in mitochondria: shutting off the reduction of cytochrome c and turning on the peroxidase activity.

Upon interaction with anionic phospholipids, particularly mitochondria-specific cardiolipin (CL), cytochrome c (cyt c) loses its tertiary structure and its peroxidase activity dramatically increases. CL-induced peroxidase activity of cyt c has been found to be important for selective CL oxidation in cells undergoing programmed death. During apoptosis, the peroxidase activity and the fraction of CL-bound cyt c markedly increase, suggesting that CL may act as a switch to regulate cyt c's mitochondrial functions. Using cyclic voltammetry and equilibrium redox titrations, we show that the redox potential of cyt c shifts negatively by 350-400 mV upon binding to CL-containing membranes. Consequently, functions of cyt c as an electron transporter and cyt c reduction by Complex III are strongly inhibited. Further, CL/cyt c complexes are not effective in scavenging superoxide anions and are not effectively reduced by ascorbate. Thus, both redox properties and functions of cyt c change upon interaction with CL in the mitochondrial membrane, diminishing cyt c's electron donor/acceptor role and stimulating its peroxidase activity.

Coupling coherence distinguishes structure sensitivity in protein electron transfer.

Quantum mechanical analysis of electron tunneling in nine thermally fluctuating cytochrome b562 derivatives reveals two distinct protein-mediated coupling limits. A structure-insensitive regime arises for redox partners coupled through dynamically averaged multiple-coupling pathways (in seven of the nine derivatives) where heme-edge coupling leads to the multiple-pathway regime. A structure-dependent limit governs redox partners coupled through a dominant pathway (in two of the nine derivatives) where axial-ligand coupling generates the single-pathway limit and slower rates. This two-regime paradigm provides a unified description of electron transfer rates in 26 ruthenium-modified heme and blue-copper proteins, as well as in numerous photosynthetic proteins.

Mechanisms of cardiolipin oxidation by cytochrome c: relevance to pro- and antiapoptotic functions of etoposide.

Execution of apoptotic program in mitochondria is associated with accumulation of cardiolipin peroxidation products required for the release of proapoptotic factors into the cytosol. This suggests that lipid antioxidants capable of inhibiting cardiolipin peroxidation may act as antiapoptotic agents. Etoposide, a widely used antitumor drug and a topoisomerase II inhibitor, is a prototypical inducer of apoptosis and, at the same time, an effective lipid radical scavenger and lipid antioxidant. Here, we demonstrate that cardiolipin oxidation during apoptosis is realized not via a random cardiolipin peroxidation mechanism but rather proceeds as a result of peroxidase reaction in a tight cytochrome c/cardiolipin complex that restrains interactions of etoposide with radical intermediates generated in the course of the reaction. Using low-temperature and ambient-temperature electron paramagnetic resonance spectroscopy of H(2)O(2)-induced protein-derived (tyrosyl) radicals and etoposide phenoxyl radicals, respectively, we established that cardiolipin peroxidation and etoposide oxidation by cytochrome c/cardiolipin complex takes place predominantly on protein-derived radicals of cytochrome c. We further show that etoposide can inhibit cytochrome c-catalyzed oxidation of cardiolipin competing with it as a peroxidase substrate. Peroxidase reaction of cytochrome c/cardiolipin complexes causes cross-linking and oligomerization of cytochrome c. With nonoxidizable tetraoleoyl-cardiolipin, the cross-linking occurs via dityrosine formation, whereas bifunctional lipid oxidation products generated from tetralinoleoyl-cardiolipin participate in the production of high molecular weight protein aggregates. Protein aggregation is effectively inhibited by etoposide. The inhibition of cardiolipin peroxidation by etoposide, however, is realized at far higher concentrations than those at which it induces apoptotic cell death. Thus, oxidation of cardiolipin by the cytochrome c/cardiolipin peroxidase complex, which is essential for apoptosis, is not inhibited by proapoptotic concentrations of the drug.

Peroxidase activity and structural transitions of cytochrome c bound to cardiolipin-containing membranes.

During apoptosis, cytochrome c (cyt c) is released from intermembrane space of mitochondria into the cytosol where it triggers the caspase-dependent machinery. We discovered that cyt c plays another critical role in early apoptosis as a cardiolipin (CL)-specific oxygenase to produce CL hydroperoxides required for release of pro-apoptotic factors [Kagan, V. E., et al. (2005) Nat. Chem. Biol. 1, 223-232]. We quantitatively characterized the activation of peroxidase activity of cyt c by CL and hydrogen peroxide. At low ionic strength and high CL/cyt c ratios, peroxidase activity of the CL/cyt c complex was increased >50 times. This catalytic activity correlated with partial unfolding of cyt c monitored by Trp(59) fluorescence and absorbance at 695 nm (Fe-S(Met(80)) band). The peroxidase activity increase preceded the loss of protein tertiary structure. Monounsaturated tetraoleoyl-CL (TOCL) induced peroxidase activity and unfolding of cyt c more effectively than saturated tetramyristoyl-CL (TMCL). TOCL/cyt c complex was found more resistant to dissociation by high salt concentration. These findings suggest that electrostatic CL/cyt c interactions are central to the initiation of the peroxidase activity, while hydrophobic interactions are involved when cyt c's tertiary structure is lost. In the presence of CL, cyt c peroxidase activity is activated at lower H(2)O(2) concentrations than for isolated cyt c molecules. This suggests that redistribution of CL in the mitochondrial membranes combined with increased production of H(2)O(2) can switch on the peroxidase activity of cyt c and CL oxidation in mitochondria-a required step in execution of apoptosis.