Interaction of cardiotoxins with membranes: a molecular modeling study

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.

The secretory phenotypes of envenomed cells: Insights into venom cytotoxicity

Snake envenomation is listed as Category A Neglected Tropical Diseases (NTD) by World Health Organization, indicates a severe public health problem. The global figures for envenomation cases are estimated to be more than 1.8 million annually. Even if the affected victims survive the envenomation, they might suffer from permanent morbidity due to local envenomation. One of the most prominent local envenomation is dermonecrosis. Dermonecrosis is a pathophysiological outcome of envenomation that often causes disability in the victims due to surgical amputations, deformities, contracture, and chronic ulceration. The key venom toxins associated with this local symptom are mainly attributed to substantial levels of enzymatic and non-enzymatic toxins as well as their possible synergistic actions. Despite so, the severity of the local tissue damage is based on macroscopic observation of the bite areas. Furthermore, limited knowledge is known about the key biomarkers involved in the pathogenesis of dermonecrosis. The current immunotherapy with antivenom is also ineffective against dermonecrosis. These local effects eventually end up as sequelae. There is also a global shortage of toxins-targeted therapeutics attributed to inadequate knowledge of the actual molecular mechanisms of cytotoxicity. This chapter discusses the characterization of secretory phenotypes of dermonecrosis as an advanced tool to indicate its severity and pathogenesis in envenomation. Altogether, the secretory phenotypes of envenomed cells and tissues represent the precise characteristics of dermonecrosis caused by venom toxins.

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.

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.

Blocking of negative charged carboxyl groups converts Naja atra neurotoxin to cardiotoxin-like protein

Naja atra cobrotoxin and cardiotoxin 3 (CTX3) exhibit neurotoxicity and cytotoxicity, respectively. In the present study, we aimed to investigate whether the carboxyl groups of cobrotoxin play a role in structural constraints, thereby preventing cobrotoxin from exhibiting cytotoxic activity. Six of the seven carboxyl groups in cobrotoxin were conjugated with semicarbazide. Measurement of circular dichroism spectra and Trp fluorescence quenching showed that the gross conformation of semicarbazide-modified cobrotoxin (SEM-cobrotoxin) and cobrotoxin differed. In sharp contrast to cobrotoxin, SEM-cobrotoxin demonstrated membrane-damaging activity and cytotoxicity, which are feature more characteristic of CTX3. Furthermore, both SEM-cobrotoxin and CTX3 induced cell death through AMPK activation. Analyses of the interaction between polydiacetylene/lipid vesicles and fluorescence-labeled lipids revealed that SEM-cobrotoxin and cobrotoxin adopted different membrane-bound states. The structural characteristics of SEM-cobrotoxin were similar to those of CTX3, including trifluoroethanol (TFE)-induced structural transformation and membrane binding-induced conformational change. Conversely, cobrotoxin was insensitive to the TFE-induced effect. Collectively, the data of this study indicate that blocking negatively charged residues confers cobrotoxin with membrane-damaging activity and cytotoxicity. The findings also suggest that the structural constraints imposed by carboxyl groups control the functional properties of snake venom α-neurotoxins during the divergent evolution of snake venom neurotoxins and cardiotoxins.

Status of Asp29 and Asp40 in the Interaction of Naja atra Cardiotoxins with Lipid Bilayers

