Computational studies on imidazole heme conformations

Study on enhancement of hemoglobin antitoxic ability modified with chromium and ruthenium

Hemoglobin is essential for carrying oxygen (O₂) in the blood. However, its ability to bind excessively to carbon monoxide (CO) makes it susceptible to CO poisoning. To reduce the risk of CO poisoning, Cr-based heme and Ru-based heme were selected from among many transition metal-based hemes based on their characteristics of adsorption conformation, binding intensity, spin multiplicity, and electronic properties. The results showed that hemoglobin modified by Cr-based heme and Ru-based heme had strong anti-CO poisoning abilities. The Cr-based heme and Ru-based heme exhibited much stronger affinity for O₂ (-190.67 kJ/mol and -143.18 kJ/mol, respectively) than Fe-based heme (-44.60 kJ/mol). Moreover, Cr-based heme and Ru-based heme exhibited much weaker affinity for CO (-121.50 kJ/mol and -120.88 kJ/mol, respectively) than their affinity for O₂, suggesting that they were less likely to cause CO poisoning. The electronic structure analysis also supported this conclusion. Additionally, molecular dynamics analysis showed that hemoglobin modified by Cr-based heme and Ru-based heme was stable. Our findings offer a novel and effective strategy for enhancing the reconstructed hemoglobin's ability to bind O₂ and reduce its potential for CO poisoning.

Common Reactivity and Properties of Heme Peroxidases: A DFT Study of Their Origin

Electronic structure calculations using the density-functional theory (DFT) have been performed to analyse the effect of water molecules and protonation on the heme group of peroxidases in different redox (ferric, ferrous, compounds I and II) and spin states. Shared geometries, spectroscopic properties at the Soret region, and the thermodynamics of peroxidases are discussed. B3LYP and M06-2X density functionals with different basis sets were employed on a common molecular model of the active site (Fe-centred porphine and proximal imidazole). Computed Gibbs free energies indicate that the corresponding aquo complexes are not thermodynamically stable, supporting the five-coordinate Fe(III) centre in native ferric peroxidases, with a water molecule located at a non-bonding distance. Protonation of the ferryl oxygen of compound II is discussed in terms of thermodynamics, Fe-O bond distances, and redox properties. It is demonstrated that this protonation is necessary to account for the experimental data, and computed Gibbs free energies reveal pK_a values of compound II about 8.5-9.0. Computation indicates that the general oxidative properties of peroxidase intermediates, as well as their reactivity towards water and protons and Soret bands, are mainly controlled by the iron porphyrin and its proximal histidine ligand.

Pyrimidinetrione-imidazoles as a Unique Structural Type of Potential Agents towards Candida Albicans: Design, Synthesis and Biological Evaluation

Substantial morbidity and mortality of fungal infections have aroused concerns all over the world, and common Candida spp. currently bring about severe systemic infections. A series of pyrimidinetrione-imidazole conjugates as potentially antifungal agents were developed. Bioassays manifested that 4-fluobenzyl pyrimidinetrione imidazole 5 f exerted favorable inhibition towards C. albicans (MIC=0.002 mM), being 6.5 folds more active than clinical antifungal drug fluconazole (MIC=0.013 mM). Preliminary mechanism research indicated that compound 5 f could not only depolarize membrane potential but also permeabilize the membrane of C. albicans. Molecular docking was operated to simulate the interaction mode between molecule 5 f and CYP51. In addition, hybrid 5 f might form 5 f-DNA supramolecular complex via intercalating into DNA. The interference of membrane and DNA might contribute to its fungicidal capacity with no obvious tendency to induce the resistance against C. albicans. Conjugate 5 f endowed good blood compatibility as well as low cytotoxicity towards HeLa and HEK-293T cells.

Heme isomers substantially affect heme's electronic structure and function

