QM/MM study of the insertion of metal ion into protoporphyrin IX by ferrochelatase

Ferrochelatase catalyzes the metallation of protoporphyrin IX in the terminal step of heme biosynthesis. Mutations in the ferrochelatase gene can lead to the disease erythropoietic porphyria. The catalyzing mechanism of ferrochelatase is still not fully understood. In this paper, we have studied the insertion of Fe²⁺ into the protoporphyrin IX ring by Bacillussubtilis ferrochelatase using combined quantum mechanical and molecular mechanics (QM/MM) calculations. Geometries were optimized at the BP86/6-31G∗ level and energies were calculated at the B3LYP/TZVP level. The overall process involves the stepwise displacement of Glu-264, His-183, and a water molecule from Fe²⁺, and the removal of two protons from the porphyrin ring. The rate-determining step is the cleavage of the bond between the oxygen atom of Glu-264 and Fe²⁺, concomitant with the formation of the first Fe-N bond. It has an energy barrier of 57 kJ mol⁻¹. The porphyrin ring is only slightly distorted in the enzyme active site. The residue Tyr-13 plays a key role for the catalytic process extracting two protons from protoporphyrin IX.     

The carbohydrate-binding site in galectin-3 is preorganized to recognize a sugarlike framework of oxygens: ultra-high-resolution structures and water dynamics

The recognition of carbohydrates by proteins is a fundamental aspect of communication within and between living cells. Understanding the molecular basis of carbohydrate-protein interactions is a prerequisite for the rational design of synthetic ligands. Here we report the high- to ultra-high-resolution crystal structures of the carbohydrate recognition domain of galectin-3 (Gal3C) in the ligand-free state (1.08 ? at 100 K, 1.25 ? at 298 K) and in complex with lactose (0.86 ?) or glycerol (0.9 ?). These structures reveal striking similarities in the positions of water and carbohydrate oxygen atoms in all three states, indicating that the binding site of Gal3C is preorganized to coordinate oxygen atoms in an arrangement that is nearly optimal for the recognition of β-galactosides. Deuterium nuclear magnetic resonance (NMR) relaxation dispersion experiments and molecular dynamics simulations demonstrate that all water molecules in the lactose-binding site exchange with bulk water on a time scale of nanoseconds or shorter. Nevertheless, molecular dynamics simulations identify transient water binding at sites that agree well with those observed by crystallography, indicating that the energy landscape of the binding site is maintained in solution. All heavy atoms of glycerol are positioned like the corresponding atoms of lactose in the Gal3C complexes. However, binding of glycerol to Gal3C is insignificant in solution at room temperature, as monitored by NMR spectroscopy or isothermal titration calorimetry under conditions where lactose binding is readily detected. These observations make a case for protein cryo-crystallography as a valuable screening method in fragment-based drug discovery and further suggest that identification of water sites might inform inhibitor design.     

Exposure to Mitochondrial Genotoxins and Dopaminergic Neurodegeneration in Caenorhabditis elegans

Neurodegeneration has been correlated with mitochondrial DNA (mtDNA) damage and exposure to environmental toxins, but causation is unclear. We investigated the ability of several known environmental genotoxins and neurotoxins to cause mtDNA damage, mtDNA depletion, and neurodegeneration in Caenorhabditis elegans. We found that paraquat, cadmium chloride and aflatoxin B₁ caused more mitochondrial than nuclear DNA damage, and paraquat and aflatoxin B₁ also caused dopaminergic neurodegeneration. 6-hydroxydopamine (6-OHDA) caused similar levels of mitochondrial and nuclear DNA damage. To further test whether the neurodegeneration could be attributed to the observed mtDNA damage, C. elegans were exposed to repeated low-dose ultraviolet C radiation (UVC) that resulted in persistent mtDNA damage; this exposure also resulted in dopaminergic neurodegeneration. Damage to GABAergic neurons and pharyngeal muscle cells was not detected. We also found that fasting at the first larval stage was protective in dopaminergic neurons against 6-OHDA-induced neurodegeneration. Finally, we found that dopaminergic neurons in C. elegans are capable of regeneration after laser surgery. Our findings are consistent with a causal role for mitochondrial DNA damage in neurodegeneration, but also support non mtDNA-mediated mechanisms.     

The importance of porphyrin distortions for the ferrochelatase reaction

