Crystal structure of cobra-venom phospholipase A2 in a complex with a transition-state analogue

Viper Venom Phospholipase A2 Database: The Structural and Functional Anatomy of a Primary Toxin in Envenomation

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.

In Silico Evaluation of Quercetin Methylated Derivatives on the Interaction with Secretory Phospholipases A2 from Crotalus durissus terrificus and Bothrops jararacussu

Quercetin derivatives have already shown their anti-inflammatory potential, inhibiting essential enzymes involved in this process. Among diverse pro-inflammatory toxins from snake venoms, phospholipase A2 is one of the most abundant in some species, such as Crotalus durissus terrificus and Bothrops jararacussu from the Viperidae family. These enzymes can induce the inflammatory process through hydrolysis at the sn-2 position of glycerophospholipids. Hence, elucidating the main residues involved in the biological effects of these macromolecules can help to identify potential compounds with inhibitory activity. In silico tools were used in this study to evaluate the potential of quercetin methylated derivatives in the inhibition of bothropstoxin I (BthTX-I) and II (BthTX-II) from Bothrops jararacussu and phospholipase A2 from Crotalus durissus terrificus. The use of a transitional analogous and two classical inhibitors of phospholipase A2 guided this work to find the role of residues involved in the phospholipid anchoring and the subsequent development of the inflammatory process. First, main cavities were studied, revealing the best regions to be inhibited by a compound. Focusing on these regions, molecular docking assays were made to show main interactions between each compound. Results reveal that analogue and inhibitors, Varespladib (Var) and p-bromophenacyl bromide (BPB), guided quercetins derivatives analysis, revealing that Leu2, Phe5, Tyr28, glycine in the calcium-binding loop, His48, Asp49 of BthTX-II and Cdtspla2 were the main residues to be inhibited. 3MQ exhibited great interaction with the active site, similar to Var results, while Q anchored better in the BthTX-II active site. However, strong interactions in the C-terminal region, highlighting His120, seem to be crucial to decreasing contacts with phospholipid and BthTX-II. Hence, quercetin derivatives anchor differently with each toxin and further in vitro and in vivo studies are essential to elucidate these data.

Catalytically Active Snake Venom PLA2 Enzymes: An Overview of Its Elusive Mechanisms of Reaction

Snake venom-secreted phospholipase A₂ (svPLA₂) enzymes, both catalytically active and inactive, are a central component in envenoming. These are responsible for disrupting the cell membrane's integrity, inducing a wide range of pharmacological effects, such as the necrosis of the bitten limb, cardiorespiratory arrest, edema, and anticoagulation. Although extensively characterized, the reaction mechanisms of enzymatic svPLA₂ are still to be thoroughly understood. This review presents and analyses the most plausible reaction mechanisms for svPLA_2, such as the "single-water mechanism" or the "assisted-water mechanism" initially proposed for the homologous human PLA₂. All of the mechanistic possibilities are characterized by a highly conserved Asp/His/water triad and a Ca²⁺ cofactor. The extraordinary increase in activity induced by binding to a lipid-water interface, known as "interfacial activation," critical for the PLA₂s activity, is also discussed. Finally, a potential catalytic mechanism for the postulated noncatalytic PLA₂-like proteins is anticipated.

Comparative phylogenomic and structural analysis of canonical secretory PLA2 and novel PLA2-like family in plants

Plant secretory phospholipase A₂ (sPLA₂) is a family of lipolytic enzymes involved in the sn-2 hydrolysis of phospholipid carboxyester bonds, characterized by the presence of a conserved PA2c domain. PLA₂ produces free fatty acids and lysophospholipids, which regulate several physiological functions, including lipid metabolism, plant growth and development, signal transduction, and response to various environmental stresses. In the present work, we have performed a comparative analysis of PA2c domain-containing genes across plants, focusing on gene distribution, phylogenetic analysis, tissue-specific expression, and homology modeling. Our data revealed the widespread occurrence of multiple sPLA₂ in most land plants and documented single sPLA₂ in multiple algal groups, indicating an ancestral origin of sPLA₂. We described a novel PA2c-containing gene family present in all plant lineages and lacking secretory peptide, which we termed PLA₂-like. Phylogenetic analysis revealed two independent clades in canonical sPLA₂ genes referred to as α and β clades, whereas PLA₂-like genes clustered independently as a third clade. Further, we have explored clade-specific gene expressions showing that while all three clades were expressed in vegetative and reproductive tissues, only sPLA₂-β and PLA₂-like members were expressed in the pollen and pollen tube. To get insight into the conservation of the gene regulatory network of sPLA₂ and PLA₂-like genes, we have analyzed the occurrence of various cis-acting promoter elements across the plant kingdom. The comparative 3D structure analysis revealed conserved and unique features within the PA2c domain for the three clades. Overall, this study will help to understand the evolutionary significance of the PA2c family and lay the foundation for future sPLA₂ and PLA₂-like characterization in plants.

Structural Dynamics and Activity of B19V VP1u during the pHs of Cell Entry and Endosomal Trafficking

Parvovirus B19 (B19V) is a human pathogen that is the causative agent of fifth disease in children. It is also known to cause hydrops in fetuses, anemia in AIDS patients, and transient aplastic crisis in patients with sickle cell disease. The unique N-terminus of Viral Protein 1 (VP1u) of parvoviruses, including B19V, exhibits phospholipase A2 (PLA2) activity, which is required for endosomal escape. Presented is the structural dynamics of B19V VP1u under conditions that mimic the pHs of cell entry and endosomal trafficking to the nucleus. Using circular dichroism spectroscopy, the receptor-binding domain of B19V VP1u is shown to exhibit an α-helical fold, whereas the PLA2 domain exhibits a probable molten globule state, both of which are pH invariant. Differential scanning calorimetry performed at endosomal pHs shows that the melting temperature (Tm) of VP1u PLA2 domain is tuned to body temperature (37 °C) at pH 7.4. In addition, PLA2 assays performed at temperatures ranging from 25-45 °C show both a temperature and pH-dependent change in activity. We hypothesize that VP1u PLA2 domain differences in Tm at differing pHs have enabled the virus to "switch on/off" the phospholipase activity during capsid trafficking. Furthermore, we propose the environment of the early endosome as the optimal condition for endosomal escape leading to B19V infection.

