Promiscuous Ability of Human Carbonic Anhydrase: QM and QM/MM Investigation of Carbon Dioxide and Carbodiimide Hydration

The mechanistic insights into different aspects of promiscuity in metalloenzymes

Enzymes are nature's ultimate machinery to catalyze complex reactions. Though enzymes are evolved to catalyze specific reactions, they also show significant promiscuity in reactions and substrate selection. Metalloenzymes contain a metal ion or metal cofactor in their active site, which is crucial in their catalytic activity. Depending on the metal and its coordination environment, the metal ion or cofactor may function as a Lewis acid or base and a redox center and thus can catalyze a plethora of natural reactions. In fact, the versatility in the oxidation state of the metal ions provides metalloenzymes with a high level of catalytic adaptability and promiscuity. In this chapter, we discuss different aspects of promiscuity in metalloenzymes by using several recent experimental and theoretical works as case studies. We start our discussion by introducing the concept of promiscuity and then we delve into the mechanistic insight into promiscuity at the molecular level.

Viral RNA Replication Suppression of SARS-CoV-2: Atomistic Insights into Inhibition Mechanisms of RdRp Machinery by ddhCTP

The nonstructural protein 12, known as RNA-dependent RNA polymerase (RdRp), is essential for both replication and repair of the viral genome. The RdRp of SARS-CoV-2 has been used as a promising candidate for drug development since the inception of the COVID-19 spread. In this work, we performed an in silico investigation on the insertion of the naturally modified pyrimidine nucleobase ddhCTP into the SARS-CoV-2 RdRp active site, in a comparative analysis with the natural one (CTP). The modification in ddhCTP involves the removal of the 3'-hydroxyl group that prevents the addition of subsequent nucleotides into the nascent strand, acting as an RNA chain terminator inhibitor. Quantum mechanical investigations helped to shed light on the mechanistic source of RdRp activity on the selected nucleobases, and comprehensive all-atom simulations provided insights about the structural rearrangements occurring in the active-site region when inorganic pyrophosphate (PPi) is formed. Subsequently, the intricate pathways for the release of PPi, the catalytic product of RdRp, were investigated using Umbrella Sampling simulations. The results are in line with the available experimental data and contribute to a more comprehensive point of view on such an important viral enzyme.

The Se-S Bond Formation in the Covalent Inhibition Mechanism of SARS-CoV-2 Main Protease by Ebselen-like Inhibitors: A Computational Study

The inhibition mechanism of the main protease (Mpro) of SARS-CoV-2 by ebselen (EBS) and its analog with a hydroxyl group at position 2 of the benzisoselenazol-3(2H)-one ring (EBS-OH) was studied by using a density functional level of theory. Preliminary molecular dynamics simulations on the apo form of Mpro were performed taking into account both the hydrogen donor and acceptor natures of the Nδ and Nε of His41, a member of the catalytic dyad. The potential energy surfaces for the formation of the Se-S covalent bond mediated by EBS and EBS-OH on Mpro are discussed in detail. The EBS-OH shows a distinctive behavior with respect to EBS in the formation of the noncovalent complex. Due to the presence of canonical H-bonds and noncanonical ones involving less electronegative atoms, such as sulfur and selenium, the influence on the energy barriers and reaction energy of the Minnesota hybrid meta-GGA functionals M06, M06-2X and M08HX, and the more recent range-separated hybrid functional wB97X were also considered. The knowledge of the inhibition mechanism of Mpro by the small protease inhibitors EBS or EBS-OH can enlarge the possibilities for designing more potent and selective inhibitor-based drugs to be used in combination with other antiviral therapies.

Conformational Change of H64 and Substrate Transportation: Insight Into a Full Picture of Enzymatic Hydration of CO2 by Carbonic Anhydrase

The enzymatic hydration of CO₂ into HCO₃ ^- by carbonic anhydrase (CA) is highly efficient and environment-friendly measure for CO₂ sequestration. Here extensive MM MD and QM/MM MD simulations were used to explore the whole enzymatic process, and a full picture of the enzymatic hydration of CO₂ by CA was achieved. Prior to CO₂ hydration, the proton transfer from the water molecule (WT1) to H64 is the rate-limiting step with the free energy barrier of 10.4 kcal/mol, which leads to the ready state with the Zn-bound OH^-. The nucleophilic attack of OH^- on CO₂ produces HCO₃ ^- with the free energy barrier of 4.4 kcal/mol and the free energy release of about 8.0 kcal/mol. Q92 as the key residue manipulates both CO₂ transportation to the active site and release of HCO₃ ^-. The unprotonated H64 in CA prefers in an inward orientation, while the outward conformation is favorable energetically for its protonated counterpart. The conformational transition of H64 between inward and outward correlates with its protonation state, which is mediated by the proton transfer and the product release. The whole enzymatic cycle has the free energy span of 10.4 kcal/mol for the initial proton transfer step and the free energy change of -6.5 kcal/mol. The mechanistic details provide a comprehensive understanding of the entire reversible conversion of CO₂ into bicarbonate and roles of key residues in chemical and nonchemical steps for the enzymatic hydration of CO₂.

