Caffeic Acid-Zinc Basic Salt/Chitosan Nanohybrid Possesses Controlled Release Properties and Exhibits In Vivo Anti-Inflammatory Activities.

Caffeic acid (CA) exhibits a myriad of biological activities including cardioprotective action, antioxidant, antitumor, anti-inflammatory, and antimicrobial properties. On the other hand, CA presents low water solubility and poor bioavailability, which have limited its use for therapeutic applications. The objective of this study was to develop a nanohybrid of zinc basic salts (ZBS) and chitosan (Ch) containing CA (ZBS-CA/Ch) and evaluate its anti-edematogenic and antioxidant activity in dextran and carrageenan-induced paw edema model. The samples were obtained by coprecipitation method and characterized by X-ray diffraction, Fourier transform infrared (FT-IR), scanning electron microscope (SEM) and UV-visible spectroscopy. The release of caffeate anions from ZBS-CA and ZBS-CA/Ch is pH-dependent and is explained by a pseudo-second order kinetics model, with a linear correlation coefficient of R2 ≥ 0.99 at pH 4.8 and 7.4. The in vivo pharmacological assays showed excellent anti-edematogenic and antioxidant action of the ZBS-CA/Ch nanoparticle with slowly releases of caffeate anions in the tissue, leading to a prolongation of CA-induced anti-edematogenic and anti-inflammatory activities, as well as improving its inhibition or sequestration antioxidant action toward reactive species. Overall, this study highlighted the importance of ZBS-CA/Ch as an optimal drug carrier.

Experimental and Computational Study of Aryl-thiosemicarbazones Inhibiting Cruzain Reveals Reversible Inhibition and a Stepwise Mechanism.

Trypanosoma cruzi is a parasite that infects about 6-7 million people worldwide, mostly in Latin America, causing Chagas disease. Cruzain, the main cysteine protease of T. cruzi, is a validated target for developing drug candidates for Chagas disease. Thiosemicarbazones are one of the most relevant warheads used in covalent inhibitors targeting cruzain. Despite its relevance, the mechanism of inhibition of cruzain by thiosemicarbazones is unknown. Here, we combined experiments and simulations to unveil the covalent inhibition mechanism of cruzain by a thiosemicarbazone-based inhibitor (compound 1). Additionally, we studied a semicarbazone (compound 2), which is structurally similar to compound 1 but does not inhibit cruzain. Assays confirmed the reversibility of inhibition by compound 1 and suggested a two-step mechanism of inhibition. The Ki was estimated to be 36.3 µM and Ki* to be 11.5 µM, suggesting the pre-covalent complex to be relevant for inhibition. Molecular dynamics simulations of compounds 1 and 2 with cruzain were used to propose putative binding modes for the ligands. One-dimensional (1D) quantum mechanics/molecular mechanics (QM/MM) potential of mean force (PMF) and gas-phase energies showed that the attack of Cys25-S- on the C═S or C═O bond yields a more stable intermediate than the attack on the C═N bond of the thiosemicarbazone/semicarbazone. Two-dimensional (2D) QM/MM PMF revealed a putative reaction mechanism for compound 1, involving the proton transfer to the ligand, followed by the Cys25-S- attack at C═S. The ΔG and energy barrier were estimated to be -1.4 and 11.7 kcal/mol, respectively. Overall, our results shed light on the inhibition mechanism of cruzain by thiosemicarbazones.

Assessment of Reversibility for Covalent Cysteine Protease Inhibitors Using Quantum Mechanics/Molecular Mechanics Free Energy Surfaces.

We have used molecular dynamics (MD) simulations with hybrid quantum mechanics/molecular mechanics (QM/MM) potentials to investigate the reaction mechanism for covalent inhibition of cathepsin K and assess the reversibility of inhibition. The computed free energy profiles suggest that a nucleophilic attack by the catalytic cysteine on the inhibitor warhead and proton transfer from the catalytic histidine occur in a concerted manner. The results indicate that the reaction is more strongly exergonic for the alkyne-based inhibitors, which bind irreversibly to cathepsin K, than for the nitrile-based inhibitor odanacatib, which binds reversibly. Gas-phase energies were also calculated for the addition of methanethiol to structural prototypes for a number of warheads of interest in cysteine protease inhibitor design in order to assess electrophilicity. The approaches presented in this study are particularly applicable to assessment of novel warheads, and computed transition state geometries can be incorporated into molecular models for covalent docking.

