Optimizing properties of translocation-enhancing transmembrane proteins

Cell membranes act as semi-permeable barriers, often restricting the entry of large or hydrophilic molecules. Nonetheless, certain amphiphilic molecules, such as antimicrobial and cell-penetrating peptides, can cross these barriers. In this study, we demonstrate that specific properties of transmembrane proteins/peptides can enhance membrane permeation of amphiphilic peptides. Using coarse-grained molecular dynamics with free-energy calculations, we identify key translocation-enhancing attributes of transmembrane proteins/peptides: a continuous hydrophilic patch, charged residues preferably in the membrane center, and aromatic hydrophobic residues. By employing both coarse-grained and atomistic simulations, complemented by experimental validation, we show that these properties not only enhance peptide translocation but also speed up lipid flip-flop. The enhanced flip-flop reinforces the idea that proteins such as scramblases and insertases not only share structural features but also operate through identical biophysical mechanisms enhancing the insertion and translocation of amphiphilic molecules. Our insights offer guidelines for the designing of translocation-enhancing proteins/peptides that could be used in medical and biotechnological applications.

Nanoparticle induced fusion of lipid membranes

Membrane fusion is crucial for infection of enveloped viruses, cellular transport, and drug delivery via liposomes. Nanoparticles can serve as fusogenic agents facilitating such membrane fusion for direct transmembrane transport. However, the underlying mechanisms of nanoparticle-induced fusion and the ideal properties of such nanoparticles remain largely unknown. Here, we used molecular dynamics simulations to investigate the efficacy of spheroidal nanoparticles with different size, prolateness, and ligand interaction strengths to enhance fusion between vesicles. By systematically varying nanoparticle properties, we identified how each parameter affects the fusion process and determined the optimal parameter range that promotes fusion. These findings provide valuable insights for the design and optimization of fusogenic nanoparticles with potential biotechnological and biomedical applications.

Valency of Ligand-Receptor Binding from Pair Potentials

Coarse grained molecular dynamics simulations have been crucial for investigating the dynamics of nanoparticle uptake by cell membranes via ligand-receptor interactions. These models have enabled researchers to evaluate the effects of nanoparticle size, shape, and ligand distribution on cellular uptake. However, when pair potentials are used to represent ligand-receptor interactions, the number of receptors interacting with one ligand, valency, may vary. We demonstrate that the curvature of a nanoparticle, strength of ligand-receptor interactions, and ligand or receptor concentration change the valency, ranging from 3.4 to 5.1 in this study. Such a change in valency can create inaccurate comparisons between nanoparticles or even result in the uptake of smaller nanoparticles than would be expected. To rectify this inconsistency, we propose the adoption of a model based on bond formation and use it to determine the extent to which previous studies may have been affected. This work recommends avoiding pair potentials for modeling ligand-receptor interactions to ensure methodological consistency in nanoparticle studies.

Insertases scramble lipids: Molecular simulations of MTCH2

Scramblases play a pivotal role in facilitating bidirectional lipid transport across cell membranes, thereby influencing lipid metabolism, membrane homeostasis, and cellular signaling. MTCH2, a mitochondrial outer membrane protein insertase, has a membrane-spanning hydrophilic groove resembling those that form the lipid transit pathway in known scramblases. Employing both coarse-grained and atomistic molecular dynamics simulations, we show that MTCH2 significantly reduces the free energy barrier for lipid movement along the groove and therefore can indeed function as a scramblase. Notably, the scrambling rate of MTCH2 in silico is similar to that of voltage-dependent anion channel (VDAC), a recently discovered scramblase of the outer mitochondrial membrane, suggesting a potential complementary physiological role for these mitochondrial proteins. Finally, our findings suggest that other insertases which possess a hydrophilic path across the membrane like MTCH2, can also function as scramblases.

Peptide translocation across asymmetric phospholipid membranes

The transport of molecules across cell membranes is vital for proper cell function and effective drug delivery. While most cell membranes naturally possess an asymmetric lipid composition, research on membrane transport predominantly uses symmetric lipid membranes. The permeation through the asymmetric membrane is then calculated as a sum of the inverse permeabilities of leaflets from symmetric bilayers. In this study, we examined how two types of amphiphilic molecules translocate across both asymmetric and symmetric membranes. Using computer simulations with both coarse-grained and atomistic force fields, we calculated the free energy profiles for the passage of model amphiphilic peptides and a lipid across various membranes. Our results consistently demonstrate that while the free energy profiles for asymmetric membranes with a small differential stress concur with symmetric ones in the region of lipid headgroups, the profiles differ around the center of the membrane. In this region, the free energy for the asymmetric membrane transitions between the profiles for two symmetric membranes. In addition, we show that peptide permeability through an asymmetric membrane cannot always be predicted from the permeabilities of the symmetric membranes. This indicates that using symmetric membranes falls short in providing an accurate depiction of peptide translocation across asymmetric membranes.

Amphipathic Helices Can Sense Both Positive and Negative Curvatures of Lipid Membranes

Curvature sensing is an essential ability of biomolecules to preferentially localize to membrane regions of a specific curvature. It has been shown that amphipathic helices (AHs), helical peptides with both hydrophilic and hydrophobic regions, could sense a positive membrane curvature. The origin of this AH sensing has been attributed to their ability to exploit lipid-packing defects that are enhanced in regions of positive curvature. In this study, we revisit an alternative framework where AHs act as sensors of local internal stress within the membrane, suggesting the possibility of an AH sensing a negative membrane curvature. Using molecular dynamics simulations, we gradually tuned the hydrophobicity of AHs, thereby adjusting their insertion depth so that the curvature preference of AHs is switched from positive to negative. This study suggests that highly hydrophobic AHs could preferentially localize proteins to regions of a negative membrane curvature.

