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

Specific Amino Acid Residues in the Three Loops of Snake Cytotoxins Determine Their Membrane Activity and Provide a Rationale for a New Classification of These Toxins

Cytotoxins (CTs) are three-finger membrane-active toxins present mainly in cobra venom. Our analysis of the available CT amino acid sequences, literature data on their membrane activity, and conformational equilibria in aqueous solution and detergent micelles allowed us to identify specific amino acid residues which interfere with CT incorporation into membranes. They include Pro9, Ser28, and Asn/Asp45 within the N-terminal, central, and C-terminal loops, respectively. There is a hierarchy in the effect of these residues on membrane activity: Pro9 > Ser28 > Asn/Asp45. Taking into account all the possible combinations of special residues, we propose to divide CTs into eight groups. Group 1 includes toxins containing all of the above residues. Their representatives demonstrated the lowest membrane activity. Group 8 combines CTs that lack these residues. For the toxins from this group, the greatest membrane activity was observed. We predict that when solely membrane activity determines the cytotoxic effects, the activity of CTs from a group with a higher number should exceed that of CTs from a group with a lower number. This classification is supported by the available data on the cytotoxicity and membranotropic properties of CTs. We hypothesize that the special amino acid residues within the loops of the CT molecule may indicate their involvement in the interaction with non-lipid targets.

Bilinear optimization of protein structure prediction: An exact approach via AB off-lattice model

Protein structure prediction (PSP) remains a central challenge in computational biology due to its inherent complexity and high dimensionality. While numerous heuristic approaches have appeared in the literature, their success varies. The AB off-lattice model, which characterizes proteins as sequences of A (hydrophobic) and B (hydrophilic) beads, presents a simplified perspective on PSP. This work presents a mathematical optimization-based methodology capitalizing on the off-lattice AB model. Dissecting the inherent non-linearities of the energy landscape of protein folding allowed for formulating the PSP as a bilinear optimization problem. This formulation was achieved by introducing auxiliary variables and constraints that encapsulate the nuanced relationship between the protein's conformational space and its energy landscape. The proposed bilinear model exhibited notable accuracy in pinpointing the global minimum energy conformations on a benchmark dataset presented by the Protein Data Bank (PDB). Compared to traditional heuristic-based methods, this bilinear approach yielded exact solutions, reducing the likelihood of local minima entrapment. This research highlights the potential of reframing the traditionally non-linear protein structure prediction problem into a bilinear optimization problem through the off-lattice AB model. Such a transformation offers a route toward methodologies that can determine the global solution, challenging current PSP paradigms. Exploration into hybrid models, merging bilinear optimization and heuristic components, might present an avenue for balancing accuracy with computational efficiency.

The application of WaterMap-guided structure-based virtual screening in novel drug discovery

INTRODUCTION

Nowadays, it is widely accepted that water molecules play a key role in binding a ligand to a molecular target. Neglecting water molecules in the process of molecular recognition was the result of several failures of the structure-based drug discovery campaigns. The application of WaterMap, in particular WaterMap-guided molecular docking, enables the reasonably accurate and quick description of the location and energetics of water molecules at the ligand-protein interface.

AREAS COVERED

In this review, the authors shortly discuss the importance of water in drug design and discovery and provide a brief overview of the computational approaches used to predict the solvent-related effects for the purposes of presenting WaterMap in the context of other available techniques and tools. A concise description of WaterMap concept is followed by the presentation of WaterMap-assisted virtual screening literature published between 2013 and 2023.

EXPERT OPINION

In recent years, WaterMap software has been extensively used to support structure-based drug design, in particular structure-based virtual screening. Indeed, it is a useful tool to rescore docking results considering water molecules in the binding pocket. Although WaterMap allows for the consideration of the dynamic behavior of water molecules in the binding site, for best accuracy, its application in conjunction with other techniques such as molecular mechanics-generalized Born surface area of FEP (Free Energy Perturbation) is recommended.

Fifty Years of Animal Toxin Research at the Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry RAS

