Synthesis,characterization, and coordination chemistry of the 2-azaphenalenyl radical

The 2-azaphenalenyl radical 2 has been synthesized and characterized by ESR spectroscopy. Variable-temperature ESR measurements were carried out on both the phenalenyl (1) and the 2-azaphenalenyl (2) radicals. The phenalenyl radical 1 has the known propensity to dimerize at temperatures below 20 degrees C, but unexpectedly less so than originally reported. The first experimental measurement of bond dissociation enthalpy for the dimerization of the phenalenyl radical 1 was obtained in CCl(4) (11.34 +/- 0.11 kcal/mol) and toluene (9.8 +/- 0.7 kcal/mol). The 2-azaphenalenyl radical 2 does not show a propensity to dimerize over the measurable temperature range (220-330 K), but does so in the presence of Cu(hfac)(2) (hfac = hexafluoroacetylacetonate). The latter complex was characterized by X-ray crystallography.     

Structure/energy correlation of bowl depth and inversion barrier in corannulene derivatives: combined experimental and quantum mechanical analysis.

Synthesis of a series of corannulene derivatives with varying bowl depths has allowed for a study correlating the structure (bowl depth) and the energy of bowl inversion. Substituents placed in the peri positions are repulsive and flatten the bowl, thus causing a decrease in the bowl inversion barrier. Conversely, annelation across the peri positions causes a deepening of the bowl, thus an increase in the bowl inversion barrier. Barriers between 8.7 and 17.3 kcal/mol have been measured, and their structures have been calculated using a variety of ab initio methods. The energy profile of an individual corannulene derivative is assumed to fit a mixed quartic/quadratic function from which an empirical correlation of bowl depth and inversion barrier that follows a quartic function is derived. Structure/energy correlations of this type speak broadly of the nature of enzymatic and catalytic activation of substrates.     

Corannulene as a Lewis base: computational modeling of protonation and lithium cation binding.

A computational modeling of the protonation of corannulene at B3LYP/6-311G(d,p)//B3LYP/6-311G(d,p) and of the binding of lithium cations to corannulene at B3LYP/6-311G(d,p)//B3LYP/6-31G(d,p) has been performed. A proton attaches preferentially to one carbon atom, forming a sigma-complex. The isomer protonated at the innermost (hub) carbon has the best total energy. Protonation at the outermost (rim) carbon and at the intermediate (bridgehead rim) carbon is less favorable by ca. 2 and 14 kcal mol(-)(1), respectively. Hydrogen-bridged isomers are transition states between the sigma-complexes; the corresponding activation energies vary from 10 to 26 kcal mol(-)(1). With an empirical correction obtained from calculations on benzene, naphthalene, and azulene, the best estimate for the proton affinity of corannulene is 203 kcal mol(-)(1). The lithium cation positions itself preferentially over a ring. There is a small energetic preference for the 6-ring over the 5-ring binding (up to 2 kcal mol(-)(1)) and of the convex face over the concave face (3-5 kcal mol(-)(1)). The Li-bridged complexes are transition states between the pi-face complexes. Movement of the Li(+) cation over either face is facile, and the activation energy does not exceed 6 kcal mol(-)(1) on the convex face and 2.2 kcal mol(-)(1) on the concave face. In contrast, the transition of Li(+) around the corannulene edge involves a high activation barrier (24 kcal mol(-)(1) with respect to the lowest energy pi-face complex). An easier concave/convex transformation and vice versa is the bowl-to-bowl inversion with an activation energy of 7-12 kcal mol(-)(1). The computed binding energy of Li(+) to corannulene is 44 kcal mol(-)(1). Calculations of the (7)Li NMR chemical shifts and nuclear independent chemical shifts (NICS) have been performed to analyze the aromaticity of the corannulene rings and its changes upon protonation.     

Density functional investigation on electron-transfer catalysis of cycloreversion of cyclobutane: radical anion mechanism