Phospholipids are imported into mitochondria by VDAC, a dimeric beta barrel scramblase

Mitochondria are double-membrane-bounded organelles that depend critically on phospholipids supplied by the endoplasmic reticulum. These lipids must cross the outer membrane to support mitochondrial function, but how they do this is unclear. We identify the Voltage Dependent Anion Channel (VDAC), an abundant outer membrane protein, as a scramblase-type lipid transporter that catalyzes lipid entry. On reconstitution into membrane vesicles, dimers of human VDAC1 and VDAC2 catalyze rapid transbilayer translocation of phospholipids by a mechanism that is unrelated to their channel activity. Coarse-grained molecular dynamics simulations of VDAC1 reveal that lipid scrambling occurs at a specific dimer interface where polar residues induce large water defects and bilayer thinning. The rate of phospholipid import into yeast mitochondria is an order of magnitude lower in the absence of VDAC homologs, indicating that VDACs provide the main pathway for lipid entry. Thus, VDAC isoforms, members of a superfamily of beta barrel proteins, moonlight as a class of phospholipid scramblases - distinct from alpha-helical scramblase proteins - that act to import lipids into mitochondria.

Insertases Scramble Lipids: Molecular Simulations of MTCH2

Scramblases play a pivotal role in facilitating bidirectional lipid transport across cell membranes, thereby influencing lipid metabolism, membrane homeostasis, and cellular signaling. MTCH2, a mitochondrial outer membrane protein insertase, has a membrane-spanning hydrophilic groove resembling those that form the lipid transit pathway in known scramblases. Employing both coarse-grained and atomistic molecular dynamics simulations, we show that MTCH2 significantly reduces the free energy barrier for lipid movement along the groove and therefore can indeed function as a scramblase. Notably, the scrambling rate of MTCH2 in silico is similar to that of VDAC, a recently discovered scramblase of the outer mitochondrial membrane, suggesting a potential complementary physiological role for these mitochondrial proteins. Finally, our findings suggest that other insertases which possess a hydrophilic path across the membrane like MTCH2, can also function as scramblases.

Magainin 2 and PGLa in bacterial membrane mimics IV: Membrane curvature and partitioning

We previously reported that the synergistically enhanced antimicrobial activity of magainin 2 (MG2a) and PGLa is related to membrane adhesion and fusion. Here, we demonstrate that equimolar mixtures of MG2a and L18W-PGLa induce positive monolayer curvature stress and sense, at the same time, positive mean and Gaussian bilayer curvatures already at low amounts of bound peptide. The combination of both abilities-membrane curvature sensing and inducing-is most likely the base for the synergistically enhanced peptide activity. In addition, our coarse-grained simulations suggest that fusion stalks are promoted by decreasing the free-energy barrier for their formation rather than by stabilizing their shape. We also interrogated peptide partitioning as a function of lipid and peptide concentration using tryptophan fluorescence spectroscopy and peptide-induced leakage of dyes from lipid vesicles. In agreement with a previous report, we find increased membrane partitioning of L18W-PGLa in the presence of MG2a. However, this effect does not prevail to lipid concentrations higher than 1 mM, above which all peptides associate with the lipid bilayers. This implies that synergistic effects of MG2a and L18W-PGLa in previously reported experiments with lipid concentrations >1 mM are due to peptide-induced membrane remodeling and not their specific membrane partitioning.

Binding of DEP domain to phospholipid membranes: More than just electrostatics

Over the past decades an extensive effort has been made to provide a more comprehensive understanding of Wnt signaling, yet many regulatory and structural aspects remain elusive. Among these, the ability of Dishevelled (DVL) protein to relocalize at the plasma membrane is a crucial step in the activation of all Wnt pathways. The membrane binding of DVL was suggested to be mediated by the preferential interaction of its C-terminal DEP domain with phosphatidic acid (PA). However, due to the scarcity and fast turnover of PA, we investigated the role on the membrane association of other more abundant phospholipids. The combined results from computational simulations and experimental measurements with various model phospholipid membranes, demonstrate that the membrane binding of DEP/DVL constructs is governed by the concerted action of generic electrostatics and finely-tuned intermolecular interactions with individual lipid species. In particular, while we confirmed the strong preference for PA lipid, we also observed a weak but non-negligible affinity for phosphatidylserine, the most abundant anionic phospholipid in the plasma membrane, and phosphatidylinositol 4,5-bisphosphate. The obtained molecular insight into DEP-membrane interaction helps to elucidate the relation between changes in the local membrane composition and the spatiotemporal localization of DVL and, possibly, other DEP-containing proteins.

Synthesis and Profiling of Highly Selective Inhibitors of Methyltransferase DOT1L Based on Carbocyclic C-Nucleosides

Histone methyltransferase DOT1L is an attractive therapeutic target for the treatment of hematological malignancies. Here, we report the design, synthesis, and profiling of new DOT1L inhibitors based on nonroutine carbocyclic C-nucleoside scaffolds. The experimentally observed SAR was found to be nontrivial as seemingly minor changes of individual substituents resulted in significant changes in the affinity to DOT1L. Molecular modeling suggested that these trends could be related to significant conformational changes of the protein upon interaction with the inhibitors. The compounds 22 and (-)-53 (MU1656), carbocyclic C-nucleoside analogues of the natural nucleoside derivative EPZ004777, and the clinical candidate EPZ5676 (pinometostat) potently and selectively inhibit DOT1L in vitro as well as in the cell. The most potent compound MU1656 was found to be more metabolically stable and significantly less toxic in vivo than pinometostat itself.

