Absorption spectra of riboflavin--a difficult case for computational chemistry

Computing accurate absorption spectra of riboflavin (RBF) has proven a difficult task for computational chemistry. Time-dependent density functional theory have herein been employed using a wide range of recent range-separated and hybrid meta functionals to investigate ultraviolet and visible spectra of RBF to determine if any progress has been made through recent developments. It is concluded that B3LYP and PBE0 perform the best throughout the entire test set. However, since all methods deviate from experimental results by at least 40 nm when computing the spectra in vacuum, two approaches to describe aqueous solution are implemented together with the MPWB1K, B3LYP, and PBE0 functionals: implicitly using integral equation formulation of the polarized continuum model (minor improvement) and explicitly through molecular dynamics (MD) simulations of the molecule embedded in a water cluster whereafter snapshots of RBF-water clusters are extracted and time-dependent density functional theory calculations performed. The resulting averaged spectra from the MD-simulated clusters show a constant blue-shift for all peaks by ∼20 nm compared to experimental data at the TD-B3LYP/6-31+G(d,p) level.

Towards echinomycin mimetics by grafting quinoxaline residues on glycophane scaffolds

Echinomycin is a natural depsipeptide, which is a bisintercalator, inserting quinoxaline units preferentially adjacent to CG base pairs of DNA. Herein the design and synthesis of echinomycin mimetics based on grafting of two quinoxaline residues onto a macrocyclic scaffold (glycophane) is addressed. Binding of the compounds to calf-thymus DNA was studied using UV-vis and steady state fluorescence spectroscopy, as well as thermal denaturation. An interesting observation was enhancement of fluorescence emission for the peptidomimetics on binding to DNA, which contrasted with observations for echinomycin. Molecular dynamics simulations were exploited to explore in more detail if bis-intercalation to DNA was possible for one of the glycophanes. Bis-intercalating echinomycin complexes with DNA were found to be stable during 20ns simulations at 298K. However, the MD simulations of a glycophane complexed with a DNA octamer displayed very different behaviour to echinomycin and its quinoxaline units were found to rapidly migrate out from the intercalation site. Release of bis-intercalation strain occurred with only one of the quinoxaline chromophores remaining intercalated throughout the simulation. The distance between the quinoxaline residues in the glycophane at the end of the MD simulation was 7.3-7.5Å, whereas in echinomycin, the distance between the residues was ∼11Å, suggesting that longer glycophane scaffolds would be required to generate bis-intercalating echinomycin mimetics.

Hydroxylation and ring-opening mechanism of an unusual flavoprotein monooxygenase, 2-methyl-3-hydroxypyridine-5-carboxylic acid oxygenase: a theoretical study

Hybrid meta-GGA density functional theory (the MPWB1K functional) was used to study the hydroxylation and ring-opening mechanism of 2-methyl-3-hydroxypyridine-5-carboxylic acid oxygenase (MHPCO). This enzyme catalyses the conversion of 2-methyl-3-hydroxypyridine-5-carboxylic acid (MHPC) to alpha-(N-acetylaminomethylene)succinic acid (AAMS), which is the essential ring-opening step in the bacterial degradation of vitamin B(6). MHPCO belongs to the flavin-containing aromatic hydroxylases family. However, MHPCO is capable of catalysing a subsequent aromatic ring-cleavage reaction to give acyclic products rather than hydroxylated aromatic ones. Our calculations show that the re-aromatisation of the hydroxylated intermediate occurs spontaneously in aqueous solution; this implies that the ring-opening process occurs inside the enzyme's active site, in which limited water is available. The instability of the hydroxylated intermediate of MHPCO is the main reason why acyclic products are formed. Previously proposed mechanisms for the ring-opening step were studied, and were shown to be less likely to occur (DeltaDeltaG(not equal298)>35 kcal mol(-1)). Two new pathways with reasonable barrier heights (DeltaDeltaG(not equal 298)<15 kcal mol(-1)) are reported herein, which are in accordance with all experimental information present to date.

Theoretical study on conformational preferences of ribose in 2-thiouridine--the role of the 2'OH group

