Proceedings of the National Academy of Sciences of the United States of America, Jan 31, 2015
Most of the main features of water oxidation in photosystem II are now well understood, including... more Most of the main features of water oxidation in photosystem II are now well understood, including the mechanism for O-O bond formation. For the intermediate S2 and S3 structures there is also nearly complete agreement between quantum chemical modeling and experiments. Given the present high degree of consensus for these structures, it is of high interest to go back to previous suggestions concerning what happens in the S2-S3 transition. Analyses of extended X-ray adsorption fine structure (EXAFS) experiments have indicated relatively large structural changes in this transition, with changes of distances sometimes larger than 0.3 Å and a change of topology. In contrast, our previous density functional theory (DFT)(B3LYP) calculations on a cluster model showed very small changes, less than 0.1 Å. It is here found that the DFT structures are also consistent with the EXAFS spectra for the S2 and S3 states within normal errors of DFT. The analysis suggests that there are severe problems ...
The heme peroxidases have a histidine group as the axial ligand of iron. This ligand forms a hydr... more The heme peroxidases have a histidine group as the axial ligand of iron. This ligand forms a hydrogen bond to an aspartate carboxylate group by the other nitrogen atom in the side chain. The aspartate is not present in the globins and it has been suggested that it gives an imidazolate character to the histidine ligand. Quantum chemical calculations have indicated that the properties of the heme site strongly depend on the position of the proton in this hydrogen bond. Therefore, we have studied the location of this proton in all intermediates in the reaction mechanism, using a set of different quantum mechanical and combined experimental and computational methods. Quantum refinements of a crystal structure of the resting FeIII state in yeast cytochrome c peroxidase show that the geometric differences of the two states are so small that it cannot be unambiguously decided where the proton is in the crystal structure. Vacuum calculations indicate that the position of the proton is sensitive to the surroundings and to the side chains of the porphyrin ring. Combined quantum and molecular mechanics (QM/MM) calculations indicate that the proton prefers to reside on the His ligand in all states in the reaction mechanism of the peroxidases. QM/MM free energy perturbations confirm these results, but reduce the energy difference between the two states to 12-44 kJ/mol.
We have developed force-field parameters for the hydrogen-abstraction transition state of aliphat... more We have developed force-field parameters for the hydrogen-abstraction transition state of aliphatic hydroxylation by cytochrome P450 using the Q2MM approach. The parametrization is based on quantum chemical (B3LYP) transition-state structures and Hessian matrices for 24 diverse substrate models (14 in the training set and 10 in the test set). The force field is intended to be applicable to any druglike molecule by the use of the general Amber force field (GAFF) for the substrates. The parameters reproduce the geometries within 0.1 Å and 1.2° for bond lengths and angles, respectively, with no significant differences between the training and test sets. The Hessian matrix is also well reproduced with a correlation coefficient of 0.99. The parametrization is performed by the ideal iterative approach of Norrby and Liljefors, which we have implemented for the Amber software.
The oxidation and dealkylation of dimethylsulfoxide (DMSO), dimethylsulfide (DMS), and trimethyla... more The oxidation and dealkylation of dimethylsulfoxide (DMSO), dimethylsulfide (DMS), and trimethylamine (TMA) by cytochrome P450 has been studied with density functional theory calculations. The results show that the oxidation reactions always occur on the doublet spin surface, whereas dealkylations can take place for both the doublet and quartet spin states. Moreover, DMS is more reactive than DMSO, and S-oxidation is more favorable than S-dealkylation, whereas N-dealkylation is more favorable than N-oxidation. This is in perfect agreement with experimental results, showing that density functional activation energies are reliable and comparable for widely different reactions with cytochrome P450.
We have estimated the activation energy for hydrogen abstraction by compound I in cytochrome P450... more 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 de...
We have studied the dynamics of water molecules in six crystal structures of four human cytochrom... more We have studied the dynamics of water molecules in six crystal structures of four human cytochromes P450, 2A6, 2C8, 2C9, and 3A4, with molecular dynamics simulations. In the crystal structures, only a few water molecules are seen and the reported sizes of the active-site cavity vary a lot. In the simulations, the cavities are completely filled with water molecules, although with approximately 20% lower density than in bulk water. The 2A6 protein differs from the other three in that it has a very small cavity with only two water molecules and no exchange with the surroundings. The other three proteins have quite big cavities, with 41 water molecules on average in 2C8 and 54-58 in 2C9 and 3A4, giving a water volume of 1500-2100 A3. The two crystal structures of 2C9 differ quite appreciably, whereas those of 3A4 are quite similar. The active-site cavity is connected to the surroundings by three to six channels, through which there is a quite frequent exchange of water molecules (one molecule is exchanged every 30-200 ps), except in 2A6. Most of the channels are observed also in the crystal structures, but two to three channels in each protein open only during the simulations. There are no water molecules close to the heme iron ion in these simulations of the high-spin ferric state (the average distance to the closest water molecule is 3.3-5 A), and there are few ordered water molecules in the active sites, none of which is conserved in all proteins.
We have estimated the activation energy for aromatic oxidation by compound I in cytochrome P450 f... more We have estimated the activation energy for aromatic oxidation by compound I in cytochrome P450 for a diverse set of 17 substrates using state-of-the-art density functional theory (B3LYP) with large basis sets. The activation energies vary from 60 to 87 kJ/mol. We then test if these results can be reproduced by computationally less demanding methods. The best methods (a B3LYP calculation of the activation energy of a methoxy-radical model or a partial least-squares model of the semiempirical AM1 bond dissociation energies and spin densities of the tetrahedral intermediate for both a hydroxyl-cation and a hydroxyl-radical model) give correlations with r(2) of 0.8 and mean absolute deviations of 3 kJ/mol. Finally, we apply these simpler methods on several sets of reactions for which experimental data are available and show that we can predict the reactive sites by combining calculations of the activation energies with the solvent-accessible surface area of each site.
We have used the density functional B3LYP method to study the effect of hydrogen bonds from the h... more We have used the density functional B3LYP method to study the effect of hydrogen bonds from the histidine ligand in various haem proteins to carboxyl groups or to the carbonyl backbone. Hydrogen bonds to carbonyl groups (encountered in globins and cytochromes, for example) have a small influence on the geometry and properties of the haem site. However, hydrogen bonds to
The homolytic cleavage of the organometallic Co-C bond in vitamin B12-dependent enzymes is accele... more 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.
We have studied the ground state of a realistic model of oxyheme with multiconfigurational second... more 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.
Proceedings of the National Academy of Sciences of the United States of America, Jan 31, 2015
Most of the main features of water oxidation in photosystem II are now well understood, including... more Most of the main features of water oxidation in photosystem II are now well understood, including the mechanism for O-O bond formation. For the intermediate S2 and S3 structures there is also nearly complete agreement between quantum chemical modeling and experiments. Given the present high degree of consensus for these structures, it is of high interest to go back to previous suggestions concerning what happens in the S2-S3 transition. Analyses of extended X-ray adsorption fine structure (EXAFS) experiments have indicated relatively large structural changes in this transition, with changes of distances sometimes larger than 0.3 Å and a change of topology. In contrast, our previous density functional theory (DFT)(B3LYP) calculations on a cluster model showed very small changes, less than 0.1 Å. It is here found that the DFT structures are also consistent with the EXAFS spectra for the S2 and S3 states within normal errors of DFT. The analysis suggests that there are severe problems ...
The heme peroxidases have a histidine group as the axial ligand of iron. This ligand forms a hydr... more The heme peroxidases have a histidine group as the axial ligand of iron. This ligand forms a hydrogen bond to an aspartate carboxylate group by the other nitrogen atom in the side chain. The aspartate is not present in the globins and it has been suggested that it gives an imidazolate character to the histidine ligand. Quantum chemical calculations have indicated that the properties of the heme site strongly depend on the position of the proton in this hydrogen bond. Therefore, we have studied the location of this proton in all intermediates in the reaction mechanism, using a set of different quantum mechanical and combined experimental and computational methods. Quantum refinements of a crystal structure of the resting FeIII state in yeast cytochrome c peroxidase show that the geometric differences of the two states are so small that it cannot be unambiguously decided where the proton is in the crystal structure. Vacuum calculations indicate that the position of the proton is sensitive to the surroundings and to the side chains of the porphyrin ring. Combined quantum and molecular mechanics (QM/MM) calculations indicate that the proton prefers to reside on the His ligand in all states in the reaction mechanism of the peroxidases. QM/MM free energy perturbations confirm these results, but reduce the energy difference between the two states to 12-44 kJ/mol.
We have developed force-field parameters for the hydrogen-abstraction transition state of aliphat... more We have developed force-field parameters for the hydrogen-abstraction transition state of aliphatic hydroxylation by cytochrome P450 using the Q2MM approach. The parametrization is based on quantum chemical (B3LYP) transition-state structures and Hessian matrices for 24 diverse substrate models (14 in the training set and 10 in the test set). The force field is intended to be applicable to any druglike molecule by the use of the general Amber force field (GAFF) for the substrates. The parameters reproduce the geometries within 0.1 Å and 1.2° for bond lengths and angles, respectively, with no significant differences between the training and test sets. The Hessian matrix is also well reproduced with a correlation coefficient of 0.99. The parametrization is performed by the ideal iterative approach of Norrby and Liljefors, which we have implemented for the Amber software.
The oxidation and dealkylation of dimethylsulfoxide (DMSO), dimethylsulfide (DMS), and trimethyla... more The oxidation and dealkylation of dimethylsulfoxide (DMSO), dimethylsulfide (DMS), and trimethylamine (TMA) by cytochrome P450 has been studied with density functional theory calculations. The results show that the oxidation reactions always occur on the doublet spin surface, whereas dealkylations can take place for both the doublet and quartet spin states. Moreover, DMS is more reactive than DMSO, and S-oxidation is more favorable than S-dealkylation, whereas N-dealkylation is more favorable than N-oxidation. This is in perfect agreement with experimental results, showing that density functional activation energies are reliable and comparable for widely different reactions with cytochrome P450.
We have estimated the activation energy for hydrogen abstraction by compound I in cytochrome P450... more 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 de...
We have studied the dynamics of water molecules in six crystal structures of four human cytochrom... more We have studied the dynamics of water molecules in six crystal structures of four human cytochromes P450, 2A6, 2C8, 2C9, and 3A4, with molecular dynamics simulations. In the crystal structures, only a few water molecules are seen and the reported sizes of the active-site cavity vary a lot. In the simulations, the cavities are completely filled with water molecules, although with approximately 20% lower density than in bulk water. The 2A6 protein differs from the other three in that it has a very small cavity with only two water molecules and no exchange with the surroundings. The other three proteins have quite big cavities, with 41 water molecules on average in 2C8 and 54-58 in 2C9 and 3A4, giving a water volume of 1500-2100 A3. The two crystal structures of 2C9 differ quite appreciably, whereas those of 3A4 are quite similar. The active-site cavity is connected to the surroundings by three to six channels, through which there is a quite frequent exchange of water molecules (one molecule is exchanged every 30-200 ps), except in 2A6. Most of the channels are observed also in the crystal structures, but two to three channels in each protein open only during the simulations. There are no water molecules close to the heme iron ion in these simulations of the high-spin ferric state (the average distance to the closest water molecule is 3.3-5 A), and there are few ordered water molecules in the active sites, none of which is conserved in all proteins.
We have estimated the activation energy for aromatic oxidation by compound I in cytochrome P450 f... more We have estimated the activation energy for aromatic oxidation by compound I in cytochrome P450 for a diverse set of 17 substrates using state-of-the-art density functional theory (B3LYP) with large basis sets. The activation energies vary from 60 to 87 kJ/mol. We then test if these results can be reproduced by computationally less demanding methods. The best methods (a B3LYP calculation of the activation energy of a methoxy-radical model or a partial least-squares model of the semiempirical AM1 bond dissociation energies and spin densities of the tetrahedral intermediate for both a hydroxyl-cation and a hydroxyl-radical model) give correlations with r(2) of 0.8 and mean absolute deviations of 3 kJ/mol. Finally, we apply these simpler methods on several sets of reactions for which experimental data are available and show that we can predict the reactive sites by combining calculations of the activation energies with the solvent-accessible surface area of each site.
We have used the density functional B3LYP method to study the effect of hydrogen bonds from the h... more We have used the density functional B3LYP method to study the effect of hydrogen bonds from the histidine ligand in various haem proteins to carboxyl groups or to the carbonyl backbone. Hydrogen bonds to carbonyl groups (encountered in globins and cytochromes, for example) have a small influence on the geometry and properties of the haem site. However, hydrogen bonds to
The homolytic cleavage of the organometallic Co-C bond in vitamin B12-dependent enzymes is accele... more 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.
We have studied the ground state of a realistic model of oxyheme with multiconfigurational second... more 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.
Uploads
Papers by Ulf Ryde