The role of membrane thickness in charged protein–lipid interactions

Charged amino acids are known to be important in controlling the actions of integral and peripheral membrane proteins and cell disrupting peptides. Atomistic molecular dynamics studies have shed much light on the mechanisms of membrane binding and translocation of charged protein groups, yet the impact of the full diversity of membrane physico-chemical properties and topologies has yet to be explored. Here we have performed a systematic study of an arginine (Arg) side chain analog moving across saturated phosphatidylcholine (PC) bilayers of variable hydrocarbon tail length from 10 to 18 carbons. For all bilayers we observe similar ion-induced defects, where Arg draws water molecules and lipid head groups into the bilayers to avoid large dehydration energy costs. The free energy profiles all exhibit sharp climbs with increasing penetration into the hydrocarbon core, with predictable shifts between bilayers of different thickness, leading to barrier reduction from 26 kcal/mol for 18 carbons to 6 kcal/mol for 10 carbons. For lipids of 10 and 12 carbons we observe narrow transmembrane pores and corresponding plateaus in the free energy profiles. Allowing for movements of the protein and side chain snorkeling, we argue that the energetic cost for burying Arg inside a thin bilayer will be small, consistent with recent experiments, also leading to a dramatic reduction in pK_a shifts for Arg. We provide evidence that Arg translocation occurs via an ion-induced defect mechanism, except in thick bilayers (of at least 18 carbons) where solubility-diffusion becomes energetically favored. Our findings shed light on the mechanisms of ion movement through membranes of varying composition, with implications for a range of charged protein–lipid interactions and the actions of cell-perturbing peptides. This article is part of a Special Issue entitled: Membrane protein structure and function.     

The role of tryptophan side chains in membrane protein anchoring and hydrophobic mismatch

Tryptophan (Trp) is abundant in membrane proteins, preferentially residing near the lipid–water interface where it is thought to play a significant anchoring role. Using a total of 3 μs of molecular dynamics simulations for a library of hydrophobic WALP-like peptides, a long poly-Leu α-helix, and the methyl-indole analog, we explore the thermodynamics of the Trp movement in membranes that governs the stability and orientation of transmembrane protein segments. We examine the dominant hydrogen-bonding interactions between the Trp and lipid carbonyl and phosphate moieties, cation–π interactions to lipid choline moieties, and elucidate the contributions to the thermodynamics that serve to localize the Trp, by ~ 4 kcal/mol, near the membrane glycerol backbone region. We show a striking similarity between the free energy to move an isolated Trp side chain to that found from a wide range of WALP peptides, suggesting that the location of this side chain is nearly independent of the host transmembrane segment. Our calculations provide quantitative measures that explain Trp's role as a modulator of responses to hydrophobic mismatch, providing a deeper understanding of how lipid composition may control a range of membrane active peptides and proteins.     

Isoform-selective inhibition of chrysin towards human cytochrome P450 1A2. Kinetics analysis,molecular docking,and molecular dynamics simulations

Our kinetics studies demonstrated that the nature product chrysin exhibited a high inhibitory affinity of 54 nM towards human cytochrome P450 1A2 and was comparable to α-naphthoflavone (49 nM), whereas it represented a moderate affinity of 5225 nM against human cytochrome P450 2C9. However, it remains unclear how this inhibitor selectively binds 1A2. To better understand the isoform selectivity of chrysin, molecular docking and molecular dynamics simulations were performed. Chrysin formed a strong H-bond with Asp313 of 1A2. The stacking interactions with Phe226 also contributed to its tight binding to 1A2. The larger and much more open active site architectures of 2C9 may explain the weaker inhibitory affinity of chrysin towards 2C9. The predicted binding free energies suggest that chrysin preferred 1A2 (ΔG_{bind, pred} = ?23.11 kcal/mol) to 2C9 (?20.41 kcal/mol). Additionally, the present work revealed that 7-hydroxy-flavone bound to 1A2 in a similar pattern as chrysin and represented a slightly less negative predicted binding free energy, which was further validated by our kinetics analysis (IC₅₀ = 240 nM). Results of the study can provide insight for designing novel isoform-selective 1A2 inhibitors.     

riDOM,a cell penetrating peptide. Interaction with phospholipid bilayers

