The importance of the active site histidine for the activity of epoxide- or aziridine-based inhibitors of cysteine proteases

Molecular dynamic simulation (MDS) and in vitro cathepsin-B inhibitory activity of decrusin angelate, ibuprofen, and thymol

Attenuation of cathepsin B (CATB) proteolytic activity and/or inhibition serves as a potential therapeutic target in cancer metastasis. Herein, we determined the specificity of FDA approved potential anti-cancer natural flavonoid decursinol angelate (DA), thymol (TH) and a propionic acid derivative ibuprofen (IB), for the inactivation of CATB. We used enzymatic assay, computational and in vitro methods for the identification of the best candidate. Out of these we found DA can inhibit CATB with lowest IC₅₀ measured after one hour of incubation using Z-Phe-Arg-4MβNA (BANA) as a substrate. Docking analysis suggested favorable interaction of DA with the catalytic site residues (GLN23, CYS26, HIS110, HIS111) of CATB (PDB Id: 1HUC) were responsible for the inhibition of its proteolytic activity. Additionally, in vitro quantification with human colorectal carcinoma (HCT 116) revealed, DA rapidly inactivates CATB as compared with commercial synthetic inhibitor CA074 with no cellular toxicity towards normal colon cells (CCD 841).

2-Sulfonylpyrimidines as Privileged Warheads for the Development of S. aureus Sortase A Inhibitors

Staphylococcus aureus is one of the most frequent causes of nosocomial and community-acquired infections, with emerging multiresistant isolates causing a significant burden to public health systems. We identified 2-sulfonylpyrimidines as a new class of potent inhibitors against S. aureus sortase A acting by covalent modification of the active site cysteine 184. Series of derivatives were synthesized to derive structure-activity relationship (SAR) with the most potent compounds displaying low micromolar K_I values. Studies on the inhibition selectivity of homologous cysteine proteases showed that 2-sulfonylpyrimidines reacted efficiently with protonated cysteine residues as found in sortase A, though surprisingly, no reaction occurred with the more nucleophilic cysteine residue from imidazolinium-thiolate dyads of cathepsin-like proteases. By means of enzymatic and chemical kinetics as well as quantum chemical calculations, it could be rationalized that the S_NAr reaction between protonated cysteine residues and 2-sulfonylpyrimidines proceeds in a concerted fashion, and the mechanism involves a ternary transition state with a conjugated base. Molecular docking and enzyme inhibition at variable pH values allowed us to hypothesize that in sortase A this base is represented by the catalytic histidine 120, which could be substantiated by QM model calculation with 4-methylimidazole as histidine analog.

Mechanisms of Proteolytic Enzymes and Their Inhibition in QM/MM Studies

Experimental evidence for enzymatic mechanisms is often scarce, and in many cases inadvertently biased by the employed methods. Thus, apparently contradictory model mechanisms can result in decade long discussions about the correct interpretation of data and the true theory behind it. However, often such opposing views turn out to be special cases of a more comprehensive and superior concept. Molecular dynamics (MD) and the more advanced molecular mechanical and quantum mechanical approach (QM/MM) provide a relatively consistent framework to treat enzymatic mechanisms, in particular, the activity of proteolytic enzymes. In line with this, computational chemistry based on experimental structures came up with studies on all major protease classes in recent years; examples of aspartic, metallo-, cysteine, serine, and threonine protease mechanisms are well founded on corresponding standards. In addition, experimental evidence from enzyme kinetics, structural research, and various other methods supports the described calculated mechanisms. One step beyond is the application of this information to the design of new and powerful inhibitors of disease-related enzymes, such as the HIV protease. In this overview, a few examples demonstrate the high potential of the QM/MM approach for sophisticated pharmaceutical compound design and supporting functions in the analysis of biomolecular structures.

Covalent Inactivation of Mycobacterium tuberculosis Isocitrate Lyase by cis-2,3-Epoxy-Succinic Acid

