Dynamic Second Harmonic Imaging of Proton Translocation Through Water Needles in Lipid Membranes.

Proton translocation through lipid membranes is a fundamental process in the field of biology. Several theoretical models have been developed and presented over the years to explain the phenomenon, yet the exact mechanism is still not well understood. Here, we show that proton translocation is directly related to membrane potential fluctuations. Using high-throughput wide-field second harmonic (SH) microscopy, we report apparently universal transmembrane potential fluctuations in lipid membrane systems. Molecular simulations and free energy calculations suggest that H+ permeation proceeds predominantly across a thin, membrane-spanning water needle and that the transient transmembrane potential drives H+ ions across the water needle. This mechanism differs from the transport of other cations that require completely open pores for transport and follows naturally from the well-known Grotthuss mechanism for proton transport in bulk water. Furthermore, SH imaging and conductivity measurements reveal that the rate of proton transport depends on the structure of the hydrophobic core of bilayer membranes.

Passive transport of Ca2+ ions through lipid bilayers imaged by widefield second harmonic microscopy.

In biology, release of Ca2+ ions in the cytosol is essential to trigger or control many cell functions. Calcium signaling acutely depends on lipid membrane permeability to Ca2+. For proper understanding of membrane permeability to Ca2+, both membrane hydration and the structure of the hydrophobic core must be taken into account. Here, we vary the hydrophobic core of bilayer membranes and observe different types of behavior in high-throughput wide-field second harmonic imaging. Ca2+ translocation is observed through mono-unsaturated (DOPC:DOPA) membranes, reduced upon the addition of cholesterol, and completely inhibited for branched (DPhPC:DPhPA) and poly-unsaturated (SLPC:SLPA) lipid membranes. We propose, using molecular dynamics simulations, that ion transport occurs through ion-induced transient pores, which requires nonequilibrium membrane restructuring. This results in different rates at different locations and suggests that the hydrophobic structure of lipids plays a much more sophisticated regulating role than previously thought.

Ion-Induced Transient Potential Fluctuations Facilitate Pore Formation and Cation Transport through Lipid Membranes.

Unassisted ion transport through lipid membranes plays a crucial role in many cell functions without which life would not be possible, yet the precise mechanism behind the process remains unknown due to its molecular complexity. Here, we demonstrate a direct link between membrane potential fluctuations and divalent ion transport. High-throughput wide-field non-resonant second harmonic (SH) microscopy of membrane water shows that membrane potential fluctuations are universally found in lipid bilayer systems. Molecular dynamics simulations reveal that such variations in membrane potential reduce the free energy cost of transient pore formation and increase the ion flux across an open pore. These transient pores can act as conduits for ion transport, which we SH image for a series of divalent cations (Cu2+, Ca2+, Ba2+, Mg2+) passing through giant unilamellar vesicle (GUV) membranes. Combining the experimental and computational results, we show that permeation through pores formed via an ion-induced electrostatic field is a viable mechanism for unassisted ion transport.

Lipid Droplets Embedded in a Model Cell Membrane Create a Phospholipid Diffusion Barrier.

Lipid droplets (LDs) are ubiquitous, cytoplasmic fat storage organelles that originate from the endoplasmic reticulum (ER) membrane. They are composed of a core of neutral lipids surrounded by a phospholipid monolayer. Proteins embedded into this monolayer membrane adopt a monotopic topology and are crucial for regulated lipid storage and consumption. A key question is, which collective properties of protein-intrinsic and lipid-mediated features determine spatio-temporal protein partitioning between phospholipid bilayer and LD monolayer membranes. To address this question, a freestanding phospholipid bilayer with physiological lipidic composition is produced using microfluidics and micrometer-sized LDs are dispersed around the bilayer that spontaneously insert into the bilayer. Using confocal microscopy, the 3D geometry of the reconstituted LDs is determined with high spatial resolution. The micrometer-sized bilayer-embedded LDs present a characteristic lens shape that obeys predictions from equilibrium wetting theory. Fluorescence recovery after photobleaching measurements reveals the existence of a phospholipid diffusion barrier at the monolayer-bilayer interface. Coarse-grained molecular dynamics simulation reveals lipid specific density distributions along the pore rim, which may rationalize the diffusion barrier. The lipid diffusion barrier between the LD covering monolayer and the bilayer may be a key phenomenon influencing protein partitioning between the ER membrane and LDs in living cells.

