Planar aggregation of the influenza viral fusion peptide alters membrane structure and hydration, promoting poration

Positive residues of the SARS-CoV-2 fusion domain are key contributors to the initiation of membrane fusion

SARS-CoV-2 is one of the most infectious viruses ever recorded. Despite a plethora of research over the last several years, the viral lifecycle is still not well understood, particularly membrane fusion. This process is initiated by the fusion domain (FD), a highly conserved stretch of amino acids consisting of a fusion peptide (FP) and fusion loop (FL), which in synergy perturbs the target cells lipid membrane to lower the energetic cost necessary for fusion. In this study, through a mutagenesis-based approach, we have investigated the basic residues within the FD (K825, K835, R847, K854) utilizing an in vitro fusion assay and 19F NMR, validated by traditional 13C 15N techniques. Alanine and charge conserving mutants revealed every basic residue plays a highly specific role within the mechanism of initiating fusion. Intriguingly, K825A led to increased fusogenecity which was found to be correlated to the number of amino acids within helix one, further implicating the role of this specific helix within the FD's fusion mechanism. This work has found basic residues to be important within the FDs fusion mechanism and highlights K825A, a specific mutation made within the FD of the SARS-CoV-2 spike protein, as requiring further investigation due to its potential to contribute to a more virulent strain of SARS-CoV-2.

Making the cut: Multiscale simulation of membrane remodeling

Biological membranes are dynamic heterogeneous materials, and their shape and organization are tightly coupled to the properties of the proteins in and around them. However, the length scales of lipid and protein dynamics are far below the size of membrane-bound organelles, much less an entire cell. Therefore, multiscale modeling approaches are often necessary to build a comprehensive picture of the interplay of these factors, and have provided critical insights into our understanding of membrane dynamics. Here, we review computational methods for studying membrane remodeling, as well as passive and active examples of protein-driven membrane remodeling. As the field advances towards the modeling of key aspects of organelles and whole cells - an increasingly accessible regime of study - we summarize here recent successes and offer comments on future trends.

Simulating asymmetric membranes using P21 periodic boundary conditions

Molecular dynamics (MD) simulations of symmetric lipid bilayers are now well established, while those of asymmetric ones are considerably less developed. This disjunction arises in part because the surface tensions of leaflets in asymmetric bilayers can differ (unlike those of symmetric ones), and there is no simple way to determine them without assumptions. This chapter describes the use of P21 periodic boundary conditions (PBC), which allow lipids to switch leaflets, to generate asymmetric bilayers under the assumption of equal chemical potentials of lipids in opposing leaflets. A series of examples, ranging from bilayers with one lipid type to those with peptides and proteins, provides a guide for the use of P21 PBC. Critical properties of asymmetric membranes, such as spontaneous curvature, are highly sensitive to differences in the leaflet surface tensions (or differential stress), and equilibration with P21 PBC substantially reduces differential stress of asymmetric bilayers assembled with surface area-based methods. Limitations of the method are discussed. Technically, the nonstandard unit cell is difficult to parallelize and to incorporate restraints. Inherently, the assumption of equal chemical potentials, and therefore the method itself, is not applicable to all target systems. Despite these limitations, it is argued that P21 simulations should be considered when designing equilibration protocols for MD studies of most asymmetric membranes.

Spontaneous curvature generation by peptides in asymmetric bilayers

Peptides and proteins play crucial roles in membrane remodeling by inducing spontaneous curvature. However, extracting spontaneous curvatures from simulations of asymmetric bilayers is challenging because differential stress (i.e., the difference of the leaflet surface tensions) arising from leaflet area strains can vary substantially among initial conditions. This study investigates peptide-induced spontaneous curvatureδc 0 p in asymmetric bilayers consisting of a single lipid type and a peptide confined to one leaflet;δc 0 p is calculated from the Helfrich equation using the first moment of the lateral pressure tensor and an alternative expression using the differential stress. It is shown that differential stress introduced during initial system generation is effectively relaxed by equilibrating using P21 periodic boundary conditions, which allows lipids to switch leaflets across cell boundaries and equalize their chemical potentials across leaflets. This procedure leads to robust estimates ofδc 0 p for the systems simulated, and is recommended when equality of chemical potentials between the leaflets is a primary consideration.

Modeling asymmetric cell membranes at all-atom resolution