This review covers briefly the work carried out at our institute (IBCh), in many cases in collaboration with other Russian and foreign laboratories, for the last 50 years. It discusses the discoveries and studies of various animal toxins, including protein and peptide neurotoxins acting on the nicotinic acetylcholine receptors (nAChRs) and on other ion channels. Among the achievements are the determination of the primary structures of the α-bungarotoxin-like three-finger toxins (TFTs), covalently bound dimeric TFTs, glycosylated cytotoxin, inhibitory cystine knot toxins (ICK), modular ICKs, and such giant molecules as latrotoxins and peptide neurotoxins from the snake, as well as from other animal venoms. For a number of toxins, spatial structures were determined, mostly by 1H-NMR spectroscopy. Using this method in combination with molecular modeling, the molecular mechanisms of the interactions of several toxins with lipid membranes were established. In more detail are presented the results of recent years, among which are the discovery of α-bungarotoxin analogs distinguishing the two binding sites in the muscle-type nAChR, long-chain α-neurotoxins interacting with α9α10 nAChRs and with GABA-A receptors, and the strong antiviral effects of dimeric phospholipases A2. A summary of the toxins obtained from arthropod venoms includes only highly cited works describing the molecules' success story, which is associated with IBCh. In marine animals, versatile toxins in terms of structure and molecular targets were discovered, and careful work on α-conotoxins differing in specificity for individual nAChR subtypes gave information about their binding sites.

Membrane-Disrupting Activity of Cobra Cytotoxins Is Determined by Configuration of the N-Terminal Loop

In aqueous solutions, cobra cytotoxins (CTX), three-finger folded proteins, exhibit conformational equilibrium between conformers with either cis or trans peptide bonds in the N-terminal loop (loop-I). The equilibrium is shifted to the cis form in toxins with a pair of adjacent Pro residues in this loop. It is known that CTX with a single Pro residue in loop-I and a cis peptide bond do not interact with lipid membranes. Thus, if a cis peptide bond is present in loop-I, as in a Pro-Pro containing CTX, this should weaken its lipid interactions and likely cytotoxic activities. To test this, we have isolated seven CTX from Naja naja and N. haje cobra venoms. Antibacterial and cytotoxic activities of these CTX, as well as their capability to induce calcein leakage from phospholipid liposomes, were evaluated. We have found that CTX with a Pro-Pro peptide bond indeed exhibit attenuated membrane-perturbing activity in model membranes and lower cytotoxic/antibacterial activity compared to their counterparts with a single Pro residue in loop-I.

Toxins' classification through Raman spectroscopy with principal component analysis

The method based on the combination of Raman spectroscopy and principal component analysis (PCA) was applied to the set of peptide and protein toxins from animal venoms and to synthetic analogues of peptides. The study demonstrated the possibility of toxin classification according to their primary and secondary structures based on Raman spectroscopy. The method described here allows discrimination of snake venom three-finger toxins from predatory marine mollusks α-conotoxins. Moreover, PCA of the Raman spectra of toxins revealed differences within the group of three-finger toxins and also within the group of conotoxins, related to their spatial structure. In particular, on the basis of the developed technique it is possible to distinguish the disulfide isomers of the same peptide toxin. The results obtained have been confirmed by bioinformatic methods. So, we have proposed a method for the rapid analysis of newly discovered venom-derived protein or peptide toxins by establishing their similarity with other already studied toxins by referring to a particular class. Taking into account a low specimen consumption by Raman spectroscopy, the proposed method could represent a first step in the study of toxins from rare and/or endangered venomous animals. The ability to distinguish configuration of disulfide bonds allows to synthesize the correct isomer of the toxin.

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.

The omega-loop of cobra cytotoxins tolerates multiple amino acid substitutions

Cobra cytotoxins (CTs), the three-fingered proteins, feature high amino acid sequence homology in the beta-strands and variations in the loop regions. We selected a pair of cytotoxins from Naja kaouthia crude venom to clarify the sequence-structure relationships. Using chromatography and mass spectroscopy, we separated and identified the mixture of cytotoxins 2 and 3, differentiated by the only Val 41/Ala 41 substitution. Here, using natural abundance ¹³C, ¹⁵N NMR-spectroscopy we performed chemical shift assignments of the signals of the both toxins in aqueous solution in the major and minor forms. Combining NOE and chemical shift data, the toxins' spatial structure was determined. Finally, we proved that the tip of the "finger"-2, or the loop-2 of cytotoxins adopts the shape of an omega-loop with a tightly-bound water molecule in its cavity. Comparison with other NMR and X-ray structures of cytotoxins possessing different amino acid sequences reveals spatial similarity in this family of proteins, including the loop-2 region, previously considered to be flexible.

