Improved Density Functional Description of the Electrochemistry and Structure−Property Descriptors of Substituted Flavins

Electronic Structure Methods for Simulating Flavin's Spectroscopy and Photophysics: Comparison of Multi-reference, TD-DFT, and Single-Reference Wave Function Methods

The use of flavins and flavoproteins in photocatalytic, sensing, and biotechnological applications has led to a growing interest in computationally modeling the excited-state electronic structure and photophysics of flavin. However, there is limited consensus regarding which computational methods are appropriate for modeling flavin's photophysics. We compare the energies of low-lying excited states of flavin computed with time-dependent density functional theory (TD-DFT), equation-of-motion coupled cluster (EOM-EE-CCSD), scaled opposite-spin configuration interaction [SOS-CIS(D)], multiconfiguration pair-density functional theory (MC-PDFT), and several multireference perturbation theory (MR-PT2) methods. In the first part, we focus on excitation energies of the first singlet excited state (S1) of five different redox and protonation states of flavin, with the goal of finding a suitable active space for MR-PT2 calculations. In the second part, we construct two sets of one-dimensional potential energy surfaces connecting the S0 and S1 equilibrium geometries (S0-S1 path) and the S1 (π,π*) and S2 (n,π*) equilibrium geometries (S1-S2 path). The first path therefore follows a Franck-Condon active mode of flavin while the second path maps crossings points between low-lying singlet and triplet states in flavin. We discuss the similarities and differences in the TD-DFT, EOM-EE-CCSD, SOS-CIS(D), MC-PDFT and MR-PT2 energy profiles along these paths. We find that (TD-)DFT methods are suitable for applications such as simulating the spectra of flavins but are inconsistent with several other methods when used for some geometry optimizations and when describing the energetics of dark (n,π*) states. MR-PT2 methods show promise for the simulation of flavin's low-lying excited states, but the selection of orbitals for the active space and the number of roots used for state averaging must be done carefully to avoid artifacts. Some properties, such as the intersystem crossing geometry and energy between the S1 (π,π*) and T2 (n,π*) states, may require additional benchmarking before they can be determined quantitatively.

Virtual Screening of Molecules via Neural Fingerprint-based Deep Learning Technique

A machine learning-based drug screening technique has been developed and optimized using convolutional neural network-derived fingerprints. The optimization of weights in the neural network-based fingerprinting technique was compared with fixed Morgan fingerprints in regard to binary classification on drug-target binding affinity. The assessment was carried out using six different target proteins using randomly chosen small molecules from the ZINC15 database for training. This new architecture proved to be more efficient in screening molecules that less favorably bind to specific targets and retaining molecules that favorably bind to it.

Scientific contribution

We have developed a new neural fingerprint-based screening model that has a significant ability to capture hits. Despite using a smaller dataset, this model is capable of mapping chemical space similar to other contemporary algorithms designed for molecular screening. The novelty of the present algorithm lies in the speed with which the models are trained and tuned before testing its predictive capabilities and hence is a significant step forward in the field of machine learning-embedded computational drug discovery.

Redox Properties of Flavin in BLUF and LOV Photoreceptor Proteins from Hybrid QM/MM Molecular Dynamics Simulation

Flavins play an important role in many oxidation and reduction processes in biological systems. For example, flavin adenine dinucleotide (FAD) and flavin mononucleotide (FMN) are common cofactors found in enzymatic proteins that use the special redox properties of these flavin molecules for their catalytic or photoactive functions. The redox potential of the flavin is strongly affected by its (protein) environment; however, the underlying molecular interactions of this effect are still unknown. Using hybrid quantum mechanics/molecular mechanics (QM/MM) simulation techniques, we have studied the redox properties of flavin in the gas phase, aqueous solution, and two different protein environments, in particular, a BLUF and a LOV photoreceptor domain. By mapping the changes in electrostatic potential and solvent structure, we gain insight into how specific polarization of the flavin by its environment tunes the reduction potential. We find also that accurate calculation of the reduction potentials of these systems by using the hybrid QM/MM approach is hampered by a too limited sampling of the counterion configurations and by artifacts at the QM/MM boundary. We make suggestions for how these issues can be overcome.