Ferrochelatase is the terminal enzyme in haem biosynthesis, i.e. the enzyme that inserts a ferrous ion into the porphyrin ring. Suggested reaction mechanisms for this enzyme involve a distortion of the porphyrin ring when it is bound to the enzyme. We have examined the energetics of such distortions using various theoretical calculations. With the density functional B3LYP method we calculate how much energy it costs to tilt one of the pyrrole rings out of the porphyrin plane for an isolated porphyrin molecule without or with a divalent metal ion in the centre of the ring. A tilt of 30 degrees costs 65-130 kJ/mol for most metal ions, but only approximately 48 kJ/mol for free-base (neutral) porphine. This indicates that once the metal is inserted, the porphyrin becomes stiffer and flatter, and therefore binds with lower affinity to a site designed to bind a distorted porphyrin. This would facilitate the release of the product from ferrochelatase. This proposal is strengthened by the fact that the only tested metal ion with a lower distortion energy than free-base porphyrin (Cd(2+)) is an inhibitor of ferrochelatase. Moreover, it costs even less energy to tilt a doubly deprotonated porphine(2-) molecule. This suggests that the protein may lower the acid constant of the pyrrole nitrogen atoms by deforming the porphyrin molecule. We have also estimated the structure of the protoporphyrin IX substrate bound to ferrochelatase using combined quantum chemical and molecular mechanics calculations. The result shows that the protein may distort the porphyrin by approximately 20 kJ/mol, leading to a distinctly non-planar structure. All four pyrrole rings are tilted out of the porphyrin mean plane (1-16 degrees ) but most towards the putative binding site of the metal ion. The predicted tilt is considerably smaller than that observed in the crystal structure of a porphyrin inhibitor.     

Metal binding to<Emphasis Type="Italic"> Bacillus subtilis</Emphasis> ferrochelatase and interaction between metal sites

Ferrochelatase, the terminal enzyme in heme biosynthesis, catalyses metal insertion into protoporphyrin IX. The location of the metal binding site with respect to the bound porphyrin substrate and the mode of metal binding are of central importance for understanding the mechanism of porphyrin metallation. In this work we demonstrate that Zn(2+), which is commonly used as substrate in assays of the ferrochelatase reaction, and Cd(2+), an inhibitor of the enzyme, bind to the invariant amino acids His183 and Glu264 and water molecules, all located within the porphyrin binding cleft. On the other hand, Mg(2+), which has been shown to bind close to the surface at 7 A from His183, was largely absent from its site. Activity measurements demonstrate that Mg(2+) has a stimulatory effect on the enzyme, lowering K(M) for Zn(2+) from 55 to 24 micro M. Changing one of the Mg(2+) binding residues, Glu272, to serine abolishes the effect of Mg(2+). It is proposed that prior to metal insertion the metal may form a sitting-atop (SAT) complex with the invariant His-Glu couple and the porphyrin. Metal binding to the Mg(2+) site may stimulate metal release from the protein ligands and its insertion into the porphyrin.     

Prediction and rationalization of the pH dependence of the activity and stability of family 11 xylanases

This paper presents a study of the pH dependence of the activity and stability of a set of family 11 xylanases for which X-ray structures are available, using the PROPKA approach. The xylanases are traditionally divided into basic and acidic xylanases, depending on whether the catalytic acid is hydrogen bonded to an Asn or Asp residue. Using X-ray structures, the predicted pH values of optimal activity of the basic xylanases are in the range of 5.2-6.9, which is in reasonable agreement with the available experimental values of 5-6.5. In the case of acidic xylanases, there are only four X-ray structures available, and using these structures, the predicted pHs of optimal activity are in the range of 4.2-5.0, compared to an observed range of 2-4.6. The influence of dynamical fluctuations of the protein structure is investigated for Bacillus agaradhaerens and Aspergillus kawachii xylanase using molecular dynamics (MD) simulations to provide snapshots from which average values can be computed. This decreases the respective predicted pH optima from 6.2-6.7 and 4.8 to 5.3 +/- 0.3 and 4.0 +/- 0.2, respectively, which are in better agreement with the observed values of 5.6 and 2, respectively. The change is primarily due to structural fluctuations of an Arg residue near the catalytic nucleophile, which lowers its pKa value compared to using the X-ray structure. The MD simulations and some X-ray structures indicate that this Arg residue can form a hydrogen bond to the catalytic base, and it is hypothesized that this hydrogen bond is stabilized by an additional hydrogen bond to another Glu residue present only in acidic xylanases. Formation of such a hydrogen bond is predicted to lower the pH optimum of A. kawachii xylanase to 2.9 +/- 0.3, which is in reasonable agreement with the observed value of 2. The predicted pH of optimal stability is in excellent agreement with the pH value at which the melting temperature (Tm) is greatest. Some correlation is observed between the pH-dependent free energy of unfolding and Tm, suggesting that the thermostability of the xylanases is partly due to a difference in residues with shifted pKa values. Thus, the thermostability of xylanases (and proteins in general) can perhaps be increased by mutations that introduce ionizable residues with pKa values significantly lower than standard values.     