A patent review on cathepsin K inhibitors to treat osteoporosis (2011 - 2021).

INTRODUCTION: Cathepsin K (CatK) is a lysosomal cysteine protease and the predominant cathepsin expressed in osteoclasts, where it degrades the bone matrix. Hence, CatK is an attractive therapeutic target related to diseases characterized by bone resorption, like osteoporosis. AREAS COVERED: This review summarizes the patent literature from 2011 to 2021 on CatK inhibitors and their potential use as new treatments for osteoporosis. The inhibitors were classified by their warheads, with the most explored nitrile-based inhibitors. Promising in vivo results have also been disclosed. EXPERT OPINION: As one of the most potent lysosomal proteins whose primary function is to mediate bone resorption, cathepsin K remains an excellent target for therapeutic intervention. Nevertheless, there is no record of any approved drug that targets CatK. The most notable cases of drug candidates targeting CatK were balicatib and odanacatib, which reached Phase II and III clinical trials, respectively, but did not enter the market. Further developments include exploring new chemical entities beyond the nitrile-based chemical space, with improved ADME and safety profiles. In addition, CatK's role in cancer immunoexpression and its involvement in the pathophysiology of osteo- and rheumatoid arthritis have raised the race to develop activity-based probes with excellent potency and selectivity.

Assessment of the PETase conformational changes induced by poly(ethylene terephthalate) binding.

Recently, a bacterium strain of Ideonella sakaiensis was identified with the uncommon ability to degrade the poly(ethylene terephthalate) (PET). The PETase from I. sakaiensis strain 201-F6 (IsPETase) catalyzes the hydrolysis of PET converting it to mono(2-hydroxyethyl) terephthalic acid (MHET), bis(2-hydroxyethyl)-TPA (BHET), and terephthalic acid (TPA). Despite the potential of this enzyme for mitigation or elimination of environmental contaminants, one of the limitations of the use of IsPETase for PET degradation is the fact that it acts only at moderate temperature due to its low thermal stability. Besides, molecular details of the main interactions of PET in the active site of IsPETase remain unclear. Herein, molecular docking and molecular dynamics (MD) simulations were applied to analyze structural changes of IsPETase induced by PET binding. Results from the essential dynamics revealed that the ß1-ß2 connecting loop is very flexible. This loop is located far from the active site of IsPETase and we suggest that it can be considered for mutagenesis to increase the thermal stability of IsPETase. The free energy landscape (FEL) demonstrates that the main change in the transition between the unbound to the bound state is associated with the ß7-α5 connecting loop, where the catalytic residue Asp206 is located. Overall, the present study provides insights into the molecular binding mechanism of PET into the IsPETase structure and a computational strategy for mapping flexible regions of this enzyme, which can be useful for the engineering of more efficient enzymes for recycling plastic polymers using biological systems.

QM/MM Study of the Fosfomycin Resistance Mechanism Involving FosB Enzyme.

Multidrug-resistant organisms contain antibiotic-modifying enzymes that facilitate resistance to a variety of antimicrobial compounds. Particularly, the fosfomycin (FOF) drug can be structurally modified by several FOF-modifying enzymes before it reaches the biological target. Among them, FosB is an enzyme that utilizes l-cysteine or bacillithiol in the presence of a divalent metal to open the epoxide ring of FOF and, consequently, inactivate the drug. Here, we have used hybrid quantum mechanics/molecular mechanics (QM/MM) and molecular dynamics (MD) simulations to explore the mechanism of the reaction involving FosB and FOF. The calculated free-energy profiles show that the cost to open the epoxide ring of FOF at the C2 atom is ∼3.0 kcal/mol higher than that at the C1 atom. Besides, our QM/MM MD results revealed the critical role of conformation change of Cys9 and Asn50 to release the drug from the active site. Overall, the present study provides insights into the mechanism of FOF-resistant proteins.

Design, synthesis and stepwise optimization of nitrile-based inhibitors of cathepsins B and L.