Magainin 2 and PGLa in bacterial membrane mimics III: Membrane fusion and disruption

We previously speculated that the synergistically enhanced antimicrobial activity of Magainin 2 and PGLa is related to membrane adhesion, fusion, and further membrane remodeling. Here we combined computer simulations with time-resolved in vitro fluorescence microscopy, cryoelectron microscopy, and small-angle X-ray scattering to interrogate such morphological and topological changes of vesicles at nanoscopic and microscopic length scales in real time. Coarse-grained simulations revealed formation of an elongated and bent fusion zone between vesicles in the presence of equimolar peptide mixtures. Vesicle adhesion and fusion were observed to occur within a few seconds by cryoelectron microscopy and corroborated by small-angle X-ray scattering measurements. The latter experiments indicated continued and time-extended structural remodeling for individual peptides or chemically linked peptide heterodimers but with different kinetics. Fluorescence microscopy further captured peptide-dependent adhesion, fusion, and occasional bursting of giant unilamellar vesicles a few seconds after peptide addition. The synergistic interactions between the peptides shorten the time response of vesicles and enhance membrane fusogenic and disruption properties of the equimolar mixture compared with the individual peptides.

Cargo Release from Nonenveloped Viruses and Virus-like Nanoparticles: Capsid Rupture or Pore Formation

Virus-like nanoparticles are protein shells similar to wild-type viruses, and both aim to deliver their content into a cell. Unfortunately, the release mechanism of their cargo/genome remains elusive. Pores on the symmetry axes were proposed to enable the slow release of the viral genome. In contrast, cryo-EM images showed that capsids of nonenveloped RNA viruses can crack open and rapidly release the genome. We combined in vitro cryo-EM observations of the genome release of three viruses with coarse-grained simulations of generic virus-like nanoparticles to investigate the cargo/genome release pathways. Simulations provided details on both slow and rapid release pathways, including the success rates of individual releases. Moreover, the simulated structures from the rapid release pathway were in agreement with the experiment. Slow release occurred when interactions between capsid subunits were long-ranged, and the cargo/genome was noncompact. In contrast, rapid release was preferred when the interaction range was short and/or the cargo/genome was compact. These findings indicate a design strategy of virus-like nanoparticles for drug delivery.

The impact of the glycan headgroup on the nanoscopic segregation of gangliosides

Gangliosides form an important class of receptor lipids containing a large oligosaccharide headgroup whose ability to self-organize within lipid membranes results in the formation of nanoscopic platforms. Despite their biological importance, the molecular basis for the nanoscopic segregation of gangliosides is not clear. In this work, we investigated the role of the ganglioside headgroup on the nanoscale organization of gangliosides. We studied the effect of the reduction in the number of sugar units of the ganglioside oligosaccharide chain on the ability of gangliosides GM1, GM2, and GM3 to spontaneously self-organize into lipid nanodomains. To reach nanoscopic resolution and to identify molecular forces that drive ganglioside segregation, we combined an experimental technique, Förster resonance energy transfer analyzed by Monte-Carlo simulations offering high lateral and trans-bilayer resolution with molecular dynamics simulations. We show that the ganglioside headgroup plays a key role in ganglioside self-assembly despite the negative charge of the sialic acid group. The nanodomains range from 7 to 120 nm in radius and are mostly composed of the surrounding bulk lipids, with gangliosides being a minor component of the nanodomains. The interactions between gangliosides are dominated by the hydrogen bonding network between the headgroups, which facilitates ganglioside clustering. The N-acetylgalactosamine sugar moiety of GM2, however, seems to impair the stability of these clusters by disrupting hydrogen bonding of neighboring sugars, which is in agreement with a broad size distribution of GM2 nanodomains. The simulations suggest that the formation of nanodomains is likely accompanied by several conformational changes in the gangliosides, which, however, have little impact on the solvent exposure of these receptor groups. Overall, this work identifies the key physicochemical factors that drive nanoscopic segregation of gangliosides.

Advances in Molecular Understanding of α-Helical Membrane-Active Peptides

Biological membranes separate the interior of cells or cellular compartments from their outer environments. This barrier function of membranes can be disrupted by membrane-active peptides, some of which can spontaneously penetrate through the membranes or open leaky transmembrane pores. However, the origin of their activity/toxicity is not sufficiently understood for the development of more potent peptides. To this day, there are no design rules that would be generally valid, and the role of individual amino acids tends to be sequence-specific.In this Account, we describe recent progress in understanding the design principles that govern the activity of membrane-active peptides. We focus on α-helical amphiphilic peptides and their ability to (1) translocate across phospholipid bilayers, (2) form transmembrane pores, or (3) act synergistically, i.e., to produce a significantly more potent effect in a mixture than the individual components.We refined the description of peptide translocation using computer simulations and demonstrated the effect of selected residues. Our simulations showed the necessity to explicitly include charged residues in the translocation description to correctly sample the membrane perturbations they can cause. Using this description, we calculated the translocation of helical peptides with and without the kink induced by the proline/glycine residue. The presence of the kink had no effect on the translocation barrier, but it decreased the peptide affinity to the membrane and reduced the peptide stability inside the membrane. Interestingly, the effects were mainly caused by the peptide's increased polarity, not the higher flexibility of the kink.Flexibility plays a crucial role in pore formation and affects distinct pore structures in different ways. The presence of a kink destabilizes barrel-stave pores, because the kink prevents the tight packing of peptides in the bundle, which is characteristic of the barrel-stave structure. In contrast, the kink facilitates the formation of toroidal pores, where the peptides are only loosely arranged and do not need to closely assemble. The exact position of the kink in the sequence further determines the preferred arrangement of peptides in the pore, i.e., an hourglass or U-shaped structure. In addition, we demonstrated that two self-associated (via termini) helical peptides could mimic the behavior of peptides with a helix-kink-helix motif.Finally, we review the recent findings on the peptide synergism of the archetypal mixture of Magainin 2 and PGLa peptides. We focused on a bacterial plasma membrane mimic that contains negatively charged lipids and lipids with negative intrinsic curvature. We showed that the synergistic action of peptides was highly dependent on the lipid composition. When the lipid composition and peptide/lipid ratios were changed, the systems exhibited more complex behavior than just the previously reported pore formation. We observed membrane adhesion, fusion, and even the formation of the sponge phase in this regime. Furthermore, enhanced adhesion/partitioning to the membrane was reported to be caused by lipid-induced peptide aggregation.In conclusion, the provided molecular insight into the complex behavior of membrane-active peptides provides clues for the design and modification of antimicrobial peptides or toxins.

