Minor electrostatic changes robustly increase VP40 membrane binding, assembly, and budding of Ebola virus matrix protein derived virus-like particles.

Ebola virus (EBOV) is a filamentous negative-sense RNA virus, which causes severe hemorrhagic fever. There are limited vaccines or therapeutics for prevention and treatment of EBOV, so it is important to get a detailed understanding of the virus lifecycle to illuminate new drug targets. EBOV encodes for the matrix protein, VP40, which regulates assembly and budding of new virions from the inner leaflet of the host cell plasma membrane (PM). In this work, we determine the effects of VP40 mutations altering electrostatics on PM interactions and subsequent budding. VP40 mutations that modify surface electrostatics affect viral assembly and budding by altering VP40 membrane-binding capabilities. Mutations that increase VP40 net positive charge by one (e.g., Gly to Arg or Asp to Ala) increase VP40 affinity for phosphatidylserine and phosphatidylinositol 4,5-bisphosphate in the host cell PM. This increased affinity enhances PM association and budding efficiency leading to more effective formation of virus-like particles. In contrast, mutations that decrease net positive charge by one (e.g., Gly to Asp) lead to a decrease in assembly and budding because of decreased interactions with the anionic PM. Taken together, our results highlight the sensitivity of slight electrostatic changes on the VP40 surface for assembly and budding. Understanding the effects of single amino acid substitutions on viral budding and assembly will be useful for explaining changes in the infectivity and virulence of different EBOV strains, VP40 variants that occur in nature, and for long-term drug discovery endeavors aimed at EBOV assembly and budding.

PI(4,5)P2 binding sites in the Ebola virus matrix protein VP40 modulate assembly and budding.

Ebola virus (EBOV) causes severe hemorrhagic fever in humans and is lethal in a large percentage of those infected. The EBOV matrix protein viral protein 40 kDa (VP40) is a peripheral binding protein that forms a shell beneath the lipid bilayer in virions and virus-like particles (VLPs). VP40 is required for virus assembly and budding from the host cell plasma membrane. VP40 is a dimer that can rearrange into oligomers at the plasma membrane interface, but it is unclear how these structures form and how they are stabilized. We therefore investigated the ability of VP40 to form stable oligomers using in vitro and cellular assays. We characterized two lysine-rich regions in the VP40 C-terminal domain (CTD) that bind phosphatidylinositol-4,5-bisphosphate (PI(4,5)P2) and play distinct roles in lipid binding and the assembly of the EBOV matrix layer. The extensive analysis of VP40 with and without lipids by hydrogen deuterium exchange mass spectrometry revealed that VP40 oligomers become extremely stable when VP40 binds PI(4,5)P2. The PI(4,5)P2-induced stability of VP40 dimers and oligomers is a critical factor in VP40 oligomerization and release of VLPs from the plasma membrane. The two lysine-rich regions of the VP40 CTD have different roles with respect to interactions with plasma membrane phosphatidylserine (PS) and PI(4,5)P2. CTD region 1 (Lys221, Lys224, and Lys225) interacts with PI(4,5)P2 more favorably than PS and is important for VP40 extent of oligomerization. In contrast, region 2 (Lys270, Lys274, Lys275, and Lys279) mediates VP40 oligomer stability via lipid interactions and has a more prominent role in release of VLPs.

AT-hook peptides bind the major and minor groove of AT-rich DNA duplexes.

The mammalian high mobility group protein AT-hook 2 (HMGA2) houses three motifs that preferentially bind short stretches of AT-rich DNA regions. These DNA binding motifs, known as 'AT-hooks', are traditionally characterized as being unstructured. Upon binding to AT-rich DNA, they form ordered assemblies. It is this disordered-to-ordered transition that has implicated HMGA2 as a protein actively involved in many biological processes, with abnormal HMGA expression linked to a variety of health problems including diabetes, obesity, and oncogenesis. In the current work, the solution binding dynamics of the three 'AT-hook' peptides (ATHPs) with AT-rich DNA hairpin substrates were studied using DNA UV melting studies, fluorescence spectroscopy, native ion mobility spectrometry-mass spectrometry (IMS-MS), solution isothermal titration calorimetry (ITC) and molecular modeling. Results showed that the ATHPs bind to the DNA to form a single, 1:1 and 2:1, 'key-locked' conformational ensemble. The molecular models showed that 1:1 and 2:1 complex formation is driven by the capacity of the ATHPs to bind to the minor and major grooves of the AT-rich DNA oligomers. Complementary solution ITC results confirmed that the 2:1 stoichiometry of ATHP: DNA is originated under native conditions in solution.