Molecular dynamics (MD) simulations are a useful tool when studying the properties of membranes as they allow for a molecular view of lipid interactions with proteins, nucleic acids, or small molecules. While model membranes are usually symmetric in their lipid composition between leaflets and include a small number of lipid components, physiological membranes are highly complex and vary in the level of asymmetry. Simulation studies have shown that changes in leaflet asymmetry can alter the properties of a membrane. It is therefore necessary to carefully build asymmetric membranes to accurately simulate membranes. This chapter carefully describes the different methods for building asymmetric membranes and the advantages/disadvantages of each method. The simplest methods involve building a membrane with either an equal number of lipids per leaflet or an equal initial surface area (SA) estimated by the area per lipid. More detailed methods include combining two symmetric membranes of equal SA or altering an asymmetric membrane and adjusting the number of lipids after equilibration to minimize an observable such as differential stress (0-DS). More complex methods that require specific simulation software are also briefly described. The challenges and assumptions are listed for each method which should help guide the researcher to choose the best method for their unique MD simulation of an asymmetric membrane.

Rigidifying of the internal dynamics of amyloid-beta fibrils generated in the presence of synaptic plasma vesicles

We investigated the changes in internal flexibility of amyloid-β1-40 (Aβ) fibrils grown in the presence of rat synaptic plasma vesicles. The fibrils are produced using a modified seeded growth protocol, in which the Aβ concentration is progressively increased at the expense of the decreased lipid to protein ratio. The morphologies of each generation are carefully assessed at several fibrils' growth time points using transmission electron microscopy. The side-chain dynamics in the fibrils is investigated using deuterium solid-state NMR measurements, with techniques spanning line shapes analysis and several NMR relaxation rates measurements. The dynamics is probed in the site-specific fashion in the hydrophobic C-terminal domain and the disordered N-terminal domain. An overall strong rigidifying effect is observed in comparison with the wild-type fibrils generated in the absence of the membranes. In particular, the overall large-scale fluctuations of the N-terminal domain are significantly reduced, and the activation energies of rotameric inter-conversion in methyl-bearing side-chains of the core (L17, L34, M35, V36), as well as the ring-flipping motions of F19 are increased, indicating a restricted core environment. Membrane-induced flexibility changes in Aβ aggregates can be important for the re-alignment of protein aggregates within the membrane, which in turn would act as a disruption pathway of the bilayers' integrity.

The underlying mechanical properties of membranes tune their ability to fuse

Membrane fusion is a ubiquitous process associated with a multitude of biological events. Although it has long been appreciated that membrane mechanics plays an important role in membrane fusion, the molecular interplay between mechanics and fusion has remained elusive. For example, although different lipids modulate membrane mechanics differently, depending on their composition, molar ratio, and complex interactions, differing lipid compositions may lead to similar mechanical properties. This raises the question of whether (i) the specific lipid composition or (ii) the average mesoscale mechanics of membranes acts as the determining factor for cellular function. Furthermore, little is known about the potential consequences of fusion on membrane disruption. Here, we use a combination of confocal microscopy, time-resolved imaging, and electroporation to shed light onto the underlying mechanical properties of membranes that regulate membrane fusion. Fusion efficiency follows a nearly universal behavior that depends on membrane fluidity parameters, such as membrane viscosity and bending rigidity, rather than on specific lipid composition. This helps explaining why the charged and fluid membranes of the inner leaflet of the plasma membrane are more fusogenic than their outer counterparts. Importantly, we show that physiological levels of cholesterol, a key component of biological membranes, has a mild effect on fusion but significantly enhances membrane mechanical stability against pore formation, suggesting that its high cellular levels buffer the membrane against disruption. The ability of membranes to efficiently fuse while preserving their integrity may have given evolutionary advantages to cells by enabling their function while preserving membrane stability.

SARS-CoV-2 fusion peptide sculpting of a membrane with insertion of charged and polar groups

The fusion peptide of SARS-CoV-2 spike is essential for infection. How this charged and hydrophobic domain occupies and affects membranes needs clarification. Its depth in zwitterionic, bilayered micelles at pH 5 (resembling late endosomes) was measured by paramagnetic NMR relaxation enhancements used to bias molecular dynamics simulations. Asp830 inserted deeply, along with Lys825 or Lys835. Protonation of Asp830 appeared to enhance agreement of simulated and NMR-measured depths. While the fusion peptide occupied a leaflet of the DMPC bilayer, the opposite leaflet invaginated with influx of water and choline head groups in around Asp830 and bilayer-inserted polar side chains. NMR-detected hydrogen exchange found corroborating hydration of the backbone of Thr827-Phe833 inserted deeply in bicelles. Pinching of the membrane at the inserted charge and the intramembrane hydration of polar groups agree with theory. Formation of corridors of hydrated, inward-turned head groups was accompanied by flip-flop of head groups. Potential roles of the defects are discussed.

Two modes of fusogenic action for influenza virus fusion peptide