Assessing the Performance of Non-Equilibrium Thermodynamic Integration in Flavodoxin Redox Potential Estimation

Flavodoxins are enzymes that contain the redox-active flavin mononucleotide (FMN) cofactor and play a crucial role in numerous biological processes, including energy conversion and electron transfer. Since the redox characteristics of flavodoxins are significantly impacted by the molecular environment of the FMN cofactor, the evaluation of the interplay between the redox properties of the flavin cofactor and its molecular surroundings in flavoproteins is a critical area of investigation for both fundamental research and technological advancements, as the electrochemical tuning of flavoproteins is necessary for optimal interaction with redox acceptor or donor molecules. In order to facilitate the rational design of biomolecular devices, it is imperative to have access to computational tools that can accurately predict the redox potential of both natural and artificial flavoproteins. In this study, we have investigated the feasibility of using non-equilibrium thermodynamic integration protocols to reliably predict the redox potential of flavodoxins. Using as a test set the wild-type flavodoxin from Clostridium Beijerinckii and eight experimentally characterized single-point mutants, we have computed their redox potential. Our results show that 75% (6 out of 8) of the calculated reaction free energies are within 1 kcal/mol of the experimental values, and none exceed an error of 2 kcal/mol, confirming that non-equilibrium thermodynamic integration is a trustworthy tool for the quantitative estimation of the redox potential of this biologically and technologically significant class of enzymes.

Polyethylene Glycol 20k. Does It Fluoresce?

Polyethylene glycol (PEG) is a polyether compound commonly used in biological research and medicine because it is biologically inert. This simple polymer exists in variable chain lengths (and molecular weights). As they are devoid of any contiguous π-system, PEGs are expected to lack fluorescence properties. However, recent studies suggested the occurrence of fluorescence properties in non-traditional fluorophores like PEGs. Herein, a thorough investigation has been conducted to explore if PEG 20k fluoresces. Results of this combined experimental and computational study suggested that although PEG 20k could exhibit "through-space" delocalization of lone pairs of electrons in aggregates/clusters, formed via intermolecular and intramolecular interactions, the actual contributor of fluorescence between 300 and 400 nm is the stabilizer molecule, i.e., 3-tert-butyl-4-hydroxyanisole present in the commercially available PEG 20k. Therefore, the reported fluorescence properties of PEG should be taken with a grain of salt, warranting further investigation.

Machine Learning for Efficient Prediction of Protein Redox Potential: The Flavoproteins Case

Determining the redox potentials of protein cofactors and how they are influenced by their molecular neighborhoods is essential for basic research and many biotechnological applications, from biosensors and biocatalysis to bioremediation and bioelectronics. The laborious determination of redox potential with current experimental technologies pushes forward the need for computational approaches that can reliably predict it. Although current computational approaches based on quantum and molecular mechanics are accurate, their large computational costs hinder their usage. In this work, we explored the possibility of using more efficient QSPR models based on machine learning (ML) for the prediction of protein redox potential, as an alternative to classical approaches. As a proof of concept, we focused on flavoproteins, one of the most important families of enzymes directly involved in redox processes. To train and test different ML models, we retrieved a dataset of flavoproteins with a known midpoint redox potential (Em) and 3D structure. The features of interest, accounting for both short- and long-range effects of the protein matrix on the flavin cofactor, have been automatically extracted from each protein PDB file. Our best ML model (XGB) has a performance error below 1 kcal/mol (∼36 mV), comparing favorably to more sophisticated computational approaches. We also provided indications on the features that mostly affect the Em value, and when possible, we rationalized them on the basis of previous studies.

Molecular properties and tautomeric equilibria of isolated flavins