Human cathepsin B (CatB) is an important biological target in cancer therapy. In this work, we performed a knowledge-based design approach and the synthesis of a new set of 19 peptide-like nitrile-based cathepsin inhibitors. Reported compounds were assayed against a panel of human cysteine proteases: CatB, CatL, CatK, and CatS. Three compounds (7h, 7i, and 7j) displayed nanomolar inhibition of CatB and selectivity over CatK and CatL. The selectivity was achieved by using the combination of a para biphenyl ring at P3, halogenated phenylalanine in P2 and Thr-O-Bz group at P1. Likewise, compounds 7i and 7j showed selective CatB inhibition among the panel of enzymes studied. We have also described a successful example of bioisosteric replacement of the amide bond for a sulfonamide one [7e â†’ 6b], where we observed an increase in affinity and selectivity for CatB while lowering the compound lipophilicity (ilogP). Our knowledge-based design approach and the respective structure-activity relationships provide insights into the specific ligand-target interactions for therapeutically relevant cathepsins.

Facile Synthesis and Metabolic Incorporation of m-DAP Bioisosteres Into Cell Walls of Live Bacteria.

Bacterial cell walls contain peptidoglycan (PG), a scaffold that provides proper rigidity to resist lysis from internal osmotic pressure and a barrier to protect cells against external stressors. It consists of repeating sugar units with a linkage to a stem peptide that becomes cross-linked by cell wall transpeptidases (TP). While synthetic PG fragments containing l-lysine in the third position on the stem peptide are easier to access, those with meso-diaminopimelic acid (m-DAP) pose a severe synthetic challenge. Herein, we describe a solid phase synthetic scheme based on widely available building blocks to assemble meso-cystine (m-CYT), which mimics key structural features of m-DAP. To demonstrate proper mimicry of m-DAP, cell wall probes were synthesized with m-CYT in place of m-DAP and evaluated for their metabolic processing in live bacterial cells. We found that m-CYT-based cell wall probes were properly processed by TPs in various bacterial species that endogenously contain m-DAP in their PG. Additionally, we have used hybrid quantum mechanical/molecular mechanical (QM/MM) and molecular dynamics (MD) simulations to explore the influence of m-DAP analogs on the PG cross-linking. The results showed that the cross-linking mechanism of transpeptidases occurred through a concerted process. We anticipate that this strategy, which is based on the use of inexpensive and commercially available building blocks, can be widely adopted to provide greater accessibility of PG mimics for m-DAP containing organisms.

Crystal structure of Leishmania mexicana cysteine protease B in complex with a high-affinity azadipeptide nitrile inhibitor.

Leishmania mexicana is an obligate intracellular protozoan parasite that causes the cutaneous form of leishmaniasis affecting South America and Mexico. The cysteine protease LmCPB is essential for the virulence of the parasite and therefore, it is an appealing target for antiparasitic therapy. A library of nitrile-based cysteine protease inhibitors was screened against LmCPB to develop a treatment of cutaneous leishmaniasis. Several compounds are sufficiently high-affinity LmCPB inhibitors to serve both as starting points for drug discovery projects and as probes for target validation. A 1.4 Å X ray crystal structure, the first to be reported for LmCPB, was determined for the complex of this enzyme covalently bound to an azadipeptide nitrile ligand. Mapping the structure-activity relationships for LmCPB inhibition revealed superadditive effects for two pairs of structural transformations. Therefore, this work advances our understanding of azadipeptidyl and dipeptidyl nitrile structure-activity relationships for LmCPB structure-based inhibitor design. We also tested the same series of inhibitors on related cysteine proteases cathepsin L and Trypanosoma cruzi cruzain. The modulation of these mammalian and protozoan proteases represents a new framework for targeting papain-like cysteine proteases.

Optimization strategy of single-digit nanomolar cross-class inhibitors of mammalian and protozoa cysteine proteases.

Cysteine proteases (CPs) are involved in a myriad of actions that include not only protein degradation, but also play an essential biological role in infectious and systemic diseases such as cancer. CPs also act as biomarkers and can be reached by active-based probes for diagnostic and mechanistic purposes that are critical in health and disease. In this paper, we present the modulation of a CP panel of parasites and mammals (Trypanosoma cruzi cruzain, LmCPB, CatK, CatL and CatS), whose inhibition by nitrile peptidomimetics allowed the identification of specificity and selectivity for a given CP. The activity cliffs identified at the CP inhibition level are useful for retrieving trends through multiple structure-activity relationships. For two of the cruzain inhibitors (10g and 4e), both enthalpy and entropy are favourable to Gibbs binding energy, thus overcoming enthalpy-entropy compensation (EEC). Group contribution of individual molecular modification through changes in enthalpy and entropy results in a separate partition on the relative differences of Gibbs binding energy (ΔΔG). Overall, this study highlights the role of CPs in polypharmacology and multi-target screening, which represents an imperative trend in the actual drug discovery effort.