The mechanism of cycloreversion of cyclobutane radical anion (c-C(4)H(8) (-)) has been investigated at the UB3LYP/6-31++G(d,p) level, and compared with those of neutral c-C(4)H(8) and c-C(4)H(8) (+) radical cation. Although both c-C(4)H(8) (-) and C(2)H(4) are shown to be Rydberg states unstable with respect to electron ejection, the activation barrier for the "rotating" cycloreversion of c-C(4)H(8) (-) (37.3 kcal/mol) is lower by about 25.2 kcal/mol than that of c-C(4)H(8), and even the intervention of tetramethylene radical anion intermediate may reduce the activation barrier for the cycloreversion of c-C(4)H(8) by about 8.4 kcal/mol, mainly due to stronger electron-deficiency of intermediate biradical species than close-shell cyclobutanes. For the cycloreversion for c-C(4)H(8) (-), side isomerization reaction may be efficiently prevented by the low kinetic stability of tetramethylene radical anion intermediate towards dissociation, just different from the radical cation case. Our theoretical results have suggested the possibility of electron-attachment catalysis of the cycloreversion of some electron-deficient substituted cyclobutanes.     

Interaction between ions and substituted buckybowls: A comprehensive computational study

Complexes formed by substituted buckybowls derived from corannulene and sumanene with sodium cation or chloride anion have been computationally studied by using a variety of methods. Best results have been obtained with the SCS‐MP2 method extrapolated to basis set limit, which reproduces the highest‐level values obtained with the MP2.X method. All bowls form stable complexes with chloride anion, with stabilities ranging from ?6 kcal/mol in the methylated corannulene derivative to ?45 kcal/mol in the CN‐substituted sumanene. The opposite trend is observed in sodium complexes, going from deeply attractive complexes with the methylated derivatives (?36 kcal/mol with sumanene derivative) to slightly repulsive ones in the CN‐substituted bowls (2 kcal/mol in the corannulene derivative). Anion complexes are stabilized by large electrostatic interactions combined with smaller though significant dispersion and induction contributions. Conversely, cation complexes are stabilized by large induction contributions capable of holding together the bowl and the cation even in cases where the electrostatic interaction is repulsive. The effect of substitution is mainly reflected on changes in the molecular electrostatic potential of the bowl and, thus, in the electrostatic contribution to the interaction. Therefore, the variations in the stability of the complexes on substitution could be roughly predicted just considering the changes in the electrostatic interaction. However, other contributions also register changes mainly as a consequence of displacements on the position of the ion at the minimum, so the accurate prediction of the stability of this kind of complexes requires going further than the electrostatic approach. © 2014 Wiley Periodicals, Inc.     

On the aromatic stabilization of corannulene and coronene

The application of set of homodesmotic reactions allowed us to estimate the aromatic stabilization energy (ASE) of corannulene and coronene. Appropriate reactions have been applied to balance syn/anti mismatches in di-, tetra- and hexamethylene substituted derivatives. Based on many different polycyclic reference structures that compensate the effect of strain in the corannulene moiety the value of ASE comes to 44.5 kcal mol(-1). Planar corannulene is more stabilized by cyclic π-electron delocalization by ca. 10.7 kcal mol(-1), as compared with a bowl-shaped system. A similar approach for coronene leads to an ASE equal to 58.4 kcal mol(-1).     

Characterizing the dimerizations of phenalenyl radicals by ab initio calculations and spectroscopy: sigma-bond formation versus resonance pi-stabilization

Electronic-structure calculations for the self-association of phenalenyl radical (P*) predict the formation of dimeric species (sigma-P2) in which both moieties are connected by a sigma-bond with rP-P approximately 1.59 A and bond dissociation enthalpy of DeltaH(D) approximately 16 kcal mol(-1). Such an unusually weak sigma-bond is related to the loss of aromatic stabilization energy of approximately 34 kcal mol(-1) per phenalenyl moiety, largely owing to rehybridization. Ab initio calculations also reveal that the corresponding (one-electron) bond between phenalenyl radical and its closed-shell cation in sigma-P2+* is unstable relative to dissociation. Time-dependent DFT computations indicate the absence of any (strongly allowed) electronic transition in the visible region of the absorption spectrum of phenalenyl sigma-dimer. Such theoretical predictions are supported by experimental (ESR and UV-NIR) spectroscopic studies, in which the availability of a series of sterically hindered phenalenyl radicals allows definitive separations of the sigma-dimerization process from interference by pi-dimerization. As such, the thermodynamic parameters (determined from the temperature dependence of the ESR signals) with DeltaH(D) = 14 kcal mol(-1) and DeltaS(D) = 52 e.u. can be assigned to the formation of the colorless sigma-dimer. Similar results are obtained for all phenalenyl derivatives (provided their substitution patterns allow sigma-bond formation) to confirm the energetic preference of sigma-dimerization over pi-dimerization.     