Flavins are employed as redox cofactors and chromophores in a plethora of flavoenzymes. Their versatility is an outcome of intrinsic molecular properties of the isoalloxazine ring modulated by the protein scaffold and surrounding solvent. Thus, an investigation of isolated flavins with high-level electronic-structure methods and with error assessment of the calculated properties will contribute to building better models of flavin reactivity. Here, we benchmarked ground-state properties such as electron affinity, gas-phase basicity, dipole moment, torsion energy, and tautomer stability for lumiflavins in all biologically relevant oxidation and charge states. Overall, multiconfigurational effects are small and chemical accuracy is achieved by coupled-cluster treatments of energetic properties. Augmented basis sets and extrapolations to the complete basis-set limit are necessary for consistent agreement with experimental energetics. Among DFT functionals tested, M06-2X shows the best performance for most properties, except gas-phase basicity, in which M06 and CAM-B3LYP perform better. Moreover, dipole moments of radical flavins show large deviations for all functionals studied. Tautomers with noncanonical protonation states are significantly populated at normal temperatures, adding to the complexity of modeling flavins. These results will guide future computational studies of flavoproteins and flavin chemistry by indicating the limitations of electronic-structure methodologies and the contributions of multiple tautomeric states.

Using structure-function relationships to understand the mechanism of phenazine-mediated extracellular electron transfer in Escherichia coli

Phenazines are redox-active nitrogen-containing heterocyclic compounds that can be produced by either bacteria or synthetic approaches. As an electron shuttles (mediators), phenazines are involved in several biological processes facilitating extracellular electron transfer (EET). Therefore, it is of great importance to understand the structural and electronic properties of phenazines that promote EET in microbial electrochemical systems. Our previous study experimentally investigated a phenazine-based library as an exogenous mediator system to facilitate EET in Escherichia coli. Herein, we combine our experimental data with density functional theory (DFT) calculations and multivariate linear regression modeling to understand the structure-function relationships in phenazine-based mediated EET. These calculations demonstrate that the computed redox properties of phenazines in lipophilic environments (e.g., cell membrane) correlate to experimental mediated current densities. Additional DFT-derived molecular properties were considered to develop a predictive model, which could be used in metabolic engineering approaches to introduce phenazines as endogenous mediators into bacteria.

Spatial effects of photosensitization on morphology of giant unilamellar vesicles

Singlet oxygen (1O2) formed through photosensitization may initiate oxidative destruction of biomembranes, however, the influence from the spatial organization of photosensitizers (PS) relative to membranes remains unclear. To clarify this issue, we loaded riboflavin 5'-(dihydrogen phosphate) monosodium (FMN-Na) as a hydrophilic PS into the lumen of halloysite nanotubes (HNTs), and attached the nanoassemblies (FMN-Na@HNTs), via Pickering effects, to the outer surfaces of giant unilamellar vesicles (GUVs) of phospholipids. We also prepared GUVs dopped with lumiflavin (LF) as a lipophilic PS having a 1O2 quantum yield comparable to FMN-Na. FMN-Na capsulated in HNT was characterized by a longer triplet excited state lifetime (12.1 μs) compared to FMN-Na free in solution (7.5 μs), and FMN-Na in both forms efficiently generated 1O2 upon illumination. The spatio-effects of PS on the photosensitized morphological changes of membranes were studied using conventional optical microscopy by monitoring GUV morphological changes. Upon light exposure (400-440 nm), the GUVs attached with FMN-Na@HNT merely experienced membrane deformation starting from the original spherical shape, ascribed to Type II photosensitization with 1O2 as oxidant. In contrast, photooxidation of LF dopped GUVs mainly led to membrane coarsening and budding assigned to Type I photosensitization. The spatial effects of PS on photosensitized morphological changes were related to the different lipid oxidation products generated through Type I and Type II photosensitized lipid oxidation.

The redox potential of flavin derivatives as a mediator in biosensors