Enhanced translocation of amphiphilic peptides across membranes by transmembrane proteins

Cell membranes are phospholipid bilayers with a large number of embedded transmembrane proteins. Some of these proteins, such as scramblases, have properties that facilitate lipid flip-flop from one membrane leaflet to another. Scramblases and similar transmembrane proteins could also affect the translocation of other amphiphilic molecules, including cell-penetrating or antimicrobial peptides. We studied the effect of transmembrane proteins on the translocation of amphiphilic peptides through the membrane. Using two very different models, we consistently demonstrate that transmembrane proteins with a hydrophilic patch enhance the translocation of amphiphilic peptides by stabilizing the peptide in the membrane. Moreover, there is an optimum amphiphilicity because the peptide could become overstabilized in the transmembrane state, in which the peptide-protein dissociation is hampered, limiting the peptide translocation. The presence of scramblases and other proteins with similar properties could be exploited for more efficient transport into cells. The described principles could also be utilized in the design of a drug-delivery system by the addition of a translocation-enhancing peptide that would integrate into the membrane.

Selecting Collective Variables and Free-Energy Methods for Peptide Translocation across Membranes

The selective permeability of cellular membranes is a crucial property for controlled transport into and out of cells. Molecules that can bypass the cellular machinery and spontaneously translocate across membranes could be used as therapeutics or drug carriers. Peptides are a prominent class of such molecules, which include natural and man-developed antimicrobial and cell-penetrating peptides. However, the necessary peptide properties for translocation remain elusive. Computer simulations could uncover these properties once we have a good collective variable (CV) that accurately describes the translocation process. Here, we developed a new CV, which includes a description of peptide insertion, local membrane deformation, and peptide internal degrees of freedom related to its charged groups. By comparison of CVs, we demonstrated that all these components are necessary for an accurate description of peptide translocation. Moreover, the advantages and disadvantages of three common methods for free-energy calculations with our CV were evaluated using the MARTINI coarse-grained model: umbrella sampling, umbrella sampling with replica exchange, and metadynamics. The developed CV leads to the reliable and effective calculation of the free energy of peptide translocation, and thus, it could be useful in the design of spontaneously translocating peptides.

Phosphorylation-induced changes in the PDZ domain of Dishevelled 3

The PDZ domain of Dishevelled 3 protein belongs to a highly abundant protein recognition motif which typically binds short C-terminal peptides. The affinity of the PDZ towards the peptides could be fine-tuned by a variety of post-translation modifications including phosphorylation. However, how phosphorylations affect the PDZ structure and its interactions with ligands remains elusive. Combining molecular dynamics simulations, NMR titration, and biological experiments, we explored the role of previously reported phosphorylation sites and their mimetics in the Dishevelled PDZ domain. Our observations suggest three major roles for phosphorylations: (1) acting as an on/off PDZ binding switch, (2) allosterically affecting the binding groove, and (3) influencing the secondary binding site. Our simulations indicated that mimetics had similar but weaker effects, and the effects of distinct sites were non-additive. This study provides insight into the Dishevelled regulation by PDZ phosphorylation. Furthermore, the observed effects could be used to elucidate the regulation mechanisms in other PDZ domains.

Capsid opening enables genome release of iflaviruses

The family Iflaviridae includes economically important viruses of the western honeybee such as deformed wing virus, slow bee paralysis virus, and sacbrood virus. Iflaviruses have nonenveloped virions and capsids organized with icosahedral symmetry. The genome release of iflaviruses can be induced in vitro by exposure to acidic pH, implying that they enter cells by endocytosis. Genome release intermediates of iflaviruses have not been structurally characterized. Here, we show that conformational changes and expansion of iflavirus RNA genomes, which are induced by acidic pH, trigger the opening of iflavirus particles. Capsids of slow bee paralysis virus and sacbrood virus crack into pieces. In contrast, capsids of deformed wing virus are more flexible and open like flowers to release their genomes. The large openings in iflavirus particles enable the fast exit of genomes from capsids, which decreases the probability of genome degradation by the RNases present in endosomes.

Effect of Helical Kink on Peptide Translocation across Phospholipid Membranes

Biological membranes present a major obstacle for the delivery of therapeutic agents into cells. Some peptides have been shown to translocate across the membrane spontaneously, and they could be thus used as drug-carriers. However, the advantageous peptide properties for the translocation remain unclear. Of particular interest is the effect of a proline-induced kink in α-helical peptides, because the kink was previously reported to both increase and decrease the antimicrobial activity. The antimicrobial activity of peptides could be related to their translocation across the membrane as is the case of the buforin 2 peptide investigated here. Using computer simulations with two independent models, we consistently showed that the presence of the kink has (1) no effect on the translocation barrier, (2) reduces the peptide affinity to the membrane, and (3) disfavors the transmembrane state. Moreover, we were able to determine that these effects are mainly caused by the peptide increased polarity, not the increased flexibility of the kink. The provided molecular understanding can be utilized for the design of cell-penetrating and drug-carrying peptides.