Inspection of heme protein structures in the protein data bank reveals four isomers of heme characterized by different relative orientations of the vinyl side chains; remarkably, all these have been reported in multiple protein structures. Density functional theory computations explain this as due to similar energy of the isomers but with a sizable (25 kJ mol^-1) barrier to interconversion arising from restricted rotation around the conjugated bonds. The four isomers, EE, EZ, ZE, and ZZ, were then investigated as 4-coordinate hemes, as 5-coordinate deoxyhemes, in 6-coordinate O₂-adducts of globins and as compound I intermediates typical of heme peroxidases. Substantial differences were observed in electronic properties relevant to heme function: notably, the spin state energy gap of O₂-heme adducts, important for fast reversible binding of O₂, depends on the isomer state, and O₂-binding enthalpies change by up to 16 kJ mol^-1; redox potentials change by up to 0.2 V depending on the isomer, and the doublet-quartet energy splitting of compound I, central to "two-state" reactivity, is affected by up to ∼15 kJ mol^-1. These effects are consistently seen with three distinct density functionals, i.e. the effects are not method-dependent. Thus, the nature of the isomer state is an important but overlooked feature of heme chemistry and function, and previous and future studies of hemes may be reconsidered in this new context.

Reprint of PSII manganese cluster: protonation of W2, O5, O4 and His337 in the S1 state explored by combined quantum chemical and electrostatic energy computations

Photosystem II (PSII) is a membrane-bound protein complex that oxidizes water to produce energized protons, which are used to built up a proton gradient across the thylakoidal membrane in the leafs of plants. This light-driven reaction is catalyzed by withdrawing electrons from the Mn₄CaO₅-cluster (Mn-cluster) in four discrete oxidation steps [S₁-(S₄/S₀)] characterized in the Kok-cycle. In order to understand in detail the proton release events and the subsequent translocation of such energized protons, the protonation pattern of the Mn-cluster need to be elucidated. The new high-resolution PSII crystal structure from Umena, Kawakami, Shen, and Kamiya is an excellent basis to make progress in solving this problem. Following our previous work on oxidation and protonation states of the Mn-cluster, in this work, quantum chemical/electrostatic calculations were performed in order to estimate the pKa of different protons of relevant groups and atoms of the Mn-cluster such as W2, O4, O5 and His337. In broad agreement with previous experimental and theoretical work, our data suggest that W2 and His337 are likely to be in hydroxyl and neutral form, respectively, O5 and O4 to be unprotonated. This article is part of a special issue entitled: photosynthesis research for sustainability: keys to produce clean energy.

PSII manganese cluster: protonation of W2, O5, O4 and His337 in the S1 state explored by combined quantum chemical and electrostatic energy computations

Photosystem II (PSII) is a membrane-bound protein complex that oxidizes water to produce energized protons, which are used to built up a proton gradient across the thylakoidal membrane in the leafs of plants. This light-driven reaction is catalyzed by withdrawing electrons from the Mn₄CaO₅-cluster (Mn-cluster) in four discrete oxidation steps [S₁-(S₄/S₀)] characterized in the Kok-cycle. In order to understand in detail the proton release events and the subsequent translocation of such energized protons, the protonation pattern of the Mn-cluster need to be elucidated. The new high-resolution PSII crystal structure from Umena, Kawakami, Shen, and Kamiya is an excellent basis to make progress in solving this problem. Following our previous work on oxidation and protonation states of the Mn-cluster, in this work, quantum chemical/electrostatic calculations were performed in order to estimate the pKa of different protons of relevant groups and atoms of the Mn-cluster such as W2, O4, O5 and His337. In broad agreement with previous experimental and theoretical work, our data suggest that W2 and His337 are likely to be in hydroxyl and neutral form, respectively, O5 and O4 to be unprotonated. This article is part of a Special Issue entitled: Photosynthesis Research for Sustainability: Keys to Produce Clean Energy.

Weak sulfur-sulfur interactions between chemically-identical atoms

Abstract Experimentally-known sulfur-sulfur distances shorter than the sum of van der Waals radii and involving two chemically-identical sulfur atoms are examined at several levels of theory (BP86/6-31G** to CCSD(T)/6-311+G**). None of the theoretical methods predict an attractive interaction from an energetic point of view, even though molecular orbitals stretching between the two sulfur atoms have been identified. Most likely, if there is indeed an attractive interaction force between chemically identical sulfur atoms, its value is comparable to the accuracy of the methods employed here — implying an attractive interaction below 1 kcal/mol. The investigation includes some simple models of 1,6,12,17-tetrathiacyclodocosa-2,4,13,15-tetrayne which was previously shown to have an S—S interaction involving two chemically-identical atoms. Attractive interactions calculated for these latter models are shown to arise from S—HC weak bonding, with the S—S interaction being again repulsive. Graphical abstract