Melittin is an amphipathic peptide which has received much attention as a model peptide for peptide–membrane interactions. It is however not suited as a transfection agent due to its cytolytic and toxicological effects. Retro-inverso-melittin, when covalently linked to the lipid 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine (riDOM), eliminates these shortcomings. The interaction of riDOM with phospholipid membranes was investigated with circular dichroism (CD) spectroscopy, dynamic light scattering, ζ-potential measurements, and high-sensitivity isothermal titration calorimetry. riDOM forms cationic nanoparticles with a diameter of ~ 13 nm which are well soluble in water and bind with high affinity to DNA and lipid membranes. When dissolved in bilayer membranes, riDOM nanoparticles dissociate and form transient pores. riDOM-induced membrane leakiness is however much reduced compared to that of authentic melittin. The secondary structure of the ri-melittin is not changed when riDOM is transferred from water to the membrane and displays a large fraction of β-structure. The ³¹P NMR spectrum of the nanoparticle is however transformed into a typical bilayer spectrum. The Gibbs free energy of riDOM binding to bilayer membranes is − 8.0 to − 10.0 kcal/mol which corresponds to the partition energy of just one fatty acyl chain. Half of the hydrophobic surface of the riDOM lipid extension with its 2 oleic acyl chains is therefore involved in a lipid–peptide interaction. This packing arrangement guarantees a good solubility of riDOM both in the aqueous and in the membrane phase. The membrane binding enthalpy is small and riDOM binding is thus entropy-driven.     

Effect of physicochemical properties of peptides from soy protein on their antimicrobial activity

Antimicrobial peptides (AMPs) kill microbial cells through insertion and damage/permeabilization of the cytoplasmic cell membranes and has applications in food safety and antibiotic replacement. Soy protein is an attractive, abundant natural source for commercial production of AMPs. In this research, explicit solvent molecular dynamics (MD) simulation was employed to investigate the effects of (i) number of total and net charges, (ii) hydrophobicity (iii) hydrophobic moment and (iv) helicity of peptides from soy protein on their ability to bind to lipid bilayer and their transmembrane aggregates to form pores. Interaction of possible AMP segments from soy protein with 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine/1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol (POPC/POPG) bilayers, a mimic of bacterial cell membrane, was investigated. Pore formation was insensitive to helicity and occurred for hydrophobicity threshold in the range of −0.3–0 kcal/mol, hydrophobic moment threshold of 0.3 kcal/mol, net charge threshold of 2. Though low hydrophobicity and high number of charges help in the formation of water channel for transmembrane aggregates, insertion of peptides with these properties requires overcome of energy barrier, as shown by potential of mean force calculations, thereby resulting in low antimicrobial activity. Experimental evaluation of antimicrobial activity of these peptides against Gram positive L. monocytogenes and Gram negative E. coli as obtained by spot-on-lawn assay was consistent with simulation results. These results should help in the development of guidelines for selection of peptides with antimicrobial activity based on their physicochemical properties.     

Towards understanding the unbound state of drug compounds: Implications for the intramolecular reorganization energy upon binding

There has been an explosion of structural information for pharmaceutical compounds bound to biological targets, but the conformations and dynamics of compounds free in solution are poorly characterized, if at all. Yet, knowledge of the unbound state is essential to understand the fundamentals of molecular recognition, including the much debated conformational intramolecular reorganization energy of a compound upon binding (ΔE_Reorg). Also, dependable observation of the unbound compounds is important for ligand-based drug discovery, e.g. with pharmacophore modelling. Here, these questions are addressed with long (⩾0.5 μs) state-of-the-art molecular dynamics (MD) simulations of 26 compounds (including 7 approved drugs) unbound in explicit solvent. These compounds were selected to be chemically diverse, with a range of flexibility, and good quality bioactive X-ray structures. The MD-simulated free compounds are compared to their bioactive structure and conformers generated with ad hoc sampling in vacuo or with implicit generalized Born (GB) aqueous solvation models. The GB conformational models clearly depart from those obtained in explicit solvent, and suffer from conformational collapse almost as severe as in vacuo. Thus, the global energy minima in vacuo or with GB are not suitable representations of the unbound state, which can instead be extensively sampled by MD simulations. Many, but not all, MD-simulated compounds displayed some structural similarity to their bioactive structure, supporting the notion of conformational pre-organization for binding. The ligand–protein complexes were also simulated in explicit solvent, to estimate ΔE_Reorg as an enthalpic difference ΔH_Reorg between the intramolecular energies in the bound and unbound states. This fresh approach yielded ΔH_Reorg values ⩽ 6 kcal/mol for 18 out of 26 compounds. For three particularly polar compounds 15 ⩽ ΔH_Reorg ⩽ 20 kcal/mol, supporting the notion that ΔH_Reorg can be substantial. Those large ΔH_Reorg values correspond to a redistribution of electrostatic interactions upon binding. Overall, the study illustrates how MD simulations offer a promising avenue to characterize the unbound state of medicinal compounds.     

