Measuring pKa of activation and pKi of inactivation for influenza hemagglutinin from kinetics of membrane fusion of virions and of HA expressing cells

Pleomorphism in Biological Units of Life: Morphological Heterogeneity in Cells Does Not Translate Uniformly to Subcellular Components

The interplay of the three-dimensional (3D) distribution of various subcellular components and their interactions are expected to control overall cellular morphology in biology. In this study, we aimed to determine whether the pleomorphy observed at the whole-cell level is being reflected by the components constituting the cells by focusing on the 3D distribution of pixel intensities at the single-cell level of the whole (cell) and its parts (the seven subcellular components of the cells-self-assemblies of smaller units). We rigorously acquired and analyzed the image data of RAW264.7 cells at the single-cell level. We report asymmetries in the spatial distribution of pixel intensities at the whole-cell and subcellular component levels along with the occurrence of alterations when pleomorphism is reduced by synchronization of the cell cycle. From our repertoire of seven subcellular components, we report ER, mitochondria, and tubulin to be independent of whole-cell apico-basal heterogeneity of optical density while nuclear, plasma membrane, lysosomal, and actin fluorescence distributions are found to contribute to the apico-basal polarity of the whole cell. While doing so, we have also developed an image analysis algorithm utilizing 2D segmentation to analyze the single cells in 3D using confocal microscopy, a technique that allows us to analyze cellular states in their native hydrated state.

(Re-)Viewing Role of Intracellular Glucose Beyond Extracellular Regulation of Glucose-Stimulated Insulin Secretion by Pancreatic Cells

For glucose-stimulated insulin secretion (GSIS) by pancreatic β-cells in animals, it is believed that ATP generated from glucose metabolism is primarily responsible. However, this ignores two well-established aspects in literature: (a) intracellular ATP generation from other sources resulting in an overall pool of ATP, regardless of the original source, and (b) that intracellular glucose transport is 10- to 100-fold higher than intracellular glucose phosphorylation in β-cells. The latter especially provides an earlier unaddressed, but highly appealing, observation pertaining to (at least transient) the presence of intracellular glucose molecules. Could these intracellular glucose molecules be responsible for the specificity of GSIS to glucose (instead of the widely believed ATP production from its metabolism)? In this work, we provide a comprehensive compilation of literature on glucose and GSIS using various cellular systems - all studies focus only on the extracellular role of glucose in GSIS. Further, we carried out a comprehensive analysis of differential gene expression in Mouse Insulinoma 6 (MIN6) cells, exposed to low and high extracellular glucose concentrations (EGC), from the existing whole transcriptome data. The expression of other genes involved in glycolysis, Krebs cycle, and electron transport chain was found to be unaffected by EGC, except Gapdh, Atp6v0a4, and Cox20. Remarkably, 3 upregulated genes (Atp6v0a4, Cacnb4, Kif11) in high EGC were identified to have an association with cellular secretion. Using glucose as a possible ligand for the 3 proteins, computational investigations were carried out (that will require future 'wet validation', both in vitro and in vivo, e.g., using primary islets and animal models). The glucose-affinity/binding scores (in kcal/mol) obtained were also compared with glucose binding scores for positive controls (GCK and GLUT2), along with negative controls (RPA1, KU70-80, POLA1, ACAA1A, POLR1A). The binding affinity scores of glucose molecules for the 3 proteins were found to be closer to positive controls. Therefore, we report the glucose binding ability of 3 secretion-related proteins and a possible direct role of intracellular glucose molecules in GSIS.

Aspects of Biological Replication and Evolution Independent of the Central Dogma: Insights from Protein-Free Vesicular Transformations and Protein-Mediated Membrane Remodeling

