Igor V. Kurnikov

Structural Requirements for Optimized Delivery, Inhibition of Oxidative Stress, and Antiapoptotic Activity of Targeted Nitroxides

Suppression of mitochondrial production of reactive oxygen species is a promising strategy against intrinsic apoptosis typical of degenerative diseases. Stable nitroxide radicals such as 4-hydroxy-2,2,6,6-tetramethyl piperidine-1-oxyl (TEMPOL) and its analogs combine several important features, including recycleability, electron acceptance from respiratory complexes, superoxide dismutase mimicry, and radical scavenging. Although successful in antioxidant protection, their effective concentrations are too high for successful in vivo applications. Recently (J Am Chem Soc 127:12460, 2005), we reported that 4-amino 2,2,6,6-tetramethyl-1-piperidinyloxy, covalently conjugated to a five-residue segment of gramicidin S (GS), was integrated into mitochondria and blocked actinomycin D (ActD)-induced superoxide generation and apoptosis. Using a model of ActD-induced apoptosis in mouse embryonic cells, we screened a library of nitroxides to explore structure-activity relationships between their antioxidant/antiapoptotic properties and chemical composition and three-dimensional (3D) structure. High hydrophobicity and effective mitochondrial integration are necessary but not sufficient for high antiapoptotic/antioxidant activity of a nitroxide conjugate. By designing conformationally preorganized peptidyl nitroxide conjugates and characterizing their 3D structure experimentally (circular dichroism and NMR) and theoretically (molecular dynamics), we established that the presence of the β-turn/β-sheet secondary structure is essential for the desired activity. Monte Carlo simulations in model lipid membranes confirmed that the conservation of the d-Phe-Pro reverse turn in hemi-GS analogs ensures the specific positioning of the nitroxide moiety at the mitochondrial membrane interface and maximizes their protective effects. These new insights into the structure-activity relationships of nitroxide-peptide and -peptide isostere conjugates are instrumental for development of new mechanism-based therapeutically effective agents.

Nitric Oxide Inhibits Peroxidase Activity of Cytochrome c· Cardiolipin Complex and Blocks Cardiolipin Oxidation

The increased production of NO during the early stages of apoptosis indicates its potential involvement in the regulation of programmed cell death through yet to be identified mechanisms. Recently, an important role for catalytically competent peroxidase form of pentacoordinate cytochrome c (cyt c) in a complex with a mitochondria-specific phospholipid, cardiolipin (CL), has been demonstrated during execution of the apoptotic program. Because the cyt c·CL complex acts as CL oxygenase and selectively oxidizes CL in apoptotic cells in a reaction dependent on the generation of protein-derived (tyrosyl) radicals, we hypothesized that binding and nitrosylation of cyt c regulates CL oxidation. Here we demonstrate by low temperature electron paramagnetic resonance spectroscopy that CL facilitated interactions of ferro- and ferri-states of cyt c with NO and NO–, respectively, to yield a mixture of penta- and hexa-coordinate nitrosylated cyt c. In the nitrosylated cyt c·CL complex, NO chemically reacted with H2O2-activated peroxidase intermediates resulting in their reduction. A dose-dependent quenching of H2O2-induced protein-derived radicals by NO donors was shown using direct electron paramagnetic resonance measurements as well as immuno-spin trapping with antibodies against protein 5,5-dimethyl-1-pyrroline N-oxide-nitrone adducts. In the presence of NO donors, H2O2-induced oligomeric forms of cyt c positively stained for 3-nitrotyrosine confirming the reactivity of NO toward tyrosyl radicals of cyt c. Interaction of NO with the cyt c·CL complex inhibited its peroxidase activity with three different substrates: CL, etoposide, and 3,3′-diaminobenzidine. Given the importance of CL oxidation in apoptosis, mass spectrometry analysis was utilized to assess the effects of NO on oxidation of 1,1′2,2′-tertalinoleoyl cardiolipin. NO effectively inhibited 1,1′2,2′-tertalinoleoyl cardiolipin oxidation catalyzed by the peroxidase activity of cyt c. Thus, NO can act as a regulator of peroxidase activity of cyt c·CL complexes.

Donor-bridge-acceptor energetics determine the distance dependence of electron tunneling in DNA