Bean peptides have higher in silico binding affinities than ezetimibe for the N-terminal domain of cholesterol receptor Niemann-Pick C1 Like-1

Niemann-Pick C1 like-1 (NPC1L1) mediates cholesterol absorption at the apical membrane of enterocytes through a yet unknown mechanism. Bean, pea, and lentil proteins are naturally hydrolyzed during digestion to produce peptides. The potential for pulse peptides to have high binding affinities for NPC1L1 has not been determined. In this study , in silico binding affinities and interactions were determined between the N-terminal domain of NPC1L1 and 14 pulse peptides (5≥ amino acids) derived through pepsin-pancreatin digestion. Peptides were docked in triplicate to the N-terminal domain using docking program AutoDock Vina, and results were compared to those of ezetimibe, a prescribed NPC1L1 inhibitor. Three black bean peptides (−7.2 to −7.0 kcal/mol) and the cowpea bean dipeptide Lys-Asp (−7.0 kcal/mol) had higher binding affinities than ezetimibe (−6.6 kcal/mol) for the N-terminal domain of NPC1L1. Lentil and pea peptides studied did not have high binding affinities. The common bean peptide Tyr-Ala-Ala-Ala-Thr (−7.2 kcal/mol), which can be produced from black or navy bean proteins, had the highest binding affinity. Ezetimibe and peptides with high binding affinities for the N-terminal domain are expected to interact at different locations of the N-terminal domain. All high affinity black bean peptides are expected to have van der Waals interactions with SER130, PHE136, and LEU236 and a conventional hydrogen bond with GLU238 of NPC1L1. Due to their high affinity for the N-terminal domain of NPC1L1, black and cowpea bean peptides produced in the digestive track have the potential to disrupt interactions between NPC1L1 and membrane proteins that lead to cholesterol absorption.     

Influenza Hemifusion Phenotype Depends on Membrane Context: Differences in Cell–Cell and Virus–Cell Fusion

Influenza viral entry into the host cell cytoplasm is accomplished by a process of membrane fusion mediated by the viral hemagglutinin protein. Hemagglutinin acts in a pH-triggered fashion, inserting a short fusion peptide into the host membrane followed by refolding of a coiled-coil structure to draw the viral envelope and host membranes together. Mutations to this fusion peptide provide an important window into viral fusion mechanisms and protein–membrane interactions. Here, we show that a well-described fusion peptide mutant, G1S, has a phenotype that depends strongly on the viral membrane context. The G1S mutant is well known to cause a “hemifusion” phenotype based on experiments in transfected cells, where cells expressing G1S hemagglutinin can undergo lipid mixing in a pH-triggered fashion similar to virus but will not support fusion pores. We compare fusion by the G1S hemagglutinin mutant expressed either in cells or in influenza virions and show that this hemifusion phenotype occurs in transfected cells but that native virions are able to support full fusion, albeit at a slower rate and 10–100 × reduced infectious titer. We explain this with a quantitative model where the G1S mutant, instead of causing an absolute block of fusion, alters the protein stoichiometry required for fusion. This change slightly slows fusion at high hemagglutinin density, as on the viral surface, but at lower hemagglutinin density produces a hemifusion phenotype. The quantitative model thus reproduces the observed virus–cell and cell–cell fusion phenotypes, yielding a unified explanation where membrane context can control the observed viral fusion phenotype.     

Inhibition of tyrosinase by gastrodin: An integrated kinetic-computational simulation analysis

We studied the inhibitory effect of gastrodin on tyrosinase using inhibition kinetics and computational simulation. Gastrodin reversibly inhibited tyrosinase in a mixed-type manner with K_i = 123.8 ± 20.2 mM. Time-interval kinetics revealed the inhibition to be a first-order process with mono- and bi-phasic components. Using AutoDock Vina, we calculated a binding energy of ?6.3 kcal/mol for gastrodin and tyrosinase, and we performed a molecular dynamics simulation of the tyrosinase–gastrodin interaction. The simulation results suggested that gastrodin interacts primarily with histidine residues in the active site. A 10-ns molecular dynamics simulation showed that one copper ion in the tyrosinase active site was responsible for the interaction with gastrodin. Our study provides insight into the inhibition of tyrosinase by the hydroxyl groups of gastrodin. A combination of inhibition kinetics and computational calculations may help to confirm the inhibitory action of gastrodin on tyrosinase and define the mechanisms of inhibition.     

