QM/MM study of the binding of H2 to MoCu CO dehydrogenase: development and applications of improved H2 van der Waals parameters

Understanding the Directed Evolution of a Natural-like Efficient Artificial Metalloenzyme

The artificial metalloenzyme containing iridium in place of iron along with four directed evolution mutations C317G, T213G, L69V, and V254L in a natural cytochrome P450 presents an important milestone in merging the extraordinary efficiency of biocatalysts with the versatility of small molecule chemical catalysts in catalyzing a new-to-nature carbene insertion reaction. This is a show-stopper enzyme, as it exhibits a catalytic efficiency similar to that of natural enzymes. Despite this remarkable discovery, there is no mechanistic and structural understanding as to why it displays extraordinary efficiency after the incorporation of the four active site mutations by directed evolution methods, which so far has been intractable to any experimental methods. In this study, we have deciphered how directed evolution mutations gradually alter the protein conformational ensemble to populate a catalytically active conformation to boost a multistep catalysis in a natural-like artificial metalloenzyme using large-scale molecular dynamics simulations, rigorous quantum chemical (QM), and multiscale quantum chemical/molecular mechanics (QM/MM) calculations. It reveals how evolution precisely positions the cofactor-substrate in an unusual but effective orientation within a reshaped active site in the catalytically active conformation stabilized by C-H···π interactions from more ordered mutated L69V and V254L residues to achieve preferential transition state stabilization compared to the ground state. This work essentially tracks down in atomistic detail the shift in the conformational ensemble of the highly active conformation from the less efficient single mutant to the most efficient quadruple mutant and offers valuable insights for designing better enzymes. The active conformation correctly reproduces the experimental barrier height and also accounts for the catalytic effect, which is in good agreement with experimental observations. Moreover, this conformation features an unusual bonding interaction in a metal-carbene species that preferentially stabilizes the rate-determining formation of an iridium porphyrin carbene intermediate to render the observed high catalytic rate acceleration. Our study provides crucial insights into the underlying rationale for directed evolution, reports the major catalytic role of nonelectrostatic interactions in enzyme catalysis different from the electrostatic model, and suggests a crucial principle toward designing enzymes with natural efficiency.

CO Oxidation Mechanism of Silver-Substituted Mo/Cu CO-Dehydrogenase - Analogies and Differences to the Native Enzyme

The aerobic oxidation of carbon monoxide to carbon dioxide is catalysed by the Mo/Cu-containing CO-dehydrogenase enzyme in the soil bacterium Oligotropha carboxidovorans, enabling the organism to grow on the small gas molecule as carbon and energy source. It was shown experimentally that silver can be substituted for copper in the active site of Mo/Cu CODH to yield a functional enzyme. In this study, we employed QM/MM calculations to investigate whether the reaction mechanism of the silver-substituted enzyme is similar to that of the native enzyme. Our results suggest that the Ag-substituted enzyme can oxidize CO and release CO2 following the same reaction steps as the native enzyme, with a computed rate-limiting step of 10.4 kcal/mol, consistent with experimental findings. Surprisingly, lower activation energies for C-O bond formation have been found in the presence of silver. Furthermore, comparison of rate constants for reduction of copper- and silver-containing enzymes suggests a discrepancy in the transition state stabilization upon silver substitution. We also evaluated the effects that differences in the water-active site interaction may exert on the overall energy profile of catalysis. Finally, the formation of a thiocarbonate intermediate along the catalytic pathway was found to be energetically unfavorable for the Ag-substituted enzyme. This finding aligns with the hypothesis proposed for the wild-type form, suggesting that the creation of such species may not be necessary for the enzymatic catalysis of CO oxidation.

How general is the effect of the bulkiness of organic ligands on the basicity of metal-organic catalysts? H2-evolving Mo oxides/sulphides as case studies

Tailoring the activity of an organometallic catalyst usually requires a targeted ligand design. Tuning the ligand bulkiness and tuning the electronic properties are popular approaches, which are somehow interdependent because substituents of different sizes within ligands can determine inter alia the occurrence of different degrees of inductive effects. Ligand basicity, in particular, turned out to be a key property for the modulation of protonation reactions occurring in vacuo at the metals in complexes bearing organophosphorus ligands; however, when the same reactions take place in a polar organic solvent, their energetics becomes dependent on the trade-off between ligand basicity and bulkiness, with the polarity of the solvent playing a key role in this regard [Bancroft et al., Inorg. Chem., 1986, 25, 3675; Rovaletti et al., J. Phys. Org. Chem., 2018, 31, e3748]. In the present contribution, we carried out molecular dynamics and density functional theory calculations on water-soluble Mo-based catalysts for proton reduction, in order to study the energetics of protonation reactions in complexes where the incipient proton binds a catalytically active ligand (i.e., an oxide or a disulphide). We considered complexes either soaked in water or in a vacuum, and featuring N-based ancillary ligands of different bulkiness (i.e. cages constituted either by pyridine or isoquinoline moieties). Our results show that the energetics of protonation events can be affected by ancillary ligand bulkiness even when the metal center does not play the role of the H⁺ acceptor. In vacuo, protonation at the O or S atom in the α position relative to the metal in complexes featuring the bulky isoquinoline-based ligand is more favored by around 10 kcal mol^-1 when compared to the case of the pyridine-based counterparts, a difference that is almost zero when the same reactions occur in water. Such an outcome is rationalized in light of the different electrostatic properties of complexes bearing ancillary ligands of different sizes. The overall picture from theory indicates that such effects of ligand bulkiness can be relevant for the design of green chemistry catalysts that undergo protonation steps in water solutions.

DFT study on the gold(I)-catalyzed cycloaddition and rearrangement reactions of allene-containing allylic silyl ether

The DFT calculation at the B3LYP/B3LYP-D3(BJ) level was carried out to explore the reaction mechanism of the synthesis of spirocyclo[4,5]decane skeleton by gold-catalyzed allenyl compounds. The more accurate energy under the CH₃CN solvent in the experiment is calculated by the single-point energy of the SMD model. Computational studies have shown that the reaction consists of three main steps: intramolecular cycloaddition of the end group carbon atoms of allenyl and vinyl groups, the semipinacol rearrangement process in which the four-membered ring is reconstructed into the five-membered ring, the elimination reaction releases the catalyst and obtains the product. The calculation results show that Zheng et al. reported that the gold-catalyzed synthesis reaction can easily occur under the experimental conditions due to its low activation free energy (12.07-15.49 kcal/mol). Furthermore, it was found that the MOMO(CH₂)₂ substituent has higher reactivity than the corresponding reactant of the phenyl substituent.