It is widely accepted that snake venom cardiotoxins (CTXs) target the plasma membranes of cells. In the present study, we investigated the role of Asp residues in the interaction of Naja atra cardiotoxin 1 (CTX1) and cardiotoxin 3 (CTX3) with phospholipid bilayers using chemical modification. CTX1 contains three Asp residues at positions 29, 40, and 57; CTX3 contains two Asp residues at positions 40 and 57. Compared to Asp29 and Asp40, Asp57 was sparingly modified with semi-carbazide, as revealed by matrix-assisted laser desorption ionization-time of flight (MALDI-TOF) mass and mass/mass analyses. Thus, semi-carbazide-modified CTX1 (SEM-CTX1) mainly contained modified Asp29 and Asp40, while SEM-CTX3 contained modified Asp40. Compared to that of native toxins, trifluoroethanol easily induced structural transition of SEM-CTX1 and SEM-CTX3, suggesting that the structural flexibility of CTXs was constrained by Asp40. Modification of Asp29 and Asp40 markedly promoted the ability of CTX1 to induce permeability of cell membranes and lipid vesicles; CTX3 and SEM-CTX3 showed similar membrane-damaging activity. Modification of Asp residues did not affect the membrane-binding capability of CTXs. Circular dichroism spectra of SEM-CTX3 and CTX3 were similar, while the gross conformation of SEM-CTX1 was distinct from that of CTX1. The interaction of CTX1 with membrane was distinctly changed by Asp modification. Collectively, our data suggest that Asp29 of CTX1 suppresses the optimization of membrane-bound conformation to a fully active state and that the function of Asp40 in the structural constraints of CTX1 and CTX3 is not important for the manifestation of membrane-perturbing activity.

Monitoring the Disulfide Bonds of Folding Isomers of Synthetic CTX A3 Polypeptide Using MS-Based Technology

Native disulfide formation is crucial to the process of disulfide-rich protein folding in vitro. As such, analysis of the disulfide bonds can be used to track the process of the folding reaction; however, the diverse structural isomers interfere with characterization due to the non-native disulfide linkages. Previously, a mass spectrometry (MS) based platform coupled with peptide demethylation and an automatic disulfide bond searching engine demonstrated the potential to screen disulfide-linked peptides for the unambiguous assignment of paired cysteine residues of toxin components in cobra venom. The developed MS-based platform was evaluated to analyze the disulfide bonds of structural isomers during the folding reaction of synthetic cardiotoxin A3 polypeptide (syn-CTX A3), an important medical component in cobra venom. Through application of this work flow, a total of 13 disulfide-linked peptides were repeatedly identified across the folding reaction, and two of them were found to contain cysteine pairings, like those found in native CTX A3. Quantitative analysis of these disulfide-linked peptides showed the occurrence of a progressive disulfide rearrangement that generates a native disulfide bond pattern on syn-CTX A3 folded protein. The formation of these syn-CTX A3 folded protein reaches a steady level in the late stage of the folding reaction. Biophysical and cell-based assays showed that the collected syn-CTX A3 folded protein have a β-sheet secondary structure and cytotoxic activity similar to that of native CTX A3. In addition, the immunization of the syn-CTX A3 folded proteins could induce neutralization antibodies against the cytotoxic activity of native CTX A3. In contrast, these structure activities were poorly observed in the other folded isomers with non-native disulfide bonds. The study highlights the ability of the developed MS platform to assay isomers with heterogeneous disulfide bonds, providing insight into the folding mechanism of the bioactive protein generation.

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.

Three-Finger Proteins from the Ly6/uPAR Family: Functional Diversity within One Structural Motif

The discovery in higher animals of proteins from the Ly6/uPAR family, which have structural homology with snake "three-finger" neurotoxins, has generated great interest in these molecules and their role in the functioning of the organism. These proteins have been found in the nervous, immune, endocrine, and reproductive systems of mammals. There are two types of the Ly6/uPAR proteins: those associated with the cell membrane by GPI-anchor and secreted ones. For some of them (Lynx1, SLURP-1, SLURP-2, Lypd6), as well as for snake α-neurotoxins, the target of action is nicotinic acetylcholine receptors, which are widely represented in the central and peripheral nervous systems, and in many other tissues, including epithelial cells and the immune system. However, the targets of most proteins from the Ly6/uPAR family and the mechanism of their action remain unknown. This review presents data on the structural and functional properties of the Ly6/uPAR proteins, which reveal a variety of functions within a single structural motif.