Biological membrane remodeling is central to living systems. In spite of serving as “containers” of whole-living systems and functioning as dynamic compartments within living systems, biological membranes still find a “blue collar” treatment compared to the “white collar” nucleic acids and proteins in biology. This may be attributable to the fact that scientific literature on biological membrane remodeling is only 50 years old compared to ~ 150 years of literature on proteins and a little less than 100 years on nucleic acids. However, recently, evidence for symbiotic origins of eukaryotic cells from data only on biological membranes was reported. This, coupled with appreciation of reproducible amphiphilic self-assemblies in aqueous environments (mimicking replication), has already initiated discussions on origins of life beyond nucleic acids and proteins. This work presents a comprehensive compilation and meta-analyses of data on self-assembly and vesicular transformations in biological membranes—starting from model membranes to establishment of Influenza Hemagglutinin-mediated membrane fusion as a prototypical remodeling system to a thorough comparison between enveloped mammalian viruses and cellular vesicles. We show that viral membrane fusion proteins, in addition to obeying “stoichiometry-driven protein folding”, have tighter compositional constraints on their amino acid occurrences than general-structured proteins, regardless of type/class. From the perspective of vesicular assemblies and biological membrane remodeling (with and without proteins) we find that cellular vesicles are quite different from viruses. Finally, we propose that in addition to pre-existing thermodynamic frameworks, kinetic considerations in de novo formation of metastable membrane structures with available “third-party” constituents (including proteins) were not only crucial for origins of life but also continue to offer morphological replication and/or functional mechanisms in modern life forms, independent of the central dogma.

New Biophysical Approaches Reveal the Dynamics and Mechanics of Type I Viral Fusion Machinery and Their Interplay with Membranes

Protein-mediated membrane fusion is a highly regulated biological process essential for cellular and organismal functions and infection by enveloped viruses. During viral entry the membrane fusion reaction is catalyzed by specialized protein machinery on the viral surface. These viral fusion proteins undergo a series of dramatic structural changes during membrane fusion where they engage, remodel, and ultimately fuse with the host membrane. The structural and dynamic nature of these conformational changes and their impact on the membranes have long-eluded characterization. Recent advances in structural and biophysical methodologies have enabled researchers to directly observe viral fusion proteins as they carry out their functions during membrane fusion. Here we review the structure and function of type I viral fusion proteins and mechanisms of protein-mediated membrane fusion. We highlight how recent technological advances and new biophysical approaches are providing unprecedented new insight into the membrane fusion reaction.

Hemagglutinin Cleavability, Acid Stability, and Temperature Dependence Optimize Influenza B Virus for Replication in Human Airways

Influenza A virus (IAV) and influenza B virus (IBV) cause yearly epidemics with significant morbidity and mortality. When zoonotic IAVs enter the human population, the viral hemagglutinin (HA) requires adaptation to achieve sustained virus transmission. In contrast, IBV has been circulating in humans, its only host, for a long period of time. Whether this entailed adaptation of IBV HA to the human airways is unknown. To address this question, we compared two seasonal IAVs (A/H1N1 and A/H3N2) and two IBVs (B/Victoria and B/Yamagata lineages) with regard to host-dependent activity of HA as the mediator of membrane fusion during viral entry. We first investigated proteolytic activation of HA by covering all type II transmembrane serine protease (TTSP) and kallikrein enzymes, many of which proved to be present in human respiratory epithelium. The IBV HA0 precursor is cleaved by a broader panel of TTSPs and activated with much higher efficiency than IAV HA0. Accordingly, knockdown of a single protease, TMPRSS2, abrogated spread of IAV but not IBV in human respiratory epithelial cells. Second, the HA fusion pH values proved similar for IBV and human-adapted IAVs (with one exception being the HA of 1918 IAV). Third, IBV HA exhibited higher expression at 33°C, a temperature required for membrane fusion by B/Victoria HA. This indicates pronounced adaptation of IBV HA to the mildly acidic pH and cooler temperature of human upper airways. These distinct and intrinsic features of IBV HA are compatible with extensive host adaptation during prolonged circulation of this respiratory virus in the human population.IMPORTANCE Influenza epidemics are caused by influenza A and influenza B viruses (IAV and IBV, respectively). IBV causes substantial disease; however, it is far less studied than IAV. While IAV originates from animal reservoirs, IBV circulates in humans only. Virus spread requires that the viral hemagglutinin (HA) is active and sufficiently stable in human airways. We resolve here how these mechanisms differ between IBV and IAV. Whereas human IAVs rely on one particular protease for HA activation, this is not the case for IBV. Superior activation of IBV by several proteases should enhance shedding of infectious particles. IBV HA exhibits acid stability and a preference for 33°C, indicating pronounced adaptation to the human upper airways, where the pH is mildly acidic and a cooler temperature exists. These adaptive features are rationalized by the long existence of IBV in humans and may have broader relevance for understanding the biology and evolution of respiratory viruses.