The two-electron reduction potential for a set of 393 flavin derivatives is presented in this article. These derivatives are substituted flavin on carbon 6, 7, 8, and 9 by coinage transition metals (Cu, Ag, and Au) and conjugated double bond hydrocarbons; and both groups are examined with and without functional groups such as OH, Cl, CH₃, COOH, and NO₂. In order to show the validity of the results, the reduction potential of human life molecules, which have experimental values, such as flavin adenine dinucleotide (FAD) and riboflavin (vitamin B₂) is calculated. The experimental value for FAD is - 0.22 V, while the obtained theoretical value is - 0.21 V, and the corresponding values for riboflavin are - 0.18 and - 0.19 V, respectively. Theoretical calculations have been carried out by DFT procedure with a 6-31+G** basis set and BLYP xc-functional for coinage transition metals substitution, and MPW1PW9 xc-functionals for conjugated double bond hydrocarbon substitution. Both xc-functionals are chosen by the DFT calibration procedure.

Investigating the Photochemistry of C7 and C8 Functionalized N(5)-Ethyl-flavinium Cation: A Computational Study

Flavins are a diverse set of compounds with a wide variety of biological and nonbiological applications. Applications of flavins receiving attention recently consist of electro- and photocatalytic oxidation of substrates for organic synthesis, bioengineered nanotechnology, and water splitting catalysts, among others. While there is vast knowledge regarding the structure-property relationships of flavins and their electrochemistry, there is much less work elucidating the structure property relationships as they pertain to flavinium photochemistry. Herein, we report the effect of molecular tailoring on the molecular properties of N(5)-ethyl-flavinium cation (Et-Fl⁺), a derivative of the biocatalytic coenzyme riboflavin, by incorporating electron withdrawing and donating groups at the C7 and C8 position of the isoalloxazine ring. The presence of electron withdrawing groups at the C8 position caused a red shift in the absorption spectrum, while the electron donating groups caused a blue shift. Functionalization at the C7 position had the opposite effect on the absorption spectrum. The effects of single substitution were relatively negated with simultaneous functionalization at both the C8 and C7 positions. Difference density plots indicate no change in the nature of the S₁ excited state, which was confirmed by optimization of the excited state geometries. The results presented in this study indicate that functionalization of the isoalloxazine unit affects the photophysical properties of N(5)-ethyl-flavinium cations.

Microscopic Picture of the Solvent Reorganization During Electron Transfer to Flavin in Water

The redox potential of molecular species is largely modulated by its molecular environment so that a change of the environment will lead to a different redox potential. However, a detailed molecular picture of reorganization of the environment upon reduction is still unclear. To unravel the details of the solvent reorganization during electron transfer, we have performed density functional theory-based molecular dynamics (DFT-MD) and hybrid quantum mechanics/molecular mechanics (QM/MM) simulations of the reduction of lumiflavin. Previously, we have calculated the reduction free energy curves of the redox half reactions of lumiflavin in water as a function of the instantaneous gap energy (ΔE) ( J. Chem. Theory Comput. 2013 , 9 , 3889 - 3899 ). In this work, we focus on finding the changes in the solvent environment that correlate with this ΔE reaction coordinate. Comparing the QM/MM simulations, in which the solvent is modeled with an empirical force field, with the (full) DFT-MD simulations, we find that the response through electronic polarization plays a significant role in the latter case. Also a small charge transfer between flavin and solvent is observed in the full DFT treatment. As a result, we find only in the case of the QM/MM model a strong correlation between ΔE and the (pairwise computed) electrostatic potential (ESP) at the flavin due to the solvent. By analyzing the contribution of the ESP at the flavin per solvent molecule, we cannot only distinguish between the different modes of hydration by solvent molecules that coordinate at the hydrophilic and hydrophobic sides of the flavin molecule but also quantify their contribution to the reorganization free energy by measuring the ESP fluctuations per solvent molecule.

Impact of fluorination on the photophysics of the flavin chromophore: a quantum chemical perspective