Conformational changes in ribose are well-known to play a significant role in biomolecular identification. The mechanism of selectivity towards C3'-endo conformation (conformer b) in ribose of 2-thiouridine has been studied using DFT (B3LYP) and MP2 methodology, together with 6-31+G(d,p) basis set. The polarity of the C2S2 bond is enhanced due to the orientation of H2' towards the S2 atoms, which leads to a difference in the corresponding bond lengths, the atomic charges and the vO2'H2' stretch vibrations in all the conformers. NBO analysis shows that charge transfer mainly occurs in the C2N3 and C2S2 orbitals. The higher stability of conformer b is attributed to its larger orbital interaction energies within the 2-thiouracil base, and total orbital interaction energies of conformer b. Our conclusion is that the distant electrostatic rather than hydrogen bonding effects between 2'OH and the S2 atoms play the dominant role in the orbital interaction, and enhance the selectivity towards the C3'-endo conformation of ribose.

Evidence of high *OH radical quenching efficiency by vitamin B6

Molecules acting as antioxidants capable of scavenging reactive oxygen species (ROS) are of the utmost importance in the living cell. The antioxidative properties of pyridoxine (vitamin B6) have recently been discovered. Previous theoretical calculations have shown a high reactivity of pyridoxine toward hydroxyl radicals, where the latter preferably abstract H from either carbon of the two methanol substituents (C8 or C9). In this study, we have explored the reactivity of pyridoxine toward further hydroxyl radicals, considering as the first step the H abstraction from either C8 or C9, also including addition reactions and cyclization. Many of the reactions display similar DeltaG, and hence, the quenching of hydroxyl radicals by pyridoxine may undergo different pathways leading to a mix of products. In addition, we observe that pyridoxine, under high hydroxyl radical concentrations, may scavenge up to eight radicals, supporting its observed high antioxidant activity.

Distinct hydroxy-radical-induced damage of 3'-uridine monophosphate in RNA: a theoretical study

