Computational analysis of plasmepsin IV bound to an allophenylnorstatine inhibitor

Molecular modelling approaches can reveal the Molecular interactions established between antimalarial targets of hemozoin pathway and the organic phytochemicals of Artocarpus species

Ayurveda, the traditional Indian medical system, has potential applications in early malaria treatment. In an in silico docking study, 50 phytochemicals from two plants Artocarpus lakoocha Roxb. (AL) And Artocarpus heterophyllus Lam. (AH), were examined for their interactions with anti-malarial proteins (PDB IDs: 3BWK, 3BPF, 1LF3). The nucleotide analogue Artemisinin, a current malaria treatment, served as a positive control. Result showed that phytochemicals from AL and AH exhibited binding affinities as high as -9.6 kcal/mol, respectively. Additionally, molecular dynamics simulation for Artocarpin: 3BPF demonstrated stable complexes over 100 ns. Notably, Artocarpin and Quercetin displayed higher binding affinities (up to -9.6 as well as -9.5 kcal/mol, respectively) compared to Artemisinin (-7.5 up to kcal/mol), have shown. Pharmacokinetic predictions indicated the compounds were likely non-carcinogenic, water-soluble and biologically safe. In-vitro analysis using β-Hematin assay supported these findings, suggesting the phytochemicals as Hemozoin pathway inhibitors.

Computational insights for predicting the binding and selectivity of peptidomimetic plasmepsin IV inhibitors against cathepsin D

Plasmepsins (Plms) are aspartic proteases involved in the degradation of human hemoglobin by P. falciparum and are essential for the survival and growth of the parasite. Therefore, Plm enzymes are reported as an important antimalarial drug target. Herein, we have applied molecular docking, molecular dynamics (MD) simulations, and binding free energy with the Linear Interaction Energy (LIE) approach to investigate the binding of peptidomimetic PlmIV inhibitors with a particular focus on understanding their selectivity against the human Asp protease cathepsin D (CatD). The residual decomposition analysis results suggest that amino acid differences in the subsite S3 of PlmIV and CatD are responsible for the higher selectivity of the 5a inhibitor. These findings yield excellent agreement with experimental binding data and provide new details regarding van der Waals and electrostatic interactions of subsite residues as well as structural properties of the PlmIV and CatD systems.

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.

Combined kinetic studies and computational analysis on kojic acid analogous as tyrosinase inhibitors

Tyrosinase is a key enzyme in melanin synthesis and widely distributed in plants and animals tissues. In mammals, this enzyme is related to pigment production, involved in wound healing, primary immune response and it can also contribute to catecholamines synthesis in the brain. Consequently, tyrosinase enzyme represents an attractive and selective target in the field of the medicine, cosmetics and bio-insecticides. In this paper, experimental kinetics and computational analysis were used to study the inhibition of tyrosinase by analogous of Kojic acid. The main interactions occurring between inhibitors-tyrosinase complexes and the influence of divalent cation (Cu2+) in enzymatic inhibition were investigated by using molecular docking, molecular dynamic simulations and electrostatic binding free energy by using the Linear Interaction Energy (LIE) method. The results showed that the electrostatic binding free energy are correlated with values of constant inhibition (r2 = 0.97).Thus, the model obtained here could contribute to future studies of this important system and, therefore, eventually facilitate development of tyrosinase inhibitors.

Computational analysis of human OGA structure in complex with PUGNAc and NAG-thiazoline derivatives

The substitution of serine and threonine residues in nucleocytoplasmic proteins with 2-acetamido-2-deoxy-β-D-glucopyranose (O-GlcNAc) residues is an essential post-translational modification found in many multicellular eukaryotes. O-glycoprotein 2-acetamino-2-deoxy-β-D-glucopyranosidase (O-GlcNAcase) hydrolyzes O-GlcNAc residues from post-translationally modified serine/threonine residues of nucleocytoplasmic protein. O-GlcNAc has been implicated in several disease states such as cancer, Alzheimer's disease, and type II diabetes. For this paper, a model of the human O-GlcNAcase (hOGA) enzyme based on the X-ray structures of bacterial Clostridium perfringens (CpNagJ) and Bacteroides thetaiotaomicrometer (BtOGA) homologues has been generated through molecular homology modeling. In addition, molecular docking, molecular dynamics (MD) simulations, and Linear Interaction Energy (LIE) were employed to determine the bind for derivatives of two potent inhibitors: O-(2-acetamido-2-deoxy-D-glucopyranosylidene)amino-N-phenylcarbamate (PUGNAc) and 1,2-dideoxy-2'-methyl-R-D-glucopyranoso-[2,1-d]-Δ2'-thiazoline (NAG-thiazoline), with hOGA. The results show that the binding free energy calculations using the Linear Interaction Energy (LIE) are correlated with inhibition constant values. Therefore, the model of the human O-GlcNAcase (hOGA) obtained here may be used as a target for rational design of new inhibitors.

