Crystal structure of cardiotoxin V from Taiwan cobra venom: pH-dependent conformational change and a novel membrane-binding motif identified in the three-finger loops of P-type cardiotoxin

Variability in the Spatial Structure of the Central Loop in Cobra Cytotoxins Revealed by X-ray Analysis and Molecular Modeling

Cobra cytotoxins (CTs) belong to the three-fingered protein family and possess membrane activity. Here, we studied cytotoxin 13 from Naja naja cobra venom (CT13Nn). For the first time, a spatial model of CT13Nn with both “water” and “membrane” conformations of the central loop (loop-2) were determined by X-ray crystallography. The “water” conformation of the loop was frequently observed. It was similar to the structure of loop-2 of numerous CTs, determined by either NMR spectroscopy in aqueous solution, or the X-ray method. The “membrane” conformation is rare one and, to date has only been observed by NMR for a single cytotoxin 1 from N. oxiana (CT1No) in detergent micelle. Both CT13Nn and CT1No are S-type CTs. Membrane-binding of these CTs probably involves an additional step—the conformational transformation of the loop-2. To confirm this suggestion, we conducted molecular dynamics simulations of both CT1No and CT13Nn in the Highly Mimetic Membrane Model of palmitoiloleoylphosphatidylglycerol, starting with their “water” NMR models. We found that the both toxins transform their “water” conformation of loop-2 into the “membrane” one during the insertion process. This supports the hypothesis that the S-type CTs, unlike their P-type counterparts, require conformational adaptation of loop-2 during interaction with lipid membranes.

Three-Finger Toxins from Brazilian Coral Snakes: From Molecular Framework to Insights in Biological Function

Studies on 3FTxs around the world are showing the amazing diversity in these proteins both in structure and function. In Brazil, we have not realized the broad variety of their amino acid sequences and probable diversified structures and targets. In this context, this work aims to conduct an in silico systematic study on available 3FTxs found in Micrurus species from Brazil. We elaborated a specific guideline for this toxin family. First, we grouped them according to their structural homologue predicted by HHPred server and further curated manually. For each group, we selected one sequence and constructed a representative structural model. By looking at conserved features and comparing with the information available in the literature for this toxin family, we managed to point to potential biological functions. In parallel, the phylogenetic relationship was estimated for our database by maximum likelihood analyses and a phylogenetic tree was constructed including the homologous 3FTx previously characterized. Our results highlighted an astonishing diversity inside this family of toxins, showing some groups with expected functional similarities to known 3FTxs, and pointing out others with potential novel roles and perhaps structures. Moreover, this classification guideline may be useful to aid future studies on these abundant toxins.

Identification of a α-helical molten globule intermediate and structural characterization of β-cardiotoxin, an all β-sheet protein isolated from the venom of Ophiophagus hannah (king cobra)

β-Cardiotoxin is a novel member of the snake venom three-finger toxin (3FTX) family. This is the first exogenous protein to antagonize β-adrenergic receptors and thereby causing reduction in heart rates (bradycardia) when administered into animals, unlike the conventional cardiotoxins as reported earlier. 3FTXs are stable all β-sheet peptides with 60-80 amino acid residues. Here, we describe the three-dimensional crystal structure of β-cardiotoxin together with the identification of a molten globule intermediate in the unfolding pathway of this protein. In spite of the overall structural similarity of this protein with conventional cardiotoxins, there are notable differences observed at the loop region and in the charge distribution on the surface, which are known to be critical for cytolytic activity of cardiotoxins. The molten globule intermediate state present in the thermal unfolding pathway of β-cardiotoxin was however not observed during the chemical denaturation of the protein. Interestingly, circular dichroism (CD) and NMR studies revealed the presence of α-helical secondary structure in the molten globule intermediate. These results point to substantial conformational plasticity of β-cardiotoxin, which might aid the protein in responding to the sometimes conflicting demands of structure, stability, and function during its biological lifetime.

The role of hydrophobic /hydrophilic balance in the activity of structurally flexible vs. rigid cytolytic polypeptides and analogs developed on their basis

INTRODUCTION

Being important representatives of various proteomes, membrane-active cationic peptides (CPs) are attractive objects as lead compounds in the design of new antibacterial, anticancer, antifungal, and antiviral molecules. Numerous CPs are found in insect and snake venoms, where many of them reveal cytolytic properties. Due to advances in omics technologies, the number of such peptides is growing dramatically. Areas covered: To understand structure-function relationships for CPs in a living cell, detailed analysis of their hydrophobic/hydrophilic properties is indispensable. We consider two structural classes of membrane-active CPs: latarcins (Ltc) from spider and cardiotoxins (CTXs) from snake venoms. While the former are void off disulfide bonds and conformationally flexible, the latter are structurally rigid and cross-linked with disulfide bonds. In order to elucidate structure-activity relationships behind their antibacterial, anticancer, and hemolytic effects, the properties of these polypeptides are considered on a side-by-side basis. Expert commentary: An ever-increasing number of venom-derived membrane-active polypeptides require new methods for identification of their functional propensities and sequence-based design of novel pharmacological substances. We address these issues considering a number of the designed peptides, based either on Ltc or CTX sequences. Experimental and computer modeling techniques required for these purposes are delineated.

Extensive Assessment of Various Computational Methods for Aspartate's pKa Shift