Stepwise Insertion of Cobra Cardiotoxin CT2 into a Lipid Bilayer Occurs as an Interplay of Protein and Membrane "Dynamic Molecular Portraits"

For many peripheral membrane-binding polypeptides(MBPs), especially β-structural ones, the precise molecular mechanisms of membrane insertion remain unclear. In most cases, only the terminal water-soluble and membrane-bound states have been elucidated, whereas potential functionally important intermediate stages are still not understood in sufficient detail. In this study, we present one of the first successful attempts to describe step-by-step embedding of the MBP cardiotoxin 2 (CT2) from cobra Naja oxiana venom into a lipid bilayer at the atomistic level. CT2 possesses a highly conservative and rigid β-structured three-finger fold shared by many other exogenous and endogenous proteins performing a wide variety of functions. The incorporation of CT2 into the lipid bilayer was analyzed via a 2 μs all-atom molecular dynamics (MD) simulation without restraints. This process was shown to occur over a number of distinct steps, while the geometry of initial membrane attachment drastically differs from that of the final equilibrated state. In the latter one, the hydrophobic platform ("bottom") formed by the tips of the three loops is deeply buried into the lipid bilayer. This agrees well with the NMR data obtained earlier for CT2 in detergent micelles. However, the bottom is too bulky to insert itself into the membrane at once. Instead, the gradual immersion of CT2 initiated by the loop-1 was observed. This initial binding stage was also demonstrated in a series of MD runs with varying starting orientations of the toxin with respect to the bilayer surface. Apart from the nonspecific long-range electrostatic attraction and hydrophobic match/mismatch factor, several specific lipid-binding sites were identified in CT2. They were shown to promote membrane insertion by engaging in strong interactions with lipid head groups, fine-tuning the toxin-membrane accommodation. We therefore propose that the toxin insertion relies on the interplay of nonspecific and specific interactions, which are determined by the "dynamic molecular portraits" of the two players, the protein and the membrane. The proposed model does not require protein oligomerization for membrane insertion and can be further employed to design MBPs with predetermined properties with regard to particular membrane targets.

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.

Naja naja oxiana Cobra Venom Cytotoxins CTI and CTII Disrupt Mitochondrial Membrane Integrity: Implications for Basic Three-Fingered Cytotoxins

Cobra venom cytotoxins are basic three-fingered, amphipathic, non-enzymatic proteins that constitute a major fraction of cobra venom. While cytotoxins cause mitochondrial dysfunction in different cell types, the mechanisms by which cytotoxins bind to mitochondria remain unknown. We analyzed the abilities of CTI and CTII, S-type and P-type cytotoxins from Naja naja oxiana respectively, to associate with isolated mitochondrial fractions or with model membranes that simulate the mitochondrial lipid environment by using a myriad of biophysical techniques. Phosphorus-31 nuclear magnetic resonance (³¹P-NMR) spectroscopy data suggest that both cytotoxins bind to isolated mitochondrial fractions and promote the formation of aberrant non-bilayer structures. We then hypothesized that CTI and CTII bind to cardiolipin (CL) to disrupt mitochondrial membranes. Collectively, ³¹P-NMR, electron paramagnetic resonance (EPR), proton NMR (¹H-NMR), deuterium NMR (²H-NMR) spectroscopy, differential scanning calorimetry, and erythrosine phosphorescence assays suggest that CTI and CTII bind to CL to generate non-bilayer structures and promote the permeabilization, dehydration and fusion of large unilamellar phosphatidylcholine (PC) liposomes enriched with CL. On the other hand, CTII but not CTI caused biophysical alterations of large unilamellar PC liposomes enriched with phosphatidylserine (PS). Mechanistically, single molecule docking simulations identified putative CL, PS and PC binding sites in CTI and CTII. While the predicted binding sites for PS and PC share a high number of interactive amino acid residues in CTI and CTII, the CL biding sites in CTII and CTI are more divergent as it contains additional interactive amino acid residues. Overall, our data suggest that cytotoxins physically associate with mitochondrial membranes by binding to CL to disrupt mitochondrial structural integrity.

Anticancer Activity of Cobra Venom Polypeptide, Cytotoxin-II, against Human Breast Adenocarcinoma Cell Line (MCF-7) via the Induction of Apoptosis

