Conformational sampling of influenza fusion peptide in membrane bilayers as a function of termini and protonation states

CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed

Since its inception nearly a half century ago, CHARMM has been playing a central role in computational biochemistry and biophysics. Commensurate with the developments in experimental research and advances in computer hardware, the range of methods and applicability of CHARMM have also grown. This review summarizes major developments that occurred after 2009 when the last review of CHARMM was published. They include the following: new faster simulation engines, accessible user interfaces for convenient workflows, and a vast array of simulation and analysis methods that encompass quantum mechanical, atomistic, and coarse-grained levels, as well as extensive coverage of force fields. In addition to providing the current snapshot of the CHARMM development, this review may serve as a starting point for exploring relevant theories and computational methods for tackling contemporary and emerging problems in biomolecular systems. CHARMM is freely available for academic and nonprofit research at https://academiccharmm.org/program.

Computational methods to study enveloped viral entry

The protein-membrane interactions that mediate viral infection occur via loosely ordered, transient assemblies, creating challenges for high-resolution structure determination. Computational methods and in particular molecular dynamics simulation have thus become important adjuncts for integrating experimental data, developing mechanistic models, and suggesting testable hypotheses regarding viral function. However, the large molecular scales of virus-host interaction also create challenges for detailed molecular simulation. For this reason, continuum membrane models have played a large historical role, although they have become less favored for high-resolution models of protein assemblies and lipid organization. Here, we review recent progress in the field, with an emphasis on the insight that has been gained using a mixture of coarse-grained and atomic-resolution molecular dynamics simulations. Based on successes and challenges to date, we suggest a multiresolution strategy that should yield the best mixture of computational efficiency and physical fidelity. This strategy may facilitate further simulations of viral entry by a broader range of viruses, helping illuminate the diversity of viral entry strategies and the essential common elements that can be targeted for antiviral therapies.

Molecular mechanisms of the influenza fusion peptide: insights from experimental and simulation studies

A key step in infections by enveloped viruses, such as influenza, is the fusion between the viral envelope and the host cell membrane, which allows the virus to insert its genetic material into the host cell and replicate. The influenza virus fusion process is promoted by hemagglutinin (HA), a glycoprotein that contains three identical monomers composed of two polypeptide chains (HA1 and HA2). Early studies on this protein revealed that HA-mediated fusion involves the insertion of the HA2 N-terminal segment into the host membrane and that this segment, known as the fusion peptide, is a key player in the fusion process. This mini-review highlights the main findings that have been obtained by experimental and computational studies on the HA fusion peptide, which give us a glimpse of its mode of action.

Hydrogen-Deuterium Exchange Supports Independent Membrane-Interfacial Fusion Peptide and Transmembrane Domains in Subunit 2 of Influenza Virus Hemagglutinin Protein, a Structured and Aqueous-Protected Connection between the Fusion Peptide and Soluble Ectodomain, and the Importance of Membrane Apposition by the Trimer-of-Hairpins Structure

The influenza virus hemagglutinin (HA) protein has HA1 and HA2 subunits, which form an initial complex. HA1's bind host cell sialic acids, which triggers endocytosis, HA1/HA2 separation, and HA2-mediated fusion between virus and endosome membranes. We report hydrogen-deuterium exchange mass spectrometry (HDX-MS) on the HA2 subunit without HA1. HA2 contains the fusion peptide (FP), soluble ectodomain (SE), transmembrane domain (TM), and endodomain. FP is a monomer by itself, while SE is a trimer of hairpins that includes an interior bundle of residue 38-105 helices, turns, and residue 154-178 strands packed antiparallel to the bundle. FP and TM extend from the same side of the SE hairpin, and fusion models often depict a FP/TM complex with membrane traversal of both domains that is important for membrane pore expansion. The HDX-MS data of this study do not support this complex and instead support independent FP and TM with respective membrane-interfacial and traversal locations. The data also show a low level of aqueous exposure of the 22-38 segment, consistent with retention of the 23-35 antiparallel β sheet observed in the initial HA1/HA2 complex. We propose the β sheet as a semirigid connector between FP and SE that enables close membrane apposition prior to fusion. The I173E mutant exhibits greater exchange for residues 22-69 and 150-191, consistent with dissociation of SE C-terminal strands from interior N-helices. Similar trends are observed for the G1E mutant as well as less exchange for G1E FP. Fusion is highly impaired with either mutant, which correlates with reduced membrane apposition and, for G1E, FP binding to SE rather than the target membrane.