The isocitrate lyases (ICL1/2) are essential enzymes of Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis. At present, no ICL1/2 inhibitors have progressed to clinical evaluation, despite extensive drug discovery efforts. Herein, we surveyed succinate analogs against ICL1 and found that dicarboxylic acids constrained in their synperiplanar conformations, such as maleic acid, comprise uncompetitive inhibitors of ICL1 and inhibit more potently than their trans-isomers. From this, we identified cis-2,3 epoxysuccinic acid (cis-EpS) as a selective, irreversible covalent inactivator of Mtb ICL1 (kinact/Kinact= (5.0 ± 1.4) × 104 M-1 s-1; Kinact = 200 ± 50 nM), the most potent inactivator of ICL1 yet characterized. Crystallographic and mass spectrometric analysis demonstrated that Cys191 of ICL1 was S-malylated by cis-EpS, and a crystallographic "snapshot" of inactivation lent insight into the chemical mechanism of this inactivation. Proteomic analysis of E. coli lysates showed that cis-EpS selectively labeled plasmid-expressed Mtb ICL1. Consistently, cis-EpS, but not its trans-isomer, inhibited the growth of Mtb under conditions in which ICL function is essential. These findings encourage the development of analogs of cis-2,3-epoxysuccinate as antituberculosis agents.

Targeting the SARS-CoV-2 main protease using FDA-approved Isavuconazonium, a P2-P3 α-ketoamide derivative and Pentagastrin: An in-silico drug discovery approach

The SARS-CoV-2 main protease (M^pro) is an attractive target towards discovery of drugs to treat COVID-19 because of its key role in virus replication. The atomic structure of M^pro in complex with an α-ketoamide inhibitor (Lig13b) is available (PDB ID:6Y2G). Using 6Y2G and the prior knowledge that protease inhibitors could eradicate COVID-19, we designed a computational study aimed at identifying FDA-approved drugs that could interact with M^pro. We searched the DrugBank and PubChem for analogs and built a virtual library containing ∼33,000 conformers. Using high-throughput virtual screening and ligand docking, we identified Isavuconazonium, a ketoamide inhibitor (α-KI) and Pentagastrin as the top three molecules (Lig13b as the benchmark) based on docking energy. The ΔG_bind of Lig13b, Isavuconazonium, α-KI, Pentagastrin was −28.1, −45.7, −44.7, −34.8 kcal/mol, respectively. Molecular dynamics simulation revealed that these ligands are stable within the M^pro active site. Binding of these ligands is driven by a variety of non-bonded interaction, including polar bonds, H-bonds, van der Waals and salt bridges. The overall conformational dynamics of the complexed-M^pro was slightly altered relative to apo-M^pro. This study demonstrates that three distinct classes molecules, Isavuconazonium (triazole), α-KI (ketoamide) and Pentagastrin (peptide) could serve as potential drugs to treat patients with COVID-19.

•

Using computational modelling, we show that SARS-CoV-2 main protease (M^pro) interacts with peptidomimetic drugs.

•

The interaction between the M^pro and the peptidomimetics is energetically favourable.

•

Isavuconazonium, a P2–P3 α-ketoamide derivative and Pentagastrin showed the tightest and most favourable binding to M^pro.

New Cysteine Protease Inhibitors: Electrophilic (Het)arenes and Unexpected Prodrug Identification for the Trypanosoma Protease Rhodesain

Electrophilic (het)arenes can undergo reactions with nucleophiles yielding π- or Meisenheimer (σ-) complexes or the products of the S_NAr addition/elimination reactions. Such building blocks have only rarely been employed for the design of enzyme inhibitors. Herein, we demonstrate the combination of a peptidic recognition sequence with such electrophilic (het)arenes to generate highly active inhibitors of disease-relevant proteases. We further elucidate an unexpected mode of action for the trypanosomal protease rhodesain using NMR spectroscopy and mass spectrometry, enzyme kinetics and various types of simulations. After hydrolysis of an ester function in the recognition sequence of a weakly active prodrug inhibitor, the liberated carboxylic acid represents a highly potent inhibitor of rhodesain (K_i = 4.0 nM). The simulations indicate that, after the cleavage of the ester, the carboxylic acid leaves the active site and re-binds to the enzyme in an orientation that allows the formation of a very stable π-complex between the catalytic dyad (Cys-25/His-162) of rhodesain and the electrophilic aromatic moiety. The reversible inhibition mode results because the S_NAr reaction, which is found in an alkaline solvent containing a low molecular weight thiol, is hindered within the enzyme due to the presence of the positively charged imidazolium ring of His-162. Comparisons between measured and calculated NMR shifts support this interpretation.