Visualization and Sequencing of Membrane Remodeling Leading to Influenza Virus Fusion

Protein-mediated membrane fusion is an essential step in many fundamental biological events, including enveloped virus infection. The nature of protein and membrane intermediates and the sequence of membrane remodeling during these essential processes remain poorly understood. Here we used cryo-electron tomography (cryo-ET) to image the interplay between influenza virus and vesicles with a range of lipid compositions. By following the population kinetics of membrane fusion intermediates imaged by cryo-ET, we found that membrane remodeling commenced with the hemagglutinin fusion protein spikes grappling onto the target membrane, followed by localized target membrane dimpling as local clusters of hemagglutinin started to undergo conformational refolding. The local dimples then transitioned to extended, tightly apposed contact zones where the two proximal membrane leaflets were in most cases indistinguishable from each other, suggesting significant dehydration and possible intermingling of the lipid head groups. Increasing the content of fusion-enhancing cholesterol or bis-monoacylglycerophosphate in the target membrane led to an increase in extended contact zone formation. Interestingly, hemifused intermediates were found to be extremely rare in the influenza virus fusion system studied here, most likely reflecting the instability of this state and its rapid conversion to postfusion complexes, which increased in population over time. By tracking the populations of fusion complexes over time, the architecture and sequence of membrane reorganization leading to efficient enveloped virus fusion were thus resolved.

IMPORTANCE

Enveloped viruses employ specialized surface proteins to mediate fusion of cellular and viral membranes that results in the formation of pores through which the viral genetic material is delivered to the cell. For influenza virus, the trimeric hemagglutinin (HA) glycoprotein spike mediates host cell attachment and membrane fusion. While structures of a subset of conformations and parts of the fusion machinery have been characterized, the nature and sequence of membrane deformations during fusion have largely eluded characterization. Building upon studies that focused on early stages of HA-mediated membrane remodeling, here cryo-electron tomography (cryo-ET) was used to image the three-dimensional organization of intact influenza virions at different stages of fusion with liposomes, leading all the way to completion of the fusion reaction. By monitoring the evolution of fusion intermediate populations over the course of acid-induced fusion, we identified the progression of membrane reorganization that leads to efficient fusion by an enveloped virus.

Environmental role in influenza virus outbreaks

The environmental drivers of influenza outbreaks are largely unknown. Despite more than 50 years of research, there are conflicting lines of evidence on the role of the environment in influenza A virus (IAV) survival, stability, and transmissibility. With the increasing and looming threat of pandemic influenza, it is important to understand these factors for early intervention and long-term control strategies. The factors that dictate the severity and spread of influenza would include the virus, natural and acquired hosts, virus-host interactions, environmental persistence, virus stability and transmissibility, and anthropogenic interventions. Virus persistence in different environments is subject to minor variations in temperature, humidity, pH, salinity, air pollution, and solar radiations. Seasonality of influenza is largely dictated by temperature and humidity, with cool-dry conditions enhancing IAV survival and transmissibility in temperate climates in high latitudes, whereas humid-rainy conditions favor outbreaks in low latitudes, as seen in tropical and subtropical zones. In mid-latitudes, semiannual outbreaks result from alternating cool-dry and humid-rainy conditions. The mechanism of virus survival in the cool-dry or humid-rainy conditions is largely determined by the presence of salts and proteins in the respiratory droplets. Social determinants of heath, including health equity, vaccine acceptance, and age-related illness, may play a role in influenza occurrence and spread.

