Biomedical subjects

Ulf Ryde

Publications and source records attributed to Ulf Ryde.

At least 19 recordsLinked to original sources

PubMed2007

Performance of density functionals for first row transition metal systems.

This article investigates the performance of five commonly used density functionals, B3LYP, BP86, PBE0, PBE, and BLYP, for studying diatomic molecules consisting of a first row transition metal bonded to H, F, Cl, Br, N, C, O, or S. Results have been compared with experiment wherever possible. Open-shell configurations are found more often in the order PBE0>B3LYP>PBE approximately BP86>BLYP. However, on average, 58 of 63 spins are correctly predicted by any functional, with only small differences. BP86 and PBE are slightly better for obtaining geometries, with errors of only 0.020 A. Hybrid functionals tend to overestimate bond lengths by a few picometers and underestimate bond strengths by favoring open shells. Nonhybrid functionals usually overestimate bond energies. All functionals exhibit similar errors in bond energies, between 42 and 53 kJmol. Late transition metals are found to be better modeled by hybrid functionals, whereas nonhybrid functionals tend to have less of a preference. There are systematic errors in predicting certain properties that could be remedied. BLYP performs the best for ionization potentials studied here, PBE0 the worst. In other cases, errors are similar. Finally, there is a clear tendency for hybrid functionals to give larger dipole moments than nonhybrid functionals. These observations may be helpful in choosing and improving existing functionals for tasks involving transition metals, and for designing new, improved functionals.

Journal Article↗

PubMed2006

Multireference ab initio calculations on reaction intermediates of the multicopper oxidases.

The multicopper oxidases (MCOs) couple the four-electron reduction of dioxygen to water with four one-electron oxidations of various substrates. Extensive spectroscopic studies have identified several intermediates in the MCO catalytic cycle, but they have not been able to settle the structures of three of the intermediates, viz. the native intermediate (NI), the peroxy intermediate (PI), and the peroxy adduct (PA). The suggested structures have been further refined and characterized by quantum mechanical/molecular mechanical (QM/MM) calculations. In this paper, we try to establish a direct link between theory and experiment, by calculating spectroscopic parameters for these intermediates using multireference wave functions from the multistate CASPT2 and MRDDCI2 methods. Thereby, we have been able to reproduce low-spin ground states (S = 0 or S = 1/2) for all the MCO intermediates, as well as a low-lying (approximately 150 cm-1) doublet state and a doublet-quartet energy gap of approximately 780 cm-1 for the NI. Moreover, we reproduce the zero-field splitting (approximately 70 cm-1) of the ground 2E state in a D3 symmetric hydroxy-bridged trinuclear Cu(II) model of the NI and obtain a quantitatively correct quartet-doublet splitting (164 cm-1) for a mu3-oxo-bridged trinuclear Cu(II) cluster. All results support the suggestion that the NI has an O2- atom in the center of the trinuclear cluster, whereas both the PI and PA have an O22- ion in the center of the cluster, in agreement with the QM/MM results and spectroscopic measurements.

Copper↗

PubMed2006

Prediction of activation energies for hydrogen abstraction by cytochrome p450.

We have estimated the activation energy for hydrogen abstraction by compound I in cytochrome P450 for a diverse set of 24 small organic substrates using state-of-the-art density functional theory (B3LYP). We then show that these results can be reproduced by computationally less demanding methods, for example, by using small organic mimics of compound I with both B3LYP and the semiempirical AM1 method (mean absolute error of 3-4 kJ/mol) or by calculating the bond dissociation energy, without relaxation of the radical (B3LYP) or estimated from three-point fit to a Morse potential (AM1; errors of 4 and 5 kJ/mol, respectively). We can assign activation energies of 74, 61, 53, 47, and 30 kJ/mol to primary carbons, secondary/tertiary carbons, carbons with adjacent sp(2) or aromatic groups, ethers/thioethers, and amines, respectively, which gives a very simple and predictive model. Finally, some of the less demanding methods are applied to study the CYP3A4 metabolism of progesterone and dextromethorphan.

Algorithms↗

PubMed2006

Ligand affinities predicted with the MM/PBSA method: dependence on the simulation method and the force field.