A 5'-Flanking C/G Pair at the Core Region Enhances the Recognition and Binding of Kaiso to Methylated DNA.

Methyl CpG binding proteins (MBPs) are transcription factors that recognize the methylated CpG sites in DNA and mediate the DNA methylation signal into various downstream cellular processes. The C2H2 zinc finger (ZF) protein, Kaiso, also an MBP, preferentially binds to two symmetrically methylated CpG sites in DNA sequences via C-terminal C2H2 ZF domains and mediates the transcription regulation process. Investigation of the molecular mechanism of the recognition of methylated DNA (meDNA) by Kaiso is important to understand how this protein reads and translates this methylation signal into downstream transcription outcomes. Despite previous studies in Kaiso-meDNA interactions, detailed structural investigations on the sequence-specific interaction of Kaiso with the meDNA sequence are still lacking. In this work, we used molecular modeling and molecular dynamics (MD) simulation-based computational approaches to investigate the recognition of various methylated DNA sequences by Kaiso. Our MD simulation results show that the Kaiso-meDNA interaction is sequence specific. The recognition of meDNA by Kaiso is enhanced in the MeECad sequence compared to the MeCG2 sequence. Compared to the 5'-flanking T/A pair in MeCG2, both MeCG2_mutCG and MeECad sequences show that a C/G base pair allows GLU535 of Kaiso to preferably recognize and bind the core mCpG site. The core mCGmCG site is crucial for the recognition process and formation of a stable complex. Our results reveal that the 5'-flanking nucleotides are also important for the enhanced binding and recognition of methylated sites.

Lipid-specific oligomerization of the Marburg virus matrix protein VP40 is regulated by two distinct interfaces for virion assembly.

Marburg virus (MARV) is a lipid-enveloped virus harboring a negative-sense RNA genome, which has caused sporadic outbreaks of viral hemorrhagic fever in sub-Saharan Africa. MARV assembles and buds from the host cell plasma membrane where MARV matrix protein (mVP40) dimers associate with anionic lipids at the plasma membrane inner leaflet and undergo a dynamic and extensive self-oligomerization into the structural matrix layer. The MARV matrix layer confers the virion filamentous shape and stability but how host lipids modulate mVP40 oligomerization is mostly unknown. Using in vitro and cellular techniques, we present a mVP40 assembly model highlighting two distinct oligomerization interfaces: the (N-terminal domain [NTD] and C-terminal domain [CTD]) in mVP40. Cellular studies of NTD and CTD oligomerization interface mutants demonstrate the importance of each interface in matrix assembly. The assembly steps include protein trafficking to the plasma membrane, homo-multimerization that induced protein enrichment, plasma membrane fluidity changes, and elongations at the plasma membrane. An ascorbate peroxidase derivative (APEX)-transmission electron microscopy method was employed to closely assess the ultrastructural localization and formation of viral particles for wildtype mVP40 and NTD and CTD oligomerization interface mutants. Taken together, these studies present a mechanistic model of mVP40 oligomerization and assembly at the plasma membrane during virion assembly that requires interactions with phosphatidylserine for NTD-NTD interactions and phosphatidylinositol-4,5-bisphosphate for proper CTD-CTD interactions. These findings have broader implications in understanding budding of lipid-enveloped viruses from the host cell plasma membrane and potential strategies to target protein-protein or lipid-protein interactions to inhibit virus budding.

Ebola virus protein VP40 binding to Sec24c for transport to the plasma membrane.