Molecular orbital calculations of ring opening of the isoelectronic cyclopropylcarbinyl radical, cyclopropoxy radical, and cyclopropylaminium radical cation series of radical clocks

Detailed molecular orbital calculations were directed to the cyclopropylcarbinyl radical (1), the cyclopropoxy radical (2), and the cyclopropylaminium radical cation (3) as well as their ring-opened products. Since a considerable amount of data are published about cyclopropylcarbinyl radicals, calculations were made for this species and related ring-opened products as a reference for 2 and 3 and their reactions. Radicals 1-3 have practical utility as "radical clocks" that can be used to time other radical reactions. Radical 3 is of further interest in photoelectron-transfer processes where the back-electron-transfer process may be suppressed by rapid ring opening. Calculations have been carried out at the UHF/6-31G*, MP4//MP2/6-31G*, DFT B3LYP/6-31G*, and CCSD(T)/cc-pVTZ//QCISD/cc-pVDZ levels. Energies are corrected to 298 K, and the barriers between species are reported in terms of Arrhenius E(a) and log A values along with differences in enthalpies, free energies, and entropies. The CCSD(T)-calculated energy barrier for ring opening of 1 is E(a) = 9.70, DeltaG* = 8.49 kcal/mol, which compares favorably to the previously calculated value of E(a) = 9.53 kcal/mol by the G2 method, but is higher than an experimental value of 7.05 kcal/mol. Our CCSD(T)-calculated E(a) value is also higher by 1.8 kcal/mol than a previously reported CBS-RAD//B3LYP/6-31G* calculation. The cyclopropoxy radical has a very small barrier to ring opening (CCSD(T), E(a) = 0.64 kcal/mol) and should be a very sensitive time clock. Of the three series studied, the cyclopropylaminium radical cation is most complex. In agreement with experimental data, bisected cyclopropylaminium radical cation is not found, but instead a ring-opened species is found. A perpendicular cyclopropylaminium radical cation (4) was found as a transition-state structure. Rotation of the 2p orbital in 4 to the bisected array results in ring opening. The minimum onset energy of photoionization of cyclopropylamine was calculated to be 201.5 kcal/mol (CCSD(T)) compared to experimental values of between about 201 and 204 kcal/mol. Calculations were made on the closely related cyclopropylcarbinyl and bicyclobutonium cations. Stabilization of the bisected cyclopropylcarbinyl conformer relative to the perpendicular species is much greater for the cations (29.1 kcal/ mol, QCISD) compared to the radicals (3.10 kcal/mol, QCISD). A search was made for analogues to the bicyclobutonium cation in the radical series 1 and 2 and the radical cation series 3. No comparable species were found. A rationale was made for some conflicting calculations involving the cyclopropylcarbinyl and bicyclobutonium cations. The order of stability of the cyclopropyl-X radicals was calculated to be X = CH2 > X = O > X = NH2+, where the latter species has no barrier for ring opening. The relative rate of ring opening for cyclopropyl-X radicals X = CH2 to X = O was calculated to be 3.1 x 10(6) s(-1) at 298 K (QCISD).     

Third-order nonlinear optical properties of open-shell supermolecular systems composed of acetylene linked phenalenyl radicals

