Electron-Transfer Induced Repair of 6-4 Photoproducts in DNA: A Computational Study

The mechanism employed by DNA photolyase to repair 6-4 photoproducts in UV-damaged DNA is explored by means of quantum chemical calculations. Considering the repair of both oxetane and azetidine lesions, it is demonstrated that reduction as well as oxidation enables a reversion reaction by creating anionic or cationic radicals that readily fragment into monomeric pyrimidines. However, on the basis of calculated reaction energies indicating that electron transfer from the enzyme to the lesion is a much more favorable process than electron transfer in the opposite direction, it is suggested that the photoenzymic repair can only occur by way of an anionic mechanism. Furthermore, it is shown that reduction of the oxetane facilitates a mechanism involving cleavage of the C-O bond followed by cleavage of the C-C bond, whereas reductive fragmentation of the azetidine may proceed with either of the intermonomeric C-N and C-C bonds cleaved as the first step. From calculations on neutral azetidine radicals, a significant increase in the free-energy barrier for the initial fragmentation step upon protonation of the carbonylic oxygens is predicted. This effect can be attributed to protonation serving to stabilize reactant complexes more than transition structures.

Identification of the brominated flame retardant 1,2-dibromo-4-(1,2-dibromoethyl)cyclohexane as an androgen agonist

To investigate androgen receptor (AR) activation by exogenous compounds, we used a combination of experimental analysis and theoretical modeling to compare a set of brominated flame retardants (BFRs) to dihydrotestosterone (DHT) with regard to ligand docking, AR binding, and AR activation in human hepatocellular liver carcinoma cells, as well as interacting energy analysis. Modeling of receptor docking was found to be a useful first step in predicting the potential to translocate to the ligand pocket of the receptor, and the computed interaction energy was found to correlate with the observed binding affinity. Flexible alignment studies of the BFR compounds demonstrated that 1,2-dibromo-4-(1,2-dibromoethyl)cyclohexane (BCH) closely overlap DHT. Combining the theoretical modeling with in vitro ligand-binding and receptor-activation assays, we show that BCH binds to and activates the human AR. The remaining BFRs did not successfully interact with the ligand pocket, were not able to replace a synthetic androgen from the receptor, and failed to activate the receptor.

The role of nucleobase carboradical and carbanion on DNA lesions: a theoretical study

DNA base release induced by H and OH radical addition to thymine and their corresponding electron adducts is studied at the DFT B3LYP/6-31+G(d,p) level in gas phase and in solution. H atom transfer after radical formation from C2' on the sugar to the C6 site on the base is shown to be prohibited for the radical species. Their corresponding electron adducts, albeit minor events in cellular systems, show excellent capabilities to proton transfer from C2' on the sugar to the C6 site on the base. The barriers for subsequent N-glycosidic bond dissociation range from 0.1 to 1.6 kcal mol(-1) at the B3LYP level and around 5 kcal mol(-1) using the BB1K functional, implying that these reactions can serve as a source to abasic sites. Analysis of bond dissociation energies show that all the reactions are exothermic, which is consistent with the changes in N-glycosidic bond lengths during the proton-transfer reactions. Bulk solvation plays a reverse influence on proton transfer and the bond rupture reactions. Molecular orbitals, NPA charges, and electron affinities are calculated to shed further light on the properties leading up to the intramolecular reactions.

Theoretical Study of the Antioxidant Properties of Pyridoxine

Molecules acting as antioxidants capable of scavenging reactive oxygen species (ROS) are of utmost importance in the living cell. The antioxidative properties of pyridoxine (vitamin B6) have recently been discovered. In this study, we have analyzed the reactivity of pyridoxine toward the ROS (.-)OH, (.-)OOH, and (.-)O(2)- at the density functional theory level (functionals B3LYP and MPW1B95). Two reaction types have been studied as follows: addition to the aromatic ring atoms and hydrogen/proton abstraction. Our results show that (.-)OH is the most reactive species, while (.-)OOH displays low reactivity and (.-)O2(-) does not react at all with pyridoxine. The most exergonic reactions are those where (.-)H is removed from the CH(2)OH groups or the ring-bound OH group and range from -33 to -39 kcal/mol. The most exergonic addition reactions occur by attacking the carbon atoms bonded to nitrogen but with an energy gain of only 6 kcal/mol.

Electronic excitation energies of Zn(i)S(i) nanoparticles