Non-canonical interactions of porphyrins in porphyrin-containing proteins

In this study we have described the noncanonical interactions between the porphyrin ring and the protein part of porphyrin-containing proteins to better understand their stabilizing role. The analysis reported in this study shows that the predominant type of non-canonical interactions at porphyrins are CH····O and CH····N interactions, with a small percentage of CH···π and noncanonical interactions involving sulfur atoms. The majority of non-canonical interactions are formed from side-chains of charged and polar amino acids, whereas backbone groups are not frequently involved. The main-chain noncanonical interactions might be slightly more linear than the side-chain interactions, and they have somewhat shorter median distances. The analysis, performed in this study, shows that about 44% of the total interactions in the dataset are involved in the formation of multiple (furcated) noncanonical interactions. The high number of porphyrin-water interactions show importance of the inclusion of solvent in protein-ligand interaction studies. Furthermore, in the present study we have observed that stabilization centers are composed predominantly from nonpolar amino acid residues. Amino acids deployed in the environment of porphyrin rings are deposited in helices and coils. The results from this study might be used for structure-based porphyrin protein prediction and as scaffolds for future porphyrin-containing protein design.

Pyridine encapsulated hyperbranched polymers as mimetic models of haeme containing proteins, that also provide interesting and unusual porphyrin-ligand geometries

This communication describes the use of non-covalent chemistry to construct recyclable porphyrin cored HBPs. The non-covalent design allows the polymeric backbone to be rescued and reused after porphyrin degradation. The steric environment within the polymeric encapsulated ligand notably affected the porphyrin coordination geometry.

Orientation and stereodynamic paths of planar monodentate ligands in square planar nickel N2S complexes

The well-established presence of histidine donors in binding sites of Ni-containing biomolecules prompts the study of orientational preference and stereodynamic nature of flat monodentate ligands (L = imidazoles, pyridine and an N-heterocyclic carbene) bound to planar N(2)SNi moieties. Square planar [N(2)SNiL](n+) complexes are accessed through bridge-splitting reactions of dimeric, thiolate-S bridged [N(2)SNi](2) complexes. The solid state molecular structures of three mononuclear products, and three monothiolate bridged dinickel complexes, reveal that the plane of the added monodentate ligand orients largely orthogonal to the N(2)SNiL square plane. Variable temperature (1)H NMR characterization of dynamic processes and ground state isomer ratios of imidazole complexes in their stopped exchange limiting spectra, readily correlate with density functional theory (DFT)-guided interpretation of Ni-L rotational activation barriers. Full DFT characterization finds Ni-L bond lengthening as well as a tetrahedral twist distortion in the transition state, reaching a maximum in the NHC complex, and relating mainly to the steric hindrance derived both from the ligand and the binding pocket. In the case of the imidazole ligands a minor electronic contribution derives from intramolecular electrostatic interactions (imidazole C-2 C-H(delta+)- - S(delta-) interaction). Computational studies find this donor-acceptor interaction is magnified in O-analogues, predicting coplanar arrangements in the ground state of N(2)ON(imid)Ni complexes.

Understanding properties of cofactors in proteins: redox potentials of synthetic cytochromes b

Haehnel et al. synthesized 399 different artificial cytochrome b (aCb) models. They consist of a template-assisted four-helix bundle with one embedded heme group. Their redox potentials were measured and cover the range from -148 to -89 mV. No crystal structures of these aCb are available. Therefore, we use the chemical composition and general structural principles to generate atomic coordinates of 31 of these aCb mutants, which are chosen to cover the whole interval of redox potentials. We start by modeling the coordinates of one aCb from scratch. Its structure remains stable after energy minimization and during molecular dynamics simulation over 2 ns. Based on this structure, coordinates of the other 30 aCb mutants are modeled. The calculated redox potentials for these 31 aCb agree within 10 mV with the experimental values in terms of root mean square deviation. Analysis of the dependence of heme redox potential on protein environment shows that the shifts in redox potentials relative to the model systems in water are due to the low-dielectric medium of the protein and the protonation states of the heme propionic acid groups, which are influenced by the surrounding amino acids. Alternatively, we perform a blind prediction of the same redox potentials using an empirical approach based on a linear scoring function and reach a similar accuracy. Both methods are useful to understand and predict heme redox potentials. Based on the modeled structure we can understand the detailed structural differences between aCb mutants that give rise to shifts in heme redox potential. On the other hand, one can explore the correlation between sequence variations and aCb redox potentials more directly and on much larger scale using the empirical prediction scheme, which--thanks to its simplicity--is much faster.