Carboxylate binding modes in zinc proteins: A theoretical study

The relative energies of different coordination modes (bidentate, monodentate, syn, and anti) of a carboxylate group bound to a zinc ion have been studied by the density functional method B3LYP with large basis sets on realistic models of the active site of several zinc proteins. In positively charged four-coordinate complexes, the mono- and bidentate coordination modes have almost the same energy (within 10 kJ/mol). However, if there are negatively charged ligands other than the carboxylate group, the monodentate binding mode is favored. In general, the energy difference between monodentate and bidentate coordination is small, 4-24 kJ/mol, and it is determined more by hydrogen-bond interactions with other ligands or second-sphere groups than by the zinc-carboxylate interaction. Similarly, the activation energy for the conversion between the two coordination modes is small, approximately 6 kJ/mol, indicating a very flat Zn-O potential surface. The energy difference between syn and anti binding modes of the monodentate carboxylate group is larger, 70-100 kJ/mol, but this figure again strongly depends on interactions with second-sphere molecules. Our results also indicate that the pK(a) of the zinc-bound water ligand in carboxypeptidase and thermolysin is 8-9.     

Combined quantum and molecular mechanics calculations on metalloproteins

The combination of quantum mechanics and molecular mechanics (QM/MM) methods is one of the most promising approaches to study the structure, function and properties of proteins. The number of QM/MM applications on metalloproteins is steadily increasing, especially studies with density functional methods on redox-active metal centres. Recent developments include new parameterised methods to treat covalent bonds between the quantum and classical systems, methods to obtain free energy from QM/MM results, and the combination of quantum chemistry and protein crystallography.     

How are hydrogen bonds modified by metal binding?

We have used density functional theory calculations to investigate how the hydrogen-bond strength is modified when a ligand is bound to a metal using over 60 model systems involving six metals and eight ligands frequently encountered in metalloproteins. We study how the hydrogen-bond geometry and energy vary with the nature of metal, the oxidation state, the coordination number, the ligand involved in the hydrogen bond, other first-sphere ligands, and different hydrogen-bond probe molecules. The results show that, in general, the hydrogen-bond strength is increased for neutral ligands and decreased for negatively charged ligands. The size of the effect is mainly determined by the net charge of the metal complex, and all effects are typically decreased when the model is solvated. In water solution, the hydrogen-bond strength can increase by up to 37 kJ/mol for neutral ligands, and that of negatively charged ligands can increase (for complexes with a negative net charge) or decrease (for positively charged complexes). If the net charge of the complex does not change, there is normally little difference between different metals or different types of complexes. The only exception is observed for sulphur-containing ligands (Met and Cys) and if the ligand is redox-active (e.g. high-valence Fe–O complexes).     

Comparison of end-point continuum-solvation methods for the calculation of protein-ligand binding free energies

We have compared the predictions of ligand‐binding affinities from several methods based on end‐point molecular dynamics simulations and continuum solvation, that is, methods related to MM/PBSA (molecular mechanics combined with Poisson–Boltzmann and surface area solvation). Two continuum‐solvation models were considered, viz., the Poisson–Boltzmann (PB) and generalised Born (GB) approaches. The nonelectrostatic energies were also obtained in two different ways, viz., either from the sum of the bonded, van der Waals, nonpolar solvation energies, and entropy terms (as in MM/PBSA), or from the scaled protein–ligand van der Waals interaction energy (as in the linear interaction energy approach, LIE). Three different approaches to calculate electrostatic energies were tested, viz., the sum of electrostatic interaction energies and polar solvation energies, obtained either from a single simulation of the complex or from three independent simulations of the complex, the free protein, and the free ligand, or the linear‐response approximation (LRA). Moreover, we investigated the effect of scaling the electrostatic interactions by an effective internal dielectric constant of the protein (?_int). All these methods were tested on the binding of seven biotin analogues to avidin and nine 3‐amidinobenzyl‐1H‐indole‐2‐carboxamide inhibitors to factor Xa. For avidin, the best results were obtained with a combination of the LIE nonelectrostatic energies with the MM+GB electrostatic energies from a single simulation, using ?_int = 4. For fXa, standard MM/GBSA, based on one simulation and using ?_int = 4–10 gave the best result. The optimum internal dielectric constant seems to be slightly higher with PB than with GB solvation. © Proteins 2012; © 2012 Wiley Periodicals, Inc.