Trifluoroethanol-Induced Changes in Activity and Conformation of Manganese-Containing Superoxide Dismutase

Superoxide dismutase (SOD, EC 1.15.1.1) plays an important role in antioxidant defense in organisms exposed to oxygen. However, there is a lack of research into the regulation of SOD activity and structural changes during folding, especially for SOD originating from extremophiles. We studied the inhibitory effects of trifluoroethanol (TFE) on the activity and conformation of manganese-containing SOD (Mn-SOD) from Thermus thermophilus. TFE decreased the degree of secondary structure of Mn-SOD, which directly resulted in enzyme inactivation and disrupted the tertiary structure of Mn-SOD. The kinetic studies showed that TFE-induced inactivation of Mn-SOD is a first-order reaction and that the regional Mn-contained active site is very stable compared to the overall structure. We further simulated the docking between Mn-SOD and TFE (binding energy for Dock 6.3, −9.68 kcal/mol) and predicted that the LEU9, TYR13, and HIS29 residues outside of the active site interact with TFE. Our results provide insight into the inactivation of Mn-SOD during unfolding in the presence of TFE and allow us to describe ligand binding via inhibition kinetics combined with computational predictions.     

Discovery of Novel Acetaldehyde Dehydrogenase 1A1(ALDH1A1) Inhibitors by Utilizing 3D-QSAR,Molecular Docking and Molecular Dynamics Simulation

Acetaldehyde dehydrogenase 1A1 is a hopeful therapeutic target to ovarian cancer. In this present work, 3D-QSAR, molecular docking and molecular dynamics(MD) simulations were implemented on a series of quinoline-based ALDH1A1 inhibitors to investigate novel acetaldehyde dehydrogenase 1A1 inhibitors as anticancer adjuvant drugs for ovarian cancer. Two reliable CoMFA(Q~2 = 0.583, R~2 = 0.967) and CoMSIA(Q~2 = 0.640, R~2 = 0.977) models of ALDH1A1 inhibitors were established. Novel ALDH1A1 inhibitors were predicted by the 3D-QSAR models. Molecular docking reveals important residues for protein-compound interactions, and the results revealed ALDH1A1 inhibitors had stronger electrostatic interaction and binding affinity with key residues of protein, such as Phe171, Val174 and Cys303. Molecular dynamics simulations further verified the results of molecular docking. The above information provided significant guidance for the design of novel ALDH1A1 inhibitors.     

Alpha-Glucosidase Folding During Urea Denaturation: Enzyme Kinetics and Computational Prediction

In this study, we investigated structural changes in alpha-glucosidase during urea denaturation. Alpha-glucosidase was inactivated by urea in a dose-dependent manner. The inactivation was a first-order reaction with a monophase process. Urea inhibited alpha-glucosidase in a mixed-type reaction. We found that an increase in the hydrophobic surface of this enzyme induced by urea resulted in aggregation caused by unstable folding intermediates. We also simulated the docking between alpha-glucosidase and urea. The docking simulation suggested that several residues, namely THR9, TRP14, LYS15, THR287, ALA289, ASP338, SER339, and TRP340, interact with urea. Our study provides insights into the alpha-glucosidase unfolding pathway and 3D structure of alpha-glucosidase.     

分子对接和动力学模拟提高嗜热蛋白酶PhpI的活力

通过分子对接和动力学模拟对嗜热蛋白酶的分子进行改造, 确定蛋白酶PH1704(PhpI)定点突变残基, 并通过分子生物学实验进行验证. 突变体K43C的蛋白酶活力提高了5.8倍. 分子动力学模拟结果表明, 经过8 ns的动力学模拟后, K43C突变体二级结构由野生型的S2片层(F11-E12-D13)变成环状结构. E12和K43均是活性位点的重要残基, 这种变化将导致活性位点的柔性增强, 有利于催化反应的发生.     

Structural and Dynamic Basis of Human Cytochrome P450 7B1: A Survey of Substrate Selectivity and Major Active Site Access Channels