Imidazole-containing (N3S)-Ni(II) complexes relating to nickel containing biomolecules

Dimeric (N(2)S)Ni complexes and the monomeric N(2)S(2) bismercaptodiazacycloheptane nickel complex, (bme-dach)Ni, serve as precursors to two N(2)-, N'-/ S- complexes where N(2) = diazacycloheptane, N' = imidazole and S = thiolate. As rare examples of nickel complexes containing a mixed thiolate/imidazole ligand set, these complexes are characterized by X-ray diffraction, UV/vis, and variable temperature (1)H NMR spectroscopies, and electrochemistry. Density functional theory computations relate the orientation of the imidazole with respect to the N(2)N'SNi square plane to the VT NMR observed fluxionality and activation parameters. The superoxide dismutase activity of the imidazole complexes was investigated by the nitroblue tetrazolium assay.

Ligand preference and orientation in b- and c-type heme-binding proteins

Hemes are often incorporated into designed proteins. The importance of the heme ligand type and its orientation is still a matter of debate. Here, heme ligands and ligand orientation were investigated using a nonredundant (87 structures) and a redundant (1503 structures) set of structures to compare and contrast design features of natural b- and c-type heme-binding proteins. Histidine is the most common ligand. Marked differences in ligation motifs between b- and c-type hemes are higher occurrence of His-Met in c-type heme binding motifs (16.4% vs. 1.4%) and higher occurrence of exchangeable, small molecules in b-type heme binding motifs (67.6% vs. 9.9%). Histidine ligands that are part of the c-type CXXCH heme-binding motif show a distinct asymmetric distribution of orientation. They tend to point between either the heme propionates or between the NA and NB heme nitrogens. Molecular mechanics calculations show that this asymmetry is due to the bonded constraints of the covalent attachment between the heme and the protein. In contrast, the orientations of b-type hemes histidine ligands are found evenly distributed with no preference. Observed histidine heme ligand orientations show no dominating influence of electrostatic interactions between the heme propionates and the ligands. Furthermore, ligands in bis-His hemes are found more frequently perpendicular rather than parallel to each other. These correlations support energetic constraints on ligands that can be used in designing proteins.

XH/pi interactions with the pi system of porphyrin ring in porphyrin-containing proteins

Searching structures of porphyrin-containing proteins from the Protein Data Bank revealed that the pi system of every porphyrin ring is involved in XH/pi interactions, with most of the porphyrins having several interactions. Both five-membered pyrrole rings and six-membered chelate rings are involved in XH/pi interactions; the number of interactions with five-membered rings is larger than the number of interactions with six-membered rings. We found interactions with C-H and N-H groups as hydrogen-atom donors; however, the number of CH/pi interactions is much larger than the number of NH/pi interactions. The amino acids involved in the interactions show a high conservation score. Our results that every porphyrin is involved in XH/pi interactions and that amino acids involved in these interactions are highly conserved demonstrate that XH/pi interactions play an important role in porphyrin-protein stability.

Ab initio molecular dynamics of heme in cytochrome c

Ab initio molecular dynamics (AIMD) calculations, based on the Car-Parrinello method, have been carried out for three models of heme c that is present in cytochrome c. Both the reduced (Fe(II)) and oxidized (Fe(III)) forms have been analyzed. The simplest models (1R and 1O, respectively) consist of a unsubstituted porphyrin (with no side chains) and two axially coordinated imidazole and ethylmethylthioether ligands. Density functional theory optimizations of these models confirm the basic electronic features and are the starting point for building more complex derivatives. AIMD simulations were performed after reaching the thermal stability at T = 300 K. The evolution of the Fe-L(ax) bond strengths is examined together with the relative rotations of the imidazole and methionine about the axial vector, which appear rather independent from each other. The next models (2R and 2O) contain side chains at the heme to better simulate the actual active site. It is observed that two adjacent propionate groups induce some important effects. The axial Fe-Sdelta bond is only weakened in 2R but is definitely cleaved in the oxidized species 2O. Also the mobility of the Im ligand seems to be reduced by the formation of a strong hydrogen bond that involves the Im Ndelta1-Hdelta1 bond and one carboxylate group. In 2O the interaction becomes so strong that a proton transfer occurs and the propionic acid is formed. Finally, the models 3 include a free N-methyl-acetamide molecule to mimic a portion of the protein backbone. This influences the orientation of carboxylate groups and limits the amount of their hydrogen bonding with the Im ligand. Residual electrostatic interactions are maintained, which are still able to modulate the dissociation of the methionine from the heme.