Effects of Membrane PEGylation on Entry and Location of Antifungal Drug Itraconazole and Their Pharmacological Implications.

Itraconazole (ITZ) is an antifungal agent used clinically to treat mycotic infections. However, its therapeutic effects are limited by low solubility in aqueous media. Liposome-based delivery systems (LDS) have been proposed as a delivery mechanism for ITZ to alleviate this problem. Furthermore, PEGylation, the inclusion in the formulation of a protective "stealth sheath" of poly(ethylene glycol) around carrier particles, is widely used to increase circulation time in the bloodstream and hence efficacy. Together, these themes highlight the importance of mechanistic and structural understanding of ITZ incorporation into liposomes both with and without PEGylation because it can provide a potential foundation for the rational design of LDS-based systems for delivery of ITZ, using alternate protective polymers or formulations. Here we have combined atomistic simulations, cryo-TEM, Langmuir film balance, and fluorescence quenching experiments to explore how ITZ interacts with both pristine and PEGylated liposomes. We found that the drug can be incorporated into conventional and PEGylated liposomes for drug concentrations up to 15 mol % without phase separation. We observed that, in addition to its protective properties, PEGylation significantly increases the stability of liposomes that host ITZ. In a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer without PEGylation, ITZ was found to reside inside the lipid bilayer between the glycerol and the double-bond regions of POPC, adopting a largely parallel orientation along the membrane surface. In a PEGylated liposome, ITZ partitions mainly to the PEG layer. The results provide a solid basis for further development of liposome-based delivery systems.

doGlycans-Tools for Preparing Carbohydrate Structures for Atomistic Simulations of Glycoproteins, Glycolipids, and Carbohydrate Polymers for GROMACS.

Carbohydrates constitute a structurally and functionally diverse group of biological molecules and macromolecules. In cells they are involved in, e.g., energy storage, signaling, and cell-cell recognition. All of these phenomena take place in atomistic scales, thus atomistic simulation would be the method of choice to explore how carbohydrates function. However, the progress in the field is limited by the lack of appropriate tools for preparing carbohydrate structures and related topology files for the simulation models. Here we present tools that fill this gap. Applications where the tools discussed in this paper are particularly useful include, among others, the preparation of structures for glycolipids, nanocellulose, and glycans linked to glycoproteins. The molecular structures and simulation files generated by the tools are compatible with GROMACS.

Structure-phenotype correlations of human CYP21A2 mutations in congenital adrenal hyperplasia.

Mutations in the cytochrome p450 (CYP)21A2 gene, which encodes the enzyme steroid 21-hydroxylase, cause the majority of cases in congenital adrenal hyperplasia, an autosomal recessive disorder. To date, more than 100 CYP21A2 mutations have been reported. These mutations can be associated either with severe salt-wasting or simple virilizing phenotypes or with milder nonclassical phenotypes. Not all CYP21A2 mutations have, however, been characterized biochemically, and the clinical consequences of these mutations remain unknown. Using the crystal structure of its bovine homolog as a template, we have constructed a humanized model of CYP21A2 to provide comprehensive structural explanations for the clinical manifestations caused by each of the known disease-causing missense mutations in CYP21A2. Mutations that affect membrane anchoring, disrupt heme and/or substrate binding, or impair stability of CYP21A2 cause complete loss of function and salt-wasting disease. In contrast, mutations altering the transmembrane region or conserved hydrophobic patches cause up to a 98% reduction in enzyme activity and simple virilizing disease. Mild nonclassical disease can result from interference in oxidoreductase interactions, salt-bridge and hydrogen-bonding networks, and nonconserved hydrophobic clusters. A simple in silico evaluation of previously uncharacterized gene mutations could, thus, potentially help predict the often diverse phenotypes of a monogenic disorder.

