Effect Of Protein Environment Research Articles

Computational Exploration of Minimum Energy Reaction Pathway of N2O Formation from Intermediate I of P450nor Using an Active Center Model.

P450nor is a heme-containing enzyme that catalyzes the conversion of nitric oxide (NO) to nitrous oxide (N2O). Its catalytic mechanism has attracted attention in chemistry, biology, and environmental engineering. The catalytic cycle of P450nor is proposed to consist of three major steps. The reaction mechanism for the last step, N2O generation, remains unknown. In this study, the reaction pathway of the N2O generation from the intermediate I was explored with the B3LYP calculations using an active center model after the examination of the validity of the model. In the validation, we compared the heme distortions between P450nor and other oxidoreductases, suggesting a small effect of protein environment on the N2O generation reaction in P450nor. We then evaluated the electrostatic environment effect of P450nor on the hydride affinity to the active site with quantum mechanics/molecular mechanics (QM/MM) calculations, confirming that the affinity was unchanged with or without the protein environment. The active center model for P450nor showed that the N2O generation process in the enzymatic reaction undergoes a reasonable barrier height without protein environment. Consequently, our findings strongly suggest that the N2O generation reaction from the intermediate I depends sorely on the intrinsic reactivity of the heme cofactor bound on cysteine residue.

Characterization of Light-Induced, Short-Lived Interacting Radicals in the Active Site of Flavoprotein Ferredoxin-NADP+ Oxidoreductase.

Radicals of flavin adenine dinucleotide (FAD), as well as tyrosine and tryptophan, are widely involved as key reactive intermediates during electron-transfer (ET) reactions in flavoproteins. Due to the high reactivity of these species and their corresponding short lifetime, characterization of these intermediates in functional processes of flavoproteins is usually challenging but can be achieved by ultrafast spectroscopic studies of light-activatable flavoproteins. In ferredoxin-NADP+ oxidoreductase from Bacillus subtilis (BsFNR), fluorescence of the FAD cofactor that very closely interacts with a neighboring tyrosine residue (Tyr50) is strongly quenched. Here we study short-lived photoproducts of this enzyme and its variants, with Tyr50 replaced by tryptophan or glycine. Using time-resolved fluorescence and absorption spectroscopies, we show that, upon the excitation of WT BsFNR, ultrafast ET from Tyr50 to the excited FAD cofactor occurs in ∼260 fs, an order of magnitude faster than the decay by charge recombination, facilitating the characterization of the reaction intermediates in the charge-separated state with respect to other recently studied systems. These studies are corroborated by experiments on the Y50W mutant protein, which yield photoproducts qualitatively similar to those observed in other tryptophan-bearing flavoproteins. By combining the experimental results with molecular dynamics simulations and quantum mechanics calculations, we investigate in detail the effects of protein environment and relaxations on the spectral properties of those radical intermediates and demonstrate that the spectral features of radical anionic FAD are highly sensitive to its environment, and in particular to the dynamics and nature of the counterions formed in the photoproducts. Altogether, comprehensive characterizations are provided for important radical intermediates that are generally involved in functional processes of flavoproteins.

Catalytic center of cytochrome c oxidase: Effects of protein environment on pKa values of cub histidine ligands

The molecular mechanism by which electron transfer (ET) is coupled to proton pumping in cytochrome oxidase is one of the main unsolved problems in biochemistry. Particularly, the nature and position of the proton-loading site is under dispute. The CuB complex has three ligated histidines, whereas only His290 and His291 are ionizable sites with the same pKa values in aqueous solution, but apparently quite different ones within the enzyme. Earlier, a model of proton pumping with the central role of His290 was proposed. Recent calculations indicate that the His291 ligand of the CuB center might play the role of the pumping element, since its protonation state depends on the oxidation state of the binuclear complex (BNC). The present electrostatic study was applied to assess the role of the protein environment on the acidity of the two histidines. Their pKa values and effects of different energy terms were evaluated to discover the nature of their diverse behavior in the enzyme. Here, a new set of pKa values for the non-standard model compounds within the BNC was applied. The enhanced results are compared with results of previous studies in the light of the plausible proton pumping mechanism. The obtained microscopic and apparent pKa values in the oxidized state of BNC are virtually the same, indicating that deprotonated form of His291 accounts for the large pKa increase of His290, since then both titratable sites on then CuB center cannot simultaneously be in the charged state. The present results support the underlined His291 pumping model.

