Noninnocent effect of axial ligand on the heme degradation process: a theoretical approach to hydrolysis pathway of verdoheme to biliverdin

Conversion of iron(II) verdoheme to iron(II) biliverdin in the presence of hydroxyl ion as a nucleophile and imidazole, pyridine, water, hydroxyl, cyanide, phenolate, chloride, thiolate and imidazolate as axial ligands was investigated using the B3LYP method and the 6-31G basis set. In the five-coordinated pathway the reactants and products are in the ground triplet state. In this path, hydroxyl ion directly attacks the macrocycle. The exothermic energy for addition of hydroxyl ion to iron(II) verdoheme with various ligands is 169.55, 166.34 and 164 kcal mol⁻¹ for water, pyridine and imidazole, energies which are around 30–60 kcal mol⁻¹ more exothermic than those for the other axial ligands used in this study. Therefore, imidazole, water and pyridine axial ligands can facilitate hydrolytic cleavage of iron(II) verdoheme to form open-chained helical iron(II) biliverdin complexes. The activation barrier for the conversion of iron(II) verdoheme hydroxyl species to the iron(II) biliverdin complex is estimated to be 5.2, 4.2, 4.35, 13.76 and 14.05 kcal mol⁻¹ for imidazole, water, cyanide, thiolate and imidazolate, respectively. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

A theoretical study of myoglobin working as a nitric oxide scavenger

The mechanism for the reaction between nitric oxide (NO) and O₂ bound to the heme iron of myoglobin (Mb), including the following isomerization to nitrate, has been investigated using hybrid density functional theory (B3LYP). Myoglobin working as a NO scavenger could be of importance, since NO reversibly inhibits the terminal enzyme in the respiration chain, cytochrome c oxidase. The concentration of NO in the cell will thus affect the respiration and thereby the synthesis of ATP. The calculations show that the reaction between NO and the heme-bound O₂ gives a peroxynitrite intermediate whose O–O bond undergoes a homolytic cleavage, forming a NO₂ radical and myoglobin in the oxo-ferryl state. The NO₂ radical then recombines with the oxo-ferryl, forming heme-bound nitrate. Nine different models have been used in the present study to examine the effect on the reaction both by the presence and the protonation state of the distal His64, and by the surroundings of the proximal His93. The barriers going from the oxy-Mb and nitric oxide reactant to the peroxynitrite intermediate and further to the oxo-ferryl and NO₂ radical are around 10 and 7 kcal/mol, respectively. Forming the product, nitrate bound to the heme iron has a barrier of less than ~7 kcal/mol. The overall reaction going from a free nitric oxide and oxy-Mb to the heme bound nitrate is exergonic by more than 30 kcal/mol.     

Theoretical study on the ground state intramolecular proton transfer (IPT) and solvation effect in two Schiff bases formed by 2-aminopyridine with 2-hydroxy-1- naphthaldehyde and 2-hydroxy salicylaldehyde

The tautomerization mechanism the isolated and monohydrated forms of two Schiff bases 1 and 2, and the effect of solvation on the proton transfer from enol-imine form to the keto-enamine form have been investigated using the B3LYP hybrid density functional method at the 6-31G** basis set level. The barrier heights for H₂O-assisted reactions are significantly lower than that of unassisted tautomerization reaction in the gas phase. Nonspecific solvent effects have also been taken into account by using the continuum model (IPCM) of four different solvent. The tautomerization energies and the potential energy barriers are decreased by increasing solvent polarity. Figure The tautomerization mechanism the isolated and monohydrated forms of two Schiff bases 1 and 2, and the effect of solvation on the proton transfer from enol-imine form to the keto-enamine form have been investigated using the B3LYP hybrid density functional method at the 6-31G** basis set level     

Reaction mechanism of the binuclear zinc enzyme glyoxalase II - A theoretical study