How the amyloid-ß peptide and membranes affect each other: an extensive simulation study.

The etiology of Alzheimer's disease is thought to be linked to interactions between amyloid-ß (Aß) and neural cell membranes, causing membrane disruption and increased ion conductance. The effects of Aß on lipid behavior have been characterized experimentally, but structural and causal details are lacking. We used atomistic molecular dynamics simulations totaling over 6 µs in simulation time to investigate the behavior of Aß(42) in zwitterionic and anionic lipid bilayers. We simulated transmembrane ß-sheets (monomer and tetramer) resulting from a global optimization study and a helical structure obtained from an NMR study. In all simulations Aß(42) remained embedded in the bilayer. It was found that the surface charge and the lipid tail type are determinants for transmembrane stability of Aß(42) with zwitterionic surfaces and unsaturated lipids promoting stability. From the considered structures, the ß-sheet tetramer is most stable as a result of interpeptide interactions. We performed an in-depth analysis of the translocation of water in the Aß(42)-bilayer systems. We observed that this process is generally fast (within a few nanoseconds) yet generally slower than in the peptide-free bilayers. It is mainly governed by the lipid type, simulation temperature and Aß(42) conformation. The rate limiting step is the permeation through the hydrophobic core, where interactions between Aß(42) and permeating H(2)O molecules slow the translocation process. The ß-sheet tetramer allows more water molecules to pass through the bilayer compared to monomeric Aß, allowing us to conclude that the experimentally observed permeabilization of membranes must be due to membrane-bound Aß oligomers, and not monomers.

Hairpin protein partitioning from the ER to lipid droplets involves major structural rearrangements.

Lipid droplet (LD) function relies on proteins partitioning between the endoplasmic reticulum (ER) phospholipid bilayer and the LD monolayer membrane to control cellular adaptation to metabolic changes. It has been proposed that these hairpin proteins integrate into both membranes in a similar monotopic topology, enabling their passive lateral diffusion during LD emergence at the ER. Here, we combine biochemical solvent-accessibility assays, electron paramagnetic resonance spectroscopy and intra-molecular crosslinking experiments with molecular dynamics simulations, and determine distinct intramembrane positionings of the ER/LD protein UBXD8 in ER bilayer and LD monolayer membranes. UBXD8 is deeply inserted into the ER bilayer with a V-shaped topology and adopts an open-shallow conformation in the LD monolayer. Major structural rearrangements are required to enable ER-to-LD partitioning. Free energy calculations suggest that such structural transition is unlikely spontaneous, indicating that ER-to-LD protein partitioning relies on more complex mechanisms than anticipated and providing regulatory means for this trans-organelle protein trafficking.

Membrane permeation induced by aggregates of human islet amyloid polypeptides.

Several neurodegenerative diseases such as Alzheimer's and Parkinson's diseases as well as nonneuropathic diseases such as type II diabetes and atrial amyloidosis are associated with aggregation of amyloid polypeptides into fibrillar structures, or plaques. In this study, we use molecular dynamics simulations to test the stability and orientation of membrane-embedded aggregates of the human islet amyloid polypeptide (hIAPP) implicated in type II diabetes. We find that in both monolayers and bilayers of dipalmitoylphosphatidylglycerol (DPPG) hIAPP trimers and tetramers remain inside the membranes and preserve their ß-sheet secondary structure. Lipid bilayer-inserted hIAPP trimers and tetramers orient inside DPPG at 60° relative to the membrane/water interface and lead to water permeation and Na(+) intrusion, consistent with ion-toxicity in islet ß-cells. In particular, hIAPP trimers form a water-filled ß-sandwich that induce water permeability comparable with channel-forming proteins, such as aquaporins and gramicidin-A. The predicted disruptive orientation is consistent with the amphiphilic properties of the hIAPP aggregates and could be probed by chiral sum frequency generation (SFG) spectroscopy, as predicted by the simulated SFG spectra.

Is Lipid Specificity Key to the Potential Antiviral Activity of Mouthwash Reagent Chlorhexidine against SARS-CoV-2?