The free energy of binding between avidin and seven biotin analogues has been calculated with the molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) method. We have studied how the force field and the method to generate geometries affect the calculated binding free energies. Four different force fields were compared, but we saw no significant difference in the results. However, mixing the force fields used for the geometry generation and energy calculations is not recommended. In the molecular dynamics simulations, explicit water molecules must be used, but the size of the simulated system and the boundary conditions are less important. In fact, nonperiodic simulations with a fixed protein outside a relatively small simulated system (18 A) seem to be a proper approach. The mean absolute error was 9-19 kJ/mol, with a standard error of 5-15 kJ/mol, which arises mainly from the entropy term.

Algorithms↗

PubMed2006

The role of axial ligands for the structure and function of chlorophylls.

We have studied the effect of axial ligation of chlorophyll and bacteriochlorophyll using density functional calculations. Eleven different axial ligands have been considered, including models of histidine, aspartate/glutamate, asparagine/glutamine, serine, tyrosine, methionine, water, the protein backbone, and phosphate. The native chlorophylls, as well as their cation and anion radical states and models of the reaction centres P680 and P700, have been studied and we have compared the geometries, binding energies, reduction potentials, and absorption spectra. Our results clearly show that the chlorophylls strongly prefer to be five-coordinate, in accordance with available crystal structures. The axial ligands decrease the reduction potentials, so they cannot explain the high potential of P680. They also redshift the Q band, but not enough to explain the occurrence of red chlorophylls. However, there is some relation between the axial ligands and their location in the various photosynthetic proteins. In particular, the intrinsic reduction potential of the second molecule in the electron transfer path is always lower than that of the third one, a feature that may prevent back-transfer of the electron.

Amino Acids↗

PubMed2006

Structural and photoluminescence properties of excited state intramolecular proton transfer capable compounds-potential emissive and electron transport materials.

Electronic factors influencing the photoluminescence properties and rates of excited state intramolecular proton transfer (ESIPT) reaction of o-hydroxy derivatives of 2,5-diphenyl-1,3,4-oxadiazole have been studied. The potential of these molecules as emissive and electron transport materials in designing improved organic light emitting diodes (OLEDs) has been studied by analyzing possible reasons for the unusually high Stokes shifts and ESIPT reaction rates. Time-dependent density functional theory (TDDFT) methods have been used to calculate the ground and excited state properties of the phototautomers that are the ESIPT reaction products. We study the relative effect of electron-withdrawing substituents on the proton-acceptor moiety and predict that the lowest ESIPT rate (1.9 x 10(11) s(-1)) is achieved with a dimethylamino substituent and that the Stokes shifts are around 11 000 cm(-1) for all three derivatives.

Journal Article↗

PubMed2006

Comparison of overlap-based models for approximating the exchange-repulsion energy.

Different ways of approximating the exchange-repulsion energy with a classical potential function have been investigated by fitting various expressions to the exact exchange-repulsion energy for a large set of molecular dimers. The expressions involve either the orbital overlap or the electron-density overlap. For comparison, the parameter-free exchange-repulsion model of the effective fragment potential (EFP) is also evaluated. The results show that exchange-repulsion energy is nearly proportional to both the orbital overlap and the density overlap. For accurate results, a distance-dependent correction is needed in both cases. If few parameters are desired, orbital overlap is superior to density overlap, but the fit to density overlap can be significantly improved by introducing more parameters. The EFP performs well, except for delocalized pi systems. However, an overlap expression with a few parameters seems to be slightly more accurate and considerably easier to approximate.

Chemistry, Physical↗

PubMed2006

Structure of reduced and oxidized manganese superoxide dismutase: a combined computational and experimental approach.

Manganese superoxide dismutases catalyze the disproportionation of the superoxide radical anion to molecular oxygen and hydrogen peroxide. Recently, atomic-resolution crystal structures of the reduced and oxidized enzymes have been reported. They show an active site with the manganese ion bound to one aspartate, three histidine residues, and a solvent molecule. In this paper, we combine crystallographic refinement with quantum mechanical methods to show that the solvent ligand is undoubtedly a water molecule in the reduced state. However, the putative oxidized structure is to a large extent reduced during data collection, so that it contains a mixture of the Mn2+ and Mn3+ structure. The crystal structures show that the Mn-bound solvent molecule accepts a hydrogen bond from the side chain of the conserved Gln-146 residue. If the solvent ligand is water, then this could lead to a steric clash, but it is avoided by the plane of water molecule forming an angle of 72 degrees to the Mn-O bond. Such a conformation is also found outside the enzyme, giving a minimal destabilization of the reduced state. We show by molecular dynamics simulations that the suggested Mn2+-H2O and Mn3+-OH- structures are stable. Moreover, we show that the superoxide substrate may bind both in the first coordination sphere of the Mn ion, opposite to the aspartate ligand, or in the second sphere, close to the conserved Tyr-34 and His-30 residues, approximately 5 A from Mn. However, the second-sphere structures are not stable in long molecular dynamics simulations. We see no difference in the coordination between the reduced and the oxidized states of the enzyme.