Electron transfer (ET) processes in DNA are of current interest because of their involvement in oxidative strand cleavage reactions and their relevance to the development of molecular electronics. Two mechanisms have been identified for ET in DNA, a single-step tunneling process and a multistep charge-hopping process. The dynamics of tunneling reactions depend on both the distance between the electron donor and acceptor and the nature of the molecular bridge separating the donor and acceptor. In the case of protein and alkane bridges, the distance dependence is not strongly dependent on the properties of the donor and acceptor. In contrast, we show here that the distance decay of DNA ET rates varies markedly with the energetics of the donor and acceptor relative to the bridge. Specifically, we find that an increase in the energy of the bridge states by 0.25 eV (1 eV = 1.602 × 10−19 J) relative to the donor and acceptor energies for photochemical oxidation of nucleotides, without changing the reaction free energy, results in an increase in the characteristic exponential distance decay constant for the ET rates from 0.71 to 1.1 Å−1. These results show that, in the small tunneling energy gap regime of DNA ET, the distance dependence is not universal; it varies strongly with the tunneling energy gap. These DNA ET reactions fill a “missing link” or transition regime between the large barrier (rapidly decaying) tunneling regime and the (slowly decaying) hopping regime in the general theory of bridge-mediated ET processes.

Ab initio based effective Hamiltonians for long‐range electron transfer: Hartree–Fock analysis

An ab initio electronic structure method is developed to describe electron transfer in large systems. The method is based on a molecular fragment effective Hamiltonian approach. The strategy pieces together results of ab initio quantum chemistry calculations on overlapping molecular segments in order to build an effective Hamiltonian that describes the long‐range electronic interactions. This is accomplished by constructing fragment effective Hamiltonians that properly describe the electronic propagation characteristics of each fragment (computed at the ab initio Hartree–Fock level in an appropriate basis set). The fragment effective Hamiltonian is projected onto the valence orbital space of each fragment, and a relatively well‐localized set of effective interactions is obtained. Combining these projected fragment Hamiltonians allows the construction of a valence effective Hamiltonian for the entire system. We find that the fragment Hamiltonian matrices constructed in this way are transferable between don...

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 H2O2-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.

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.

Pathways, pathway tubes, pathway docking, and propagators in electron transfer proteins

The simplest views of long-range electron transfer utilize flat one-dimensional barrier tunneling models, neglecting structural details of the protein medium. The pathway model of protein electron transfer reintroduces structure by distinguishing between covalent bonds, hydrogen bonds, and van der Waals contacts. These three kinds of interactions in a tunneling pathway each have distinctive decay factors associated with them. The distribution and arrangement of these bonded and nonbonded contacts in a folded protein varies tremendously between structures, adding a richness to the tunneling problem that is absent in simpler views. We review the pathway model and the predictions that it makes for protein electron transfer rates in small proteins, docked proteins, and the photosynthetic reactions center. We also review the formulation of the protein electron transfer problem as an effective two-level system. New multi-pathway approaches and improved electronic Hamiltonians are described briefly as well.

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 toxins 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.

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.

On the importance of accounting for nuclear quantum effects in ab initio calibrated force fields in biological simulations

Significance In molecular modeling the motion of nuclei, especially hydrogen, cannot be described using the laws of classical mechanics. The importance of nuclear quantum effects has long been appreciated by the ab initio molecular dynamics and by the water simulation communities. However, the vast majority of simulations of biological systems performed at ambient conditions treat atomic motion classically. Even in the new-generation force fields parameterized from quantum mechanics these effects are thought to be minor compared with other inaccuracies at room temperature and pressure. We show that a force field in excellent agreement with quantum mechanical energies and forces will not produce acceptably inaccurate predictions at ambient conditions unless the nuclear motion and interaction are accounted for in the simulation. In many important processes in chemistry, physics, and biology the nuclear degrees of freedom cannot be described using the laws of classical mechanics. At the same time, the vast majority of molecular simulations that employ wide-coverage force fields treat atomic motion classically. In light of the increasing desire for and accelerated development of quantum mechanics (QM)-parameterized interaction models, we reexamine whether the classical treatment is sufficient for a simple but crucial chemical species: alkanes. We show that when using an interaction model or force field in excellent agreement with the “gold standard” QM data, even very basic simulated properties of liquid alkanes, such as densities and heats of vaporization, deviate significantly from experimental values. Inclusion of nuclear quantum effects via techniques that treat nuclear degrees of freedom using the laws of classical mechanics brings the simulated properties much closer to reality.