C-H···O Hydrogen Bonding. The Prototypical Methane-Formaldehyde System: A Critical Assessment.

Distinguishing the functionality of C-H···O hydrogen bonds (HBs) remains challenging, because their properties are difficult to quantify reliably. Herein, we present a study of the model methane-formaldehyde complex (MFC). Six stationary points on the MFC potential energy surface (PES) were obtained at the CCSD(T)/ANO2 level. The CCSDT(Q)/CBS interaction energies of the conformers range from only -1.12 kcal mol-1 to -0.33 kcal mol-1, denoting a very flat PES. Notably, only the lowest energy stationary point (MFC1) corresponds to a genuine minimum, whereas all other stationary points-including the previously studied ideal case of ae(C-H···O) = 180°-exhibit some degree of freedom that leads to MFC1. Despite the flat PES, we clearly see that the HB properties of MFC1 align with those of the prototypical water dimer O-H···O HB. Each HB property generally becomes less prominent in the higher-energy conformers. Only the MFC1 conformer prominently exhibits (1) elongated C-H donor bonds, (2) attractive C-H···O═C interactions, (3) n(O) → σ*(C-H) hyperconjugation, (4) critical points in the electron density from Bader's method and from the noncovalent interactions method, (5) positively charged donor hydrogen, and (6) downfield NMR chemical shifts and nonzero 2J(CM-HM···OF) coupling constants. Based on this research, some issues merit further study. The flat PES hinders reliable determinations of the HB-induced shifts of the C-H stretches; a similarly difficult challenge is observed for the experiment. The role of charge transfer in HBs remains an intriguing open question, although our BLW and NBO computations suggest that it is relevant to the C-H···O HB geometries. These issues notwithstanding, the prominence of the HB properties in MFC1 serves as clear evidence that the MFC is predominantly bound by a C-H···O HB.

A Dynamic Equilibrium of Three Hydrogen-Bond Conformers Explains the NMR Spectrum of the Active Site of Photoactive Yellow Protein.

A theoretical study on the NMR shifts of the hydrogen bond network around the chromophore, para-coumaric acid (pCA), of photoactive yellow protein (PYP) is presented. Previous discrepancies between theoretical and experimental studies are resolved by our findings of a previously unknown rapid conformational exchange near the active site of PYP. This exchange caused by the rotation of Thr50 takes place in the ground state of PYP's active site and results in three effectively energetically equal conformations characterized by the formation of new hydrogen bonds, all of which contribute to the overall NMR signals of the investigated protons. In light of these findings, we are able to successfully explain the experimental results and provide valuable insight into the behavior of PYP in solution. We further investigated related PYP mutants (T50V, E46Q, and Y42F), and found the same conformational exchange in E46Q and Y42F to be responsible for the experimentally observed NMR and UV/vis spectra.

Unraveling the Base Excision Repair Mechanism of Human DNA Glycosylase.

Human DNA glycosylase, hOGG1, is known to perform DNA repair by cleaving oxidized guanine (8OG) from the DNA. Despite numerous experimental and theoretical investigations, the underlying selective molecular mechanism has remained a mystery. Here we present a mechanism that explains how hOGG1's catalytic pocket is able to host an undamaged guanine base, but is not able to cleave it from the DNA. Using linear-scaling quantum mechanics/molecular mechanics (QM/MM) techniques with more than 500 atoms in the QM-region, we have investigated previously proposed mechanisms that all rely on protonating the 8OG nucleobase. We have found that the repair mechanisms propagated in the literature to this date are not capable of differentiating between the G and 8OG nucleobase. Besides this nonselectivity, they also involve reaction barriers that are too high, hence rendering the corresponding reaction intermediates inaccessible. Instead, we present a completely different reaction mechanism, where hOGG1 initially targets the ribose moiety of the substrate and cleaves the glycosidic bond at the very last stage. Our ribose-protonated repair mechanism is not only energetically more preferable, but also explains the selectivity utilized by hOGG1 to block processing a guanine base.

A base-independent repair mechanism for DNA glycosylase--no discrimination within the active site.

The ubiquitous occurrence of DNA damages renders its repair machinery a crucial requirement for the genomic stability and the survival of living organisms. Deficiencies in DNA repair can lead to carcinogenesis, Alzheimer, or Diabetes II, where increased amounts of oxidized DNA bases have been found in patients. Despite the highest mutation frequency among oxidized DNA bases, the base-excision repair process of oxidized and ring-opened guanine, FapydG (2,6-diamino-4-hydroxy-5-formamidopyrimidine), remained unclear since it is difficult to study experimentally. We use newly-developed linear-scaling quantum-chemical methods (QM) allowing us to include up to 700 QM-atoms and achieving size convergence. Instead of the widely assumed base-protonated pathway we find a ribose-protonated repair mechanism which explains experimental observations and shows strong evidence for a base-independent repair process. Our results also imply that discrimination must occur during recognition, prior to the binding within the active site.

