Maciej Haranczyk

Electrostatic basis for the unidirectionality of the primary proton transfer in cytochrome c oxidase

Gaining detailed understanding of the energetics of the proton-pumping process in cytochrome c oxidase (CcO) is one of the challenges of modern biophysics. Despite promising mechanistic proposals, most works have not related the activation barriers of the different assumed steps to the protein structure, and there has not been a physically consistent model that reproduced the barriers needed to create a working pump. This work reevaluates the activation barriers for the primary proton transfer (PT) steps by calculations that reflect all relevant free energy contributions, including the electrostatic energies of the generated charges, the energies of water insertion, and large structural rearrangements of the donor and acceptor. The calculations have reproduced barriers that account for the directionality and sequence of events in the primary PT in CcO. It has also been found that the PT from Glu-286 (E) to the propionate of heme a3 (Prd) provides a gate for an initial back leakage from the high pH side of the membrane. Interestingly, the rotation of E that brings it closer to Prd appears to provide a way for blocking competing pathways in the primary PT. Our study elucidates and quantifies the nature of the control of the directionality in the primary PT in CcO and provides instructive insight into the role of the water molecules in biological PT, showing that “bridges” of several water molecules in hydrophobic regions present a problem (rather than a solution) that is minimized in the primary PT.

Comparison of some representative density functional theory and wave function theory methods for the studies of amino acids

Energies of different conformers of 22 amino acid molecules and their protonated and deprotonated species were calculated by some density functional theory (DFT; SVWN, B3LYP, B3PW91, MPWB1K, BHandHLYP) and wave function theory (WFT; HF, MP2) methods with the 6‐311++G(d,p) basis set to obtain the relative conformer energies, vertical electron detachment energies, deprotonation energies, and proton affinities. Taking the CCSD/6‐311++G(d,p) results as the references, the performances of the tested DFT and WFT methods for amino acids with various intramolecular hydrogen bonds were determined. The BHandHLYP method was the best overall performer among the tested DFT methods, and its accuracy was even better than that of the more expensive MP2 method. The computational dependencies of the five DFT methods and the HF and MP2 methods on the basis sets were further examined with the 6‐31G(d,p), 6‐311++G(d,p), aug‐cc‐pVDZ, 6‐311++G(2df,p), and aug‐cc‐pVTZ basis sets. The differences between the small and large basis set results have decreased quickly for the hybrid generalized gradient approximation (GGA) methods. The basis set convergence of the MP2 results has been, however, very slow. Considering both the cost and the accuracy, the BHandHLYP functional with the 6‐311++G(d,p) basis set is the best choice for the amino acid systems that are rich in hydrogen bonds.

Electron-Driven Acid-Base Chemistry : Proton Transfer from Hydrogen Chloride to Ammonia

In contrast to widely familiar acid-base behavior in solution, single molecules of NH3 and HCl do not react to form the ionic salt, NH+4Cl–, in isolation. We applied anion photoelectron spectroscopy and ab initio theory to investigate the interaction of an excess electron with the hydrogen-bonded complex NH3···HCl. Our results show that an excess electron induces this complex to form the ionic salt. We propose a mechanism that proceeds through a dipole-bound state to form the negative ion of ionic ammonium chloride, a species that can also be characterized as a deformed Rydberg radical, NH4, polarized by a chloride anion, Cl–.

Catalysis models for the enolization of the acetaldehyde radical cation

Abstract The catalysis of the enolization of the prototype aldehyde radical cation in the gas phase using methanol or acetaldehyde as the catalyst was studied theoretically. For the computations the CBS-Q model was used, in combination with RHF/DZP optimized geometries. Various reaction path models were tested, including the proton-transport catalysis (PTC) model, the Spectator model and the Quid pro Quo (QpQ) model. For the latter model there are several possibilities, depending on the choice made for the catalyst hydrogen to be exchanged. For some reaction steps no transition states exist. Instead, there are intersections of two potential energy surfaces corresponding to different localizations of the radical site in the ion–molecule complex. For acetaldehyde both the PTC and QpQ (all variants) models are shown to be feasible. However, the barrier heights involved in these models are very different and the corresponding reaction rates are also expected to differ widely. As a consequence different reaction pathways may be applicable for reactions taking place in different time frames.

Solvation free energies of molecules. The most stable anionic tautomers of uracil

Anionic states of nucleic acid bases are suspected to play a role in the radiation damage processes of DNA. Our recent studies suggested that the excess electron attachment to the nucleic acid bases can stabilize some rare tautomers, i.e. imine-enamine tautomers and other tautomers with a proton being transferred from nitrogen sites to carbon sites (with respect to the canonical tautomer). So far, these new anionic tautomers have been characterized by the gas-phase electronic structure calculations and photoelectron spectroscopy experiments. In the current contribution we explore the effect of water solvation on the stability of the new anionic tautomers of uracil. The accurate free energies of solvation are calculated in a two step approach. The major contribution was calculated using the classical free-energy perturbation adiabatic-charging approach, where it is assumed that the solvated molecule has the charge distribution given by the polarizable continuum model. In the second step the free energy of solvation is refined by taking into account the real, average solvent charge distribution. This is done using our accelerated QM/MM simulations, where the QM energy of the solute is calculated in the mean potential averaged over many MD steps. We found that in water solution three of the recently identified anionic tautomers are 6.5-3.6 kcal mol(-1) more stable than the anion of the canonical tautomer.

STABLE VALENCE ANIONS OF NUCLEIC ACID BASES AND DNA STRAND BREAKS INDUCED BY LOW ENERGY ELECTRONS