An in silico approach unveils the potential of antiviral compounds in preclinical and clinical trials as SARS-CoV-2 omicron inhibitors

The increased transmissibility and highly infectious nature of the new variant of concern (VOC) that is severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron and lack of effective therapy need the rapid discovery of therapeutic antivirals against it. The present investigation aimed to identify antiviral compounds that would be effective against SARS-CoV-2 Omicron. In this study, molecular docking experiments were carried out using the recently reported experimental structure of omicron spike protein in complex with human angiotensin-converting enzyme 2 (ACE2) and various antivirals in preclinical and clinical trial studies. Out of 36 tested compounds, Abemaciclib, Dasatinib and Spiperone are the three top-ranked molecules which scored binding energies of ?10.08 kcal/mol, ?10.06 kcal/mol and ?9.54 kcal/mol respectively. Phe338, Asp339, and Asp364 are crucial omicron receptor residues involved in hydrogen bond interactions, while other residues were mostly involved in hydrophobic interactions with the lead molecules. The identified lead compounds also scored well in terms of drug-likeness. Molecular dynamics (MD) simulation, essential dynamics (ED) and entropic analysis indicate the ability of these molecules to modulate the activity of omicron spike protein. Therefore, Abemaciclib, Dasatinib and Spiperone are likely to be viable drug-candidate molecules that can block the interaction between the omicron spike protein and the host cellular receptor ACE2. Though our findings are compelling, more research into these molecules is needed before they can be employed as drugs to treat SARS-CoV-2 omicron infections.     

Ethyl isovalerate synthesis using Candida rugosa lipase immobilized on silica nanoparticles prepared in nonionic reverse micelles

Uniform and monodispersed silica nanoparticles were synthesized with a mean diameter of 100 ± 20 nm as analyzed by Transmission Electron Microscopy (TEM). Glutaraldehyde was used as a coupling agent for efficient binding of the lipase onto the silica nanoparticles. For the hydrolysis of pNPP at pH 7.2, the activation energy within 25–40 °C for free and immobilized lipase was 7.8 and 1.25 KJ/mol, respectively. The V_max and K_m of immobilized lipase at 25 °C for pNPP hydrolysis were found to be 212 μmol/min/mg and 0.3 mM, whereas those for free lipase were 26.17 μmol/min and 1.427 mM, respectively. The lower activation energy of immobilized lipase in comparison to free lipase suggests a change in conformation of the enzyme leading to a requirement for lower energy on the surface of the nanoparticles. A better yield (7 fold higher) of ethyl isovalerate was observed using lipase immobilized onto silica nanoparticles in comparison to free lipase.     

Comparison of lipid-dependent bilayer insertion of pHLIP and its P20G variant

The ability of the pH-Low Insertion Peptide (pHLIP) to insert into lipid membranes in a transbilayer conformation makes it an important tool for targeting acidic diseased tissues. pHLIP can also serve as a model template for thermodynamic studies of membrane insertion. We use intrinsic fluorescence and circular dichroism spectroscopy to examine the effect of replacing pHLIP's central proline on the pH-triggered lipid-dependent conformational switching of the peptide. We find that the P20G variant (pHLIP-P20G) has a higher helical propensity than the native pHLIP (pHLIP-WT), in both water:organic solvent mixtures and in the presence of lipid bilayers. Spectral shifts of tryptophan fluorescence reveal that with both pHLIP-WT and pHLIP-P20G, the deeply penetrating interfacial form (traditionally called State II) is populated only in pure phosphocholine bilayers. The presence of either anionic lipids or phosphatidylethanolamine leads to a much shallower penetration of the peptide (referred to here as State II^S, for “shallow”). This novel state can be differentiated from soluble state by a reduction in accessibility of tryptophans to acrylamide and by FRET to vesicles doped with Dansyl-PE, but not by a spectral shift in fluorescence emission. FRET experiments indicate free energies for interfacial partitioning range from 6.2 to 6.8 kcal/mol and are marginally more favorable for pHLIP-P20G. The effective pKa for the insertion of both peptides depends on the lipid composition, but is always higher for pHLIP-P20G than for pHLIP-WT by approximately one pH unit, which corresponds to a difference of 1.3 kcal/mol in free energy of protonation favoring insertion of pHLIP-P20G.     