Architecture of a nascent viral fusion pore

Enveloped viruses use specialized protein machinery to fuse the viral membrane with that of the host cell during cell invasion. In influenza virus, hundreds of copies of the haemagglutinin (HA) fusion glycoprotein project from the virus surface. Despite intensive study of HA and its fusion activity, the protein's modus operandi in manipulating viral and target membranes to catalyse their fusion is poorly understood. Here, the three-dimensional architecture of influenza virus-liposome complexes at pH 5.5 was investigated by electron cryo-tomography. Tomographic reconstructions show that early stages of membrane remodeling take place in a target membrane-centric manner, progressing from punctate dimples, to the formation of a pinched liposomal funnel that may impinge on the apparently unperturbed viral envelope. The results suggest that the M1 matrix layer serves as an endoskeleton for the virus and a foundation for HA during membrane fusion. Fluorescence spectroscopy monitoring fusion between liposomes and virions shows that leakage of liposome contents takes place more rapidly than lipid mixing at pH 5.5. The relation of 'leaky' fusion to the observed prefusion structures is discussed.

Avian influenza virus (H5N1); effects of physico-chemical factors on its survival

Present study was performed to determine the effects of physical and chemical agents on infective potential of highly pathogenic avian influenza (HPAI) H5N1 (local strain) virus recently isolated in Pakistan during 2006 outbreak. H5N1 virus having titer 10^8.3 ELD₅₀/ml was mixed with sterilized peptone water to get final dilution of 4HA units and then exposed to physical (temperature, pH and ultraviolet light) and chemical (formalin, phenol crystals, iodine crystals, CID 20, virkon^®-S, zeptin 10%, KEPCIDE 300, KEPCIDE 400, lifebuoy, surf excel and caustic soda) agents. Harvested amnio-allantoic fluid (AAF) from embryonated chicken eggs inoculated with H5N1 treated virus (0.2 ml/egg) was subjected to haemagglutination (HA) and haemagglutination inhibition (HI) tests. H5N1 virus lost infectivity after 30 min at 56°C, after 1 day at 28°C but remained viable for more than 100 days at 4°C. Acidic pH (1, 3) and basic pH (11, 13) were virucidal after 6 h contact time; however virus retained infectivity at pH 5 (18 h), 7 and 9 (more than 24 h). UV light was proved ineffectual in inactivating virus completely even after 60 min. Soap (lifebuoy^®), detergent (surf excel^®) and alkali (caustic soda) destroyed infectivity after 5 min at 0.1, 0.2 and 0.3% dilution. All commercially available disinfectants inactivated virus at recommended concentrations. Results of present study would be helpful in implementing bio-security measures at farms/hatcheries levels in the wake of avian influenza virus (AIV) outbreak.

Single-particle kinetics of influenza virus membrane fusion

Membrane fusion is an essential step during entry of enveloped viruses into cells. Conventional fusion assays are generally limited to observation of ensembles of multiple fusion events, confounding more detailed analysis of the sequence of the molecular steps involved. We have developed an in vitro, two-color fluorescence assay to monitor kinetics of single virus particles fusing with a target bilayer on an essentially fluid support. Analysis of lipid- and content-mixing trajectories on a particle-by-particle basis provides evidence for multiple, long-lived kinetic intermediates leading to hemifusion, followed by a single, rate-limiting step to pore formation. We interpret the series of intermediates preceding hemifusion as a result of the requirement that multiple copies of the trimeric hemagglutinin fusion protein be activated to initiate the fusion process.

Inactivation of avian influenza viruses by chemical agents and physical conditions: a review