Charged N-terminus of Influenza Fusion Peptide Facilitates Membrane Fusion

Cleavage of hemagglutinin precursor (HA0) by cellular proteases results in the formation of two subunits, HA1 and HA2. The N-terminal fragment of HA2, named a fusion peptide (HAfp), possess a charged, amine N-terminus. It has been shown that the N-terminus of HAfp stabilizes the structure of a helical hairpin observed for a 23-amino acid long peptide (HAfp1-23), whose larger activity than HAfp1-20 has been demonstrated recently. In this paper, we analyze the effect of N-terminal charge on peptide-mediated fusion efficiency and conformation changes at the membrane interface by comparison with the corresponding N-acetylated peptides of 20- and 23-amino acid lengths. We found that higher fusogenic activities of peptides with unmodified amino termini correlates with their ability to form helical hairpin structures oriented perpendicularly to the membrane plane. Molecular dynamics simulations showed that acetylated peptides adopt open and surface-bound conformation more often, which induced less disorder of the phospholipid chains, as compared to species with unmodified amino termini.

Influence of membrane composition on the binding and folding of a membrane lytic peptide from the non-enveloped flock house virus

Using a combination of coarse-grained and atomistic molecular dynamics simulations we have investigated the membrane binding and folding properties of the membrane lytic peptide of Flock House virus (FHV). FHV is an animal virus and an excellent model system for studying cell entry mechanisms in non-enveloped viruses. FHV undergoes a maturation event where the 44 C-terminal amino acids are cleaved from the major capsid protein, forming the membrane lytic (γ) peptides. Under acidic conditions, γ is released from the capsid interior allowing the peptides to bind and disrupt membranes. The first 21 N-terminal residues of γ, termed γ₁, have been resolved in the FHV capsid structure and γ₁ has been the subject of in vitro studies. γ₁ is structurally dynamic as it adopts helical secondary structure inside the capsid and on membranes, but it is disordered in solution. In vitro studies have shown the binding free energies to POPC or POPG membranes are nearly equivalent, but binding to POPC is enthalpically driven, while POPG binding is entropically driven. Through coarse-grained and multiple microsecond all-atom simulations the membrane binding and folding properties of γ₁ are investigated against homogeneous and heterogeneous bilayers to elucidate the dependence of the microenvironment on the structural properties of γ₁. Our studies provide a rationale for the thermodynamic data and suggest binding of γ₁ to POPG bilayers occurs in a disordered state, but γ₁ must adopt a helical conformation when binding POPC bilayers.

Heterogeneous dielectric generalized Born model with a van der Waals term provides improved association energetics of membrane-embedded transmembrane helices

The heterogeneous dielectric generalized Born (HDGB) implicit membrane formalism is extended by the addition of a van der Waals dispersion term to better describe the nonpolar components of the free energy of solvation. The new model, termed HDGBvdW, improves the energy estimates in the hydrophobic interior of the membrane, where polar and charged species are rarely found and nonpolar interactions become significant. The implicit van der Waals term for the membrane environment extends the model from Gallicchio et al. (J. Comput. Chem. 2004, 25, 479) by combining separate contributions from each of the membrane components. The HDGBvdW model is validated with a series of test cases ranging from membrane insertion and pair association profiles of amino acid side chain analogs and transmembrane helices. Overall, the HDGBvdW model leads to increased agreement with explicit membrane simulation results and experimental data. © 2016 Wiley Periodicals, Inc.

Three conserved C-terminal residues of influenza fusion peptide alter its behavior at the membrane interface

The N-terminal fragment of the viral hemagglutinin HA2 subunit is termed a fusion peptide (HAfp). The 23-amino acid peptide (HAfp1-23) contains three C-terminal W21-Y22-G23 residues which are highly conserved among serotypes of influenza A and has been shown to form a tight helical hairpin very distinct from the boomerang structure of HAfp1-20. We studied the effect of peptide length on fusion properties, structural dynamics, and binding to the membrane interface. We developed a novel fusion visualization assay based on FLIM microscopy on giant unilamellar vesicles (GUV). By means of molecular dynamics simulations and spectroscopic measurements, we show that the presence of the three C-terminal W21-Y22-G23 residues promotes the hairpin formation, which orients perpendicularly to the membrane plane and induces more disorder in the surrounding lipids than the less structured HAfp1-20. Moreover, we report cholesterol-enriched domain formation induced exclusively by the longer fusion peptide.