RNA strand scission and base release in 3'-uridine monophosphate (UMP), induced by OH radical addition to uracil, is studied at the DFT B3LYP/6-31+G(d,p) level in the gas phase and in solution. In particular, the mechanism of hydrogen-atom transfer subsequent to radical formation, from C2' on the sugar to the C6 site on the base, is explored. The barriers of (C2'-)H2'(a) abstraction by the C6 radical site range from 11.2 to 20.0 kcal mol(-1) in the gas phase and 14.1 to 21.0 kcal mol(-1) in aqueous solution, indicating that the local surrounding governs the hydrogen-abstraction reaction in a stereoselective way. The calculated N1-C1' (N1-glycosidic bond) and beta-phosphate bond strengths show that homolytic and heterolytic bond-breaking processes are largely favored in each case, respectively. The barrier for beta-phosphate bond rupture is approximately 3.2-4.0 kcal mol(-1) and is preferred by 8-12 kcal mol(-1) over N1-glycosidic bond cleavage in both the gas phase and solution. The beta-phosphate bond-rupture reactions are exothermal in the gas phase and solution, whereas N1-C1' bond-rupture reactions require both solvation and thermal corrections at 298 K to be energetically favored. The presence of the ribose 2'-OH group and its formation of low-barrier hydrogen bonds with oxygen atoms of the 3'-phosphate linkage are highly important for hydrogen transfer and the subsequent bond-breakage reactions.

Theoretical study of the phototoxicity of naproxen and the active form of nabumetone

Density functional theory using the hybrid functional B3LYP has been employed in order to study the mechanisms of photoinduced decomposition of the closely related nonsteroidal anti-inflammatory drugs naproxen (NP) and 6-methoxy-2-naphthylacetic acid (MNAA; the active form of nabumetone). The photochemical properties and computed energies of various species obtained in this study show that both drugs dominate in their deprotonated forms at physiological pH. The deprotonated acids are unable to decarboxylate from their excited singlets; instead, they decarboxylate from their first excited triplet states with high efficiency, overcoming energy barriers less than 3 and 1 kcal/mol for MNAA and NP, respectively. The ultraviolet and visible spectra of the neutral, deprotonated, and decarboxylated moieties of MNAA and NP are more-or-less similar but with higher probabilites (oscillator strength) for the latter. This fact, as well as the higher reactivity of NP, is explained in terms of the electron-donating effect of the additional methyl group present in NP. Singlet oxygen, superoxide radical anion, and corresponding peroxyl radical species are expected to be formed in different steps throughout the proposed photodegradation pathways of both drugs, which give rise to their effects on biomolecules, for example, lipid peroxidation.

Deamination of the radical cation of the base moiety of 2'-deoxycytidine: a theoretical study

Five pathways leading to the deamination of cytosine (to uracil) after formation of its deprotonated radical cation are investigated in the gas phase, at the UB3LYP/6-311G(d,p) level of theory, and in bulk aqueous solvent. The most favorable pathway involves hydrogen-atom transfer from a water molecule to the N3 nitrogen of the deprotonated radical cation, followed by addition of the resulting hydroxyl radical to the C4 carbon of the cytosine derivative. Following protonation of the amino group (N4), the C4--N4 bond is broken with elimination of the NH3+(. ) and formation of a protonated uracil. The rate-determining step of this mechanism is hydrogen-atom transfer from a water molecule to the cytosine derivative. The associated free energy barrier is 70.2 kJ mol(-1).

Theoretical studies of damage to 3'-uridine monophosphate induced by electron attachment

Low-energy electrons (LEE) are well known to induce nucleic acid damage. However, the damage mechanisms related to charge state and structural features remain to be explored in detail. In the present work, we have investigated the N1-glycosidic and C3'-O(P) bond ruptures of 3'-UMP (UMP=uridine monophosphate) and the protonated form 3'-UMPH with -1 and zero charge, respectively, based on hybrid density functional theory (DFT) B3 LYP together with the 6-31+G(d,p) basis set. The glycosidic bond breakage reactions of the 3'UMP and 3'UMPH electron adducts are exothermic in both cases, with barrier heights of 19-20 kcal mol(-1) upon inclusion of bulk solvation. The effects of the charge state on the phosphate group are marginal, but the C2'-OH group destabilizes the transition structure of glycosidic bond rupture of 3'-UMPH in the gas phase by approximately 5.0 kcal mol(-1). This is in contrast with the C3'-O(P) bond ruptures induced by LEE in which the charge state on the phosphate influences the barrier heights and reaction energies considerably. The barrier towards C3'-O(P) bond dissociation in the 3'UMP electron adduct is higher in the gas phase than the one corresponding to glycosidic bond rupture and is dramatically influenced by the C2'-OH group and bulk salvation, which decreases the barrier to 14.7 kcal mol(-1). For the C3'-O(P) bond rupture of the 3'UMPH electron adduct, the reaction is exothermic and the barrier is even lower, 8.2 kcal mol(-1), which is in agreement with recent results for 3'-dTMPH and 5'-dTMPH (dTMPH=deoxythymidine monophosphate). Both the Mulliken atomic charges and unpaired-spin distribution play significant roles in the reactions.

*H atom and *OH radical reactions with 5-methylcytosine

The reactions between either a hydrogen atom or a hydroxyl radical and 5-methylcytosine (5-MeCyt) are studied by using the hybrid kinetic energy meta-GGA functional MPW1B95. *H atom and *OH radical addition to positions C5 and C6 of 5-MeCyt, or *OH radical induced H-abstraction from the C5 methyl group, are explored. All systems are optimized in bulk solvent. The data presented show that the barriers to reaction are very low: ca. 7 kcal/mol for the *H atom additions and 1 kcal/mol for the reactions involving the *OH radical. Thermodynamically, the two C6 radical adducts and the *H-abstraction product are the most stable ones. The proton hyperfine coupling constants (HFCC), computed at the IEFPCM/MPW1B95/6-311++G(2d,2p) level, agree well with B3LYP results and available experimental and theoretical data on related thymine and cytosine radicals.

A triplet mechanism for the formation of cyclobutane pyrimidine dimers in UV-irradiated DNA

The reaction pathways for the photochemical formation of cyclobutane thymine dimers in DNA are explored using hybrid density functional theory techniques. It is concluded that the thymine-thymine [2 + 2] cycloaddition displays favorable energy barriers and reaction energies in both the triplet and the singlet excited states. The stepwise cycloaddition in the triplet excited state involves the initial formation of a diradical followed by ring closure via singlet-triplet interaction. The triplet mechanism is thus completely different from the concerted singlet state cycloaddition processes. The key geometric features and electron spin densities are also discussed. Bulk solvation has a major effect by reducing the barriers and increasing the diradical stabilities. The present results provide a rationale for the faster cycloreaction observed in the singlet excited states than in the triplet excited states.

Hydrolysis process of the second generation platinum-based anticancer drug cis-amminedichlorocyclohexylamineplatinum(II)

The hydrolysis process of the anticancer drug cis-amminedichlorocyclohexylamineplatinum(II) (JM118 or cis-[PtCl2(NH3)cyclohexylamine]) and the influence of solvent models therein have been studied using hybrid density functional theory (B3LYP). The aquation reactions leading to the activated drug forms a key step for the reaction with the target DNA. In this study, the stepwise hydrolysis, cis-[PtCl2(NH3)cyclohexylamine] + 2 H2O --> cis-[Pt(NH3)cyclohexylamine(OH2)2]2+ + 2 Cl- was explored, using three different models. Implicit solvent effects were incorporated through polarized continuum models. The stationary points on the potential energy surfaces for the first and second hydrolysis steps, proceeding via a general S(N)2 pathway, were fully optimized and characterized. It was found that the explicit solvent effects originating from the inclusion of extra water molecules into the system are significantly stronger than those arising from the bulk aqueous medium, especially for the second aquation step, emphasizing the use of appropriate models for these types of problems. In comparison with previous work on the parent compound cisplatin, a slower rate of hydrolysis is determined for the first (rate determining) reaction. The results furthermore imply that the doubly aquated form of JM118 will be the main DNA binding form of the drug. The results provide detailed energy profiles for the mechanism of hydrolysis of JM118, which may assist in understanding the reaction mechanism of the drug with the DNA target and in the design of novel Pt-containing anticancer drugs.

Metal Ion Dependent Adhesion Sites in Integrins: A Combined DFT and QMC Study on Mn2+

The theoretical study of relative energies of different spin states of Mn2+ has been carried out for the isolated cation and for structures in which the cation is coordinated to ligands that represent the first coordination shell in a protein environment that contains a metal ion dependent adhesion site (MIDAS, found in the ligand binding domain of protein LFA-1). The calculations determine whether the ligand field generated by a prototype protein environment affects the relative energies between high, intermediate, and low spin states. Geometry optimizations and vibrational frequency calculations were carried out at the B3LYP/SKBJ+* level of theory. Single point calculations were performed at the B3LYP/6-311++G(2df,2p) and diffusion monte carlo (DMC) levels for the refinement of the electronic energies. These calculations reveal important differences in the relative energies between high/low spin complexes obtained by B3LYP and DMC and show that although both DFT and DMC show similar trends, a higher level method such as DMC is necessary for a quantitative description of the interactions between Mn2+ and its natural ligands. (G)s of acetate-type ligand binding reactions were calculated that show that the higher the spin of the manganese complex, the lower the affinity for the ligand.

Resonance structures of the amide bond: the advantages of planarity

Delocalization indexes based on magnitudes derived from electron-pair densities are demonstrated to be useful indicators of electron resonance in amides. These indexes, based on the integration of the two-electron density matrix over the atomic basins defined through the zero-flux condition, have been calculated for a series of amides at the B3LYP/6-31+G* level of theory. These quantities, which can be viewed as a measure of the sharing of electrons between atoms, behave in concordance with the traditional resonance model, even though they are integrated in Bader atomic basins. Thus, the use of these quantities overcomes contradictory results from analyses of atomic charges, yet keeps the theoretical appeal of using nonarbitrary atomic partitions and unambiguously defined functions such as densities and pair densities. Moreover, for a large data set consisting of 24 amides plus their corresponding rotational transition states, a linear relation was found between the rotational barrier for the amide and the delocalization index between the nitrogen and oxygen atoms, indicating that this parameter can be used as an ideal physical-chemical indicator of the electron resonance in amides.

Theoretical Study of the Reaction of Vitamin B6 with1O2

Singlet oxygen is known to cause oxidative stress in cells, leading to severe damage (e.g., lipid peroxidation, membrane degradation, mutagenic alterations to DNA, protein misfunctionality). Recently, pyridoxine has been discovered to be capable of quenching singlet oxygen, however, the mechanism of this reaction remains essentially unknown. In this work, we have investigated four sets of reactions: 1) 1,3-addition to a double bond connected to a hydrogen-carrying group, resulting in the formation of allylic hydroperoxides; 2) [pi2+pi2] 1,2-cycloaddition to an isolated double bond, resulting in the formation of 1,2-peroxides; 3) 1,4-cycloaddition to a system containing at least two conjugated double bonds, resulting in the formation of the so-called 1,4-peroxides; 4) 1,4-addition to phenols and naphthols with the formation of hydroperoxide ketones. Thermodynamically, reaction 4 and the 6(9), 3(8), and 5(8) cases of reaction 1 are the most exergonic ones, with energies ranging from -16 to -18 kcal mol(-1). Furthermore, reaction 4 shows the lowest barrier through the reaction path, and is predicted to be the preferred mechanism for the pyridoxine + singlet-oxygen reaction, which is in agreement with previous experimental results.