Effect of helical kink in antimicrobial peptides on membrane pore formation

Every cell is protected by a semipermeable membrane. Peptides with the right properties, for example Antimicrobial peptides (AMPs), can disrupt this protective barrier by formation of leaky pores. Unfortunately, matching peptide properties with their ability to selectively form pores in bacterial membranes remains elusive. In particular, the proline/glycine kink in helical peptides was reported to both increase and decrease antimicrobial activity. We used computer simulations and fluorescence experiments to show that a kink in helices affects the formation of membrane pores by stabilizing toroidal pores but disrupting barrel-stave pores. The position of the proline/glycine kink in the sequence further controls the specific structure of toroidal pore. Moreover, we demonstrate that two helical peptides can form a kink-like connection with similar behavior as one long helical peptide with a kink. The provided molecular-level insight can be utilized for design and modification of pore-forming antibacterial peptides or toxins.

Magainin 2 and PGLa in Bacterial Membrane Mimics II: Membrane Fusion and Sponge Phase Formation

We studied the synergistic mechanism of equimolar mixtures of magainin 2 (MG2a) and PGLa in phosphatidylethanolamine/phosphatidylglycerol mimics of Gram-negative cytoplasmic membranes. In a preceding article of this series, we reported on the early onset of parallel heterodimer formation of the two antimicrobial peptides already at low concentrations and the resulting defect formation in the membranes. Here, we focus on the structures of the peptide-lipid aggregates occurring in the synergistic regime at elevated peptide concentrations. Using a combination of calorimetric, scattering, electron microscopic, and in silico techniques, we demonstrate that the two peptides, even if applied individually, transform originally large unilamellar vesicles into multilamellar vesicles with a collapsed interbilayer spacing resulting from peptide-induced adhesion. Interestingly, the adhesion does not lead to a peptide-induced lipid separation of charged and charge-neutral species. In addition to this behavior, equimolar mixtures of MG2a and PGLa formed surface-aligned fibril-like structures, which induced adhesion zones between the membranes and the formation of transient fusion stalks in molecular dynamics simulations and a coexisting sponge phase observed by small-angle x-ray scattering. The previously reported increased leakage of lipid vesicles of identical composition in the presence of MG2a/PGLa mixtures is therefore related to a peptide-induced cross-linking of bilayers.

Effect of membrane composition on DivIVA-membrane interaction

DivIVA is a crucial membrane-binding protein that helps to localize other proteins to negatively curved membranes at cellular poles and division septa in Gram-positive bacteria. The N-terminal domain of DivIVA is responsible for membrane binding. However, to which lipids the domain binds or how it recognizes the membrane negative curvature remains elusive. Using computer simulations, we demonstrate that the N-terminal domain of Streptomyces coelicolor DivIVA adsorbs to membranes with affinity and orientation dependent on the lipid composition. The domain interacts non-specifically with lipid phosphates via its arginine-rich tip and the strongest interaction is with cardiolipin. Moreover, we observed a specific attraction between a negatively charged side patch of the domain and ethanolamine lipids, which addition caused the change of the domain orientation from perpendicular to parallel alignment to the membrane plane. Similar but less electrostatically dependent behavior was observed for the N-terminal domain of Bacillus subtilis. The domain propensity for lipids which prefer negatively curved membranes could be a mechanism for the cellular localization of DivIVA protein.

Magainin 2 and PGLa in Bacterial Membrane Mimics I: Peptide-Peptide and Lipid-Peptide Interactions

We addressed the onset of synergistic activity of the two well-studied antimicrobial peptides magainin 2 (MG2a) and PGLa using lipid-only mimics of Gram-negative cytoplasmic membranes. Specifically, we coupled a joint analysis of small-angle x-ray and neutron scattering experiments on fully hydrated lipid vesicles in the presence of MG2a and L18W-PGLa to all-atom and coarse-grained molecular dynamics simulations. In agreement with previous studies, both peptides, as well as their equimolar mixture, were found to remain upon adsorption in a surface-aligned topology and to induce significant membrane perturbation, as evidenced by membrane thinning and hydrocarbon order parameter changes in the vicinity of the inserted peptide. These effects were particularly pronounced for the so-called synergistic mixture of 1:1 (mol/mol) L18W-PGLa/MG2a and cannot be accounted for by a linear combination of the membrane perturbations of two peptides individually. Our data are consistent with the formation of parallel heterodimers at concentrations below a synergistic increase of dye leakage from vesicles. Our simulations further show that the heterodimers interact via salt bridges and hydrophobic forces, which apparently makes them more stable than putatively formed antiparallel L18W-PGLa and MG2a homodimers. Moreover, dimerization of L18W-PGLa and MG2a leads to a relocation of the peptides within the lipid headgroup region as compared to the individual peptides. The early onset of dimerization of L18W-PGLa and MG2a at low peptide concentrations consequently appears to be key to their synergistic dye-releasing activity from lipid vesicles at high concentrations.

Force Field Comparison of GM1 in a DOPC Bilayer Validated with AFM and FRET Experiments