A series of computational methods for pK_a shift prediction are extensively tested on a set of benchmark protein systems, aiming at identifying pitfalls and evaluating their performance on high variants. Including 19 ASP residues in 10 protein systems, the benchmark set consists of both residues with highly shifted pK_a values as well as those varying little from the reference value, with an experimental RMS free energy differences of 2.49 kcal/mol with respect to blocked amino acid, namely the RMS pK_a shift being 1.82 pK_a units. The constant pH molecular dynamics (MD), alchemical methods, PROPKA3.1, and multiconformation continuum electrostatics give RMSDs of 1.52, 2.58, 1.37, and 3.52 pK_a units, respectively, on the benchmark set. The empirical scoring method is the most accurate one with extremely low computational cost, and the pH-dependent model is also able to provide accurate results, while the accuracy of MD sampling incorporating alchemical free energy simulation is prohibited by convergence achievement and the performance of conformational search incorporating multiconformation continuum electrostatics is bad. Former research works did not define statistical uncertainty with care and yielded the questionable conclusion that alchemical methods perform well in most benchmarks. In this work the traditional alchemical methods are thoroughly tested for high variants. We also performed the first application of nonequilibrium alchemical methods to the pK_a cases.

Ringhalexin from Hemachatus haemachatus: A novel inhibitor of extrinsic tenase complex

Anticoagulant therapy is used for the prevention and treatment of thromboembolic disorders. Blood coagulation is initiated by the interaction of factor VIIa (FVIIa) with membrane-bound tissue factor (TF) to form the extrinsic tenase complex which activates FX to FXa. Thus, it is an important target for the development of novel anticoagulants. Here, we report the isolation and characterization of a novel anticoagulant ringhalexin from the venom of Hemachatus haemachatus (African Ringhals Cobra). Amino acid sequence of the protein indicates that it belongs to the three-finger toxin family and exhibits 94% identity to an uncharacterized Neurotoxin-like protein NTL2 from Naja atra. Ringhalexin inhibited FX activation by extrinsic tenase complex with an IC50 of 123.8 ± 9.54 nM. It is a mixed-type inhibitor with the kinetic constants, Ki and Ki' of 84.25 ± 3.53 nM and 152.5 ± 11.32 nM, respectively. Ringhalexin also exhibits a weak, irreversible neurotoxicity on chick biventer cervicis muscle preparations. Subsequently, the three-dimensional structure of ringhalexin was determined at 2.95 Å resolution. This study for the first time reports the structure of an anticoagulant three-finger toxin. Thus, ringhalexin is a potent inhibitor of the FX activation by extrinsic tenase complex and a weak, irreversible neurotoxin.

DNA Aptamers against Taiwan Banded Krait α-Bungarotoxin Recognize Taiwan Cobra Cardiotoxins

Bungarus multicinctus α-bungarotoxin (α-Bgt) and Naja atra cardiotoxins (CTXs) share a common structural scaffold, and their tertiary structures adopt three-fingered loop motifs. Four DNA aptamers against α-Bgt have been reported previously. Given that the binding of aptamers with targeted proteins depends on structural complementarity, in this study, we investigated whether DNA aptamers against α-Bgt could also recognize CTXs. It was found that N.atra cardiotoxin 3 (CTX3) reduced the electrophoretic mobility of aptamers against α-Bgt. Analysis of the changes in the fluorescence intensity of carboxyfluorescein-labeled aptamers upon binding toxin molecules revealed that CTX3 and α-Bgt could bind the tested aptamers. Moreover, the aptamers inhibited the membrane-damaging activity and cytotoxicity of CTX3. In addition to CTX3, other N. atra CTX isotoxins also bound to the aptamer against α-Bgt. Taken together, our data indicate that aptamers against α-Bgt show cross-reactivity with CTXs. The findings that aptamers against α-Bgt also suppress the biological activities of CTX3 highlight the potential utility of aptamers in regard to the broad inhibition of snake venom three-fingered proteins.

Fusogenicity of Naja naja atra cardiotoxin-like basic protein on sphingomyelin vesicles containing oxidized phosphatidylcholine and cholesterol

This study investigated the effect of oxidized phosphatidylcholine (oxPC) and cholesterol (Chol) on Naja naja atra cardiotoxin-like basic protein (CLBP)-induced fusion and leakage in sphingomyelin (SM) vesicles. Compared with those on PC/SM/Chol vesicles, CLBP showed a lower activity to induce membrane permeability but a higher fusogenicity on oxPC/SM/Chol vesicles. A reduction in inner-leaflet fusion elucidated that CLBP fusogenicity was not in parallel to its membrane-leakage activity on oxPC/SM/Chol vesicles. The lipid domain formed by Chol and SM supported CLBP fusogenicity on oxPC/SM/Chol vesicles, while oxPC altered the interacted mode of CLBP with oxPC/SM/Chol vesicles as evidenced by Fourier transform infrared spectra analyses and colorimetric phospholipid/polydiacetylene membrane assay. Although CLBP showed similar binding affinity with PC/SM/Chol and oxPC/SM/Chol vesicles, the binding capability of CLBP with PC/SM/Chol and oxPC/SM/Chol vesicles was affected differently by NaCl. This emphasized that CLBP adopted different membrane interaction modes upon binding with PC/SM/Chol and oxPC/SM/Chol vesicles. CLBP induced fusion in vesicles containing oxPC bearing the aldehyde group, and aldehyde scavenger methoxyamine abrogated the CLBP ability to induce oxPC/SM/Chol fusion. Taken together, our data indicate that Chol and oxPC bearing aldehyde group alter the CLBP membrane-binding mode, leading to fusogenicity promotion while reducing the membrane-damaging activity of CLBP.

Identification and structural characterization of a new three-finger toxin hemachatoxin from Hemachatus haemachatus venom