10-Methylisoalloxazine (MIA) and its fluorinated derivatives (6-9F-MIA) were investigated by means of quantum chemistry, looking into the influence of fluorination on fluorescence, absorption and inter-system crossing (ISC) in vacuum and in aqueous solution.

Cyclic Changes in Active Site Polarization and Dynamics Drive the 'Ping-pong' Kinetics in NRH:Quinone Oxidoreductase 2: An Insight from QM/MM Simulations

Quinone reductases belong to the family of flavin-dependent oxidoreductases. With the redox active cofactor, flavin adenine dinucleotide, quinone reductases are known to utilize a 'ping-pong' kinetic mechanism during catalysis in which a hydride is bounced back and forth between flavin and its two substrates. However, the continuation of this catalytic cycle requires product displacement steps, where the product of one redox half-cycle is displaced by the substrate of the next half-cycle. Using improved hybrid quantum mechanical/molecular mechanical simulations, both the catalytic hydride transfer and the product displacement reactions were studied in NRH:quinone oxidoreductase 2. Initially, the self-consistent charge-density functional tight binding theory was used to describe flavin ring and the substrate atoms, while embedded in the molecular mechanically-treated solvated active site. Then, for each step of the catalytic cycle, a further improvement of energetics was made using density functional theory-based corrections. The present study showcases an integrated interplay of solvation, protonation, and protein matrix-induced polarization as the driving force behind the thermodynamic wheel of the 'ping-pong' kinetics. Reported here is the first-principles model of the 'ping-pong' kinetics that portrays how cyclic changes in the active site polarization and dynamics govern the oscillatory hydride transfer and product displacement in this enzyme.

The first water coordination sphere of lanthanide(iii) motexafins (Ln-Motex2+, Ln = La, Gd, Lu) and its effects on structures, reduction potentials and UV-vis absorption spectra. Theoretical studies

An explicit treatment of strongly bound water molecules is mandatory to calculate correct UV-vis absorption spectra of lanthanide(iii) motexafins.

Insight into the kinetics and thermodynamics of the hydride transfer reactions between quinones and lumiflavin: a density functional theory study

The kinetics and equilibrium of the hydride transfer reaction between lumiflavin and a number of substituted quinones was studied using density functional theory. The impact of electron withdrawing/donating substituents on the redox potentials of quinones was studied. In addition, the role of these substituents on the kinetics of the hydride transfer reaction with lumiflavin was investigated in detail under the transition state (TS) theory assumption. The hydride transfer reactions were found to be more favorable for an electron-withdrawing substituent. The activation barrier exhibited a quadratic relationship with the driving force of these reactions as derived under the formalism of modified Marcus theory. The present study found a significant extent of electron delocalization in the TS that is stabilized by enhanced electrostatic, polarization, and exchange interactions. Analysis of geometry, bond-orders, and energetics revealed a predominant parallel (Leffler-Hammond) effect on the TS. Closer scrutiny reveals that electron-withdrawing substituents, although located on the acceptor ring, reduce the N-H bond order of the donor fragment in the precursor complex. Carried out in the gas-phase, this is the first ever report of a theoretical study of flavin's hydride transfer reactions with quinones, providing an unfiltered view of the electronic effect on the nuclear reorganization of donor-acceptor complexes.

Single Hydrogen Bond Donation from Flavin N5 to Proximal Asparagine Ensures FAD Reduction in DNA Photolyase

The spread of the absorbance of the stable FADH(•) radical (300-700 nm) allows CPD photolyase to highly efficiently form FADH(-), making it functional for DNA repair. In this study, FTIR spectroscopy detected a strong hydrogen bond, from FAD N5-H to the carbonyl group of the Asn378 side chain, that is modulated by the redox state of FAD. The observed characteristic frequency shifts were reproduced in quantum-mechanical models of the flavin binding site, which were then employed to elucidate redox tuning governed by Asn378. We demonstrate that enhanced hydrogen bonding of the Asn378 side chain with the FADH(•) radical increases thermodynamic stabilization of the radical state, and further ensures kinetic stabilization and accumulation of the fully reduced FADH(-) state.

