Natural polycationic membrane-active peptides typically lack disulfide bonds and exhibit fusion, cell-penetrating, antimicrobial activities. They are mostly unordered in solution, but adopt a helical structure, when bound to phospholipid membranes. Structurally different are cardiotoxins (or cytotoxins, CTs) from cobra venom. They are fully β- structured molecules, characterized by the three-finger fold (TFF). Affinity of CTs to lipid bilayer was shown to depend on amino acid sequence in the tips of the three loops. In the present review, CT-membrane interactions are analyzed through the prism of data on binding of the toxins to phospholipid liposomes and detergent micelles, as well as their structural and computational studies in membrane mimicking environments. We assess different hydrophobicity scales to compare membrane partitioning of various CTs and their membrane effects. A comparison of hydrophobic/hydrophilic properties of CTs and linear polycationic peptides provides a key to their biological activity and creates a fundamental basis for rational design of new membrane-interacting compounds, including new promising drugs. For instance, from the viewpoint of the data obtained on model lipid membranes, cytotoxic activity of CTs against cancer cells is discussed.

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.

Cytotoxins (or cardiotoxins; CTs) are toxins from cobra venom characterized by the three-finger (TF) fold. CTs are on average 60-residue-long peptides, possessing as many as 4 disulfide bonds. The elements of antiparallel β-structure take origin from the hydrophobic core formed by the disulfides. The β-strands adopt the shape of the three loops, giving the name of the fold. While neurotoxins (NTs) - also TF proteins from snake venom - exert their effect through specific interactions with protein receptors, no specific protein target has been found for CTs. Unlike NTs, CTs are amphiphilic and cytotoxic against a variety of cells, including cancer ones. Thus, the hypothesis that the activity of CTs is caused by their interactions with lipid membranes is currently central. To understand molecular basis behind variations in toxicities of CTs highly homologous in their sequences, detailed knowledge of their structure and dynamics is required. The present review summarizes experimental and computational data on the spatial organization of CTs and their dynamics in various environments (aqueous solution, membranous milieus).

Membrane and membrane-active peptides and proteins play a crucial role in numerous cell processes, such as signaling, ion conductance, fusion, and others. Many of them act as highly specific and efficient drugs or drug targets, and, therefore, attract growing interest of medicinal chemists. Because of experimental difficulties with characterization of their spatial structure and mode of membrane binding, essential attention is given now to molecular modeling techniques. During the last years an important progress has been achieved in molecular dynamics (MD) and Monte Carlo (MC) simulations of peptides and proteins with explicit and/or implicit theoretical models of membranes. The first ones allow atomic-resolution studies of peptides behavior on the membrane-water interfaces. Models with implicit consideration of membrane are of a special interest because of their computational efficiency and ability to account for principal trends in protein-lipid interactions. In this approximation, the bilayer is usually treated as continuum whose properties vary along the membrane thickness, and membrane insertion is simulated using either MC or MD methods. This review surveys recent applications of both types of lipid bilayer models in computer simulations of a wide variety of peptides and proteins with different biological activities. Theoretical background of the membrane models is considered with examples of their applications to biologically relevant problems. The emphasis of the review is made on recent MC and MD computations, on structural and/or functional information, which may be obtained via molecular modeling. The approximations and shortcomings of the models, along with their perspectives in design of new membrane active drugs, are discussed.

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.