Understanding the Effect of Single Cysteine Mutations on Gold Nanoclusters as Studied by Spectroscopy and Density Functional Theory Modeling.

Fluorescent metal nanoclusters have generated considerable excitement in nanobiotechnology, particularly in the applications of biolabeling, targeted delivery, and biological sensing. The present work is an experimental and computational study that aims to understand the effects of protein environment on the synthesis and electronic properties of gold nanoclusters. MPT63, a drug target of Mycobacterium tuberculosis, was used as the template protein to synthesize, for the first time, gold nanoclusters at a low micromolar concentration of the protein. Two single cysteine mutants of MPT63, namely, MPT63Gly20Cys (mutant I) and MPT63Gly40Cys (mutant II) were employed for this study. The experimental results show that cysteine residues positioned in two different regions of the protein induce varying electronic states of the nanoclusters depending on the surrounding amino acids. A mixture of five-atom and eight-atom clusters was generated for each mutant, and the former was found to be predominant in both cases. Computational studies, including density functional theory (DFT), frontier molecular orbital (FMO), and natural bond orbital (NBO) calculations, validated the experimental observations. The as-prepared protein-stabilized nanoclusters were found to have applications in the imaging of live cells.

Protein environmental effects on iron-sulfur clusters: A set of rules for constructing computational models for inner and outer coordination spheres.

The structural properties and reactivity of iron-sulfur proteins are greatly affected by interactions between the prosthetic groups and the surrounding amino acid residues. Thus, quantum chemical investigations of the structure and properties of protein-bound iron-sulfur clusters can be severely limited by truncation of computational models. The aim of this study was to identify, a priori, significant interactions that must be included in a quantum chemical model. Using the [2Fe-2S] accessory cluster of the FeFe-hydrogenase as a demonstrative example with rich electronic structural features, the electrostatic and covalent effects of the surrounding side chains, charged groups, and backbone moieties were systematically mapped through density functional theoretical calculations. Electron affinities, spin density differences, and delocalization indexes from the quantum theory of atoms in molecules were used to evaluate the importance of each interaction. Case studies for hydrogen bonding and charged side-chain interactions were used to develop selection rules regarding the significance of a given protein environmental effect. A set of general rules is proposed for constructing quantum chemical models for iron-sulfur active sites that capture all significant interactions from the protein environment. This methodology was applied to our previously used models in galactose oxidase and the 6Fe-cluster of FeFe-hydrogenase. © 2016 Wiley Periodicals, Inc.

Linkage Isomerization in Nitrosothiols (RSNOs): The X-ray Crystal Structure of an S-nitrosocysteine and DFT Analysis of its Metastable MS1 and MS2 Isomers

Density functional theory calculations reveal that the previously unreported nitrosothiol MS1 and MS2 linkage isomers are metastable species with energies similar to those determined for other experimentally observed XNO linkage isomers.Nitrosation of protein cysteines to yield S-nitrosocysteine (cySNO) is an important component of nitric oxide biology. However, high-resolution structural data of cySNO in the absence of protein environment effects is lacking, as is information on cySNO linkage isomers. We report the ordered X-ray crystal structure of a cySNO compound. We employ density functional theory calculations to probe the ground-state and linkage isomers of cySNO, CH3SNO, and CF3SNO. The HOMO of CH3SNO contains an intramolecular CH···O interaction that helps rationalize the stability of the cis conformation; this interaction is absent in CF3SNO, which favors the trans conformation. We show that the isonitroso MS1 (RSON) and side-on MS2 (RS(η2-ON)) linkage isomers are metastable species, with MS1 and MS2 energies similar to those of other experimentally observed nitroso linkage isomers.