Computing Metal-Binding Proteins for Therapeutic Benefit

Over one third of biomolecules rely on metal ions to exert their cellular functions. Metal ions can play a structural role by stabilizing the structure of biomolecules, a functional role by promoting a wide variety of biochemical reactions, and a regulatory role by acting as messengers upon binding to proteins regulating cellular metal-homeostasis. These diverse roles in biology ascribe critical implications to metal-binding proteins in the onset of many diseases. Hence, it is of utmost importance to exhaustively unlock the different mechanistic facets of metal-binding proteins and to harness this knowledge to rationally devise novel therapeutic strategies to prevent or cure pathological states associated with metal-dependent cellular dysfunctions. In this compendium, we illustrate how the use of a computational arsenal based on docking, classical, and quantum-classical molecular dynamics simulations can contribute to extricate the minutiae of the catalytic, transport, and inhibition mechanisms of metal-binding proteins at the atomic level. This knowledge represents a fertile ground and an essential prerequisite for selectively targeting metal-binding proteins with small-molecule inhibitors aiming to (i) abrogate deregulated metal-dependent (mis)functions or (ii) leverage metal-dyshomeostasis to selectively trigger harmful cells death.

Frontiers of metal-coordinating drug design

Introduction: The occurrence of metal ions in biomolecules is required to exert vital cellular functions. Metal-containing biomolecules can be modulated by small-molecule inhibitors targeting their metal-moiety. As well, the discovery of cisplatin ushered the rational discovery of metal-containing-drugs. The use of both drug types exploiting metal-ligand interactions is well established to treat distinct pathologies. Therefore, characterizing and leveraging metal-coordinating drugs is a pivotal, yet challenging, part of medicinal chemistry.Area covered: Atomic-level simulations are increasingly employed to overcome the challenges met by traditional drug-discovery approaches and to complement wet-lab experiments in elucidating the mechanisms of drugs' action. Multiscale simulations, allow deciphering the mechanism of metal-binding inhibitors and metallo-containing-drugs, enabling a reliable description of metal-complexes in their biological environment. In this compendium, the authors review selected applications exploiting the metal-ligand interactions by focusing on understanding the mechanism and design of (i) inhibitors targeting iron and zinc-enzymes, and (ii) ruthenium and gold-based anticancer agents targeting the nucleosome and aquaporin protein, respectively.Expert opinion: The showcased applications exemplify the current role and the potential of atomic-level simulations and reveal how their synergic use with experiments can contribute to uncover fundamental mechanistic facets and exploit metal-ligand interactions in medicinal chemistry.

On the Catalytic Activity of the Engineered Coiled-Coil Heptamer Mimicking the Hydrolase Enzymes: Insights from a Computational Study

Recently major advances were gained on the designed proteins aimed to generate biomolecular mimics of proteases. Although such enzyme-like catalysts must still suffer refinements for improving the catalytic activity, at the moment, they represent a good example of artificial enzymes to be tested in different fields. Herein, a de novo designed homo-heptameric peptide assembly (CC-Hept) where the esterase activity towards p-nitro-phenylacetate was obtained for introduction of the catalytic triad (Cys-His-Glu) into the hydrophobic matrix, is the object of the present combined molecular dynamics and quantum mechanics/molecular mechanics investigation. Constant pH Molecular Dynamics simulations on the apoform of CC-Hept suggested that the Cys residues are present in the protonated form. Molecular dynamics (MD) simulations of the enzyme-substrate complex evidenced the attitude of the enzyme-like system to retain water molecules, necessary in the hydrolytic reaction, in correspondence of the active site, represented by the Cys-His-Glu triad on each of the seven chains, without significant structural perturbations. A detailed reaction mechanism of esterase activity of CC-Hept-Cys-His-Glu was investigated on the basis of the quantum mechanics/molecular mechanics calculations employing a large quantum mechanical (QM) region of the active site. The proposed mechanism is consistent with available esterases kinetics and structural data. The roles of the active site residues were also evaluated. The deacylation phase emerged as the rate-determining step, in agreement with esterase activity of other natural proteases.

