1
|
de Oliveira ALN, Lacerda MT, Ramos MJ, Fernandes PA. Viper Venom Phospholipase A2 Database: The Structural and Functional Anatomy of a Primary Toxin in Envenomation. Toxins (Basel) 2024; 16:71. [PMID: 38393149 PMCID: PMC10893444 DOI: 10.3390/toxins16020071] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2023] [Revised: 01/04/2024] [Accepted: 01/13/2024] [Indexed: 02/25/2024] Open
Abstract
Viper venom phospholipase A2 enzymes (vvPLA2s) and phospholipase A2-like (PLA2-like) proteins are two of the principal toxins in viper venom that are responsible for the severe myotoxic and neurotoxic effects caused by snakebite envenoming, among other pathologies. As snakebite envenoming is the deadliest neglected tropical disease, a complete understanding of these proteins' properties and their mechanisms of action is urgently needed. Therefore, we created a database comprising information on the holo-form, cofactor-bound 3D structure of 217 vvPLA2 and PLA2-like proteins in their physiologic environment, as well as 79 membrane-bound viper species from 24 genera, which we have made available to the scientific community to accelerate the development of new anti-snakebite drugs. In addition, the analysis of the sequenced, 3D structure of the database proteins reveals essential aspects of the anatomy of the proteins, their toxicity mechanisms, and the conserved binding site areas that may anchor universal interspecific inhibitors. Moreover, it pinpoints hypotheses for the molecular origin of the myotoxicity of the PLA2-like proteins. Altogether, this study provides an understanding of the diversity of these toxins and how they are conserved, and it indicates how to develop broad, interspecies, efficient small-molecule inhibitors to target the toxin's many mechanisms of action.
Collapse
Affiliation(s)
| | | | | | - Pedro A. Fernandes
- Requimte-Faculty of Sciences, University of Porto, Rua do Campo Alegre s/n, 4169-000 Porto, Portugal; (A.L.N.d.O.); (M.T.L.); (M.J.R.)
| |
Collapse
|
2
|
S S, V S, I JM, P VM, P LK, Nair AS, R SP, Oommen OV. In silico screening of the phytochemicals present in Clitoria ternatea L. as the inhibitors of snake venom phospholipase A 2 (PLA 2). J Biomol Struct Dyn 2023; 41:7874-7883. [PMID: 36153001 DOI: 10.1080/07391102.2022.2126889] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2021] [Accepted: 09/15/2022] [Indexed: 10/14/2022]
Abstract
Millions of people suffer from snake bite envenomation, and its management is a challenge, even today. Medicinal plants have attracted the researcher's attention for their outstanding advantages in treating many diseases, including snake venom poisoning. Clitoria ternatea L, is a plant popularly known for its various pharmacological effects especially, anti-snake venom property. However, the molecular mechanism behind this is poorly understood. It is reported that snake venom PLA2 is an extensively studied toxic factor. This study is meant to screen the compound's capability to act as inhibitors of the Daboia russelli snake venom PLA2 through molecular docking and dynamics studies. Our results show that among the 27 compounds taken for the study, only Kaempferol showed good interaction profile with the conserved catalytic active site residues, His48 and Asp49. The pharmacophore features of the compound also demonstrate its exact fitting at the binding pocket. Further RMSD, RMSF, Rg, and hydrogen bond analysis confirmed the stable binding of Kaempferol with PLA2 through molecular dynamic simulations for 100 ns. In addition, the MM/PBSA binding free energy calculation of the complex was also affirming the docking results. The binding free energy (BFE) of Kaempferolis better than the reference compound. ADME and Lipinski's rule of five reveals its drug like properties.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Suveena S
- Centre for Venom Informatics, Department of Computational Biology & Bioinformatics, University of Kerala, Trivandrum, Kerala, India
| | - Saraswathy V
- Centre for Venom Informatics, Department of Computational Biology & Bioinformatics, University of Kerala, Trivandrum, Kerala, India
| | - Junaida M I
- Centre for Venom Informatics, Department of Computational Biology & Bioinformatics, University of Kerala, Trivandrum, Kerala, India
| | - Vinod M P
- Centre for Venom Informatics, Department of Computational Biology & Bioinformatics, University of Kerala, Trivandrum, Kerala, India
| | - Laladhas K P
- Department of Zoology, St.Stephen's College, Kollam, Kerala, India
| | - Achuthsankar S Nair
- Centre for Venom Informatics, Department of Computational Biology & Bioinformatics, University of Kerala, Trivandrum, Kerala, India
| | - Sudhakaran P R
- Centre for Venom Informatics, Department of Computational Biology & Bioinformatics, University of Kerala, Trivandrum, Kerala, India
| | - Oommen V Oommen
- Centre for Venom Informatics, Department of Computational Biology & Bioinformatics, University of Kerala, Trivandrum, Kerala, India
| |
Collapse
|
3
|
Sharifpour S, Fakhraee S, Behjatmanesh-Ardakani R. Insights into the mechanism of inhibition of phospholipase A2 by resveratrol: An extensive molecular dynamics simulation and binding free energy calculation. J Mol Graph Model 2020; 100:107649. [PMID: 32739638 DOI: 10.1016/j.jmgm.2020.107649] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2020] [Revised: 06/02/2020] [Accepted: 06/03/2020] [Indexed: 10/23/2022]
Abstract
Phospholipase A2 (PLA2) is one of the enzymes involved in the development of cardiovascular diseases, vascular inflammation, risk of heart attacks, and strokes. This enzyme is responsible for catalyzing the hydrolytic cleavage of ester bonds of phospholipids in the biological pathway of inflammation. To prevent the undesired hydrolysis of phospholipids, the catalytic activity of PLA2 needs to be blocked. Resveratrol is a plant-derived polyphenol inhibitor, proven to have anti-inflammatory properties. However, there is still substantial ambiguity about its inhibitory function. The present study uncovers a detailed molecular mechanism behind the resveratrol action in inhibition of PLA2, by applying and comparing two 200-ns molecular dynamics simulations. The results of structural analyses revealed that the binding of resveratrol to PLA2 reduces the content of β-sheets and increases a 5-helix to PLA2 structure, producing more folding and stability in protein. In the active site, the resveratrol is placed between the N-terminal α-helix and the newly formed 5-helix through the hydrophobic interactions with ILE19 and LEU3 residues, as well as the hydrogen bond interactions. These interactions play the role of a network at the entrance of the enzyme active site and prevent the penetration of water molecules into the PLA2 cavity. A high occupancy hydrogen bonding has been identified between SER23 of the protein and hydroxyl group of resveratrol. Furthermore, the estimation of binding free energy verified the binding affinity of resveratrol is thermodynamically sufficient to be stably bounded to PLA2. It also proved that the van der Waals interactions, particularly hydrophobic interactions, have the most significant role in PLA2-resveratrol binding and stability. Overall, our results provide useful information on the stepwise mechanism of the inhibition of PLA2 enzyme by resveratrol, as a target for improving the pharmacological applications.
Collapse
Affiliation(s)
- Sajedeh Sharifpour
- Department of Chemistry, Payame Noor University, 19395-3697, Tehran, Iran
| | - Sara Fakhraee
- Department of Chemistry, Payame Noor University, 19395-3697, Tehran, Iran.
| | | |
Collapse
|
4
|
Vindas J, Carrera Y, Lomonte B, Gutiérrez JM, Calvete JJ, Sanz L, Fernández J. A novel pentameric phospholipase A2 myotoxin (PophPLA2) from the venom of the pit viper Porthidium ophryomegas. Int J Biol Macromol 2018; 118:1-8. [DOI: 10.1016/j.ijbiomac.2018.06.028] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2018] [Revised: 06/03/2018] [Accepted: 06/06/2018] [Indexed: 01/07/2023]
|
5
|
Sivaramakrishnan V, Ilamathi M, Ghosh KS, Sathish S, Gowda TV, Vishwanath BS, Rangappa KS, Dhananjaya BL. Virtual analysis of structurally diverse synthetic analogs as inhibitors of snake venom secretory phospholipase A2. J Mol Recognit 2015. [PMID: 26218369 DOI: 10.1002/jmr.2492] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
Due to the toxic pathophysiological role of snake venom phospholipase A2 (PLA2 ), its compelling limitations to anti-venom therapy in humans and the need for alternative therapy foster considerable pharmacological interest towards search of PLA2 specific inhibitors. In this study, an integrated approach involving homology modeling, molecular dynamics and molecular docking studies on VRV-PL-V (Vipera russellii venom phospholipase A2 fraction-V) belonging to Group II-B secretory PLA2 from Daboia russelli pulchella is carried out in order to study the structure-based inhibitor design. The accuracy of the model was validated using multiple computational approaches. The molecular docking study of this protein was undertaken using different classes of experimentally proven, structurally diverse synthetic inhibitors of secretory PLA2 whose selection is based on IC50 value that ranges from 25 μM to 100 μM. Estimation of protein-ligand contacts by docking analysis sheds light on the importance of His 47 and Asp 48 within the VRV-PL-V binding pocket as key residue for hydrogen bond interaction with ligands. Our virtual analysis revealed that compounds with different scaffold binds to the same active site region. ADME analysis was also further performed to filter and identify the best potential specific inhibitor against VRV-PL-V. Additionally, the e-pharmacophore was generated for the best potential specific inhibitor against VRV-PL-V and reported here. The present study should therefore play a guiding role in the experimental design of VRV-PL-V inhibitors that may provide better therapeutic molecular models for PLA2 recognition and anti-ophidian activity.
Collapse
Affiliation(s)
- V Sivaramakrishnan
- Cardiomyocyte Toxicity and Oncology Research Lab, Department of Bioinformatics, School of Chemical and Biotechnology, SASTRA University, Thanjavur, 613401, India
| | - M Ilamathi
- Cardiomyocyte Toxicity and Oncology Research Lab, Department of Bioinformatics, School of Chemical and Biotechnology, SASTRA University, Thanjavur, 613401, India
| | - K S Ghosh
- Department of Chemistry, National Institute of Technology (NIT), Hamirpur, 177 005, India
| | - S Sathish
- Department of Studies in Biochemistry, University of Mysore, Mysore, 570006, India
| | - T V Gowda
- Department of Studies in Biochemistry, University of Mysore, Mysore, 570006, India
| | - B S Vishwanath
- Department of Studies in Biochemistry, University of Mysore, Mysore, 570006, India
| | - K S Rangappa
- Department of Chemistry, University of Mysore, Mysore, 570006, India
| | - B L Dhananjaya
- Toxinology/Toxicology and Drug Discovery Unit, Center for Emerging Technologies, Jain University, Jakkasandra post, Ramanagara, 562112, India
| |
Collapse
|
6
|
Shukla PK, Gautam L, Sinha M, Kaur P, Sharma S, Singh TP. Structures and binding studies of the complexes of phospholipase A2 with five inhibitors. BIOCHIMICA ET BIOPHYSICA ACTA 2015; 1854:269-77. [PMID: 25541253 DOI: 10.1016/j.bbapap.2014.12.017] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/12/2014] [Revised: 12/01/2014] [Accepted: 12/16/2014] [Indexed: 11/25/2022]
Abstract
Phospholipase A2 (PLA2) catalyzes the hydrolysis of phospholipids into arachidonic acid and lysophospholipids. Arachidonic acid is used as a substrate in the next step of the multistep pathway leading to the production of eicosanoids. The eicosanoids, in extremely low concentrations, are required in a number of physiological processes. However, the increase in their concentrations above the essential physiological requirements leads to various inflammatory conditions. In order to prevent the unwanted rise in the concentrations of eicosanoids, the actions of PLA2 and other enzymes of the pathway need to be blocked. We report here the structures of five complexes of group IIA PLA2 from Daboia russelli pulchella with tightly binding inhibitors, (i) p-coumaric acid, (ii) resveratrol, (iii) spermidine, (iv) corticosterone and (v) gramine derivative. The binding studies using fluorescence spectroscopy and surface plasmon resonance techniques for the interactions of PLA2 with the above five compounds showed high binding affinities with values of dissociation constants (KD) ranging from 3.7×10(-8) M to 2.1×10(-9) M. The structure determinations of the complexes of PLA2 with the above five compounds showed that all the compounds bound to PLA2 in the substrate binding cleft. The protein residues that contributed to the interactions with these compounds included Leu2, Leu3, Phe5, Gly6, Ile9, Ala18, Ile19, Trp22, Ser23, Cys29, Gly30, Cys45, His48, Asp49 and Phe106. The positions of side chains of several residues including Leu2, Leu3, Ile19, Trp31, Lys69, Ser70 and Arg72 got significantly shifted while the positions of active site residues, His48, Asp49, Tyr52 and Asp99 were unperturbed.
Collapse
Affiliation(s)
- Prakash Kumar Shukla
- Department of Biophysics, All India Institute of Medical Sciences, New Delhi, India
| | - Lovely Gautam
- Department of Biophysics, All India Institute of Medical Sciences, New Delhi, India
| | - Mau Sinha
- Department of Biophysics, All India Institute of Medical Sciences, New Delhi, India
| | - Punit Kaur
- Department of Biophysics, All India Institute of Medical Sciences, New Delhi, India
| | - Sujata Sharma
- Department of Biophysics, All India Institute of Medical Sciences, New Delhi, India
| | - Tej P Singh
- Department of Biophysics, All India Institute of Medical Sciences, New Delhi, India.
| |
Collapse
|
7
|
Kang TS, Georgieva D, Genov N, Murakami MT, Sinha M, Kumar RP, Kaur P, Kumar S, Dey S, Sharma S, Vrielink A, Betzel C, Takeda S, Arni RK, Singh TP, Kini RM. Enzymatic toxins from snake venom: structural characterization and mechanism of catalysis. FEBS J 2011; 278:4544-76. [PMID: 21470368 DOI: 10.1111/j.1742-4658.2011.08115.x] [Citation(s) in RCA: 185] [Impact Index Per Article: 14.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
Abstract
Snake venoms are cocktails of enzymes and non-enzymatic proteins used for both the immobilization and digestion of prey. The most common snake venom enzymes include acetylcholinesterases, l-amino acid oxidases, serine proteinases, metalloproteinases and phospholipases A(2) . Higher catalytic efficiency, thermal stability and resistance to proteolysis make these enzymes attractive models for biochemists, enzymologists and structural biologists. Here, we review the structures of these enzymes and describe their structure-based mechanisms of catalysis and inhibition. Some of the enzymes exist as protein complexes in the venom. Thus we also discuss the functional role of non-enzymatic subunits and the pharmacological effects of such protein complexes. The structures of inhibitor-enzyme complexes provide ideal platforms for the design of potent inhibitors which are useful in the development of prototypes and lead compounds with potential therapeutic applications.
Collapse
Affiliation(s)
- Tse Siang Kang
- Department of Pharmacy, National University of Singapore, Singapore
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
8
|
Structural analysis of a group III Glu62-phospholipase A2 from the scorpion, Mesobuthus tamulus: Targeting and reversible inhibition by native peptides. Int J Biol Macromol 2011; 48:423-31. [PMID: 21238479 DOI: 10.1016/j.ijbiomac.2011.01.004] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2010] [Revised: 12/22/2010] [Accepted: 01/04/2011] [Indexed: 11/22/2022]
Abstract
Group III phospholipase A(2) enzyme transcript from the Mesobuthus tamulus (Indian red scorpion) codes for three distinct products that include a large enzymatic subunit, a pentameric peptide and a small non-enzymatic subunit. The structures of these two subunits were modeled based on their sequence identity to bee venom PLA(2) and the partial sequence of MU2 adaptin subunit of AP2 clathrin adaptor, respectively. The enzymatic subunit comprises of three helices, the calcium binding loop and a substrate binding hydrophobic channel where the structure is stabilized by four disulfide bonds. The active site of the enzyme shows a catalytic histidine residue. Interestingly, there is a conservative mutation of the conserved aspartic acid, a classical participant of catalysis in this enzyme family, to glutamic acid. However, the side chain oxygen atoms of this glutamate are oriented away from the catalytic histidine implicating the non-participation of this residue in stabilizing the tautomeric conformation of the histidine. The acidic non-enzymatic subunit comprises of extensive hydrophobic residues with a conformation of an anti-parallel β-sheets making it ideal for tissue specific targeting. The native pentapeptide with the sequence Alanine-Arginine-Serine-Alanine-Arginine was docked to the enzymatic subunit. The peptide ligand occupies the hydrophobic cavity and makes a plethora of interactions with the residues in the channel, including a hydrogen bond with the crucial catalytic histidine and coordinate bond with the calcium ion. This ligand has a binding constant (K(D)) of 1.5μM. This makes the ligand a potential reversible inhibitor, ideal to prevent the enzyme from interacting with non-specific molecules enroute to the target. The enzyme-ligand complex also provides a model to understand the stereochemistry required for the design of more potent drug molecules against such enzyme drug targets.
Collapse
|
9
|
Hariprasad G, Kumar M, Kaur P, Singh TP, Kumar RP. Human group III PLA2 as a drug target: structural analysis and inhibitor binding studies. Int J Biol Macromol 2010; 47:496-501. [PMID: 20654644 DOI: 10.1016/j.ijbiomac.2010.07.004] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2010] [Revised: 07/13/2010] [Accepted: 07/13/2010] [Indexed: 11/24/2022]
Abstract
Group III phospholipase A(2) is a known mediator of inflammation, atherosclerosis and cancer in mammals. This enzyme, therefore, is a potential drug target. The availability of the human group III phospholipase A(2) (hIIIPLA(2)) amino acid sequence offers an opportunity to study its structural features by modeling. The monomeric hIII PLA(2) model is based on the 44% identity it has with the bee venom PLA(2), the only known representative structure of this group. The overall structure comprises of three α-helices, a β-wing and the calcium binding loop which is present at the N-terminus of the enzyme. However, the unique structural features of hIIIPLA(2) in comparison to the other well known group I/II PLA(2)s are: (1) the replacement of the 'conserved' tyrosine residue by phenylalanine at position 87 in the active site; (2) a decrease in the volume of the substrate binding hydrophobic channel and (3) presence of a C-terminal extension which has a close proximity to the third helix. Docking studies of the enzyme with small molecules gives a detailed insight into the participating residues of the enzyme and also the possible type of interactions with the drug molecules. The ligand molecules have binding affinities predicted to range from micromolar to nanomolar range, thereby making them either potential lead molecules or potent drugs. This analysis paves the way for possible therapeutic applications in pathological states caused by this enzyme.
Collapse
Affiliation(s)
- Gururao Hariprasad
- Department of Biophysics, All India Institute of Medical Sciences, New Delhi 110029, India.
| | | | | | | | | |
Collapse
|
10
|
Nonantibiotic properties of tetracyclines: structural basis for inhibition of secretory phospholipase A2. J Mol Biol 2010; 398:83-96. [PMID: 20211188 DOI: 10.1016/j.jmb.2010.02.049] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2009] [Revised: 02/26/2010] [Accepted: 02/28/2010] [Indexed: 11/21/2022]
Abstract
Secretory phospholipase A(2) is involved in inflammatory processes and was previously shown to be inhibited by lipophilic tetracyclines such as minocycline (minoTc) and doxycycline. Lipophilic tetracyclines might be a new lead compound for the design of specific inhibitors of secretory phospholipase A(2), which play a crucial role in inflammatory processes. Our X-ray crystal structure analysis at 1.65 A resolution of the minoTc complex of phospholipase A(2) (PLA(2)) of the Indian cobra (Naja naja naja) is the first example of nonantibiotic tetracycline interactions with a protein. MinoTc interferes with the conformation of the active-site Ca(2+)-binding loop, preventing Ca(2)(+) binding, and shields the active site from substrate entrance, resulting in inhibition of the enzyme. MinoTc binding to PLA(2) is dominated by hydrophobic interactions quite different from antibiotic recognition of tetracyclines by proteins or the ribosome. The affinity of minoTc for PLA(2) was determined by surface plasmon resonance, resulting in a dissociation constant K(d)=1.8 x 10(-)(4) M.
Collapse
|