Coproheme decarboxylases (ChdCs) catalyze the final step in heme b biosynthesis of monoderm and some diderm bacteria. In this reaction, coproheme is converted to heme b via monovinyl monopropionate deuteroheme (MMD) in two consecutive decarboxylation steps. In Firmicutes decarboxylation of propionates 2 and 4 of coproheme depend on hydrogen peroxide and the presence of a catalytic tyrosine. Here we demonstrate that ChdCs from Actinobacteria are unique in using a histidine (H118 in ChdC from Corynebacterium diphtheriae, CdChdC) as a distal base in addition to the redox-active tyrosine (Y135). We present the X-ray crystal structures of coproheme-CdChdC and MMD-CdChdC, which clearly show (i) differences in the active site architecture between Firmicutes and Actinobacteria and (ii) rotation of the redox-active reaction intermediate (MMD) after formation of the vinyl group at position 2. Distal H118 is shown to catalyze the heterolytic cleavage of hydrogen peroxide (k app = (4.90 ± 1.25) × 104 M–1 s–1). The resulting Compound I is rapidly converted to a catalytically active Compound I* (oxoiron­(IV) Y135•) that initiates the radical decarboxylation reactions. As a consequence of the more efficient Compound I formation, actinobacterial ChdCs exhibit a higher catalytic efficiency in comparison to representatives from Firmicutes. On the basis of the kinetic data of wild-type CdChdC and the variants H118A, Y135A, and H118A/Y135A together with high-resolution crystal structures and molecular dynamics simulations, we present a molecular mechanism for the hydrogen peroxide dependent conversion of coproheme via MMD to heme b and discuss differences between ChdCs from Actinobacteria and Firmicutes.

Alchemically derived free energies are artifacted when the perturbed moiety has a nonzero net charge. The source of the artifacts lies in the effective treatment of the electrostatic interactions within and between the perturbed atoms and remaining (partial) charges in the simulated system. To treat the electrostatic interactions effectively, lattice‐summation (LS) methods or cutoff schemes in combination with a reaction‐field contribution are usually employed. Both methods render the charging component of the calculated free energies sensitive to essential parameters of the system like the cutoff radius or the box side lengths. Here, we discuss the results of three previously published studies of ligand binding. These studies presented estimates of binding free energies that were artifacted due to the charged nature of the ligands. We show that the size of the artifacts can be efficiently calculated and raw simulation data can be corrected. We compare the corrected results with experimental estimates and nonartifacted estimates from path‐sampling methods. Although the employed correction scheme involves computationally demanding continuum‐electrostatics calculations, we show that the correction estimate can be deduced from a small sample of configurations rather than from the entire ensemble. This observation makes the calculations of correction terms feasible for complex biological systems. To show the general applicability of the proposed procedure, we also present results where the correction scheme was used to correct independent free energies obtained from simulations employing a cutoff scheme or LS electrostatics. In this work, we give practical guidelines on how to apply the appropriate corrections easily.

Currently, two different methods dominate the field of biomolecular free-energy calculations for the prediction of binding affinities. Pathway methods are frequently used for large ligands that bind on the surface of a host, such as protein–protein complexes. Alchemical methods, on the other hand, are preferably applied for small ligands that bind to deeply buried binding sites. The latter methods are also widely known to be heavily artifacted by the representation of electrostatic energies in periodic simulation boxes, in particular, when net-charge changes are involved. Different methods have been described to deal with these artifacts, including postsimulation correction schemes and instantaneous correction schemes (e.g., co-alchemical perturbation of ions). Here, we use very simple test systems to show that instantaneous correction schemes with no change in the system net charge lower the artifacts but do not eliminate them. Furthermore, we show that free energies from pathway methods suffer from the same artifacts.

Monoderm bacteria utilize coproheme decarboxylases (ChdCs) to generate heme b by a stepwise decarboxylation of two propionate groups of iron coproporphyrin III (coproheme), forming two vinyl groups. This work focuses on actinobacterial ChdC from Corynebacterium diphtheriae (CdChdC) to elucidate the hydrogen peroxide-mediated decarboxylation of coproheme via monovinyl monopropionyl deuteroheme (MMD) to heme b, with the principal aim being to understand the reorientation mechanism of MMD during turnover. Wild-type CdChdC and variants, namely H118A, H118F, and A207E, were studied by resonance Raman and ultraviolet-visible spectroscopy, mass spectrometry, and molecular dynamics simulations. As actinobacterial ChdCs use a histidine (H118) as a distal base, we studied the H118A and H118F variants to elucidate the effect of 1) the elimination of the proton acceptor and 2) steric constraints within the active site. The A207E variant mimics the proximal H-bonding network found in chlorite dismutases. This mutation potentially increases the rigidity of the proximal site and might impair the rotation of the reaction intermediate MMD. We found that both wild-type CdChdC and the variant H118A convert coproheme mainly to heme b upon titration with H 2 O 2 . Interestingly, the variant A207E mostly accumulates MMD along with small amounts of heme b, whereas H118F is unable to produce heme b and accumulates only MMD. Together with molecular dynamics simulations, the spectroscopic data provide insight into the reaction mechanism and the mode of reorientation of MMD, i.e., a rotation in the active site versus a release and rebinding.

Hybrid quantum mechanics/molecular mechanics (QM/MM) simulations have advanced the field of computational chemistry tremendously. However, they require the partitioning of a system into two different regions that are treated at different levels of theory, which can cause artifacts at the interface. Furthermore, they are still limited by high computational costs of quantum chemical calculations. In this work, we develop the buffer region neural network (BuRNN), an alternative approach to existing QM/MM schemes, which introduces a buffer region that experiences full electronic polarization by the inner QM region to minimize artifacts. The interactions between the QM and the buffer region are described by deep neural networks (NNs), which leads to the high computational efficiency of this hybrid NN/MM scheme while retaining quantum chemical accuracy. We demonstrate the BuRNN approach by performing NN/MM simulations of the hexa-aqua iron complex.