Assembly of Influenza Hemagglutinin Fusion Peptides in a Phospholipid Bilayer by Coarse-grained Computer Simulations

Membrane fusion is critical to eukaryotic cellular function and crucial to the entry of enveloped viruses such as influenza and human immunodeficiency virus. Influenza viral entry in the host cell is mediated by a 20–23 amino acid long sequence, called the fusion peptide (FP). Recently, possible structures for the fusion peptide (ranging from an inverted V shaped α-helical structure to an α-helical hairpin, or to a complete α-helix) and their implication in the membrane fusion initiation have been proposed. Despite the large number of studies devoted to the structure of the FP, the mechanism of action of this peptide remains unclear with several mechanisms having been suggested, including the induction of local disorder, promoting membrane curvature, and/or altering local membrane composition. In recent years, several research groups have employed atomistic and/or coarse-grained molecular dynamics (MD) simulations to investigate the matter. In all previous works, the behavior of a single FP monomer was studied, while in this manuscript, we use a simplified model of a tripeptide (TP) monomer of FP (TFP) instead of a single FP monomer because each Influenza Hemagglutinin contains three FP molecules in the biological system. In this manuscript we report findings targeted at understanding the fusogenic properties and the collective behavior of these trimers of FP peptides on a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine model membrane. Here we show how the TFP monomers self-assemble into differently sized oligomers in the presence of the membrane. We measure the perturbation to the structure of the phospholipid membrane caused by the presence of these TFP oligomers. Our work (i) shows how self-assembly of TFP in the presence of the membrane induces non negligible deformation to the membrane and (ii) could be a useful starting point to stimulate discussion and further work targeted to fusion pore formation.

Closed and Semiclosed Interhelical Structures in Membrane vs Closed and Open Structures in Detergent for the Influenza Virus Hemagglutinin Fusion Peptide and Correlation of Hydrophobic Surface Area with Fusion Catalysis

The ∼25 N-terminal "HAfp" residues of the HA2 subunit of the influenza virus hemagglutinin protein are critical for fusion between the viral and endosomal membranes at low pH. Earlier studies of HAfp in detergent support (1) N-helix/turn/C-helix structure at pH 5 with open interhelical geometry and N-helix/turn/C-coil structure at pH 7; or (2) N-helix/turn/C-helix at both pHs with closed interhelical geometry. These different structures led to very different models of HAfp membrane location and different models of catalysis of membrane fusion by HAfp. In this study, the interhelical geometry of membrane-associated HAfp is probed by solid-state NMR. The data are well-fitted to a population mixture of closed and semiclosed structures. The two structures have similar interhelical geometries and are planar with hydrophobic and hydrophilic faces. The different structures of HAfp in detergent vs membrane could be due to the differences in interaction with the curved micelle vs flat membrane with better geometric matching between the closed and semiclosed structures and the membrane. The higher fusogenicity of longer sequences and low pH is correlated with hydrophobic surface area and consequent increased membrane perturbation.

Capturing Spontaneous Membrane Insertion of the Influenza Virus Hemagglutinin Fusion Peptide

Hemagglutinin (HA) is a protein located on the surface of the influenza virus that mediates viral fusion to the host cellular membrane. During the fusion process the HA fusion peptide (HAfp), formed by the first 23 N-terminal residues of HA and structurally characterized by two alpha helices (Helix A and Helix B) tightly packed in a hairpin-like arrangement, is the only part of the virus in direct contact with the host membrane. After encountering the host cell, HAfp is believed to insert into the membrane, thereby initiating the fusion of the viral and host membranes. Detailed characterization of the interactions between the HAfp and cellular membrane is therefore of high relevance to the mechanism of viral entry into the host cell. Employing HMMM membrane representation with enhanced lipid mobility, we have performed a large set of independent simulations of unbiased membrane binding of HAfp. We have been able to capture spontaneous binding and insertion of HAfp consistently in nearly all the simulations. A reproducible membrane-bound configuration emerges from these simulations, despite employing a diverse set of initial configurations. Extension of several of the simulations into full membrane systems confirms the stability of the membrane-bound form obtained from HMMM binding simulations. The resulting model allows for the characterization of important interactions between the peptide and the membrane and the details of the binding process of the peptide for the first time. Upon membrane binding, Helix A inserts much deeper into the membrane than Helix B, suggesting that the former is responsible for hydrophobic anchoring of the peptide into the membrane. Helix B, in contrast, is found to establish major amphipathic interactions at the interfacial region thereby contributing to binding strength of HAfp.