Novel Naja atra cardiotoxin 1 (CTX-1) derived antimicrobial peptides with broad spectrum activity

Naja atra subsp. atra cardiotoxin 1 (CTX-1), produced by Chinese cobra snakes, belonging to Elapidae family, is included in the three-finger toxin family and exerts high cytotoxicity and antimicrobial activity too. Using as template mainly the tip and the subsequent β-strand of the first "finger" of this toxin, different sequences of 20 amino acids linear peptides have been designed in order to avoid toxic effects but to maintain or even strengthen the partial antimicrobial activity already seen for the complete toxin. As a result, the sequence NCP-0 (Naja Cardiotoxin Peptide-0) was designed as ancestor and subsequently 4 other variant sequences of NCP-0 were developed. These synthesized variant sequences have shown microbicidal activity towards a panel of reference and field strains of Gram-positive and Gram-negative bacteria. The sequence named NCP-3, and its variants NCP-3a and NCP-3b, have shown the best antimicrobial activity, together with low cytotoxicity against eukaryotic cells and low hemolytic activity. Bactericidal activity has been demonstrated by minimum bactericidal concentration (MBC) assay at values below 10 μg/ml for most of the tested bacterial strains. This potent antimicrobial activity was confirmed even for unicellular fungi Candida albicans, Candida glabrata and Malassezia pachydermatis (MBC 50-6.3 μg/ml), and against the fast-growing mycobacteria Mycobacterium smegmatis and Mycobacterium fortuitum. Moreover, NCP-3 has shown virucidal activity on Bovine Herpesvirus 1 (BoHV1) belonging to Herpesviridae family. The bactericidal activity is maintained even in a high salt concentration medium (125 and 250 mM NaCl) and phosphate buffer with 20% Mueller Hinton (MH) medium against E. coli, methicillin resistant Staphylococcus aureus (MRSA) and Pseudomonas aeruginosa reference strains. Considering these in vitro obtained data, the search for active sequences within proteins presenting an intrinsic microbicidal activity could provide a new way for discovering a large number of novel and promising antimicrobial peptides families.

Towards universal approach for bacterial production of three-finger Ly6/uPAR proteins: Case study of cytotoxin I from cobra N. oxiana

Cytotoxins or cardiotoxins is a group of polycationic toxins from cobra venom belonging to the 'three-finger' protein superfamily (Ly6/uPAR family) which includes small β-structural proteins (60-90 residues) with high disulfide bond content (4-5 disulfides). Due to a high cytotoxic activity for cancer cells, cytotoxins are considered as potential anticancer agents. Development of the high-throughput production methods is required for the prospective applications of cytotoxins. Here, efficient approach for bacterial production of recombinant analogue of cytotoxin I from N. oxiana containing additional N-terminal Met-residue (rCTX1) was developed. rCTX1 was produced in the form of E. coli inclusion bodies. Refolding in optimized conditions provided ∼6 mg of correctly folded protein from 1 L of bacterial culture. Cytotoxicity of rCTX1 for C6 rat glioma cells was found to be similar to the activity of wild type CTX1. The milligram quantities of ¹³C,¹⁵N-labeled rCTX1 were obtained. NMR study confirmed the similarity of the spatial structures of recombinant and wild-type toxins. Additional Met residue does not perturb the overall structure of the three-finger core. The analysis of available data for different Ly6/uPAR proteins of snake and human origin revealed that efficiency of their folding in vitro is correlated with the number of proline residues in the third loop and the surface area of hydrophobic residues buried within the protein interior. The obtained data indicate that hydrophobic core is important for the folding of proteins with high disulfide bond content. Developed expression method opens new possibilities for structure-function studies of CTX1 and other related three-finger proteins.

Effect of externally applied jidesheng anti-venom on skin and soft-tissue necrosis after Chinese cobra bite: a retrospective study

OBJECTIVE