The recent outbreaks of avian influenza (AI) worldwide have highlighted the difficulties in controlling this disease both in developed and in developing countries. Biosecurity is considered the most important tool to prevent and control AI. In certain areas of the world, AI has become endemic and the recent outbreaks in Europe and Africa show that the epidemiological situation is evolving in an unprecedented way. The consequences of this situation are economic losses to the poultry industry, food security issues in developing countries and a serious threat to human health, due to the direct consequences of AI infection in humans, and more alarmingly due to the risk of the generation of a new pandemic virus from the animal reservoir. In this paper, the physical and chemical methods of inactivating AI viruses are reviewed, with particular emphasis on the practicalities of using such methods in the poultry industry.

The relevance of salt bridges for the stability of the influenza virus hemagglutinin

Hemagglutinin (HA) of influenza virus undergoes an irreversible conformational change at acidic pH, mediating viral fusion with the host endosomal membrane. To unravel the molecular basis of the pH-dependent stability of HA, we demonstrate by mutagenesis of the prototype HA of virus strain X31 (H3 subtype) that salt bridges, especially a tetrad salt bridge within the monomers, are crucial for folding and stability of the trimeric ectodomain. This complex (tetrad) salt bridge is highly conserved among influenza virus subtypes. Introducing additional sites of electrostatic attraction between monomers in the distal region enhanced the stability of ectodomain at low pH mimicking the natural variant H2 subtype. We propose that distinct salt bridges in the distal domain may contribute to the enhanced stability of HA of natural virus variants. This hypothesis may provide clues to understanding adaptations of virus strains (for example, avian influenza viruses) in order to preserve stability of the protein in the host-specific environment.

Paramyxovirus fusion: real-time measurement of parainfluenza virus 5 virus-cell fusion

Although cell-cell fusion assays are useful surrogate methods for studying virus fusion, differences between cell-cell and virus-cell fusion exist. To examine paramyxovirus fusion in real time, we labeled viruses with fluorescent lipid probes and monitored virus-cell fusion by fluorimetry. Two parainfluenza virus 5 (PIV5) isolates (W3A and SER) and PIV5 containing mutations within the fusion protein (F) were studied. Fusion was specific and temperature-dependent. Compared to many low pH-dependent viruses, the kinetics of PIV5 fusion was slow, approaching completion within several minutes. As predicted from cell-cell fusion assays, virus containing an F protein with an extended cytoplasmic tail (rSV5 F551) had reduced fusion compared to wild-type virus (W3A). In contrast, virus-cell fusion for SER occurred at near wild-type levels, despite the fact that this isolate exhibits a severely reduced cell-cell fusion phenotype. These results support the notion that virus-cell and cell-cell fusion have significant differences.

Stoichiometry of envelope glycoprotein trimers in the entry of human immunodeficiency virus type 1

The human immunodeficiency virus type 1 (HIV-1) envelope glycoproteins (Envs) function as a trimer, mediating virus entry by promoting the fusion of the viral and target cell membranes. HIV-1 Env trimers induce membrane fusion through a pH-independent pathway driven by the interaction between an Env trimer and its cellular receptors, CD4 and CCR5/CXCR4. We studied viruses with mixed heterotrimers of wild-type and dominant-negative Envs to determine the number (T) of Env trimers required for HIV-1 entry. To our surprise, we found that a single Env trimer is capable of supporting HIV-1 entry; i.e., T = 1. A similar approach was applied to investigate the entry stoichiometry of envelope glycoproteins from amphotropic murine leukemia virus (A-MLV), avian sarcoma/leukosis virus type A (ASLV-A), and influenza A virus. When pseudotyped on HIV-1 virions, the A-MLV and ASLV-A Envs also exhibit a T = 1 entry stoichiometry. In contrast, eight to nine influenza A virus hemagglutinin trimers function cooperatively to achieve membrane fusion and virus entry, using a pH-dependent pathway. The different entry requirements for cooperativity among Env trimers for retroviruses and influenza A virus may influence viral strategies for replication and evasion of the immune system.

Exact kinetic analysis of passive transport across a polarized confluent MDCK cell monolayer modeled as a single barrier