Snake venoms are rich sources of biologically active proteins and polypeptides. Three-finger toxins are non-enzymatic proteins present in elapid (cobras, kraits, mambas and sea snakes) and colubrid venoms. These proteins contain four conserved disulfide bonds in the core to maintain the three-finger folds. Although all three-finger toxins have similar fold, their biological activities are different. A new three-finger toxin (hemachatoxin) was isolated from Hemachatus haemachatus (Ringhals cobra) venom. Its amino acid sequence was elucidated, and crystal structure was determined at 2.43 Å resolution. The overall fold is similar to other three-finger toxins. The structure and sequence analysis revealed that the fold is maintained by four highly conserved disulfide bonds. It exhibited highest similarity to particularly P-type cardiotoxins that are known to associate and perturb the membrane surface with their lipid binding sites. Also, the increased B value of hemachotoxin loop II suggests that loop II is flexible and may remain flexible until its interaction with membrane phospholipids. Based on the analysis, we predict hemachatoxin to be cardiotoxic/cytotoxic and our future experiments will be directed to characterize the activity of hemachatoxin.

Snake cytotoxins bind to membranes via interactions with phosphatidylserine head groups of lipids

The major representatives of Elapidae snake venom, cytotoxins (CTs), share similar three-fingered fold and exert diverse range of biological activities against various cell types. CT-induced cell death starts from the membrane recognition process, whose molecular details remain unclear. It is known, however, that the presence of anionic lipids in cell membranes is one of the important factors determining CT-membrane binding. In this work, we therefore investigated specific interactions between one of the most abundant of such lipids, phosphatidylserine (PS), and CT 4 of Naja kaouthia using a combined, experimental and modeling, approach. It was shown that incorporation of PS into zwitterionic liposomes greatly increased the membrane-damaging activity of CT 4 measured by the release of the liposome-entrapped calcein fluorescent dye. The CT-induced leakage rate depends on the PS concentration with a maximum at approximately 20% PS. Interestingly, the effects observed for PS were much more pronounced than those measured for another anionic lipid, sulfatide. To delineate the potential PS binding sites on CT 4 and estimate their relative affinities, a series of computer simulations was performed for the systems containing the head group of PS and different spatial models of CT 4 in aqueous solution and in an implicit membrane. This was done using an original hybrid computational protocol implementing docking, Monte Carlo and molecular dynamics simulations. As a result, at least three putative PS-binding sites with different affinities to PS molecule were delineated. Being located in different parts of the CT molecule, these anion-binding sites can potentially facilitate and modulate the multi-step process of the toxin insertion into lipid bilayers. This feature together with the diverse binding affinities of the sites to a wide variety of anionic targets on the membrane surface appears to be functionally meaningful and may adjust CT action against different types of cells.

Biomimetic affinity purification of cardiotoxin and its pharmacological effects on the nervous system

Cobra venom is a very precious natural resource. The traditional method for purification of cardiotoxin from cobra venom is a multi-step, high cost, and low recovery procedure. By molecular modeling and docking with SYBYL software, we designed and synthesized an affinity ligand, m-aminobenzoic acid, for high efficiency purification of this therapeutically useful Chinese cobra venom cardiotoxin. The one-step recovery of cardiotoxin reached 64% and the purity reached 92% upon purification. The binding capacity of this synthetic ligand was 9.1 mg cardiotoxin/g moist weight gel and the affinity constant for cardiotoxin was 5.5 x 10(3) M(-1). Unlike a natural affinity ligand, this synthetic ligand is highly stable, and has great potential for industrial scale production of cardiotoxin. In addition, we examined the effects of cardiotoxin on the nervous system in a mouse model. Results showed that cardiotoxin could maintain analgesic effects for 120 min with a dose of less than 0.06 mg/kg (2.8% of the LD(50)). Administration of 0.12 mg/kg cardiotoxin could improve scopolamine impairments of memory in mice. These results suggest that cardiotoxin may be a potential drug for nervous system diseases.

Purification and characterization of a novel antinociceptive toxin from Cobra venom (Naja naja atra)

Snake venoms have demonstrated antinociceptive activity, and certain isolated neurotoxins have demonstrated significant analgesia in animal models. Here we report a novel analgesic toxin which was isolated from Naja naja atra and was given the name 'najanalgesin'. The LD(50) of the crude venom and najanalgesin were 0.89mg/kg and 2.69mg/kg, respectively. We used the writhing test and hot plate test to evaluate the antinociceptive properties of the crude venom and najanalgesin after intraperitoneal (ip) administration. The analgesic mechanism of najanalgesin was also studied. The response latency time was significantly prolonged in the hot plate test after ip administration of the crude venom of Naja naja atra (0.111-0.445mg/kg) in a dose-dependent manner. Najanalgesin (1mg/kg) elicited almost the same antinociceptive effect as that of the crude venom of Naja naja atra at the dose of 0.445mg/kg and remained for 6h after intraperitoneal injection, shown by hot plate test. The percentage of increase in the latency time for the venom and the najanalgesin 3h after drug administration was 96.2% and 112%, respectively. The number of writhes decreased to almost 1/3, 1/6, and 1/12 of the NS (physiological saline) group after intraperitoneal administration of najanalgesin at 0.25, 0.5, and 1.0mg/kg, respectively. Pretreatment with atropine (1mg/kg) or naloxone (3mg/kg) blocked the antinociception of najanalgesin in the hot plate test. Based on the sequence information, najanalgesin is found to be highly homologous with the conventional CTXs (cardiotoxins). To our knowledge, no study had previously reported that a toxin which was homologous with CTXs possessed the antinociceptive activity. Thus, this is the first report that the antinociceptive effect induced by najanalgesin is mediated by cholinergic and opioidergic mechanisms.

The mechanism of cytotoxicity by Naja naja atra cardiotoxin 3 is physically distant from its membrane-damaging effect

In order to dissect out whether multiple activities of cardiotoxins (CTXs) are connected, to some extent, with each other, studies on reduced and S-carboxyamidomethylated (Rcam) Naja naja atra CTX3 were carried out in the present study. Although both CTX3 and Rcam-CTX3 induced apoptotic death of PC-3 cells as evidenced by propodium iodide/annexin V double staining, degradation of procaspases and DNA fragmentation, the cytotoxicity of Rcam-CTX3 was mostly 100-fold lower than that noted with native toxin. However, Rcam-CTX3 retained approximately 38% of the membrane-damaging activity of native toxin as revealed by the decrease in calcein self-quenching from phospholipid vesicles. These results are likely to reflect that the mechanism of cytotoxicity by CTX3 is not heavily dependent on its membrane-perturbing effect, and suggest that the structural elements within CTX3 responsible for the two activities are probably separated.