Cytochrome P450 (CYP) 7B1 is a steroid cytochrome P450 7α‐hydroxylase that has been linked directly with bile salt synthesis and hereditary spastic paraplegia type 5 (SPG5). The enzyme provides the primary metabolic route for neurosteroids dehydroepiandrosterone (DHEA), cholesterol derivatives 25‐hydroxycholesterol (25‐HOChol), and other steroids such as 5α‐androstane‐3β,17β‐diol (anediol), and 5α‐androstene‐3β,17β‐diol (enediol). A series of investigations including homology modeling, molecular dynamics (MD), and automatic docking, combined with the results of previous experimental site‐directed mutagenesis studies and access channels analysis, have identified the structural features relevant to the substrate selectivity of CYP7B1. The results clearly identify the dominant access channels and critical residues responsible for ligand binding. Both binding free energy analysis and total interaction energy analysis are consistent with the experimental conclusion that 25‐HOChol is the best substrate. According to 20 ns MD simulations, the Phe cluster residues that lie above the active site, particularly Phe489, are proposed to merge the active site with the adjacent channel to the surface and accommodate substrate binding in a reasonable orientation. The investigation of CYP7B1–substrate binding modes provides detailed insights into the poorly understood structural features of human CYP7B1 at the atomic level, and will be valuable information for drug development and protein engineering.     

Structural,Biochemical, and Computational Studies Reveal the Mechanism of Selective Aldehyde Dehydrogenase 1A1 Inhibition by Cytotoxic Duocarmycin Analogues

Analogues of the natural product duocarmycin bearing an indole moiety were shown to bind aldehyde dehydrogenase 1A1 (ALDH1A1) in addition to DNA, while derivatives without the indole solely addressed the ALDH1A1 protein. The molecular mechanism of selective ALDH1A1 inhibition by duocarmycin analogues was unraveled through cocrystallization, mutational studies, and molecular dynamics simulations. The structure of the complex shows the compound embedded in a hydrophobic pocket, where it is stabilized by several crucial π‐stacking and van der Waals interactions. This binding mode positions the cyclopropyl electrophile for nucleophilic attack by the noncatalytic residue Cys302, thereby resulting in covalent attachment, steric occlusion of the active site, and inhibition of catalysis. The selectivity of duocarmycin analogues for ALDH1A1 is unique, since only minor alterations in the sequence of closely related protein isoforms restrict compound accessibility.     

Structural,Biochemical, and Computational Studies Reveal the Mechanism of Selective Aldehyde Dehydrogenase 1A1 Inhibition by Cytotoxic Duocarmycin Analogues

Analogues of the natural product duocarmycin bearing an indole moiety were shown to bind aldehyde dehydrogenase 1A1 (ALDH1A1) in addition to DNA, while derivatives without the indole solely addressed the ALDH1A1 protein. The molecular mechanism of selective ALDH1A1 inhibition by duocarmycin analogues was unraveled through cocrystallization, mutational studies, and molecular dynamics simulations. The structure of the complex shows the compound embedded in a hydrophobic pocket, where it is stabilized by several crucial π‐stacking and van der Waals interactions. This binding mode positions the cyclopropyl electrophile for nucleophilic attack by the noncatalytic residue Cys302, thereby resulting in covalent attachment, steric occlusion of the active site, and inhibition of catalysis. The selectivity of duocarmycin analogues for ALDH1A1 is unique, since only minor alterations in the sequence of closely related protein isoforms restrict compound accessibility.     

Chemical- and Thermal-Induced Unfolding of Leishmania donovani Ribose-5-Phosphate Isomerase B: a Single-Tryptophan Protein