Thiol versus hydroxamate as zinc binding group in HDAC inhibition: An ab initio QM/MM molecular dynamics study.

Zinc-dependent histone deacetylases (HDACs) play a critical role in transcriptional repression and gene silencing, and are among the most attractive targets for the development of new therapeutics against cancer and various other diseases. Two HDAC inhibitors have been approved by FDA as anti-cancer drugs: one is SAHA whose hydroxamate is directly bound to zinc, the other is FK228 whose active form may use thiol as the zinc binding group. In spite of extensive studies, it remains to be ambiguous regarding how thiol and hydroxamate are bound to the zinc active site of HDACs. In this work, our computational approaches center on Born-Oppenheimer ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics with umbrella sampling, which allow for modeling of the zinc active site with reasonable accuracy while properly including dynamics and effects of protein environment. Meanwhile, an improved short-long effective function (SLEF2) to describe non-bonded interactions between zinc and other atoms has been employed in initial MM equilibrations. Our ab initio QM/MM MD simulations have confirmed that hydroxamate is neutral when it is bound to HDAC8, and found that thiol is deprotonated when directly bound to zinc in the HDAC active site. By comparing thiol and hydroxamate, our results elucidated the differences in their binding environment in the HDAC active sites, and emphasized the importance of the linker design to achieve more specific binding toward class IIa HDACs.

The Allosteric Regulation of Axial/Rhombic Population in a “Type 1” Copper Site: Multi-Edge X-ray Absorption Spectroscopic and Density Functional Studies of Pseudoazurin

Abstract The co-existence of “axial” and “rhombic” coordination environments has been demonstrated in a “Type 1” copper site of Pseudoazurin. This observation opens up previously not considered interpretations for the relationship between geometry and electronic structure of the four coordinate copper site. The Met16 variants of pseudoazurin were considered as model systems for investigating the effect of weak interactions from the second coordination sphere. The correlation between geometric and electronic structures of “Type 1” copper site was evaluated by the multi-edge (Cu K-edge and S K-edge) X-ray absorption spectroscopy (XAS) of Met16 variants of pseudoazurin. The co-existing axial and rhombic sites in pseudoazurin were characterized by Cu–ligand distances, effective nuclear charge, and Cu–S(Cys) covalency from XAS. The XAS results were correlated with DFT calculations for investigating the effect of protein environment from the inner-sphere and beyond around the Cu site. The combined experimental and theoretical results support the presence of a close correlation between outer sphere environment and inner sphere coordination environment. This is achieved in pseudoazurin by a previously undisclosed allosteric effect that involves a rearrangement of the protein tertiary structure.

Electronic circular dichroism of fluorescent proteins: a computational study.

The electronic circular dichroism (ECD) properties of the green fluorescent protein and other fluorescent proteins have been calculated with density functional theory. The influence of different embedding models on the ECD signal of the chromophore has been investigated by modeling the protein environment by the polarizable continuum model (QM/PCM), by the polarizable embedding model (PE-QM/MM), by treating the minimal environment quantum mechanically at the same footing as the chromophore (QM/QM), and by adding the remaining part of the protein by means of PCM (QM/QM/PCM). The rotatory strength is found to be more sensitive than the oscillatory strength to changes in the geometry of the chromophore and its surroundings and to the type of embedding model used. In general, explicit embedding of the surrounding protein (PE-QM/MM or QM/QM) induces an increase in the rotatory strength of the chromophore. Explicit inclusion of the whole protein through polarizable embedding is found to be an affordable embedding model that gives the correct sign of the rotatory strength for all fluorescent proteins. PCM is useful as a first approximation to protein environment effects, but as a rule seems to underestimate the rotatory strength.

Molecular-Level Insight into the Spectral Tuning Mechanism of the DsRed Chromophore.

