An Extended Polyanion Activation Surface in Insulin Degrading Enzyme

Insulin degrading enzyme (IDE) is believed to be the major enzyme that metabolizes insulin and has been implicated in the degradation of a number of other bioactive peptides, including amyloid beta peptide (Aβ), glucagon, amylin, and atrial natriuretic peptide. IDE is activated toward some substrates by both peptides and polyanions/anions, possibly representing an important control mechanism and a potential therapeutic target. A binding site for the polyanion ATP has previously been defined crystallographically, but mutagenesis studies suggest that other polyanion binding modes likely exist on the same extended surface that forms one wall of the substrate-binding chamber. Here we use a computational approach to define three potential ATP binding sites and mutagenesis and kinetic studies to confirm the relevance of these sites. Mutations were made at four positively charged residues (Arg 429, Arg 431, Arg 847, Lys 898) within the polyanion-binding region, converting them to polar or hydrophobic residues. We find that mutations in all three ATP binding sites strongly decrease the degree of activation by ATP and can lower basal activity and cooperativity. Computational analysis suggests conformational changes that result from polyanion binding as well as from mutating residues involved in polyanion binding. These findings indicate the presence of multiple polyanion binding modes and suggest the anion-binding surface plays an important conformational role in controlling IDE activity

Molecular mechanism of activation of human musk receptors OR5AN1 and OR1A1 by (R)-muscone and diverse other musk-smelling compounds