Chlorhexidine (CHX), a popular antibacterial drug, is widely used for oral health. Emerging pieces of evidence suggest that commercially available chlorhexidine mouthwash formulations are effective in suppressing the spread of SARS-CoV-2, possibly through destabilization of the viral lipid envelope. CHX is known for its membrane-active properties; however, the molecular mechanism revealing how it damages the viral lipid envelope is yet to be understood. Here we used extensive conventional and umbrella sampling simulations to quantify the effects of CHX on model membranes mimicking the composition of the SARS-CoV-2 outer lipid membrane as well as the host plasma membrane. Our results show that the lipid composition and physical properties of the membrane play an important role in binding and insertion, with CHX binding favorably to the viral membrane over the plasma membrane. Among the simulated lipids, CHX preferentially binds to anionic lipids, PS and PI, which are more concentrated in the viral membrane. The deeper and stable binding of CHX to the viral membrane results in more pronounced swelling of the membrane laterally with a thinning of the bilayer. The overall free energies of pore formation are strongly reduced for the viral membrane compared to the plasma membrane; however, CHX has a larger concentration-dependent effect on free energies of pore formation in the plasma membrane than the viral membrane. The results indicate that CHX is less toxic to the human plasma membrane at low concentrations. Our simulations reveal that CHX facilitates pore formation by the combination of thinning the membrane and accumulation at the water defect. This study provides insights into the mechanism underlying the anti-SARS-CoV-2 potency of CHX, supporting its potential for application as an effective and safe oral rinse agent for preventing viral transmission.

Free energies of membrane stalk formation from a lipidomics perspective.

Many biological membranes are asymmetric and exhibit complex lipid composition, comprising hundreds of distinct chemical species. Identifying the biological function and advantage of this complexity is a central goal of membrane biology. Here, we study how membrane complexity controls the energetics of the first steps of membrane fusions, that is, the formation of a stalk. We first present a computationally efficient method for simulating thermodynamically reversible pathways of stalk formation at coarse-grained resolution. The method reveals that the inner leaflet of a typical plasma membrane is far more fusogenic than the outer leaflet, which is likely an adaptation to evolutionary pressure. To rationalize these findings by the distinct lipid compositions, we computed ~200 free energies of stalk formation in membranes with different lipid head groups, tail lengths, tail unsaturations, and sterol content. In summary, the simulations reveal a drastic influence of the lipid composition on stalk formation and a comprehensive fusogenicity map of many biologically relevant lipid classes.

Cooperative Effects of an Antifungal Moiety and DMSO on Pore Formation over Lipid Membranes Revealed by Free Energy Calculations.

Itraconazole is a triazole drug widely used in the treatment of fungal infections, and it is in clinical trials for treatment of several cancers. However, the drug suffers from poor solubility, while experiments have shown that itraconazole delivery in liposome nanocarriers improves both circulation half-life and tissue distribution. The drug release mechanism from the nanocarrier is still unknown, and it depends on several factors including membrane stability against defect formation. In this work, we used molecular dynamics simulations and potential of mean force (PMF) calculations to quantify the influence of itraconazole on pore formation over lipid membranes, and we compared the effect by itraconazole with a pore-stabilizing effect by the organic solvent dimethyl sulfoxide (DMSO). According to the PMFs, both itraconazole and DMSO greatly reduce the free energy of pore formation, by up to ∼20 kJ mol-1. However, whereas large concentrations of itraconazole of 8 mol % (relative to lipid) were required, only small concentrations of a few mole % DMSO (relative to water) were sufficient to stabilize pores. In addition, itraconazole and DMSO facilitate pore formation by different mechanisms. Whereas itraconazole predominantly aids the formation of a partial defect with a locally thinned membrane, DMSO mainly stabilizes a transmembrane water needle by shielding it from the hydrophobic core. Notably, the two distinct mechanisms act cooperatively upon adding both itraconazole and DMSO to the membrane, as revealed by an additional reduction of the pore free energy. Overall, our simulations reveal molecular mechanisms and free energies of membrane pore formation by small molecules. We suggest that the stabilization of a locally thinned membrane as well as the shielding of a transmembrane water needle from the hydrophobic membrane core may be a general mechanism by which amphiphilic molecules facilitate pore formation over lipid membranes at sufficient concentrations.