Structural and chemical insights into the covalent-allosteric inhibition of the protein kinase Akt

Structure-based driven synthesis and biological evaluation provide innovative novel covalent-allosteric Akt inhibitors.

The Ser/Thr kinase Akt (Protein Kinase B/PKB) is a master switch in cellular signal transduction pathways. Its downstream signaling influences cell proliferation, cell growth, and apoptosis, rendering Akt a prominent drug target. The unique activation mechanism of Akt involves a change of the relative orientation of its N-terminal pleckstrin homology (PH) and the kinase domain and makes this kinase suitable for highly specific allosteric modulation. Here we present a unique set of crystal structures of covalent-allosteric interdomain inhibitors in complex with full-length Akt and report the structure-based design, synthesis, biological and pharmacological evaluation of a focused library of these innovative inhibitors.

A proton transfer network that generates deprotonated tyrosine is a key to producing reactive oxygen species in phototoxic KillerRed protein

KillerRed is the first genetically encoded photosensitizer that can induce cytotoxicity upon light exposure. Nevertheless, its phototoxicity is still lower than that of chemical photosensitizers, and the efforts to further develop KillerRed variants with enhanced phototoxicity have been impeded because the mechanism by which it generates cytotoxic reactive oxygen species (ROS) has remained elusive. To shed light on this issue, we employ quantum mechanics/molecular mechanics (QM/MM) modeling with statistical free energy analysis to examine the photo-induced electron transfer reaction occurring in KillerRed. We identify a deprotonated tyrosine residue (Tyr110) as an electron donor and further show that adjacent glutamate and serine residues play essential roles in deprotonating Tyr110. We also show that water mediation is important in the proton transfer and that protein fluctuations importantly govern the fate of the excited system. We provide clues about why KillerRed can only exhibit a low ROS yield and suggest future directions of mutagenesis toward an enhanced phototoxicity.

Modeling covalent-modifier drugs

In this review, we present a summary of how computer modeling has been used in the development of covalent-modifier drugs. Covalent-modifier drugs bind by forming a chemical bond with their target. This covalent binding can improve the selectivity of the drug for a target with complementary reactivity and result in increased binding affinities due to the strength of the covalent bond formed. In some cases, this results in irreversible inhibition of the target, but some targeted covalent inhibitor (TCI) drugs bind covalently but reversibly. Computer modeling is widely used in drug discovery, but different computational methods must be used to model covalent modifiers because of the chemical bonds formed. Structural and bioinformatic analysis has identified sites of modification that could yield selectivity for a chosen target. Docking methods, which are used to rank binding poses of large sets of inhibitors, have been augmented to support the formation of protein-ligand bonds and are now capable of predicting the binding pose of covalent modifiers accurately. The pK_a's of amino acids can be calculated in order to assess their reactivity towards electrophiles. QM/MM methods have been used to model the reaction mechanisms of covalent modification. The continued development of these tools will allow computation to aid in the development of new covalent-modifier drugs. This article is part of a Special Issue entitled: Biophysics in Canada, edited by Lewis Kay, John Baenziger, Albert Berghuis and Peter Tieleman.

Evidence for substrate binding-induced zwitterion formation in the catalytic Cys-His dyad of the SARS-CoV main protease

The coronavirus main protease (M(pro)) represents an attractive drug target for antiviral therapy of coronavirus (CoV) infections, including severe acute respiratory syndrome (SARS). The SARS-CoV M(pro) and related CoV proteases have several distinct features, such as an uncharged Cys-His catalytic dyad embedded in a chymotrypsin-like protease fold, that clearly separate these enzymes from archetypical cysteine proteases. To further characterize the catalytic system of CoV main proteases and to obtain information about improved inhibitors, we performed comprehensive simulations of the proton-transfer reactions in the SARS-CoV M(pro) active site that lead to the Cys(-)/His(+) zwitterionic state required for efficient proteolytic activity. Our simulations, comprising the free enzyme as well as substrate-enzyme and inhibitor-enzyme complexes, lead us to predict that zwitterion formation is fostered by substrate binding but not inhibitor binding. This indicates that M(pro) employs a substrate-induced catalytic mechanism that further enhances its substrate specificity. Our computational data are in line with available experimental results, such as X-ray geometries, measured pKa values, mutagenesis experiments, and the measured differences between the kinetic parameters of substrates and inhibitors. The data also provide an atomistic picture of the formerly postulated electrostatic trigger involved in SARS-CoV M(pro) activity. Finally, they provide information on how a specific microenvironment may finely tune the activity of M(pro) toward specific viral protein substrates, which is known to be required for efficient viral replication. Our simulations also indicate that the low inhibition potencies of known covalently interacting inhibitors may, at least in part, be attributed to insufficient fostering of the proton-transfer reaction. These findings suggest ways to achieve improved inhibitors.