Quantum Chemical Studies on Stability and Chemical Activities in Calcium Ion Bound Calmodulin Loops

Quantum chemical (QC) calculations for macromolecules require truncation of the molecule, highlighting the portion of interest due to heavy computation cost. As a result, an estimation of the effects of truncation is important to interpret the energy spectrum of such calculations. We perform density functional theory based QC calculations on calcium ion bound EF-hand loops of Calmodulin isolated from the crystal structure in an implicit solvent. We find that the terminal contributions of neutral capping are negligible across the entire ground-state energy spectrum. The coordination energy range and the nature of hybridization of the coordination state molecular orbitals remain qualitatively similar across these loops. While the HOMO and LUMO of loops in the N-terminal domain are dominated by the acidic aspartates, and the polar/hydrophobic residues, respectively, these levels of the C-terminal domain loops show strong localized electron density on the phenyl rings of the tyrosines. The Fukui index calculation identifies the hydroxyl oxygen in the phenyl ring of Y99 as a potent nucleophile. Our analysis indicates a general way of interpreting the electronic energy spectra to understand stability and functions of large biomolecules where the truncation of the molecule and, hence, the terminal capping effects are inevitable.

Effect of stacking interactions on the thermodynamics and kinetics of lumiflavin: a study with improved density functionals and density functional tight-binding protocol

The π-π stacking interaction between lumiflavin and a number of π-electron-rich molecules has been studied by density functional theory using several new-generation density functionals. Six known lumiflavin-aromatic adducts were used and the models were evaluated by comparing the geometry and energetics with experimental results. The study found that dispersion-corrected and hybrid functionals with larger (>50%) Hartree-Fock exchanges produced superior results in modeling thermodynamic characteristics of these complexes. The functional producing the best energetics for these model systems was used to study the stacking interactions of lumiflavin with biologically relevant aromatic groups. Additionally, the reduction of flavin-in the presence of both a hydride donor and a nondonor π-electronic system was also studied. Weak interactions were observed in the stacked lumiflavin complexes of benzene, phenol, and indole, mimicking phenyl alanine, tryptophan, and tyrosine side chains, respectively, of an enzyme. The stacked complex of naphthalene and flavin showed little change in flavin's redox potential indicating insignificant effect on the thermodynamics of the hydride transfer reaction. In contrast, the hydride transfer reaction with the hydride donor N-methyl nicotinamide tells a different story, as the transition state was found to be strongly impacted by the stacking interactions. A comparison of performance between the density functional theory (DFT) and the computationally less expensive dispersion-corrected self-consistent density functional tight-binding (SCC-DFTB-D) theory revealed that the latter produces consistent energetics for this hydride transfer reaction and additional DFT-computed perturbative corrections could significantly improve these results.

Computational electrochemistry: prediction of liquid-phase reduction potentials

This article reviews recent developments and applications in the area of computational electrochemistry. Our focus is on predicting the reduction potentials of electron transfer and other electrochemical reactions and half-reactions in both aqueous and nonaqueous solutions. Topics covered include various computational protocols that combine quantum mechanical electronic structure methods (such as density functional theory) with implicit-solvent models, explicit-solvent protocols that employ Monte Carlo or molecular dynamics simulations (for example, Car-Parrinello molecular dynamics using the grand canonical ensemble formalism), and the Marcus theory of electronic charge transfer. We also review computational approaches based on empirical relationships between molecular and electronic structure and electron transfer reactivity. The scope of the implicit-solvent protocols is emphasized, and the present status of the theory and future directions are outlined.

Probing protonation sites of isolated flavins using IR spectroscopy: from lumichrome to the cofactor flavin mononucleotide