We present a detailed study of the protein environmental effects on the one- and two-photon absorption (1PA and 2PA, respectively) properties of the S0-S1 transition in the DsRed protein using the polarizable embedding density functional theory formalism. We find that steric factors and chromophore-protein interactions act in concert to enhance the 2PA activity inside the protein while adversely blue-shifting the 1PA maximum. A two-state model reveals that the 2PA intensity gain is primarily governed by the increased change in the permanent dipole moment between the ground and the excited states acquired inside the protein. Our results indicate that this mainly is attributable to counter-directional contributions stemming from Lys163 and the conserved Arg95 with the former additionally identified as a key residue in the color tuning mechanism. The results provide new insight into the tuning mechanism of DsRed and suggest a possible strategy for simultaneous improvement of its 1PA and 2PA properties.

Theoretical studies of electrostatic effect of protein environment on electronic structures and chemical indices of the active site of oxygenated and deoxygenated hemerythrin

AbstractWe have studied the electrostatic effect of protein environment on electronic structures and chemical indices of the active site of oxygenated and deoxygenated hemerythrin (oxyHr and deoxyHr) using hybrid density functional calculations. We found that the UBHandHLYP method can be regarded as a reliable method for the nature of chemical bonds of Hr even in the protein. The electrostatic effects of the protein environment on oxyHr and deoxyHr greatly improve the magnetic couplings and provide the polarized electronic structures. The electrostatic effects on the active site of oxyHr hardly affect the shape and symmetry of singly occupied natural orbitals but contribute to the enhanced orbital interaction between the active core of Hr and oxygen molecule. This strong orbital interaction leads to easy binding of oxygen molecule to Hr for oxygen transport. Thus, our computation suggests that protein environment regulates the electronic structures of the active site to intensify the intrinsic ability of the active site and to express biological functions such as oxygen transport, energy conversion, and electron transfer. © 2012 Wiley Periodicals, Inc.

The Effect of Protein Environment on Photoexcitation Properties of Retinal

Retinal is the photon absorbing chromophore of rhodopsin and other visual pigments, enabling the vertebrate vision process. The effects of the protein environment on the primary photoexcitation process of retinal were studied by time-dependent density functional theory (TDDFT) and the algebraic diagrammatic construction through second order (ADC(2)) combined with our recently introduced reduction of virtual space (RVS) approximation method. The calculations were performed on large full quantum chemical cluster models of the bluecone (BC) and rhodopsin (Rh) pigments with 165-171 atoms. Absorption wavelengths of 441 and 491 nm were obtained at the B3LYP level of theory for the respective models, which agree well with the experimental values of 414 and 498 nm. Electrostatic rather than structural strain effects were shown to dominate the spectral tuning properties of the surrounding protein. The Schiff base retinal and a neighboring Glu-113 residue were found to have comparable proton affinities in the ground state of the BC model, whereas in the excited state, the proton affinity of the Schiff base is 5.9 kcal/mol (0.26 eV) higher. For the ground and excited states of the Rh model, the proton affinity of the Schiff base is 3.2 kcal/mol (0.14 eV) and 7.9 kcal/mol (0.34 eV) higher than for Glu-113, respectively. The protein environment was found to enhance the bond length alternation (BLA) of the retinyl chain and blueshift the first absorption maxima of the protonated Schiff base in the BC and Rh models relative to the chromophore in the gas phase. The protein environment was also found to decrease the intensity of the second excited state, thus improving the quantum yield of the photoexcitation process. Relaxation of the BC model on the excited state potential energy surface led to a vanishing BLA around the isomerization center of the conjugated retinyl chain, rendering the retinal accessible for cis-trans isomerization. The energy of the relaxed excited state was found to be 30 kcal/mol (1.3 eV) above the minimum ground state energy, and might be related to the transition state of the thermal activation process.

Theoretical Studies of Chromophore Maturation in the Wild-Type Green Fluorescent Protein: ONIOM(DFT:MM) Investigation of the Mechanism of Cyclization