QM/MM investigations of organic chemistry oriented questions

About 35 years after its first suggestion, QM/MM became the standard theoretical approach to investigate enzymatic structures and processes. The success is due to the ability of QM/MM to provide an accurate atomistic picture of enzymes and related processes. This picture can even be turned into a movie if nuclei-dynamics is taken into account to describe enzymatic processes. In the field of organic chemistry, QM/MM methods are used to a much lesser extent although almost all relevant processes happen in condensed matter or are influenced by complicated interactions between substrate and catalyst. There is less importance for theoretical organic chemistry since the influence of nonpolar solvents is rather weak and the effect of polar solvents can often be accurately described by continuum approaches. Catalytic processes (homogeneous and heterogeneous) can often be reduced to truncated model systems, which are so small that pure quantum-mechanical approaches can be employed. However, since QM/MM becomes more and more efficient due to the success in software and hardware developments, it is more and more used in theoretical organic chemistry to study effects which result from the molecular nature of the environment. It is shown by many examples discussed in this review that the influence can be tremendous, even for nonpolar reactions. The importance of environmental effects in theoretical spectroscopy was already known. Due to its benefits, QM/MM can be expected to experience ongoing growth for the next decade.In the present chapter we give an overview of QM/MM developments and their importance in theoretical organic chemistry, and review applications which give impressions of the possibilities and the importance of the relevant effects. Since there is already a bunch of excellent reviews dealing with QM/MM, we will discuss fundamental ingredients and developments of QM/MM very briefly with a focus on very recent progress. For the applications we follow a similar strategy.

Benchmark Study for the Cysteine-Histidine Proton Transfer Reaction in a Protein Environment: Gas Phase, COSMO, QM/MM Approaches

Proton transfer reactions are of crucial interest for the investigation of proteins. We have investigated the accuracy of commonly used quantum chemical methods for the description of proton transfer reactions in different environments (gas phase, COSMO, QM/MM) using the proton transfer between the catalytic dyad residues cysteine 145 and histidine 41 of SARS coronavirus main protease as a case study. The test includes thermodynamic, kinetic, and structural properties. The study comprises computationally demanding ab initio approaches (HF, CC2, MP2, SCS-CC2, SCS-MP2, CCSD(T)), popular density functional theories (BLYP, B3LYP, M06-2X), and semiempirical methods (MNDO/d, AM1, RM1, PM3, PM6). The approximated coupled cluster approach LCCSD(T) is taken as a reference method. We find that the robustness of the tested methods with respect to the environment correlates well with the level of theory. As an example HF, CC2, MP2, and their SCS variants show similar errors for gas phase, COSMO, or QM/MM computations. In contrast for semiempirical methods, the errors strongly diversify if one goes from gas phase to COSMO or QM/MM. Particular problems are observed for the recent semiempirical methods PM6 and RM1, which show the best performance for gas phase calculations but possess larger errors in conjunction with COSMO. Finally, a combination of SCS-MP2 and B3LYP or M06-2X allows reliable estimates about remaining errors.

Mechanism of falcipain-2 inhibition by α,β-unsaturated benzo[1,4]diazepin-2-one methyl ester

Falcipain-2 (FP-2) is a papain-family cysteine protease of Plasmodium falciparum whose primary function is to degrade the host red cell hemoglobin, within the food vacuole, in order to provide free amino acids for parasite protein synthesis. Additionally it promotes host cell rupture by cleaving the skeletal proteins of the erythrocyte membrane. Therefore, the inhibition of FP-2 represents a promising target in the search of novel anti-malarial drugs. A potent FP-2 inhibitor, characterized by the presence in its structure of the 1,4-benzodiazepine scaffold and an α,β-unsaturated methyl ester moiety capable to react with the Cys42 thiol group located in the active site of FP-2, has been recently reported in literature. In order to study in depth the inhibition mechanism triggered by this interesting compound, we carried out, through ONIOM hybrid calculations, a computational investigation of the processes occurring when the inhibitor targets the enzyme and eventually leads to an irreversible covalent Michael adduct. Each step of the reaction mechanism has been accurately characterized and a detailed description of each possible intermediate and transition state along the pathway has been reported.