Linear interaction energy: method and applications in drug design

A broad range of computational methods exist for the estimation of ligand-protein binding affinities. In this chapter we will provide a guide to the linear interaction energy (LIE) method for binding free energy calculations, focusing on the drug design problem. The method is implemented in combination with molecular dynamics (MD) sampling of relevant conformations of the ligands and complexes under consideration. The detailed procedure for MD sampling is followed by key notes in order to properly analyze such sampling and obtain sufficiently accurate estimations of ligand-binding affinities.

Computational perspectives into plasmepsins structure-function relationship: implications to inhibitors design

The development of efficient and selective antimalariais remains a challenge for the pharmaceutical industry. The aspartic proteases plasmepsins, whose inhibition leads to parasite death, are classified as targets for the design of potent drugs. Combinatorial synthesis is currently being used to generate inhibitor libraries for these enzymes, and together with computational methodologies have been demonstrated capable for the selection of lead compounds. The high structural flexibility of plasmepsins, revealed by their X-ray structures and molecular dynamics simulations, made even more complicated the prediction of putative binding modes, and therefore, the use of common computational tools, like docking and free-energy calculations. In this review, we revised the computational strategies utilized so far, for the structure-function relationship studies concerning the plasmepsin family, with special focus on the recent advances in the improvement of the linear interaction estimation (LIE) method, which is one of the most successful methodologies in the evaluation of plasmepsin-inhibitor binding affinity.

Crystal structures of the free and inhibited forms of plasmepsin I (PMI) from Plasmodium falciparum

Plasmepsin I (PMI) is one of the four vacuolar pepsin-like proteases responsible for hemoglobin degradation by the malarial parasite Plasmodium falciparum, and the only one with no crystal structure reported to date. Due to substantial functional redundancy of these enzymes, lack of inhibition of even a single plasmepsin can defeat efforts in creating effective antiparasitic agents. We have now solved crystal structures of the recombinant PMI as apoenzyme and in complex with the potent peptidic inhibitor, KNI-10006, at the resolution of 2.4 and 3.1Å, respectively. The apoenzyme crystallized in the orthorhombic space group P2(1)2(1)2(1) with two molecules in the asymmetric unit and the structure has been refined to the final R-factor of 20.7%. The KNI-10006 bound enzyme crystallized in the tetragonal space group P4(3) with four molecules in the asymmetric unit and the structure has been refined to the final R-factor of 21.1%. In the PMI-KNI-10006 complex, the inhibitors were bound identically to all four enzyme molecules, with the opposite directionality of the main chain of KNI-10006 relative to the direction of the enzyme substrates. Such a mode of binding of inhibitors containing an allophenylnorstatine-dimethylthioproline insert in the P1-P1' positions, previously reported in a complex with PMIV, demonstrates the importance of satisfying the requirements for the proper positioning of the functional groups in the mechanism-based inhibitors towards the catalytic machinery of aspartic proteases, as opposed to binding driven solely by the specificity of the individual enzymes. A comparison of the structure of the PMI-KNI-10006 complex with the structures of other vacuolar plasmepsins identified the important differences between them and may help in the design of specific inhibitors targeting the individual enzymes.

Computational analysis of aspartic protease plasmepsin II complexed with EH58 inhibitor: a QM/MM MD study