Plethodontid modulating factor, a hypervariable salamander courtship pheromone in the three-finger protein superfamily

The soluble members of the three-finger protein superfamily all share a relatively simple 'three-finger' structure, yet perform radically different functions. Plethodontid modulating factor (PMF), a pheromone protein produced by the lungless salamander, Plethodon shermani, is a new and unusual member of this group. It affects female receptivity when delivered to the female's nares during courtship. As with other plethodontid pheromone genes, PMF is hyperexpressed in a specialized male mental (chin) gland. Unlike other plethodontid pheromone genes, however, PMF is also expressed at low levels in the skin, liver, intestine and kidneys of both sexes. The PMF sequences obtained from all tissue types were highly variable, with 103 unique haplotypes identified which averaged 35% sequence dissimilarity (range 1-60%) at the protein level. Despite this variation, however, all PMF sequences contained a conserved approximately 20-amino-acid secretion signal sequence and a pattern of eight cysteines that is also found in cytotoxins and short neurotoxins from snake venoms, as well as xenoxins from Xenopus. Although they share a common cysteine pattern, PMF isoforms differ from other three-finger proteins in: (a) amino-acid composition outside of the conserved motif; (b) length of the three distinguishing 'fingers'; (c) net charge at neutral pH. Whereas most three-finger proteins have a net positive charge at pH 7.0, PMF has a high net negative charge at neutral pH (pI range of most PMFs 3.5-4.0). Sequence comparisons suggest that PMF belongs to a distinct multigene subfamily within the three-finger protein superfamily.

Purification, partial characterization, crystallization and preliminary X-ray diffraction of a novel cardiotoxin-like basic protein from Naja naja atra (South Anhui) venom

A novel cardiotoxin-like basic protein was isolated from the venom of the Chinese cobra (Naja naja atra) from the south of Anhui in China. The protein inhibits the expression of vascular endothelial growth factor and basic fibroblast growth factor in human lung cancer cell line H1299 and induces the haemolysis of rabbit erythrocytes under low-lecithin conditions. After a two-step chromatographic purification, the resultant 7 kDa protein was crystallized by the hanging-drop vapour-diffusion method at room temperature. A complete data set was collected to 2.35 A resolution using an in-house X-ray diffraction system. The crystal belongs to space group P4(1)2(1)2, with unit-cell parameters a = b = 43.2, c = 147.9 A. There are two molecules in the crystallographic asymmetric unit.

Purification and characterization of Ophiophagus hannah cytotoxin-like proteins

Three cytotoxin-like proteins from the venom of Ophiophagus hannah were isolated by a combination of ion exchange chromatography and reverse phase HPLC. Amino acid sequence analysis revealed that these proteins all consisted of 63 amino acids and shared approximate 50% and 56% sequence identity with Naja naja atra cardiotoxins and cardiotoxin-like basic proteins (CLBPs), respectively. CD spectra revealed that their secondary structure was dominated with beta-sheet as those noted with cardiotoxins and CLBPs. O. hannah cytotoxin-like protein exhibited a cell-lytic activity on SK-N-SH cells, but its activity was more weak than that noted for N. naja atra cardiotoxin 3. Alternatively, apoptotic cell death was induced by the addition of N. naja atra CLBP. Based on the sequence information with the toxin molecules, the functional residues and regions related to the differential activity with O. hannah cytotoxin-like protein, cardiotoxin and CLBP are discussed.

Chemically accurate protein structures: validation of protein NMR structures by comparison of measured and predicted pKa values

