Integrated Computational Study of the Cu-Catalyzed Hydration of Alkenes in Water Solvent and into the Context of an Artificial Metallohydratase

Rationalization of a Streptavidin Based Enantioselective Artificial Suzukiase: An Integrative Computational Approach

An Artificial Metalloenzyme (ArM) built employing the streptavidin-biotin technology has been used for the enantioselective synthesis of binaphthyls by means of asymmetric Suzuki-Miyaura cross-coupling reactions. Despite its success, it remains a challenge to understand how the length of the biotin cofactors or the introduction of mutations to streptavidin leads the preferential synthesis of one atropisomer over the other. In this study, we apply an integrated computational modeling approach, including DFT calculations, protein-ligand dockings and molecular dynamics to rationalize the impact of mutations and length of the biotion cofactor on the enantioselectivities of the biaryl product. The results unravel that the enantiomeric differences found experimentally can be rationalized by the disposition of the first intermediate, coming from the oxidative addition step, and the entrance of the second substrate. The work also showcases the difficulties facing to control the enantioselection when engineering ArM to catalyze enantioselective Suzuki-Miyaura couplings and how the combination of DFT calculations, molecular dockings and MD simulations can be used to rationalize artificial metalloenzymes.

The Impact of DFT Functional, Cluster Model Size, and Implicit Solvation on the Structural Description of Single-Metal-Mediated DNA Phosphodiester Bond Cleavage: The Case Study of APE1

Phosphodiester bond hydrolysis in nucleic acids is a ubiquitous reaction that can be facilitated by enzymes called nucleases, which often use metal ions to achieve catalytic function. While a two-metal-mediated pathway has been well established for many enzymes, there is growing support that some enzymes require only one metal for the catalytic step. Using human apurinic/apyrimidinic endonuclease (APE1) as a prototypical example and cluster models, this study clarifies the impact of DFT functional, cluster model size, and implicit solvation on single-metal-mediated phosphodiester bond cleavage and provides insight into how to efficiently model this chemistry. Initially, a model containing 69 atoms built from a high-resolution X-ray crystal structure is used to explore the reaction pathway mapped by a range of DFT functionals and basis sets, which provides support for the use of standard functionals (M06-2X and B3LYP-D3) to study this reaction. Subsequently, systematically increasing the model size to 185 atoms by including additional amino acids and altering residue truncation points highlights that small models containing only a few amino acids or β carbon truncation points introduce model strains and lead to incorrect metal coordination. Indeed, a model that contains all key residues (general base and acid, residues that stabilize the substrate, and amino acids that maintain the metal coordination) is required for an accurate structural depiction of the one-metal-mediated phosphodiester bond hydrolysis by APE1, which results in 185 atoms. The additional inclusion of the broader enzyme environment through continuum solvation models has negligible effects. The insights gained in the present work can be used to direct future computational studies of other one-metal-dependent nucleases to provide a greater understanding of how nature achieves this difficult chemistry.

Copper-Containing Artificial Polyenzymes as a Clickase for Bioorthogonal Chemistry

Artificial polyenzymes (ArPoly) are tailored combinations of universal protein scaffolds and polymers newly proposed as promising alternatives to natural enzymes to expand the biocatalyst toolbox. The concept of ArPoly has been continuously extended to metal-containing ArPoly to overcome the drawbacks faced by conventional artificial metalloenzymes. Herein, we present a sustainable route to synthesize a novel water-soluble metalloenzyme for copper-catalyzed azide-alkyne cycloadditions in water with remarkable selectivity. In this case, synthetic l-proline monomers were polymerized onto bovine serum albumen in an aqueous medium via copper-mediated "grafting-from" atom-transfer radical polymerization, resulting in protein-polymer-copper conjugates named ArPoly_clickase. The copper in ArPoly_clickase plays pivotal bifunctional roles, not only as the catalyst for polymerization but also as the coordinated active site for alkyne-azide click catalysis. ArPoly_clickase showcases high efficiency, substrate generality, regioselectivity, and ease of product separation for "click chemistry" in water. Notably, ArPoly_clickase displays good biocompatibility without imposing copper toxicity on living cells, which offers the prospect for the upcoming bioorthogonal chemistry.

Rationalizing the Decavanadate(V) and Oxidovanadium(IV) Binding to G-Actin and the Competition with Decaniobate(V) and ATP