Crystal structure and molecular modeling study of N-carbamoylsarcosine amidase Ta0454 from Thermoplasma acidophilum

A crystal structure of the putative N-carbamoylsarcosine amidase (CSHase) Ta0454 from Thermoplasma acidophilum was solved by single-wavelength anomalous diffraction and refined at a resolution of 2.35 Å. CSHases are involved in the degradation of creatinine. Ta0454 shares a similar fold and a highly conserved C-D-K catalytic triad (Cys123, Asp9, and Lys90) with the structures of three cysteine hydrolases (PDB codes 1NBA, 1IM5, and 2H0R). Molecular dynamics (MD) simulations of Ta0454/N-carbamoylsarcosine and Ta0454/pyrazinamide complexes were performed to determine the structural basis of the substrate binding pattern for each ligand. Based on the MD-simulated trajectories, the MM/PBSA method predicts binding free energies of ?24.5 and ?17.1 kcal/mol for the two systems, respectively. The predicted binding free energies suggest that Ta0454 is selective for N-carbamoylsarcosine over pyrazinamide, and zinc ions play an important role in the favorable substrate bound states.     

Effect of d-amino acids on collagen fibrillar assembly and stability: Experimental and modelling studies

The effect of selected d-amino acids (d-AAs) on collagen with 1-ethyl-3-(3-dimethylamino propyl)carbodiimide (EDC)/N-hydroxysuccinimide (NHS) initiated crosslinking is evaluated by using experimental and modelling tools. The experimental results suggest that d-Lysine (d-Lys) plays a pivotal role in the self-assembly and conformation of collagen fibrils than d-Alanine (d-Ala) and d-Glutamic acid (d-Glu). The SDS-PAGE, absorption spectrum and viscosity measurements indicate significant differences in the d-Lys crosslinked collagen when compared to other systems. The CD spectra show an increase in the peak intensity at 220 nm in the presence of d-Lys, which could be due to increase in propensity of the structure to form a triple helix. Modelling studies indicated that d-Lys bind with collagen-like peptide (CLP) through multiple H-bonding and hydrophobic interactions. d-Lys has the lowest binding energy (−4.2 kcal/mol, indicating strongest interactions) when compared to d-Ala and d-Glu (−3.6 and −3.7 kcal/mol, respectively). Orientational changes in the collagenase on CLP-d-Lys are observed which may decrease its accessibility to degradation and stabilise CLP against the action of the former. d-Lys has the lowest binding energy and improved fibrillar-assembly and staggered alignment without the undesired structural stiffness and aggregations. The information derived from the present study could help in designing heterochiral collagen-based biomaterial.     

Conformation and Trimer Association of the Transmembrane Domain of the Parainfluenza Virus Fusion Protein in Lipid Bilayers from Solid-State NMR: Insights into the Sequence Determinants of Trimer Structure and Fusion Activity

Enveloped viruses enter cells by using their fusion proteins to merge the virus lipid envelope and the cell membrane. While crystal structures of the water-soluble ectodomains of many viral fusion proteins have been determined, the structure and assembly of the C-terminal transmembrane domain (TMD) remains poorly understood. Here we use solid-state NMR to determine the backbone conformation and oligomeric structure of the TMD of the parainfluenza virus 5 fusion protein. ¹³C chemical shifts indicate that the central leucine-rich segment of the TMD is α-helical in POPC/cholesterol membranes and POPE membranes, while the Ile- and Val-rich termini shift to the β-strand conformation in the POPE membrane. Importantly, lipid mixing assays indicate that the TMD is more fusogenic in the POPE membrane than in the POPC/cholesterol membrane, indicating that the β-strand conformation is important for fusion by inducing membrane curvature. Incorporation of para-fluorinated Phe at three positions of the α-helical core allowed us to measure interhelical distances using ¹⁹F spin diffusion NMR. The data indicate that, at peptide:lipid molar ratios of ~ 1:15, the TMD forms a trimeric helical bundle with inter-helical distances of 8.2–8.4 Å for L493F and L504F and 10.5 Å for L500F. These data provide high-resolution evidence of trimer formation of a viral fusion protein TMD in phospholipid bilayers, and indicate that the parainfluenza virus 5 fusion protein TMD harbors two functions: the central α-helical core is the trimerization unit of the protein, while the two termini are responsible for inducing membrane curvature by transitioning to a β-sheet conformation.     