A new method is presented for evaluating the quality of protein structures obtained by NMR. This method exploits the dependence between measurable chemical properties of a protein, namely pKa values of acidic residues, and protein structure. The accurate and fast empirical computational method employed by the PROPKA program ( http://www.propka.chem.uiowa.edu) allows the user to test the ability of a given structure to reproduce known pKa values, which in turn can be used as a criterion for the selection of more accurate structures. We demonstrate the feasibility of this novel idea for a series of proteins for which both NMR and X-ray structures, as well as pKa values of all ionizable residues, have been determined. For the 17 NMR ensembles used in this study, this criterion is shown effective in the elimination of a large number of NMR structure ensemble members.

Very fast empirical prediction and rationalization of protein pKa values

A very fast empirical method is presented for structure-based protein pKa prediction and rationalization. The desolvation effects and intra-protein interactions, which cause variations in pKa values of protein ionizable groups, are empirically related to the positions and chemical nature of the groups proximate to the pKa sites. A computer program is written to automatically predict pKa values based on these empirical relationships within a couple of seconds. Unusual pKa values at buried active sites, which are among the most interesting protein pKa values, are predicted very well with the empirical method. A test on 233 carboxyl, 12 cysteine, 45 histidine, and 24 lysine pKa values in various proteins shows a root-mean-square deviation (RMSD) of 0.89 from experimental values. Removal of the 29 pKa values that are upper or lower limits results in an RMSD = 0.79 for the remaining 285 pKa values.

Interaction of three-finger toxins with phospholipid membranes: comparison of S- and P-type cytotoxins

The CTs (cytotoxins) I and II are positively charged three-finger folded proteins from venom of Naja oxiana (the Central Asian cobra). They belong to S- and P-type respectively based on Ser-28 and Pro-30 residues within a putative phospholipid bilayer binding site. Previously, we investigated the interaction of CTII with multilamellar liposomes of dipalmitoylphosphatidylglycerol by wide-line (31)P-NMR spectroscopy. To compare interactions of these proteins with phospholipids, we investigated the interaction of CTI with the multilamellar liposomes of dipalmitoylphosphatidylglycerol analogously. The effect of CTI on the chemical shielding anisotropy and deformation of the liposomes in the magnetic field was determined at different temperatures and lipid/protein ratios. It was found that both the proteins do not affect lipid organization in the gel state. In the liquid crystalline state of the bilayer they disturb lipid packing. To get insight into the interactions of the toxins with membranes, Monte Carlo simulations of CTI and CTII in the presence of the bilayer membrane were performed. It was found that both the toxins penetrate into the bilayer with the tips of all the three loops. However, the free-energy gain on membrane insertion of CTI is smaller (by approximately 7 kcal/mol; 1 kcal identical with 4.184 kJ) when compared with CTII, because of the lower hydrophobicity of the membrane-binding site of CTI. These results clearly demonstrate that the P-type cytotoxins interact with membranes stronger than those of the S-type, although the mode of the membrane insertion is similar for both the types.

Cobra cardiotoxin and phospholipase A2 as GAG-binding toxins: on the path from structure to cardiotoxicity and inflammation

Glycosaminoglycans (GAGs) represent the sulfated carbohydrate moieties of proteoglycans which occur abundantly in tissues of the cardiovascular system. Many proteins bind specifically to GAGs and perform an important role in inflammation, cell proliferation, and blood coagulation processes. Recently, in vitro GAG-binding studies of cardiotoxins (CTXs) and basic phospholipase A(2) (PLA(2)) from cobra venom established the toxins as two new families of GAG-binding proteins. In particular, discontinuous basic residues in beta-sheet CTXs may form a cationic cradle suitable for heparin binding, as in the case of fibronectin module III-13. The binding specificity of beta-sheet proteins to different GAGs can be further enhanced by involving other cationic clusters near the flexible loop of the molecule. Since the three-dimensional structures of many CTXs and PLA(2) are available, these two toxins may serve as models for the elucidation of the molecular recognition of GAG-binding proteins and also as polypeptide templates for further improvement of the binding specificity suitable for future biomedical application. Research along the line of GAG-guided toxicity of cobra venom components may help us to understand the functional role of GAGs and the action mechanism of cobra venom components in the cardiovascular system.

Structural basis of membrane-induced cardiotoxin A3 oligomerization

Cobra cardiotoxins (CTXs) have previously been shown to induce membrane fusion of vesicles formed by phospholipids such as cardiolipin or sphingomyelin. CTX can also form a pore in membrane bilayers containing a anionic lipid such as phosphatidylserine or phosphatidylglycerol. Herein, we show that the interaction of CTX with negatively charged lipids causes CTX dimerization, an important intermediate for the eventual oligomerization of CTX during the CTX-induced fusion and pore formation process. The structural basis of the lipid-induced oligomerization of CTX A3, a major CTX from Naja atra, is then illustrated by the crystal structure of CTX A3 in complex with SDS; SDS likely mimics anionic lipids of the membrane under micelle conditions at 1.9-A resolution. The crystal packing reveals distinct SDS-free and SDS-rich regions; in the latter two types of interconnecting CTX A3 dimers, D1 and D2, and several SDS molecules can be identified to stabilize D1 and D2 by simultaneously interacting with residues at each dimer interface. When the three CTXSDS complexes in the asymmetric unit are overlaid, the orientation of CTX A3 monomers relative to the SDS molecules in the crystal is strikingly similar to that of the toxin with respect to model membranes as determined by NMR and Fourier transform infrared methods. These results not only illustrate how lipid-induced CTX dimer formation may be transformed into oligomers either as inverted micelles of fusion intermediates or as membrane pore of anionic lipid bilayers but also underscore a potential role for SDS in x-ray diffraction study of protein-membrane interactions in the future.

Interaction of the P-type cardiotoxin with phospholipid membranes

The cardiotoxin (cytotoxin II, or CTII) isolated from cobra snake (Naja oxiana) venom is a 60-residue basic membrane-active protein featuring three-finger beta sheet fold. To assess possible modes of CTII/membrane interaction 31P- and 1H-NMR spectroscopy was used to study binding of the toxin and its effect onto multilamellar vesicles (MLV) composed of either zwitterionic or anionic phospholipid, dipalmitoylglycerophosphocholine (Pam2Gro-PCho) or dipalmitoylglycerophosphoglycerol (Pam2Gro-PGro), respectively. The analysis of 1H-NMR linewidths of the toxin and 31P-NMR spectral lineshapes of the phospholipid as a function of temperature, lipid-to-protein ratios, and pH values showed that at least three distinct modes of CTII interaction with membranes exist: (a) nonpenetrating mode; in the gel state of the negatively charged MLV the toxin is bound to the surface electrostatically; the binding to Pam2Gro-PCho membranes was not observed; (b) penetrating mode; hydrophobic interactions develop due to penetration of the toxin into Pam2Gro-PGro membranes in the liquid-crystalline state; it is presumed that in this mode CTII is located at the membrane/water interface deepening the side-chains of hydrophobic residues at the tips of the loops 1-3 down to the boundary between the glycerol and acyl regions of the bilayer; (c) the penetrating mode gives way to isotropic phase, stoichiometrically well-defined CTII/phospholipid complexes at CTII/lipid ratio exceeding a threshold value which was found to depend at physiological pH values upon ionization of the imidazole ring of His31. Biological implications of the observed modes of the toxin-membrane interactions are discussed.

Interaction of cardiotoxins with membranes: a molecular modeling study

Incorporation of beta-sheet proteins into membrane is studied theoretically for the first time, and the results are validated by the direct experimental data. Using Monte Carlo simulations with implicit membrane, we explore spatial structure, energetics, polarity, and mode of insertion of two cardiotoxins with different membrane-destabilizing activity. Both proteins, classified as P- and S-type cardiotoxins, are found to retain the overall "three-finger" fold interacting with membrane core and lipid/water interface by the tips of the "fingers" (loops). The insertion critically depends upon the structure, hydrophobicity, and electrostatics of certain regions. The simulations reveal apparently distinct binding modes for S- and P-type cardiotoxins via the first loop or through all three loops, respectively. This rationalizes an earlier empirical classification of cardiotoxins into S- and P-type, and provides a basis for the analysis of experimental data on their membrane affinities. Accomplished with our previous simulations of membrane alpha-helices, the computational method may be used to study partitioning of proteins with diverse folds into lipid bilayers.

Dynamic characterization of the water binding loop in the P-type cardiotoxin: implication for the role of the bound water molecule

Recent studies of cobra P-type cardiotoxins (CTXs) have shown that the water-binding loop (loop II) plays a crucial role in toxin binding to biological membranes and in their cytotoxicity. To understand the role of bound water in the loop, the structure and dynamics of the major P-type CTX from Taiwan cobra, CTX A3, were determined by a comprehensive NMR analysis involving (1)H NOESY/ROESY, (13)C[1)H]NOE/T(1) relaxation, and (17)O triple-quantum filtered NMR. A single water molecule was found to be tightly hydrogen bonded to the NH of Met26 with a correlation time (5-7 ns) approaching the isotropic tumbling time (3.8-4.5 ns) of the CTX A3 molecule. Surprisingly, despite the relatively long residence time (ca. 5 ns to 100 micros), the bound water molecule of CTX A3 is located within a dynamic (order parameter S(2) approximately 0.7) and solvent accessible loop. Comparison among several P-type CTXs suggests that proline residues in the consensus sequence of MxAxPxVPV should play an important role in the formation of the water binding loop. It is proposed that the exchange rate of the bound water may play a role in regulating the lipid binding mode of amphiphilic CTX molecules near membrane surfaces.