Infrared spectra of the isolated protonated flavin molecules lumichrome, lumiflavin, riboflavin (vitamin B2), and the biologically important cofactor flavin mononucleotide are measured in the fingerprint region (600-1850 cm(-1)) by means of IR multiple-photon dissociation (IRMPD) spectroscopy. Using density functional theory calculations, the geometries, relative energies, and linear IR absorption spectra of several low-energy isomers are calculated. Comparison of the calculated IR spectra with the measured IRMPD spectra reveals that the N10 substituent on the isoalloxazine ring influences the protonation site of the flavin. Lumichrome, with a hydrogen substituent, is only stable as the N1-protonated tautomer and protonates at N5 of the pyrazine ring. The presence of the ribityl unit in riboflavin leads to protonation at N1 of the pyrimidinedione moiety, and methyl substitution in lumiflavin stabilizes the tautomer that is protonated at O2. In contrast, flavin mononucleotide exists as both the O2- and N1-protonated tautomers. The frequencies and relative intensities of the two C=O stretch vibrations in protonated flavins serve as reliable indicators for their protonation site.

Acidity constants of lumiflavin from first principles molecular dynamics simulations

DFT-based molecular dynamics simulations predict the acidity of lumiflavin in different redox states.

Density functional theory calculations of the redox potentials of actinide(VI)/actinide(V) couple in water

The measured redox potential of an actinide at an electrode surface involves the transfer of a single electron from the electrode surface on to the actinide center. Before electron transfer takes place, the complexing ligands and molecules of solvation need to become structurally arranged such that the electron transfer is at its most favorable. Following the electron transfer, there is further rearrangement to obtain the minimum energy structure for the reduced state. As such, there are three parts to the total energy cycle required to take the complex from its ground state oxidized form to its ground state reduced form. The first part of the energy comes from the structural rearrangement and solvation energies of the actinide species before the electron transfer or charge transfer process; the second part, the energy of the electron transfer; the third part, the energy required to reorganize the ligands and molecules of solvation around the reduced species. The time resolution of electrochemical techniques such as cyclic voltammetry is inadequate to determine to what extent bond and solvation rearrangement occurs before or after electron transfer; only for a couple to be classed as reversible is it fast in terms of the experimental time. Consequently, the partitioning of the energy theoretically is of importance to obtain good experimental agreement. Here we investigate the magnitude of the instantaneous charge transfer through calculating the fast one electron reduction energies of AnO2(H2O)n(2+), where An = U, Np, and Pu, for n = 4-6, in solution without inclusion of the structural optimization energy of the reduced form. These calculations have been performed using a number of DFT functionals, including the recently developed functionals of Zhao and Truhlar. The results obtained for calculated electron affinities in the aqueous phase for the AnO2(H2O)5(2+/+) couples are within 0.04 V of accepted experimental redox potentials, nearly an order of magnitude improvement on previous calculated standard potentials E(0) values, obtained using both DFT and high level multireference approaches.

Manipulation of microbial extracellular electron transfer by changing molecular structure of phenazine-type redox mediators

Phenazines, as a type of electron shuttle, are involved in various biological processes to facilitate microbial energy metabolism and electron transfer. They constitute a large group of nitrogen-containing heterocyclic compounds, which can be produced by a diverse range of bacteria or by artificial synthesis. They vary significantly in their properties, depending mainly on the nature and position of substitutent group. Thus, it is of great interest to find out the most favorable substituent type and molecular structure of phenazines for electron transfer routes. Here, the impacts of the substituent group on the reduction potentials of phenazine-type redox mediators in aqueous solution were investigated by quantum chemical calculations, and the calculation results were further validated with experimental data. The results show that the reaction free energy was substantially affected by the location of substituent groups on the phenazine molecule and the protonated water clusters. For the main proton addition process, the phenazines substituted with electron-donating groups and those with electron-withdrawing groups interacted with different protonated water clusters, attributed to the proximity effect of water molecules on proton transfer. Thus, high energy conversion efficiency could be achieved by controlling electron flow route with appropriate substituted phenazines to reduce the biological energy acquisition. This study provides useful information for designing efficient redox mediators to promote electron transfer between microbes and terminal acceptors, which are essential to bioenergy recovery from wastes and environmental bioremediation.