Hydrogen Bonding↗

PubMed2006

Implicit versus explicit solvent in free energy calculations of enzyme catalysis: Methyl transfer catalyzed by catechol O-methyltransferase.

We compare free energy calculations for the methyl transfer reaction catalyzed by catechol O-methyltransferase using the quantum mechanical/molecular mechanical free energy method with implicit and explicit solvents. An analogous methylation reaction in a solution is also studied. For the explicit solvent model, we use the three-point transferable intermolecular potential model, and for the implicit model, we use the generalized Born molecular volume model as implemented in CHARMM. We find that activation and reaction free energies calculated with the two models are very similar, despite some structural differences that exist. A significant change in the polarization of the environment occurs as the reaction proceeds. This is more pronounced for the reaction in a solution than for the enzymatic reaction. For the enzymatic reaction, most of the changes take place in the protein rather than in the solvent, and, hence, the benefit of having an instantaneous relaxation of the solvent degrees of freedom is less pronounced for the enzymatic reaction than for the reaction in a solution. This is a likely reason why energies of the enzyme reaction are less sensitive to the choice of atomic radii than are energies of the reaction in a solution.

Catalysis↗

PubMed2006

Structures of the high-valent metal-ion haem-oxygen intermediates in peroxidases, oxygenases and catalases.

Peroxidases, oxygenases and catalases have similar high-valent metal-ion intermediates in their respective reaction cycles. In this review, haem-based examples will be discussed. The intermediates of the haem-containing enzymes have been extensively studied for many years by different spectroscopic methods like UV-Vis, EPR (electron paramagnetic resonance), resonance Raman, Mössbauer and MCD (magnetic circular dichroism). The first crystal structure of one of these high-valent intermediates was on cytochrome c peroxidase in 1987. Since then, structures have appeared for catalases in 1996, 2002, 2003, putatively for cytochrome P450 in 2000, for myoglobin in 2002, for horseradish peroxidase in 2002 and for cytochrome c peroxidase again in 1994 and 2003. This review will focus on the most recent structural investigations for the different intermediates of these proteins. The structures of these intermediates will also be viewed in light of quantum mechanical (QM) calculations on haem models. In particular quantum refinement, which is a combination of QM calculations and crystallography, will be discussed. Only small structural changes accompany the generation of these intermediates. The crystal structures show that the compound I state, with a so called pi-cation radical on the haem group, has a relatively short iron-oxygen bond (1.67-1.76A) in agreement with a double-bond character, while the compound II state or the compound I state with a radical on an amino acid residue have a relatively long iron-oxygen bond (1.86-1.92A) in agreement with a single-bond character where the oxygen-atom is protonated.

Catalase↗

PubMed2006

The reaction mechanism of iron and manganese superoxide dismutases studied by theoretical calculations.

We have studied the detailed reaction mechanism of iron and manganese superoxide dismutase with density functional calculations on realistic active-site models, with large basis sets and including solvation, zero-point, and thermal effects. The results indicate that the conversion of O2- to O2 follows an associative mechanism, with O2- directly binding to the metal, followed by the protonation of the metal-bound hydroxide ion, and the dissociation of 3O2. All these reaction steps are exergonic. Likewise, we suggest that the conversion of O2- to H2O2 follows an at least a partly second-sphere pathway. There are small differences in the preferred oxidation and spin states, as well as in the geometries, of Fe and Mn, but these differences have little influence on the energetics, and therefore on the reaction mechanism of the two types of superoxide dismutases. For example, the two metals have very similar reduction potentials in the active-site models, although they differ by 0.7 V in water solution. The reaction mechanisms and spin states seem to have been designed to avoid spin conversions or to facilitate them by employing nearly degenerate spin states.

Binding Sites↗

PubMed2006

The structures of frataxin oligomers reveal the mechanism for the delivery and detoxification of iron.