We acknowledge support from NSF (CHE-1265679) and NIH (5R01DC014423 subaward) (EB), NIH (5R01 DC014423) (HM), the European Reasearch Council (ERC) and the Engineering Science Research Council (EPSRC) (DO'H), FAPESP and CNPq (RAC), the Chinese Scholarship Council (CSC) for studentship support (MY), National Science Foundation (31070972) (HZ), Science and Technology Commission of Shanghai Municipality (16ZR1418300) (HZ), the Shanghai Eastern Scholar Program (J50201) (HZ). VSB thanks NIH grant 1R01GM106121-01A1 and computational time from NERSC.Understanding olfaction at the molecular level is challenging due to the lack of crystallographic models of odorant receptors (ORs). To better understand the molecular mechanism of OR activation, we focused on chiral (R)-muscone and other musk smelling odorants due to their great importance and widespread use in perfumery and traditional medicine, as well as environmental concerns associated with bioaccumulation of musks with estrogenic/antiestrogenic properties.  We experimentally and computationally examined the activation of human receptors OR5AN1 and OR1A1, recently identified as specifically responding to musk compounds.  OR5AN1 responds at nanomolar concentrations to musk ketone and robustly to macrocyclic sulfoxides and fluorine-substituted macrocyclic ketones; OR1A1 responds only to nitromusks. Structural models of OR5AN1 and OR1A1 based on quantum mechanics/molecular mechanics (QM/MM) hybrid methods were validated through direct comparisons with activation profiles from site-directed mutagenesis experiments and analysis of binding energies for 35 musk-related odorants.  The experimentally found chiral selectivity of OR5AN1 to (R)- over (S)-muscone was also computationally confirmed for muscone and fluorinated (R)-muscone analogs.  Structural models show that OR5AN1, highly responsive to nitromusks over macrocyclic musks, stabilizes odorants by hydrogen bonding to Tyr260 of transmembrane a-helix 6 and hydrophobic interactions with surrounding aromatic residues Phe105, Phe194 and, Phe207.  The binding of OR1A1 to nitromusks is stabilized by hydrogen bonding to Tyr258 along with hydrophobic interactions with surrounding aromatic residues Tyr251 and Phe206.  Hydrophobic/nonpolar and hydrogen bonding interactions contribute, respectively, 77% and 13% to the odorant binding affinities, as shown by an atom-based quantitative structure-activity relationship model.PostprintPeer reviewe

Disulfide Bridges Remain Intact while Native Insulin Converts into Amyloid Fibrils

Amyloid fibrils are β-sheet-rich protein aggregates commonly found in the organs and tissues of patients with various amyloid-associated diseases. Understanding the structural organization of amyloid fibrils can be beneficial for the search of drugs to successfully treat diseases associated with protein misfolding. The structure of insulin fibrils was characterized by deep ultraviolet resonance Raman (DUVRR) and Nuclear Magnetic Resonance (NMR) spectroscopy combined with hydrogen-deuterium exchange. The compositions of the fibril core and unordered parts were determined at single amino acid residue resolution. All three disulfide bonds of native insulin remained intact during the aggregation process, withstanding scrambling. Three out of four tyrosine residues were packed into the fibril core, and another aromatic amino acid, phenylalanine, was located in the unordered parts of insulin fibrils. In addition, using all-atom MD simulations, the disulfide bonds were confirmed to remain intact in the insulin dimer, which mimics the fibrillar form of insulin

Integrated Theoretical and Computational Approaches to Study Hydrolysis and Aggregation of Biomolecules

In this thesis, we have focused on the following two projects: (1) mechanisms of peptide hydrolysis, and (2) aggregation of biomolecules. A plethora of theoretical and computational techniques; quantum mechanics, molecular dynamics, and molecular docking have been employed to accomplish the goals of these projects. The selective hydrolytic cleavage of the extremely stable peptide bond (-(O=)C-NH-) of peptides and proteins by enzymes and their analogues is required in a wide range of biological, biotechnological and industrial applications. In the first project, we have elucidated the mechanisms of peptide hydrolysis by non-metallic (beta secretase (BACE)) and metallic (insulin degrading enzyme (IDE) and neprilysin (NEP)) enzymes and the synthetic analogues of metallo-enzymes. In addition, mechanisms of activation of IDE through small molecules and alteration in substrate specificity of NEP have been investigated. The aggregation of biomolecules has been implicated in a large number of neurological disorders. In the second project, the aggregation mechanisms of insulin and Alzheimer`s amyloid beta (Aβ) peptide have been studied. The interactions between the Aβ peptide and Ru complexes have also been explored. In another project, we have investigated the mechanism oxygen activation by Cu-containing amine oxidases (CAOs) by computing kinetic isotope effect (KIEs)

Comparison of Clinically Approved Molecules on SARS-CoV-2 Drug Target Proteins: A Molecular Docking Study

The new type of coronavirus, SARS-CoV-2 has affected more than 6.3 million people worldwide. Since the first day the virus has been spotted in Wuhan, China, there are numerous drug design studies conducted all over the globe. Most of these studies target the receptor-binding domain of spike protein of SASR-CoV-2, which is known to bind human ACE2 receptor and SARS-CoV-2 main protease, vital for the virus’ replication. However, there might be a third target, human furin protease, which cleaves the virus’ S1-S2 domains taking active role in its entry into the host cell. In this study we docked five clinically used drug molecules, favipiravir, hydroxychloroquine, remdesivir, lopinavir, and ritonavir onto three target proteins, receptor binding domain of SARS-CoV-2 spike protein, SARS-CoV-2 main protease, and human furin protease. Results of molecular docking simulations revealed that human furin protease might be targeted against COVID-19. Remdesivir, a nucleic acid derivative, strongly bound to the active site of this protease, suggesting this molecule can be used as a template for designing novel furin protease inhibitorsto fight with the disease. Protein-drug interactions revealed at the molecular level in this study can pave the way for better drug design for each specific target.</p

Elucidation of insulin degrading enzyme catalyzed site specific hydrolytic cleavage of amyloid beta peptide: a comparative density functional theory study

In this B3LYP study, the catalytic mechanisms for the hydrolysis of the three different peptide bonds (Lys28-Gly29, Phe19-Phe20, and His14-Gln15) of Alzheimer amyloid beta (Abeta) peptide by insulin-degrading enzyme (IDE) have been elucidated. For all these peptides, the nature of the substrate was found to influence the structure of the active enzyme-substrate complex. The catalytic mechanism is proposed to proceed through the following three steps: (1) activation of the metal-bound water molecule, (2) formation of the gem-diol intermediate, and (3) cleavage of the peptide bond. With the computed barrier of 14.3, 18.8, and 22.3 kcal/mol for the Lys28-Gly29, Phe19-Phe20, and His14-Gln15 substrates, respectively, the process of water activation was found to be the rate-determining step for all three substrates. The computed energetics show that IDE is the most efficient in hydrolyzing the Lys28-Gly29 (basic polar-neutral nonpolar) peptide bond followed by the Phe19-Phe20 (neutral nonpolar-neutral nonpolar) and His14-Gln15 (basic polar-neutral polar) bonds of the Abeta substrate

Which one among aspartyl protease, metallopeptidase, and artificial metallopeptidase is the most efficient catalyst in peptide hydrolysis?

In this comparative DFT study, the hydrolysis of a peptide bond (Phe1-Phe2) by the following three types of catalysts has been studied: (1) beta-secretase (BACE2), (2) matrix metalloproteinase (MMP) and insulin degrading enzyme (IDE), and (3) [Pd(H(2)O)(4)](2+) (I(MPC)) and [Pd(2)(mu-OH)([18]aneN(6))](3+) (I(DPC)). The computed energetics predict that among these catalysts, the Zn(2+) metal center containing MMP is the most efficient in catalyzing this reaction. The two active site aspartate residues containing BACE2 catalyze this reaction with 5.0 kcal/mol higher barrier than MMP. The substitution of a His ligand with Glu in the metal center of MMP generates the active site of IDE that catalyzes the reaction with a 6.9 kcal/mol higher barrier than MMP. Both artificial peptidases I(MPC) and I(DPC) catalyze this reaction with significantly high barriers of 35.4 and 31.0 kcal/mol, respectively. The computed energetics of all the catalysts are in line with the available experimental and theoretical data

Mechanism of peptide hydrolysis by co-catalytic metal centers containing leucine aminopeptidase enzyme: a DFT approach

In this density functional theory study, reaction mechanisms of a co-catalytic binuclear metal center (Zn1–Zn2) containing enzyme leucine aminopeptidase for two different metal bridging nucleophiles (H2O and –OH) have been investigated. In addition, the effects of the substrate (l-leucine-p-nitroanilide → l-leucyl-p-anisidine) and metal (Zn1 → Mg and Zn2 → Co, i.e., Mg1–Zn2 and Mg1–Co2 variants) substitutions on the energetics of the mechanism have been investigated. The general acid/base mechanism utilizing a bicarbonate ion followed by this enzyme is divided into two steps: (1) the formation of the gem-diolate intermediate, and (2) the cleavage of the peptide bond. With the computed barrier of 17.8 kcal/mol, the mechanism utilizing a hydroxyl nucleophile was found to be in excellent agreement with the experimentally measured barrier of 18.7 kcal/mol. The rate-limiting step for reaction with l-leucine-p-nitroanilide is the cleavage of the peptide bond with a barrier of 17.8 kcal/mol. However, for l-leucyl-p-anisidine all steps of the mechanism were found to occur with similar barriers (18.0–19.0 kcal/mol). For the metallovariants, cleavage of the peptide bond occurs in the rate-limiting step with barriers of 17.8, 18.0, and 24.2 kcal/mol for the Zn1–Zn2, Mg1–Zn2, and Mg1–Co2 enzymes, respectively. The nature of the metal ion was found to affect only the creation of the gem-diolate intermediate, and after that all three enzymes follow essentially the same energetics. The results reported in this study have elucidated specific roles of both metal centers, the nucleophile, indirect ligands, and substrates in the catalytic functioning of this important class of binuclear metallopeptidases

Novel N-(1-thia-4-azaspiro[4.5]decan-4-yl)carboxamide derivatives as potent and selective influenza virus fusion inhibitors

Hemagglutinin is the surface protein of the influenza virus that mediates both binding and penetration of the virus into host cells. We here report on the synthesis and structure-activity relationship of some novel N-(1-thia-4-azaspiro[4.5]decan-4-yl)-carboxamide compounds carrying the 5-chloro-2-methoxybenzamide structure, designed as influenza virus fusion inhibitors. The carboxamides (1a-h, 2a-h) have a similar backbone structure as the fusion inhibitors that we reported on previously. Compounds 2b and 2d displayed inhibitory activity against influenza A/H3N2 virus replication (average antiviral EC50 : 2.1 µM for 2b and 3.4 µM for 2d). Data obtained in the hemolysis inhibition assay supported that these compounds act as inhibitors of the influenza virus hemagglutinin-mediated fusion process.status: publishe

Unraveling the photoluminescence response of light-switching ruthenium(II) complexes bound to amyloid-β

Photoluminescent molecules are widely used for real-time monitoring of peptide aggregation. In this Article, we detail both experimental and computational modeling to elucidate the interaction between [Ru(bpy)2dppz](2+) and amyloid-β (Aβ(1-40)) aggregates. The transition from monomeric to fibrillar Aβ is of interest in the study of Alzheimer's disease. Concentration-dependent experiments allowed the determination of a dissociation constant of 2.1 μM, while Job plots provided a binding stoichiometry of 2.6 Aβ monomers per [Ru(bpy)2dppz](2+). Our computational approach that combines molecular docking (both rigid and flexible) and all-atom molecular dynamics (MD) simulations predicts that the hydrophobic cleft between Val18 and Phe20 is a plausible binding site, which could also explain the increase in photoluminescence of [Ru(bpy)2dppz](2+) upon binding. This binding site is parallel to the fibril axis, in marked contrast to the binding site of these complexes in DNA (perpendicular to the DNA axis). Other binding sites may exist at the edges of the Aβ fibril, but they are actually of low abundance in an Aβ fibril several micrometers long. The assignment of the binding site was confirmed by binding studies in an Aβ fragment (Aβ(25-35)) that lacked the amino acids necessary to form the binding site. The agreement between the experimental and computational work is remarkable and provides a general model that can be used for studying the interaction of amyloid-binding molecules to Aβ