The great physiological relevance of glycolipids is being increasingly recognized, and glycolipid interactions have been shown to be central to cell-cell recognition, neuronal plasticity, protein-ligand recognition, and other important processes. However, detailed molecular-level understanding of these processes remains to be fully resolved. Molecular dynamics simulations could reveal the details of the glycolipid interactions, but the results may be influenced by the choice of the employed force field. Here, we have compared the behavior and properties of GM1, a common, biologically important glycolipid, using the CHARMM36, OPLS, GROMOS, and Amber99-GLYCAM06 (in bilayers comprising SLIPIDS and LIPID14 lipids) force fields in bilayers comprising 1,2-dioleoyl-sn-glycero-3-phosphocholine lipids and compared the results to atomic force microscopy and fluorescence resonance energy transfer experiments. We found discrepancies within the GM1 behavior displayed between the investigated force fields. Based on a direct comparison with complementary experimental results derived from fluorescence and AFM measurements, we recommend using the Amber99-GLYCAM force field in bilayers comprising LIPID14 or SLIPIDS lipids followed by CHARMM36 and OPLS force fields in simulations. The GROMOS force field is not recommended for reproducing the properties of the GM1 head group.

Synergism of Antimicrobial Frog Peptides Couples to Membrane Intrinsic Curvature Strain

Mixtures of the frog peptides magainin 2 and PGLa are well-known for their pronounced synergistic killing of Gram-negative bacteria. We aimed to gain insight into the underlying biophysical mechanism by interrogating the permeabilizing efficacies of the peptides as a function of stored membrane curvature strain. For Gram-negative bacterial-inner-membrane mimics, synergism was only observed when the anionic bilayers exhibited significant negative intrinsic curvatures imposed by monounsaturated phosphatidylethanolamine. In contrast, the peptides and their mixtures did not exhibit significant activities in charge-neutral mammalian mimics, including those with negative curvature, which is consistent with the requirement of charge-mediated peptide binding to the membrane. Our experimental findings are supported by computer simulations showing a significant decrease of the peptide-insertion free energy in membranes upon shifting intrinsic curvatures toward more positive values. The physiological relevance of our model studies is corroborated by a remarkable agreement with the peptide’s synergistic activity in Escherichia coli. We propose that synergism is related to a lowering of a membrane-curvature-strain-mediated free-energy barrier by PGLa that assists membrane insertion of magainin 2, and not by strict pairwise interactions of the two peptides as suggested previously.

Dishevelled-3 conformation dynamics analyzed by FRET-based biosensors reveals a key role of casein kinase 1

Dishevelled (DVL) is the key component of the Wnt signaling pathway. Currently, DVL conformational dynamics under native conditions is unknown. To overcome this limitation, we develop the Fluorescein Arsenical Hairpin Binder- (FlAsH-) based FRET in vivo approach to study DVL conformation in living cells. Using this single-cell FRET approach, we demonstrate that (i) Wnt ligands induce open DVL conformation, (ii) DVL variants that are predominantly open, show more even subcellular localization and more efficient membrane recruitment by Frizzled (FZD) and (iii) Casein kinase 1 ɛ (CK1ɛ) has a key regulatory function in DVL conformational dynamics. In silico modeling and in vitro biophysical methods explain how CK1ɛ-specific phosphorylation events control DVL conformations via modulation of the PDZ domain and its interaction with DVL C-terminus. In summary, our study describes an experimental tool for DVL conformational sampling in living cells and elucidates the essential regulatory role of CK1ɛ in DVL conformational dynamics.

Enterovirus particles expel capsid pentamers to enable genome release

Viruses from the genus Enterovirus are important human pathogens. Receptor binding or exposure to acidic pH in endosomes converts enterovirus particles to an activated state that is required for genome release. However, the mechanism of enterovirus uncoating is not well understood. Here, we use cryo-electron microscopy to visualize virions of human echovirus 18 in the process of genome release. We discover that the exit of the RNA from the particle of echovirus 18 results in a loss of one, two, or three adjacent capsid-protein pentamers. The opening in the capsid, which is more than 120 Å in diameter, enables the release of the genome without the need to unwind its putative double-stranded RNA segments. We also detect capsids lacking pentamers during genome release from echovirus 30. Thus, our findings uncover a mechanism of enterovirus genome release that could become target for antiviral drugs.

Genome release of enteroviruses relies on exposure to acidic pH, but the mechanism of uncoating remains unclear. Here, Buchta et al. show that echovirus 18 loses one to three adjacent capsid-protein pentamers, resulting in an opening of more than 120 Å for genome release.

Optimal Hydrophobicity and Reorientation of Amphiphilic Peptides Translocating through Membrane

Cell-penetrating and some antimicrobial peptides can translocate across lipid bilayers without disrupting the membrane structure. However, the molecular properties required for efficient translocation are not fully understood. We employed the Metropolis Monte Carlo method together with coarse-grained models to systematically investigate free-energy landscapes associated with the translocation of secondary amphiphilic peptides. We studied α-helical peptides with different length, amphiphilicity, and distribution of hydrophobic content and found a common translocation path consisting of adsorption, tilting, and insertion. In the adsorbed state, the peptides are parallel to the membrane plane, whereas, in the inserted state, the peptides are perpendicular to the membrane. Our simulations demonstrate that, for all tested peptides, there is an optimal ratio of hydrophilic/hydrophobic content at which the peptides cross the membrane the easiest. Moreover, we show that the hydrophobicity of peptide termini has an important effect on the translocation barrier. These results provide general guidance to optimize peptides for use as carriers of molecular cargos or as therapeutics themselves.

Aggregate Size Dependence of Amyloid Adsorption onto Charged Interfaces

Amyloid aggregates are associated with a range of human neurodegenerative disorders, and it has been shown that neurotoxicity is dependent on aggregate size. Combining molecular simulation with analytical theory, a predictive model is proposed for the adsorption of amyloid aggregates onto oppositely charged surfaces, where the interaction is governed by an interplay between electrostatic attraction and entropic repulsion. Predictions are experimentally validated against quartz crystal microbalance-dissipation experiments of amyloid beta peptides and fragmented fibrils in the presence of a supported lipid bilayer. Assuming amyloids as rigid, elongated particles, we observe nonmonotonic trends for the extent of adsorption with respect to aggregate size and preferential adsorption of smaller aggregates over larger ones. Our findings describe a general phenomenon with implications for stiff polyions and rodlike particles that are electrostatically attracted to a surface.