Mechanistic study of the reaction of thiol-containing enzymes with alpha,beta-unsaturated carbonyl substrates by computation and chemoassays

We investigated the reactions between substituted α,β‐unsaturated carbonyl compounds (Michael systems) and thiols by computations as well as chemoassays. The results give insight into variations in the underlying mechanisms as a function of the substitution pattern. This is of interest for the mechanisms of inhibition of the SARS coronavirus main protease (SARS‐CoV M^pro) by etacrynic acid derivatives as well as for the excess toxicity of substituted α,β‐unsaturated carbonyl compounds. This study compares possible reaction courses including 1,4‐addition followed by a ketonization step, and underscores the importance of a base‐catalyzed step for the reactivity of thiol groups in enzymes. Phenyl and methyl substituents at the Michael system decrease the reactivity of the electrophilic compound, but chlorophenyl substituents partly recover the reactivity. Computations also indicate that electron‐pushing substituents lead to a change in the reaction mechanism. The conformation of the Michael system is also found to significantly influence reactivity: the s‐cis conformation leads to higher reactivity than the s‐trans conformation. The computed data explain the trends in measured inhibition potencies of substituted α,β‐unsaturated carbonyl compounds and of reaction rates in chemical assays. They also indicate that the reversibility of inhibition does not stand in contrast to the formation of a new covalent bond between inhibitor and protease.

Falcipain-2 inhibitors

Malaria, particularly that one caused by Plasmodium falciparum, remains a serious health problem in Africa, South America, and many parts of Asia where it afflicts about 500 million people and is responsible for the death of more than one million children each year. The main reasons for the persistence of malaria are the emergence of resistance to common antimalarial drugs, inadequate control of mosquito vectors, and the lack of effective vaccines. Therefore, the identification and characterization of new targets for antimalarial chemotherapy are of urgent priority. This review is focused on inhibitors of falcipain-2, a cysteine protease from P. falciparum, which represents one of the most promising targets for antimalarial drug design. Falcipain-2 is a key enzyme in the life cycle of P. falciparum since it degrades hemoglobin, at the early trophozoite stage, and cleaves ankyrin and protein 4.1, the cytoskeletal elements vital to the stability of red cell membrane, at the schizont stage. The main classes of falcipain-2 inhibitors are peptides or peptidomimetics bearing the most popular pharmacophores of cysteine protease inhibitors, such as vinyl sulfones, halomethyl ketones, and aldehydes. Furthermore, many other chemotypes have been identified as inhibitors of falcipain-2, such as isoquinolines, thiosemicarbazones, and chalcones. These inhibitors represent all classes, which, to the best of our knowledge, have been disclosed in journal articles to date.

Challenging a paradigm: theoretical calculations of the protonation state of the Cys25-His159 catalytic diad in free papain

A central mechanistic paradigm of cysteine proteases is that the His-Cys catalytic diad forms an ion-pair NH(+)/S(-) already in the catalytically active free enzyme. Most molecular modeling studies of cysteine proteases refer to this paradigm as their starting point. Nevertheless, several recent kinetics and X-ray crystallography studies of viral and bacterial cysteine proteases depart from the ion-pair mechanism, suggesting general base catalysis. We challenge the postulate of the ion-pair formation in free papain. Applying our QM/SCRF(VS) molecular modeling approach, we analyzed all protonation states of the catalytic diad in free papain and its SMe derivative, comparing the predicted and experimental pK(a) data. We conclude that the His-Cys catalytic diad in free papain is fully protonated, NH(+)/SH. The experimental pK(a) = 8.62 of His159 imidazole in free papain, obtained by NMR-controlled titration and originally interpreted as the NH(+)/S(-) <==> N/S(-) NH(+)/S(-) <==> N/S(-) equilibrium, is now assigned to the NH(+)/SH <==> N/SH NH(+)/SH <==> N/SH equilibrium.