To evaluate the effects of Jidesheng anti-venom used externally for skin and soft-tissue necrosis from Chinese cobra bite.

METHODS

A retrospective review was performed according to the clinical data recorded from January 2002 to December 2012. A total of 126 patients (116 females and 10 males) with skin and soft-tissue necrosis due to Chinese cobra bite were divided into two groups. The control group was treated externally with 40% glyceride magnesium sulfate (n = 52), and the treatment group was given Jidesheng anti-venom externally (n = 74). The data collected included maximum local necrotic area of skin and soft tissues, de-tumescence time, healing time, and skin-grafting rate.

RESULTS

There were no significant differences in gender, age, and disease condition between the control and treatment groups (P > 0.05). No statistically significant difference was found in de-tumescence time between the two groups (P > 0.05). The maximum local necrotic area of skin and soft tissues was (19.9 +/- 7.3) cm2 in the treatment group, while it was (23.3 +/- 6.4) cm2 in the control group. The healing time of the treatment group was shorter than that of the control group [(32.1 +/- 3.7) vs (34.4 +/- 4.5) days)]. The skin-grafting rate in the treatment group was lower than that of the control group (10.81% vs 25.00%). There were statistically significant differences in maximum local necrotic area of skin and soft tissues, healing time, and skin-grafting rate between the control and treatment groups (all P < 0.05).

CONCLUSION

External application of Jidesheng anti-venom may help to promote wound healing and reduce the skin-grafting rate in cases of skin and soft-tissue necrosis due to Chinese cobra bite.

Novel lynx spider toxin shares common molecular architecture with defense peptides from frog skin

A unique 30-residue cationic peptide oxyopinin 4a (Oxt 4a) was identified in the venom of the lynx spider Oxyopes takobius (Oxyopidae). Oxt 4a contains a single N-terminally located disulfide bond, Cys4-Cys10, and is structurally different from any spider toxin studied so far. According to NMR findings, the peptide is disordered in water, but assumes a peculiar torpedo-like structure in detergent micelles. It features a C-terminal amphipathic α-helical segment (body; residues 12-25) and an N-terminal disulfide-stabilized loop (head; residues 1-11), and has an unusually high density of positive charge in the head region. Synthetic Oxt 4a was produced and shown to possess strong and broad-spectrum cytolytic and antimicrobial activity. cDNA cloning showed that the peptide is synthesized in the form of a conventional prepropeptide with an acidic prosequence. Unlike other arachnid toxins, Oxt 4a exhibits striking similarity with defense peptides from the skin of ranid frogs that contain the so-called Rana-box motif (a C-terminal disulfide-enclosed loop). Parallelism or convergence is apparent on several levels: the structure, function and biosynthesis of a lynx spider toxin are mirrored by those of Rana-box peptides from frogs.

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.

Coarse-grained molecular dynamics simulations of cobra cytotoxin A3 interactions with a lipid bilayer: penetration of loops into membranes

Cobra cytotoxins, which are small three-looped proteins composed of approximately 60 amino acid residues, primarily act by destroying the bilayer membranes of cells and artificial vesicles. However, the molecular mechanism governing this process is not yet completely understood. We used coarse-grained molecular dynamics (CGMD) simulations to study the mechanism underlying the penetration of cardiotoxin A3 (CTX A3), the major toxic component of Naja atra (Chinese cobra) venom, into a hydrated 1-palmitoyl-2-oleoyl-1-sn-3-phosphatidylcholine (POPC) lipid bilayer. We performed CGMD simulations for three different conformations of the cobra cytotoxin-the tail, lying, and harrow conformations. The results of our simulations indicate that two of these, the tail and lying conformations, did not penetrate the bilayer system. Further, for the harrow conformation, loops 2 and 3 played important roles in penetration of CTX A3 into the bilayer system.

Anionic lipids: determinants of binding cytotoxins from snake venom on the surface of cell membranes

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.