Effect of including torsional parameters for histidine-metal interactions in classical force fields for metalloproteins.

Classical force-field parameters of the metal site of metalloproteins usually comprise only the partial charges of the involved atoms, as well as the bond-stretching and bending parameters of the metal-ligand interactions. Although for certain metal ligands such as histidine residues, the torsional motions at the metal site play an important role for the dynamics of the protein, no such terms have been considered to be crucial in the parametrization of the force fields, and they have therefore been omitted in the parametrization. In this work, we have optimized AMBER-compatible force-field parameters for the reduced state of the metal site of copper, zinc superoxide dismutase (SOD1) and assessed the effect of including torsional parameters for the histidine-metal interactions in molecular dynamics simulations. On the basis of the obtained results, we recommend that torsion parameters of the metal site are included when processes at the metal site are investigated or when free-energy calculations are performed. As the torsion parameters mainly affect the structure of the metal site, other kinds of structural studies can be performed without considering the torsional parameters of the metal site.

Ribose-protonated DNA base excision repair: a combined theoretical and experimental study.

Living organisms protect the genome against external influences by recognizing and repairing damaged DNA. A common source of gene mutation is the oxidized guanine, which undergoes base excision repair through cleavage of the glycosidic bond between the ribose and the nucleobase of the lesion. We unravel the repair mechanism utilized by bacterial glycosylase, MutM, using quantum-chemical calculations involving more than 1000 atoms of the catalytic site. In contrast to the base-protonated pathway currently favored in the literature, we show that the initial protonation of the lesion's ribose paves the way for an almost barrier-free glycosidic cleavage. The combination of theoretical and experimental data provides further insight into the selectivity and discrimination of MutM's binding site toward various substrates.

Deamination, oxidation, and C-C bond cleavage reactivity of 5-hydroxymethylcytosine, 5-formylcytosine, and 5-carboxycytosine.

Three new cytosine derived DNA modifications, 5-hydroxymethyl-2'-deoxycytidine (hmdC), 5-formyl-2'-deoxycytidine (fdC) and 5-carboxy-2'-deoxycytidine (cadC) were recently discovered in mammalian DNA, particularly in stem cell DNA. Their function is currently not clear, but it is assumed that in stem cells they might be intermediates of an active demethylation process. This process may involve base excision repair, C-C bond cleaving reactions or deamination of hmdC to 5-hydroxymethyl-2'-deoxyuridine (hmdU). Here we report chemical studies that enlighten the chemical reactivity of the new cytosine nucleobases. We investigated their sensitivity toward oxidation and deamination and we studied the C-C bond cleaving reactivity of hmdC, fdC, and cadC in the absence and presence of thiols as biologically relevant (organo)catalysts. We show that hmdC is in comparison to mdC rapidly oxidized to fdC already in the presence of air. In contrast, deamination reactions were found to occur only to a minor extent. The C-C bond cleavage reactions require the presence of high concentration of thiols and are acid catalyzed. While hmdC dehydroxymethylates very slowly, fdC and especially cadC react considerably faster to dC. Thiols are active site residues in many DNA modifiying enzymes indicating that such enzymes could play a role in an alternative active DNA demethylation mechanism via deformylation of fdC or decarboxylation of cadC. Quantum-chemical calculations support the catalytic influence of a thiol on the C-C bond cleavage.

Why BLUF photoreceptors with roseoflavin cofactors lose their biological functionality.

The photophysics of roseoflavin in three different environments is investigated by using ab initio and quantum mechanics/molecular mechanics methods. Intramolecular charge transfer is shown to be responsible for the quenching of the fluorescence in the gas phase, and in the water environment. However, for the roseoflavin incorporated into the blue light using flavin (BLUF) protein environment (substituting the native flavin) no such deactivation is found. The conical intersection between the locally excited state of the chromophore and the charge transfer state involving the tyrosine residue, which in the native BLUF domain is responsible for initiating the photocycle, is missing for the roseoflavin substituted protein. This explains the experimental observations of the lack of any photocycle, and the loss of the biological function of the BLUF photoreceptor reported earlier.

Intermolecular interactions in photodamaged DNA from density functional theory symmetry-adapted perturbation theory.