Self-assembled clusters of patchy rod-like molecules

The design of complex self-assembled structures remains a challenging task due to the intricate relationship between the properties of the building blocks and the final morphology of the aggregates. Here, we report such a relationship for rod-like particles with one or two attractive patches based on a combination of computer simulations and analytical theory. We investigated the formation of finite aggregates under various conditions and constructed structure diagrams, which can be used to determine and extrapolate the system composition. The size of the clusters is mainly determined by the size of the attractive patches and their geometrical arrangement. We showed that it is challenging to obtain clusters with more than four particles in high yields, and more complex building blocks or additional molecules would need to be used. Moreover, the particles with patch sizes close to the structure boundaries can switch between the aggregate state by a small change in conditions. These findings can be useful for the development of self-assembling building blocks and for the understanding of protein folds of coiled coils under various conditions.

Design of Multivalent Inhibitors for Preventing Cellular Uptake

Cellular entry, the first crucial step of viral infection, can be inhibited by molecules adsorbed on the virus surface. However, apart from using stronger affinity, little is known about the properties of such inhibitors that could increase their effectiveness. Our simulations showed that multivalent inhibitors can be designed to be much more efficient than their monovalent counterparts. For example, for our particular simulation model, a single multivalent inhibitor spanning 5 to 6 binding sites is enough to prevent the uptake compared to the required 1/3 of all the receptor binding sites needed to be blocked by monovalent inhibitors. Interestingly, multivalent inhibitors are more efficient in inhibiting the uptake not only due to their increased affinity but mainly due to the co-localization of the inhibited receptor binding sites at the virion's surface. Furthermore, we show that Janus-like inhibitors do not induce virus aggregation. Our findings may be generalized to other uptake processes including bacteria and drug-delivery.

GM1 Ganglioside Inhibits β-Amyloid Oligomerization Induced by Sphingomyelin

β-Amyloid (Aβ) oligomers are neurotoxic and implicated in Alzheimer's disease. Neuronal plasma membranes may mediate formation of Aβ oligomers in vivo. Membrane components sphingomyelin and GM1 have been shown to promote aggregation of Aβ; however, these studies were performed under extreme, non-physiological conditions. We demonstrate that physiological levels of GM1 , organized in nanodomains do not seed oligomerization of Aβ40 monomers. We show that sphingomyelin triggers oligomerization of Aβ40 and that GM1 is counteractive thus preventing oligomerization. We propose a molecular explanation that is supported by all-atom molecular dynamics simulations. The preventive role of GM1 in the oligomerization of Aβ40 suggests that decreasing levels of GM1 in the brain, for example, due to aging, could reduce protection against Aβ oligomerization and contribute to the onset of Alzheimer's disease.

Influence of ligand distribution on uptake efficiency

Cellular uptake is a crucial process in nanomedicine and drug-delivery; however, the factors that affect its efficiency/speed are not well understood. We report computer simulations on passive uptake via receptor-mediated endocytosis of nanoparticle coated with ligands. In particular, we study how the distribution of ligands on the nanoparticle surface influences the uptake rate. The speed of membrane wrapping and uptake was found to be the fastest for nanoparticles with homogeneous ligand distributions, where ligands are spread evenly on the surface. We show that the diffusion of the ligands on the nanoparticle can hinder its uptake, since upon the interaction with the membrane the ligand distribution becomes extremely inhomogeneous, with a large ligand-free patch. When the ligand-free-area was more than 20% of the surface, we did not observe uptake within the scale of our simulations.

Stability of bicelles: a simulation study

Aqueous mixtures of long-tailed lipids (e.g., dimyristoylphosphatidylcholine - DMPC) and detergents can sometimes form membrane disks called bicelles. Bicelles have found applications as an embedding medium for membrane proteins in the context of NMR studies and protein crystallization. However, the parameters that determine the thermodynamic stability of bicelles are not well understood. Here we report a coarse-grained simulation study of the relationship between lipid-aggregate morphology and the composition and temperature of the surfactant mixture. In agreement with experiments, we find that bicellar mixtures are destabilized at higher temperatures and detergents are present at membrane edges as well as in flat membranes with a strong preference for the edges. In addition, our results suggest that the free-energy difference between bicelles and the perforated lamellar phase is typically very small for molecules without intrinsic curvature and charge. Cone shaped surfactant molecules tend to favor the formation of bicelles; however, none of the systems that we have studied provide unambiguous evidence for the existence of thermodynamically stable bicelles in mixtures of uncharged lipids with long and short tails. We speculate that small changes in the properties of the system (charge, dopants) may make bicelles thermodynamically stable.

Simulations suggest possible novel membrane pore structure

Amphiphilic proteins and peptides can induce the formation of stable and metastable pores in membranes. Using coarse-grained simulations, we have studied the factors that affect structure of peptide-stabilized pores. Our simulations are able to reproduce the formation of the well-known barrel-stave or toroidal pores, but in addition, we find evidence for a novel "double-belt" pore structure: in this structure the peptides that coat the membrane pore are oriented parallel to the membrane plane. To check the predictions of our coarse-grained model, we have performed more detailed simulations, using the MARTINI force field. These simulations show that the double-belt structure is stable up to at least the microsecond time scale.

Intracellular release of endocytosed nanoparticles upon a change of ligand-receptor interaction