Role of heparan sulfates and glycosphingolipids in the pore formation of basic polypeptides of cobra cardiotoxin

Cobra venom contains cardiotoxins (CTXs) that induce tissue necrosis and systolic heart arrest in bitten victims. CTX-induced membrane pore formation is one of the major mechanisms responsible for the venom's designated cytotoxicity. This chapter examines how glycoconjugates such as heparan sulfates (HS) and glycosphingolipids, located respectively in the extracellular matrix and lipid bilayers of the cell membranes, facilitate CTX pore formation. Evidences for HS-facilitated cell surface retention and glycosphingolipid-facilitated membrane bilayer insertion of CTX are reviewed. We suggest that similar physical steps could play a role in the mediation of other pore forming toxins (PFT). The membrane pores formed by PFT are expected to have limited lifetime on biological cell surface as a result of membrane dynamics during endocytosis and/or rearrangement of lipid rafts.

Steered molecular dynamics simulations of cobra cytotoxin interaction with zwitterionic lipid bilayer: no penetration of loop tips into membranes

Cobra cytotoxins, small proteins of three-fingered toxin family, unspecifically damage membranes in different cells and artificial vesicles. However, the molecular mechanism of this damage is not yet completely understood. We used steered molecular dynamics simulations to study the interaction of cardiotoxin A3 from Naja atra cobra venom with hydrated 1-palmitoyl-2-oleoyl-1-sn-3-phosphatidylcholine (POPC) bilayer. The studied system included one cytotoxin molecule, 64 lipid molecules (32 molecules in each monolayer) and 2500 water molecules. It was found that the toxin interacted with zwitterionic bilayer formed by POPC. During first nanosecond of simulation the toxin molecule was oriented toward membrane surface by loops' basement including cytotoxin regions Cys14-Asn19 and Cys38-Ser46. This orientation was stable enough and was not changed during next 6 ns of simulation. The obtained data suggest that cytotoxin molecule cannot penetrate into membrane composed of zwitterionic lipids without some auxiliary interaction.

Monte Carlo simulations of tBid association with the mitochondrial outer membrane

Bid, a BH3-only pro-apoptopic member of the BCL-2 protein family, regulates cell death at the level of mitochondrial cytochrome c efflux. Bid consists of 8 alpha-helices (H1-H8, respectively) and is soluble cytosolic protein in its native state. Proteolysis of the N-terminus (encompassing H1 and H2) of Bid by caspase 8 in apoptosis yields activated "tBid" (truncated Bid), which translocates to the mitochondria and induces the efflux of cytochrome c. The release of cytochrome c from mitochondria to the cytosol constitutes a critical control point in apoptosis that is regulated by interaction of tBid protein with mitochondrial membrane. tBid displays structural homology to channel-forming bacterial toxins, such as colicins or transmembrane domain of diphtheria toxin. By analogy, it has been hypothesized that tBid would unfold and insert into the lipid bilayer of the mitochondria outer membrane (MOM) upon membrane association. However, it has been shown recently that unlike colicins and the transmembrane domain of diphtheria toxin, tBid binds to the lipid bilayer maintaining alpha-helical conformation of its helices without adopting a transmembrane orientation by them. Here, the mechanism of the association of tBid with the model membrane mimicking the mitochondrial membrane is studied by Monte Carlo simulations, taking into account the underlying energetics. A novel two-stage hierarchical simulation protocol combining coarse-grained discretization of conformational space with subsequent refinements was applied which was able to generate the protein conformation and its location in the membrane using modest computational resources. The simulations show that starting from NMR-established conformation in the solution, the protein associates with the membrane without adopting the transmembrane orientation. The configuration (conformation and location) of tBid providing the lowest free energy for the system protein/membrane/solvent has been obtained. The simulations reveal that tBid upon association with the membrane undergoes significant conformational changes primarily due to rotations within the loops between helices H4 and H5, H6 and H7, H7 and H8. It is established that in the membrane-bound state of tBid-monomer helices H3 and H5 have the locations exposed to the solution, helices H6 and H8 are partly buried and helices H4 and H7 are buried into the membrane at shallow depth. The average orientation of tBid bound to the membrane in the most stable configuration reported here is in satisfactory agreement with the evaluations obtained by indirect experimental means.

