The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipidprotein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Infection of human cells by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV2) relies on its binding to a specific receptor and subsequent fusion of the viral and host cell membranes. The fusion peptide (FP), a short peptide segment in the spike protein, plays a central role in the initial penetration of the virus into the host cell membrane, followed by the fusion of the two membranes. Here, we use an array of molecular dynamics simulations that take advantage of the highly mobile membrane mimetic model to investigate the interaction of the SARS-CoV2 FP with a lipid bilayer representing mammalian cellular membranes at an atomic level and to characterize the membrane-bound form of the peptide. Six independent systems were generated by changing the initial positioning and orientation of the FP with respect to the membrane, and each system was simulated in five independent replicas, each for 300 ns. In 73% of the simulations, the FP reaches a stable, membrane-bound configuration, in which the peptide deeply penetrated into the membrane. Clustering of the results reveals three major membrane-binding modes (binding modes 1-3), in which binding mode 1 populates over half of the data points. Taking into account the sequence conservation among the viral FPs and the results of mutagenesis studies establishing the role of specific residues in the helical portion of the FP in membrane association, the significant depth of penetration of the whole peptide, and the dense population of the respective cluster, we propose that the most deeply inserted membrane-bound form (binding mode 1) represents more closely the biologically relevant form. Analysis of FP-lipid interactions shows the involvement of specific residues, previously described as the ''fusion-active core residues,'' in membrane binding. Taken together, the results shed light on a key step involved in SARS-CoV2 infection, with potential implications in designing novel inhibitors.

Eukaryotic cell homeostasis requires transfer of cellular components among organelles and relies on membrane fusion catalyzed by SNARE proteins. Inactive SNARE bundles are reactivated by hexameric N-ethylmaleimide-sensitive factor, vesicle-fusing ATPase (Sec18/NSF)-driven disassembly that enables a new round of membrane fusion. We previously found that phosphatidic acid (PA) binds Sec18 and thereby sequesters it from SNAREs and that PA dephosphorylation dissociates Sec18 from the membrane, allowing it to engage SNARE complexes. We now report that PA also induces conformational changes in Sec18 protomers and that hexameric Sec18 cannot bind PA membranes. Molecular dynamics (MD) analyses revealed that the D1 and D2 domains of Sec18 contain PA-binding sites and that the residues needed for PA binding are masked in hexameric Sec18. Importantly, these simulations also disclosed that a major conformational change occurs in the linker region between the D1 and D2 domains, which is distinct from the conformational changes that occur in hexameric Sec18 during SNARE priming. Together, these findings indicate that PA regulates Sec18 function by altering its architecture and stabilizing membrane-bound Sec18 protomers.Membrane fusion is necessary for all eukaryotes to effectively transport cellular components between organelles. The trafficking and fusion of vesicles is carried out through a series of events that are highly conserved across eukarya (1). Although many proteins that drive the process may differ between eukaryotic species, they all perform similar roles allowing compartment contact, bilayer fusion, and luminal content mixing (2). The final stage of membrane fusion, and luminal content mixing, is catalyzed by SNARE 3 proteins. Each participating membrane contributes either an R-SNARE or three Q-SNARE coils that wrap around each other to form a parallel four-helical trans-SNARE complex that brings membranes into close apposition. The formation of such complexes releases free energy that is transmitted to the membranes to trigger fusion. Once fusion occurs and membranes are merged, the four-helical SNARE bundle, now a cis-SNARE complex, is inactive and requires disassembly to undergo a new round of fusion.The disassembly of cis-SNAREs, also known as Priming, is carried out by the AAA ϩ protein Sec18/NSF and its adaptor protein Sec17/␣-SNAP (3) (Fig. 1A). Current models suggest that NSF primes cis-SNAREs through a "loaded spring" mechanism triggered by cis-SNARE recognition and ATP hydrolysis (4). NSF binds to cis-SNAREs with the help of ␣-SNAP to form what is known as the 20S complex (5-8). Although NSF was originally isolated as a trimer or tetramer, it can only prime SNAREs as a homohexamer that surrounds the cis-SNAREs and ␣-SNAP proteins to form the 20S particle (9 -11). Association with cis-SNARE-␣-SNAP complexes triggers ATP hydrolysis, which leads to a large conformational change in the protein, with the major change occurring at the N terminus where it folds back over the D1-D2 rings (8). This generates eno...

No abstract

Infection of human cells by the SARS-CoV2 relies on its binding to a specific receptor and subsequent fusion of the viral and host cell membranes. The fusion peptide (FP), a short peptide segment in the spike protein, plays a central role in the initial penetration of the virus into the host cell membrane, followed by the fusion of the two membranes. Here, we use an array of molecular dynamics (MD) simulations taking advantage of the Highly Mobile Membrane Mimetic (HMMM) model, to investigate the interaction of the SARS-CoV2 FP with a lipid bilayer representing mammalian cellular membranes at an atomic level, and to characterize the membrane-bound form of the peptide. Six independent systems were generated by changing the initial positioning and orientation of the FP with respect to the membrane, and each system was simulated in five independent replicas. In 60% of the simulations, the FP reaches a stable, membrane-bound configuration where the peptide deeply penetrated into the membrane. Clustering of the results reveals two major membrane binding modes, the helix-binding mode and the loop-binding mode. Taken into account the sequence conservation among the viral FPs and the results of mutagenesis studies establishing the role of specific residues in the helical portion of the FP in membrane association, we propose that the helix-binding mode represents more closely the biologically relevant form. In the helix-binding mode, the helix is stabilized in an oblique angle with respect to the membrane with its N-terminus tilted towards the membrane core. Analysis of the FP-lipid interactions shows the involvement of specific residues of the helix in membrane binding previously described as the fusion active core residues. Taken together, the results shed light on a key step involved in SARS-CoV2 infection with potential implications in designing novel inhibitors.