The third-order nonlinear optical (NLO) properties, at the molecular level, the static second hyperpolarizabilities, γ, of supermolecular systems composed of phenalenyl and pyrene rings linked by acetylene units are investigated by employing the long-range corrected spin-unrestricted density functional theory, LC-UBLYP, method. The phenalenyl based superethylene, superallyl, and superbutadiene in their lowest spin states have intermediate diradical characters and exhibit larger γ values than the closed-shell pyrene based superpolyene systems. The introduction of a positive charge into the phenalenyl based superallyl radical changes the sign of γ and enhances its amplitude by a factor of 35. Although such sign inversion is also observed in the allyl radical and cation systems in their ground state equilibrium geometries, the relative amplitude of γ is much different, that is, |γ(regular allyl cation)/γ(regular allyl radical)| = 0.61 versus |γ(phenalenyl based superallyl cation)/γ(phenalenyl based superallyl radical)| = 35. In contrast, the model ethylene, allyl radical/cation, and butadiene systems with stretched carbon-carbon bond lengths (2.0 ?), having intermediate diradical characters, exhibit similar γ features to those of the phenalenyl based superpolyene systems. This exemplifies that the size dependence of γ as well as its sign change by introducing a positive charge on the phenalenyl based superpolyene systems originate from their intermediate diradical characters. In addition, the change from the lowest to the highest π-electron spin states significantly reduces the γ amplitudes of the neutral phenalenyl based superpolyene systems. For phenalenyl based superallyl cation, the sign inversion of γ (from negative to positive) is observed upon switching between the singlet and triplet states, which is predicted to be associated with a modification of the balance between the positive and negative contributions to γ. The present study paves the way toward designing a variety of open-shell NLO supermolecular systems composed of phenalenyl radical building blocks.     

Electronic Factors Influencing the Decarboxylation of beta-Keto Acids. A Model Enzyme Study