Self-assembly molecular dynamics simulations shed light into the interaction of the influenza fusion Peptide with a membrane bilayer

Influenza virus is one of the most devastating human pathogens. In order to infect host cells, this virus fuses its membrane with the host membrane in a process mediated by the glycoprotein hemagglutinin. During fusion, the N-terminal region of hemagglutinin, which is known as the fusion peptide (FP), inserts into the host membrane, promoting lipid mixing between the viral and host membranes. Therefore, this peptide plays a key role in the fusion process, but the exact mechanism by which it promotes lipid mixing is still unclear. To shed light into this matter, we performed molecular dynamics (MD) simulations of the influenza FP in different environments (water, dodecylphosphocholine (DPC) micelles, and a dimyristoylphosphatidylcholine (DMPC) membrane). While in pure water the peptide lost its initial secondary structure, in simulations performed in the presence of DPC micelles it remained stable, in agreement with previous experimental observations. In simulations performed in the presence of a preassembled DMPC bilayer, the peptide became unstructured and was unable to insert into the membrane as a result of technical limitations of the method used. To overcome this problem, we used a self-assembly strategy, assembling the membrane together with the peptide. These simulations revealed that the peptide can adopt a membrane-spanning conformation, which had not been predicted by previous MD simulation studies. The peptide insertion had a strong effect on the membrane, lowering the bilayer thickness, disordering nearby lipids, and promoting lipid tail protrusion. These results contribute to a better understanding of the role of the FP in the fusion process.

Interactions of amino acid side-chain analogs within membrane environments

The interactions among four amino acid analog pairs (Asn, Ser, Phe, and Val) within the membrane environment were investigated using umbrella sampling molecular dynamics simulations. The results confirm generally expected qualitative trends of preferential association of polar compounds inside the membrane vs preferential interaction of hydrophobic compounds outside the membrane. Furthermore, correlations between amino acid interactions, membrane insertion, and membrane deformations are discussed and a detailed analysis of pair interaction energies is presented. A comparison of the energetics obtained from explicit lipid simulations with those from implicit membrane models reveals significant deviations and an improved parametrization of the heterogeneous dielectric generalized Born implicit model is provided that partially corrects for deficiencies in the implicit membrane model when compared with the new reference data from this study.

Transferring the PRIMO Coarse-Grained Force Field to the Membrane Environment: Simulations of Membrane Proteins and Helix-Helix Association

An extension of the recently developed PRIMO coarse-grained force field to membrane environments, PRIMO-M, is described. The membrane environment is modeled with the heterogeneous dielectric generalized Born (HDGB) methodology that simply replaces the standard generalized Born model in PRIMO without further parametrization. The resulting model was validated by comparing amino acid insertion free energy profiles and application in molecular dynamics simulations of membrane proteins and membrane-interacting peptides. Membrane proteins with 148–661 amino acids show stable root-mean-squared-deviations (RMSD) between 2 and 4 Å for most systems. Transmembrane helical peptides maintain helical shape and exhibit tilt angles in good agreement with experimental or other simulation data. The association of two glycophorin A (GpA) helices was simulated using replica exchange molecular dynamics simulations yielding the correct dimer structure with a crossing angle in agreement with previous studies. Finally, conformational sampling of the influenza fusion peptide also generates structures in agreement with previous studies. Overall, these findings suggest that PRIMO-M can be used to study membrane bound peptides and proteins and validates the transferable nature of the PRIMO coarse-grained force field.

Structure and dynamics of a fusion peptide helical hairpin on the membrane surface: comparison of molecular simulations and NMR