Assessment of the Cruzain Cysteine Protease Reversible and Irreversible Covalent Inhibition Mechanism.

Reversible and irreversible covalent ligands are advanced cysteine protease inhibitors in the drug development pipeline. K777 is an irreversible inhibitor of cruzain, a necessary enzyme for the survival of the Trypanosoma cruzi (T. cruzi) parasite, the causative agent of Chagas disease. Despite their importance, irreversible covalent inhibitors are still often avoided due to the risk of adverse effects. Herein, we replaced the K777 vinyl sulfone group with a nitrile moiety to obtain a reversible covalent inhibitor (Neq0682) of cysteine protease. Then, we used advanced experimental and computational techniques to explore details of the inhibition mechanism of cruzain by reversible and irreversible inhibitors. The isothermal titration calorimetry (ITC) analysis shows that inhibition of cruzain by an irreversible inhibitor is thermodynamically more favorable than by a reversible one. The hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) and Molecular Dynamics (MD) simulations were used to explore the mechanism of the reaction inhibition of cruzain by K777 and Neq0682. The calculated free energy profiles show that the Cys25 nucleophilic attack and His162 proton transfer occur in a single step for a reversible inhibitor and two steps for an irreversible covalent inhibitor. The hybrid QM/MM calculated free energies for the inhibition reaction correspond to -26.7 and -5.9 kcal mol-1 for K777 and Neq0682 at the MP2/MM level, respectively. These results indicate that the ΔG of the reaction is very negative for the process involving K777, consequently, the covalent adduct cannot revert to a noncovalent protein-ligand complex, and its binding tends to be irreversible. Overall, the present study provides insights into a covalent inhibition mechanism of cysteine proteases.

Evaluating QM/MM Free Energy Surfaces for Ranking Cysteine Protease Covalent Inhibitors.

One tactic for cysteine protease inhibition is to form a covalent bond between an electrophilic atom of the inhibitor and the thiol of the catalytic cysteine. In this study, we evaluate the reaction free energy obtained from a hybrid quantum mechanical/molecular mechanical (QM/MM) free energy profile as a predictor of affinity for reversible, covalent inhibitors of rhodesain. We demonstrate that the reaction free energy calculated with the PM6/MM potential is in agreement with the experimental data and suggest that the free energy profile for covalent bond formation in a protein environment may be a useful tool for the inhibitor design.

On the intrinsic reactivity of highly potent trypanocidal cruzain inhibitors.

The cysteine protease cruzipain is considered to be a validated target for therapeutic intervention in the treatment of Chagas disease. Hence, peptidomimetic cruzipain inhibitors having a reactive group (known as warhead) are subject to continuous studies to discover novel antichagasic compounds. Here, we evaluated how different warheads for a set of structurally similar related compounds could inhibit the activity of cruzipain and, ultimately, their trypanocidal effect. We first investigated in silico the intrinsic reactivity of these compounds by applying the Fukui index to correlate it with the enzymatic affinity. Then, we evaluated their potency against T. cruzi (Y and Tulahuen strains), which revealed the reversible cruzain inhibitor Neq0656 as a better trypanocidal agent (ECY.strain 50 = 0.1 µM; SI = 58.4) than the current drug benznidazole (ECY.strain 50 = 5.1 µM; SI > 19.6). We also measured the half-life time by HPLC analysis of three lead compounds in the presence of glutathione and cysteine to experimentally assess their intrinsic reactivity. Results clearly illustrated the reactivity trend for the warheads (azanitrile > aldehyde > nitrile), where the aldehyde displayed an intermediate intrinsic reactivity. Therefore, the aldehyde bearing peptidomimetic compounds should be subject for in-depth evaluation in the drug discovery process.

Targeting Peptidyl-prolyl Cis-trans Isomerase NIMA-interacting 1: A Structure-based Virtual Screening Approach to Find Novel Inhibitors.