Characterization of New Allergens from the Venom of the European Paper Wasp Polistes dominula

Discriminating Polistes dominula and Vespula spp. venom allergy is of growing importance worldwide, as systemic reactions to either species’ sting can lead to severe outcomes. Administering the correct allergen-specific immunotherapy is therefore a prerequisite to ensure the safety and health of venom-allergic patients. Component-resolved diagnostics of Hymenoptera venom allergy might be improved by adding additional allergens to the diagnostic allergen panel. Therefore, three potential new allergens from P. dominula venom—immune responsive protein 30 (IRP30), vascular endothelial growth factor C (VEGF C) and phospholipase A2 (PLA2)—were cloned, recombinantly produced and biochemically characterized. Sera sIgE titers of Hymenoptera venom-allergic patients were measured in vitro to assess the allergenicity and potential cross-reactivity of the venom proteins. IRP30 and VEGF C were classified as minor allergens, as sensitization rates lay around 20–40%. About 50% of P. dominula venom-allergic patients had measurable sIgE titers directed against PLA2 from P. dominula venom. Interestingly, PLA2 was unable to activate basophils of allergic patients, questioning its role in the context of clinically relevant sensitization. Although the obtained results hint to a questionable benefit of the characterized P. dominula venom proteins for improved diagnosis of venom-allergic patients, they can contribute to a deeper understanding of the molecular mechanisms of Hymenoptera venoms and to the identification of factors that determine the allergenic potential of proteins.

Mann T, Bastard K, F. Scott K, Church WB. Structural and Functional Aspects of Targeting the Secreted Human Group IIA Phospholipase A2

Human group IIA secretory phospholipase A₂ (hGIIA) promotes the proliferation of cancer cells, making it a compelling therapeutic target, but it is also significant in other inflammatory conditions. Consequently, suitable inhibitors of hGIIA have always been sought. The activation of phospholipases A₂ and the catalysis of glycerophospholipid substrates generally leads to the release of fatty acids such as arachidonic acid (AA) and lysophospholipid, which are then converted to mediator compounds, including prostaglandins, leukotrienes, and the platelet-activating factor. However, this ability of hGIIA to provide AA is not a complete explanation of its biological role in inflammation, as it has now been shown that it also exerts proinflammatory effects by a catalysis-independent mechanism. This mechanism is likely to be highly dependent on key specific molecular interactions, and the full mechanistic descriptions of this remain elusive. The current candidates for the protein partners that may mediate this catalysis-independent mechanism are also introduced in this review. A key discovery has been that selective inhibition of the catalysis-independent activity of hGIIA is achieved with cyclised derivatives of a pentapeptide, FLSYK, derived from the primary sequence of hGIIA. The effects of hGIIA on cell function appear to vary depending on the pathology studied, and so its mechanism of action is complex and context-dependent. This review is comprehensive and covers the most recent developments in the understanding of the many facets of hGIIA function and inhibition and the insight they provide into their clinical application for disease treatment. A cyclic analogue of FLSYK, c2, the most potent analogue known, has now been taken into clinical trials targeting advanced prostate cancer.

Phospholipase A2 way to hydrolysis: Dint formation, hydrophobic mismatch, and lipid exclusion

Phospholipase A2 (PLA2) exerts a wide range of biological effects and attracts a lot of attention of researchers. Two sites are involved in manifestation of PLA2 enzymatic activity: catalytic site responsible for substrate binding and fatty acid cleavage from the sn-2 position of a glycerophospholipid, and interface binding site (IBS) responsible for the protein binding to lipid membrane. IBS is formed by positively charged and hydrophobic amino acids on the outer surface of the protein molecule. Understanding the mechanism of PLA2 interaction with the lipid membrane is the most challenging step in biochemistry of this enzyme. We used a combination of experimental and computer simulation techniques to clarify molecular details of bee venom PLA2 interaction with lipid bilayers formed by palmitoyloleoylphosphatidylcholine or dipalmitoylphosphatidylcholine. We found that after initial enzyme contact with the membrane, a network of hydrogen bonds was formed. This led to deformation of the interacting leaflet and dint formation. The bilayer response to the deformation depended on its phase state. In a gel-phase bilayer, diffusion of lipids is restricted therefore chain melting occurred in both leaflets of the bilayer. In the case of a fluid-phase bilayer, lateral diffusion is possible, and lipid polar head groups were excluded from the contact area. As a result, the bilayer became thinner and a large hydrophobic area was formed. We assume that relative ability of a bilayer to come through lipid redistribution process defines the rate of initial stages of the catalysis.

An Integrated Proteomic and Transcriptomic Analysis Reveals the Venom Complexity of the Bullet Ant Paraponera clavata

A critical hurdle in ant venom proteomic investigations is the lack of databases to comprehensively and specifically identify the sequence and function of venom proteins and peptides. To resolve this, we used venom gland transcriptomics to generate a sequence database that was used to assign the tandem mass spectrometry (MS) fragmentation spectra of venom peptides and proteins to specific transcripts. This was performed alongside a shotgun liquid chromatography-mass spectrometry (LC-MS/MS) analysis of the venom to confirm that these assigned transcripts were expressed as proteins. Through the combined transcriptomic and proteomic investigation of Paraponera clavata venom, we identified four times the number of proteins previously identified using 2D-PAGE alone. In addition to this, by mining the transcriptomic data, we identified several novel peptide sequences for future pharmacological investigations, some of which conform with inhibitor cysteine knot motifs. These types of peptides have the potential to be developed into pharmaceutical or bioinsecticide peptides.

Water-Intake and Water-Molecule Paths to the Active Site of Secretory Phospholipase A2 Studied Using MD Simulations and the Tracking Tool AQUA-DUCT