Defects in the mitochondrial protein frataxin are responsible for Friedreich ataxia, a neurodegenerative and cardiac disease that affects 1:40,000 children. Here, we present the crystal structures of the iron-free and iron-loaded frataxin trimers, and a single-particle electron microscopy reconstruction of a 24 subunit oligomer. The structures reveal fundamental aspects of the frataxin mechanism. The trimer has a central channel in which one atom of iron binds. Two conformations of the channel with different metal-binding affinities suggest that a gating mechanism controls whether the bound iron is delivered to other proteins or transferred to detoxification sites. The trimer constitutes the basic structural unit of the 24 subunit oligomer. The architecture of this oligomer and several features of the trimer structure demonstrate striking similarities to the iron-storage protein ferritin. The data reveal how stepwise assembly provides frataxin with the structural flexibility to perform two functions: metal delivery and detoxification.

Amino Acid Sequence↗

PubMed2005

A combined quantum and molecular mechanical study of the O2 reductive cleavage in the catalytic cycle of multicopper oxidases.

The four-electron reduction of dioxygen to water in multicopper oxidases takes place in a trinuclear copper cluster, which is linked to a mononuclear blue copper site, where the substrates are oxidized. Recently, several intermediates in the catalytic cycle have been spectroscopically characterized, and two possible structural models have been suggested for both the peroxy and native intermediates. In this study, these spectroscopic results are complemented by hybrid quantum and molecular mechanical (QM/MM) calculations, taking advantage of recently available crystal structures with a full complement of copper ions. Thereby, we obtain optimized molecular structures for all of the experimentally studied intermediates involved in the reductive cleavage of the O(2) molecule and energy profiles for individual reaction steps. This allows identification of the experimentally observed intermediates and further insight into the reaction mechanism that is probably relevant for the whole class of multicopper oxidases. We suggest that the peroxy intermediate contains an O(2)(2-) ion, in which one oxygen atom bridges the type 2 copper ion and one of the type 3 copper ions, whereas the other one coordinates to the other type 3 copper ion. One-electron reduction of this intermediate triggers the cleavage of the O-O bond, which involves the uptake of a proton. The product of this cleavage is the observed native intermediate, which we suggest to contain a O(2-) ion coordinated to all three of the copper ions in the center of the cluster.

Algorithms↗

PubMed2005

How the Co-C bond is cleaved in coenzyme B12 enzymes: a theoretical study.

The homolytic cleavage of the organometallic Co-C bond in vitamin B12-dependent enzymes is accelerated by a factor of approximately 10(12) in the protein compared to that of the isolated cofactor in aqueous solution. To understand this much debated effect, we have studied the Co-C bond cleavage in the enzyme glutamate mutase with combined quantum and molecular mechanics methods. We show that the calculated bond dissociation energy (BDE) of the Co-C bond in adenosyl cobalamin is reduced by 135 kJ/mol in the enzyme. This catalytic effect can be divided into four terms. First, the adenosine radical is kept within 4.2 angstroms of the Co ion in the enzyme, which decreases the BDE by 20 kJ/mol. Second, the surrounding enzyme stabilizes the dissociated state by 42 kJ/mol using electrostatic and van der Waals interactions. Third, the protein itself is stabilized by 11 kJ/mol in the dissociated state. Finally, the coenzyme is geometrically distorted by the protein, and this distortion is 61 kJ/mol larger in the Co(III) state. This deformation of the coenzyme is caused mainly by steric interactions, and it is especially the ribose moiety and the Co-C5'-C4' angle that are distorted. Without the polar ribose group, the catalytic effect is much smaller, e.g. only 42 kJ/mol for methyl cobalamin. The deformation of the coenzyme is caused mainly by the substrate, a side chain of the coenzyme itself, and a few residues around the adenosine part of the coenzyme.

Carbon↗

PubMed2005

Quantum mechanical free energy barrier for an enzymatic reaction.

We discuss problems related to in silico studies of enzymes and show that accurate and converged free energy changes for complex chemical reactions can be computed if a method based on a thermodynamic cycle is employed. The method combines the sampling speed of molecular mechanics with the accuracy of a high-level quantum mechanics method. We use the method to compute the free energy barrier for a methyl transfer reaction catalyzed by the enzyme catechol O-methyltransferase at the level of density functional theory. The surrounding protein and solvent are found to have a profound effect on the reaction, and we show that energies can be extrapolated easily from one basis set and exchange-correlation functional to another. Using this procedure we calculate a barrier of 69 kJ/mol, in excellent agreement with the experimental value of 75 kJ/mol.