Ribose-5-phosphate isomerase B (RpiB), a crucial enzyme of pentose phosphate pathway, was proposed to be a potential drug target for visceral leishmaniasis. In this study, we have analyzed the biophysical properties of Leishmania donovani RpiB (LdRpiB) enzyme to gain insight into its unfolding pathway under various chemical and thermal denaturation conditions by using fluorescence and CD spectroscopy. LdRpiB inactivation precedes the structural transition at lower concentrations of both urea and guanidine hydrochloride (GdHCl). 8-Anilinonapthalene 1-sulfonic (ANS) binding experiments revealed the presence of molten globule intermediate at 1.5 M GdHCl and a nonnative intermediate state at 6-M urea concentration. Acrylamide quenching experiments further validated the above findings, as solvent accessibility of tryptophan residues increased with increase in GdHCl and urea concentration. The recombinant LdRpiB was completely unfolded at 6 M GdHCl, whereas the enzyme molecule was resistant to complete unfolding even at 8-M urea concentration. The GdHCl- and urea-mediated unfolding involves a three-state transition process. Thermal-induced denaturation revealed complete loss of enzyme activity at 65 °C with only 20 % secondary structure loss. The formation of the well-ordered β-sheet structures of amyloid fibrils was observed after 55 °C which increased linearly till 85 °C as detected by thioflavin T dye. This study depicts the stability of the enzyme in the presence of chemical and thermal denaturants and stability-activity relationship of the enzyme. The presence of the intermediate states may have major implications in the way the enzyme binds to its natural ligand under various conditions. Also, the present study provides insights into the properties of intermediate entities of this important enzyme.     

Allostery Inhibition of BACE1 by Psychotic and Meroterpenoid Drugs in Alzheimer’s Disease Therapy

In over a century since its discovery, Alzheimer’s disease (AD) has continued to be a global health concern due to its incurable nature and overwhelming increase among older people. In this paper, we give an overview of the efforts of researchers towards identifying potent BACE1 exosite-binding antibodies and allosteric inhibitors. Herein, we apply computer-aided drug design (CADD) methods to unravel the interactions of some proposed psychotic and meroterpenoid BACE1 allosteric site inhibitors. This study is aimed at validating the allosteric potentials of these selected compounds targeted at BACE1 inhibition. Molecular docking, molecular dynamic (MD) simulations, and post-MD analyses are carried out on these selected compounds, which have been experimentally proven to exhibit allosteric inhibition on BACE1. The SwissDock software enabled us to identify more than five druggable pockets on the BACE1 structural surface using docking. Besides the active site region, a melatonin derivative (compound 1) previously proposed as a BACE1 allostery inhibitor showed appreciable stability at eight different subsites on BACE1. Refinement with molecular dynamic (MD) simulations shows that the identified non-catalytic sites are potential allostery sites for compound 1. The allostery and binding mechanism of the selected potent inhibitors show that the smaller the molecule, the easier the attachment to several enzyme regions. This finding hereby establishes that most of these selected compounds failed to exhibit strong allosteric binding with BACE1 except for compound 1. We hereby suggest that further studies and additional identification/validation of other BACE1 allosteric compounds be done. Furthermore, this additional allosteric site investigation will help in reducing the associated challenges with designing BACE1 inhibitors while exploring the opportunities in the design of allosteric BACE1 inhibitors.     

HIV蛋白酶抑制剂——利迪链菌素的分子对接研究

用分子对接的方法, 对利迪链菌素的抗HIV蛋白酶活性进行了研究. 为了更准确地反映利迪链菌素分子与酶蛋白结合的情况, 充分考虑受体活性部位的柔性, 采用了FlexX(初步对接)和Flexidock(精确对接)分两步将配体与受体进行对接. 在初步对接中, 设计了不同的受体活性部位来考察是否有结合水分子参与抑制剂与酶的结合. 对一种作用方式已知的非肽类HIV蛋白酶抑制剂Aha006进行的对接研究显示, 分子模拟的结果与实际情况吻合得较好, 证明了本文所采用的方法的可靠性. 利迪链菌素与蛋白酶活性部位的对接结果显示, 配体分子与受体之间的结合没有结合水分子的参与, 两者通过5对氢键作用结合成为稳定的复合物. 利迪链菌素占据结合腔, 覆盖了蛋白酶的活性三联体Asp25-Thr26-Gly27, 从而起到抑制其生物活性的作用.     

Influence of protonation on substrate and inhibitor interactions at the active site of human monoamine oxidase-A