Purpose

Breast cancer is a significant health problem worldwide, accounting for a quarter of all cancer diagnoses in women. Current strategies for breast cancer treatment are not fully effective, and there is substantial interest in the identification of novel anticancer agents especially from natural products including toxins. Cytotoxins are polypeptides found in the venom of cobras and have various physiological effects. In the present study, the anticancer potential of cytotoxin-II against the human breast adenocarcinoma cell line (MCF-7) was investigated.

Methods

The cytotoxic effects of cytotoxin-II were determined by morphological analysis and 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide assay. The mode and mechanism of cell death were investigated via acridine orange/ethidium bromide (AO/EtBr) double staining, flow cytometric analysis of cell death, detection of mitochondrial membrane potential, measurement of intracellular reactive oxygen species (ROS), annexin V/propidium iodide staining, and caspase-9 activity assays.

Results

The half maximal inhibitory concentration (IC₅₀) of cytotoxin-II in MCF-7 cells was 4.18±1.23 µg/mL, while the value for cisplatin was approximately 28.02±1.87 µg/mL. Morphological analysis and AO/EtBr double staining showed typical manifestations of apoptotic cell death (in doses lower than 8 µg/mL). Dose- and time-dependent ROS generation, loss of mitochondrial membrane potential, caspase-9 activation, and cell cycle arrest were observed in their respective tests.

Conclusion

In conclusion, cytotoxin-II has potent anticancer effects in the MCF-7 cell line, which are induced via the intrinsic pathways of apoptosis. Based on these findings, cytotoxin-II is a suitable choice for breast cancer treatment.

NOT THAT RIGID MIDGETS AND NOT SO FLEXIBLE GIANTS: ON THE ABUNDANCE AND ROLES OF INTRINSIC DISORDER IN SHORT AND LONG PROTEINS

Intrinsically disordered proteins or proteins with disordered regions are very common in nature. These proteins have numerous biological functions which are complementary to the biological activities of traditional ordered proteins. A noticeable difference in the amino acid sequences encoding long and short disordered regions was found and this difference was used in the development of length-dependent predictors of intrinsic disorder. In this study, we analyze the scaling of intrinsic disorder in eukaryotic proteins and investigate the presence of length-dependent functions attributed to proteins containing long disordered regions.

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.

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.

Structural Difference between Group I and Group II Cobra Cardiotoxins: X-ray, NMR, and CD Analysis of the Effect ofcis-Proline Conformation on Three-Fingered Toxins†

Natural homologues of cobra cardiotoxins (CTXs) were classified into two structural subclasses of group I and II based on the amino acid sequence and circular dichroism analysis, but the exact differences in their three-dimensional structures and biological significance remain elusive. We show by circular dichroism, NMR spectroscopic, and X-ray crystallographic analyses of a newly purified group I CTX A6 from eastern Taiwan cobra (Naja atra) venoms that its loop I conformation adopts a type VIa turn with a cis peptide bond located between two proline residues of PPxY. A similar "banana-twisted" conformation can be observed in other group I CTXs and also in cyclolinopeptide A and its analogues. By binding to the membrane environment, group I CTX undergoes a conformational change to adopt a more extended hydrophobic domain with beta-sheet twisting closer to the one adopted by group II CTX. This result resolves a discrepancy in the CTX structural difference reported previously between solution as well as crystal state and shows that, in addition to the hydrophobicity, the exact loop I conformation also plays an important role in CTX-membrane interaction. Potential protein targets of group I CTXs after cell internalization are also discussed on the basis of the determined loop I conformation.

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.

Spatial structure of zervamicin IIB bound to DPC micelles: implications for voltage-gating

