IMSJ International Marine Science Journal 2643-0282Open Access PubUnited States10.14302/issn.2643-0282.imsj-18-2448IMSJ-18-2448research-articleIn Silico Inhibition of Essential Candida albicans Proteins by Arenicin, a Marine Antifungal PeptideLíviado Carmo Silva1AndreiaMarcelino Barbosa2IsabelePereira Tannous1ThaynaraGonzaga Santos1JulianaSantana de Curcio1LucasNojosa Oliveira1MarielleGarcia Silva1AmandaAlves de Oliveira1RaisaMelo Lima1KleberSantiago Freitas e Silva1*Biological Sciences Institute, Federal University of Goiás, BrazilReplicon, Pontifical Catholic University of Goiás, BrazilCorresponding authorCesarAmaral1DNA Diagnostic Laboratory, University of the State of Rio de Janeiro, Brazil
Kleber Santiago Freitas e Silva, KSF, Biological Sciences Institute, Federal University of Goiás, Brazil, Address: Campus Samambaia, Goiânia - GO, Phone: 74710-310 (fax): 5562 32820744, Email: smallbinho@hotmail.com
The authors have declared that no competing interests exist.
Fungal infections increased substantially in the last years, becoming a relevant public health problem. Many of these infections account for high rates of morbidity and mortality. The emergence of resistant fungal clinical isolates have also motivate studies to find new antifungal therapies. Candida albicans is an oportunistic pathogen and affects a great number of immunocompromised patients worldwide. The marine ecosystem has been considered a rich source of bioactive metabolites due to the complexity and originality of its structures. Proteins and peptides from marine organisms have been shown to have antiviral, anti-inflammatory, antimalarial, anticancer, antimicrobial and antifungal properties. Arenicins are antimicrobial peptides isolated from the marine lugworm Arenicola marina with 21 amino acid residues in a β-hairpin structure. Dihydrofolate reductase, exo-b-(1,3)-glucanase and sterol 14α-demethylase are essential C. albincas enzymes that take part in DNA, cell wall and membrane metabolism, respectively. The present study evaluates the interaction of arenicin with important enzymes of C. albicans related to cell wall, ergosterol and DNA metabolism in order to elucidate possible molecular targets. We showed through an in silico approach, that a single compound from a marine worm (A. marina), can bind to three C. albicans essential proteins. The interaction occurs in regions inside the active site or at least near, with amino acid residues evaluated as hot spots. Arenicin is a new promising antifugal drug. The next step is to investigate protein-protein interactions performed by DHFR, EBG and CYP51 and assess whether arenicin is able to disrupt essential interaction or not.
Fungal infections increased substantially in the last years, becoming a relevant public health problem. Many of these infections account for high rates of morbidity and mortality. Candidiasis, an opportunistic fungal infection is caused by fungi belonging to the Candida genus 1. Although the antifungal agents used in clinical treatments of mycosis have sometimes proved efficient, they are restricted to few targets in fungal cells with high toxicity to humans 2. In addition, the problem increases with the emergence of resistant clinical isolates 34, which makes the search for new antifungal therapies extremely relevant.
Since the discovery of spongouridine, molecular model of the first antiviral used as a therapeutic resource for Acquired Immunodeficiency Syndrome (AIDS), the biodiversity of the marine ecosystem has been considered a rich source of bioactive metabolites due to the complexity and originality of its structures 5. In this scenario, peptides from marine organisms have been shown to have antiviral, anti-inflammatory, antimalarial, anticancer, antimicrobial and antifungal properties 678.
Arenicins are antimicrobial peptides isolated from the marine lugworm Arenicola marina with 21 amino acid residues in a β-hairpin structure Figure 1. In vitro and in vivo studies have shown that arenicins exhibit a very potent bactericidal and antifungal effect against multi-resistant gram-negative bacteria and Candida albicans, respectively 9. In bacterial, Arenicin interacts with membrane lipids promoting membrane permeabilization and detachment besides affecting cytoplasm metabolism 10. However, the interaction of arenicin with important antifungal targets from Candida albicans is unknown.
The arenicin structure. The conformational structure of the peptide folds into an antiparallel beta-sheet comprising two strands, which are important for its antimicrobial activity. Arenicin has a high affinity for lipids and disrupt fungal and other microorganisms membrane.
Dihydrofolate reductase (DHFR) catalyzes the reduction of dihydrofolic acid to tetrahydrofolic acid 11. The enzyme acts on the synthesis of purines, thymidylic acid, amino acids and the regulation of tetrahydrofolate levels in the cell, participating in proliferation and growth 12. The protein exo-b-(1,3)-glucanase (EBG) takes part in cell expansion and differentiation through cell wall remodeling. EBG is also secreted during pathogenesis 13. Finally, CYP51 (sterol 14α-demethylase) is a major drug target in the cytochrome P450 superfamily. This enzyme is related to ergosterol and membranes metabolism 14. All the enzymes described above are important to survivability and pathogenesis of Candida species. The identification of novel therapeutics can drive a more efficient treatment of immunocompromised patients that are susceptible to fungal diseases.
Molecular anchoring is an important computational tool in drug design and targeting. The purpose of the ligand-protein anchor is to predict the best interaction orientation of a linker with a protein of known three-dimensional structure 15. Thus, in this study we verified the interaction of arenicin with important enzymes of C. albicans related to cell wall, ergosterol and DNA metabolism in order to elucidate possible molecular targets.
Material and Methods
All the 3-D structures used in the analysis are available in the PDB (protein databank; https://www.rcsb.org/). KBDOCK server were used in order to assess protein domains and possible interaction between protein domains 16. The protein docking was performed by ClusPro 17. We used PyMol (https://pymol.org) for the visualization of the interface of interaction, the visualization of hot spots and for creating the figures presented in this manuscript .The hot spots in the proteins under study were identified by KFC2. The server offers an automated analysis of a protein complex interface. The server analyses the structural environment around amino acid residues and checks for already known hot spots environments determined experimentally. The hot spot prediction is based on characteristics regarding conformation specificity (K-FADE) and biochemical features such as hydrophobicity (K-CON) 1819.
Results and DiscussionThe Arenicin Properties are Defined by its Amino Acid Composition
Arginine and valine account for 54% of the amino acid residues present in the arenicin structure Figure 2 and they account for the properties presented by the peptide. The large amount of arginine residues has a role in maintaining the conformational state of the peptide and its overall charge 202122. Moreover, arginine residues may interact with active sites of proteins 232425 that bind to phosphorylated substrates through hydrogen bonds. The large amount of valine residues helps the peptide folds into its β-hairpin structure Figure 326 and this is related to the presence of two non-hydrogen substituent attached to the valine C-beta carbon. In addition, valine plays a role in substrate recognition of hydrophobic ligands such as lipids 2728, which could explain the high affinity the peptide has for membranes.
The arenicin arginine and valine residues. Arginine and valine amino acids correspond to 54% of residues in the arenicin structure. They contribute to the conformational state of the peptide, the overall charge and affinity for ligands. Red – arginine; blue – valine; green – other amino acid residues.
Dihydrofolate reductase (DHFR) interaction with arenicin. DHFR is a homodimeric protein and the inhibitor marine compound arenicin binds to a loop within the active site of the protein. Arenicin anchors inside the active site and alter slightly the conformational structure of DHFR. Green – DHFR; red – arenicin.
Arenicin Binds to Residues in the Active Site of Dihydrofolate Reductase
The main feature of the conformational structure of DHFR is a central region comprised by beta-pleated sheets 29. The active site is surrounded by a domain that contain loops and regulate the interaction pattern with partner proteins and compounds 30. An efficient strategy for inhibiting DHFR from C. albicans is to target essential amino acid residues at the active site of the enzyme 313233. Arenicin interacts with a loop inside the active site of DHFR Figure 4. There are five main amino acid residues (Ile-Pro-Gln-Lys-Phe) that contribute considerably for the stability of the interaction between DHFR and the arenicin peptide and they have been shown to link to other inhibitor compounds as well 33.
Arenicin anchors inside the active site of dihydrofolate reductase. The loop in the active site of DHFR (pink) contains the five amino acid residues Ile-Pro-Gln-Lys-Phe. The binding of the inhibitor compound arenicin within the active site of the enzyme might hinder the efficiency of interaction between DHFR and partner proteins and consequently reduces its activity
Arenicin Binds to Hot Spots Residues Near the Active Site of EBG
EBG folds into 8-barrel structure with both alfa and beta sheets, a typical conformation presented by the glycosyl hydrolase family. The active site is in a pocket within the protein structure 34. The pocket is energetically maintained by hydrogen-bonded interactions 35. Figure 5 shows the surface of EBG and the place where arenicin invades EBG towards the active site. We found seven hot spot residues on the EBG surface. They energetically contribute to the binding between the protein and the marine inhibitor. The main hot spot residues from EBG interacting with the inhibitor peptide are near residues belonging to the active site Figure 6. Arenicin bind to hot spot residues and changes the conformational structure of amino acids surrounding the active site of the protein, consequently inhibiting its activity.
Surface overview of EBG and arenicin interaction. Arenicin (red) binds to EBG and invades a pocket region of the protein that has amino acid residues belonging to the active site. Blue – EBG; red – arenicin; white – interface of interaction between the inhibitor peptide arenicin and EBG; green – hot spot residue.
Arenicin and the pocket region of EBG. Arenicin (red) binds to EBG and invades a pocket region of the protein that has amino acid residues belonging to the active site (pink). Green amino acids are hot spots on EBG and the yellow amino acid residue was classified as a hot spot and belongs to the active site of EBG, an important residue for inhibiting assays.
Arenicin is Surround by Hot Spots Residues in a Cleft of CYP51
CYP51 is required for ergosterol biosynthesis 36. It has been extensively investigated as the main target for azole antimicrobial drugs 373839. We found 11 hot spot residues on the C. albicans CYP51 interface of interaction with the peptide arenicin Figure 7. The necessity of finding new target to fungal diseases is real, since azoles, although having the ability to inhibit CYP51, are toxic to mammals cells and lead to severe side effects in patients 40414243. The marine peptide arenicin is held near the active site of the enzyme by the hydrophobicity in that area. Inhibition of CYP51 leads to disruption of the membrane and cell death. The marine peptide fits in a cleft of the protein, at the entrance of the active site, and is surrounded by hot spot residues, interacting with them via hydrogen bonds.
Arenicin surrounded by hot spot residues on CYP51 structure. Arenicin (red) binds to CYP51 (pink) and interacts with hot spot residues. Green amino acids are hot spots on CYP51, the conformation of the region of interaction changes, affecting the ability of the protein to perform its function.
Concluding Remarks
The number of deaths of immunocompromised patients caused by opportunistic fungi, including C. albicans, has increased significantly 44. The therapeutics used against such infections are mostly based on azoles and other drugs that develop severe side effects in those patients. Here, we showed through an in silico approach, that a single compound from a marine worm (A. marina), can bind to three C. albicans essential proteins. The interaction occurs in regions inside the active site or at least near, with amino acid residues evaluated as hot spots. Arenicin is a new promising antifugal drug. The next step is to investigate protein-protein interactions performed by DHFR, EBG and CYP51 and assess whether arenicin is able to disrupt essential interaction or not.
BongominFGagoSOladeleR ODenningD WGlobal and Multi-National Prevalence of Fungal Diseases-Estimate Precision18
de outubro de 2017;3(4)
PerfectJ RThe antifungal pipeline: a reality check. Nat Rev Drug Discov. setembro de 2017;16(9):60316RSNBrilhanteARodrigues PH dedeAlencar LPRielloG BRibeiroJ FdeOliveira JSEvidence of Fluconazole-Resistant Candida Species in Tortoises and Sea Turtles. Mycopathologia. dezembro de2015180PrigentGAït-AmmarNLevesqueEFekkarACostaJ-MSEl AnbassiEchinocandin Resistance in Candida Species Isolates from Liver Transplant Recipients. Antimicrob Agents Chemother2017LiTDingTLiJMedicinal Purposes: Bioactive Metabolites from Marine-derived Organisms. Mini Rev Med Chem. 26 de setembro de2017LeeYPhatCHongS-CStructural diversity of marine cyclic peptides and their molecular mechanisms for anticancer, antibacterial, antifungal, and other clinical applications. Peptides. setembro de 2017;95:94105Ruiz-RuizFMancera-AndradeE IHMNIqbalMarine-Derived Bioactive Peptides for Biomedical Sectors: A Review. Protein Pept Lett24210917GogineniVHamannM TMarine natural product peptides with therapeutic potential: Chemistry, biosynthesis, and pharmacology1862181196WangY-HDongH-HZhaoFWangJYanFJiangY-YThe synthesis and synergistic antifungal effects of chalcones against drug resistant Candida albicans26133098102AndräJJakovkinIGrötzingerJHechtOKrasnosdembskayaA DGoldmannTStructure and mode of action of the antimicrobial peptide arenicin410111322ChenM JShimadaTMoultonA DHarrisonMNienhuisA WIntronless human dihydrofolate reductase genes are derived from processed RNA moleculesProc Natl Acad Sci USA. dezembro de 1982;79(23):74359UrlaubGChasinL AIsolation of Chinese hamster cell mutants deficient in dihydrofolate reductase activityProc Natl Acad Sci USA. julho de 1980;77(7):421620CôtéFHahnM GOligosaccharins: structures and signal transduction. Plant Mol Biol. dezembro de 1994;26(5):1379411MarichalPGorrensJLaurijssensLVermuytenKCVan HoveLAccumulation of 3-ketosteroids induced by itraconazole in azole-resistant clinical Candida albicans isolates. Antimicrob Agents Chemother. novembro de 1999;43(11):266370FerreiraL GDosSantos RNOlivaGAndricopuloA DMolecular docking and structure-based drug design strategies20713384421GhoorahA WDevignesM-DSmaïl-TabboneMRitchieD WClassification and Exploration of 3D Protein Domain Interactions Using Kbdock. Methods Mol Biol201691105KozakovDHallD RXiaBPorterK APadhornyDYuehCThe ClusPro web server for protein-protein docking. Nat Protoc. fevereiro de 2017;12(2):25578DarnellS JLeGaultLMitchellJ CKFC Server: interactive forecasting of protein interaction hot spots. Nucleic Acids Res1ode julho de 2008;36(Web Server265269DarnellS JPageDMitchellJ CAn automated decision-tree approach to predicting protein interaction hot spotsProteins. 1ode setembro de68481323LipskaA GSieradzanA KKrupaPMozolewskaM AD’AuriaSLiwoAStudies of conformational changes of an arginine-binding protein from Thermotoga maritima in the presence and absence of ligand via molecular dynamics simulations with the coarse-grained UNRES force field21364VostrikovSoller KJHaK NGopinathTVegliaGEffects of naturally occurring arginine 14 deletion on phospholamban conformational dynamics and membrane interactions. Biochim Biophys Acta. janeiro de 2015;1848(1 Pt B):31522WenLChenYLiaoJZhengXYinZPreferential interactions between protein and arginine: effects of arginine on tertiary conformational and colloidal stability of protein solution478275361BiswasAShuklaARSKVijayanJeyakanthanJSekarKCrystal structures of an archaeal thymidylate kinase from Sulfolobus tokodaii provide insights into the role of a conserved active site Arginine residue197323649QinMSongHDaiXChenYGuoZTwo active site arginines are critical determinants of substrate binding and catalysis in MenD, a thiamine-dependent enzyme in menaquinone biosynthesis201819FordwourO BWolthersK RActive site arginine controls the stereochemistry of hydride transfer in cyclohexanone monooxygenase. Arch Biochem Biophys. 1ode outubro de 2018;659:4756KakudaKYamaguchiK-IKuwataKHondaRA valine-to-lysine substitution at position 210 induces structural conversion of prion protein into a β-sheet rich oligomer. Biochem Biophys Res Commun. 16 de outubro de2018MasbuchinA NRohmanM SPutriJ FCahyaningtyasMWidodo null. 279(Val→Phe) Polymorphism of lipoprotein-associated phospholipase A(2) resulted in changes of folding kinetics and recognition to substrate. Comput Biol Chem. dezembro de 2015;59 Pt A:MoreauMTakahashiHSariM-ABoucherJ-LSagamiIShimizuTImportance of valine 567 in substrate recognition and oxidation by neuronal nitric oxide synthase98712009MatthewsD AAldenR ABolinJ TFreerS THamlinRXuongNDihydrofolate reductase: x-ray structure of the binary complex with methotrexate19743024525OsborneM JSchnellJBenkovicS JDysonH JWrightP EBackbone dynamics in dihydrofolate reductase complexes: role of loop flexibility in the catalytic mechanism4033984659LiuJBolstadD BSmithA EPriestleyN DWrightD LAndersonA CStructure-guided development of efficacious antifungal agents targeting Candida glabrata dihydrofolate reductase. Chem Biol. 22 de setembro de 2008;15(9):9906LiuJBolstadD BSmithA EPriestleyN DWrightD LAndersonA CProbing the active site of Candida glabrata dihydrofolate reductase with high resolution crystal structures and the synthesis of new inhibitors. Chem Biol Drug Des. janeiro de 2009;73(1):6274PaulsenJ LBendelS DAndersonA CCrystal structures of Candida albicans dihydrofolate reductase bound to propargyl-linked antifolates reveal the flexibility of active site loop residues critical for ligand potency and selectivity. Chem Biol Drug Des. outubro de 2011;78(4):50512CutfieldS MDaviesG JMurshudovGAndersonB FMoodyP CSullivanP AThe structure of the exo-beta-(1,3)-glucanase from Candida albicans in native and bound forms: relationship between a pocket and groove in family 5 glycosyl hydrolases294377183SakonJAdneyW SHimmelM EThomasS RKarplusP ACrystal structure of thermostable family 5 endocellulase E1 from Acidothermus cellulolyticus in complex with cellotetraose. Biochemistry. 20 de agosto de 1996;35(33):1064860KellyS LLambD CJacksonC JWarrilowA GKellyD EThe biodiversity of microbial cytochromes P450. Adv Microb Physiol200313186AGSWarrilowMartelC MParkerJ EMeloNLambD CNesW DAzole binding properties of Candida albicans sterol 14-alpha demethylase (CaCYP51). Antimicrob Agents Chemother. outubro de 2010;54(10):423545ParkerJ EMerkammMManningN JPomponDKellyS LKellyD EDifferential azole antifungal efficacies contrasted using a Saccharomyces cerevisiae strain humanized for sterol 14 alpha-demethylase at the homologous locus. Antimicrob Agents Chemother. outubro de 2008;52(10):3597603LambD CKellyD EWatermanM RStromstedtMRozmanDKellyS LCharacteristics of the heterologously expressed human lanosterol 14 alpha-demethylase (other names: P45014DM, CYP51, P45051) and inhibition of the purified human and Candida albicans CYP51 with azole antifungal agents15975563KorashyH MShayeganpourABrocksD RAOSEl-KadiInduction of cytochrome P450 1A1 by ketoconazole and itraconazole but not fluconazole in murine and human hepatoma cell lines. Toxicol Sci. maio de 2007;97(1):3243TaxvigCVinggaardA MHassUAxelstadMMetzdorffSNellemannCEndocrine-disrupting properties in vivo of widely used azole fungicides3121707PontAWilliamsP LAzharSReitzR EBochraCSmithE RKetoconazole blocks testosterone synthesis. Arch Intern Med. novembro de 1982;142(12):213740PontAWilliamsP LLooseD SFeldmanDReitzR EBochraCKetoconazole blocks adrenal steroid synthesis. Ann Intern Med. setembro de 1982;97(3):3702LowC-YRotsteinCEmerging fungal infections in immunocompromised patients. F1000 Med Rep [Internet]. 1ode julho de 2011 [citado 29 de outubro de 2018];3. Disponível em: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3155160/