Outbreaks of the Ebola virus (EBOV) continue to occur and while a vaccine and treatment are now available, there remains a dearth of options for those who become sick with EBOV disease. An understanding at the atomic and molecular level of the various steps in the EBOV replication cycle can provide molecular targets for disrupting the virus. An important step in the EBOV replication cycle is the transport of EBOV structural matrix VP40 protein molecules to the plasma membrane inner leaflet, which involves VP40 binding to the host cell's Sec24c protein. Though some VP40 residues involved in the binding are known, the molecular details of VP40-Sec24c binding are not known. We use various molecular computational techniques to investigate the molecular details of how EBOV VP40 binds with the Sec24c complex of the ESCRT-I pathway. We employed different docking programs to identify the VP40-binding site on Sec24c and then performed molecular dynamics simulations to determine the atomic details and binding interactions of the complex. We also investigated how the inter-protein interactions of the complex are affected upon mutations of VP40 amino acids in the Sec24c-binding region. Our results provide a molecular basis for understanding previous coimmunoprecipitation experimental studies. In addition, we found that VP40 can bind to a site on Sec24c that can also bind Sec23 and suggests that VP40 may use the COPII transport mechanism in a manner that may not need the Sec23 protein in order for VP40 to be transported to the plasma membrane.

Mutation-induced changes in the receptor-binding interface of the SARS-CoV-2 Delta variant B.1.617.2 and implications for immune evasion.

Following the initial surges of the Alpha (B.1.1.7) and the Beta (B.1.351) variants, a more infectious Delta variant (B.1.617.2) is now surging, further deepening the health crises caused by the pandemic. The sharp rise in cases attributed to the Delta variant has made it especially disturbing and is a variant of concern. Fortunately, current vaccines offer protection against known variants of concern, including the Delta variant. However, the Delta variant has exhibited some ability to dodge the immune system as it is found that neutralizing antibodies from prior infections or vaccines are less receptive to binding with the Delta spike protein. Here, we investigated the structural changes caused by the mutations in the Delta variant's receptor-binding interface and explored the effects on binding with the ACE2 receptor as well as with neutralizing antibodies. We find that the receptor-binding ß-loop-ß motif adopts an altered but stable conformation causing separation in some of the antibody binding epitopes. Our study shows reduced binding of neutralizing antibodies and provides a possible mechanism for the immune evasion exhibited by the Delta variant.

Structural insights into the repair mechanism of AGT for methyl-induced DNA damage.

Methylation induced DNA base-pairing damage is one of the major causes of cancer. O6-alkylguanine-DNA alkyltransferase (AGT) is considered a demethylation agent of the methylated DNA. Structural investigations with thermodynamic properties of the AGT-DNA complex are still lacking. In this report, we modeled two catalytic states of AGT-DNA interactions and an AGT-DNA covalent complex and explored structural features using molecular dynamics (MD) simulations. We utilized the umbrella sampling method to investigate the changes in the free energy of the interactions in two different AGT-DNA catalytic states, one with methylated GUA in DNA and the other with methylated CYS145 in AGT. These non-covalent complexes represent the pre- and post-repair complexes. Therefore, our study encompasses the process of recognition, complex formation, and separation of the AGT and the damaged (methylated) DNA base. We believe that the use of parameters for the amino acid and nucleotide modifications and for the protein-DNA covalent bond will allow investigations of the DNA repair mechanism as well as the exploration of cancer therapeutics targeting the AGT-DNA complexes at various functional states as well as explorations via stabilization of the complex.

A cationic, C-terminal patch and structural rearrangements in Ebola virus matrix VP40 protein control its interactions with phosphatidylserine.

Ebola virus (EBOV) is a filamentous lipid-enveloped virus that causes hemorrhagic fever with a high fatality rate. Viral protein 40 (VP40) is the major EBOV matrix protein and regulates viral budding from the plasma membrane. VP40 is a transformer/morpheein that can structurally rearrange its native homodimer into either a hexameric filament that facilitates viral budding or an RNA-binding octameric ring that regulates viral transcription. VP40 associates with plasma-membrane lipids such as phosphatidylserine (PS), and this association is critical to budding from the host cell. However, it is poorly understood how different VP40 structures interact with PS, what essential residues are involved in this association, and whether VP40 has true selectivity for PS among different glycerophospholipid headgroups. In this study, we used lipid-binding assays, MD simulations, and cellular imaging to investigate the molecular basis of VP40-PS interactions and to determine whether different VP40 structures (i.e. monomer, dimer, and octamer) can interact with PS-containing membranes. Results from quantitative analysis indicated that VP40 associates with PS vesicles via a cationic patch in the C-terminal domain (Lys224, 225 and Lys274, 275). Substitutions of these residues with alanine reduced PS-vesicle binding by >40-fold and abrogated VP40 localization to the plasma membrane. Dimeric VP40 had 2-fold greater affinity for PS-containing membranes than the monomer, whereas binding of the VP40 octameric ring was reduced by nearly 10-fold. Taken together, these results suggest the different VP40 structures known to form in the viral life cycle harbor different affinities for PS-containing membranes.

