Maréchal, J.-D.

ACS Catal. 2014, 4, 833-842, 10.1021/cs400921n

We present a computational study that combines protein–ligand docking, quantum mechanical, and quantum mechanical/molecular mechanical calculations to scrutinize the mechanistic behavior of the first artificial enzyme able to enantioselectively reduce cyclic imines. We applied a novel strategy that allows the characterization of transition state structures in the protein host and their associated reaction paths. Of the most striking results of our investigation is the identification of major conformational differences between the transition state geometries of the lowest energy paths leading to (R)- and (S)-reduction products. The molecular features of (R)- and (S)-transition states highlight distinctive patterns of hydrophobic and polar complementarities between the substrate and the binding site. These differences lead to an activation energy gap that stands in very good agreement with the experimentally determined enantioselectivity. This study sheds light on the mechanism by which transfer hydrogenases operate and illustrates how the change of environment (from homogeneous solution conditions to the asymmetric protein frame) affect the reactivity of the organometallic cofactor. It provides novel insights on the complexity in integrating unnatural organometallic compounds into biological scaffolds. The modeling strategy that we pursued, based on the generation of “pseudo transition state” structures, is computationally efficient and suitable for the discovery and optimization of artificial enzymes. Alternatively, this approach can be applied on systems for which a large conformational sampling is needed to identify relevant transition states.

Maréchal, J.-D.; Roelfes, G.

Chem. Sci. 2017, 8, 7228-7235, 10.1039/C7SC03477F

The design of artificial metalloenzymes is a challenging, yet ultimately highly rewarding objective because of the potential for accessing new-to-nature reactions. One of the main challenges is identifying catalytically active substrate–metal cofactor–host geometries. The advent of expanded genetic code methods for the in vivo incorporation of non-canonical metal-binding amino acids into proteins allow to address an important aspect of this challenge: the creation of a stable, well-defined metal-binding site. Here, we report a designed artificial metallohydratase, based on the transcriptional repressor lactococcal multidrug resistance regulator (LmrR), in which the non-canonical amino acid (2,2′-bipyridin-5yl)alanine is used to bind the catalytic Cu(II) ion. Starting from a set of empirical pre-conditions, a combination of cluster model calculations (QM), protein–ligand docking and molecular dynamics simulations was used to propose metallohydratase variants, that were experimentally verified. The agreement observed between the computationally predicted and experimentally observed catalysis results demonstrates the power of the artificial metalloenzyme design approach presented here.

Maréchal, J.-D.; Ward, T.R.

Angew. Chem. Int. Ed. 2018, 57, 1863-1868, 10.1002/anie.201711016

Artificial metalloenzymes, resulting from incorporation of a metal cofactor within a host protein, have received increasing attention in the last decade. The directed evolution is presented of an artificial transfer hydrogenase (ATHase) based on the biotin‐streptavidin technology using a straightforward procedure allowing screening in cell‐free extracts. Two streptavidin isoforms were yielded with improved catalytic activity and selectivity for the reduction of cyclic imines. The evolved ATHases were stable under biphasic catalytic conditions. The X‐ray structure analysis reveals that introducing bulky residues within the active site results in flexibility changes of the cofactor, thus increasing exposure of the metal to the protein surface and leading to a reversal of enantioselectivity. This hypothesis was confirmed by a multiscale approach based mostly on molecular dynamics and protein–ligand dockings.

Maréchal, J.-D.

Acc. Chem. Res. 2020, 53, 896-905, 10.1021/acs.accounts.0c00031

Artificial metalloenzymes (ArMs) are obtained by inserting homogeneous catalysts into biological scaffolds and are among the most promising strategies in the quest for new-to-nature biocatalysts. The quality of their design strongly depends on how three partners interact: the biological host, the “artificial cofactor,” and the substrate. However, structural characterization of functional artificial metalloenzymes by X-ray or NMR is often partial, elusive, or absent. How the cofactor binds to the protein, how the receptor reorganizes upon the binding of the cofactor and the substrate, and which are the binding mode(s) of the substrate for the reaction to proceed are key questions that are frequently unresolved yet crucial for ArM design. Such questions may eventually be solved by molecular modeling but require a step change beyond the current state-of-the-art methodologies. Here, we summarize our efforts in the study of ArMs, presenting both the development of computational strategies and their application. We first focus on our integrative computational framework that incorporates a variety of methods such as protein–ligand docking, classical molecular dynamics (MD), and pure quantum mechanical (QM) methods, which, when properly combined, are able to depict questions that range from host–cofactor binding predictions to simulations of entire catalytic mechanisms. We also pay particular attention to the protein–ligand docking strategies that we have developed to accurately predict the binding of transition metal-containing molecules to proteins. While this aspect is fundamental to many bioinorganic fields beyond ArMs, it has been disregarded from the molecular modeling landscape until very recently. Next we describe how to apply this computational framework to particular ArMs including systems previously characterized experimentally as well as others where computation served to guide the design. We start with the prediction of the interactions between homogeneous catalysts and biological hosts. Protein–ligand docking is pivotal at that stage, but it needs to be combined with QM/MM or MD approaches when the binding of the cofactor implies significant conformational changes of the protein or involve changes of the electronic state of the metal. Then, we summarize molecular modeling studies aimed at identifying cofactor–substrate arrangements inside the ArM active pocket that are consistent with its reactivity. These calculations stand on “Theozyme”-like dockings, MD-refined or not, which provide molecular rationale of the catalytic profiles of the artificial systems. In the third section, we present case studies to decode the entire catalytic mechanism of two ArMs: (1) an iridium based asymmetric transfer hydrogenase obtained by insertion of Noyori’s catalyst into streptavidin and (2) a metallohydrolase achieved by including a receptor. Transition states, second coordination sphere effects, as well as motions of the cofactors are identified as drivers of the enantiomeric profiles. Finally, we report computer-aided designs of ArMs to guide experiments toward chemical and mutational changes that improve their activity and/or enantioselective profiles and expand toward future directions.