Secretory phospholipases A₂ (sPLA₂s) are a subclass of enzymes that catalyze the hydrolysis at the sn-2 position of glycerophospholipids, producing free fatty acids and lysophospholipids. In this study, different phospholipids with structural modifications close to the scissile sn-2 ester bond were studied to determine the effect of the structural changes on the formation of the Michaelis-Menten complex and the water entry/exit pathways using molecular dynamics simulations and the computational tracking tool AQUA-DUCT. Structural modifications include methylation, dehydrogenation, and polarization close to the sn-2 scissile bond. We found that all water molecules reaching the active site of sPLA₂-IIA pass by the aromatic residues Phe⁵ and Tyr⁵¹ and enter the active site through an active-site cleft. The relative amount of water available for the enzymatic reaction of the different phospholipid-sPLA₂ complexes was determined together with the distance between key atoms in the catalytic machinery. The results showed that (Z)-unsaturated phospholipid is a good substrate for sPLA₂-IIA. The computational results are in good agreement with previously reported experimental data on the ability of sPLA₂-IIA to hydrolyze liposomes made from the different phospholipids, and the results provide new insights into the necessary active-site solvation of the Michaelis-Menten complex and can pave the road for rational design in engineering applications.

Neutralization of a bothropic PLA2-like protein by caftaric acid, a novel potent inhibitor of ophidian myotoxicity

Envenoming by snakebite is an important global health issue that has received little attention, leading the World Health Organization to naming it as neglected tropical disease. Several snakebites present serious local symptoms manifested on victims that may not be efficiently neutralized by serum therapy. Phospholipase A₂-like (PLA₂-like) toxins are present in Viperidae venoms and are responsible for local myotoxic activity. Herein, we investigated the association between BthTX-I toxin and caftaric acid (CFT), a molecule present in plants. CFT neutralized neuromuscular blocking and muscle-damaging activities promoted by BthTX-I. Calorimetric and light-scattering assays demonstrated that CFT inhibitor interacted with dimeric BthTX-I. Bioinformatics simulations indicated that CFT inhibitor binds to the toxin's hydrophobic channel (HCh). According to the current myotoxic mechanism, three different regions of PLA₂-like toxins have specific tasks: protein allosteric activation (HCh), membrane dockage (MDoS), and membrane rupture (MDiS). We propose CFT inhibitor interferes with the allosteric activation, which is related to the conformation change leading to the exposure/alignment of MDoS/MDiS region. This is the first report of a PLA₂-like toxin fully inhibited by a compound that interacts only with its HCh region. Thus, CFT is a novel candidate to complement serum therapy and improve the treatment of snakebite.

Intravascular hemolysis induced by phospholipases A2 from the venom of the Eastern coral snake, Micrurus fulvius: Functional profiles of hemolytic and non-hemolytic isoforms

A unique feature of the venom of Micrurus fulvius (Eastern coral snake) is its ability to induce severe intravascular hemolysis in particular species, such as dogs or mice. This effect was previously shown to be induced by distinct phospholipase A₂ (PLA₂) isoforms which cause direct hemolysis in vitro, an uncommon finding for such enzymes. The functional profiles of PLA₂-17, a direct hemolytic enzyme, and PLA₂-12, a co-existing venom isoform lacking such effect, were compared. The enzymes differed not only in their ability to cause intravascular hemolysis: PLA₂-17 additionally displayed lethal, myotoxic, and anticoagulant actions, whereas PLA₂-12 lacked these effects. PLA₂-12 was much more active in hydrolyzing a monodisperse synthetic substrate than PLA₂-17, but the catalytic activity of latter was notably higher on a micellar substrate, or towards pure phospholipid artificial monolayers under controlled lateral pressures. Interestingly, PLA₂-17 could hydrolyze substrate at a pressure of 20 mN m^-1, in contrast to PLA₂-12 or the non-toxic pancreatic PLA₂. This suggests important differences in the monolayer penetrating power, which could be related to differences in toxicity. Comparative examination of primary structures and predicted three-dimensional folding of PLA₂-12 and PLA₂-17, revealed that differences concentrate in their N-terminal and central regions, leading to variations of the surface properties at the membrane interacting interface. PLA₂-17 presents a less basic interfacial surface than PLA₂-12, but more bulky aromatic residues, which could be associated to its higher membrane-penetrating strength. Altogether, these structural and functional comparative observations suggest that the ability of PLA₂s to penetrate substrate interfaces could be a major determinant of toxicity, perhaps more important than protein surface charge.

A two-helix motif positions the lysophosphatidic acid acyltransferase active site for catalysis within the membrane bilayer

Phosphatidic acid (PA), the central intermediate in membrane phospholipid synthesis, is generated by two acyltransferases in a pathway conserved in all life forms. The second step in this pathway is catalyzed by 1-acyl-sn-glycerol-3-phosphate acyltransferase, called PlsC in bacteria. Here we present the crystal structure of PlsC from Thermotoga maritima, revealing an unusual hydrophobic/aromatic N-terminal two-helix motif linked to an acyltransferase αβ-domain that contains the catalytic HX₄D motif. PlsC dictates the acyl chain composition of the 2-position of phospholipids, and the acyl chain selectivity 'ruler' is an appropriately placed and closed hydrophobic tunnel. We confirmed this by site-directed mutagenesis and membrane composition analysis of Escherichia coli cells that expressed mutant PlsC. Molecular dynamics (MD) simulations showed that the two-helix motif represents a novel substructure that firmly anchors the protein to one leaflet of the membrane. This binding mode allows the PlsC active site to acylate lysophospholipids within the membrane bilayer by using soluble acyl donors.

The detoxifying effects of structural elements of persimmon tannin on Chinese cobra phospholipase A 2 correlated with their structural disturbing effects well

The effects of persimmon tannin (PT) characteristic structural elements on Naja atra phospholipase A₂ (PLA₂)-induced lethality, myotoxicity, and hemolysis in mice models were determined. In addition, methods including surface plasmon resonance, dynamic light scattering, and Fourier transform infrared spectroscopy were explored to uncover the possible detoxifying mechanisms of PT on snake venom PLA₂. Our results revealed that PT characteristic elements (EGCG, ECG, A-type EGCG dimer, and A-type ECG dimer) could neutralize the lethality, myotoxicity, and hemolysis of PLA₂. Moreover, the detoxifying effects of the four structural elements correlated with their structural disturbing effects well. Our results proved that A-type EGCG dimer and A-type ECG dimer may be structural requirements for the detoxifying effects of PT. We propose that the high affinity of A-type EGCG dimer and A-type ECG dimer for PLA₂ and the considerable spatial structural disturbance of PLA₂ induced by the dimers may be responsible for their antilethality, antimyotoxicity, and antihemolysis on Chinese cobra PLA₂in vivo.