During passive endocytosis, nanosized particles are initially encapsulated by a membrane separating it from the cytosol. Yet, in many applications the nanoparticles need to be in direct contact with the cytosol in order to be active. We report a simulation study that elucidates the physical mechanisms by which such nanoparticles can shed their bilayer coating. We find that nanoparticle release can be readily achieved by a pH-induced lowering of the attraction between nanoparticle and membrane only if the nanoparticle is either very small or nonspherical. Interestingly, we find that in the case of large spherical nanoparticles, the reduction of attraction needs to be accompanied by exerting an additional tension on the membrane (e.g., via nanoparticle expansion) to achieve release. We expect these findings will contribute to the rational design of drug delivery strategies via nanoparticles.

Connecting macroscopic observables and microscopic assembly events in amyloid formation using coarse grained simulations

The pre-fibrillar stages of amyloid formation have been implicated in cellular toxicity, but have proved to be challenging to study directly in experiments and simulations. Rational strategies to suppress the formation of toxic amyloid oligomers require a better understanding of the mechanisms by which they are generated. We report Dynamical Monte Carlo simulations that allow us to study the early stages of amyloid formation. We use a generic, coarse-grained model of an amyloidogenic peptide that has two internal states: the first one representing the soluble random coil structure and the second one the -sheet conformation. We find that this system exhibits a propensity towards fibrillar self-assembly following the formation of a critical nucleus. Our calculations establish connections between the early nucleation events and the kinetic information available in the later stages of the aggregation process that are commonly probed in experiments. We analyze the kinetic behaviour in our simulations within the framework of the theory of classical nucleated polymerisation, and are able to connect the structural events at the early stages in amyloid growth with the resulting macroscopic observables such as the effective nucleus size. Furthermore, the free-energy landscapes that emerge from these simulations allow us to identify pertinent properties of the monomeric state that could be targeted to suppress oligomer formation.

A number of normally soluble proteins can form amyloid structures in a process associated with neurodegenerative diseases such as Alzheimer's and Parkinson's diseases. Mature amyloid structures consist of large fibrils containing thousands of individual proteins aggregated into linear nanostructures; there is increasing evidence, however, that the toxic species responsible for neurodegeneration are not the mature fibrils themselves but rather lower molecular weight precursors commonly known as amyloid oligomers. Unfortunately, these early oligomers are commonly thermodynamically unstable and of nanometer scale dimensions, factors which make them highly challenging to probe in detail in experiments. We have used computer simulations of a model inspired by Alzheimer's Abeta peptide to investigate the early stages of protein aggregation. The results that we obtain were shown to fit Oosawa's polymerization theory, a finding which allows us to provide a connection between the microscopic molecular parameters and macroscopic growth. One crucial parameter is size of the nucleus, i.e. the basic oligomer existing at origin of the formation of each fiber. We have revealed a path for the formation of this nucleus and validate its size by several methods. Our results provide fundamental information for influencing the early stages of amyloid formation in a rational manner.

Aqueous solutions at the interface with phospholipid bilayers

In a sense, life is defined by membranes, because they delineate the barrier between the living cell and its surroundings. Membranes are also essential for regulating the machinery of life throughout many interfaces within the cell's interior. A large number of experimental, computational, and theoretical studies have demonstrated how the properties of water and ionic aqueous solutions change due to the vicinity of membranes and, in turn, how the properties of membranes depend on the presence of aqueous solutions. Consequently, understanding the character of aqueous solutions at their interface with biological membranes is critical to research progress on many fronts. The importance of incorporating a molecular-level description of water into the study of biomembrane surfaces was demonstrated by an examination of the interaction between phospholipid bilayers that can serve as model biological membranes. The results showed that, in addition to well-known forces, such as van der Waals and screened Coulomb, one has to consider a repulsion force due to the removal of water between surfaces. It was also known that physicochemical properties of biological membranes are strongly influenced by the specific character of the ions in the surrounding aqueous solutions because of the observation that different anions produce different effects on muscle twitch tension. In this Account, we describe the interaction of pure water, and also of aqueous ionic solutions, with model membranes. We show that a symbiosis of experimental and computational work over the past few years has resulted in substantial progress in the field. We now better understand the origin of the hydration force, the structural properties of water at the interface with phospholipid bilayers, and the influence of phospholipid headgroups on the dynamics of water. We also improved our knowledge of the ion-specific effect, which is observed at the interface of the phospholipid bilayer and aqueous solution, and its connection with the Hofmeister series. Nevertheless, despite substantial progress, many issues remain unresolved. Thus, for example, we still cannot satisfactorily explain the force of interaction between phospholipid bilayers immersed in aqueous solutions of NaI. Although we try to address many issues here, the scope of the discussion is limited and does not cover such important topics as the influence of ionic solutions on phases of bilayers, the influence of salts on the properties of Langmuir monolayers containing lipid molecules, or the influence of aqueous solutions on bilayers containing mixtures of lipids. We anticipate that the future application of more powerful experimental techniques, in combination with more advanced computational hardware, software, and theory, will produce molecular-level information about these important topics and, more broadly, will further illuminate our understanding of interfaces between aqueous solutions and biological membranes.

Receptor-mediated endocytosis of nanoparticles of various shapes

Cellular uptake through endocytosis is crucial for drug delivery and nanomedicine. However, the conditions under which passive endocytosis (i.e., not ATP driven) takes place are not well understood. We report MD simulations of the passive uptake of ligand-coated nanoparticles with varying size, shape, coverage, and membrane-binding strength. We find that the efficiency of passive endocytosis is higher for spherocylindrical particles than for spheres and that endocytosis is suppressed for particles with sharp edges.