The conserved N-terminal residues of the HA2 subunit of influenza hemagglutinin (fusion peptide) are essential for membrane fusion and viral entry. Recent NMR studies showed that the 23-residue fusion peptide forms a helical hairpin that undergoes rocking motion relative to the membrane surface on a nanosecond time scale. To compare with NMR and to obtain a detailed molecular picture of the peptide-membrane interaction, we performed molecular dynamics simulations of the fusion peptide in explicit dimyristoylphosphatidylcholine and with the IMM1 implicit membrane model. To account for low and neutral pH conditions, simulations were performed with acidic groups (E11 and D19) protonated and unprotonated, respectively. The hairpin structure was stable in the simulations, with the N-terminal helix buried more deeply into the hydrophobic membrane interior than the C-terminal helix. Interactions between the tryptophans in the fusion peptide and phospholipid residues contribute to peptide orientation. Higher flexibility of the hairpin was observed in the implicit membrane simulations. Internal correlation functions of backbone N-H vectors were fit to the extended Lipari-Szabo model-free approach to obtain order parameters and correlation times. Good agreement with the NMR results was obtained for orientational fluctuations around the hairpin axis (rotation), but those around the perpendicular axis (tilting) were more limited in the simulations than inferred from the NMR experiments.

Plasticity and conformational equilibria of influenza fusion peptides in model lipid bilayers

Membrane fusion is critical to eukaryotic cellular function and crucial to the entry of enveloped viruses such as influenza and human immunodeficiency virus. Influenza viral entry in the host cell is mediated by a 20-23 amino acid long sequence, called the fusion peptide. In the last years, possible structures for the fusion peptide and their implication in the membrane fusion initiation have been proposed; these ranging from an inverted V shaped α-helical structure to an α-helical hairpin, or to a complete α-helix. Here we develop a coarse grained approach to describe effectively the plasticity of the fusion peptide and the explored conformational states. We describe also a trimeric assembly for the fusion peptide and analyse the explored states in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine model membrane. For the single fusion peptide systems the kink angle observed experimentally for the V shaped structure shows a strong correlation with the orientation of the fusion peptide within the lipid bilayer. The trimeric fusion peptide model also experiences different conformational states and represents a more realistic model for the anchoring mechanism of one influenza haemagglutinin molecule. This article is part of a Special Issue entitled: Viral Membrane Proteins - Channels for Cellular Networking.

Dynamic Heterogeneous Dielectric Generalized Born (DHDGB): An implicit membrane model with a dynamically varying bilayer thickness

An extension to the heterogeneous dielectric generalized Born (HDGB) implicit membrane formalism is presented to allow dynamic membrane deformations in response to membrane-inserted biomolecules during molecular dynamic simulations. The flexible membrane is implemented through additional degrees of freedom that represent the membrane deformation at the contact points of a membrane-inserted solute with the membrane. The extra degrees of freedom determine the dielectric and non-polar solvation free energy profiles that are used to obtain the solvation free energy in the presence of the membrane and are used to calculate membrane deformation free energies according to an elastic membrane model. With the dynamic HDGB (DHDGB) model the membrane is able to deform in response to the insertion of charged molecules thereby avoiding the overestimation of insertion free energies with static membrane models. The DHDGB model also allows the membrane to respond to the insertion of membrane-spanning solutes with hydrophobic mismatch. The model is tested with the membrane insertion of amino acid side chain analogs, arginine-containing helices, the WALP23 peptide, and the gramicidin A channel.

Detection of closed influenza virus hemagglutinin fusion peptide structures in membranes by backbone (13)CO- (15)N rotational-echo double-resonance solid-state NMR

The influenza virus fusion peptide is the N-terminal ~20 residues of the HA2 subunit of the hemagglutinin protein and this peptide plays a key role in the fusion of the viral and endosomal membranes during initial infection of a cell. The fusion peptide adopts N-helix/turn/C-helix structure in both detergent and membranes with reports of both open and closed interhelical topologies. In the present study, backbone (13)CO-(15)N REDOR solid-state NMR was applied to the membrane-associated fusion peptide to detect the distribution of interhelical distances. The data clearly showed a large fraction of closed and semi-closed topologies and were best-fitted to a mixture of two structures that do not exchange. One of the earlier open structural models may have incorrect G13 dihedral angles derived from TALOS analysis of experimentally correct (13)C shifts.

Role of conformational sampling of Ser16 and Thr17-phosphorylated phospholamban in interactions with SERCA