Plasmepsin (PM) II is one of four enzymes in the food vacuole of Plasmodium falciparum. It has become an attractive target for combating malaria through research regarding its importance in the P. falciparum metabolism and life cycle, making it the target of choice for structure-based drug design. This paper reports the results of hybrid quantum mechanics / molecular mechanics (QM/MM) molecular dynamics (MD) simulations employed to study the details of the interactions established between PM II and N-(3-{(2-benzo[1, 3]dioxol-5-yl-ethyl)[3-(1-methyl-3-oxo-1,3-dihydro-isoindol-2-yl) propionyl]-amino}-1-benzyl-2-(hydroxyl-propyl)-4-benzyloxy-3,5dimethoxy-benzamide (EH58), a well-known potent inhibitor for this enzyme. Electrostatic binding free energy and energy terms decomposition have been computed for PM II complexed with the EH58 inhibitor. The results reveal that there is a strong interaction between Asp34, Val78, Ser79, Tyr192 and Asp214 residues and the EH58 inhibitor. In addition, we have computed the potential of the mean force (PMF) profile in order to assign the protonation state of the two catalytic aspartates in PM II-EH58 complex. The results indicate that the protonation of Asp214 favors a stable active site structure, which is consistent with our electrostatic binding free energy calculation and with previous published works.

Insight into estrogenicity of phytoestrogens using in silico simulation

Phytoestrogens, including miroestrol and deoxymiroestrol, have the ability to act through competition with estrogen for binding to the estrogen receptor (ER). Here, we utilize manual ligand docking followed by molecular dynamics simulations and binding free energy calculations with the linear interaction energy method to predict the binding modes and the binding affinities of phytoestrogens on the ligand binding domain of ER (ERalpha-LBD). The calculations brought about the good correlation between the calculated binding free energy and the bioassays. Furthermore, consideration of Lennard-Jones and Coulomb interaction energies of miroestrol and deoxymiroestrol on ERalpha-LBD provided the information to develop the phytoestrogen derivatives as the preferred drug for ER positive breast cancer treatment.

Antimalarial activity enhancement in hydroxymethylcarbonyl (HMC) isostere-based dipeptidomimetics targeting malarial aspartic protease plasmepsin

Plasmepsin (Plm) is a potential target for new antimalarial drugs, but most reported Plm inhibitors have relatively low antimalarial activities. We synthesized a series of dipeptide-type HIV protease inhibitors, which contain an allophenylnorstatine-dimethylthioproline scaffold to exhibit potent inhibitory activities against Plm II. Their activities against Plasmodium falciparum in the infected erythrocyte assay were largely different from those against the target enzyme. To improve the antimalarial activity of peptidomimetic Plm inhibitors, we attached substituents on a structure of the highly potent Plm inhibitor KNI-10006. Among the derivatives, we identified alkylamino compounds such as 44 (KNI-10283) and 47 (KNI-10538) with more than 15-fold enhanced antimalarial activity, to the sub-micromolar level, maintaining their potent Plm II inhibitory activity and low cytotoxicity. These results suggest that auxiliary substituents on a specific basic group contribute to deliver the inhibitors to the target Plm.

A Model for How Ribosomal Release Factors Induce Peptidyl-tRNA Cleavage in Termination of Protein Synthesis

A major unresolved question in messenger RNA translation is how ribosomal release factors terminate protein synthesis. Class 1 release factors decode stop codons and trigger hydrolysis of the bond between the nascent polypeptide and tRNA some 75 A away from the decoding site. While the gross features of the release factor-ribosome interaction have been revealed by low-resolution crystal structures, there is no information on the atomic level at either the decoding or peptidyl transfer center. We used extensive computer simulations, constrained by experimental data, to predict how bacterial release factors induce peptide dissociation from the ribosome. A distinct structural solution is presented for how the methylated Gln residue of the universally conserved GGQ release factor motif inserts into the ribosomal A site and promotes rapid reaction with the peptidyl-tRNA substrate. This model explains key mutation experiments and shows that the ribosomal peptidyl transfer center catalyzes its two chemical reactions by a common mechanism.

The protonation state of the catalytic aspartates in plasmepsin II

Assigning the correct protonation state to the catalytic residues is essential for a realistic modelling of an enzyme's active site. Plasmepsins are pharmaceutically relevant aspartic proteases involved in haemoglobin degradation by Plasmodium spp. In aspartic proteases, one of the two catalytic aspartates is protonated, while the other is negatively charged. Here, multiple explicit-water molecular dynamics simulations of plasmepsin II, uncomplexed and with a hydroxypropylamine peptidomimetic inhibitor, indicate that protonation of Asp214 favours a stable active site structure. Moreover, the protonation state of the catalytic aspartate has a strong influence on a linear chain of hydrogen bonds with the adjacent side chains.

Additional interaction of allophenylnorstatine-containing tripeptidomimetics with malarial aspartic protease plasmepsin II

Based on a highly potent allophenylnorstatine-containing inhibitor, KNI-10006, against the plasmepsins of Plasmodium falciparum, we synthesized a series of tripeptide-type compounds with various N-terminal moieties and evaluated their inhibitory activities against plasmepsin II. Certain phenylacetyl derivatives exhibited extremely high affinity with K(i) values of less than 0.1n M suggesting successful hydrophobic interactions.