Effects of OH Radical Addition on Proton Transfer in the Guanine−Cytosine Base Pair

Double proton transfer (PT) reactions in guanine-cytosine OH radical adducts are studied by the hybrid density functional B3LYP approach. Concerted and stepwise proton-transfer processes are explored between N1(H) on guanine (G) and N3 on cytosine (C), and between N4(H) on C and O6 on G. All systems except GC6OH display a concerted mechanism. 8OHGC has the highest dissociation energy and is 1.2 kcal/mol more stable than the nonradical GC base pair. The origin of the interactions are investigated through the estimation of intrinsic acid-basic properties of the *OH-X monomer (X = G or C). Solvent effects play a significant role in reducing the dissociation energy. The reactions including *OH-C adducts have significantly lower PT barriers than both the nonradical GC pair and the *OH-G adducts. All reactions are endothermic, with the GC6OH --> GC6OHPT reaction has the lowest reaction energy (4.6 kcal/mol). In accordance with earlier results, the estimated NBO charges show that the G moiety carries a slight negative charge (and C a corresponding positive one) in each adduct. The formation of a partial ion pair may be a potential factor leading to the PT reactions being thermodynamically unfavored.

Theoretical Modeling of Hydroxyl-Radical-Induced Lipid Peroxidation Reactions