The glyoxalase system catalyzes the conversion of toxic methylglyoxal to nontoxic d-lactic acid using glutathione (GSH) as a coenzyme. Glyoxalase II (GlxII) is a binuclear Zn enzyme that catalyzes the second step of this conversion, namely the hydrolysis of S-d-lactoylglutathione, which is the product of the Glyoxalase I (GlxI) reaction. In this paper we use density functional theory method to investigate the reaction mechanism of GlxII. A model of the active site is constructed on the basis of the X-ray crystal structure of the native enzyme. Stationary points along the reaction pathway are optimized and the potential energy surface for the reaction is calculated. The calculations give strong support to the previously proposed mechanism. It is found that the bridging hydroxide is capable of performing nucleophilic attack at the substrate carbonyl to form a tetrahedral intermediate. This step is followed by a proton transfer from the bridging oxygen to Asp58 and finally C-S bond cleavage. The roles of the two zinc ions in the reaction mechanism are analyzed. Zn2 is found to stabilize the charge of tetrahedral intermediate thereby lowering the barrier for the nucleophilic attack, while Zn1 stabilizes the charge of the thiolate product, thereby facilitating the C-S bond cleavage. Finally, the energies involved in the product release and active-site regeneration are estimated and a new possible mechanism is suggested.     

A theoretical study of zinc(II) interactions with amino acid models and peptide fragments

Density functional theory calculations have been employed to study the interaction between the Zn²⁺ ion and some standard amino acid models. The highest affinities towards the Zn²⁺ ion are predicted for serine, cysteine, and histidine. Relatively high affinities are reported also for proline and glutamate/aspartate residues. It was found that the zinc complexes with cysteine adopt a tetrahedral conformation. Conversely, complexes with one or two histidine moieties remain in an octahedral geometry, while those with three or more histidine groups adopt a square-planar geometry. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Novel derivatives of ansa-titanocenes procured from 6-phenylfulvene: a combined experimental and theoretical study

The previously prepared trans-[(1,2-diphenyl-1,2-dicyclopentadienyl)ethanediyl] titanium(IV) dichloride, [1,2-(Ph)₂C₂H₂{η⁵-C₅H₄}₂]Ti(Cl)₂, was synthesised using an alternative procedure, from which its crystal structure was determined. Using this compound, a variety of other ansa-titanocene derivatives were synthesised by replacement of the chloride ligands with selected substituents. Thus RTi(X)(Y) systems where R=1,2-(Ph)₂C₂H₂η⁵-C₅H₄₂; X=Y=CH₃; X=CH₃, Y=Cl; X=Y=NCS; X=Y=NCO; X=Y=OPh and (X/Y)=O have been synthesised and characterised. DFT calculations were performed on the complexes trans-[(1,2-diphenyl-1,2-dicyclopentadienyl)-ethanediyl] titanium(IV) dichloride, bis-(6,6-diphenylfulvene)titanium and bis-(6,6-diphenylfulvene)iron. This demonstrated the role that the metal centre plays in their formation, generating either an ansa-metallocene, a ‘tucked in’ fulvene complex or a metallocene coordinating fulvene anions.     

Ligation of Aza bases to the AgF<Subscript>2</Subscript> molecule: a theoretical study

We have analyzed the electronic structure and chemical bonding for molecular adducts of the Ag(II)F₂ molecule with various aza Lewis bases including ammonia, nitriles, secondary amines, and their derivatives exhibiting various degrees of fluorination. Density functional theory calculations indicate that a progressive shift occurs of the spin density from the Ag center towards the coordinating nitrogen atoms of aza ligands, as the ligation energy increases. Chemistry of Ag(II) might be extended with little effort beyond the known aza connections, to include nitriles, perfluorinated nitriles and perfluorinated amines. Figure Properties of a variety of novel adducts of the AgF₂ molecule with two aza bases (L), possible precursors of the AgF₂L₂ extended solids, were assessed by the DFT calculations Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. This work is dedicated to memory of Wojciech Ochmański, unforgettable person, good-hearted man, whose craftsmanship in work was second-to-none.     

A gas-phase study of the gold-catalyzed coupling of alkynes and alcohols

Mass spectrometric methods are used for an investigation of the C-O bond formation between alkynes and alcohols mediated by cationic gold(I) complexes. In the condensed-phase catalysis, an efficient coupling to the corresponding enol ethers has been observed. In the gas phase, however, kinetic and entropic restrictions are too large for both bare Au⁺ as well as (CH₃)₃PAu⁺ to permit bond formation in strictly bimolecular collisions occurring at thermal energies. The experimental findings are complemented by a theoretical study of the Au⁺/C₂H₂/CH₃OH system using density functional theory.     

A combined theoretical and spectroscopic study of 4,6-di-O-acetyl-2,3-dideoxy-d-erythro-hex-2-enopyranosyl sulfamide: a novel glycosyl carbonic anhydrase IX inhibitor