Matsuo, T.

Chem. - Asian J. 2011, 6, 2491-2499, 10.1002/asia.201100107

H64D myoglobin mutant was reconstituted with two different types of synthetic hemes that have aromatic rings and a carboxylate‐based cluster attached to the terminus of one or both of the heme‐propionate moieties, thereby forming a “single‐winged cofactor” and “double‐winged cofactor,” respectively. The reconstituted mutant myoglobins have smaller Km values with respect to 2‐methoxyphenol oxidation activity relative to the parent mutant with native heme. This suggests that the attached moiety functions as a substrate‐binding domain. However, the kcat value of the mutant myoglobin with the double‐winged cofactor is much lower than that of the mutant with the native heme. In contrast, the mutant reconstituted with the single‐winged cofactor has a larger kcat value, thereby resulting in overall catalytic activity that is essentially equivalent to that of the native horseradish peroxidase. Enhanced peroxygenase activity was also observed for the mutant myoglobin with the single‐winged cofactor, thus indicating that introduction of an artificial substrate‐binding domain at only one of the heme propionates in the H64D mutant is the optimal engineering strategy for improving the peroxidase activity of myoglobin.

Matsuo, T.

Tetrahedron Lett. 2019, 60, 151226, 10.1016/j.tetlet.2019.151226

Organic synthesis using biocatalysts has been developed over many years and is still a prominent area of research. In this context, various hybrid biocatalysts composed of a synthetic metal complex catalyst and a protein scaffold (i.e. artificial metalloenzymes) have been constructed. One of the most recent research areas in biocatalysts-mediated synthesis is CC bond/cleavage, the most important type of reaction in organic chemistry. Some of the artificial enzymes were applied to in-cell reactions as well as in vitro systems. The effects of the structural fluctuation in biomacromolecules on their functions have also been realized. This review article includes recent research examples of artificial metalloenzymes used to CC bond formation/cleavage. As a perspective, we also focus on how we apply protein dynamics factor for the creation of new generation artificial metalloenzymes.

Maréchal, J.-D.; Ward, T.R.

J. Am. Chem. Soc. 2014, 136, 15676-15683, 10.1021/ja508258t

An artificial imine reductase results upon incorporation of a biotinylated Cp*Ir moiety (Cp* = C5Me5–) within homotetrameric streptavidin (Sav) (referred to as Cp*Ir(Biot-p-L)Cl] ⊂ Sav). Mutation of S112 reveals a marked effect of the Ir/streptavidin ratio on both the saturation kinetics as well as the enantioselectivity for the production of salsolidine. For [Cp*Ir(Biot-p-L)Cl] ⊂ S112A Sav, both the reaction rate and the selectivity (up to 96% ee (R)-salsolidine, kcat 14–4 min–1 vs [Ir], KM 65–370 mM) decrease upon fully saturating all biotin binding sites (the ee varying between 96% ee and 45% ee R). In contrast, for [Cp*Ir(Biot-p-L)Cl] ⊂ S112K Sav, both the rate and the selectivity remain nearly constant upon varying the Ir/streptavidin ratio [up to 78% ee (S)-salsolidine, kcat 2.6 min–1, KM 95 mM]. X-ray analysis complemented with docking studies highlight a marked preference of the S112A and S112K Sav mutants for the SIr and RIr enantiomeric forms of the cofactor, respectively. Combining both docking and saturation kinetic studies led to the formulation of an enantioselection mechanism relying on an “induced lock-and-key” hypothesis: the host protein dictates the configuration of the biotinylated Ir-cofactor which, in turn, by and large determines the enantioselectivity of the imine reductase.

Maréchal, J.-D.

ACS Catal. 2015, 5, 2469-2480, 10.1021/acscatal.5b00010

The development of artificial enzymes aims at expanding the scope of biocatalysis. Over recent years, artificial metalloenzymes based on the insertion of homogeneous catalysts in biomolecules have received an increasing amount of attention. Rational or pseudorational design of these composites is a challenging task because of the complexity of the identification of efficient complementarities among the cofactor, the substrate, and the biological partner. Molecular modeling represents an interesting alternative to help in this task. However, little attention has been paid to this field so far. In this manuscript, we aim at reviewing our efforts in developing strategies efficient to computationally drive the design of artificial metalloenzymes. From protein–ligand dockings to multiscale approaches, we intend to demonstrate that modeling could be useful at the different steps of the design. This Perspective ultimately aims at providing computational chemists with illustration of the applications of their tools for artificial metalloenzymes and convincing enzyme designers of the capabilities, qualitative and quantitative, of computational methodologies.