Coralsnake Venomics: Analyses of Venom Gland Transcriptomes and Proteomes of Six Brazilian Taxa

Venom gland transcriptomes and proteomes of six Micrurus taxa (M. corallinus, M. lemniscatus carvalhoi, M. lemniscatus lemniscatus, M. paraensis, M. spixii spixii, and M. surinamensis) were investigated, providing the most comprehensive, quantitative data on Micrurus venom composition to date, and more than tripling the number of Micrurus venom protein sequences previously available. The six venomes differ dramatically. All are dominated by 2-6 toxin classes that account for 91-99% of the toxin transcripts. The M. s. spixii venome is compositionally the simplest. In it, three-finger toxins (3FTxs) and phospholipases A₂ (PLA₂s) comprise >99% of the toxin transcripts, which include only four additional toxin families at levels ≥0.1%. Micrurus l. lemniscatus venom is the most complex, with at least 17 toxin families. However, in each venome, multiple structural subclasses of 3FTXs and PLA₂s are present. These almost certainly differ in pharmacology as well. All venoms also contain phospholipase B and vascular endothelial growth factors. Minor components (0.1-2.0%) are found in all venoms except that of M. s. spixii. Other toxin families are present in all six venoms at trace levels (<0.005%). Minor and trace venom components differ in each venom. Numerous novel toxin chemistries include 3FTxs with previously unknown 8- and 10-cysteine arrangements, resulting in new 3D structures and target specificities. 9-cysteine toxins raise the possibility of covalent, homodimeric 3FTxs or heterodimeric toxins with unknown pharmacologies. Probable muscarinic sequences may be reptile-specific homologs that promote hypotension via vascular mAChRs. The first complete sequences are presented for 3FTxs putatively responsible for liberating glutamate from rat brain synaptosomes. Micrurus C-type lectin-like proteins may have 6-9 cysteine residues and may be monomers, or homo- or heterodimers of unknown pharmacology. Novel KSPIs, 3× longer than any seen previously, appear to have arisen in three species by gene duplication and fusion. Four species have transcripts homologous to the nociceptive toxin, (MitTx) α-subunit, but all six species had homologs to the β-subunit. The first non-neurotoxic, non-catalytic elapid phospholipase A₂s are reported. All are probably myonecrotic. Phylogenetic analysis indicates that the six taxa diverged 15-35 million years ago and that they split from their last common ancestor with Old World elapines nearly 55 million years ago. Given their early diversification, many cryptic micrurine taxa are anticipated.

[Ants: a chemical library of anticancer molecules]

Animal venoms are complex mixtures containing simple organic molecules, proteins, peptides, and other bioactive elements with extraordinary biological properties associated with their ability to act on a number of molecular receptors in the process of incapacitating their target organisms. In such a context, arthropod venoms are invaluable sources of bioactive substances, with therapeutic interest but the limited availability of some venom such as those from ants, has restricted the potential that these biomolecules could represent. We investigated for the first time transcriptomic expression from the ant species Tetramorium bicarinatum. Four hundred randomly selected clones from cDNA libraries were sequenced and a total of 374 expressed sequence tags (ESTs) were generated. Based on the results of BLAST searches, these sequences were clustered and assembled into 269 contigs. About 72% (269) of these matched BLASTx hits with an interesting diversity and unusual abundance of cellular transcripts (48%) related to gene and protein expression reflecting the specialization of this tissue. In addition, transcripts encoding transposases were relatively highly expressed (14%). It may be that transposable elements are present and that their presence accounts for some of the variation in venom toxins. About twenty per cent of the ESTs were categorized as putative toxins, the major part represented by allergens (48% of the total venom toxins) such as pilosulin 5, sol i 3 and Myp p I and II. Several contigs encoding enzymes, including zinc-metalloproteases (17%) that are likely involved in the processing and activation of venom proteins/peptides, were also identified from the library. In addition, a number of sequences (8%) had no significant similarity to any known sequence which indicates a potential source of for the discovery of new toxins. In order to provide a global insight on the transcripts expressed in the venom gland of the Brazilian ant species Tetramorium bicarinatum and to unveil the potential of their products, high-throughput expressed sequence tags were generated using Illumina paired-end sequencing technology. A total of 212 371 758 pairs of quality-filtered, 100-base-pair Illumina reads were obtained. The de novo assemblies yielded 36 042 contigs for which 27 873 have at least one predicted ORF among which 59.77% produce significant hits in the available databases. The investigation of the read mapping toxin class revealed and confirmed a high diversification with the major part consistent with the classical hymenopteran venom protein signature represented by venom allergen (33.3%) followed by a diverse toxin-expression profile including several distinct isoforms of phospholipase A1 and A2, venom serine protease, hyaluronidase, protease inhibitor and secapin. Moreover, our results revealed for the first time the presence of toxin-like peptides that have been previously identified from unrelated venomous animals such as waprin-like (snakes) and agatoxins (spiders and conus). These studies provide a first insight of the gene expression scenario of the venom gland of T. bicarinatum which might contribute to acquiring a more comprehensive view about the origin and functional diversity of venom proteins of this ant. Based on such results, we conducted cytotoxic tests from the crude venom of T.bicarinatum ant and reported toxic effect on tumoral cells lines from one of the fifth of the most frequently occurring cancers with a 3-year survival rate of only 30%. In such a context, new therapeutic strategies are essential and the discovery of new molecules in ant venom could be one possible avenue. Thus our project aims to characterize, from the crude venom of T.bicarinatum, the molecule(s) which have potential anti-cancerous toxicity as well as their mechanisms of action.

Heterogeneous behavior of metalloproteins toward metal ion binding and selectivity: insights from molecular dynamics studies