Computer simulations of protein-membrane systems.

The interactions between proteins and membranes play critical roles in signal transduction, cell motility, and transport, and they are involved in many types of diseases. Molecular dynamics (MD) simulations have greatly contributed to our understanding of protein-membrane interactions, promoted by a dramatic development of MD-related software, increasingly accurate force fields, and available computer power. In this chapter, we present available methods for studying protein-membrane systems with MD simulations, including an overview about the various all-atom and coarse-grained force fields for lipids, and useful software for membrane simulation setup and analysis. A large set of case studies is discussed.

Cholesterol Reduces Partitioning of Antifungal Drug Itraconazole into Lipid Bilayers.

Cholesterol plays a crucial role in modulating the physicochemical properties of biomembranes, both increasing mechanical strength and decreasing permeability. Cholesterol is also a common component of vesicle-based delivery systems, including liposome-based drug delivery systems (LDSs). However, its effect on the partitioning of drug molecules to lipid membranes is very poorly recognized. Herein, we performed a combined experimental/computational study of the potential for the use of the LDS formulation for the delivery of the antifungal drug itraconazole (ITZ). We consider the addition of cholesterol to the lipid membrane. Since ITZ is only weakly soluble in water, its bioavailability is limited. Use of an LDS has thus been proposed. We studied lipid membranes composed of cholesterol, 1-palmitoyl-2-oleoyl-sn-glycerol-3-phosphocholine (POPC), and ITZ using a combination of computational molecular dynamics (MD) simulations of lipid bilayers and Brewster angle microscopy (BAM) experiments of monolayers. Both experimental and computational results show separation of cholesterol and ITZ. Cholesterol has a strong preference to orient parallel to the bilayer normal. However, ITZ, a long and relatively rigid molecule with weakly hydrophilic groups along the backbone, predominantly locates below the interface between the hydrocarbon chain region and the polar region of the membrane, with its backbone oriented parallel to the membrane surface; the orthogonal orientation in the membrane could be the cause of the observed separation. In addition, fluorescence measurements demonstrated that the affinity of ITZ for the lipid membrane is decreased by the presence of cholesterol, which is thus probably not a suitable formulation component of an LDS designed for ITZ delivery.

Role of Oxidized Gly25, Gly29, and Gly33 Residues on the Interactions of Aß1-42 with Lipid Membranes.

Oxidative stress is known to play an important role in the pathogenesis of Alzheimer's disease. Moreover, it is becoming increasingly evident that the plasma membrane of neurons plays a role in modulating the aggregation and toxicity of Alzheimer's amyloid-ß peptide (Aß). In this study, the combined and interdependent effects of oxidation and membrane interactions on the 42 residues long Aß isoform are investigated using molecular simulations. Hamiltonian replica exchange molecular dynamics simulations are utilized to elucidate the impact of selected oxidized glycine residues of Aß42 on the interactions of the peptide with a model membrane comprised of 70% POPC, 25% cholesterol, and 5% of the ganglioside GM1. The main findings are that, independent of the oxidation state, Aß prefers binding to GM1 over POPC, which is further enhanced by the oxidation of Gly29 and Gly33 and reduced the formation of ß-sheet. Our results suggest that the differences observed in Aß42 conformations and its interaction with a lipid bilayer upon oxidation originate from the position of the oxidized Gly residue with respect to the hydrophobic sequence of Aß42 involving the Gly29-XXX-Gly33-XXX-Gly37 motif and from specific interactions between the peptide and the terminal sugar groups of GM1.

Cryo-EM structure of the complete and ligand-saturated insulin receptor ectodomain.

