Matthew D. Liptak

Absolute pKa Determinations for Substituted Phenols

The CBS-QB3 method was used to calculate the gas-phase free energy difference between 20 phenols and their respective anions, and the CPCM continuum solvation method was applied to calculate the free energy differences of solvation for the phenols and their anions. The CPCM solvation calculations were performed on both gas-phase and solvent-phase optimized structures. Absolute pK(a) calculations with solvated phase optimized structures for the CPCM calculations yielded standard deviations and root-mean-square errors of less than 0.4 pK(a) unit. This study is the most accurate absolute determination of the pK(a) values of phenols, and is among the most accurate of any such calculations for any group of compounds. The ability to make accurate predictions of pK(a) values using a coherent, well-defined approach, without external approximations or fitting to experimental data, is of general importance to the chemical community. The solvated phase optimized structures of the anions are absolutely critical to obtain this level of accuracy, and yield a more realistic charge separation between the negatively charged oxygen and the ring system of the phenoxide anions.

Accurate relative pKa calculations for carboxylic acids using complete basis set and Gaussian-n models combined with continuum solvation methods

The complete basis set methods CBS-4, CBS-QB3, and CBS-APNO, and the Gaussian methods G2 and G3 were used to calculate the gas phase energy differences between six different carboxylic acids and their respective anions. Two different continuum methods, SM5.42R and CPCM, were used to calculate the free energy differences of solvation for the acids and their anions. Relative pKa values were calculated for each acid using one of the acids as a reference point. The CBS-QB3 and CBS-APNO gas phase calculations, combined with the CPCM/HF/6-31+G(d)//HF/6-31G(d) or CPCM/HF/6-31+G(d)//HF/6-31+G(d) continuum solvation calculations on the lowest energy gas phase conformer, and with the conformationally averaged values, give results accurate to 12 pKa unit.

Suppression of Kasha's rule as a mechanism for fluorescent molecular rotors and aggregation-induced emission

Although there are some proposed explanations for aggregation-induced emission, a phenomenon with applications that range from biosensors to organic light-emitting diodes, current understanding of the quantum-mechanical origin of this photophysical behaviour is limited. To address this issue, we assessed the emission properties of a series of BF2–hydrazone-based dyes as a function of solvent viscosity. These molecules turned out to be highly efficient fluorescent molecular rotors. This property, in addition to them being aggregation-induced emission luminogens, enabled us to probe deeper into their emission mechanism. Time-dependent density functional theory calculations and experimental results showed that the emission is not from the S1 state, as predicted from Kashas rule, but from a higher energy (>S1) state. Furthermore, we found that suppression of internal conversion to the dark S1 state by restricting the rotor rotation enhances fluorescence, which leads to the proposal that suppression of Kashas rule is the photophysical mechanism responsible for emission in both viscous solution and the solid state. A family of fluorescent molecular rotors has been developed and their mechanism for emission understood. It has been observed that, although most fluorescent molecules emit from their lowest energy excited state, S1 (in accordance with Kashas rule), BODIHY dyes do not. Furthermore, their fluorescence is enhanced through restricted rotor rotation, which suppresses internal conversion to the dark S1 state.

NMR and DFT Investigation of Heme Ruffling: Functional Implications for Cytochrome c

Out-of-plane (OOP) deformations of the heme cofactor are found in numerous heme-containing proteins and the type of deformation tends to be conserved within functionally related classes of heme proteins. We demonstrate correlations between the heme ruffling OOP deformation and the (13)C and (1)H nuclear magnetic resonance (NMR) hyperfine shifts of heme aided by density functional theory (DFT) calculations. The degree of ruffling in the heme cofactor of Hydrogenobacter thermophilus cytochrome c(552) has been modified by a single amino acid mutation in the second coordination sphere of the cofactor. The (13)C and (1)H resonances of the cofactor have been assigned using one- and two-dimensional NMR spectroscopy aided by selective (13)C-enrichment of the heme. DFT has been used to predict the NMR hyperfine shifts and electron paramagnetic resonance (EPR) g-tensor at several points along the ruffling deformation coordinate. The DFT-predicted NMR and EPR parameters agree with the experimental observations, confirming that an accurate theoretical model of the electronic structure and its response to ruffling has been established. As the degree of ruffling increases, the heme methyl (1)H resonances move upfield while the heme methyl and meso (13)C resonances move downfield. These changes are a consequence of altered overlap of the Fe 3d and porphyrin pi orbitals, which destabilizes all three occupied Fe 3d-based molecular orbitals and decreases the positive and negative spin density on the beta-pyrrole and meso carbons, respectively. Consequently, the heme ruffling deformation decreases the electronic coupling of the cofactor with external redox partners and lowers the reduction potential of heme.