Spatial Structure and Activity Mechanism of a Novel Spider Antimicrobial Peptide,

Latarcins (Ltc), linear peptides (ca. 25 amino acid long) isolated from the venom of the Lachesana tarabaevi spider, exhibit a broad-spectrum antibacterial activity, most likely acting on the bacterial plasmatic membrane. We study the structure-activity relationships in the series of these compounds. At the first stage, we investigated the spatial structure of one of the peptides, Ltc2a, and its mode of membrane perturbation. This was done by a combination of experimental and theoretical methods. The approach includes (i) structural study of the peptide by CD spectroscopy in phospholipid liposomes and by (1)H NMR in detergent micelles, (ii) determination of the effect on the liposomes by a dye leakage fluorescent assay and (31)P NMR spectroscopy, (iii) refinement of the NMR-derived spatial structure via Monte Carlo simulations in an implicit water-octanol slab, and (iv) calculation of the molecular hydrophobicity potential. The molecule of Ltc2a was found to consist of two helical regions (residues 3-9 and 13-21) connected via a poorly ordered fragment. The effect of the peptide on the liposomes suggests the carpet mechanism of the membrane deterioration. This is also supported by the analysis of hydrophobic/hydrophilic characteristics of Ltc2a and homologous antimicrobial peptides. These peptides exhibiting a helix-hinge-helix structural motif are characterized by a distinct and feebly marked amphiphilicity of their N- and C-terminal helices, respectively, and by a hydrophobicity gradient along the peptide chain. The approach we suggested may be useful in studying not only other latarcins but also a wider class of membrane-active peptides.

Implicit solvent simulations of peptide interactions with anionic lipid membranes

A recently developed implicit membrane model (IMM1) is supplemented with a Gouy-Chapman term describing counterion-screened electrostatic interactions of a solute with negatively charged membrane lipids. The new model is tested on peptides that bind to anionic membranes. Pentalysine binds just outside the plane of negative charge, whereas Lys-Phe peptides insert their aromatic rings into the hydrophobic core. Melittin and magainin 2 bind more strongly to anionic than to neutral membranes and in both cases insert their hydrophobic residues into the hydrocarbon core. The third domain of Antennapedia homeodomain (penetratin) binds as an alpha-helix in the headgroup region. Cardiotoxin II binds strongly to anionic membranes but marginally to neutral ones. In all cases, the location and configuration of the peptides are consistent with experimental data, and the effective energy changes upon binding compare favorably with experimental binding free energies. The model opens the way to exploring the effect of membrane charge on the location, conformation, and dynamics of a large variety of biologically active peptides on membranes.

Association of transmembrane helices: what determines assembling of a dimer?

Self-association of two hydrophobic alpha-helices is studied via unrestrained Monte Carlo (MC) simulations in a hydrophobic slab described by an effective potential. The system under study represents two transmembrane (TM) segments of human glycophorin A (GpA), which form homo-dimers in membranes. The influence of TM electrostatic potential, thickness and hydrophobicity degree of lipid bilayer is investigated. It is shown that the membrane environment stabilizes alpha-helical conformation of GpA monomers, induces their TM insertion and facilitates inter-helical contacts. Head-to-head orientation of the helices is promoted by the voltage difference across the membrane. Subsequent "fine-tuned" assembling of the dimer is mediated by van der Waals interactions. Only the models of dimer, calculated in a hydrophobic slab with applied voltage agree with experimental data, while simulations in vacuo or without TM voltage fail to give reasonable results. The moderate structural heterogeneity of GpA dimers (existence of several groups of states with close energies) is proposed to reflect their equilibrium dynamics in membrane-mimics. The calculations performed for GpA mutants G83A and G86L permit rationalization of mutagenesis data for them. The results of Monte Carlo simulations critically depend on the parameters of the membrane model: adequate description of helix association is achieved in the water-cyclohexane-water system with the membrane thickness 30-34 A, while in membranes with different hydrophobicities and thickness unrealistic conformations of the dimer are found. The computational approach permits efficient prediction of TM helical oligomers based solely on the sequences of interacting peptides.