Membrane binding motif of the P-type cardiotoxin

Carditoxins (CTXs) from cobra snake venoms, the basic 60-62 residue all-beta sheet polypeptides, are known to bind to and impair the function of cell membranes. To assess the membrane induced conformation and orientation of CTXs, the interaction of the P-type cardiotoxin II from Naja oxiana snake venom (CTII) with perdeuterated dodecylphosphocholine (DPC) was studied using ( 1 )H-NMR spectroscopy and diffusion measurements. Under conditions where the toxin formed a well-defined complex with DPC, the spatial structure of CTII with respect to the presence of tightly bound water molecules in loop II, was calculated using the torsion angle dynamics program DYANA. The structure was found to be similar, except for subtle changes in the tips of all three loops, to the previously described "major" form of CTII in aqueous solution illustrated by the "trans" configuration of the Val7-Pro8 peptide bond. No "minor" form with the "cis" configuration of the above bond was found in the micelle-bound state. The broadening of the CTII backbone proton signals by 5, 16-doxylstearate relaxation probes, together with modeling based on the spatial structure of CTII, indicated a periphery mode of binding of the toxin molecule to the micelle and revealed its micelle interacting domain. The latter includes a hydrophobic region of CTII within the extremities of loops I and III (residues 5-11, 46-50), the basement of loop II (residues 24-29,31-37) and the belt of polar residues encircling these loops (lysines 4,5,12,23,50, serines 11,46, histidine 31, arginine 36). It is suggested that this structural motif and the mode of binding can be realized during interaction of CTXs with lipid and biological membranes.

The multiplicity of cardiotoxins from Naja naja atra (Taiwan cobra) venom

Four novel cardiotoxins were isolated from Naja naja atra (Taiwan cobra) venom by successive separation on a SP-Sephadex C-25 column and a reverse phase column. Amino acid sequences of the cardiotoxins were determined by Edman degradation and carboxypeptidase digestion. It shows that these cardiotoxins comprise 60 amino acid residues. Comparative analyses on the amino acid sequences of cardiotoxins from the venoms of N. naja atra and other Naja species indicated that amino acid substitutions of cardiotoxin isoforms frequently occurred at positions 7-11, 27-32 and 45-47. The hypervariable segments encoded by the second and third exon of cardiotoxin genes are located at or near the tips of loop structure of cardiotoxin molecules. These results, together with the suggestions that the residues at the tips of cardiotoxins' loop structure were involved in the manifestation of the biological activities of cardiotoxins, reflect that the preferential mutations may contribute to alterations in the function of cardiotoxin molecules. Analysis on the secondary structure of pre-mRNAs of N. naja atra cardiotoxin 4 gene and N. naja sputatrix cardiotoxin 3 gene has shown that the hypervariable regions of the exon 2 pertain to form intra-exon pairings and are not involved in the formation of intron-exon pairings. Since the pairings of splice sites and gene architecture were supposed to be associated with intron-exon recognition, it is likely that the preferred loci of mutations occurring with the evolution of cardiotoxin genes would not affect the processing of cardiotoxin precursors.

Structure-function studies on Taiwan cobra long neurotoxin homolog

A novel long neurotoxin homolog was purified from Naja naja atra (Taiwan cobra) venom using the combination of ion exchange chromatography and reverse phase high performance liquid chromatography. The determined protein sequence was essentially the same as that deduced from the cDNA amplified by reverse transcriptase-polymerase chain reaction. The long neurotoxin homolog exhibited an activity that inhibited acetylcholine-induced muscle contractions, as with N. naja atra cobrotoxin. The degree of inhibition caused by the addition of long neurotoxin homolog was approximately 70% of that observed with the addition of cobrotoxin. Unlike the well-known short and long neurotoxins, this neurotoxin homolog contained two additional cysteine residues forming a disulfide linkage in the N-terminal region. Circular dichroism measurement and computer models of the neurotoxin reveal that its secondary structure was not abundant in beta-sheet as noted with short and long neurotoxins. This less ordered structure may be associated with the lower activity noted with the long neurotoxin homolog. Together with the finding that the known long neurotoxin homologs exclusively appear in the venoms of the Naja and Bungarus genera, the long neurotoxin homologs should represent an evolutionary branch from the long and short neurotoxins in the Elapidae family.