Although substrate conversion mediated by human monoaminooxidase (hMAO) has been associated with the deprotonated state of their amine moiety, data regarding the influence of protonation on substrate binding at the active site are scarce. Thus, in order to assess protonation influence, steered molecular dynamics (SMD) runs were carried out. These simulations revealed that the protonated form of the substrate serotonin (5-HT) exhibited stronger interactions at the protein surface compared to the neutral form. The latter displayed stronger interactions in the active site cavity. These observations support the possible role of the deprotonated form in substrate conversion. Multigrid docking studies carried out to rationalize the role of 5-HT protonation in other sites besides the active site indicated two energetically favored docking sites for the protonated form of 5-HT on the enzyme surface. These sites seem to be interconnected with the substrate/inhibitor cavity, as revealed by the tunnels observed by means of CAVER program. pK(a) calculations in the surface loci pointed to Glu32?, Asp32?, His???, and Asp132 as candidates for a possible in situ deprotonation step. Docking analysis of a group of inhibitors (structurally related to substrates) showed further interactions with the same two docking access sites. Interestingly, the protonated/deprotonated amine moiety of almost all compounds attained different docking poses in the active site, none of them oriented to the flavin moiety, thus producing a more variable and less productive orientations to act as substrates. Our results highlight the role of deprotonation in facilitating substrate conversion and also might reflect the necessity of inhibitor molecules to adopt specific orientations to achieve enzyme inhibition.     

Link between allosteric signal transduction and functional dynamics in a multisubunit enzyme: S-adenosylhomocysteine hydrolase

S-adenosylhomocysteine hydrolase (SAHH), a cellular enzyme that plays a key role in methylation reactions including those required for maturation of viral mRNA, is an important drug target in the discovery of antiviral agents. While targeting the active site is a straightforward strategy of enzyme inhibition, evidence of allosteric modulation of active site in many enzymes underscores the molecular origin of signal transduction. Information of co-evolving sequences in SAHH family and the key residues for functional dynamics that can be identified using native topology of the enzyme provide glimpses into how the allosteric signaling network, dispersed over the molecular structure, coordinates intra- and intersubunit conformational dynamics. To study the link between the allosteric communication and functional dynamics of SAHHs, we performed Brownian dynamics simulations by building a coarse-grained model based on the holo and ligand-bound structures. The simulations of ligand-induced transition revealed that the signal of intrasubunit closure dynamics is transmitted to form intersubunit contacts, which in turn invoke a precise alignment of active site, followed by the dimer-dimer rotation that compacts the whole tetrameric structure. Further analyses of SAHH dynamics associated with ligand binding provided evidence of both induced fit and population shift mechanisms and also showed that the transition-state ensemble is akin to the ligand-bound state. Besides the formation of enzyme-ligand contacts at the active site, the allosteric couplings from the residues distal to the active site are vital to the enzymatic function.     

Aspects of the mechanism of catalysis of glucose oxidase: A docking, molecular mechanics and quantum chemical study

The complex structure of glucose oxidase (GOX) with the substrate glucose was determined using a docking algorithm and subsequent molecular dynamics simulations. Semiempirical quantum chemical calculations were used to investigate the role of the enzyme and FAD co-enzyme in the catalytic oxidation of glucose. On the basis of a small active site model, substrate binding residues were determined and heats of formation were computed for the enzyme substrate complex and different potential products of the reductive half reaction. The influence of the protein environment on the active site model was estimated with a point charge model using a mixed QM/MM method. Solvent effects were estimated with a continuum model. Possible modes of action are presented in relation to experimental data and discussed with respect to related enzymes. The calculations indicate that the redox reaction of GOX differs from the corresponding reaction of free flavins as a consequence of the protein environment. One of the active site histidines is involved in substrate binding and stabilization of potential intermediates, whereas the second histidine is a proton acceptor. The former one, being conserved in a series of oxidoreductases, is also involved in the stabilization of a C4a-hydroperoxy dihydroflavin in the course of the oxidative half reaction.     

UV radiation effects on flavocytochrome b2 in dilute aqueous solution

A dilute aqueous solution of flavocytochrome b2 when exposed to inactivating doses of UV radiation at 280 nm underwent reversible loss in activity both under aerated and deaerated conditions. The active site as well as the substrate binding sites were found to be modified in the irradiated enzyme. Irradiation of the enzyme in the UV-C range resulted in partial unfolding of the polypeptide framework. Destruction and/or modification of both tryptophan and tyrosine residues as well as heme moieties took place. Preliminary laser flash photolysis studies suggest that the initial photo-ionization takes place with tryptophan and tyrosine residues with the formation of excited states and radicals, and then rapid transfer of electrons takes place to histidyl and cystinyl sites which might have eventually been altered in the process.     