Phosphorylation of phospholamban (PLB) at Ser16 and/ or Thr17 is believed to release its inhibitory effect on sarcoplasmic reticulum calcium ATPase. Ser16 phosphorylation of PLB has been suggested to cause a conformational change that alters the interaction between the enzyme and protein. Using computer simulations, the conformational sampling of Ser16 phosphorylated PLB in implicit membrane environment is compared here with the unphosphorylated PLB system to investigate these conformational changes. The results suggest that conformational changes in the cytoplasmic domain of PLB upon phosphorylation at Ser16 increase the likelihood of unfavorable interactions with SERCA in the E2 state prompting a conformational switch of SERCA from E2 to E1. Phosphorylation of PLB at Thr17 on the other hand does not appear to affect interactions with SERCA significantly suggesting that the mechanism of releasing the inhibitory effect is different between Thr17 phosphorylated and Ser16 phosphorylated PLB.

The influenza fusion peptide adopts a flexible flat V conformation in membranes

Knowledge about the influenza fusion peptide (FP) membrane insertion mode is crucial for understanding its fusogenic mechanism. NMR and electron paramagnetic resonance experiments showed that in micelles, the FP inserted as a fixed-angle inverted V. In membranes, however, it was shown to insert as a straight α-helix (by molecular-dynamics simulations) and to adopt multiple kinked conformations (by solid-state NMR). In this work we performed explicit-solvent molecular-dynamics simulations of the influenza FP, and its F9A and W14A mutants, in POPC membranes. The Hα1 chemical shifts predicted from the molecular-dynamics structures are in excellent agreement with the experimental values obtained for the three peptides. The peptide orientation and conformations observed from the simulations lead to a flexible flat-V model in which the peptide lies almost flat on the membrane surface and alternates between kinked and straight-helix conformations.

Helical hairpin structure of influenza hemagglutinin fusion peptide stabilized by charge-dipole interactions between the N-terminal amino group and the second helix

The fusion domain of the influenza coat protein hemagglutinin HA2, bound to dodecyl phosphocholine micelles, was recently shown to adopt a structure consisting of two antiparallel α-helices, packed in an exceptionally tight hairpin configuration. Four interhelical H(α) to C═O aliphatic H-bonds were identified as factors stabilizing this fold. Here, we report evidence for an additional stabilizing force: a strong charge-dipole interaction between the N-terminal Gly(1) amino group and the dipole moment of helix 2. pH titration of the amino-terminal (15)N resonance, using a methylene-TROSY-based 3D NMR experiment, and observation of Gly(1 13)C' show a strongly elevated pK = 8.8, considerably higher than expected for an N-terminal amino group in a lipophilic environment. Chemical shifts of three C-terminal carbonyl carbons of helix 2 titrate with the protonation state of Gly(1)-N, indicative of a close proximity between the N-terminal amino group and the axis of helix 2, providing an optimal charge-dipole stabilization of the antiparallel hairpin fold. pK values of the side-chain carboxylate groups of Glu(11) and Asp(19) are higher by about 1 and 0.5 unit, respectively, than commonly seen for solvent-exposed side chains in water-soluble proteins, indicative of dielectric constants of ε = ∼30 (Glu(11)) and ∼60 (Asp(19)), placing these groups in the headgroup region of the phospholipid micelle.

Effect of flanking residues on the conformational sampling of the internal fusion peptide from Ebola virus

Fusion peptides mediate viral and host-cell membrane fusion during viral entry. The monomeric form of the internal fusion peptide from Ebola virus was studied in membrane bilayer and water environments with computer simulations using replica exchange sampling and an implicit solvent description of the environment. Wild-type Ebola fusion peptide (EFP), the W8A mutant form, and an extended construct with flanking residues were examined. It was found that the monomeric form of wild-type EFP adopts coil-helix-coil structure with a short helix from residues 8 to 11 mostly sampling orientations parallel to the membrane surface. W8A mutation disrupts the helicity in the N-terminal region of the peptide and leads to a preference for slightly oblique orientation relative to the membrane surface. The addition of flanking residues also alters the fusion peptide conformation with either a helix-break-helix structure or extended N and C-termini and reduced membrane insertion. In water, the fusion peptide is found to adopt structures with low helicity.

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.