About one-third of the existing proteins require metal ions as cofactors for their catalytic activities and structural complexities. While many of them bind only to a specific metal, others bind to multiple (different) metal ions. However, the exact mechanism of their metal preference has not been deduced to clarity. In this study, we used molecular dynamics (MD) simulations to investigate whether a cognate metal (bound to the structure) can be replaced with other similar metal ions. We have chosen seven different proteins (phospholipase A2, sucrose phosphatase, pyrazinamidase, cysteine dioxygenase (CDO), plastocyanin, monoclonal anti-CD4 antibody Q425, and synaptotagmin 1 C2B domain) bound to seven different divalent metal ions (Ca(2+), Mg(2+), Zn(2+), Fe(2+), Cu(2+), Ba(2+), and Sr(2+), respectively). In total, 49 MD simulations each of 50 ns were performed and each trajectory was analyzed independently. Results demonstrate that in some cases, cognate metal ions can be exchanged with similar metal ions. On the contrary, some proteins show binding affinity specifically to their cognate metal ions. Surprisingly, two proteins CDO and plastocyanin which are known to bind Fe(2+) and Cu(2+), respectively, do not exhibit binding affinity to any metal ion. Furthermore, the study reveals that in some cases, the active site topology remains rigid even without cognate metals, whereas, some require them for their active site stability. Thus, it will be interesting to experimentally verify the accuracy of these observations obtained computationally. Moreover, the study can help in designing novel active sites for proteins to sequester metal ions particularly of toxic nature.

Cobra venom proteome and glycome determined from individual snakes of Naja atra reveal medically important dynamic range and systematic geographic variation

Recent progress in snake venomics has shed much light on the intra-species variation among the toxins from different geographical regions and has provided important information for better snakebite management. Most previous reports on snake venomics were based on venoms pooled from different snakes. In this study, we present the proteomic and glycomic profiles of venoms from individual Naja atra snakes. The results reveal wide dynamic range of three-finger toxins. Systematic classification based on cardiotoxin (CTX-) profiles of A2/A4 and A6, respectively, allowed the identification of two putative subspecies of Taiwan cobra from the eastern and western regions. We also identified four major N-glycan moieties on cobra snake venom metalloproteinase on the bi-antennary glycan core. ELISA showed that these glycoproteins (<3%) could elicit much higher antibody response in antiserum when compared to other high-abundance cobra venom toxins such as small molecular weight CTXs (~60%). By removing these high-molecular weight glycoproteins from the immunogen, we demonstrated better protection than that achieved with conventional crude venom immunization in mice challenged by crude venom. We conclude that both intra-species and inter-individual variations of proteomic and glycomic profiles of snake venomics should be considered to provide better antivenomic approach for snakebite management.

BIOLOGICAL SIGNIFICANCE

Based on the proteomic and glycomic profiles of venoms obtained from individual snakes, we demonstrated a surprisingly wide dynamic range and geographical variation of three-finger toxins in cobra venomics. This provides a reasonable explanation for the variable neutralization effects of antivenom treatment on victims suffering from cobra snakebite and suggests a simple and economic method to produce potent antivenom with better efficacy. Since two major venomic profiles with distinct dynamic ranges were observed for Taiwan cobra venoms isolated from the eastern and western regions, the current venomic profile should be used as a quality control for future production of antivenom in clinical applications.

De Novo sequencing and transcriptome analysis for Tetramorium bicarinatum: a comprehensive venom gland transcriptome analysis from an ant species

BACKGROUND

Arthropod venoms are invaluable sources of bioactive substances with biotechnological application. The limited availability of some venoms, such as those from ants, has restricted the knowledge about the composition and the potential that these biomolecules could represent. In order to provide a global insight on the transcripts expressed in the venom gland of the Brazilian ant species Tetramorium bicarinatum and to unveil the potential of its products, high-throughput approach using Illumina technology has been applied to analyze the genes expressed in active venom glands of this ant species.

RESULTS

A total of 212,371,758 pairs of quality-filtered, 100-base-pair Illumina reads were obtained. The de novo assemblies yielded 36,042 contigs for which 27,873 have at least one predicted ORF among which 59.77% produce significant hits in the available databases. The investigation of the reads mapping toxin class revealed a high diversification with the major part consistent with the classical hymenopteran venom protein signature represented by venom allergen (33.3%), followed by a diverse toxin-expression profile including several distinct isoforms of phospholipase A1 and A2, venom serine protease, hyaluronidase, protease inhibitor and secapin. Moreover, our results revealed for the first time the presence of toxin-like peptides that have been previously identified from unrelated venomous animals such as waprin-like (snakes) and agatoxins (spiders and conus).The non-toxin transcripts were mainly represented by contigs involved in protein folding and translation, consistent with the protein-secretory function of the venom gland tissue. Finally, about 40% of the generated contigs have no hits in the databases with 25% of the predicted peptides bearing signal peptide emphasizing the potential of the investigation of these sequences as source of new molecules. Among these contigs, six putative novel peptides that show homologies with previously identified antimicrobial peptides were identified.

CONCLUSIONS

To the best of our knowledge, this work reports the first large-scale analysis of genes transcribed by the venomous gland of the ant species T. bicarinatum and helps with the identification of Hymenoptera toxin arsenal. In addition, results from this study demonstrate that de novo transcriptome assembly allows useful venom gene expression analysis in a species lacking a genome sequence database.

A facile assay to monitor secretory phospholipase A₂ using 8-anilino-1-naphthalenesulfonic acid

Secretory phospholipases A2 (sPLA2s) are present in snake venoms, serum, and biological fluids of patients with various inflammatory, autoimmune and allergic disorders. Lipid mediators in the inflammatory processes have potential value for controlling phospholipid metabolism through sPLA2 inhibition. Thus, it demands the need for screening of potential leads for sPLA2 inhibition. To date, sPLA2 activity has been assayed using expensive radioactive or chromogenic substrates, thereby limiting a large number of assays. In this study, a simple and sensitive NanoDrop assay was developed using non-fluorogenic and non-chromogenic phospholipid substrate 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and 8-anilino-1-naphthalenesulfonic acid (ANS) as interfacial hydrophobic probe. The modified assay required a 10ng concentration of sPLA2. ANS, as a strong anion, binds predominantly to cationic group of choline head of DMPC through ion pair formation, imparting hydrophobicity and lipophilicity and resulting in an increase in fluorescence. Triton X-100 imparts correct geometrical space during sPLA2 catalyzing DMPC, releasing lysophospholipid and acidic myristoyl acid, which in turn alters the hydrophobic environment prevailing around ANS-DMPC, which leads to weakening of the electrostatic ion pair interaction between DMPC and ANS ensuing decrease in fluorescence. These characteristic fluorescence changes between DMPC and ANS in response to sPLA2 catalysis are well documented and validated in this study.