The novel 4,6-di-O-acetyl-2,3-dideoxy-d-erythro-hex-2-enopyranosyl sulfamide, which exhibits selectivity for inhibiting isoform IX of carbonic anhydrase as overexpressed in many tumors, has been investigated from a combined theoretical and spectroscopic point of view. The conformational study of the compound shows that the α-anomeric form is more stable than the β-anomeric form from a thermodynamic point of view after including solvent effects. This fact suggests that the synthesis reaction could take place mainly under thermodynamic control as the main experimental product is the α-anomeric form of the sulfamide. Calculated α/β ratio is about 95:5, in excellent agreement with experimental data. Optimized geometries of the α-anomeric form agree quite well with crystallographic data. The inclusion of a solvent has negligible effects on the conformations. A detailed analysis of some geometric parameters shed light into the conformational behavior of the sulfamide in terms of both exo- and endo-anomeric effects and antiperiplanar relationships. Natural bond orbital calculations confirm those findings. Several intramolecular hydrogen bonds, characterized through the Atoms-in-Molecules theory, were found in the stable conformers. They, however, seem to play no relevant role in determining the relative stability of α conformers with respect to the β ones. Calculated ¹H and ¹³C NMR chemical shifts support previous findings concerning configuration and conformation assignments of the title sulfamide. The IR spectrum of the compound is recorded and reported for the first time and the assignment of some of the most important bands is accomplished with the aid of calculated harmonic vibrational frequencies.     

A theoretical study of spin states in Ni-S<Subscript>4</Subscript> complexes and models of the [NiFe] hydrogenase active site

We have applied density functional theory, using both pure (BP86) and hybrid (B3LYP and B3LYP*) functionals, to investigate structural parameters and reaction energies for nickel(II)-sulfur coordination compounds, as well as for small cluster models of the Ni-SI and Ni-R redox state of [NiFe] hydrogenases. Results obtained investigating experimentally well-characterized complexes show that BP86 is well suited to describe the structural features of this class of compounds. However, the singlet-triplet energy splitting and even the computed ground state are strongly dependent on the applied functional. Results for the cluster models of [NiFe] hydrogenases lead to the conclusion that in the reduced protein structures characterized by X-ray diffraction a hydride bridges the two metal centres. The energy splitting of the singlet and triplet states in Ni-R and Ni-SI models is calculated to be very small and may be overcome at room temperature to allow a spin crossover. Moreover, the relative stability of the Ni-SI and Ni-R structures adopted in the present investigation is fully compatible with their involvement in the reversible heterolytic cleavage of H(2).     

Structural properties of sH hydrate: a DFT study of anisotropy and equation of state

ABSTRACT

Structure-H (sH) hydrate is one of the canonical gas hydrates with significant potential applications and scarce characterised material properties despite the wide knowledge available on other gas hydrates. In this work we characterise some of the important physical properties of this hydrate at the atomistic level using Density Functional Theory. Two exchange-correlation functionals (revPBE and DRSLL) were used to simulate six sH hydrate systems encapsulating neohexane and different help gas molecules. The important role of dispersion forces is quantified. The density and isothermal bulk modulus of sH hydrate are higher when dispersion interactions are considered. The presence of those interactions imposes a direct relationship between the hydrate density and its bulk modulus, while their absence reveals the bulk modulus dependency on hydrogen bond density. Anisotropy is a distinguishing feature of this hydrate in distinction to nearly isotropic sI and sII hydrates. Structure-H hydrate experiences a compressional anisotropy in which the a-lattice and the c-lattice constants respond differently to applied pressure showing less compressibility along the c-axis. This compressional anisotropy was found dependant on the chemistry of help gas molecules. Taken together, these property characterisation results and analysis are a significant and novel contribution to the material physics of sH hydrates.     

Density functional theoretical study of the structure, bonding and electrochemistry of N2S2 nitrosyl complexes of iron and cobalt

Density functional calculations are used to study the structure, bonding and properties of the N₂S₂ nitrosyl complexes (bme-pda)M(NO) and (bme-dach)M(NO) M = Fe, Co. The complexes studied in this work have been previously shown to exhibit unusual structural and electrochemical properties. Calculated optimal geometries and vibrational frequencies are compared to experimental data and found to be in excellent agreement. A very small energy barrier for rotation of the NO ligand relative to the S-donor atoms of the N₂S₂ ligand is calculated. The similarity in reduction potentials of Fe and Co complexes is investigated using molecular orbital analysis and calculation of absolute reduction potentials. It is concluded that reduction of the Fe complexes may produce triplet-state products. Such products are determined to be of lower energy than singlet-state products, and also consistent with the observed reduction potentials of the Fe and Co complexes.     