Origin of the reactivity differences of substituted aziridines: CN vs CC bond breakages

Aziridines are broadly used as starting materials for various chemical syntheses, and the underlying reactions (CN vs CC bond breaking accompanied by an attack of a nucleophile or a dipolarophile) are strongly influenced by the substitution pattern. The present study investigates reaction courses of possible ring-opening reactions accompanied by the attack of a nucleophile for different substitution patterns of the aziridine. Information is obtained through the computation of the underlying potential energy surfaces and reaction paths. The results provide insight into the mechanisms of different ring-opening reactions and explain how the kinetics and thermodynamics of the reaction are influenced by substituents. This allows predicting substitution patterns that steer the reaction course to either CN or CC bond cleavage.

Environmental effects on charge densities of biologically active molecules: do molecule crystal environments indeed approximate protein surroundings?

In the present paper, we investigate whether crystal and enzyme environments influence the electron density (ED) of active compounds in a similar manner. This supposition is essential for high-resolution X-ray studies, which use the EDs obtained from crystals of the pure active compound as approximations for the ED of the active compound in its complex with the target enzyme. The EDs of such complexes determine the molecular recognition process between the targeted enzyme and active compound and are, hence, extremely useful tools for rational drug design. The approximation of such EDs by data obtained from crystals of the pure active compound is needed since high-resolution X-ray experiments of the target-ligand complexes are still extremely demanding. Quantum mechanical/molecular mechanical (QM/MM) and pure QM calculations are employed to determine the EDs of two inhibitors, the reversible trans-4-(aminomethyl)cyclohexane-1-carboxylic acid (AMCHA) and the irreversible E64c in four different environments (the enzyme-inhibitor complex, crystals of the pure compounds, a continuum solvation model, and the gas phase). Our investigation shows that the environment inside of the crystal of the pure active compound generally influences the ED of an active compound in a very similar way as the enzyme surrounding in the complex between the active compound and target enzyme. However, this does not hold any more if the geometrical arrangement of the inhibitor in the enzyme differs significantly from that in the crystal. While EDs computed for gas-phase environments deviate strongly from those in crystal and protein surroundings, polar solvent environments provide rather similar electron distributions. Thus, such continuum solvation models are very well suited to compute density databases which are to be employed for the determination of the ED of macromolecules.

QM/MM methods for biomolecular systems

Combined quantum-mechanics/molecular-mechanics (QM/MM) approaches have become the method of choice for modeling reactions in biomolecular systems. Quantum-mechanical (QM) methods are required for describing chemical reactions and other electronic processes, such as charge transfer or electronic excitation. However, QM methods are restricted to systems of up to a few hundred atoms. However, the size and conformational complexity of biopolymers calls for methods capable of treating up to several 100,000 atoms and allowing for simulations over time scales of tens of nanoseconds. This is achieved by highly efficient, force-field-based molecular mechanics (MM) methods. Thus to model large biomolecules the logical approach is to combine the two techniques and to use a QM method for the chemically active region (e.g., substrates and co-factors in an enzymatic reaction) and an MM treatment for the surroundings (e.g., protein and solvent). The resulting schemes are commonly referred to as combined or hybrid QM/MM methods. They enable the modeling of reactive biomolecular systems at a reasonable computational effort while providing the necessary accuracy.

Atomistic insights into the inhibition of cysteine proteases: first QM/MM calculations clarifying the stereoselectivity of epoxide-based inhibitors

Due to their important role in many diseases, cysteine proteases represent new promising drug targets. An important class of cysteine-protease inhibitors is derived from the naturally occurring compound E64, possessing an epoxysuccinyl moiety as warhead. Experimental studies show stereoselectivity concerning the inhibition potency, e.g., a trans-configured epoxide ring is essential for inhibition, and furthermore, in most cases, the ( S, S)-configured inhibitors have a higher inhibition potency than their ( R, R)-counterparts. However, the underlying effects are not fully understood. In this work, such effects are investigated by classical molecular dynamics simulations and combined quantum mechanics/molecular modeling (QM/MM) calculations for the E64c-cathepsin B complex. Our computations reveal that the hydrogen bonding network between the enzyme and the E64c (or its derivatives) determines the stereoselectivity of the subsequent ring opening reaction by governing the distance between the attacking thiolate and the attacked C2 atom of the epoxide ring. For the ( S, S)-configuration, a strong network can be realized which enables a close contact between the reacting centers, so that the irreversible step becomes very efficient. The ( R, S)-configuration ( cis-configuration) can only form networks in which the two reacting centers are so far away from each other that the irreversible step can hardly happen. The ( R, R)-configuration is in between, less optimal than the ( S, S)-configuration but much better than the ( R, S)-configuration. Exceptions where the ( R, R)-configurations shows higher potency than the ( S, S) ones are also explained.