The availability of a gene encoding green fluorescence immediately stimulates interest in the puzzle of autocatalytic formation of the green fluorescent protein (GFP) chromophore. Numerous experimental and theoretical studies have indicated that cyclization is the first and most important step in the maturation process of the GFP. In our previous paper based on cluster models [J. Phys. Chem. B2010, 114, 9698-9705], two possible mechanisms have been investigated with the conclusion that the backbone condensation initiated by deprotonation of the Gly67 amide nitrogen is easier than deprotonation of the Tyr66 α-carbon. However, the impact of the protein environment on the reaction mechanism remains to be explored. In this paper, we investigated the two possible mechanisms with inclusion of protein environmental effects by using molecular dynamics (MD) and combined quantum mechanics/molecular mechanics (QM/MM) calculations. Our calculations reveal no hydrogen bonding network that would facilitate deprotonation of the amide nitrogen of Gly67, although it is the lower energy pathway in the cluster model system. Contrastingly, there is a hydrogen bonding network between Tyr66 α-carbon and Glu222, which is in good agreement with X-ray data. The ONIOM studies show that proton transfer from Tyr66 α-carbon to Glu222 is a long-distance charge transfer process. The charge distribution of the MM region has a significant perturbation to the wave function for the QM region, with the QM energy for the proton transfer product being increased under the influence of the electrostatic protein environment. The barrier for the rate-limiting step in cyclization is quite high, about 40.0 kcal/mol in the case of ONIOM-EE.

The Effects of Protein Environment and Dispersion on the Formation of Ferric-Superoxide Species in myo-Inositol Oxygenase (MIOX): A Combined ONIOM(DFT:MM) and Energy Decomposition Analysis

The catalytic reaction of myo-inositol oxygenase, a nonheme diiron enzyme, is initiated by the binding of an O(2) molecule to the ferrous center of a mixed-valence Fe(II)Fe(III) intermediate. This generates a (superoxo)Fe(III)Fe(III) reactive species that abstracts a hydrogen atom from the myo-inositol substrate. To understand the effects of protein environment and intracluster dispersion on this O(2)-binding process, we undertook a combined ONIOM(B3LYP:AMBER) and energy decomposition analysis. The interaction energy between the active site and the thousands of atoms present in the protein environment was decomposed into electrostatic, van der Waals (vdW) and polarization terms. These terms were further decomposed into contributions from individual amino acid residues. The dispersion effect, which is not adequately accounted for by the B3LYP method, was estimated in an empirical manner. The results show that the electrostatic, vdW, and polarization effects slightly enhance the O(2) binding process. The dispersion effect enhances O(2) binding more significantly than these effects. Despite these stabilizing effects, the entropy effect disfavors O(2) binding, making the process almost thermoneutral.

Optical properties of flavin mononucleotide: A QM/MM study of protein environment effects

We investigate the optical properties of flavin mononucleotide (FMN) as chromophore of a bacterial protein. FMN is studied both in the oxidized and anionic state, by combining hybrid quantum mechanics/molecular mechanics (QM/MM) dynamics with time-dependent density functional theory (TDDFT). Our first-principles calculations show that the inclusion of the protein biological environment is crucial to reproduce the experimental optical spectra.

Effect of Protein Environment on Electronically Excited and Ionized States of the Green Fluorescent Protein Chromophore

The effect of the protein environment on the electronic structure of the green fluorescent protein (GFP) chromophore is investigated by QM/MM (quantum mechanics/molecular mechanics) calculations. The protein has very small effect on the excitation energy of the bright absorbing and the lowest triplet states of the anionic GFP chromophore, deprotonated 4-hydroxybenzylidene-2,3-dimethylimidazolinone (HBDI) anion, however, it increases vertical detachment energy from 2.5 eV (gas-phase deprotonated HBDI anion) to 5.0 eV (solvated protein). We also investigated possible existence of the charge-transfer-to-solvent (CTTS) states associated with the GFP chromophore. Although precursors of such states appear in cluster calculations, a tightly packed structure of the protein prevents the formation of the CTTS states in this system. Motivated by a recently discovered new type of photoconversion, oxidative redding, we characterized the redox properties of GFP. The computed standard reduction potential of the anionic form of GFP is 0.47 V (for the GFP(•) + 1e → GFP(-) reaction), and the reduction potential at physiological conditions (pH 7, T = 25 °C) is 0.06 V.