BACKGROUND: Peptidyl-prolyl cis-trans isomerase NIMA-interacting 1 (Pin1) is an enzyme that isomerizes phosphorylated serine or threonine motifs adjacent to proline residues. Pin1 has important roles in several cellular signaling pathways, consequently impacting the development of multiple types of cancers. METHODS: Based on the previously reported inhibitory activity of pentacyclic triterpenoids isolated from the gum resin of Boswellia genus against Pin1, we designed a computational experiment using molecular docking, pharmacophore filtering, and structural clustering allied to molecular dynamics (MD) simulations and binding free energy calculations to explore the inhibitory activity of new triterpenoids against Pin1 structure. RESULTS: Here, we report different computational evidence that triterpenoids from neem (Azadirachta indica A. Juss), such as 6-deacetylnimbinene, 6-Oacetylnimbandiol, and nimbolide, replicate the binding mode of the Pin1 substrate peptide, interacting with high affinity with the binding site and thus destabilizing the Pin1 structure. CONCLUSIONS: Our results are supported by experimental data, and provide interesting structural insights into their molecular mechanism of action, indicating that their structural scaffolds could be used as a start point to develop new inhibitors against Pin1.

In silico identification of natural products with anticancer activity using a chemo-structural database of Brazilian biodiversity.

Cancer is one of the leading causes of death worldwide, and the number of patients has only increased each year, despite the considerable efforts and investments in scientific research. Since natural products (NPs) may serve as suitable sources for drug development, the cytotoxicity against cancer cells of 2221 compounds from the Nuclei of Bioassays, Ecophysiology, and Biosynthesis of Natural Products Database (NuBBEDB) was predicted using CDRUG algorithm. Molecular modeling, chemoinformatics, and chemometric tools were then used to analyze the structural and physicochemical properties of these compounds. We compared the positive NPs with FDA-approved anticancer drugs and predicted the molecular targets involved in the anticancer activity. In the present study, 46 families comprising potential anticancer compounds and at least 19 molecular targets involved in oncogenesis. To the best of our knowledge, this is the first large-scale study conducted to evaluate the potentiality of NPs sourced from Brazilian biodiversity as anticancer agents, using in silico approaches. Our results provided interesting insights about the mechanism of action of these compounds, and also suggested that their structural diversity may aid structure-based optimization strategies for developing novel drugs for cancer therapy.

Understanding the Mechanism of Direct Activation of AMP-Kinase: Toward a Fine Allosteric Tuning of the Kinase Activity.

Mammalian AMP-activated protein kinase (AMPK) is a Ser/Thr protein kinase with a key role as a sensor in cellular energy homeostasis. It has a major role in numerous metabolic disorders, such as type 2 diabetes, obesity, and cancer, and hence it has gained progressive interest as a potential therapeutic target. AMPK is a heterotrimeric enzyme composed by an α-catalytic subunit and two regulatory subunits, ß and γ. It is regulated by several mechanisms, including indirect activators such as metformin and direct activators such as compound A-769662. The crystal structure of AMPK bound to A-769662 has been recently reported, suggesting a hypothetical allosteric mechanism of AMPK activation assisted by phosphorylated Ser108 at the ß-subunit. Here, we have studied the direct activation mechanism of A-769662 by means of molecular dynamics simulations, suggesting that the activator may act as a glue, coupling the dynamical motion of the ß-subunit and the N-terminal domain of the α-subunit, and assisting the preorganization of the ATP-binding site. This is achieved through the formation of an allosteric network that connects the activator and ATP-binding sites, particularly through key interactions formed between αAsp88 and ßArg83 and between ßpSer108 and αLys29. Overall, these studies shed light into key mechanistic determinants of the allosteric regulation of this cellular energy sensor, and pave the way for the fine-tuning of the rational design of direct activators of this cellular energy sensor.

Experimental study and computational modelling of cruzain cysteine protease inhibition by dipeptidyl nitriles.

Chagas disease affects millions of people in Latin America. This disease is caused by the protozoan parasite Trypanossoma cruzi. The cysteine protease cruzain is a key enzyme for the survival and propagation of this parasite lifecycle. Nitrile-based inhibitors are efficient inhibitors of cruzain that bind by forming a covalent bond with this enzyme. Here, three nitrile-based inhibitors dubbed Neq0409, Neq0410 and Neq0570 were synthesized, and the thermodynamic profile of the bimolecular interaction with cruzain was determined using isothermal titration calorimetry (ITC). The result suggests the inhibition process is enthalpy driven, with a detrimental contribution of entropy. In addition, we have used hybrid Quantum Mechanical/Molecular Mechanical (QM/MM) and Molecular Dynamics (MD) simulations to investigate the reaction mechanism of reversible covalent modification of cruzain by Neq0409, Neq0410 and Neq0570. The computed free energy profile shows that the nucleophilic attack of Cys25 on the carbon C1 of inhibitiors and the proton transfer from His162 to N1 of the dipeptidyl nitrile inhibitor take place in a single step. The calculated free energy of the inhibiton reaction is in agreement with covalent experimental binding. Altogether, the results reported here suggests that nitrile-based inhibitors are good candidates for the development of reversible covalent inhibitors of cruzain and other cysteine proteases.