Interaction of characteristic structural elements of persimmon tannin with Chinese cobra PLA2

To more fully understand the mechanism by which persimmon tannin (PT) inhibited phospholipase A2 (PLA2) and the structural requirements of PT for the inhibition, the interactions between PLA2 and seven characteristic structural elements of PT including epigallocatechin-3-gallate (EGCG), myricetin, epicatechin-3-gallate (ECG), epicatechin-3-gallate-(4β → 8, 2β → O → 7)-epicatechin-3-gallate (A-type ECG dimer), epigallocatechin-3-gallate-(4β → 8, 2β → O → 7)-epigallocatechin-3-gallate (A-type EGCG dimer), epicatechin-(4β → 8, 2β → O → 7)-epicatechin (A-type EC dimer) and epicatechin-(4β → 8)-epicatechin (B-type EC dimer) were studied by enzymatic and spectroscopic methods. Molecular docking was also used to explore the possible residues involved in the interactions. The results revealed that A-type EGCG dimer and A-type ECG dimer showed higher inhibitory effects on the catalytic activity of PLA2 than monomers and B-type dimer. They induced greater conformational changes in PLA2 than other structural elements. In addition, molecular docking studies revealed that expect for lysine residues, other residues such as Trp18, Try27, Gly29, His47 and Tyr63 were involved in the interactions. We propose that A-type EGCG and ECG dimer units may be structural requirements for the interaction between PT and PLA2. Our data provide an additional structural basis for anti-PLA2 activity of persimmon tannin.

Structure and dynamics of adeno-associated virus serotype 1 VP1-unique N-terminal domain and its role in capsid trafficking

The importance of the phospholipase A2 domain located within the unique N terminus of the capsid viral protein VP1 (VP1u) in parvovirus infection has been reported. This study used computational methods to characterize the VP1 sequence for adeno-associated virus (AAV) serotypes 1 to 12 and circular dichroism and electron microscopy to monitor conformational changes in the AAV1 capsid induced by temperature and the pHs encountered during trafficking through the endocytic pathway. Circular dichroism was also used to monitor conformational changes in AAV6 capsids assembled from VP2 and VP3 or VP1, VP2, and VP3 at pH 7.5. VP1u was predicted (computationally) and confirmed (in solution) to be structurally ordered. This VP domain was observed to undergo a reversible pH-induced unfolding/refolding process, a loss/gain of α-helical structure, which did not disrupt the capsid integrity and is likely facilitated by its difference in isoelectric point compared to the other VP sequences assembling the capsid. This study is the first to physically document conformational changes in the VP1u region that likely facilitate its externalization from the capsid interior during infection and establishes the order of events in the escape of the AAV capsid from the endosome en route to the nucleus.

Using hydrogen/deuterium exchange mass spectrometry to define the specific interactions of the phospholipase A2 superfamily with lipid substrates, inhibitors, and membranes

The phospholipase A(2) (PLA(2)) superfamily consists of 16 groups and many subgroups and constitutes a diverse set of enzymes that have a common catalytic activity due to convergent evolution. However, different PLA(2) types have unique three-dimensional structures and catalytic residues as well as specific tissue localization and distinct biological functions. Understanding how the different PLA(2) enzymes associate with phospholipid membranes, specific phospholipid substrate molecules, and inhibitors on a molecular basis has advanced in recent years due to the introduction of hydrogen/deuterium exchange mass spectrometry. Its theory, practical considerations, and application to understanding PLA(2)/membrane interactions are addressed.

Assay of phospholipases A(2) and their inhibitors by kinetic analysis in the scooting mode

Several cellular processes are regulated by interfacial catalysis on biomembrane surfaces. Phospholipases A₂ (PLA₂) are interesting not only as prototypes for interfacial catalysis, but also because they mobilize precursors for the biosynthesis of eicosanoids and platelet activating factor, and these agents ultimately control a wide range of secretory and inflammatory processes. Since PLA₂ carry out their catalytic function at membrane surfaces, the kinetics of these enzymes depends on what the enzyme ‘sees’ at the interface, and thus the observed rate is profoundly influenced by the organization and dynamics of the lipidwater interface (‘quality of the interface’). In this review we elaborate the advantages of monitoring interfacial catalysis in the scooting mode, that is, under the conditions where the enzyme remains bound to vesicles for several thousand catalytic turnover cycles. Such a highly processive catalytic turnover in the scooting mode is useful for a rigorous and quantitative characterization of the kinetics of interfacial catalysis. This analysis is now extended to provide insights into designing strategy for PLA₂ assays and screens for their inhibitors.

Therapeutic application of natural inhibitors against snake venom phospholipase A(2)

Natural inhibitors occupy an important place in the potential to neutralize the toxic effects caused by snake venom proteins and enzymes. It has been well recognized for several years that animal sera, some of the plant and marine extracts are the most potent in neutralizing snake venom phospholipase A(2) (svPLA(2)). The implication of this review to update the latest research work which has been accomplished with svPLA(2) inhibitors from various natural sources like animal, marine organisms presents a compilation of research in this field over the past decade and revisiting the previous research report including those found in plants. In addition to that the bioactive compounds/inhibitor molecules from diverse sources like aristolochic alkaloid, flavonoids and neoflavonoids from plants, hydrocarbones -2, 4 dimethyl hexane, 2 methylnonane, and 2, 6 dimethyl heptane obtained from traditional medicinal plants Tragia involucrata (Euphorbiaceae) member of natural products involved for the inhibitory potential of phospholipase A(2) (PLA(2)) enzymes in vitro and also decrease both oedema induced by snake venom as well as human synovial fluid PLA(2). Besides marine natural products that inhibit PLA(2) are manoalide and its derivatives such as scalaradial and related compounds, pseudopterosins and vidalols, tetracylne from synthetic chemicals etc. There is an overview of the role of PLA(2) in inflammation that provides a rationale for seeking inhibitors of PLA(2) as anti-inflammatory agents. However, more studies should be considered to evaluate antivenom efficiency of sera and other agents against a variety of snake venoms found in various parts of the world. The implications of these new groups of svPLA(2) toxin inhibitors in the context of our current understanding of snake biology as well as in the development of new novel antivenoms therapeutics agents in the efficient treatment of snake envenomations are discussed.