The intermolecular interactions of the photodamaged cyclobutane pyrimidine dimer (CPD) lesion with adjacent nucleobases in the native intrahelical DNA double strand are investigated at the level of density functional theory symmetry-adapted perturbation theory (DFT-SAPT) and compared to the original (or repaired) case with pyrimidines (TpT) instead of CPD. The CPD aggregation is on average destabilized by about 6 kcal mol(-1) relative to that involving TpT. The effect of destabilization is asymmetric, that is, it involves a single H-bonding (Watson-Crick (WC) type) base-pair interaction.

Theoretical study on the repair mechanism of the (6-4) photolesion by the (6-4) photolyase.

UV irradiation of DNA can lead to the formation of mutagenic (6-4) pyrimidine-pyrimidone photolesions. The (6-4) photolyases are the enzymes responsible for the photoinduced repair of such lesions. On the basis of the recently published crystal structure of the (6-4) photolyase bound to DNA [Maul et al. 2008] and employing quantum mechanics/molecular mechanics techniques, a repair mechanism is proposed, which involves two photoexcitations. The flavin chromophore, initially being in its reduced anionic form, is photoexcited and donates an electron to the (6-4) form of the photolesion. The photolesion is then protonated by the neighboring histidine residue and forms a radical intermediate. The latter undergoes a series of energy stabilizing hydrogen-bonding rearrangements before the electron back transfer to the flavin semiquinone. The resulting structure corresponds to the oxetane intermediate, long thought to be formed upon DNA-enzyme binding. A second photoexcitation of the flavin promotes another electron transfer to the oxetane. Proton donation from the same histidine residue allows for the splitting of the four-membered ring, hence opening an efficient pathway to the final repaired form. The repair of the lesion by a single photoexcitation was shown not to be feasible.

A QM/MM study on the fast photocycle of blue light using flavin photoreceptors in their light-adapted/active form.

Based on QM/MM calculations using a combination of time-dependent Hartree-Fock and coupled cluster response theory a mechanism is proposed for the photocycle of blue light using flavin (BLUF) domains in the signaling/light adapted conformation. In analogy to the dark-adapted form, a charge transfer state drives proton transfer from the highly conserved tyrosine residue to the flavin chromophore. The latter step is mediated by the adjacent glutamine residue, which, in the light adapted conformation, exists as its imidic tautomer. However, before the proton transfer is even halfway completed, a conical intersection seam between the charge transfer and ground state is reached. Two channels for the decay back to the initial light-adapted conformation are open, a rapid one leading directly through the funnel of the conical intersection, bypassing the formation of the biradical intermediate, and a slower one via the biradical intermediate. The mechanism as proposed here: (i) explains the very rapid photocycle; and (ii) confirms the concept of photoirreversibility, both of which have been experimentally observed for BLUF domains in their light-adapted conformations.

A conclusive mechanism of the photoinduced reaction cascade in blue light using flavin photoreceptors.

On the basis of extensive first-principle calculations within the framework of quantum mechanics/molecular mechanics (QM/MM), a conclusive mechanism for the formation of the signaling state of blue light using flavin (BLUF) domain proteins is proposed which is compatible with the experimental data presently available. Time-dependent density functional, as well as advanced coupled cluster response theory was employed for the QM part in order to describe the relevant excited states. One of the key residues involved in the mechanism is the glutamine adjacent to the flavin chromophore. The reaction cascade, triggered by the initial photoexcitation of the flavin chromophore, involves isomerization of this residue but no rotation as assumed previously. The fact that only the environment, but not the flavin chromophore by itself, is chemically transformed along the individual steps of the mechanism is unique for biological photoreceptors. The final isomer of the glutamine tautomer, i.e., the imidic acid, is further stabilized by the interchange of a methionine residue in the binding pocket with a tryptophan residue. The flip of these two residues might be the trigger for the large conformational change of this protein which is consequently transmitted as the signal to the biological environment.

On the photophysics of artificial blue-light photoreceptors: an ab initio study on a flavin-based dye dyad at the level of coupled-cluster response theory.

The photophysical behavior of a phenothiazine-phenyl-isoalloxazine dye dyad, a model system for blue-light photoreceptors functioning on the basis of photoinduced electron transfer, was investigated by employing a combination of time-dependent density functional and coupled-cluster response theory. A conical intersection between a "bright" locally excited and a "dark" charge-transfer state was found in the low-energy region of the corresponding potential energy surfaces. We propose that, for the solvated dyad, this conical intersection is responsible for the experimentally observed fast fluorescence quenching in that system.