Prediction of amphipathic in-plane membrane anchors in monotopic proteins using a SVM classifier

Background

Membrane proteins are estimated to represent about 25% of open reading frames in fully sequenced genomes. However, the experimental study of proteins remains difficult. Considerable efforts have thus been made to develop prediction methods. Most of these were conceived to detect transmembrane helices in polytopic proteins. Alternatively, a membrane protein can be monotopic and anchored via an amphipathic helix inserted in a parallel way to the membrane interface, so-called in-plane membrane (IPM) anchors. This type of membrane anchor is still poorly understood and no suitable prediction method is currently available.

Results

We report here the "AmphipaSeeK" method developed to predict IPM anchors. It uses a set of 21 reported examples of IPM anchored proteins. The method is based on a pattern recognition Support Vector Machine with a dedicated kernel.

Conclusion

AmphipaSeeK was shown to be highly specific, in contrast with classically used methods (e.g. hydrophobic moment). Additionally, it has been able to retrieve IPM anchors in naively tested sets of transmembrane proteins (e.g. PagP). AmphipaSeek and the list of the 21 IPM anchored proteins is available on NPS@, our protein sequence analysis server.

Modelling of proteins in membranes

This review describes some recent theories and simulations of mesoscopic and microscopic models of lipid membranes with embedded or attached proteins. We summarize results supporting our understanding of phenomena for which the activities of proteins in membranes are expected to be significantly affected by the lipid environment. Theoretical predictions are pointed out, and compared to experimental findings, if available. Among others, the following phenomena are discussed: interactions of interfacially adsorbed peptides, pore-forming amphipathic peptides, adsorption of charged proteins onto oppositely charged lipid membranes, lipid-induced tilting of proteins embedded in lipid bilayers, protein-induced bilayer deformations, protein insertion and assembly, and lipid-controlled functioning of membrane proteins.

Molecular dynamics simulation of the enhancement of cobra cardiotoxin and E6 protein binding on mixed self-assembled monolayer molecules

Molecular dynamics simulations are performed on n-alkinethiol self-assembled monolayers (SAMs) and their mixture on a gold surface so that the orientations of the binding of cobra cardiotoxin and E6 protein molecules can be selected using the mixing ratio of CH3-terminated SAMs with different chain lengths. The simulations suggest that a SAM surface with different mixing ratios may provide a possible platform for aligning protein molecules with a desired orientation and for enhancing the binding energy of the protein on the designed surface.

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.

Interactions of hydrophobic peptides with lipid bilayers: Monte Carlo simulations with M2delta

We introduce here a novel Monte Carlo simulation method for studying the interactions of hydrophobic peptides with lipid membranes. Each of the peptide's amino acids is represented as two interaction sites: one corresponding to the backbone alpha-carbon and the other to the side chain, with the membrane represented as a hydrophobic profile. Peptide conformations and locations in the membrane and changes in the membrane width are sampled using the Metropolis criterion, taking into account the underlying energetics. Using this method we investigate the interactions between the hydrophobic peptide M2delta and a model membrane. The simulations show that starting from an extended conformation in the aqueous phase, the peptide first adsorbs onto the membrane surface, while acquiring an ordered helical structure. This is followed by formation of a helical-hairpin and insertion into the membrane. The observed path is in agreement with contemporary understanding of peptide insertion into biological membranes. Two stable orientations of membrane-associated M2delta were obtained: transmembrane (TM) and surface, and the value of the water-to-membrane transfer free energy of each of them is in agreement with calculations and measurements on similar cases. M2delta is most stable in the TM orientation, where it assumes a helical conformation with a tilt of 14 degrees between the helix principal axis and the membrane normal. The peptide conformation agrees well with the experimental data; average root-mean-square deviations of 2.1 A compared to nuclear magnetic resonance structures obtained in detergent micelles and supported lipid bilayers. The average orientation of the peptide in the membrane in the most stable configurations reported here, and in particular the value of the tilt angle, are in excellent agreement with the ones calculated using the continuum-solvent model and the ones observed in the nuclear magnetic resonance studies. This suggests that the method may be used to predict the three-dimensional structure of TM peptides.

