Putative membrane lytic sites of P-type and S-type cardiotoxins from snake venoms as probed by all-atom molecular dynamics simulations.

Cardiotoxins (CTXs) belonging to the three-finger toxin superfamily of snake venoms are one of principal toxic components and the protein toxins exhibit membrane lytic activities when the venoms are injected into victims. In the present study, complex formations between CTX VI (a P-type CTX from Naja atra) and CTX1 (an S-type CTX from Naja naja) on zwitterionic POPC bilayers (a major lipid component of cell membranes) have been studied in near physiological conditions for a total dynamic time scale of 1.35 μs using all-atom molecular dynamics (MD) simulations. Comprehensive analyses of the MD data revealed that residues such as Leu1, Lys2, Tyr11, Lys31, Asp57 and Arg58 of CTX VI, and Ala16, Lys30 and Arg58 of CTX1 were crucial for establishing interactions with the POPC bilayer. Moreover, loop I, along with globular head and loop II of CTX VI, and loop II of CTX1 were found to be the structural regions chiefly governing complex formation of the respective proteins with POPC. Rationalizations for the differential binding modes of CTXs and implications of the findings for designing small molecular inhibitors to the toxins are also discussed. Graphical Abstract Binding modes of a P-type CTX and an S-type CTX towards the POPC bilayer.