Computational Simulations Highlight the IL2Rα Binding Potential of Polyphenol Stilbenes from Fenugreek

Widely used in global households, fenugreek is well known for its culinary and medicinal uses. The various reported medicinal properties of fenugreek are by virtue of the different natural phytochemicals present in it. Regarded as a promising target, interleukin 2 receptor subunit alpha (IL2Rα) has been shown to influence immune responses. In the present research, using in silico techniques, we have demonstrated the potential IL2Rα binding properties of three polyphenol stilbenes (desoxyrhaponticin, rhaponticin, rhapontigenin) from fenugreek. As the first step, molecular docking was performed to assess the binding potential of the fenugreek phytochemicals with IL2Rα. All three phytochemicals demonstrated interactions with active site residues. To confirm the reliability of our molecular docking results, 100 ns molecular dynamics simulations studies were undertaken. As discerned by the RMSD and RMSF analyses, IL2Rα in complex with the desoxyrhaponticin, rhaponticin, and rhapontigenin indicated stability. The RMSD analysis of the phytochemicals alone also demonstrated no significant structural changes. Based on the stable molecular interactions and comparatively slightly better MM/PBSA binding free energy, rhaponticin seems promising. Additionally, ADMET analysis performed for the stilbenes indicated that all of them obey the ADMET rules. Our computational study thus supports further in vitro IL2Rα binding studies on these stilbenes, especially rhaponticin.     

Molecular modeling studies on the urease active site and the enzyme-catalyzed urea hydrolysis

Summary These studies are an attempt to gain better insight into the pharmacophore requirements of urease. On the basis of published information on this enzyme (EXAFS, amino acid sequence, essential groups at the active site) a hypothetical nickel-tripeptide complex, as preliminary substitute for the urease active site was modeled using computer-aided molecular modeling techniques. The results suggest two alternative docking modes of urea and reaction intermediates, corresponding to two different reaction mechanisms. Both binding modes are compatible with the docking of known potent inhibitors such as selected hydroxamic acids and phosphorodiamides. The results can be used to help in the design of new potential inhibitors of urease.     

Picosecond dynamics of ligand interconversion in the primary docking site of heme proteins

We have used femtosecond IR spectroscopy to probe interconversion dynamics of ligand in the primary docking site of heme proteins under physiological conditions. The docking site, fashioned with highly conserved amino acid residues, modulates ligand-binding activity by mediating the passage of ligand to and from the active binding site. Ligands in two states of the docking site interconvert on the picosecond time scale, and the rates are about 4 times slower in hemoglobin than that in myoglobin. The accurate interconversion rates on the time scale readily accessible by MD simulations can be used to refine computer simulations, which could in turn provide a detailed mechanistic picture of ligand binding in heme proteins.     

Molecular docking of cathepsin L inhibitors in the binding site of papain

The papain/CLIK-148 coordinate system was employed as a model to study the interactions of a nonpeptide thiocarbazate inhibitor of cathepsin L ( 1). This small molecule inhibitor, a thiol ester containing a diacyl hydrazine functionality and one stereogenic center, was most active as the S-enantiomer, with an IC 50 of 56 nM; the R-enantiomer ( 2) displayed only weak activity (33 microM). Correspondingly, molecular docking studies with Extra Precision Glide revealed a correlation between score and biological activity for the two thiocarbazate enantiomers when a structural water was preserved. The molecular interactions between 1 and papain were very similar to the interactions observed for CLIK-148 ( 3a and 3b) with papain, especially with regard to the hydrogen-bonding and lipophilic interactions of the ligands with conserved residues in the catalytic binding site. Subsequent docking of virtual compounds in the binding site led to the identification of a more potent inhibitor ( 5), with an IC 50 of 7.0 nM. These docking studies revealed that favorable energy scores and correspondingly favorable biological activities could be realized when the virtual compound design included occupation of the S2, S3, and S1' subsites by hydrophobic and aromatic functionalities of the ligand, and at least three hydrogen bonding contacts between the ligand and the conserved binding site residues of the protein.