Knowledge of the passive permeability coefficient for new drugs is useful for estimating the fraction absorbed across the gastrointestinal tract. The commonly used approximate formula for the passive permeability coefficient is based on the initial rate of permeation across cell monolayers, requires measurement during the linear phase of permeation, and is not applicable when there is significant back flux of compound or mass balance problem. To develop a rigorous equation that can be used at any time point, i.e., that is valid outside of the linear phase, the mass action equations were integrated for a standard single barrier model of passive permeability. The simple analytical solution found also allows correction for both loss of drug (e.g., due to binding and/or hydrolysis) and sampling volume loss for multiple time point experiments. To test this equation, we measured the passive permeation of three well characterized drugs (amprenavir, quinidine, and loperamide) across confluent monolayers of MDCKII-hMDR1 cells. The potent P-glycoprotein inhibitor GF120918 was used to inhibit P-glycoprotein activity, so only passive permeability was determined. Dramatically different time-dependent behavior was observed for the three compounds, with loperamide showing significant loss of compound, and loperamide and quinidine causing plasma membrane modifications over time. The simple and exact equation for the permeability coefficient developed here works from start of transport to equilibrium, being valid when the commonly used approximate equation may not be. Thus, the exact equation is safer to use in any context, even for single time point estimates in high-throughput permeability assays.

The elementary mass action rate constants of P-gp transport for a confluent monolayer of MDCKII-hMDR1 cells

The human multi-drug resistance membrane transporter, P-glycoprotein, or P-gp, has been extensively studied due to its importance to human health and disease. Thus far, the kinetic analysis of P-gp transport has been limited to steady-state Michaelis-Menten approaches or to compartmental models, neither of which can prove molecular mechanisms. Determination of the elementary kinetic rate constants of transport will be essential to understanding how P-gp works. The experimental system we use is a confluent monolayer of MDCKII-hMDR1 cells that overexpress P-gp. It is a physiologically relevant model system, and transport is measured without biochemical manipulations of P-gp. The Michaelis-Menten mass action reaction is used to model P-gp transport. Without imposing the steady-state assumptions, this reaction depends upon several parameters that must be simultaneously fitted. An exhaustive fitting of transport data to find all possible parameter vectors that best fit the data was accomplished with a reasonable computation time using a hierarchical algorithm. For three P-gp substrates (amprenavir, loperamide, and quinidine), we have successfully fitted the elementary rate constants, i.e., drug association to P-gp from the apical membrane inner monolayer, drug dissociation back into the apical membrane inner monolayer, and drug efflux from P-gp into the apical chamber, as well as the density of efflux active P-gp. All three drugs had overlapping ranges for the efflux active P-gp, which was a benchmark for the validity of the fitting process. One novel finding was that the association to P-gp appears to be rate-limited solely by drug lateral diffusion within the inner monolayer of the plasma membrane for all three drugs. This would be expected if P-gp structure were open to the lipids of the apical membrane inner monolayer, as has been suggested by recent structural studies. The fitted kinetic parameters show how P-gp efflux of a wide range of xenobiotics has been maximized.

Kinetically differentiating influenza hemagglutinin fusion and hemifusion machines

Membrane fusion mediated by influenza virus hemagglutinin (HA) yields different phenotypes depending on the surface density of activated HAs. A key question is whether different phenotypes arise from different fusion machines or whether different numbers of identical fusion machines yield different probabilistic outcomes. If fusion were simply a less probable event than hemifusion, requiring a larger number of identical fusion machines to occur first, then two predictions can be made. First, fusion should have a shorter average delay time than hemifusion, since there are more machines. Second, fusion should have a longer execution time of lipid mixing after it begins than hemifusion, since the full event cannot be faster than the partial event. Using a new automated video microscopy technique, we simultaneously monitored many HA-expressing cells fusing with erythrocytes and identified individual cell pairs with either full or only partial redistribution of fluorescent lipids. The full lipid mixing phenotype also showed contents mixing, i.e., fusion. Kinetic screening of the digitized fluorescence data showed that the execution of lipid mixing after the onset is faster for fusion than hemifusion. We found no correlation between the delay times before the onset of lipid mixing and the final fusion phenotype. We also found that the execution time for fusion was faster than that for hemifusion. Thus, we provide the first experimental evidence for fusion and hemifusion arising from different machines.