A theoretical study of the mechanism of decarboxylation of beta-keto acids is described. A cyclic transition structure was found with essentially complete proton transfer from the carboxylic acid to the beta-carbonyl group. The activation barrier for decarboxylation of formylacetic acid is predicted to be 28.6 kcal/mol (MP4SDTQ/6-31G//MP2/6-31G) while loss of CO(2) from its anion exhibits a barrier of only 20.6 kcal/mol (MP4SDTQ/6-31+G//MP2/6-31+G). Barrier heights of decarboxylation of malonic acid and alpha,alpha-dimethylacetoacetic acid are predicted to be 33.2 and 26.7 kcal/mol, respectively. Model enzyme studies using a thio methyl ester of malonate anion suggests that the role of malonyl-CoA is to afford a polarizable sulfur atom to stabilize the developing enolate anion in the transition structure for decarboxylation. Adjacent positively charged ammonium ions are also observed to stabilize the loss of CO(2) from a carboxylate anion by through-bond Coulombic stabilization of the transition structure.     

Conformational analysis of 1-butyl-3-methylimidazolium by CCSD(T) level ab initio calculations: effects of neighboring anions

Conformational energies for the butyl group of 1-butyl-3-methylimidazolium (bmim) were calculated by high-level ab initio methods. Estimated relative energies for the TT, GT and G'T rotamers of an isolated bmim cation at the CCSD(T)/cc-pVTZ level are 0.0 -0.02 and -0.50 kcal/mol, respectively. The close contact of a Cl anion to theC(2)-H of imidazolium considerably increases the relative stability of the GT rotamer. Estimated relative energies for the three rotamers of the [bmim]Cl complex, in which the Cl anion exists close to the C(2)-H, are 0.0, -1.61 and -0.25 kcal/mol, respectively. The GT rotamer is favored by the strong attractive electrostatic interaction between the bmim cation and Cl anion. The C(2)-H group in the GT rotamer has a larger positive charge compared with those in the TT and G'T rotamers. The contact of a Br anion to the C(2)-H also stabilizes the GT rotamer. The effects of the Cl anion close to the C(4)-Hand C(5)-Hare small. The anion effects suggest that the GT rotamer is the most stable in ionic liquids. The positive charge on imidazolium ring does not largely change the conformational energies. Estimated relative energies for the three rotamers of N-butylimidazole (0.0, -0.29 and -0.75 kcal/mol, respectively) are not largely different from those for isolated bmim. Calculated MP2/cc-pVTZ level torsional potential for the C im-N im-C-C bond has a minimum when the torsional angle is close to 90 degrees. Coplanar conformation is not a stable structure. Calculated torsional barrier height between the two nonplanar minima is less than 1 kcal/mol.     

Methyl and silyl mesolytic dissociations in the radical cations and radical anions of but-1-ene,allylsilane, hexa-1,3-diene,and penta-2,4-dienylsilane. CAS-MCSCF and coupled cluster theoretical study

Methyl or silyl dissociation in the CH(2)=CHCH(2)-XH(3) (a-XH(3)(*)(+)) and CH(2)=CHCH=CHCH(2)-XH(3) (p-XH(3)(*) (+)) radical cations (X = C, Si) yields a(+) or p(+) and XH(3)(*). Similarly, the radical anions a-CH(3)(*) (-) and p-CH(3)(*) (-) give the pi-delocalized anion and CH(3)(*) preferentially. In contrast, a-SiH(3)(*) (-) and p-SiH(3)(*-) prefer to dissociate into the pi-delocalized radical and silide. All reactions are endoergic: by 43-50 kcal mol(-)(1) in the radical cations, and easier to some extent in the radical anions, that require 29-33 (X = C) and 13-14 kcal mol(-)(1) (X = Si). The fragmentation energy profiles do not present significant barriers for the backward process in the case of the radical cations. All radical anions exhibit an energy maximum along the dissociation pathway, but the barrier is lower than the dissociation limit. Fragmentation is "activated" more in the anions than in the cations with respect to homolysis in the corresponding neutrals (that requires 72-81 kcal mol(-)(1)). Wave function analysis indicates that the C-X bond cleavage in the hydrocarbon radical ions, although formally comparable to a homolytic process, is at variance with this model, due to the spin recoupling of one of the two C-X bond electrons with the originally unpaired electron. This is basically true also for the silyl-substituted radical anions, in which the initial more delocalized charge distribution might suggest some heterolytic character of the bond cleavage.     

Computational insights into the mechanism of radical generation in B12-dependent methylmalonyl-CoA mutase

ONIOM calculations have provided novel insights into the mechanism of homolytic Co-C5' bond cleavage in the 5'-deoxyadenosylcobalamin cofactor catalyzed by methylmalonyl-CoA mutase. We have shown that it is a stepwise process in which conformational changes in the 5'-deoxyadenosine moiety precede the actual homolysis step. In the transition state structure for homolysis, the Co-C5' bond elongates by approximately 0.5 Angstroms from the value found in the substrate-bound reactant complex. The overall barrier to homolysis is approximately 10 kcal/mol, and the radical products are approximately 2.5 kcal/mol less stable than the initial ternary complex of enzyme, substrate, and cofactor. The movement of the deoxyadenosine moiety during the homolysis step positions the resulting 5'-deoxyadenosyl radical for the subsequent hydrogen atom transfer from the substrate, methylmalonyl-CoA.     

Is the protein surrounding the active site critical for hydrogen peroxide reduction by selenoprotein glutathione peroxidase? An ONIOM study

In this ONIOM(QM:MM) study, we evaluate the role of the protein surroundings in the mechanism of H2O2 reduction catalyzed by the glutathione peroxidase enzyme, using the whole monomer (3113 atoms in 196 amino acid residues) as a model. A new optimization scheme that allows the full optimization of transition states for large systems has been utilized. It was found that in the presence of the surrounding protein the optimized active site structure bears a closer resemblance to the one in the X-ray structure than that without the surrounding protein. H2O2 reduction occurs through a two-step mechanism. In the first step, the selenolate anion (E-Se(-)) formation occurs with a barrier of 16.4 kcal/mol and is endothermic by 12.0 kcal/mol. The Gln83 residue plays the key role of the proton abstractor, which is in line with the experimental suggestion. In the second step, the O-O bond is cleaved, and selenenic acid (R-Se-OH) and a water molecule are formed. The calculated barrier for this process is 6.0 kcal/mol, and it is exothermic by 80.9 kcal/mol. The overall barrier of 18.0 kcal/mol for H2O2 reduction is in reasonable agreement with the experimentally measured barrier of 14.9 kcal/mol. The protein surroundings has been calculated to exert a net effect of only 0.70 kcal/mol (in comparison to the "active site only" model including solvent effects) on the overall barrier, which is most likely due to the active site being located at the enzyme surface.     

On the possibility of the detachment of hydrogen as a result of electron capture by a Broensted center on zeolites

The structure and reactivity of the radical anion center in ZSM-5 zeolite were studied by the density functional theory method. It was shown that the interaction of the hydrogen zeolite form with adsorbed olefins and aromatic hydrocarbons could be accompanied by electron transfer from hydrocarbon molecules to the Broensted acid center with the formation of a radical anion fragment. The radical anion fragment formed is unstable, which contributes to the probability of the exothermic process (ΔE = ?21 kcal/mol) of atomic hydrogen detachment with the activation energy not exceeding 10–12 kcal/mol. The atomic hydrogen split off can initiate hydrocarbon transformations, and such a radical anion center can play the role of a catalytic activity carrier on acid zeolites.     

Understanding the planar tetracoordinate carbon atom: spiropentadiene dication

The atoms in molecule theory shows that the spiropentadiene dication has a planar tetracoordinate carbon (ptC) atom stabilized mainly through the sigma bonds and this atom has a negative charge. The bonds to the ptC atom have less covalent character than the central carbon from neutral spiropentadiene. The total positive charge is spread along the structure skeleton. The analysis of the potential energy surface shows that the dication spiropentadiene has a 2.3 kcal/mol activation barrier for ring opening.     

The possibility of the decomposition of 2'-deoxyribose moiety of thymidine induced by the low energy electron attachment

To evaluate the possibility of the decomposition of 2-deoxyribose moiety of thymidine induced by low energy electrons (LEE) attachment, the transition states and the energy barriers of the bond breaking processes of the ribose of the nucleoside have been studied theoretically by applying the density functional theory with the double zeta basis sets (DZP++). The energy barriers for the breakage of the C-C bonds (C(1')-C(2'), C(2')-C(3'), C(3')-C(4'), and C(4')-C(5')) of the ribose group of the radical anion of thymidine are found to be high (ca. 42-57 kcal/mol). The total energies of the C-C bond-broken products are significantly higher than that of the radical anion dT(*-). The decomposition of dT(*-) through the C-C bond rupture is unlikely to take place. The rupture of the C(1')-O(4') bond of dT(*-) needs an activation energy as low as 10.4 kcal/mol. However, the reversed reaction (C(1')-O(4') bond formation) needs the activation energy low as 0.3 kcal/mol. Therefore, the intermediate product LM1(C1')-(O4') is unlikely to be stable and the C(1')-O(4') bond-broken is not favored. The activation energy of the C(4')-O(4') bond rupture process amounts to 20.5 kcal/mol. The total energy of the C(4')-O(4') bond broken product is about 6.5 kcal/mol lower than that of the reactant dT(*-). The subsequent N1-glycosidic bond breaking process is found to have a very low energy barrier. Therefore, the LEE-induced base release through the C(4')-O(4') bond rupture might be a possible pathway.     

McLafferty rearrangement of the radical cations of butanal and 3-fluorobutanal: a theoretical investigation of the concerted and stepwise mechanisms

The stepwise and concerted pathways for the McLafferty rearrangement of the radical cations of butanal (Bu(+)) and 3-fluorobutanal (3F-Bu(+)) are investigated with density functional theory (DFT) and ab initio methods in conjunction with the 6-311+G(d,p) basis set. A concerted transition structure (TS) for Bu(+), (H), is located with a Gibbs barrier height of 37.7 kcal/mol as computed with CCSD(T)//BHandHLYP. Three pathways for the stepwise rearrangement of Bu(+) have been located, which are all found to involve different complexes. The barrier height for the H(gamma) transfer is found to be 2.2 kcal/mol, while the two most favorable TSs for the C(alpha)-C(beta) cleavage are located 8.9 and 9.2 kcal/mol higher. The energies of the 3F-Bu(+) system have been calculated with the promising hybrid meta-GGA MPWKCIS1K functional of DFT. Interestingly, the fluorine substitution yields a barrier height of only 20.5 kcal/mol for the concerted TS, (3F-H). A smaller computed dipole moment, 12.1 D, for (3F-H) compared with 103.2 D for (H) might explain the stabilization of the substituted TS. The H(gamma) transfer, with a barrier height of 4.9 kcal/mol, is found to be rate-determining for the stepwise McLafferty rearrangement of 3F-Bu(+), in contrast to the unsubstituted case. By inspection of the spin and charge distributions of the stationary points, it is noted that the bond cleavages in the concerted rearrangements are mainly of heterolytic nature, while those in the stepwise channels are found to be homolytic.     

Electron attachment-induced DNA single strand breaks: C3'-O3' sigma-bond breaking of pyrimidine nucleotides predominates

A detailed understanding of DNA strand breaks induced by low energy electrons (LEE) is of crucial importance for the advancement of many areas of molecular biology and medicine. To elucidate the mechanism of DNA strand breaks by LEEs, theoretical investigations of the electron attachment-induced C3'-O3' sigma-bond breaking of the pyrimidine nucleotides have been performed. Calculations of 2'-deoxycytidine-3'-monophosphate and 2'-deoxythymidine-3'-monophosphate in their protonated form (denoted as 3'-dCMPH and 3'-dTMPH) have been carried out with the reliably calibrated B3LYP/DZP++ theoretical approach. Our results demonstrate that the transfer of the negative charge from the pi*-orbital of the radical anion of pyrimidines to the DNA backbone does not pass through the N1-glycosidic bond. Instead, the migration of the excessive negative charge through the atomic orbital overlap between the C6 of pyrimidine and the C3' of ribose most likely represents a pathway that subsequently leads to the strand breaks. The proposed mechanism of the LEE-induced single strand breaks in DNA assumes that the formation of the base-centered radical anions is the first step in this process. Subsequently, these electronically stable radical anions may undergo either C-O bond breaking or N-glycosidic bond rupture. The present investigation of 3'-dCMPH and 3'-dTMPH yields an energy barrier of 6.2-7.1 kcal/mol for the C3'-O3' sigma-bond cleavage. This is much lower than the energy barriers required for the C5'-O5' sigma-bond and the N1-glycosidic bond break. Therefore, we conclude that the C3'-O3' sigma-bond rupture dominates the LEE-induced single strand breaks of DNA.     

Computational investigation of adsorption of molecular hydrogen on lithium-doped corannulene

Density functional theory and classical molecular dynamics simulations are used to investigate the prospect of lithium-doped corannulene as adsorbent material for H(2) gas. Potential energy surface scans at the level of B3LYP/6-311G(d,p) show an enhanced interaction of molecular hydrogen with lithium-atom-doped corannulene complexes with respect to that found in undoped corannulene. MP2(FC)/6-31G(d,p) optimizations of 4H(2)-(Li(2)-C(20)H(10)) yield H(2) binding energies of -1.48 kcal/mol for the H(2)-Li interaction and -0.92 kcal/mol for the H(2)-C interaction, whereas values of -0.94 and -0.83 kcal/mol were reported (J. Phys. Chem. B 2006, 110, 7688-7694) for physisorption of H(2) on the concave and the convex side of corannulene using MP2(full)/6-31G(d), respectively. Classical molecular dynamics simulations predict hydrogen uptakes in Li-doped corannulene assemblies that are significantly enhanced with respect to that found in undoped molecules, and the hydrogen uptake ability is dependent on the concentration of lithium dopant. For the Li(6)-C(20)H(10) complex, a hydrogen uptake of 4.58 wt % at 300 K and 230 bar is obtained when the adsorbent molecules are arranged in stack configurations separated by 6.5 A, and with interlayer distances of 10 A, hydrogen uptake reaches 6.5 wt % at 300 K and 215 bar.