Effective energy function for proteins in lipid membranes

A simple extension of the EEF1 energy function to heterogeneous membrane-aqueous media is proposed. The extension consists of (a) development of solvation parameters for a nonpolar phase using experimental data for the transfer of amino acid side-chains from water to cyclohexane, (b) introduction of a heterogeneous membrane-aqueous system by making the reference solvation free energy of each atom dependent on the vertical coordinate, (c) a modification of the distance-dependent dielectric model to account for reduced screening of electrostatic interactions in the membrane, and (d) an adjustment of the EEF1 aqueous model in light of recent calculations of the potential of mean force between amino acid side-chains in water. The electrostatic model is adjusted to match experimental observations for polyalanine, polyleucine, and the glycophorin A dimer. The resulting energy function (IMM1) reproduces the preference of Trp and Tyr for the membrane interface, gives reasonable energies of insertion into or adsorption onto a membrane, and allows stable 1-ns MD simulations of the glycophorin A dimer. We find that the lowest-energy orientation of melittin in bilayers varies, depending on the thickness of the hydrocarbon layer.

Peripheral binding mode and penetration depth of cobra cardiotoxin on phospholipid membranes as studied by a combined FTIR and computer simulation approach

Cobra cardiotoxin, a cytotoxic beta-sheet basic polypeptide, is known to cause membrane leakage in many cells including human erythrocytes. Herein, we demonstrate that the major cobra cardiotoxin from Naja atra, CTX A3, can cause leakage of vesicle contents in phosphatidylglycerol (PG) and phosphatidylserine containing, but not in pure phosphatidylcholine (PC), membrane bilayers. By the combined polarized attenuated total reflection infrared spectroscopy and computer simulation studies, CTX A3 is shown to peripherally bind to both zwitterionic and anionic monolayers in a similar edgewise manner with a tilted angle of approximately 48 +/- 20 degrees between the beta-sheet plane of the CTX molecule and the normal of the membrane surface. The average surface area expansion induced by CTX A3 binding to the PG monolayer, however, is two times larger than that of the PC monolayer as determined by the Langmuir minitrough method. Interaction energy considerations of CTX A3 on neutral and negatively charged membrane surfaces suggests that the electrostatic interaction between anionic lipid and cationic CTXs plays a role in modulating the penetration depth of CTX molecules on the initial peripheral binding mode and reveals a pathway leading to the formation of an inserted mode in negatively charged membrane bilayers.

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.

Monte Carlo simulations of voltage-driven translocation of a signal sequence

Transmembrane potentials play important but poorly understood roles in many biological processes, including signal sequence-mediated protein translocation across bacterial membranes. In this study we applied Monte Carlo techniques to simulate the way the potential acts on a signal sequence. The simulations demonstrate that in the absence of a potential the signal sequence prefers insertion in both helical hairpin and transmembrane alpha-helical conformations. However, in the presence of a potential only the transmembrane alpha-helical conformation is the state of lowest energy for the signal sequence. From these results it is concluded that the membrane potential stabilizes the transmembrane orientation of a signal sequence, explaining the membrane potential dependence of preprotein translocation.