The Mycobacterium tuberculosis secreted protein Rv0203 transfers heme to membrane proteins MmpL3 and MmpL11.

Background: A novel Mycobacterium tuberculosis heme acquisition has recently been discovered. Results: The membrane protein MmpL11 is required for efficient heme uptake, and secreted Rv0203 may transfer heme to extracellular domains of both MmpL3 and MmpL11. Conclusion: MmpL3 and MmpL11 are potential heme transporters, whereby heme is transported into the cytosol. Significance: This work enhances our understanding of Mycobacterium tuberculosis heme uptake. Mycobacterium tuberculosis is the causative agent of tuberculosis, which is becoming an increasingly global public health problem due to the rise of drug-resistant strains. While residing in the human host, M. tuberculosis needs to acquire iron for its survival. M. tuberculosis has two iron uptake mechanisms, one that utilizes non-heme iron and another that taps into the vast host heme-iron pool. To date, proteins known to be involved in mycobacterial heme uptake are Rv0203, MmpL3, and MmpL11. Whereas Rv0203 transports heme across the bacterial periplasm or scavenges heme from host heme proteins, MmpL3 and MmpL11 are thought to transport heme across the membrane. In this work, we characterize the heme-binding properties of the predicted extracellular soluble E1 domains of both MmpL3 and MmpL11 utilizing absorption, electron paramagnetic resonance, and magnetic circular dichroism spectroscopic methods. Furthermore, we demonstrate that Rv0203 transfers heme to both MmpL3-E1 and MmpL11-E1 domains at a rate faster than passive heme dissociation from Rv0203. This work elucidates a key step in the mycobacterial uptake of heme, and it may be useful in the development of anti-tuberculosis drugs targeting this pathway.

The Proapoptotic G41S Mutation to Human Cytochrome c Alters the Heme Electronic Structure and Increases the Electron Self-Exchange Rate.

The naturally occurring G41S mutation to human (Hs) cytochrome (cyt) c enhances apoptotic activity based upon previous in vitro and in vivo studies, but the molecular mechanism underlying this enhancement remains unknown. Here, X-ray crystallography, nuclear magnetic resonance (NMR) spectroscopy, and density functional theory (DFT) calculations have been used to identify the structural and electronic differences between wild-type (WT) and G41S Hs cyt c. S41 is part of the hydrogen bonding network for propionate 7 of heme pyrrole ring A in the X-ray structure of G41S Hs cyt c and, compared to WT, G41S Hs cyt c has increased spin density on pyrrole ring C and a faster electron self-exchange rate. DFT calculations illustrate an electronic mechanism where structural changes near ring A can result in electronic changes at ring C. Since ring C is part of the solvent-exposed protein surface, we propose that this heme electronic structure change may ultimately be responsible for the enhanced proapoptotic activity of G41S Hs cyt c.

Comparison of CBS-QB3, CBS-APNO, G2, and G3 thermochemical predictions with experiment for formation of ionic clusters of hydronium and hydroxide ions complexed with water

The GAUSSIAN 2, GAUSSIAN 3, complete basis set-QB3, and complete basis set-APNO methods have been used to calculate DeltaH( composite function) and DeltaG( composite function) values for ionic clusters of hydronium and hydroxide ions complexed with water. Results for the clusters H3O+(H2O)n and OH-(H2O)n, where n=1-4 are reported in this paper, and compared against experimental values contained in the National Institutes of Standards and Technology (NIST) database. Agreement with experiment is excellent for the three ab initio methods for formation of these clusters. The high accuracy of these methods makes them reliable for calculating energetics for the formation of ionic clusters containing water. In addition this allows them to serve as a valuable check on the accuracy of experimental data reported in the NIST database, and makes them useful tools for addressing unresolved issues in atmospheric chemistry.