Membrane pore formation and ion selectivity of the Ebola virus delta peptide.

The Ebola virus (EBOV) is a filamentous lipid-enveloped virus that causes severe hemorrhagic fever with a high fatality rate in humans. The EBOV encodes a glycoprotein that when cleaved, produces the delta peptide. Experimental evidence suggests that the delta peptide functions as a viroporin that enhances virus particle release through the host cell membrane. However, the viroporin forming mechanism of the delta peptide is still not well understood. Guided by experimental information, we have computationally investigated the pore formation by different oligomers of the delta peptide. We have performed all-atom molecular dynamics (MD) simulations in an explicit membrane environment to investigate the pore-forming mechanism and stability of the pores. Our results suggest that the delta peptide forms stable pentameric pores. In addition, the pore is selective with respect to chloride ions, and the disulfide bond formed between Cys-29 and Cys-38 in the C-terminal of the peptide is essential for the pore stabilization and ion permeation. Our study provides helpful information on the pore-forming mechanism of filovirus delta peptides and such structural information can be important in designing and developing molecular modulators that target the delta peptide pore and disrupt the pathology of the Ebola virus.

Molecular mechanisms of pore formation and membrane disruption by the antimicrobial lantibiotic peptide Mutacin 1140.

The emergence of antibiotic-resistance is a major concern to global human health and identification of novel antibiotics is critical to mitigate the threat. Mutacin 1140 (MU1140) is a promising antimicrobial lanthipeptide and is effective against Gram-positive bacteria. Like nisin, MU1140 targets and sequesters lipid II and interferes with its function, which results in the inhibition of bacterial cell wall synthesis, and leads to bacteria cell lysis. MU1140 contains a structurally similar thioether cage for binding the lipid II pyrophosphate as for nisin. In addition to lipid II binding, nisin is known to form membrane pores. Membrane pore formation and membrane disruption is a common mode of action for many antimicrobial peptides, including gallidermin, a lantibiotic peptide with similar structural features as MU1140. However, whether and how MU1140 and its variants can form permeable membrane pores remains to be demonstrated. In this work, we explored the potential mechanisms of membrane pore formation by performing molecular simulations of the MU1140-lipid II complex in the bacterial membrane. Our results suggest that MU1140-lipid II complexes are able to form water permeating membrane pores. We find that a single chain of MU1140 complexed with lipid II in the transmembrane region can permeate water molecules across the membrane via a single-file water transport mechanism. The ordering of the water molecules in the single-file chain region as well as the diffusion behavior is similar to those observed in other biological water channels. Multiple complexes of MU1140-lipid II in the membrane showed enhanced permeability for the water molecules, as well as a noticeable membrane distortion and lipid relocation, suggesting that a higher concentration of MU1140 assembly in the membrane can cause significant disruption of the bacterial membrane. These investigations provide an atomistic level insight into a novel mode of action for MU1140 that can be exploited to develop optimized peptide variants with improved antimicrobial properties.

Detection of lipid-induced structural changes of the Marburg virus matrix protein VP40 using hydrogen/deuterium exchange-mass spectrometry.

Marburg virus (MARV) is a lipid-enveloped virus from the Filoviridae family containing a negative sense RNA genome. One of the seven MARV genes encodes the matrix protein VP40, which forms a matrix layer beneath the plasma membrane inner leaflet to facilitate budding from the host cell. MARV VP40 (mVP40) has been shown to be a dimeric peripheral protein with a broad and flat basic surface that can associate with anionic phospholipids such as phosphatidylserine. Although a number of mVP40 cationic residues have been shown to facilitate binding to membranes containing anionic lipids, much less is known on how mVP40 assembles to form the matrix layer following membrane binding. Here we have used hydrogen/deuterium exchange (HDX) mass spectrometry to determine the solvent accessibility of mVP40 residues in the absence and presence of phosphatidylserine and phosphatidylinositol 4,5-bisphosphate. HDX analysis demonstrates that two basic loops in the mVP40 C-terminal domain make important contributions to anionic membrane binding and also reveals a potential oligomerization interface in the C-terminal domain as well as a conserved oligomerization interface in the mVP40 N-terminal domain. Lipid binding assays confirm the role of the two basic patches elucidated with HD/X measurements, whereas molecular dynamics simulations and membrane insertion measurements complement these studies to demonstrate that mVP40 does not appreciably insert into the hydrocarbon region of anionic membranes in contrast to the matrix protein from Ebola virus. Taken together, we propose a model by which association of the mVP40 dimer with the anionic plasma membrane facilitates assembly of mVP40 oligomers.