Computational study of bindings of HL9, a nonapeptide fragment of human lysozyme, to HIV-1 fusion protein gp41

HL9 is a nonapeptide fragment of human lysozyme which has been shown to have anti-HIV-1 activity in nanomolar concentration. This study aims to explain this inhibitory activity by using molecular dynamics (MD) simulation, focusing on the ectodomain of gp41, the envelope glycoprotein of HIV-1 crucial to membrane fusion. It was found that in HL9, two Trp residues separated by two others occupy the conserved hydrophobic pocket on gp41 and thus inhibit fusion in dominant-negative manner. Detailed HL9-gp41 binding interactions and free energies of binding were obtained through MD simulation and solvated interaction energies (SIE) calculation, giving a binding free energy of −8.25 kcal/mol which is in close agreement with the experimental value of −9.96 kcal/mol. Since C-helical region (C34) of gp41 also has two Trp residues separated by two others, this arrangement may be generalised and used to scan peptide library and to find those having similar manner of inhibition.     

Stapled peptide-based membrane fusion inhibitors of hepatitis C virus

The strategy of peptide stapling was used to develop new molecules to inhibit the hepatitis C virus infection via disrupting the binding of HCV envelope glycoprotein E2 with human cell surface protein CD81. The peptide sequence was designed based on the large extra-cellular loop of CD81 with known importance in the HCV E2 binding interaction. Our results showed that the stapled peptides exhibited significantly higher α-helicity and proteolytic stability as compared to their linear peptide counterpart. The optimal compound was found to have an EC₅₀ value of ca. 17–39 μM against different HCV subtypes and represented a new HCV membrane fusion inhibitor.     

Docking based virtual screening and molecular dynamics study to identify potential monoacylglycerol lipase inhibitors

Monoacylglycerol lipase (MAGL) is one of the key enzymes of the endocannabinoid system (ECS). It hydrolyzes one of the major endocannabinoid, 2-arachidonoylglycerol (2-AG), an endogenous full agonist at G protein coupled cannabinoid receptors CB1 and CB2. Numerous studies showed that MGL inhibitors are potentially useful for the treatment of pain, inflammation, cancer and CNS disorders. These provocative findings suggested that pharmacological inhibition of MAGL function may confer significant therapeutic benefits. In this study, we presented hybrid ligand and structure-based approaches to obtain a novel set of virtual leads as MAGL inhibitors. The constraints used in this study, were Glide score, binding free energy estimates and ADME properties to screen the ZINC database, containing approximately 21 million compounds. A total of seven virtual hits were obtained, which showed significant binding affinity towards MAGL protein. Ligand, ZINC24092691 was employed in complex form with the protein MAGL, for molecular dynamics simulation study, because of its excellent glide score, binding free energy and ADME properties. The RMSD of ZINC24092691 was observed to stay at 0.1 nm (1 Å) in most of the trajectories, which further confirmed its ability to inhibit the protein MAGL. The hits were then evaluated for their ability to inhibit human MAGL. The compound ZINC24092691 displayed the noteworthy inhibitory activity reducing MAGL activity to 21.15% at 100 nM concentration, with an IC₅₀ value of 10 nM.     

Synthesis and biological activity of oxadiazole and triazolothiadiazole derivatives as tyrosinase inhibitors

A series of 16 oxadiazole and triazolothiadiazole derivatives were designed, synthesized and evaluated as mushroom tyrosinase inhibitors. Five derivatives were found to display high inhibition on the tyrosinase activity ranging from 0.87 to 1.49 μM. Compound 5 exhibited highest tyrosinase inhibitory activity with an IC₅₀ value of 0.87 ± 0.16 μM. The in silico protein–ligand docking using autodock 4.1 was successfully performed on compound 5 with significant binding energy value of ?5.58 kcal/mol. The docking results also showed that the tyrosinase inhibition might be due to the metal chelating effect by the presence of thione functionality in compounds 1–5. Further studies revealed that the presence of hydrophobic group such as cycloamine derivatives played a major role in the inhibition. Piperazine moiety in compound 5 appeared to be involved in an extensive hydrophobic contact and a 2.9 Å hydrogen bonding with residue Glu 182 in the active site.