Measurement of Heme Ruffling Changes in MhuD Using UV–vis Spectroscopy

For decades it has been known that an out-of-plane ruffling distortion of heme perturbs its UV-vis absorption (Abs) spectrum, but whether increased ruffling induces a red or blue shift of the Soret band has remained a topic of debate. This debate has been resolved by the spectroscopic and computational characterization of Mycobacterium tuberculosis MhuD presented here, an enzyme that converts heme, oxygen, and reducing equivalents to nonheme iron and mycobilin. W66F and W66A MhuD have been characterized using (1)H nuclear magnetic resonance, Abs, and magnetic circular dichroism spectroscopies, and the data have been used to develop an experimentally validated theoretical model of ruffled, ferric heme. The PBE density functional theory (DFT) model that has been developed accurately reproduces the observed spectral changes from wild type enzyme, and the underlying quantum mechanical origins of these ruffling-induced changes were revealed by analyzing the PBE DFT description of the electronic structure. Small amounts of heme ruffling have no influence on the energy of the Q-band and blue-shift the Soret band due to symmetry-allowed mixing of the Fe 3dxy and porphyrin a2u orbitals. Larger amounts of ruffling red-shift both the Q and Soret bands due to disruption of π-bonding within the porphyrin ring.

Hydrogen bond donation to the heme distal ligand of Staphylococcus aureus IsdG tunes the electronic structure

Staphylococcus aureus IsdG catalyzes the final step of staphylococcal iron acquisition from host hemoglobin, whereby host-derived heme is converted to iron and organic products. The Asn7 distal pocket residue is known to be critical for enzyme activity, but the influence of this residue on the substrate electronic structure was unknown prior to this work. Here, an optical spectroscopic and density functional theory characterization of azide- and cyanide-inhibited wild type and N7A IsdG is presented. Magnetic circular dichroism data demonstrate that Asn7 perturbs the electronic structure of azide-inhibited, but not cyanide-inhibited, IsdG. As the iron-ligating α-atom of azide, but not cyanide, can act as a hydrogen bond acceptor, these data indicate that the terminal amide of Asn7 is a hydrogen bond donor to the α-atom of a distal ligand to heme in IsdG. Circular dichroism characterization of azide- and cyanide-inhibited forms of WT and N7A IsdG strongly suggests that the Asn7···N3 hydrogen bond influences the orientation of a distal azide ligand with respect to the heme substrate. Specifically, density functional theory calculations suggest that Asn7···N3 hydrogen bond donation causes the azide ligand to rotate about an axis perpendicular to the porphyrin plane and weakens the π-donor strength of the azide ligand. This lowers the energies of the Fe 3dxz and 3dyz orbitals, mixes Fe 3dxy and porphyrin a2u character into the singly-occupied molecular orbital, and results in spin delocalization onto the heme meso carbons. These discoveries have important implications for the mechanism of heme oxygenation catalyzed by IsdG.

12:Computational Studies of Bioorganometallic Enzymes and Cofactors

Because of their complex geometric and electronic structures, the active sites and cofactors of bioorganometallic enzymes, which are characterized by their metal–carbon bonds, pose a major challenge for computational chemists. However, recent progress in computer technology and theoretical chemistry, along with insights gained from mechanistic, spectroscopic, and X-ray crystallographic studies, have established an excellent foundation for the successful completion of computational studies aimed at elucidating the electronic structures and catalytic cycles of these species. This chapter briefly reviews the most popular computational approaches employed in theoretical studies of bioorganometallic species and summarizes important information obtained from computational studies of (i) the enzymatic formation and cleavage of the Co–C bond of coenzyme B12; (ii) the catalytic cycle of methyl-coenzyme M reductase and its nickel-containing cofactor F430; (iii) the polynuclear active-site clusters of the bifunctional enzyme carbon monoxide dehydrogenase/acetyl-coenzyme A synthase; and (iv) the magnetic properties of the active-site cluster of Fe-only hydrogenases.