Membrane association and localization dynamics of the Ebola virus matrix protein VP40.

The Ebola virus matrix protein VP40 is a major structural protein that provides the scaffolding for new Ebola virus particles. For this, VP40 is first trafficked to the lower leaflet of the plasma membrane (PM) in its dimeric form. Once associated with the PM, the VP40 dimers undergo structural rearrangements and oligomerize into hexamers and filaments that make up the virus matrix. Therefore, association of the VP40 dimers and their stabilization at the PM is a crucial step in the Ebola life-cycle. To understand the molecular details of the VP40 dimer-PM interactions, we investigated the dimer association with the inner leaflet of the PM using detailed all-atom molecular dynamics (MD) simulations. The formation of the dimer-PM complex is facilitated by the interactions of the VP40 lysine residues and the anionic lipids POPS, POPI, and PIP2 in the PM. In contrast, the dimer fails to associate with a membrane without POPS, POPI, or PIP2 lipids. We explored the mechanisms of the association and identified important residues and lipids involved in localization and stabilization of VP40 dimers at the PM. MD simulations elucidate the role of a C-terminal α-helix alignment parallel to the lipid bilayer surface as well as the creation of membrane defects that allow partial insertion of the hydrophobic residue V276 into the membrane to further stabilize the VP40 dimer-PM complex. Understanding the mechanisms of the VP40 dimer-PM association that facilitate oligomerization can be important for potentially targeting the VP40 for small molecules that can interfere with the virus life-cycle.

Graphene-VP40 interactions and potential disruption of the Ebola virus matrix filaments.

Ebola virus infections cause hemorrhagic fever that often results in very high fatality rates. In addition to exploring vaccines, development of drugs is also essential for treating the disease and preventing the spread of the infection. The Ebola virus matrix protein VP40 exists in various conformational and oligomeric forms and is a potential pharmacological target for disrupting the virus life-cycle. Here we explored graphene-VP40 interactions using molecular dynamics simulations and graphene pelleting assays. We found that graphene sheets associate strongly with VP40 at various interfaces. We also found that the graphene is able to disrupt the C-terminal domain (CTD-CTD) interface of VP40 hexamers. This VP40 hexamer-hexamer interface is crucial in forming the Ebola viral matrix and disruption of this interface may provide a method to use graphene or similar nanoparticle based solutions as a disinfectant that can significantly reduce the spread of the disease and prevent an Ebola epidemic.

Cooperative structural transitions in amyloid-like aggregation.

Amyloid fibril aggregation is associated with several horrific diseases such as Alzheimer's, Creutzfeld-Jacob, diabetes, Parkinson's, and others. Although proteins that undergo aggregation vary widely in their primary structure, they all produce a cross-ß motif with the proteins in ß-strand conformations perpendicular to the fibril axis. The process of amyloid aggregation involves forming myriad different metastable intermediate aggregates. To better understand the molecular basis of the protein structural transitions and aggregation, we report on molecular dynamics (MD) computational studies on the formation of amyloid protofibrillar structures in the small model protein ccß, which undergoes many of the structural transitions of the larger, naturally occurring amyloid forming proteins. Two different structural transition processes involving hydrogen bonds are observed for aggregation into fibrils: the breaking of intrachain hydrogen bonds to allow ß-hairpin proteins to straighten, and the subsequent formation of interchain H-bonds during aggregation into amyloid fibrils. For our MD simulations, we found that the temperature dependence of these two different structural transition processes results in the existence of a temperature window that the ccß protein experiences during the process of forming protofibrillar structures. This temperature dependence allows us to investigate the dynamics on a molecular level. We report on the thermodynamics and cooperativity of the transformations. The structural transitions that occurred in a specific temperature window for ccß in our investigations may also occur in other amyloid forming proteins but with biochemical parameters controlling the dynamics rather than temperature.