The excitation energies of small ZnS nanoclusters characterized in previous studies have been calculated using TDDFT. The relativistic pseudopotentials of Stevens et al have been used, including Zn 4s(2) electrons and S 3s(2) and 3p(4) electrons as valence electrons. Results obtained with these pseudopotentials are compared to those obtained considering also Zn 3s(2)3p(6)3d(10) electrons in the valence part, and demonstrated to be consistent. The results show that spheroid-like bubble structures have absorption energies in the range of 5-5.3 eV for small sizes, which decreases to 5 eV with increasing particle size.

pH-Dependent Electronic and Spectroscopic Properties of Pyridoxine (Vitamin B6)

The key electronic and spectroscopic properties of vitamin B(6) (pyridoxine) and some of its main charged and protonated/deprotonated species are explored using hybrid density functional theory (DFT) methods including polarized solvation models. It is found that the dominant species at low pH is the N(1)-protonated form and, at high pH, the O(3)(')-deprotonated compound. Computed and experimental UV-spectra for these species (experimental spectra recorded at pH 1.7 and 11.1, respectively) show a very close resemblance. At pH 4.3, the protonated species dominates, but with onset of the zwitterionic oxo form which is also the dominant species at neutral pH. The computational studies furthermore show that neither a polarized continuum model of the polar aqueous solvent or explicit hydrogen bonding through additional water molecules are sufficient to describe accurately the spectrum at physiological pH. Instead, Na(+) and Cl(-) counterions were required to give a blue-shift of approximately 0.15 eV.

Permeability of psoralen derivatives in lipid membranes

Molecular dynamics simulations have been performed to explore the distribution and translocation of a set of furocoumarins (psoralen derivatives) inside saturated and partially unsaturated lipid membranes. Within the simulations, strong accumulation of the photodynamic drugs is observed near the polar headgroup region, although the populations also extend out into the membrane/water interface as well as to the membrane center. The computed transverse (Dz) diffusion coefficients are in the range 0.01-0.03 x 10(-5) cm2 s(-1)-significantly slower than those reported for small molecules like water, ethane, and ammonia-and are related to the low mobility inside the polar headgroup region. Trimethylpsoralen (TMP) has a very low free energy barrier to transversion, only approximately 10 kJ/mol, whereas 5- and 8-methoxy psoralens (5-MOP, 8-MOP) have the largest barriers of the compounds studied-between 25 and 40 kJ/mol. Upper bounds to the permeation coefficients, obtained by integrating the resistance profiles across the bilayers, range from 5.2 x 10(-8) cm s(-1) for TMP to 4.1 x 10(-12) cm s(-1) for 5-MOP. The current simulations explain the high level of furocoumarin-lipid membrane complexes found in experimental studies of albino Wistar rats exposed to topical application of 8-MOP, and points to the possibility of membrane photodamage as a viable mechanism in psoralen ultraviolet-A treatment.

Photoreaction of skin-sensitizing trimethyl psoralen with lipid membrane models

The concerted photoreaction between trimethyl psoralen (TMP) and isomeric model compounds of oleic acid methyl ester (OAME) and eadilic acid methyl ester (EAME) are explored using density functional theory. The S(1) surfaces all reveal large barriers (25-31 kcal/mol) to reach the decay channels, whereas the S(2) surfaces differ considerably between the furan and pyrone side adducts. For the pyrone side adducts, a small barrier to reach the intersection between the S(1) and S(2) states is found, followed by a small second barrier on the S(1) surface to the low-lying decay channel. For the furan-side adducts, no such intersection between S(1) and S(2) is seen, which thus prevents these products from being formed in high yields. The photoinduced cycloaddition is more favorable between OAME and the pyrone side double bond of TMP than for any of the other systems, which agrees with the experimental findings that this is formed in the highest yield. The computed UV absorption spectra of TMP and the TMP(pyrone)-OAME cycloadduct agree well with the experimental spectra.

Sandwich Complexes Based on the “All-Metal” Al42− Aromatic Ring

We report on novel sandwichlike structures [Al(4)MAl(4)](q-) (q=0-2 and M=Ti, V, Cr, Zr, Nb, Mo, Hf, Ta and W) based on the recently synthesized all-metal aromatic Al(4)(2-) square ring. The sandwichlike structures have two aromatic tetraaluminum square rings which trap a transition-metal cation from either the first, second, or third row. The stability of the anionic sandwichlike complexes towards electron detachment is discussed, and addition of alkali cations is found to stabilize the 2- charged complexes, preventing spontaneous electron detachment. Once the sandwichlike complexes are formed, the Al(4)(2-) square properties remain nearly unchanged; this fact strongly supports the hypothesis that in these complexes the Al(4)(2-) square rings remain aromatic.

Protein-bound chromophores astaxanthin and phytochromobilin: excited state quantum chemical studies

We present an overview of excited state quantum chemical calculations aimed at elucidating controversial issues regarding the photochemistry of the protein-bound chromophores astaxanthin and phytochromobilin. In particular, we show how the application of time-dependent density functional theory and other single-reference quantum chemical excited state methods have contributed to shed new light on the origin of the >0.5 eV bathochromic shift of the electronic absorption by the carotenoid astaxanthin in the protein macromolecular complex crustacyanin, and the mechanism for C15-Z,syn --> C15-E,anti isomerization of the tetrapyrrole phytochromobilin that underlies the photoactivation of the plant photoreceptor phytochrome. Within the approximation that exciton coupling is neglected, the calculations on astaxanthin provide support for the notion that the bathochromic shift, which is responsible for the slate-blue coloration of lobster shell, is due to polarization rather than a conformational change of the chromophore in the protein-bound state. Furthermore, the polarization is attributed to a hydrogen-bonded protonated histidine residue. The calculations on phytochromobilin, in turn, suggest that a stepwise C15-Z,syn --> C15-E,syn (photochemical), C15-E,syn --> C15-E,anti (thermal) mechanism is much more favorable than a concerted, fully photochemical mechanism, and that neutral forms of the chromophore are much less likely to photoisomerize than the parent, protonated form. Accordingly, the calculations indirectly support the view that the photoactivation of phytochrome does not involve a proton transfer from the chromophore to the surrounding protein.

Dehydration of ribonucleotides catalyzed by ribonucleotide reductase: the role of the enzyme

This article focuses on the second step of the catalytic mechanism for the reduction of ribonucleotides catalyzed by the enzyme Ribonucleotide Reductase (RNR). This step corresponds to the protonation/elimination of the substrate's C-2' hydroxyl group. Protonation is accomplished by the neighbor Cys-225, leading to the formation of one water molecule. This is a very relevant step since most of the known inhibitors of this enzyme, which are already used in the fight against certain forms of cancer, are 2'-substituted substrate analogs. Even though some theoretical studies have been performed in the past, they have modeled the enzyme with minimal gas-phase models, basically represented by a part of the side chain of the relevant amino acids, disconnected from the protein backbone. This procedure resulted in a limited accuracy in the position and/or orientation of the participating residues, which can result in erroneous energetics and even mistakes in the choice of the correct mechanism for this step. To overcome these limitations we have used a very large model, including a whole R1 model with 733 residues plus the substrate and 10 A thick shell of water molecules, instead of the minimal gas-phase models used in previous works. The ONIOM method was employed to deal with such a large system. This model can efficiently account for the restrained mobility of the reactive residues, as well as the long-range enzyme-substrate interactions. The results gave additional information about this step, which previous small models could not provide, allowing a much clearer evaluation of the role of the enzyme. The interaction energy between the enzyme and the substrate along the reaction coordinate and the substrate steric strain energy have been obtained. The conclusion was that the barrier obtained with the present model was very similar to the one previously determined with minimal gas-phase models. Therefore, the role of the enzyme in this step was concluded to be mainly entropic, rather than energetic, by placing the substrate and the two reactive residues in a position that allows for the highly favorable concerted trimolecular reaction, and to protect the enzyme radical from the solvent.

Photodegradation of Substituted Stilbene Compounds: What Colors Aging Paper Yellow?

Photodegradation of lignin is one of the major postprocessing problems in paper production, as this renders yellowing of the paper and reduced paper quality. In this study, we have explored the photochemical properties of substituted stilbene derivatives believed to be key chromophores in the photodegradation of lignin derived from cinnamyl alcohol. In particular, the present work focuses on the computation of UV/vis electronic absorption spectra for different methoxylated stilbenes and their proposed photodegradation products. All calculations were performed using the time-dependent formalism of density functional theory (TD-DFT) and the B3LYP hybrid functional. It is concluded that the methodology employed is capable of reproducing not only the overall spectra, but also subtle features owing to the effects of different substitution patterns. For the strongly absorbing first excited singlet state (HOMO --> LUMO excitation) of the methoxylated stilbenes, the calculated transition energies are, albeit somewhat fortuitously, in excellent agreement with experimental data. The light-induced yellowing indirectly caused by the presence of stilbenes can be rationalized in terms of the absorption spectra of the resulting photodegraded o-quinones, for which distinct transitions in the 420-500 nm region of the visible spectrum lacking prior to degradation are observed.

Theoretical Study of Cisplatin Binding to DNA: The Importance of Initial Complex Stabilization

The first and second substitution reactions between activated (hydrolyzed) cisplatin, Pt(NH3)2(H2O)2(2+), and purine bases guanine and adenine are explored using the B3LYP hybrid functional, IEF-PCM solvation models, and large basis sets. The computed free energy barrier for the first substitution is 19.5 kcal/mol for guanine (exptl value = 18.3 kcal/mol) and 24.0 kcal/mol for adenine. The observed predominance toward guanine in the first substitution is explained in terms of significantly larger stabilization energy for the initially formed complex, compared with adenine, in combination with favored kinetics, and represents a revised view of the proposed mechanism for cisplatin binding to DNA. For the second substitution, the computed barrier for Pt(NH3)2G2(2+) head-to-head formation is 22.5 kcal/mol, in very good agreement with experimental data for adduct closure (23.4 kcal/mol). Again, a higher stability in complexation with G over A is ascribed as the main contributing factor favoring G over A substitution. The calculations provide a first explanation for the predominance of 1,2-d(GpG) over 1,2-d(ApG) intrastrand didentate adducts, and the origin of the 5'-3' direction specificity of the 1,2-d(ApG) adducts.

Ribonucleotide activation by enzyme ribonucleotide reductase: understanding the role of the enzyme

This article focuses on the first step of the catalytic mechanism for the reduction of ribonucleotides catalyzed by the enzyme Ribonucleotide Reductase (RNR). This corresponds to the activation of the substrate. In this work a large model of the active site region involving 130 atoms was used instead of the minimal gas phase models used in previous works. The ONIOM method was employed to deal with such a large system. The results gave additional information, which previous small models could not provide, allowing a much clearer evaluation of the role of the enzyme in this step. Enzyme-substrate interaction energies, specific transition state stabilization, and substrate steric strain energies were obtained. It was concluded that the transition state is stabilized in 4.0 kcal/mol by specific enzyme-substrate interactions. However, this stabilization is cancelled by the cost in conformational energy for the enzyme to adopt the transition state geometry; the overall result is that the enzyme machinery does not lead to a rate enhancement in this step. It was also found that the substrate binds to the active site with almost no steric strain, emphasizing the complementarity and specificity of the RNR active site for nucleotide binding. The main role of the enzyme at the very beginning of the catalytic cycle was concluded to be to impose stereospecifity upon substrate activation and to protect the enzyme radical from the solvent, rather than to be an reaction rate enhancement.

On the formation of cyclobutane pyrimidine dimers in UV-irradiated DNA: why are thymines more reactive?

The reaction pathways for thermal and photochemical formation of cyclobutane pyrimidine dimers in DNA are explored using density functional theory techniques. Although it is found that the thermal [2 + 2] cycloadditions of thymine + thymine (T + T --> T x T), cytosine + cytosine (C + C --> C x C) and cytosine + thymine (C + T --> C x T) all are similarly unfavorable in terms of energy barriers and reaction energies, the excited-state energy curves associated with the corresponding photochemical cycloadditions display differences that--in line with experimental findings--unanimously point to the predominance of T x T in UV-irradiated DNA. It is shown that the photocycloaddition of thymines is facilitated by the fact that the S1 state of the corresponding reactant complex lies comparatively high in energy. Moreover, at a nuclear configuration coinciding with the ground-state transition structure, the excited-state energy curve displays an absolute minimum only for the T + T system. Finally, the T + T system is also associated with the most favorable excited-state energy barriers and has the smallest S2-S0 energy gap at the ground-state transition structure.

Theoretical studies of the cross-linking mechanisms between cytosine and tyrosine

DNA-protein cross-linking is one of the many DNA lesions mediated by hydroxyl radicals, the most damaging among the reactive oxygen species in biological systems. Density functional theory methods are employed to investigate the complex reaction mechanisms of the formation of cytosine-tyrosine cross-links as observed in gamma-irradiated aqueous solutions of cytosine and tyrosine, as well as in gamma-irradiated nucleohistone. The majority of the radical addition mechanisms considered are found to have significant barriers and therefore to be thermodynamically unfavorable for the formation of the initial cross-linked product. Our calculated reaction potential energy surfaces suggest that a feasible complete mechanism consists of radical combination forming the initial cross-linked product, a hydrogen shuffle within the initial cross-linked product, and an acid-catalyzed dehydration reaction. Water and hydrogen-bonding interactions are suggested to play a key role in catalyzing the hydrogen-transfer step of the reaction.