First homology model of Plasmodium falciparum glucose-6-phosphate dehydrogenase: Discovery of selective substrate analog-based inhibitors as novel antimalarial agents.

In Plasmodium falciparum the bifunctional enzyme glucose-6-phosphate dehydrogenase‒6-phosphogluconolactonase (PfG6PD‒6PGL) is involved in the catalysis of the first reaction of the pentose phosphate pathway. Since this enzyme has a key role in parasite development, its unique structure represents a potential target for the discovery of antimalarial drugs. Here we describe the first 3D structural model of the G6PD domain of PfG6PD‒6PGL. Compared to the human enzyme (hG6PD), the 3D model has enabled the identification of a key difference in the substrate-binding site, which involves the replacement of Arg365 in hG6PD by Asp750 in PfG6PD. In a prospective validation of the model, this critical change has been exploited to rationally design a novel family of substrate analog-based inhibitors that can display the necessary selectivity towards PfG6PD. A series of glucose derivatives featuring an α-methoxy group at the anomeric position and different side chains at position 6 bearing distinct basic functionalities has been synthesized, and their PfG6PD and hG6PD inhibitory activities and their toxicity against parasite and mammalian cells have been assessed. Several compounds displayed micromolar affinity (Ki up to 23 µM), favorable selectivity (up to > 26-fold), and low cytotoxicity. Phenotypic assays with P. falciparum cultures revealed high micromolar IC50 values, likely as a result of poor internalization of the compounds in the parasite cell. Overall, these results endorse confidence to the 3D model of PfG6PD, paving the way for the use of target-based drug design approaches in antimalarial drug discovery studies around this promising target.

Computational analyses of interactions between ALK-5 and bioactive ligands: insights for the design of potential anticancer agents.

Activin Receptor-Like Kinase 5 (ALK-5) is related to some types of cancer, such as breast, lung, and pancreas. In this study, we have used molecular docking, molecular dynamics simulations, and free energy calculations in order to explore key interactions between ALK-5 and six bioactive ligands with different ranges of biological activity. The motivation of this work is the lack of crystal structure for inhibitor-protein complexes for this set of ligands. The understanding of the molecular structure and the protein-ligand interaction could give support for the development of new drugs against cancer. The results show that the calculated binding free energy using MM-GBSA, MM-PBSA, and SIE is correlated with experimental data with r2 = 0.88, 0.80, and 0.94, respectively, which indicates that the calculated binding free energy is in excellent agreement with experimental data. In addition, the results demonstrate that H bonds with Lys232, Glu245, Tyr249, His283, Asp351, and one structural water molecule play an important role for the inhibition of ALK-5. Overall, we discussed the main interactions between ALK-5 and six inhibitors that may be used as starting points for designing new molecules to the treatment of cancer.

Catalysis by solvation rather than the desolvation effect: exploring the catalytic efficiency of SAM-dependent chlorinase.

Chlorinase SalL halogenate S-adenosyl-l-methionine (SAM) reacts with chloride to generate 5'-chloro-5'-deoxyadenosine and l-methionine through a nucleophilic substitution mechanism. Although it is known that chlorinase enhances the rate of reaction by a factor of 1.2 × 1017 fold, it is not entirely clear how this is accomplished. The search for the origin of the catalysis of chlorinase and other enzymes has led to a desolvation hypothesis. In the present work, we have used well defined computational simulations in order to evaluate the origin of the catalytic efficiency of chlorinase. The results demonstrate that the catalytic effect of chlorinase is associated with the fact that Cl- is "solvated" by the protein more than by the reference solution reaction, which is not in accordance with proposed catalysis by desolvation. It is found that chlorinase SalL active sites provide electrostatic stabilization of the transition state which is the origin of its catalytic effect.