The experimental data collected over the past 15 years on the interaction of decavanadate(V) (V₁₀O₂₈^6-; V₁₀), a polyoxometalate (POM) with promising anticancer and antibacterial action, with G-actin, were rationalized by using several computational approaches (docking, density functional theory (DFT), and molecular dynamics (MD)). Moreover, a comparison with the isostructural and more stable decaniobate(V) (Nb₁₀O₂₈^6-; Nb₁₀) was carried out. Four binding sites were identified, named α, β, γ, and δ, the site α being the catalytic nucleotide site located in the cleft of the enzyme at the interface of the subdomains II and IV. It was observed that the site α is preferred by V₁₀, whereas Nb₁₀ is more stable at the site β; this indicates that, differently from other proteins, G-actin could contemporaneously bind the two POMs, whose action would be synergistic. Both decavanadate and decaniobate induce conformational rearrangements in G-actin, larger for V₁₀ than Nb₁₀. Moreover, the binding mode of oxidovanadium(IV) ion, V^IVO²⁺, formed upon the reduction of decavanadate(V) by the -SH groups of accessible cysteine residues, is also found in the catalytic site α with (His161, Asp154) coordination; this adduct overlaps significantly with the region where ATP is bound, accounting for the competition between V₁₀ and its reduction product V^IVO²⁺ with ATP, as previously observed by EPR spectroscopy. Finally, the competition with ATP was rationalized: since decavanadate prefers the nucleotide site α, Ca²⁺-ATP displaces V₁₀ from this site, while the competition is less important for Nb₁₀ because this POM shows a higher affinity for β than for site α. A relevant consequence of this paper is that other metallodrug-protein systems, in the absence or presence of eventual inhibitors and/or competition with molecules of the organism, could be studied with the same approach, suggesting important elements for an explanation of the biological data and a rational drug design.

Unlocking the Full Evolutionary Potential of Artificial Metalloenzymes Through Direct Metal-Protein Coordination : A review of recent advances for catalyst development

Generation of artificial metalloenzymes (ArMs) has gained much inspiration from the general understanding of natural metalloenzymes. Over the last decade, a multitude of methods generating transition metal-protein hybrids have been developed and many of these new-to-nature constructs catalyse reactions previously reserved for the realm of synthetic chemistry. This perspective will focus on ArMs incorporating 4d and 5d transition metals. It aims to summarise the significant advances made to date and asks whether there are chemical strategies, used in nature to optimise metal catalysts, that have yet to be fully recognised in the synthetic enzyme world, particularly whether artificial enzymes produced to date fully take advantage of the structural and energetic context provided by the protein. Further, the argument is put forward that, based on precedence, in the majority of naturally evolved metalloenzymes the direct coordination bonding between the metal and the protein scaffold is integral to catalysis. Therefore, the protein can attenuate metal activity by positioning ligand atoms in the form of amino acids, as well as making non-covalent contributions to catalysis, through intermolecular interactions that pre-organise substrates and stabilise transition states. This highlights the often neglected but crucial element of natural systems that is the energetic contribution towards activating metal centres through protein fold energy. Finally, general principles needed for a different approach to the formation of ArMs are set out, utilising direct coordination inspired by the activation of an organometallic cofactor upon protein binding. This methodology, observed in nature, delivers true interdependence between metal and protein. When combined with the ability to efficiently evolve enzymes, new problems in catalysis could be addressed in a faster and more specific manner than with simpler small molecule catalysts.

Designer metalloenzymes for synthetic biology: Enzyme hybrids for catalysis

Combining organometallics and biology has generated broad interest from scientists working on applications from in situ drug release to biocatalysis. Engineered enzymes and biohybrid catalysts (also referred to as artificial enzymes) have introduced a wide range of abiotic chemistry into biocatalysis. Predominantly, this work has concentrated on using these catalysts for single step in vitro reactions. However, the promise of using these hybrid catalysts in vivo and combining them with synthetic biology and metabolic engineering is vast. This report will briefly review recent advances in artificial metalloenzyme design, followed by summarising recent studies that have looked at the use of these hybrid catalysts in vivo and in enzymatic cascades, therefore exploring their potential for synthetic biology.

Molecular Modeling for Artificial Metalloenzyme Design and Optimization

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.