The entry of influenza virus into the host cell requires fusion of its lipid envelope with the host membrane. It is catalysed by viral hemagglutinin protein, whose fragments called fusion peptides become inserted into the target bilayer and initiate its merging with the viral membrane. Isolated fusion peptides are already capable of inducing lipid mixing between liposomes. Years of studies indicate that upon membrane binding they form bend helical structure whose degree of opening fluctuates between tightly closed hairpin and an extended boomerang. The actual way in which they initiate fusion remains elusive. In this work we employ atomistic simulations of wild type and fusion inactive W14A mutant of influenza fusion peptides confined between two closely apposed lipid bilayers. We characterise peptide induced membrane perturbation and determine the potential of mean force for the formation of the first fusion intermediate, an interbilayer lipid bridge called stalk. Our results demonstrate two routes through which the peptides can lower free energy barrier towards fusion. The first one assumes peptides capability to adopt transmembrane configuration which subsequently promotes the creation of a stalk-hole complex. The second involves surface bound peptide configuration and proceeds owing to its ability to stabilise stalk by fitting into the region of extreme negative membrane curvature resulting from its formation. In both cases, the active peptide conformation corresponds to tight helical hairpin, whereas extended boomerang geometry appears to be unable to provide favourable thermodynamic effect. The latter observation offers plausible explanation for long known inactivity of boomerang-stabilising W14A mutation.

pH-dependent endocytosis mechanisms for influenza A and SARS-coronavirus

The ongoing SARS-CoV-2 pandemic and the influenza epidemics have revived the interest in understanding how these highly contagious enveloped viruses respond to alterations in the physicochemical properties of their microenvironment. By understanding the mechanisms and conditions by which viruses exploit the pH environment of the host cell during endocytosis, we can gain a better understanding of how they respond to pH-regulated anti-viral therapies but also pH-induced changes in extracellular environments. This review provides a detailed explanation of the pH-dependent viral structural changes preceding and initiating viral disassembly during endocytosis for influenza A (IAV) and SARS coronaviruses. Drawing upon extensive literature from the last few decades and latest research, I analyze and compare the circumstances in which IAV and SARS-coronavirus can undertake endocytotic pathways that are pH-dependent. While there are similarities in the pH-regulated patterns leading to fusion, the mechanisms and pH activation differ. In terms of fusion activity, the measured activation pH values for IAV, across all subtypes and species, vary between approximately 5.0 to 6.0, while SARS-coronavirus necessitates a lower pH of 6.0 or less. The main difference between the pH-dependent endocytic pathways is that the SARS-coronavirus, unlike IAV, require the presence of specific pH-sensitive enzymes (cathepsin L) during endosomal transport. Conversely, the conformational changes in the IAV virus under acidic conditions in endosomes occur due to the specific envelope glycoprotein residues and envelope protein ion channels (viroporins) getting protonated by H+ ions. Despite extensive research over several decades, comprehending the pH-triggered conformational alterations of viruses still poses a significant challenge. The precise mechanisms of protonation mechanisms of certain during endosomal transport for both viruses remain incompletely understood. In absence of evidence, further research is needed.

CHARMM-GUI Membrane Builder: Past, Current, and Future Developments and Applications

Molecular dynamics simulations of membranes and membrane proteins serve as computational microscopes, revealing coordinated events at the membrane interface. As G protein-coupled receptors, ion channels, transporters, and membrane-bound enzymes are important drug targets, understanding their drug binding and action mechanisms in a realistic membrane becomes critical. Advances in materials science and physical chemistry further demand an atomistic understanding of lipid domains and interactions between materials and membranes. Despite a wide range of membrane simulation studies, generating a complex membrane assembly remains challenging. Here, we review the capability of CHARMM-GUI Membrane Builder in the context of emerging research demands, as well as the application examples from the CHARMM-GUI user community, including membrane biophysics, membrane protein drug-binding and dynamics, protein-lipid interactions, and nano-bio interface. We also provide our perspective on future Membrane Builder development.

Initiation and evolution of pores formed by influenza fusion peptides probed by lysolipid inclusion

The fusion peptide (FP) domain is necessary for the fusogenic activity of spike proteins in a variety of enveloped viruses, allowing the virus to infect the host cell, and is the only part of the protein that interacts directly with the target membrane lipid tails during fusion. There are consistent findings of poration by this domain in experimental model membrane systems, and, in certain conditions, the isolated FPs can generate pores. Here, we use molecular dynamics simulations to investigate the specifics of how these FP-induced pores form in membranes with different compositions of lysolipid and POPC. The simulations show that pores form spontaneously at high lysolipid concentrations via hybrid intermediates, where FP aggregates in the cis leaflet tilt to form a funnel-like structure that spans the leaflet and locally reduces the hydrophobic thickness that must be traversed by water to form a pore. By restraining a single FP within an FP aggregate to this tilted conformation, pores can be formed in lower-lysolipid-content membranes, including pure POPC, on the 100-ns timescale, much more rapidly than in unbiased simulations in bilayers with the same composition. The pore formation pathway is similar to the spontaneous formation in high lysolipid concentrations. Depending on the membrane composition, the pores can be metastable (as seen in POPC) or lead to membrane rupture.