A theoretical study on H-bonding interactions in maleic acid: calculated 17O, 1H NMR parameters and QTAIM analysis

A theoretical study was carried out to examine intra- and intermolecular hydrogen bond (HB) properties in crystalline maleic acid (MA). We investigated geometries, ¹⁷O and ¹H nuclear magnetic resonance parameters of various MA clusters by means of M06 and B3LYP functionals using recently developed Jensen's polarisation-consistent basis sets, pcJ-n (n = 0, 1, 2, 3). Our results reveal that the calculated chemical shift isotropy, δ_iso, at the sites of ¹⁷O and ¹H nuclei depends markedly on the size of the basis set. Overall, convergence of the pcJ-n series is rather similar for both B3LYP and M06 functionals. An increase of δ_iso(¹⁷O) and δ_iso(¹H) in going from the pcJ-1 to the pcJ-3 basis set is a typical feature of the (MA)_1–3 clusters. The quantum theory of atoms in molecules (QTAIM) and energy decomposition analyses were also used to elucidate the interaction characteristics in the MA H-bonded network. According to QTAIM results obtained, it is concluded that strong HBs are more covalent in nature and weak HBs are mainly electrostatic interactions.     

Synthesis and theoretical analysis of the structures and bonding in five new neutral square planar bis[2,4-di(aryl)-1,3,5-triazapentadienato]nickel(II) complexes

Solvothermal reactions in methanol of nickel acetate tetrahydrate, Ni(OAc)₂ · 4H₂O, with benzonitrile derivatives NC(C₆H₄)X, where X is one of the electron withdrawing substituents -CN, -NO₂, or -CF₃, located at the m- or p-positions relative to -CN, yield complexes of the general formula Ni{HNC(R)-NC(R)-NH}₂. More specifically, 3-nitrobenzonitrile, 4-nitrobenzonitrile, 1,3-dicyanobenzene, 1,4-dicyanobenzene, and ααα-trifluoro-p-toluonitrile are found to react with Ni(OAc)₂ · 4H₂O to yield Ni{HNC(R)-NC(R)-NH}₂, where R = 3-(NO₂)C₆H₄, 4-(NO₂)C₆H₄, 3-(CN)C₆H₄, 4-(CN)C₆H₄, or 4-(CF₃)C₆H₄, respectively. Analogous reactions of nitriles lacking electron withdrawing groups do not occur under similar conditions. Solid-state structures have been determined for the complexes with p-NO₂, p-CN, and p-CF₃ substituents on the phenyl rings. In addition, we describe density functional theory (DFT) and natural bonding orbital theory (NBO) studies on a simplified analog of these compounds, aimed at understanding their molecular bonding. It is shown that the new compounds for which solid-state structures have been determined are model examples of coordination compounds containing robust ω-bonds.     

DFT: a dynamic study of the interaction of ethanol and methanol with platinum

The interaction between alcohol molecules and platinum (Pt) was studied using molecular dynamics (MD; Born-Oppenheimer method). Alcohol molecules like ethanol and methanol present a similar molecular structure, with a methyl group (CH₃) at one end and a fragment of hydroxyl (OH) at the other. This fact generates two orientations that are considered in the interaction with Pt. The MD calculation results for these two orientations indicate a preferential orientation due to energy interactions. A plausible reaction mechanism that takes into account the interaction between Pt and alcohol is presented. The charge transference obtained from the Pt–alcohol interaction was also analyzed. The energy for the two orientations was calculated by indicating the preferential orientation. The methyl and hydroxyl groups are involved in heterolytic breakage of hydrogen bonds, joined to a carbon atom in the former and to an oxygen atom in the latter; however, the methyl group reaction seems to be the most important.     

Structural and electronic properties of diazonium functionalized (4, 4) single walled carbon nanotube: an ab initio study

In this work, ab initio density functional theory (DFT) calculations are performed to study the structural and electronic properties of diazonium reagent functionalized (4, 4) single-walled carbon nanotube (SWCNT). We find the aryl group covalently bonds with SWCNT and prefers to be perpendicular to the side wall of nanotube. It has a rotational barrier of 0.35 eV around the formed aryl-tube bond axis and should be thermodynamically stable at room temperature. Additionally, new peaks appeared around the Fermi energy in the density of state (DOS) due to the weak band dispersion. Increasing of the coverage of the functional group will result in significant upshift of the Fermi level.     