Computational Study of the pKa Values of Potential Catalytic Residues in the Active Site of Monoamine Oxidase B

Monoamine oxidase (MAO), which exists in two isozymic forms, MAO A and MAO B, is an important flavoenzyme responsible for the metabolism of amine neurotransmitters such as dopamine, serotonin, and norepinephrine. Despite extensive research effort, neither the catalytic nor the inhibition mechanisms of MAO have been completely understood. There has also been dispute with regard to the protonation state of the substrate upon entering the active site, as well as the identity of residues that are important for the initial deprotonation of irreversible acetylenic inhibitors, in accordance with the recently proposed mechanism. Therefore, in order to investigate features essential for the modes of action of MAO, we have calculated pKa values of three relevant tyrosine residues in the MAO B active site, with and without dopamine bound as the substrate (as well as the pKa of the dopamine itself in the active site). The calculated pKa values for Tyr188, Tyr398, and Tyr435 in the complex are found to be shifted upward to 13.0, 13.7, and 14.7, respectively, relative to 10.1 in aqueous solution, ruling out the likelihood that they are viable proton acceptors. The altered tyrosine pKa values could be rationalized as an interplay of two opposing effects: insertion of positively charged bulky dopamine that lowers tyrosine pKa values, and subsequent removal of water molecules from the active site that elevates tyrosine pKa values, in which the latter prevails. Additionally, the pKa value of the bound dopamine (8.8) is practically unchanged compared to the corresponding value in aqueous solution (8.9), as would be expected from a charged amine placed in a hydrophobic active site consisting of aromatic moieties. We also observed potentially favorable cation-π interactions between the -NH3(+) group on dopamine and aromatic moieties, which provide a stabilizing effect to the charged fragment. Thus, we offer here theoretical evidence that the amine is most likely to be present in the active site in its protonated form, which is similar to the conclusion from experimental studies of MAO A (Jones et al. J. Neural Trans.2007, 114, 707-712). However, the free energy cost of transferring the proton from the substrate to the bulk solvent is only 1.9 kcal mol(-1), leaving open the possibility that the amine enters the chemical step in its neutral form. In conjunction with additional experimental and computational work, the data presented here should lead toward a deeper understanding of mechanisms of the catalytic activity and irreversible inhibition of MAO B, which can allow for the design of novel and improved MAO B inhibitors.

Interplay of flavin's redox states and protein dynamics: an insight from QM/MM simulations of dihydronicotinamide riboside quinone oxidoreductase 2

Dihydronicotinamide riboside quinone oxidoreductase 2 is known to catalyze a two-electron reduction of quinone to hydroquinone using its cofactor, flavin adenine dinucleotide. Using quantum mechanical/molecular mechanical simulations, we have computed the reorganization free energies of the electron and proton transfer processes of flavin in the free state as well as when it is bound in the active site of the enzyme. The calculated energetics for electron transfer processes demonstrate that the enzyme active site lowers the reorganization energy for the redox process as compared to the enzyme-free aqueous state. This is most apparent in the two electron reduction step, which eliminates the possibility of flavosemiquinone generation. In addition, essential dynamics study of the simulated motions revealed spectacular changes in the principal components of atomic fluctuations upon reduction of flavin. This alteration of active site dynamics provides an insight into the "ping-pong" kinetics exhibited by the enzyme upon a change in the redox state of the enzyme-bound flavin. A charge perturbation analysis provides further support that the observed change in dynamics is correlated with the change in energetics due to the altered electrostatic interactions between the flavin ring and the active site residues. This study shows that the effect of electrostatic preorganization goes beyond the chemical catalysis as it strongly impacts the postcatalytic intrinsic protein dynamics.