On the origin of the stabilization of the zwitterionic resting state of cysteine proteases: a theoretical study

Papain-like cysteine proteases are ubiquitous proteolytic enzymes. The protonated His199/deprotonated Cys29 ion pair (cathepsin B numbering) in the active site is essential for their proper functioning. The presence of this ion pair stands in contrast to the corresponding intrinsic residue p K a values, indicating a strong influence of the enzyme environment. In the present work we show by molecular dynamics simulations on quantum mechanical/molecular mechanical (QM/MM) potentials that the ion pair is stabilized by a complex hydrogen bond network which comprises several amino acids situated in the active site of the enzyme and 2-4 water molecules. QM/MM reaction path computations for the proton transfer from His199 to the thiolate of the Cys29 moiety indicate that the ion pair is about 32-36 kJ mol (-1) more stable than the neutral form if the whole hydrogen bonding network is active. Without any hydrogen bonding network the ion pair is predicted to be significantly less stable than the neutral form. QM/MM charge deletion analysis and QM model calculations are used to quantify the stabilizing effect of the active-site residues and the L1 helix in favor of the zwitterionic form. The active-site water molecules contribute about 30 kJ mol (-1) to the overall stabilization. Disruption of the hydrogen bonding network upon substrate binding is expected to enhance the nucleophilic reactivity of the thiolate.

Atomistic insights into the inhibition of cysteine proteases: first QM/MM calculations clarifying the regiospecificity and the inhibition potency of epoxide- and aziridine-based inhibitors

Epoxides and aziridines are important building blocks for inhibitors of cysteine proteases which are promising drug targets for many diseases. In spite of the large amount of experimental data concerning inhibition potency, structure-activity relationships, and structural arrangements of enzyme-inhibitor complexes, little is known about the basic principles which connect the substitution pattern with the resulting activities. To shed some light on this issue which is essential for the rational design of improved compounds, we have studied the inhibition processes theoretically for various inhibitors using quantum mechanical/molecular mechanical hybrid approaches and classical molecular dynamics simulations. The careful analysis of the computational results allows insight into the interactions which govern the regio- and stereospecificity of the interactions. Known structure-activity relationships are rationalized in terms of the same interactions that determine the measured pH dependencies. Inconsistencies in existing X-ray structures are resolved through comparison with the computed structures, which leads to a reassessment of the factors that control the inhibition potency. Similarities and differences in the mode of action of epoxide- and aziridine-based inhibitors are elucidated. Finally the small reaction barriers computed for the irreversible step in E64 analogues call into question the commonly accepted two-step model of inhibition since the second, irreversible step is predicted to be so fast that suitably oriented enzyme-inhibitor complexes will react rather than dissociate and equilibrate.

Conformational analysis of arginine in gas phase—A strategy for scanning the potential energy surface effectively

The determination of all possible low-lying energy conformers of flexible molecules is of fundamental interest for various applications. It necessitates a reliable conformational search that is able to detect all important minimum structures and calculates the energies on an adequate level of theory. This work presents a strategy to identify low-energy conformers using arginine as an example by means of a force-field based conformational search in combination with high-level geometry optimizations (RI-MP2/TZVPP+). The methods used for various stages in the conformational search strategy are shown and various pitfalls are discussed. We can show that electronic energies calculated on a DFT level of theory with standard exchange-correlation functionals strongly underestimate the intramolecular stabilization resulting from stacked orientations of the guanidine and carbonyl moiety of arginine due to the deficiency of DFT to describe dispersion effects. In this case by usage of electron correlation methods, low energy conformers comprising stacked arrangements that are counterintuitive become favorable.