The last decade has witnessed immense advances in our understanding of the effects of ionizing radiation on biological systems. As the genetic information carrier in biological systems, DNA is the most important species which is prone to damage by high energy photons. Ionizing radiations destroy DNA indirectly by forming low energy electrons (LEEs) as secondary products of the interaction between ionizing radiation and water. An understanding of the mechanism that leads to the formation of single and double strand breaks may be important in guiding the further development of anticancer radiation therapy. In this article we demonstrate the likely involvement of stable nucleobases anions in the formation of DNA strand breaks – a concept which the radiation research community has not focused on so far. In Section refch21:sec21.1 we discuss the current status of studies related to the interaction between DNA and LEEs. The next section is devoted to the description of proton transfer induced by electron attachment to the complexes between nucleobases and various proton donors – a process leading to the strong stabilization of nucleobases anions. Then, we review our results concerning the anionic binary complexes of nucleobases with particular emphasize on the GC and AT systems. Next, the possible consequences of interactions between DNA and proteins in the context of electron attachment are briefly discussed. Further, we focus on existing proposal of single strand break formation in DNA. Ultimately, open questions as well perspectives of studies on electron induced DNA damage are discussed

Discovery of Most Stable Structures of Neutral and Anionic Phenylalanine through Automated Scanning of Tautomeric and Conformational Spaces.

We have developed a software tool for combinatorial generation of tautomers and conformers of small molecules. We have demonstrated it by performing a systematic search for the most stable structures of neutral and anionic phenylalanine (Phe) using electronic structure methods. For the neutral canonical tautomer we found out that the conformers with and without the intramolecular (O)H···NH2 hydrogen bond are similarly stable, within the error bars of our method. A unique IR signature of the conformer without the hydrogen bond has been identified. We also considered anions of Phe, both valence type and dipole-bound. We have found out that tautomers resulting from proton transfer from the carboxylic OH to the phenyl ring do support valence anions that are vertically strongly bound, with electron vertical detachment energies (VDE) in a range of 3.2-3.5 eV. The most stable conformer of these valence anions remains adiabatically unbound with respect to the canonical neutral by only 2.17 kcal/mol at the CCSD(T)/aug-cc-pVDZ level. On the basis of our past experience with valence anions of nucleic acid bases, we suggest that the valence anions of Phe identified in this report can be observed experimentally. The most stable conformer of canonical Phe is characterized by an adiabatic electron affinity of 53 meV (a dipole-bound state).

Structure and Energetics of Clustered Damage Sites

Abstract Miller, J. H., Aceves-Gaona, A., Ernst, M. B., Haranczyk, M., Gutowski, M., Vorpagel, E. R. and Dupuis, M. Structure and Energetics of Clustered Damage Sites. Radiat. Res. 164, 582–585 (2005). Quantum calculations on duplex DNA trimers were used to model the changes in structure, hydrogen bonding, stacking properties, and electrostatic potential induced by oxidized purine bases and abasic (AP) sites. Results for oxidized purine bases were consistent with experimental data that show small structural and energetic perturbations induced by isolated 8-oxoguanine (8oG). Watson-Crick base pairing was preserved, and no major distortions of the backbone were induced. The thermal destabilization of DNA induced by 8oG was comparable to the energy of a single hydrogen bond. In contrast, AP sites caused substantial distortions of the DNA backbone that were accompanied by relocation of counterions. The loss of Watson-Crick hydrogen bonds in AP sites had the potential to destabilize DNA by 10–20 kcal/mol (0.4–0.8 eV); however, new inter- and intrastrand hydrogen bonds formed after removal of a nucleic acid base that significantly affected the energy of AP sites and introduced a strong dependence on sequence context. Quantum calculations on small DNA fragments provided starting conformations and force-field parameters for classical molecular dynamics simulations of radiation-induced single-strand breaks that most often combine hydrolysis of a phosphate-oxygen (P-O) bond with an AP site and fully or partially degraded sugar ring. P-O bond hydrolysis increased the freedom in backbone torsion angles, which allowed the broken strand to compress and partially fill the hole in the DNA created by the AP site. Results for strand breaks with a 3′phosphoglycolate were similar to those with phosphate end groups.

Structure and singly occupied molecular orbital analysis of anionic tautomers of guanine

Recently, we reported the discovery of adiabatically bound anions of guanine that might be involved in the processes of DNA damage by low‐energy electrons and in charge transfer through DNA. These anions correspond to some tautomers that have been ignored thus far. They were identified using a hybrid quantum mechanical–combinatorial approach in which an energy‐based screening was performed on the library of 499 tautomers with their relative energies calculated with quantum chemistry methods. In the current study, we analyze the adiabatically bound anions of guanine in two aspects: (1) the geometries and excess electron distributions are analyzed and compared with anions of the most stable neutrals to identify the sources of stability; (2) the chemical space of guanine tautomers is explored to verify if these new tautomers are contained in a particular subspace of the tautomeric space. The first task involves the development of novel approaches—the quantum chemical data like electron density, orbital, and information on its bonding/antibonding character are coded into holograms and analyzed using chemoinformatics techniques. The second task is completed using substructure analysis and clustering techniques performed on molecules represented by 2D fingerprints. The major conclusion is that the high stability of adiabatically bound anions originates from the bonding character of the π orbital occupied by the excess electron. This compensates for the antibonding character that usually causes significant buckling of the ring. Also, the excess electron is more homogenously distributed over both rings than in the case of anions of the most stable neutral species. In terms of 2D substructure, the most stable anionic tautomers generally have additional hydrogen atoms at C8 and/or C2 and they do not have hydrogen atoms attached to C4, C5, and C6. They also form an “island of stability” in the tautomeric space of guanine.