The cytotoxic properties of cytotoxins (CTs) from snake venom are mediated by their interaction with the cell membrane. The hydrophobic pattern containing the tips of loops I-III and flanked by polar residues is known to be a membrane-binding motif of CTs. However, this is not enough to explain the difference in activity among various CTs which are similar in sequence and in 3D structure. The mechanism of further CT-membrane interaction leading to pore formation and cell death still remains unknown. Published experimental data on the specific interaction between CT and low molecular weight anionic components (sulphatide) of the bilayer point to the existence of corresponding ligand binding sites on the surface of toxin molecules. In this work we study the membrane-lytic properties of CT I, CT II (Naja oxiana), and CT 4 (Naja kaouthia), which belong to different structural and functional types (P-and S-type) of CTs, by measuring the intensity of a fluorescent dye, calcein released from liposomes containing a phosphatidylserine (PS) lipid as an anionic component. Using molecular docking simulations, we find and characterize three sites in CT molecules that can potentially bind the PS polar head. Based on the data obtained, we suggest a hypothesis that CTs can specifically interact with one or more of the anionic lipids (in particular, with PS) contained in the membrane, thus facilitating the interaction between CTs and the lipid bilayer of a cell membrane.

Cobra cytotoxins (CTs) belong to the three-fingered protein family. They are classified into S- and P-types, the latter exhibiting higher membrane-perturbing capacity. In this work, we investigated the interaction of CTs with phospholipid bilayers, using coarse-grained (CG) and full-atom (FA) molecular dynamics (MD). The object of this work is a CT of an S-type, cytotoxin I (CT1) from N.oxiana venom. Its spatial structure in aqueous solution and in the micelles of dodecylphosphocholine (DPC) were determined by H-NMR spectroscopy. Then, via CG- and FA MD-computations, we evaluated partitioning of CT1 molecule into palmitoyloleoylphosphatidylcholine (POPC) membrane, using the toxin spatial models, obtained either in aqueous solution, or detergent micelle. The latter model exhibits minimal structural changes upon partitioning into the membrane, while the former deviates from the starting conformation, loosing the tightly bound water molecule in the loop-2. These data show that the structural changes elicited by CT1 molecule upon incorporation into DPC micelle take place likely in the lipid membrane, although the mode of the interaction of this toxin with DPC micelle (with the tips of the all three loops) is different from its mode in POPC membrane (primarily with the tip of the loop-1 and both the tips of the loop-1 and loop-2).

Cardiotoxins (CTs) from snake venoms are a family of homologous highly basic proteins that have extended hydrophobic patterns on their molecular surfaces. CTs are folded into three β-structured loops stabilized by four disulfide bridges. Being well-structured in aqueous solution, most of these proteins are membrane-active, although the exact molecular mechanisms of CT-induced cell damage are still poorly understood. To elucidate the structure-function relationships in CTs, a detailed knowledge of their spatial organization and local conformational dynamics is required. Protein domain motions can be either derived from a set of experimental structures or generated via molecular dynamics (MD). At the same time, traditional clustering algorithms in the Cartesian coordinate space often fail to properly take into account the local large-scale dihedral angle transitions that occur in MD simulations. This is because such perturbations are usually offset by changes in the neighboring dihedrals, thus preserving the overall protein fold. States with a "locally perturbed" backbone were found in experimental 3D models of some globular proteins and have been shown to be functionally meaningful. In this work, the possibility of large-scale dihedral angle transitions in the course of long-term MD in explicit water was explored for three CTs with different membrane activities: CT 1, 2 (Naja oxiana) and CT A3 (Naja atra). Analysis of the MD-derived distributions of backbone torsion angles revealed several important common and specific features in the structural/dynamic behavior of these proteins. First, large-amplitude transitions were detected in some residues located in the functionally important loop I region. The K5/L6 pair of residues was found to induce a perturbation of the hydrophobic patterns on the molecular surface of CTs-reversible breaking of a large nonpolar zone ("bottom") into two smaller ones and their subsequent association. Second, the characteristic sizes of these patterns perfectly coincided with the dimensions of the nonpolar zones on the surfaces of model two-component (zwitterionic/anionic) membranes. Taken together with experimental data on the CT-induced leakage of fluorescent dye from such membranes, these results allowed us to formulate a two-stage mechanism of CT-membrane binding. The principal finding of this study is that even local conformational dynamics of CTs can seriously affect their functional activity via a tuning of the membrane binding site - specific "hot spots" (like the K5/L6 pair) in the protein structure.