The OH-radical-induced mechanism of lipid peroxidation, involving hydrogen abstraction followed by O2 addition, is explored using the kinetically corrected hybrid density functional MPWB1K in conjunction with the MG3S basis set and a polarized continuum model to mimic the membrane interior. Using a small nonadiene model of linoleic acid, it is found that hydrogen abstraction preferentially occurs at the mono-allylic methylene groups at the ends of the conjugated segment rather than at the central bis-allylic carbon, in disagreement with experimental data. Using a full linoleic acid, however, abstraction is correctly predicted to occur at the central carbon, giving a pentadienyl radical. The Gibbs free energy for abstraction at the central C11 is approximately 8 kcal/mol, compared to 9 kcal/mol at the end points (giving an allyl radical). Subsequent oxygen addition will occur at one of the terminal atoms of the pentadienyl radical fragment, giving a localized peroxy radical and a conjugated butadiene fragment, but is associated with rather high free energy barriers and low exergonicity at the CPCM-MPWB1K/MG3S level. The ZPE-corrected potential energy surfaces obtained without solvent effects, on the other hand, display considerably lower barriers and more exergonic reactions.

Theoretical Assessment of Naphazoline Redoxchemistry and Photochemistry

The imidazoline derivative naphazoline (2-(1-naphtylmethyl)-2-imidazoline) is an alpha2-adrenergic agonist used as non-prescription eye and nasal preparations. Besides its functionality in generating vascoconstriction and decongestion in the patient, the toxicity, ROS generating capability, and recently also possible antioxidant capacity of the compound have been reported in the literature. In the current work the structural and electronic features of the drug are explored, using computational chemical tools. Electron affinities, ionization potentials, and excitation energies are reported, as well as charge and spin distributions of various forms of the drug. The difference in photochemical behavior between the protonated and unprotonated (basic) species is explained by the molecular orbital distributions, allowing for efficient excitation quenching in the basic structure but clear naphthalene to imidazolene charge transfer upon HOMO--> LUMO excitation in the protonated form, enabling larger intersystem crossing capability to the imidazole localized excited triplet and a resulting higher singlet oxygen quantum yield.

Mechanism of Photoinduced Decomposition of Ketoprofen

UV-induced decarboxylation of the NSAID ketoprofen, followed by activation of molecular oxygen or formation of a decarboxylated peroxide adduct, is explored using computational quantum chemistry. The excited energy surfaces reveal that the neutral species will not decarboxylate, whereas the deprotonated acid decarboxylates spontaneously in the triplet state, and with an associated 3-5 kcal/mol barrier from several low-lying excited singlet states. The observed long lifetimes of the decarboxylated anion is explained in terms of the high stability of the triplet benzoyl ethyl species with protonated carbonylic oxygen, from which there is no obvious decay channel. Mechanisms for the generation of singlet oxygen and superoxide are discussed in detail. Addition of molecular oxygen to give the corresponding peroxyl radical capable of initiating propagating lipid peroxidation reactions is also explored. The computed data explains all features of the observed experimental observations made to date on the photodegradation of ketoprofen.