On the oxidation of alkyl and aryl sulfides by [(Me3TACN)MnO(OH)2]: A density functional study

Density functional theory suggests that the formal 2-electron oxidation of sulfides, RR′S, to sulfoxides, by the model Mn^VO catalyst, [(TACN)Mn^V O(OH)₂]⁺, proceeds in two quite distinct 1-electron steps. Transfer of the first electron is barrierless and generates a sulfur radical cation, antiferromagnetically coupled to a Mn^IV centre via a covalent μ-oxo bridge. The second electron-transfer step is accompanied by migration of the oxygen atom to the sulfur centre, and is rate-determining. The absence of a barrier in the first step, where a sulfur radical is formed, means that the presence of electron-donating or withdrawing substituents on the sulfide has only a minor impact on the rate of reaction.     

Investigation of intermolecular hydrogen bond interactions in crystalline l-Cysteine by DFT calculations of the oxygen-17, nitrogen-14, and hydrogen-2 EFG tensors and AIM analysis

A systematic computational study is carried out to investigate hydrogen bond (HB) interactions in the real crystalline structures of l-Cysteine at 30 and 298 K by density functional theory (DFT) calculations of electric field gradient (EFG) tensors at the sites of O-17, N-14, and H-2 nuclei. One-molecule (monomer) and nine-molecule (cluster) models of l-Cysteine are created by available crystal coordinates at both temperatures and the EFG tensors are calculated for both models to indicate the effect of HB interactions on the tensors. The calculated EFG tensors at the level of B3LYP and B3PW91 DFT methods and 6-311++G?? and cc-pVTZ basis sets are converted to those experimentally measurable nuclear quadrupole resonance (NQR) parameters i.e. quadrupole coupling constants (qcc) and asymmetry parameters (η_Q). The evaluated NQR parameters reveal that the EFG tensors of ¹⁷O, ¹⁴N, and ²H are influenced and show particular trends from monomer to the target molecule in the cluster due to the contribution of target molecule to classic N–H…O, and non-classic S–H…O and S–H…S types of HB interactions. On the other hand, atoms in molecules (AIM) analyses confirm the presence of HB interactions and rationalize the observed EFG trends. The results indicate different contribution of various nuclei to HB interactions in the cluster where O2 and N1 have major contributions. The EFG tensors as well as AIM analysis at the H6 site show that the N1-H6…O2 HB undergoes a significant change from 30 to 298 K where changes in other N–H…O interactions are almost negligible. There is a good agreement between the calculated ¹⁴N NQR parameters and reported experimental data.     

Insights into a highly conserved network of hydrogen bonds in the agonist binding site of nicotinic acetylcholine receptors: A structural and theoretical study

Structural and theoretical studies on the geometrical features of a hydrogen‐bond network occurring in the binding site of nicotinic acetylcholine receptors (nAChRs) and composed of interconnected WxPD (Trp‐x‐Pro‐Asp) and SWyz (Ser‐Trp‐yz) sequences from loops A and B, respectively, have been carried out. Multiple sequence alignments using as template the sequence of the apoform of Aplysia californica acetylcholine binding protein (Ac‐AChBP) show the strict conservation of serine and tryptophan residues of the loop B SWyz sequence. Considering a sample of 19 high resolution AChBP structures, the strong conformational preferences of the key tryptophan residue has been pointing out, whatever the form, free or bounded, of AChBP. The geometry of the motif hydrogen‐bond network has been characterized through the analyses of seven distances. The robustness of the various hydrogen‐bond interactions is pointed out, the one involving the aspartate carboxylate group and the serine residue being the shortest of the network. The role of a cooperative effect involving a NH(His145)…OH (Ser142) hydrogen bond is highlighted. Density functional theory calculations on several simplified models based on the motif hydrogen‐bond network allow probing the importance of the various hydrogen‐bond interactions. The removal of the Ser142 hydroxyl group induces strong structural rearrangements, in agreement with the structural observations. Molecular electrostatic potential calculations on model systems highlight the importance of a cooperative effect in the whole hydrogen‐bond network. More precisely, the key role of the Ser142 hydroxyl group, involved in several hydrogen bonds, is underlined. Proteins 2014; 82:2303–2317. © 2014 Wiley Periodicals, Inc.