Two forms of cytotoxin II (cardiotoxin) from Naja naja oxiana in aqueous solution: spatial structures with tightly bound water molecules

1H-NMR spectroscopy data, such as NOE intraprotein and (bound water)/protein contacts, 3J coupling constants and deuterium exchange rates were used to determine the in-solution spatial structure of cytotoxin II from Naja naja oxiana snake venom (CTII). Exploiting information from two 1H-NMR spectral components, shown to be due to cis/trans isomerization of the Val7-Pro8 peptide bond, spatial structures of CTII minor and major forms (1 : 6) were calculated using the torsion angle dynamics algorithm of the DYANA program and then energy refined using the FANTOM program. Each form, major and minor, is represented by 20 resulting conformers, demonstrating mean backbone rmsd values of 0.51 and 0.71 A, respectively. Two forms of CTII preserve the structural skeleton as three large loops, including two beta-sheets with bend regions, and demonstrate structural differences at loop I, where cis/trans isomerization occurs. The CTII side-chain distribution constitutes hydrophilic and hydrophobic belts around the protein, alternating in the trend of the three main loops. Because of the Omega-shaped backbone, formed in participation with two bound water molecules, the tip of loop II bridges the tips of loops I and III. This ensures the continuity of the largest hydrophobic belt, formed with the residues of these tips. Comparison revealed pronounced differences in the spatial organization of the tips of the three main loops between CTII and previous structures of homologous cytotoxins (cardiotoxins) in solution.

Identification of specific sites involved in ligand binding by photoaffinity labeling of the receptor for the urokinase-type plasminogen activator. Residues located at equivalent positions in uPAR domains I and III participate in the assembly of a composite ligand-binding site

Plasminogen activation by the urokinase-type plasminogen activator (uPA) is facilitated in the presence of cells expressing the glycolipid-anchored high-affinity receptor for uPA (denoted uPAR). Structures involved in the interaction between human uPAR and a decamer peptide antagonist of uPA binding (SLNFSQYLWS) were previously tagged by specific site-directed photoaffinity labeling [Ploug, M., Ostergaard, S., Hansen, L. B. L., Holm, A., and Dano, K. (1998) Biochemistry 37, 3612-3622]. Replacement of the key functional residues Phe4 and Trp9 with either benzophenone or (trifluoromethyl)aryldiazirine rendered this peptide antagonist photoactivatable, and as a consequence, it incorporated covalently upon photolysis into either uPAR domain I or domain III depending on the actual position of the photophore in the sequence. The residues of uPAR specifically targeted by photoaffinity labeling were identified by matrix-assisted laser desorption mass spectrometry, NH2-terminal sequence analysis, and amino acid composition analysis after enzymatic fragmentation and HPLC purification. According to these data, the formation of the receptor-ligand complex positions Phe4 of the peptide antagonist very close to Arg53 and Leu66 in uPAR domain I and Trp9 of the antagonist in the vicinity of His251 in uPAR domain III. The gross molecular arrangement of the deduced receptor-ligand interface provides a rational structural basis for the observed requirement for the intact multidomain state of uPAR for achieving high-affinity ligand binding, since according to this model ligand binding must rely on a close spatial proximity of uPAR domains I and III. In addition, these data suggest that the assembly of the composite ligand binding site in uPAR may resemble the homophilic interdomain dimerization of kappa-bungarotoxin, a structural homologue of the Ly-6/uPAR domain family.

Effect of D57N mutation on membrane activity and molecular unfolding of cobra cardiotoxin

Cobra cardiotoxins (CTXs) are able to adopt a three-fingered beta-strand structure with continuous hydrophobic patch that is capable of interacting with zwitterionic phospholipid bilayer. In addition to the four disulfide bonds that form the rigid core of CTXs, Asp57 near the C-terminus interacts electrostatically with Lys2 near the N-terminus (Chiang et al. 1996. Biochemistry. 35:9177-9186). We indicate herein, using circular dichroism and the time-resolved polarized tryptophan fluorescence measurement, that Asp57 to Asn57 (D57N) mutation perturbs the structure of CTX molecules at neutral pH. The structural stability of the D57N mutant was found to be lower, as evidenced by the reduced effective concentration of the 2,2,2-trifluoethanol (TFE)-induced beta-sheet to alpha-helix transition. Interestingly, the single mutation also allows a greater degree of molecular unfolding, because the rotational correlation time of the TFE-induced unfolding intermediate is larger for the D57N mutant. It is suggested that the electrostatic interaction between N- and C-termini also contributes to the formation of the functionally important continuous hydrophobic stretch on the distant end of CTX molecules, because both the binding to anilinonaphthalene fluorescent probe and the interaction with phospholipid bilayer were also reduced for D57N mutant. The result emphasizes the importance of the hydrophobic amino acid residues near the tip of loop 3 as a continuous part of the three-fingered beta-strand CTX molecule and indicates how a distant electrostatic interaction might be involved. It is also implicated that electrostatic interaction plays a role in expanding the radius of gyration of the folding/unfolding intermediate of proteins.

Glycosaminoglycans bind to homologous cardiotoxins with different specificity