The Ebola virus protein VP40 hexamer enhances the clustering of PI(4,5)P2 lipids in the plasma membrane.

The Ebola virus is a lipid-enveloped virus that obtains its lipid coat from the plasma membrane of the host cell it infects during the budding process. The Ebola virus protein VP40 localizes to the inner leaflet of the plasma membrane and forms the viral matrix, which provides the major structure for the Ebola virus particles. VP40 is initially a dimer that rearranges to a hexameric structure that mediates budding. VP40 hexamers and larger filaments have been shown to be stabilized by PI(4,5)P2 in the plasma membrane inner leaflet. Reduction in the plasma membrane levels of PI(4,5)P2 significantly reduce formation of VP40 oligomers and virus-like particles. We investigated the lipid-protein interactions in VP40 hexamers at the plasma membrane. We quantified lipid-lipid self-clustering by calculating the fractional interaction matrix and found that the VP40 hexamer significantly enhances the PI(4,5)P2 clustering. The radial pair distribution functions suggest a strong interaction between PI(4,5)P2 and the VP40 hexamer. The cationic Lys side chains are found to mediate the PIP2 clustering around the protein, with cholesterol filling the space between the interacting PIP2 molecules. These computational studies support recent experimental data and provide new insights into the mechanisms by which VP40 assembles at the plasma membrane inner leaflet, alters membrane curvature, and forms new virus-like particles.

Structures of the kinetically trapped i-motif DNA intermediates.

In the present work, the conformational dynamics and folding pathways of i-motif DNA were studied in solution and in the gas-phase as a function of the solution pH conditions using circular dichroism (CD), photoacoustic calorimetry analysis (PAC), trapped ion mobility spectrometry-mass spectrometry (TIMS-MS), and molecular dynamics (MD). Solution studies showed at thermodynamic equilibrium the existence of a two-state folding mechanism, whereas during the pH = 7.0 → 4.5 transition a fast and slow phase (ΔHfast + ΔHslow = 43 ± 7 kcal mol-1) with a volume change associated with the formation of hemiprotonated cytosine base pairs and concomitant collapse of the i-motif oligonucleotide into a compact conformation were observed. TIMS-MS experiments showed that gas-phase, kinetically trapped i-motif DNA intermediates produced by nanoESI are preserved, with relative abundances depending on the solution pH conditions. In particular, a folded i-motif DNA structure was observed in nanoESI-TIMS-MS for low charge states in both positive and negative ion mode (e.g., z = ±3 to ±5) at low pH conditions. As solution pH increases, the cytosine neutralization leads to the loss of cytosine-cytosine+ (C·CH+) base pairing in the CCC strands and in those conditions we observe partially unfolded i-motif DNA conformations in nanoESI-TIMS-MS for higher charge states (e.g., z = -6 to -9). Collisional induced activation prior to TIMS-MS showed the existence of multiple local free energy minima, associated with the i-motif DNA unfolding at z = -6 charge state. For the first time, candidate gas-phase structures are proposed based on mobility measurements of the i-motif DNA unfolding pathway. Moreover, the inspection of partially unfolded i-motif DNA structures (z = -7 and z = -8 charge states) showed that the presence of inner cations may or may not induce conformational changes in the gas-phase. For example, incorporation of ammonium adducts does not lead to major conformational changes while sodium adducts may lead to the formation of sodium mediated bonds between two negatively charged sides inducing the stabilization towards more compact structures in new local, free energy minima in the gas-phase.

Minor changes in electrostatics robustly increase VP40 membrane binding, assembly, and budding of Ebola virus matrix protein derived virus-like particles.