Multireference and Multiconfiguration Ab Initio Methods in Heme-Related Systems: What Have We Learned So Far?

This work reviews the recent applications of ab initio multireference/multiconfiguration (MR/MC) electronic structure methods to heme-related systems, involving tetra-, penta-, and hexa-coordinate species, as well as the high-valent iron-oxo species. The current accuracy of these methods in the various systems is discussed, with special attention to potential sources of systematic errors. Thus, the review summarizes and tries to rationalize the key elements of MR/MC calculations, namely, the choice of the employed active space, especially the so-called double-shell effect that has already been recognized to be important in transition-metal-containing systems, and the impact of these elements on the spin-state energetics of heme species, as well as on the bonding mechanism of small molecules to the heme. It is shown that expansion of the MC wave function into one based on localized orbitals provides a compact and insightful view on some otherwise complex electronic structures. The effects of protein environment on the MR/MC results are summarized for the few available quantum mechanical/molecular mechanical (QM/MM) studies. Comparisons with corresponding DFT results are also made wherever available. Potential future directions are proposed.

Spectroscopic Comparison of Photogenerated Tryptophan Radicals in Azurin: Effects of Local Environment and Structure

Tryptophan radicals play a significant role in mediating biological electron transfer. We report the photogeneration of a long-lived, neutral tryptophan radical (Az48W*) from the native residue tryptophan-48 in the hydrophobic core of azurin. The optical absorption, electron paramagnetic resonance, and resonance Raman spectra strongly support the formation of a neutral radical, and the data are consistent with direct electron transfer between tryptophan and the copper(II) center. Spectra of the long-lived Az48W* species are compared to those of a previously studied, solvent-exposed radical at position 108 to identify signatures of tryptophan radicals that are sensitive to the local environment. The absorption maxima of Az48W* display an approximately 23 nm hypsochromic shift in the nonpolar environment. The majority of the resonance Raman frequencies are downshifted by approximately 7 cm(-1) relative to the solvent-exposed radical, and large changes in intensity are observed for some modes. The resonance Raman excitation profiles for Az48W* exhibit distinct maxima within the absorption envelope. Electron paramagnetic resonance spectroscopy yields spectra with partially resolved lines caused by hyperfine couplings; the differences between the coupling constants for the buried and solvent-exposed radical are primarily caused by variations in structure. The insights gained by electronic, vibrational, and magnetic resonance spectroscopy enhance our fundamental understanding of the effects of protein environment on radical properties. Hypotheses for the proton transfer pathway within azurin and a deprotonation rate of approximately 5 x 10(6) s(-1) are proposed.

The Chromophore Structure of the Cyanobacterial Phytochrome Cph1 As Predicted by Time-Dependent Density Functional Theory

The UV-vis absorption spectra of the photoreceptor chromophores biliverdin (BV) in the ZZZssa conformation and the phycocyanobilin (PCB) with conformations ZZZssa and ZZZasa have been investigated by means of time-dependent density functional theory (TD-DFT) with a polarized continuum model. The three systems are studied in different conditions to include protonation, solvation- and protein-environmental effects on gas phase and available X-ray structures. The crystal structures of BV in bacteriophytochrome of Deinococcus radiodurans and PCB in C-Phycocyanin serve to calibrate the performance of the TD-DFT method and allow estimating the spectral shifts created when gas phase structures instead of a proper environment are used. In contrast, the structure of PCB in the cyanobacterial phytochrome Cphl is unknown. The excellent agreement of the theoretical spectrum with experimentally recorded data for the PCB in the cyanobacterial phytochrome Cph1 strongly supports a semicyclic ZZZssa structure, similar to that found for the BV chromophore.

2P354 Theoretical analysis of the electronic asymmetry of the special pair in the photosynthetic reaction center : Effect of protein environment(43. Photosynthesis,Poster Session,Abstract,Meeting Program of EABS & BSJ 2006)

2P354 Theoretical analysis of the electronic asymmetry of the special pair in the photosynthetic reaction center : Effect of protein environment(43. Photosynthesis,Poster Session,Abstract,Meeting Program of EABS & BSJ 2006)