Structural basis for functional diversity of animal toxins

The diversity of biological functions that are exerted by toxins from snake and scorpion venoms is associated with a limited number of structural frameworks. At present, one predominant basic fold has been observed among scorpion toxins whereas six folds have been found among snake toxins. Most toxin folds have the capacity to accept multiple insertions, deletions and mutations and to exert various recognition functions. We suggest that such folds may serve as guides to engineer new protein functions.

Surface pressure-dependent interactions of secretory phospholipase A2 with zwitterionic phospholipid membranes

The hydrolytic activity of secretory phospholipase A(2) (PLA(2)) is regulated by many factors, including the physical state of substrate aggregates and the chemical nature of phospholipid molecules. In order to achieve strong binding of PLA(2) on its substrates, many previous works have used anionic lipid dispersion to characterize the orientation and penetration depth of PLA(2) molecules on membrane surfaces. In this study, we applied monolayer technique with controllable surface area to investigate the PLA(2)s of Taiwan cobra venom and bee venom on zwitterionic phophatidylcholine monolayers and demonstrated an optimum hydrolytic activity at a surface pressure of 18 and 24 mN/m, respectively. By combining polarized attenuated total reflection Fourier-transform infrared spectroscopy and monolayer-binding experiments, we found that the amount of membrane-bound PLA(2) decreased markedly as the surface pressure of the monolayer was increased. Interestingly, the insertion area of the PLA(2)s decreased to near zero as the surface pressure increased to the optimum pressure for hydrolytic activity. On the basis of the measured infrared dichroic ratio, the orientation of the PLA(2)s bound to zwitterionic membranes was similar to that observed on a negatively charged membrane and was independent of the surface pressure. Our findings suggest that both PLA(2)s were located on the membrane surface rather than penetrating the membrane bilayer and that the deeply inserted mode is not a favorable condition for the hydrolysis of phospholipids in zwitterionic phospholipid membranes. The results are discussed in terms of the easy access of catalytic water for the PLA(2) activity and the mobilization of its substrate and product to facilitate the catalytic process.

Secretory phospholipase A2 activity toward diverse substrates

We have studied secretory phospholipase A(2)-IIA (sPLA(2)) activity toward different phospholipid analogues by performing biophysical characterizations and molecular dynamics simulations. The phospholipids were natural substrates, triple alkyl phospholipids, a prodrug anticancer etherlipid, and an inverted ester. The latter were included to study head group-enzyme interactions. Our simulation results show that the lipids are optimally placed into the binding cleft and that water molecules can freely reach the active site through a well-defined pathway; both are indicative that these substrates are efficiently hydrolyzed, which is in good agreement with our experimental data. The phospholipid analogue with three alkyl side chains forms aggregates of different shapes with no well-defined sizes due to its cone-shape structure. Phosphatidylglycerol and phosphatidylcholine head groups interact with specific charged residues, but relatively large fluctuations are observed, suggesting that these interactions are not necessarily important for stabilizing substrate binding to the enzyme.

Light controls phospholipase A2alpha and beta gene expression in Citrus sinensis

The low-molecular weight secretory phospholipase A2alpha (CssPLA2alpha) and beta (CsPLA2beta) cloned in this study exhibited diurnal rhythmicity in leaf tissue of Citrus sinensis. Only CssPLA2alpha displayed distinct diurnal patterns in fruit tissues. CssPLA2alpha and CsPLA2beta diurnal expression exhibited periods of approximately 24 h; CssPLA2alpha amplitude averaged 990-fold in the leaf blades from field-grown trees, whereas CsPLA2beta amplitude averaged 6.4-fold. Diurnal oscillation of CssPLA2alpha and CsPLA2beta gene expression in the growth chamber experiments was markedly dampened 24 h after transfer to continuous light or dark conditions. CssPLA2alpha and CsPLA2beta expressions were redundantly mediated by blue, green, red and red/far-red light, but blue light was a major factor affecting CssPLA2alpha and CsPLA2beta expression. Total and low molecular weight CsPLA2 enzyme activity closely followed diurnal changes in CssPLA2alpha transcript expression in leaf blades of seedlings treated with low intensity blue light (24 micromol m(-2) s(-1)). Compared with CssPLA2alpha basal expression, CsPLA2beta expression was at least 10-fold higher. Diurnal fluctuation and light regulation of PLA2 gene expression and enzyme activity in citrus leaf and fruit tissues suggests that accompanying diurnal changes in lipophilic second messengers participate in the regulation of physiological processes associated with phospholipase A2 action.

Differential gene expression profiles in the venom gland/sac of Eumenes pomiformis (Hymenoptera: Eumenidae)

To search for novel transcripts encoding biologically active venom components, a subtractive cDNA library specific to the venom gland and sac (gland/sac) of a solitary hunting wasp species, Eumenes pomiformis Fabricius (1781), was constructed by suppression subtractive hybridization. A total of 541 expressed sequence tags (ESTs) were clustered and assembled into 102 contigs (31 multiple sequences and 71 singletons). In total, 37 cDNAs were found in the library via BLASTx searching and manual annotation. Eight contigs (337 ESTs) encoding short venom peptides (10 to 16 amino acids) occupied 62% of the library. The deduced amino acid sequence (78 amino acids) of a novel venom peptide transcript shared sequence similarity with trypsin inhibitors and dendrotoxin-like venom peptides known to be K(+) channel blockers, implying that this novel peptide may play a role in the paralysis of prey. In addition to phospholipase A2 and hyaluronidase, which are known to be the main components of wasp venoms, several transcripts encoding enzymes, including three metallopeptidases and a decarboxylase likely involved in the processing and activation of venomous proteins, peptides, amines, and neurotransmitters, were also isolated from the library. The presence of a transcript encoding a putative insulin/insulin-like peptide binding protein suggests that solitary hunting wasps use their venom to control their prey, leading to larval growth cessation. The abundance of these venom components in the venom gland/sac and in the alimentary canal was confirmed by quantitative real-time PCR. Discovery of venom gland/sac-specific transcripts should promote further studies on biologically active components in the venom of solitary hunting wasps.