Glucose homeostasis and growth essentially depend on the hormone insulin engaging its receptor. Despite biochemical and structural advances, a fundamental contradiction has persisted in the current understanding of insulin ligand-receptor interactions. While biochemistry predicts two distinct insulin binding sites, 1 and 2, recent structural analyses have resolved only site 1. Using a combined approach of cryo-EM and atomistic molecular dynamics simulation, we present the structure of the entire dimeric insulin receptor ectodomain saturated with four insulin molecules. Complementing the previously described insulin-site 1 interaction, we present the first view of insulin bound to the discrete insulin receptor site 2. Insulin binding stabilizes the receptor ectodomain in a T-shaped conformation wherein the membrane-proximal domains converge and contact each other. These findings expand the current models of insulin binding to its receptor and of its regulation. In summary, we provide the structural basis for a comprehensive description of ligand-receptor interactions that ultimately will inform new approaches to structure-based drug design.

Complexity of seemingly simple lipid nanodiscs.

Lipid nanodiscs are macromolecular assemblies, where a scaffold protein is wrapped around a nanosized disc of a lipid bilayer, thus protecting the hydrocarbon chains at the disc edges from unfavorable interactions with water. These nanostructures have numerous applications in, e.g., nanotechnology and pharmaceutics, and in investigations of membrane proteins. Here, we present results based on atomistic molecular dynamics simulations combined with electron paramagnetic spectroscopy measurements on the structure and dynamics of lipids in single-component nanodiscs. Our data highlight the existence of three distinctly different lipid fractions: central lipids residing in the center of a nanodisc, boundary lipids in direct contact with a scaffold protein, and intermediate lipids between these two regions. The central lipids are highly ordered and characterized by slow diffusion. In this part of the nanodisc, the membrane is the thickest and characterized by a gel-like or liquid-ordered phase, having features common to cholesterol-rich membranes. The boundary lipids in direct contact with the scaffold protein turned out to be less ordered and characterized by faster diffusion, and they remained in the liquid-disordered phase even at temperatures that were somewhat below the main phase transition temperature (Tm). The enthalpies associated with the central-boundary and central-intermediate transitions were similar to those observed for lipids going through the main phase transition. Overall, the study reveals lipid nanodiscs to be characterized by a complex internal structure, which is expected to influence membrane proteins placed in nanodiscs.

Mechanism of homodimeric cytokine receptor activation and dysregulation by oncogenic mutations.

Homodimeric class I cytokine receptors are assumed to exist as preformed dimers that are activated by ligand-induced conformational changes. We quantified the dimerization of three prototypic class I cytokine receptors in the plasma membrane of living cells by single-molecule fluorescence microscopy. Spatial and spatiotemporal correlation of individual receptor subunits showed ligand-induced dimerization and revealed that the associated Janus kinase 2 (JAK2) dimerizes through its pseudokinase domain. Oncogenic receptor and hyperactive JAK2 mutants promoted ligand-independent dimerization, highlighting the formation of receptor dimers as the switch responsible for signal activation. Atomistic modeling and molecular dynamics simulations based on a detailed energetic analysis of the interactions involved in dimerization yielded a mechanistic blueprint for homodimeric class I cytokine receptor activation and its dysregulation by individual mutations.

MPI-LIT: a literature-curated dataset of microbial binary protein--protein interactions.

UNLABELLED: Prokaryotic protein-protein interactions are underrepresented in currently available databases. Here, we describe a 'gold standard' dataset (MPI-LIT) focusing on microbial binary protein-protein interactions and associated experimental evidence that we have manually curated from 813 abstracts and full texts that were selected from an initial set of 36 852 abstracts. The MPI-LIT dataset comprises 1237 experimental descriptions that describe a non-redundant set of 746 interactions of which 659 (88%) are not reported in public databases. To estimate the curation quality, we compared our dataset with a union of microbial interaction data from IntAct, DIP, BIND and MINT. Among common abstracts, we achieve a sensitivity of up to 66% for interactions and 75% for experimental methods. Compared with these other datasets, MPI-LIT has the lowest fraction of interaction experiments per abstract (0.9) and the highest coverage of strains (92) and scientific articles (813). We compared methods that evaluate functional interactions among proteins (such as genomic context or co-expression) which are implemented in the STRING database. Most of these methods discriminate well between functionally relevant protein interactions (MPI-LIT) and high-throughput data. AVAILABILITY: http://www.jcvi.org/mpidb/interaction.php?dbsource=MPI-LIT. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.