Ebola virus (EBOV) is a filamentous negative-sense RNA virus which causes severe hemorrhagic fever. There are limited vaccines or therapeutics for prevention and treatment of EBOV, so it is important to get a detailed understanding of the virus lifecycle to illuminate new drug targets. EBOV encodes for the matrix protein, VP40, which regulates assembly and budding of new virions from the inner leaflet of the host cell plasma membrane (PM). In this work we determine the effects of VP40 mutations altering electrostatics on PM interactions and subsequent budding. VP40 mutations that modify surface electrostatics affect viral assembly and budding by altering VP40 membrane binding capabilities. Mutations that increase VP40 net positive charge by one (e.g., Gly to Arg or Asp to Ala) increase VP40 affinity for phosphatidylserine (PS) and PI(4,5)P2 in the host cell PM. This increased affinity enhances PM association and budding efficiency leading to more effective formation of virus-like particles (VLPs). In contrast, mutations that decrease net positive charge by one (e.g., Gly to Asp) lead to a decrease in assembly and budding because of decreased interactions with the anionic PM. Taken together our results highlight the sensitivity of slight electrostatic changes on the VP40 surface for assembly and budding. Understanding the effects of single amino acid substitutions on viral budding and assembly will be useful for explaining changes in the infectivity and virulence of different EBOV strains, VP40 variants that occur in nature, and for long-term drug discovery endeavors aimed at EBOV assembly and budding.

Role of phosphatidic acid lipids on plasma membrane association of the Ebola virus matrix protein VP40.

The Ebola virus matrix protein VP40 is responsible for the formation of the viral matrix by localizing at the inner leaflet of the human plasma membrane (PM). Various lipid types, including PI(4,5)P2 (i.e. PIP2) and phosphatidylserine (PS), play active roles in this process. Specifically, the negatively charged headgroups of both PIP2 and PS interact with the basic residues of VP40 and stabilize it at the membrane surface, allowing for eventual egress. Phosphatidic acid (PA), resulting from the enzyme phospholipase D (PLD), is also known to play an active role in viral development. In this work, we performed a biophysical and computational analysis to investigate the effects of the presence of PA on the membrane localization and association of VP40. We used coarse-grained molecular dynamics simulations to quantify VP40 hexamer interactions with the inner leaflet of the PM. Analysis of the local distribution of lipids shows enhanced lipid clustering when PA is abundant in the membrane. We observed that PA lipids have a similar role to that of PS lipids in VP40 association due to the geometry and charge. Complementary experiments performed in cell culture demonstrate competition between VP40 and a canonical PA-binding protein for the PM. Also, inhibition of PA synthesis reduced the detectable budding of virus-like particles. These computational and experimental results provide new insights into the early stages of Ebola virus budding and the role that PA lipids have on the VP40-PM association.

Computational and experimental identification of keystone interactions in Ebola virus matrix protein VP40 dimer formation.

The Ebola virus (EBOV) is a lipid-enveloped virus with a negative sense RNA genome that can cause severe and often fatal viral hemorrhagic fever. The assembly and budding of EBOV is regulated by the matrix protein, VP40, which is a peripheral protein that associates with anionic lipids at the inner leaflet of the plasma membrane. VP40 is sufficient to form virus-like particles (VLPs) from cells, which are nearly indistinguishable from authentic virions. Due to the restrictions of studying EBOV in BSL-4 facilities, VP40 has served as a surrogate in cellular studies to examine the EBOV assembly and budding process from the host cell plasma membrane. VP40 is a dimer where inhibition of dimer formation halts budding and formation of new VLPs as well as VP40 localization to the plasma membrane inner leaflet. To better understand VP40 dimer stability and critical amino acids to VP40 dimer formation, we integrated computational approaches with experimental validation. Site saturation/alanine scanning calculation, combined with molecular mechanics-based generalized Born with Poisson-Boltzmann surface area (MM-GB/PBSA) method and molecular dynamics simulations were used to predict the energetic contribution of amino acids to VP40 dimer stability and the hydrogen bonding network across the dimer interface. These studies revealed several previously unknown interactions and critical residues predicted to impact VP40 dimer formation. In vitro and cellular studies were then pursued for a subset of VP40 mutations demonstrating reduction in dimer formation (in vitro) or plasma membrane localization (in cells). Together, the computational and experimental approaches revealed critical residues for VP40 dimer stability in an alpha-helical interface (between residues 106-117) as well as in a loop region (between residues 52-61) below this alpha-helical region. This study sheds light on the structural origins of VP40 dimer formation and may inform the design of a small molecule that can disrupt VP40 dimer stability.