Electron delocalization mediates the metal-dependent capacity for CH/pi interactions of acetylacetonato chelates

CH/pi interactions between the coordinated acetylacetonato ligand and phenyl rings were analyzed in the crystal structures from the Cambridge Structural Database and by quantum chemical calculations. The acetylacetonato ligand may engage in two types of interactions: it can be hydrogen atom donor or acceptor. The analysis of crystal structures and calculations show that interactions with the acetylacetonato ligand acting as hydrogen atom donor depend on the metal in an acetylacetonato chelate ring; the chelate rings with soft metals make stronger interactions. The same trend was not observed in the interactions where the acetylacetonato chelate ring acts as the hydrogen atom acceptor.

Molecular statistics of cytochrome c: structural plasticity and molecular environment

Nuclear magnetic resonance experiments performed on yeast mitochondrial cytochrome c (Cytc), a paradigmatic electron transfer protein, reveal that the two oxidation states have similar structures, but different mobility: despite the few structural differences compared with the reduced form, the oxidized form displays a larger unfolding propensity. Molecular dynamics simulations performed on both NMR reduced and NMR oxidized forms show that the reduced form has a larger solvent-accessible surface area (SASA). Starting from this observation, a molecular statistical approach was then applied in order to correlate the molecular surface to molecular mobility. Simulations started from biased initial conditions corresponding to different molecular sizes were combined with the maximal constrained entropy method. The NMR structure of oxidized Cytc is more suited to expose a smaller SASA than the NMR structure of the reduced form, but the accessible conformational landscape at 300 K around the NMR oxidized structure is flatter than for the NMR reduced structure. Protein configurations of smaller SASA and size display larger plasticity when they resemble the NMR oxidized structure, whereas they are more rigid when they resemble the NMR reduced structure. Implications of the results for the protein properties during its functional process are discussed.

Factors influencing the energetics of electron and proton transfers in proteins. What can be learned from calculations

A protein structure should provide the information needed to understand its observed properties. Significant progress has been made in developing accurate calculations of acid/base and oxidation/reduction reactions in proteins. Current methods and their strengths and weaknesses are discussed. The distribution and calculated ionization states in a survey of proteins is described, showing that a significant minority of acidic and basic residues are buried in the protein and that most of these remain ionized. The electrochemistry of heme and quinones are considered. Proton transfers in bacteriorhodopsin and coupled electron and proton transfers in photosynthetic reaction centers, 5-coordinate heme binding proteins and cytochrome c oxidase are highlighted as systems where calculations have provided insight into the reaction mechanism.

One-Electron Reduction Potential for Oxygen- and Sulfur-Centered Organic Radicals in Protic and Aprotic Solvents

We estimated one-electron reduction potentials of redox-active organic molecules for a spectrum of eight different functional groups (phenoxyl, p-benzoquinone, phenylthiyl, p-benzodithiyl, carboxyl, benzoyloxyl, carbthiyl, and benzoylthiyl) in protic (water) and aprotic (acetonitrile, N,N-dimethylacetamide) solvents. Electron affinities (EA) were evaluated in a vacuum with high level quantum chemical methods using Gaussian3-MP2 (G3MP2) and Becke 3 Lee, Yang, and Parr functional B3LYP with aug-cc-pVTZ basis set. To evaluate one-electron redox potentials, gas-phase free energies were combined with solvation energies obtained in a two-step computational approach. First, atomic partial charges were determined in a vacuum by the quantum chemical method B3LYP/6-31G(d,p). Second, solvation energies were determined, solving the Poisson equation with these atomic partial charges. Redox potentials computed this way, compared to experimental data for the 21 considered organic compounds in different solvents, yielded overall root-mean-square deviations of 0.058 and 0.131 V using G3MP2 or B3LYP to compute electronic energies, respectively, while B3LYP/6-31G(d,p) was used to compute solvation energies.