Substrate efflux propensity plays a key role in the specificity of secretory A-type phospholipases

To better understand the principles underlying the substrate specificity of A-type phospholipases (PLAs), a high throughput mass spectrometric assay was employed to study the effect of acyl chain length and unsaturation of phospholipids on their rate of hydrolysis by three different secretory PLAs in micelles and vesicle bilayers. With micelles, each enzyme responded differently to substrate acyl chain unsaturation and double bond position, probably reflecting differences in the accommodative properties of their substrate binding sites. Experiments with saturated acyl positional isomers indicated that the length of the sn2 chain was more critical than that of the sn1 chain, suggesting tighter association of the former with the enzyme. Only the first 9-10 carbons of the sn2 acyl chain seem to interact intimately with the active site. Strikingly, no discrimination between positional isomers was observed with vesicles, and the rate of hydrolysis decreased far more with increasing chain length than with micelles, suggesting that translocation of the phospholipid substrate to the active site is rate-limiting with bilayers. Supporting this conclusion, acyl chain structure affected hydrolysis and spontaneous intervesicle transfer, which correlates with lipid efflux propensity, analogously. We conclude that substrate efflux propensity plays a more important role in the specificity of secretory PLA(2)s than commonly thought and could also be a key attribute in phospholipid homeostasis in which (unknown) PLA(2)s are key players.

Relationship between antibody 2F5 neutralization of HIV-1 and hydrophobicity of its heavy chain third complementarity-determining region

The membrane-proximal external region (MPER) of the HIV-1 gp41 transmembrane glycoprotein is the target of the broadly neutralizing antibody 2F5. Prior studies have suggested a two-component mechanism for 2F5-mediated neutralization involving both structure-specific recognition of a gp41 protein epitope and nonspecific interaction with the viral lipid membrane. Here, we mutationally alter a hydrophobic patch on the third complementarity-determining region of the heavy chain (CDR H3) of the 2F5 antibody and assess the abilities of altered 2F5 variants to bind gp41 and to neutralize diverse strains of HIV-1. CDR H3 alterations had little effect on the affinity of 2F5 variants for a peptide corresponding to its gp41 epitope. In contrast, strong effects and a high degree of correlation (P < 0.0001) were found between virus neutralization and CDR H3 hydrophobicity, as defined by predicted free energies of transfer from water to a lipid bilayer interface or to octanol. The effect of CDR H3 hydrophobicity on neutralization was independent of isolate sensitivity to 2F5, and CDR H3 variants with tryptophan substitutions were able to neutralize HIV-1 approximately 10-fold more potently than unmodified 2F5. A threshold was observed for increased hydrophobicity of the 2F5 CDR H3 loop beyond which effects on 2F5-mediated neutralization leveled off. Together, the results provide a more complete understanding of the 2F5 mechanism of HIV-1 neutralization and indicate ways to enhance the potency of MPER-directed antibodies.

Differential gene expression profiles in the venom gland/sac of Orancistrocerus drewseni (Hymenoptera: Eumenidae)

To determine differential gene expression profiles in the venom gland and sac (gland/sac) of a solitary hunting wasp species, Orancistrocerus drewseni Saussure (1857), a subtractive cDNA library was constructed by suppression subtractive hybridization. A total of 498 expressed sequence tags (EST) were clustered and assembled into 205 contigs (94 multiple sequences and 111 singletons). About 65% (134) of the contigs had matched BLASTx hits (E< or =10(-4)). Among these, 115 contigs had similarity to proteins with assigned molecular function in the Gene Ontology database, and most of them (112 contigs, 83%) were homologous to genes from Hymenoptera, particularly to Apis mellifera (98 contigs). The contigs encoding hyaluronidase and phospholipase A2, known to be main components of wasp venoms, were found in high frequencies (27 and 4%, respectively, as judged by the number of ESTs) in the gene ontology category of catalytic activity. Full-length open reading frames of hyaluronidase and phospholipase A2 were characterized and their abundance in the venom gland/sac was confirmed by quantitative real-time PCR. Several contigs encoding enzymes, including zinc-metallopeptidases that are likely involved in the processing and activation of venomous proteins or peptides, were also identified from the library. Discovery of venom gland/sac-specific genes should promote further studies on biologically active components in the venom of O. drewseni.

Crystal structure of a class XIB phospholipase A2 (PLA2): rice (oryza sativa) isoform-2 pla2 and an octanoate complex

Phospholipase A(2) catalyzes the specific hydrolysis of the sn-2 acyl bond of various glycerophospholipids, producing fatty acids and lysophospholipids. Phospholipase A(2)s (PLA(2)s) constitute a large superfamily of enzymes whose products are important for a multitude of signal transduction processes, lipid mediator release, lipid metabolism, development, plant stress responses, and host defense. The crystal structure of rice (Oryza sativa) isoform 2 phospholipase A(2) has been determined to 2.0 A resolution using sulfur SAD phasing, and shows that the class XIb phospholipases have a unique structure compared with other secreted PLA(2)s. The N-terminal half of the chain contains mainly loop structure, including the conserved Ca(2+)-binding loop, but starts with a short 3(10)-helix and also includes two short anti-parallel beta-strands. The C-terminal half is folded into three anti-parallel alpha-helices, of which the two first are also present in other secreted PLA(2)s and contain the conserved catalytic histidine and calcium liganding aspartate residues. The structure is stabilized by six disulfide bonds. The water structure around the calcium ion binding site suggests the involvement of a second water molecule in the mechanism for hydrolysis, the water-assisted calcium-coordinate oxyanion mechanism. The octanoate molecule in the complex structure is bound in a hydrophobic pocket, which extends to the likely membrane interface and is proposed to model the binding of the product fatty acid. Due to the differences in structure, the suggested surface for binding to the membrane has a different morphology in the rice PLA(2) compared with other phospholipases.