Architecture of the influenza hemagglutinin membrane fusion site

The mechanism of influenza hemagglutinin (HA) mediated membrane fusion has been intensively studied for over 20 years after the bromelain-released ectodomain of HA at neutral pH was first crystallized. Nearly 10 years ago, the low-pH-induced "spring coiled" conformational change of HA was predicted from peptide chemistry and confirmed by crystallography. Other work has yielded a wealth of knowledge on the observed changes in HA fusion/hemifusion phenotypes as a function of site-specific mutations of HA, or added amphipathic molecules or particular IgGs. It is becoming clear that the conformational changes predicted by the crystallography are necessary to cause fusion and that interfering with these changes can block fusion or reduce it to hemifusion. What is not known is how the conformational changes cause fusion. In particular, while it is generally agreed that fusion requires an aggregate of HAs, how the aggregate may act to transduce the energy of the HA conformational changes to creating the initial fusion defect is not known. We have used a comprehensive mass action kinetic model of HA-mediated fusion to carry out a "meta-analysis" of several key data sets, using HA-expressing cells and using virions. The consensus result of these detailed kinetic studies was that the fusion site of influenza hemagglutinin (HA) is an aggregate with at least eight HAs. The high-energy conformational change of only two of these HAs within the aggregate permits the formation of the first fusion pore. This "8 and 2" result was required to best fit all the data. We review these studies and how this kinetic result can guide and constrain HA fusion models. The kinetic analysis suggests that the sequence of fusion intermediates starts with protein control and ends with lipid control, which makes sense. While curvature intermediates, e.g. the lipid stalk, are almost certainly within the fusion sequence, the "8 and 2" result does not suggest that they are the first step after HA aggregation. The stabilized hydrophobic defect model we have proposed as a precursor to the lipid stalk can form and is consistent with the "8 and 2" result.

Molecular dynamics simulation of the evolution of hydrophobic defects in one monolayer of a phosphatidylcholine bilayer: relevance for membrane fusion mechanisms

The spontaneous formation of the phospholipid bilayer underlies the permeability barrier function of the biological membrane. Tears or defects that expose water to the acyl chains are spontaneously healed by lipid lateral diffusion. However, mechanical barriers, e.g., protein aggregates held in place, could sustain hydrophobic defects. Such defects have been postulated to occur in processes such as membrane fusion. This gives rise to a new question in bilayer structure: What do the lipids do in the absence of lipid lateral diffusion to minimize the free energy of a hydrophobic defect? As a first step to understand this rather fundamental question about bilayer structure, we performed molecular dynamic simulations of up to 10 ns of a planar bilayer from which lipids have been deleted randomly from one monolayer. In one set of simulations, approximately one-half of the lipids in the defect monolayer were restrained to form a mechanical barrier. In the second set, lipids were free to diffuse around. The question was simply whether the defects caused by removing a lipid would aggregate together, forming a large hydrophobic cavity, or whether the membrane would adjust in another way. When there are no mechanical barriers, the lipids in the defect monolayer simply spread out and thin with little effect on the other intact monolayer. In the presence of a mechanical barrier, the behavior of the lipids depends on the size of the defect. When 3 of 64 lipids are removed, the remaining lipids adjust the lower one-half of their chains, but the headgroup structure changes little and the intact monolayer is unaffected. When 6 to 12 lipids are removed, the defect monolayer thins, lipid disorder increases, and lipids from the intact monolayer move toward the defect monolayer. Whereas this is a highly simplified model of a fusion site, this engagement of the intact monolayer into the fusion defect is strikingly consistent with recent results for influenza hemagglutinin mediated fusion.