Zervamicin IIB is a 16-amino acid peptaibol that forms voltage-dependent ion channels with multilevel conductance states in planar lipid bilayers and vesicular systems. The spatial structure of zervamicin IIB bound to dodecylphosphocholine micelles was studied by nuclear magnetic resonance spectroscopy. The set of 20 structures obtained has a bent helical conformation with a mean backbone root mean square deviation value of approximately 0.2 A and resembles the structure in isotropic solvents (Balashova et al., 2000. NMR structure of the channel-former zervamicin IIB in isotropic solvents. FEBS Lett 466:333-336). The N-terminus represents an alpha-helix, whereas the C-terminal part has a mixed 3(10)/alpha(R) hydrogen-bond pattern. In the anisotropic micelle environment, the bending angle on Hyp10 (23 degrees) is smaller than that (47 degrees) in isotropic solvents. In the NOESY (Nuclear Overhauser Effect Spectroscopy) spectra, the characteristic attenuation of the peptide signals by 5- and 16-doxylstearate relaxation probes indicates a peripheral mode of the peptaibol binding to the micelle with the N-terminus immersed slightly deeper into micelle interior. Analysis of the surface hydrophobicity reveals that the zervamicin IIB helix is amphiphilic and well suited to formation of a tetrameric transmembrane bundle, according to the barrel-stave mechanism. The results are discussed in a context of voltage-driven peptaibol insertion into membrane.

Heparin binding stabilizes the membrane-bound form of cobra cardiotoxin

It has been shown previously that the long chain fragments of heparin bind to the beta-strand cationic belt of the three-finger cobra cardiotoxin (or cytotoxin, CTX) and hence enhance its penetration into phospholipid monolayer under physiological ionic conditions. By taking lysophosphatidylcholine (LPC) micelles as a membrane model, we have shown by (1)H NMR study that the binding of heparin-derived hexasaccharide (Hep-6) to CTX at the beta-strand region can induce conformational changes of CTX near its membrane binding loops and promote the binding activity of CTX toward LPC. The Fourier-transform infrared spectra and NMR nuclear Overhauser effect of Hep-6.CTX and CTX.LPC complex in aqueous buffer also supplemented the aforementioned observation. Thus, the detected conformational change may presumably be the result of structural coupling between the connecting loops and its beta-strands. This is the first documentation of results showing how the association of hydrophilic carbohydrate molecules with amphiphilic proteins can promote hydrophobic protein-lipid interaction via the stabilization of its membrane-bound form. A similar mechanism involving tripartite interactions of heparin, protein, and lipid molecules may be operative near the extracellular matrix of cell membranes.

NMR structure of bucandin, a neurotoxin from the venom of the Malayan krait (Bungarus candidus)

A high-resolution solution structure of bucandin, a neurotoxin from Malayan krait (Bungarus candidus), was determined by (1)H-NMR spectroscopy and molecular dynamics. The average backbone root-mean-square deviation for the 20 calculated structures and the mean structure is 0.47 A (1 A=0.1 nm) for all residues and 0.24 A for the well-defined region that spans residues 23-58. Secondary-structural elements include two antiparallel beta-sheets characterized by two and four strands. According to recent X-ray analysis, bucandin adopts a typical three-finger loop motif and yet it has some peculiar characteristics that set it apart from other common alpha-neurotoxins. The presence of a fourth strand in the second antiparallel beta-sheet had not been observed before in three-finger toxins, and this feature was well represented in the NMR structure. Although the overall fold of the NMR structure is similar to that of the X-ray crystal structure, there are significant differences between the two structures that have implications for the pharmacological action of the toxin. These include the extent of the beta-sheets, the conformation of the region spanning residues 42-49 and the orientation of some side chains. In comparison with the X-ray structure, the NMR structure shows that the hydrophobic side chains of Trp(27) and Trp(36) are stacked together and are orientated towards the tip of the middle loop. The NMR study also showed that the two-stranded beta-sheet incorporated in the first loop, as defined by residues 1-22, and the C-terminus from Asn(59), is probably flexible relative to the rest of the molecule. On the basis of the dispositions of the hydrophobic and hydrophilic side chains, the structure of bucandin is clearly different from those of cytotoxins.

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.

NMR structure of the channel-former zervamicin IIB in isotropic solvents

Spatial structure of the membrane channel-forming hexadecapeptide, zervamicin IIB, was studied by NMR spectroscopy in mixed solvents of different polarity ranging from CDCl3/CD3OH (9:1, v/v) to CD3OH/H2O (1:1, v/v). The results show that in all solvents used the peptide has a very similar structure that is a bent amphiphilic helix with a mean backbone root mean square deviation (rmsd) value of ca. 0.3 A. Side chains of Trp1, Ile2, Gln3, Ile5 and Thr6 are mobile. The results are discussed in relation to the validity of the obtained structure to serve as a building block of zervamicin IIB ion channels.