We have reported that some cardiotoxins (CTXs), homologous basic polypeptides of cobra venom, bind strongly to heparin. Herein we show that CTXs from spitting cobra venom bind avidly to chondroitin-6-sulfate (CS6) and dermatan sulfate (DS), the glycosaminoglycans (GAGs) abounding in the cornea and elsewhere. We compared the binding strength of Tgamma, a major component of spitting cobra, Naja nigricollis, venom with that of CTX A3, a major component of Naja atra venom to various GAGs including CS6, chondroitin-4-sulfate (CS4), DS, keratan sulfate (KS), hyaluronan (HYA), and heparin. The binding strength of Tgamma followed the order CS6 > KS > HYA > DS > CS4 > heparin, whereas that of CTX A3 was heparin > KS > CS4 > DS > CS6 > HYA. The binding specificity displayed by different CTXs toward GAGs is impressive, given the high homology among CTXs and among GAGs. Strong binding of Tgamma to CS6, rather than to the highly anionic and versatile cousin, heparin, implies specific interaction with CS6. Heparin, at high concentration, displaced CS6 from CS6-Tgamma and CS6-A3 complexes. We also show that corneal CS/DS likely allow Tgamma to bind to corneal epithelium. CTXs of spitting cobra venom are known to cause corneal opacity and/or blindness. Taken together with these observations, our results suggest that corneal CS/DS play a role in the action of CTX in the eye. Most importantly, the present results establish CTXs as cationic, readily available, avidly binding ligands of CS/DS.

Photoaffinity labeling of the human receptor for urokinase-type plasminogen activator using a decapeptide antagonist. Evidence for a composite ligand-binding site and a short interdomain separation

Binding of urokinase-type plasminogen activator (uPA) to its cellular receptor (uPAR) renders the cell surface a favored site for plasminogen activation. Recently, a 15-mer peptide antagonist of the uPA-uPAR interaction, with an IC50 value of 10 nM, was identified using phage display technology [Goodson, R. J., Doyle, M. V., Kaufman, S. E., and Rosenberg, S. (1994) Proc. Natl. Acad. Sci. 91, 7129-7133]. In the present study, the molecular aspects of the interaction between this peptide and uPAR have been investigated. We have characterized the real-time receptor binding kinetics for the antagonist using surface plasmon resonance and identified critical residues by alanine replacements. The minimal peptide antagonist thus derived (SLNFSQYLWS) was rendered photoactivatable by replacing residues important for uPAR binding with photochemically active derivatives of phenylalanine containing either (trifluoromethyl)diazirine or benzophenone. These peptides incorporated covalently into purified soluble uPAR upon photoactivation, and this was inhibited by preincubation with receptor binding derivatives of uPA. The intact three-domain structure of uPAR was essential for efficient photoaffinity labeling. Proteolytic domain mapping using chymotrypsin revealed a specific labeling of both uPAR domain I and domains II + III dependent on the position of the photoprobe in the antagonist. On the basis of these studies, we propose the existence of a composite ligand binding site in uPAR combined of residues located in distinct structural domains. According to this model, a close spatial proximity between uPAR domain I and either domains II or III in intact uPAR is required for the assembly of this composite binding site. Since the receptor binding properties of the peptide antagonist closely mimic those of uPA itself, these two ligands presumably share coincident binding site in uPAR.

Genomic structures of cardiotoxin 4 and cobrotoxin from Naja naja atra (Taiwan cobra)

Two genomic DNAs with the size of 2.3 kb and 2.4 kb, which were isolated from the liver of Naja naja atra (Taiwan cobra), encoded the precursors of cardiotoxin 4 and cobrotoxin, respectively. Both genes shared virtually identical overall organization with three exons separated by two introns, which were inserted in the similar positions of the gene's coding regions. Moreover, their nucleotide sequences shared approximately 84.2% identity. This result reveals the evolutionary relationship between cardiotoxin and cobrotoxin. The exon/intron structures of cardiotoxin 4 and cobrotoxin genes were similar to that reported for erabutoxin c gene, a neurotoxin genomic DNA from a sea snake (Laticauda semifasciata). However, in contrast to the finding that the intron 2 of these genes had a similar size, a notable variation with the size of intron 1 was observed (1233 bp, 1269 bp and 197 bp for cardiotoxin 4, cobrotoxin and erabutoxin c genes, respectively). The different size with intron 1 is due to the middle region at the first intron of cardiotoxin 4 and cobrotoxin genes, which encoded small nucleolar RNA (snoRNA), being absent in that of erabutoxin c gene. These results, together with the finding of the potential mobility of snoRNA genes during evolution, suggest that intron insertions or deletions of snoRNA genes occur with the evolutionary divergence of snake neurotoxins and cardiotoxins.

Analysis of binding of cobra cardiotoxins to heparin reveals a new beta-sheet heparin-binding structural motif

Heparin and heparan sulfate have recently been shown to bind to snake cardiotoxin (CTX) and to potentiate its penetration into phospholipid monolayer under physiological ionic conditions. Herein we analyze the heparin-binding domain of CTX using 10 CTXs from Taiwan and African cobra venom. We also performed computer modeling to obtain more information of the binding at molecular level. The results provide a molecular model for interaction of CTX-heparin complex where the cationic belt of the conserved residues on the concave surface of three finger beta-sheet polypeptides initiates ionic interaction with heparin-like molecules followed by specific binding of Lys residues near the tip of loop 2 of CTX. The dissociation constants of CTXs differ by as much as 4 orders of magnitude, ranging from approximately 140 microM for toxin gamma to approximately 20 nM for CTX M3, depending on the presence of Lys residues near the tip of loop 2. High affinity heparin binding becomes possible due to the presence of Arg-28, Lys-33, or the so-called consensus heparin binding sequence of XKKXXXKRX near the tip of the loop. The well defined three-finger loop structure of CTX provides an interesting template for the design of high affinity heparin-binding polypeptides with beta-sheet structure. The finding that several cobra CTXs and phospholipase A2 bind to heparin with different affinity may provide information on the synergistic action of the two venom proteins.