The Catalytic Mechanism of Human Transketolase

We have computationally determined the catalytic mechanism of human transketolase (hTK) using a cluster model approach and density functional theory calculations. We were able to determine all the relevant structures, bringing solid evidences to the proposed experimental mechanism, and to add important detail to the structure of the transition states and the energy profile associated with catalysis. Furthermore, we have established the existence of a crucial intermediate of the catalytic cycle, in agreement with experiments. The calculated data brought new insights to hTK's catalytic mechanism, providing free-energy values for the chemical reaction, as well as adding atomistic detail to the experimental mechanism.

Principles Governing Catalytic Activity of Self-Assembled Short Peptides

Molecular self-assembly provides a chemical strategy for the synthesis of nanostructures by using the principles of nature, and peptides serve as the promising building blocks to construct adaptable molecular architectures. Recently, a series of heptapeptides with alternative hydrophobic and hydrophilic residues were reported to form amyloid-like structures, which are capable of catalyzing acyl ester hydrolysis with remarkable efficiency. However, information remains elusive about the atomic structures of the fibrils. What is the origin of the sequence-dependent catalytic activity? How is the ester hydrolysis catalyzed by the fibrils? In this work, the atomic structures of the aggregates were determined by using molecular modeling and further validated by solid-state NMR experiments, where the fibril with high activity adopts twisted parallel configuration within each layer, and the one with low activity is in flat antiparallel configuration. The polymorphism originates from the interactions between different regions of the building block peptides, where the delicate balance between rigidity and flexibility plays an important role. We further show that the p-nitrophenylacetate ( pNPA) hydrolysis reactions catalyzed by two different fibrils follow a similar mechanism, and the difference in microenvironment at the active site between the natural enzyme and the present self-assembled fibrils should account for the discrepancy in catalytic activities. The present work provides understanding of the structure and function of self-assembled fibrils formed with short peptides at an atomic level and thus sheds new insight on designing aggregates with better functions.

Kinetics of CO2 diffusion in human carbonic anhydrase: a study using molecular dynamics simulations and the Markov-state model

Molecular dynamics (MD) simulations, in combination with the Markov-state model (MSM), were applied to probe CO₂ diffusion from an aqueous solution into the active site of human carbonic anhydrase II (hCA-II), an enzyme useful for enhanced CO₂ capture and utilization. The diffusion process in the hydrophobic pocket of hCA-II was illustrated in terms of a two-dimensional free-energy landscape. We found that CO₂ diffusion in hCA-II is a rate-limiting step in the CO₂ diffusion-binding-reaction process. The equilibrium distribution of CO₂ shows its preferential accumulation within a hydrophobic domain in the protein core region. An analysis of the committors and reactive fluxes indicates that the main pathway for CO₂ diffusion into the active site of hCA-II is through a binding pocket where residue Gln¹³⁶ contributes to the maximal flux. The simulation results offer a new perspective on the CO₂ hydration kinetics and useful insights toward the development of novel biochemical processes for more efficient CO₂ sequestration and utilization.

Mechanistic Explanation of the Weak Carbonic Anhydrase's Esterase Activity

In order to elucidate the elementary mechanism of the promiscuous esterase activity of human carbonic anhydrase (h-CA), we present an accurate theoretical investigation on the hydrolysis of fully-acetylated d-glucose functionalized as sulfamate. This h-CA's inhibitor is of potential relevance in cancer therapy. The study has been performed within the framework of three-layer ONIOM (QM-high:QM'-medium:MM-low) hybrid approach. The computations revealed that the hydrolysis process is not energetically favored, in agreement with the observed weak carbonic anhydrase's esterase activity.

The working mechanism of the β-carbonic anhydrase degrading carbonyl sulphide (COSase): a theoretical study

In order to give insights into the working mechanism of the novel characterized enzyme carbonyl sulphide hydrolase (COSase), which efficiently converts COS into H2S and CO2, we have performed a detailed theoretical investigation using the framework of density functional theory (using B3LYP and M06 exchange-correlation functionals) by the cluster model approach. In the final part of the reaction the metal ion is unable to form a pentacoordinated species. The B3LYP-D3 and M06 potential energy surfaces have a very similar shape. The elucidation of the catalytic reduction of COS is important in view of its role in environmental chemistry.