Catechol O-Methyltransferase↗

PubMed2005

Reaction mechanism of porphyrin metallation studied by theoretical methods.

We have studied the reaction mechanism for the insertion of Mg2+ and Fe2+ into a porphyrin ring with density functional calculations with large basis set and including solvation, zero-point and thermal effects. We have followed the reaction from the outer-sphere complex, in which the metal is coordinated with six water molecules and the porphyrin is doubly protonated, until the metal ion is inserted into the deprotonated porphyrin ring with only one water ligand remaining. This reaction involves the stepwise displacement of five water molecules and the removal of two protons from the porphyrin ring. In addition, a step seems to be necessary in which a porphyrin pyrrolenine nitrogen atom changes its interaction from a hydrogen bond to a metal-bound solvent molecule to a direct coordination to the metal ion. If the protons are taken up by a neutral imidazole molecule, the deprotonation reactions are exothermic with minimal barriers. However, with a water molecule as an acceptor, they are endothermic. The ligand exchange reactions were approximately thermoneutral (+/-20 kJ mol(-1), with one exception) with barriers of up to 72 kJ mol(-1) for Mg and 51 kJ mol(-1) for Fe. For Mg, the highest barrier was found for the formation of the first bond to the porphyrin ring. For Fe, a higher barrier was found for the formation of the second bond to the porphyrin ring, but this barrier is probably lower in solution. No evidence was found for an initial pre-equilibrium between a planar and a distorted porphyrin ring. Instead, the porphyrin becomes more and more distorted as the number of metal-porphyrin bonds increase (by up to 191 kJ mol(-1)). This strain is released when the porphyrin becomes deprotonated and the metal moves into the ring plane. Implications of these findings for the chelatase enzymes are discussed.

Iron↗

PubMed2005

NMR structure determination of proteins supplemented by quantum chemical calculations: detailed structure of the Ca2+ sites in the EGF34 fragment of protein S.

We present and test two methods to use quantum chemical calculations to improve standard protein structure refinement by molecular dynamics simulations restrained to experimental NMR data. In the first, we replace the molecular mechanics force field (employed in standard refinement to supplement experimental data) for a site of interest by quantum chemical calculations. This way, we obtain an accurate description of the site, even if a molecular-mechanics force field does not exist for this site, or if there is little experimental information about the site. Moreover, the site may change its bonding during the refinement, which often is the case for metal sites. The second method is to extract a molecular mechanics potential for the site of interest from a quantum chemical geometry optimisation and frequency calculation. We apply both methods to the two Ca2+ sites in the epidermal growth factor-like domains 3 and 4 in the vitamin K-dependent protein S and compare them to various methods to treat these sites in standard refinement. We show that both methods perform well and have their advantages and disadvantages. We also show that the glutamate Ca2+ ligand is unlikely to bind in a bidentate mode, in contrast to the crystal structure of an EGF domain of factor IX.

Binding Sites↗

PubMed2005

O2-binding to heme: electronic structure and spectrum of oxyheme, studied by multiconfigurational methods.

We have studied the ground state of a realistic model of oxyheme with multiconfigurational second-order perturbation theory (CASPT2). Our results show that the ground-state electronic structure is strongly multiconfigurational in character. Thus, the wavefunction is a mixture of many different configurations, of which the three most important ones are approximately 1FeII-1O2 (70%), FeIV-2O2(2-) (12%) and 3FeII-3O2 (3%). Thus, the wavefunction is dominated by closed-shell configurations, as suggested by Pauling, whereas the Weiss 2FeIII-2O2- configuration is not encountered among the 10 most important configurations. However, many other states are also important for this multiconfigurational wavefunction. Moreover, the traditional view is based on an oversimplified picture of the atomic-orbital contributions to the molecular orbitals. Thus, the population analysis indicates that all five iron orbitals are significantly occupied (by 0.5-2.0 electrons) and that the total occupation is most similar to the 3FeII-3O2 picture. The net charge on O2 is small, -0.20 e. Thus, it is quite meaningless to discuss which is the best valence-bond description of this inherently multiconfigurational system. Finally, we have calculated the eleven lowest ligand-field excited states of oxyheme and assigned the experimental spectrum of oxyhemoglobin with an average error of 0.24 eV.

Electrons↗