Search

In the erythrocyte membrane, the interactions between glycophorin A (GPA) and Band 3 are associated strongly with the biological function of the membrane and several blood disorders. In this work, using coarse-grained molecular-dynamics simulations, we systematically investigate the effects of cholesterol and phosphatidylinositol-4,5-bisphosphate (PIP2) on the interactions of GPA with Band 3 in the model erythrocyte membranes. We examine the dynamics of the interactions of GPA with Band 3 in different lipid bilayers on the microsecond time scale and calculate the binding free energy between GPA and Band 3. The results indicate that cholesterols thermodynamically favor the binding of GPA to Band 3 by increasing the thickness of the lipid bilayer and by producing an effective attraction between the proteins due to the depletion effect. Cholesterols also slow the kinetics of the binding of GPA to Band 3 by reducing the lateral mobility of the lipids and proteins and may influence the binding sites between the proteins. The anionic PIP2 lipids prefer binding to the surface of the proteins through electrostatic attraction between the PIP2 headgroup and the positively charged residues on the protein surface. Ions in the solvent facilitate PIP2 aggregation, which promotes the binding of GPA to Band 3.

Simvastatin (Zocor) is one of the most prescribed drugs for reducing high cholesterol. Although simvastatin is ingested in its inactive lactone form, it is converted to its active dihydroxyheptanoate form by carboxylesterases in the liver. The dihydroxyheptanoate form can also be converted back to its original lactone form. Unfortunately, some of the side effects associated with the intake of simvastatin and other lipophilic statins at higher doses include statin-associated myopathy (SAM) and, in more severe cases, kidney failure. While the cause of SAM is unknown, it is hypothesized that these side effects are dependent on the localization of statins in lipid bilayers and their impact on bilayer properties. In this work, we carry out all-atom molecular dynamics simulations on both the lactone and dihydroxyheptanoate forms of simvastatin (termed "SN" and "SA", respectively) with a pure 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid bilayer and a POPC/cholesterol (30 mol %) binary mixture as membrane models. Additional simulations were carried out with multiple simvastatin molecules to mimic in vitro conditions that produced pleiotropic effects. Both SN and SA spontaneously diffused into the lipid bilayer, and a longer simulation time of 4 μs was needed for the complete incorporation of multiple SAs into the bilayer. By constructing potential mean force and electron density profiles, we find that SN localizes deeper within the hydrophobic interior of the bilayer and that SA has a greater tendency to form hydrogen-bonding interactions with neighboring water molecules and lipid headgroups. For the pure POPC bilayer, both SN and SA increase membrane order, while membrane fluidity increases for the POPC/cholesterol bilayer.

The lipid regulation of mammalian ion channel function has emerged as a fundamental mechanism in the control of electrical signalling and transport specificity in various cell types. In this work, we combine molecular dynamics simulations, mutagenesis, and electrophysiology to provide mechanistic insights into how lipophilic molecules (ceramide-sphingolipid probe) alter gating kinetics and K+ currents of hERG1. We show that the sphingolipid probe induced a significant left shift of activation voltage, faster deactivation rates, and current blockade comparable to traditional hERG1 blockers. Microseconds-long MD simulations followed by experimental mutagenesis elucidated ceramide specific binding locations at the interface between the pore and voltage sensing domains. This region constitutes a unique crevice present in mammalian channels with a non-swapped topology. The combined experimental and simulation data provide evidence for ceramide-induced allosteric modulation of the channel by a conformational selection mechanism., The lipid regulation of mammalian ion channel function has emerged as a fundamental mechanism in the control of electrical signalling and transport specificity. Here, the authors combine molecular dynamics simulations, mutagenesis, and electrophysiology to provide mechanistic insights into how lipophilic molecules alter gating kinetics and K+ currents of hERG1.

G protein-coupled receptors (GPCRs) are membrane-bound proteins that depend on their lipid environment to carry out their physiological function. Combined efforts from many theoretical and experimental studies on the lipid-protein interaction profile of several GPCRs hint at an intricate relationship of these receptors with their surrounding membrane environment, with several lipids emerging as particularly important. Using coarse-grained molecular dynamics simulations, we explore the lipid-protein interaction profiles of 28 different GPCRs, spanning different levels of classification and conformational states and totaling to 1 ms of simulation time. We find a close relationship with lipids for all GPCRs simulated, in particular, cholesterol and phosphatidylinositol phosphate (PIP) lipids, but the number, location, and estimated strength of these interactions is dependent on the specific GPCR as well as its conformational state. Although both cholesterol and PIP lipids bind specifically to GPCRs, they utilize distinct mechanisms. Interactions with PIP lipids are mediated by charge-charge interactions with intracellular loop residues and stabilized by one or both of the transmembrane helices linked by the loop. Interactions with cholesterol, on the other hand, are mediated by a hydrophobic environment, usually made up of residues from more than one helix, capable of accommodating its ring structure and stabilized by interactions with aromatic and charged/polar residues. Cholesterol binding to GPCRs occurs in a small number of sites, some of which (like the binding site on the extracellular side of transmembrane 6/7) are shared among many class A GPCRs. Combined with a thorough investigation of the local membrane structure, our results provide a detailed picture of GPCR-lipid interactions. Additionally, we provide an accompanying website to interactively explore the lipid-protein interaction profile of all GPCRs simulated to facilitate analysis and comparison of our data.

Deletion of phenylalanine 508 (F508del), in the cystic fibrosis transmembrane conductance regulator (CFTR) anion channel, is the most common cause of cystic fibrosis (CF). F508 is located on nucleotide-binding domain 1 (NBD1) in contact with cytosolic extensions of transmembrane helices, in particular intracellular loop 4 (ICL4). We carried out a mutagenesis scan of ICL4 by introducing five or six second-site mutations at eleven positions in cis with F508del, and quantifying changes in membrane proximity and ion-channel function of CFTR. The scan strongly validated the effectiveness of R1070W at rescuing F508del defects. Molecular dynamics simulations highlighted two features characterizing the ICL4/NBD1 interface of F508del/R1070W-CFTR: flexibility, with frequent transient formation of interdomain hydrogen bonds, and loosely stacked aromatic sidechains, (F1068, R1070W, and F1074, mimicking F1068, F508 and F1074 in wild-type CFTR). F508del-CFTR had a distorted aromatic stack, with F1068 displaced towards space vacated by F508. In F508del/R1070F-CFTR, which largely retained F508del defects, R1070F could not form hydrogen bonds, and the interface was less flexible. Other ICL4 second-site mutations which partially rescued F508del-CFTR are F1068M and F1074M. Methionine side chains allow hydrophobic interactions without the steric rigidity of aromatic rings, possibly conferring flexibility to accommodate the absence of F508 and retain a dynamic interface. Finally, two mutations identified in a yeast scan (A141S and R1097T, on adjacent transmembrane helices linked to ICL1 and ICL4) also partially rescued F508del-CFTR function. These studies highlight the importance of hydrophobic interactions and conformational flexibility at the ICL4/NBD1 interface, advancing understanding of the structural underpinning of F508del dysfunction.

Members of the ATP-binding cassette (ABC) transporter superfamily translocate a broad spectrum of chemically diverse substrates. While their eponymous ATP-binding cassette in the nucleotide-binding domains (NBDs) is highly conserved, their transmembrane domains (TMDs) forming the translocation pathway exhibit distinct folds and topologies, suggesting that during evolution, the ancient motor domains were combined with different transmembrane mechanical systems to orchestrate a variety of cellular processes. In recent years, it has become increasingly evident that the distinct TMD folds are best suited to categorize the multitude of ABC transporters. We therefore propose a new ABC transporter classification that currently comprises seven different types based on structural homology in the TMDs.

Specific and nonspecific protein-lipid interactions in cell membranes have important roles in an abundance of biological functions. We have used coarse-grained (CG) molecular dynamics (MD) simulations to assess lipid distributions and cholesterol flipping dynamics around surfaces in a model asymmetric plasma membrane containing one of six structurally distinct entities: aquaporin-1 (AQP1), the bacterial β-barrel outer membrane proteins OmpF and OmpX, KcsA potassium channel, WALP23 peptide, and a carbon nanotube (CNT). Our findings revealed varied lipid partitioning and cholesterol flipping times around the different solutes, and putative cholesterol binding sites in AQP1 and KcsA. The results suggest that protein-lipid interactions can be highly variable, and that surface-dependant lipid profiles are effectively manifested in CG simulations with the Martini force field.

The Ser/Thr protein kinase ataxia telangiectasia mutated (ATM) plays an important role in the DNA damage response, signaling in response to redox signals, the control of metabolic processes, and mitochondrial homeostasis. ATM localizes to the nucleus and at the plasma membrane, mitochondria, peroxisomes, and other cytoplasmic vesicular structures. It has been shown that the C-terminal FATC domain of human ATM (hATMfatc) can interact with a range of membrane mimetics and may thereby act as a membrane-anchoring unit. Here, NMR structural and N-15 relaxation data, NMR data using spin-labeled micelles, and MD simulations of micelle-associated hATMfatc revealed that it binds the micelle by a dynamic assembly of three helices with many residues of hATMfatc located in the headgroup region. We observed that none of the three helices penetrates the micelle deeply or makes significant tertiary contacts to the other helices. NMR-monitored interaction experiments with hATMfatc variants in which two conserved aromatic residues (Phe(3049) and Trp(3052)) were either individually or both replaced by alanine disclosed that the double substitution does not abrogate the interaction with micelles and bicelles at the high concentrations at which these aggregates are typically used, but impairs interactions with small unilamellar vesicles, usually used at much lower lipid concentrations and considered a better mimetic for natural membranes. We conclude that the observed dynamic structure of micelle-associated hATMfatc may enable it to interact with differently composed membranes or membrane-associated interaction partners and thereby regulate ATM's kinase activity. Moreover, the FATC domain of ATM may function as a membrane-anchoring unit for other biomolecules.

The cardiac KCNQ1 channel is a promising anti-arrhythmic target. Yazdi et al. report on how PUFAs interact with two binding sites in KCNQ1 to trigger channel activation. These findings further our mechanistic understanding of how to modulate KCNQ1 activity., Polyunsaturated fatty acids (PUFAs), but not saturated fatty acids, modulate ion channels such as the cardiac KCNQ1 channel, although the mechanism is not completely understood. Using both simulations and experiments, we find that PUFAs interact directly with the KCNQ1 channel via two different binding sites: one at the voltage sensor and one at the pore. These two amphiphilic binding pockets stabilize the negatively charged PUFA head group by electrostatic interactions with R218, R221, and K316, while the hydrophobic PUFA tail is selectively stabilized by cassettes of hydrophobic residues. The rigid saturated tail of stearic acid prevents close contacts with KCNQ1. By contrast, the mobile tail of PUFA linoleic acid can be accommodated in the crevice of the hydrophobic cassette, a defining feature of PUFA selectivity in KCNQ1. In addition, we identify Y268 as a critical PUFA anchor point underlying fatty acid selectivity. Combined, this study provides molecular models of direct interactions between PUFAs and KCNQ1 and identifies selectivity mechanisms. Long term, this understanding may open new avenues for drug development based on PUFA mechanisms.

The human ether-a-go-go-related gene (hERG) encodes the voltage gated potassium channel (KCNH2 or Kv11.1, commonly known as hERG). This channel plays a pivotal role in the stability of phase 3 repolarization of the cardiac action potential. Although a high-resolution cryo-EM structure is available for its depolarized (open) state, the structure surprisingly did not feature many functionally important interactions established by previous biochemical and electrophysiology experiments. Using Molecular Dynamics Flexible Fitting (MDFF), we refined the structure and recovered the missing functionally relevant salt bridges in hERG in its depolarized state. We also performed electrophysiology experiments to confirm the functional relevance of a novel salt bridge predicted by our refinement protocol. Our work shows how refinement of a high-resolution cryo-EM structure helps to bridge the existing gap between the structure and function in the voltage-sensing domain (VSD) of hERG.Statement of SignificanceCryo-EM has emerged as a major breakthrough technique in structural biology of membrane proteins. However, even high-resolution Cryo-EM structures contain poor side chain conformations and interatomic clashes. A high-resolution cryo-EM structure of hERG1 has been solved in the depolarized (open) state. The state captured by Cryo-EM surprisingly did not feature many functionally important interactions established by previous experiments. Molecular Dynamics Flexible Fitting (MDFF) used to enable refinement of the hERG1 channel structure in complex membrane environment re-establishing key functional interactions in the voltage sensing domain.

The hydrophobic surfactant proteins, SP-B and SP-C, promote rapid adsorption by the surfactant lipids to the surface of the liquid that lines the alveolar air sacks of the lungs. To gain insights into the mechanisms of their function, we used X-ray diffuse scattering (XDS) and molecular dynamics (MD) simulations to determine the location of SP-B and SP-C within phospholipid bilayers. Initial samples contained the surfactant lipids from extracted calf surfactant with increasing doses of the proteins. XDS located protein density near the phospholipid headgroup and in the hydrocarbon core, presumed to be SP-B and SP-C, respectively. Measurements on dioleoylphosphatidylcholine (DOPC) with the proteins produced similar results. MD simulations of the proteins with DOPC provided molecular detail and allowed direct comparison of the experimental and simulated results. Simulations used conformations of SP-B based on other members of the saposin-like family, which form either open or closed V-shaped structures. For SP-C, the amino acid sequence suggests a partial α-helix. Simulations fit best with measurements of XDS for closed SP-B, which occurred at the membrane surface, and SP-C oriented along the hydrophobic interior. Our results provide the most definitive evidence yet concerning the location and orientation of the hydrophobic surfactant proteins.

Band 3 is an anion exchanger for chloride/bicarbonate in the plasma membrane of erythrocytes. The function of Band 3 may be influenced by the interactions between Band 3 and lipids or lipid domains in the plasma membrane. In this work, using atomistic molecular dynamics simulation, we investigate the interactions between Band 3 and nanosized lipid domains in ternary lipid bilayers composed of saturated lipids, unsaturated lipids, and cholesterol. The simulations show asymmetric interactions between Band 3 and lipid nanodomains in the two leaflets of a neutral lipid bilayer with lower cholesterol concentration. With an increase in cholesterol concentration in the bilayer, cholesterol affects the interactions between Band 3 and lipid domains by deforming the structure of the protein. Additionally, the anionic lipids, which prefer to bind to some specific sites of Band 3, also affect the interactions between Band 3 and lipid domains. This work provides some new insight into understanding the distribution of Band 3 in the plasma membrane of erythrocytes as well as its anion exchange function.

Anthrax toxin action requires triggering of natural endocytic transport mechanisms whereby the binding component of the toxin forms channels (PA63) within endosomal limiting and intraluminal vesicle membranes to deliver the toxin's enzymatic components into the cytosol. Membrane lipid composition varies at different stages of anthrax toxin internalization, with intraluminal vesicle membranes containing ~70% of anionic bis(monoacylglycero)phosphate lipid. Using model bilayer measurements, we show that membrane lipids can have a strong effect on the anthrax toxin channel properties, including the channel-forming activity, voltage-gating, conductance, selectivity, and enzymatic factor binding. Interestingly, the highest PA63 insertion rate was observed in bis(monoacylglycero)phosphate membranes. The molecular dynamics simulation data show that the conformational properties of the channel are different in bis(monoacylglycero)phosphate compared to PC, PE, and PS lipids. The anthrax toxin protein/lipid bilayer system can be advanced as a novel robust model to directly investigate lipid influence on membrane protein properties and protein/protein interactions.

The activity of CTP:phosphocholine cytidylyltransferase (CCT), a key enzyme in phosphatidylcholine synthesis, is regulated by reversible interactions of a lipid-inducible amphipathic helix (domain M) with membrane phospholipids. When dissociated from membranes, a portion of the M domain functions as an auto-inhibitory (AI) element to suppress catalysis. The AI helix from each subunit binds to a pair of α helices (αE) that extend from the base of the catalytic dimer to create a four-helix bundle. The bound AI helices make intimate contact with loop L2, housing a key catalytic residue, Lys122. The impacts of the AI helix on active-site dynamics and positioning of Lys122 are unknown. Extensive MD simulations with and without the AI helix revealed that backbone carbonyl oxygens at the point of contact between the AI helix and loop L2 can entrap the Lys122 side chain, effectively competing with the substrate, CTP. In silico, removal of the AI helices dramatically increased αE dynamics at a predicted break in the middle of these helices, enabling them to splay apart and forge new contacts with loop L2. In vitro cross-linking confirmed the reorganization of the αE element upon membrane binding of the AI helix. Moreover, when αE bending was prevented by disulfide engineering, CCT activation by membrane binding was thwarted. These findings suggest a novel two-part auto-inhibitory mechanism for CCT involving capture of Lys122 and restraint of the pliable αE helices. We propose that membrane binding enables bending of the αE helices, bringing the active site closer to the membrane surface.

P-glycoprotein contributes to multidrug resistance by exporting a broad range of substrates across the cell membrane. Using molecular dynamics simulations, Barreto-Ojeda et al. identify key lipid-binding sites and reveal lipid access pathways toward the cavity of the transporter., P-glycoprotein (P-gp) exports a broad range of dissimilar compounds, including drugs, lipids, and lipid-like molecules. Because of its substrate promiscuity, P-gp is a key player in the development of cancer multidrug resistance. Although P-gp is one of the most studied ABC transporters, the mechanism by which its substrates access the cavity remains unclear. In this study, we perform coarse-grained molecular dynamics simulations to explore possible lipid access pathways in the inward-facing conformation of P-gp embedded in bilayers of different lipid compositions. In the inward-facing orientation, only lipids from the lower leaflet access the cavity of the transporter. We identify positively charged residues at the portals of P-gp that favor lipid entrance to the cavity, as well as lipid-binding sites at the portals and within the cavity, which is in good agreement with previous experimental studies. This work includes several examples of lipid pathways for phosphatidylcholine and phosphatidylethanolamine lipids that help elucidate the molecular mechanism of lipid binding in P-gp.

Cardiac troponin C (cTnC) is the regulatory protein that initiates cardiac contraction in response to Ca2+. TnC binding Ca2+ initiates a cascade of protein–protein interactions that begins with the opening of the N-terminal domain of cTnC, followed by cTnC binding the troponin I switch peptide (TnISW). We have evaluated, through isothermal titration calorimetry and molecular-dynamics simulation, the effect of several clinically relevant mutations (A8V, L29Q, A31S, L48Q, Q50R, and C84Y) on the Ca2+ affinity, structural dynamics, and calculated interaction strengths between cTnC and each of Ca2+ and TnISW. Surprisingly the Ca2+ affinity measured by isothermal titration calorimetry was only significantly affected by half of these mutations including L48Q, which had a 10-fold higher affinity than WT, and the Q50R and C84Y mutants, each of which had affinities 3-fold higher than wild type. This suggests that Ca2+ affinity of the N-terminal domain of cTnC in isolation is insufficient to explain the pathogenicity of these mutations. Molecular-dynamics simulation was used to evaluate the effects of these mutations on Ca2+ binding, structural dynamics, and TnI interaction independently. Many of the mutations had a pronounced effect on the balance between the open and closed conformations of the TnC molecule, which provides an indirect mechanism for their pathogenic properties. Our data demonstrate that the structural dynamics of the cTnC molecule are key in determining myofilament Ca2+ sensitivity. Our data further suggest that modulation of the structural dynamics is the underlying molecular mechanism for many disease mutations that are far from the regulatory Ca2+-binding site of cTnC.

Crystallographic data of the dimeric and octameric forms of fragaceatoxin C (FraC) suggested the key role of a small hydrophobic protein-protein interaction surface for actinoporins oligomerization and pore formation in membranes. However, site-directed mutagenesis studies supporting this hypothesis for others actinoporins are still lacking. Here, we demonstrate that disrupting the key hydrophobic interaction between V60 and F163 (FraC numbering scheme) in the oligomerization interface of FraC, equinatoxin II (EqtII) and sticholysin II (StII) impairs the pore formation activity of these proteins. Our results allow for the extension of the importance of FraC protein-protein interactions in the stabilization of the oligomeric intermediates of StII and EqtII pointing out that all of these proteins follow a similar pathway of membrane disruption. These findings support the hybrid pore proposal as the universal model of actinoporins pore formation. Moreover, we reinforce the relevance of dimer formation, which appears to be a functional intermediate in the assembly pathway of some different pore-forming proteins. This article is protected by copyright. All rights reserved.

The dynamics of phosphocholine and maltoside micelles, detergents frequently used for membrane protein structure determination, were investigated using electron paramagnetic resonance of spin probes doped into the micelles. Specifically, phosphocholines are frequently used detergents in NMR studies, and maltosides are frequently used in x-ray crystallography structure determination. Beyond the structural and electrostatic differences, this study aimed to determine whether there are differences in the local chain dynamics (i.e., fluidity). The nitroxide probe rotational dynamics in longer chain detergents is more restricted than in shorter chain detergents, and maltoside micelles are more restricted than phosphocholine micelles. Furthermore, the micelle microviscosity can be modulated with mixtures, as demonstrated with mixtures of 3-[(3-cholamidopropyl)dimethylammonio]-1-propanesulfonate with n-dodecylphosphocholine, n-tetradecylphosphocholine, n-decyl-β-D-maltoside, or n-dodecyl-β-D-maltoside. These results indicate that observed differences in membrane protein stability in these detergents could be due to fluidity in addition to the already determined structural differences.

Lung surfactant lines the gas-exchange interface in the lungs and reduces the surface tension, which is necessary for breathing. Lung surfactant consists mainly of lipids with a small amount of proteins and forms a monolayer at the air-water interface connected to bilayer reservoirs. Lung surfactant function involves transfer of material between the monolayer and bilayers during the breathing cycle. Lipids and proteins are organized laterally in the monolayer; selected species are possibly preferentially transferred to bilayers. The complex 3D structure of lung surfactant and the exact roles of lipid organization and proteins remain important goals for research. We review recent simulation studies on the properties of lipid monolayers, monolayers with phase coexistence, monolayer-bilayer transformations, lipid-protein interactions, and effects of nanoparticles on lung surfactant. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Cholesterol is the most abundant molecule in the plasma membrane of mammals. Its distribution across the two membrane leaflets is critical for understanding how cells work. Cholesterol trans-bilayer motion (flip-flop) is a key process influencing its distribution in membranes. Despite extensive investigations, the rate of cholesterol flip-flop and its dependence on the lateral heterogeneity of membranes remain uncertain. In this work, we used atomistic molecular dynamics simulations to sample spontaneous cholesterol flip-flop events in a DPPC:DOPC:cholesterol mixture with heterogeneous lateral distribution of lipids. In addition to an overall flip-flop rate at the time scale of sub-milliseconds, we identified a significant impact of local environment on flip-flop rate. We discuss the atomistic details of the flip-flop events observed in our simulations.

Membrane curvature controls the spatial organization and activity of cells. Lipid sorting in cell membranes can be explained by matching lipid molecular shape to regions of different curvatures. A molecular view of curvature-induced lipid sorting is obtained using coarse-grained molecular dynamics. A model membrane consisting of an asymmetric bilayer of multiple lipid species is simulated. Curvature is induced by pulling a tether, that is, a bilayer nanotube, from a flat membrane. Pulling is performed both from the inner and outer leaflets, corresponding to directions in and out of the cell. Redistribution of different lipid types between the tether and the bilayer is observed, leading to spatial variations in the composition of both leaflets, and, in turn, softening of the tether. Depending on the direction of pulling, the lipid distributions and the tether properties differ. Formation of a tether from the planar membrane thus induces lipid sorting without phase separation.

Cell membranes contain hundreds of different proteins and lipids in an asymmetric arrangement. Our current understanding of the detailed organization of cell membranes remains rather elusive, because of the challenge to study fluctuating nanoscale assemblies of lipids and proteins with the required spatiotemporal resolution. Here, we use molecular dynamics simulations to characterize the lipid environment of 10 different membrane proteins. To provide a realistic lipid environment, the proteins are embedded in a model plasma membrane, where more than 60 lipid species are represented, asymmetrically distributed between the leaflets. The simulations detail how each protein modulates its local lipid environment in a unique way, through enrichment or depletion of specific lipid components, resulting in thickness and curvature gradients. Our results provide a molecular glimpse of the complexity of lipid–protein interactions, with potentially far-reaching implications for our understanding of the overall organization of real cell membranes., Computer simulations show how 10 different membrane proteins representing major protein families uniquely shape the membrane, which may be an important organizing principle for cell membranes.

The cystic fibrosis transmembrane conductance regulator (CFTR) is an ion channel that regulates the flow of anions across epithelia. Mutations in CFTR cause cystic fibrosis. CFTR belongs to the ATP-Binding Cassette (ABC) transporter superfamily, and gating is controlled by phosphorylation and ATP binding and hydrolysis. Recent ATP-free and ATP-bound structures of zebrafish CFTR revealed an unwound segment of transmembrane helix (TM) 8, which appears to be a unique feature of CFTR not present in other ABC transporter structures. Here, by means of 1 μs long molecular dynamics simulations, we investigate the interactions formed by this TM8 segment with nearby helices, in both ATP-free and ATP-bound states. We highlight the structural role of TM8 in maintaining the functional architecture of the pore and we describe a distinct membrane defect that is found near TM8 only in the ATP-free state. The results of the MD simulations are discussed in the context of the gating mechanism of CFTR.

Cholesterol is a key component of eukaryotic membranes, but its role in cellular biology in general and in lipid rafts in particular remains controversial. Model membranes are used extensively to determine the phase behavior of ternary mixtures of cholesterol, a saturated lipid, and an unsaturated lipid with liquid-ordered and liquid-disordered phase coexistence. Despite many different experiments that determine lipid-phase diagrams, we lack an understanding of the molecular-level driving forces for liquid phase coexistence in bilayers with cholesterol. Here, we use atomistic molecular dynamics computer simulations to address the driving forces for phase coexistence in ternary lipid mixtures. Domain formation is directly observed in a long-timescale simulation of a mixture of 1,2-distearoyl-sn-glycero-3-phosphocholine, unsaturated 1,2-dilinoleoyl-sn-glycero-3-phosphocholine, and cholesterol. Free-energy calculations for the exchange of the saturated and unsaturated lipids between the ordered and disordered phases give insight into the mixing behavior. We show that a large energetic contribution to domain formation is favorable enthalpic interactions of the saturated lipid in the ordered phase. This favorable energy for forming an ordered, cholesterol-rich phase is opposed by a large unfavorable entropy. Martini coarse-grained simulations capture the unfavorable free energy of mixing but do not reproduce the entropic contribution because of the reduced representation of the phospholipid tails. Phospholipid tails and their degree of unsaturation are key energetic contributors to lipid phase separation.

The detailed organization of cellular membranes remains rather elusive. Based on large-scale molecular dynamics simulations, we provide a high-resolution view of the lipid organization of a plasma membrane at an unprecedented level of complexity. Our plasma membrane model consists of 63 different lipid species, combining 14 types of headgroups and 11 types of tails asymmetrically distributed across the two leaflets, closely mimicking an idealized mammalian plasma membrane. We observe an enrichment of cholesterol in the outer leaflet and a general non-ideal lateral mixing of the different lipid species. Transient domains with liquid-ordered character form and disappear on the microsecond time scale. These domains are coupled across the two membrane leaflets. In the outer leaflet, distinct nanodomains consisting of gangliosides are observed. Phosphoinositides show preferential clustering in the inner leaflet. Our data provide a key view on the lateral organization of lipids in one of life's fundamental structures, the cell membrane.

Collapse of homogeneous lipid monolayers is known to proceed via wrinkling/buckling, followed by folding into bilayers in water. For heterogeneous monolayers with phase coexistence, the mechanism of collapse remains unclear. Here, we investigated collapse of lipid monolayers with coexisting liquid-liquid and liquid-solid domains using molecular dynamics simulations. The MARTINI coarse-grained model was employed to simulate monolayers of ∼80 nm in lateral dimension for 10–25 μs. The monolayer minimum surface tension decreased in the presence of solid domains, especially if they percolated. Liquid-ordered domains facilitated monolayer collapse due to the spontaneous curvature induced at a high cholesterol concentration. Upon collapse, bilayer folds formed in the liquid (disordered) phase; curved domains shifted the nucleation sites toward the phase boundary. The liquid (disordered) phase was preferentially transferred into bilayers, in agreement with the squeeze-out hypothesis. As a result, the composition and phase distribution were altered in the monolayer in equilibrium with bilayers compared to a flat monolayer at the same surface tension. The composition and phase behavior of the bilayers depended on the degree of monolayer compression. The monolayer-bilayer connection region was enriched in unsaturated lipids. Percolation of solid domains slowed down monolayer collapse by several orders of magnitude. These results are important for understanding the mechanism of two-to-three-dimensional transformations in heterogeneous thin films and the role of lateral organization in biological membranes. The study is directly relevant for the function of lung surfactant, and can explain the role of nanodomains in its surface activity and inhibition by an increased cholesterol concentration.

Bioenergetic membranes are the key cellular structures responsible for coupled energy-conversion processes, which supply ATP and important metabolites to the cell. Here, we report the first 100-million atom-scale model of an entire photosynthetic organelle, a chromatophore membrane vesicle from a purple bacterium, which reveals the rate-determining steps of membrane-mediated energy conversion. Molecular dynamics simulations of this bioenergetic organelle elucidate how the network of bioenergetic proteins influences membrane curvature and demonstrates the impact of thermal disorder on photosynthetic excitation transfer. Brownian dynamics simulations of the quinone and cytochrome c(2) charge carriers within the chromatophore interior probe the mechanisms of nanoscale charge transport under various pH and salinity conditions. Reproducing phenotypic properties from atomistic details, a rate-kinetic model evinces that low-light adaptations of the bacterium emerge as a spontaneous outcome of optimizing the balance between the chromatophore’s structural integrity and robust energy conversion. Put together, the hybrid structure determination and systems-level modeling of the chromatophore, in conjunction with optical spectroscopy, illuminate the chemical and organizational design principles of biological membranes that foster energy storage and transduction in living cells. Parallels are drawn with the more universal mitochondrial bioenergetic machinery, from whence molecular-scale insights on the mechanism of cellular aging are inferred. This endeavor made feasible through the advent of petascale supercomputers, paves the way to first-principles modeling of whole living cells.

Small-angle scattering is a powerful technique that can probe the structure of lipid bilayers on the nanometer scale. Retrieving the real space structure of lipid bilayers from the scattering intensity can be a challenging task, as their fluid nature results in a liquid-like scattering pattern which is hard to interpret. The standard approach to this problem is to describe the bilayer structure as a sum of density distributions of separate components of the lipid molecule and then to fit the parameters of the distributions against experimental data. The accuracy of the density-based analysis is partially limited by the choice of the functions used to describe component distributions, especially in the case of multi-component bilayers. The number of parameters in the model is balanced by the need for an accurate description of the underlying bilayer structure and the risk of overfitting the data. Here, we present an alternative method for the interpretation of small-angle scattering intensity data for lipid bilayers. The method is based on restrained ensemble molecular dynamics simulations that allow direct incorporation of the scattering data into the simulations in the form of a restraining potential. This approach combines the information implicitly contained in the simulation force field with structural data from the scattering intensity and is free from prior assumptions regarding the bilayer structure.

Cellular membranes separate distinct aqueous compartments, but can be breached by transient hydrophilic pores. A large energetic cost prevents pore formation, which is largely dependent on the composition and structure of the lipid bilayer. The softness of bilayers and the disordered structure of pores make their characterization difficult. We use molecular-dynamics simulations with atomistic detail to study the thermodynamics, kinetics, and mechanism of pore formation and closure in DLPC, DMPC, and DPPC bilayers, with pore formation free energies of 17, 45, and 78 kJ/mol, respectively. By using atomistic computer simulations, we are able to determine not only the free energy for pore formation, but also the enthalpy and entropy, which yields what is believed to be significant new insights in the molecular driving forces behind membrane defects. The free energy cost for pore formation is due to a large unfavorable entropic contribution and a favorable change in enthalpy. Changes in hydrogen bonding patterns occur, with increased lipid-water interactions, and fewer water-water hydrogen bonds, but the total number of overall hydrogen bonds is constant. Equilibrium pore formation is directly observed in the thin DLPC lipid bilayer. Multiple long timescale simulations of pore closure are used to predict pore lifetimes. Our results are important for biological applications, including the activity of antimicrobial peptides and a better understanding of membrane protein folding, and improve our understanding of the fundamental physicochemical nature of membranes.

Gangliosides are glycolipids in which an oligosaccharide headgroup containing one or more sialic acids is connected to a ceramide. Gangliosides reside in the outer leaflet of the plasma membrane and play a crucial role in various physiological processes such as cell signal transduction and neuronal differentiation by modulating structures and functions of membrane proteins. Because the detailed behavior of gangliosides and protein-ganglioside interactions are poorly known, we investigated the interactions between the gangliosides GM1 and GM3 and the proteins aquaporin (AQP1) and WALP23 using equilibrium molecular dynamics simulations and potential of mean force calculations at both coarse-grained (CG) and atomistic levels. In atomistic simulations, on the basis of the GROMOS force field, ganglioside aggregation appears to be a result of the balance between hydrogen bond interactions and steric hindrance of the headgroups. GM3 clusters are slightly larger and more ordered than GM1 clusters due to the smaller headgroup of GM3. The different structures of GM1 and GM3 clusters from atomistic simulations are not observed at the CG level based on the Martini model, implying a difference in driving forces for ganglioside interactions in atomistic and CG simulations. For protein-ganglioside interactions, in the atomistic simulations, GM1 lipids bind to specific sites on the AQP1 surface, whereas they are depleted from WALP23. In the CG simulations, the ganglioside binding sites on the AQP1 surface are similar, but ganglioside aggregation and protein-ganglioside interactions are more prevalent than in the atomistic simulations. Using the polarizable Martini water model, results were closer to the atomistic simulations. Although experimental data for validation is lacking, we proposed modified Martini parameters for gangliosides to more closely mimic the sizes and structures of ganglioside clusters observed at the atomistic level.

Cell membranes contain multiple lipid and protein components having heterogeneous in-plane (lateral) distribution. Nanoscale rafts are believed to play an important functional role, but their phase state—domains of coexisting phases or composition fluctuations—is unknown. As a step toward understanding lateral organization of cell membranes, we investigate the difference between nanoscale domains of coexisting phases and composition fluctuations in lipid bilayers. We simulate model lipid bilayers with the MARTINI coarse-grained force field on length scales of tens of nanometers and timescales of tens of microseconds. We use a binary and a ternary mixture: a saturated and an unsaturated lipid, or a saturated lipid, an unsaturated lipid, and cholesterol, respectively. In these mixtures, the phase behavior can be tuned from a mixed state to a coexistence of a liquid-crystalline and a gel, or a liquid-ordered and a liquid-disordered phase. Transition from a two-phase to a one-phase state is achieved by raising the temperature and adding a hybrid lipid (with a saturated and an unsaturated chain). We analyze the evolution of bilayer properties along this transition: domains of two phases transform to fluctuations with local ordering and compositional demixing. Nanoscale domains and fluctuations differ in several properties, including interleaflet overlap and boundary length. Hybrid lipids show no enrichment at the boundary, but decrease the difference between the coexisting phases by ordering the disordered phase, which could explain their role in cell membranes.

Due to antimicrobial resistance, the development of new drugs to combat bacterial and fungal infections is an important area of research. Nature uses short, charged, and amphipathic peptides for antimicrobial defense, many of which disrupt the lipid membrane in addition to other possible targets inside the cell. Computer simulations have revealed atomistic details for the interactions of antimicrobial peptides and cell-penetrating peptides with lipid bilayers. Strong interactions between the polar interface and the charged peptides can induce bilayer deformations - including membrane rupture and peptide stabilization of a hydrophilic pore. Here, we performed microsecond-long simulations of the antimicrobial peptide CM15 in a POPC bilayer expecting to observe pore formation (based on previous molecular dynamics simulations). We show that caution is needed when interpreting results of equilibrium peptide-membrane simulations, given the length of time single trajectories can dwell in local energy minima for 100's of ns to microseconds. While we did record significant membrane perturbations from the CM15 peptide, pores were not observed. We explain this discrepancy by computing the free energy for pore formation with different force fields. Our results show a large difference in the free energy barrier (ca. 40 kJ/mol) against pore formation predicted by the different force fields that would result in orders of magnitude differences in the simulation time required to observe spontaneous pore formation. This explains why previous simulations using the Berger lipid parameters reported pores induced by charged peptides, while with CHARMM based models pores were not observed in our long time-scale simulations. We reconcile some of the differences in the distance dependent free energies by shifting the free energy profiles to account for thickness differences between force fields. The shifted curves show that all the models describe small defects in lipid bilayers in a consistent manner, suggesting a common physical basis.

Helix A of bacteriorhodopsin was simulated by using molecular dynamics both in isolation and as part of the complete protein. A POPC lipid bilayer and an octane monolayer were used as model membranes. Comparison of various systems showed octane to be a good alternative to lipid bilayer membranes providing fast equilibration, increased sampling, and decreased computational cost. Similarly single-helix simulations were found to capture some of the details of the full-protein simulations. In particular, aromatic side chains were found to anchor in identical conformations in all simulations, regardless of the absence of a lipid-water interface or the remainder of the protein. Simulations displayed a remarkable robustness with respect to simulation parameters and system set up.

Zebrafish, as a model for teleost fish, have two paralogous troponin C (TnC) genes that are expressed in the heart differentially in response to temperature acclimation. Upon Ca(2+) binding, TnC changes conformation and exposes a hydrophobic patch that interacts with troponin I and initiates cardiac muscle contraction. Teleost-specific TnC paralogs have not yet been functionally characterized. In this study we have modeled the structures of the paralogs using molecular dynamics simulations at 18°C and 28°C and calculated the different Ca(2+)-binding properties between the teleost cardiac (cTnC or TnC1a) and slow-skeletal (ssTnC or TnC1b) paralogs through potential-of-mean-force calculations. These values are compared with thermodynamic binding properties obtained through isothermal titration calorimetry (ITC). The modeled structures of each of the paralogs are similar at each temperature, with the exception of helix C, which flanks the Ca(2+) binding site; this region is also home to paralog-specific sequence substitutions that we predict have an influence on protein function. The short timescale of the potential-of-mean-force calculation precludes the inclusion of the conformational change on the ΔG of Ca(2+) interaction, whereas the ITC analysis includes the Ca(2+) binding and conformational change of the TnC molecule. ITC analysis has revealed that ssTnC has higher Ca(2+) affinity than cTnC for Ca(2+) overall, whereas each of the paralogs has increased affinity at 28°C compared to 18°C. Microsecond-timescale simulations have calculated that the cTnC paralog transitions from the closed to the open state more readily than the ssTnC paralog, an unfavorable transition that would decrease the ITC-derived Ca(2+) affinity while simultaneously increasing the Ca(2+) sensitivity of the myofilament. We propose that the preferential expression of cTnC at lower temperatures increases myofilament Ca(2+) sensitivity by this mechanism, despite the lower Ca(2+) affinity that we have measured by ITC.

Most promising gene therapy approaches use siRNA. Sophisticated delivery systems are required to protect and deliver siRNA effectively to the target tissues. Lipid nanoparticles (LNPs) containing ionizable amino lipids with pKa of 7 or lower, e.g. DLinKC2DMA (XTC2) or DLinMC3DMA (MC3), are currently the leading delivery systems for siRNA in liver-diseases, including cancer. Despite a high encapsulation efficiency of siRNA-LNPs, currently, only 1% of the siRNA in internalized LNPs is released in the cytoplasm while the balance is transferred to the lysosome and degraded. Therefore, further optimization can be done on drug release from endosomes.Using molecular dynamics simulation and NMR experiments we aim to better understand interactions between LNP's cationic lipids and endosomal anionic lipids, which have been proposed to destabilize endosomal membranes, allowing drug release. We investigate bilayers composed of XTC2 and DSPS lipids as ionizable amino lipid and endosomal anionic lipid model, respectively, for a variety of mixing ratios and environment conditions, e.g. salt concentration, pH, and temperature. Lipid order parameters are calculated computationally and experimentally to validate our computational models. Further analysis on simulated systems provides us atomistic details about the structure of and interactions in the system.Our results will assist in rational design of new generations of LNPs with improved release efficiencies.

Biological membranes are crucial as they define essential compartments for cells. Cell membranes are heterogeneous mixtures of membrane proteins and lipids. Their structure and function are fundamentally dependent on lipid-lipid and lipid-protein interactions which play a crucial role in regulating cell protein functions and are involved in many diseases when altered. However, the exact properties and the relation among the distinct components that form the biological membrane are not completely understood due the limitations of experimental methods in studying the lipid-protein interactions in living cells. Lately, computer simulations have become a promising tool to understand and clarify experimental results. The goal of our research is to identify specific lipid-protein interactions of biological relevance and work towards large-scale modeling of realistic biological membranes. In this study, we placed ten main classes of eukaryotic membrane proteins with different ratios in a prototypical plasma membrane model, containing various lipid types found in the plasma membranes of eukaryotic cells. The protein-membrane system was prepared using Martini force field. Using molecular dynamic simulations, we intend to investigate how proteins sort lipids, and to identify specific binding partners and lipids of interest for more detailed studies.

Glycolipids, which are composed of oligosaccharide headgroups and glycerol or sphingosine backbones, are important components of plasma membranes. Gangliosides are glycolipids in which an oligosaccharide headgroup containing one or more sialic acids is connected to ceramide. The negatively charged headgroups of gangliosides extend out of the plasma membrane extracellular side and play a crucial role in cell-cell recognition. Since the behavior of gangliosides and their interactions with proteins are poorly known, we conducted molecular dynamics simulations at both the coarse-grained and atomistic levels. We investigated the interactions between the gangliosides GM1 and GM3 and the proteins aquaporin and WALP23. GM1 and GM3 aggregated in clusters in POPC bilayers containing 17% gangliosides. These clusters pack around the protein and peptide in the bilayer. The cluster size is much smaller in polarizable MARTINI water than in standard MARTINI water. In atomistic simulations, the large GM1 and GM3 clusters found in standard MARTINI water also fell apart to clusters of smaller size. Gangliosides and ganglioside-protein interaction free energies in POPC bilayers obtained with polarizable MARTINI water are roughly half of those with standard MARTINI water. Corresponding free energies in atomistic simulations revealed different free energy landscapes and even weaker interactions.

Lipid nanoparticles (LNP) containing ionizable cationic lipids are the leading systems for enabling therapeutic applications of siRNA; however, the structure of these systems has not been defined. Here we examine the structure of LNP siRNA systems containing DLinKC2-DMA(an ionizable cationic lipid), phospholipid, cholesterol and a polyethylene glycol (PEG) lipid formed using a rapid microfluidic mixing process. Techniques employed include cryo-transmission electron microscopy, (31)P NMR, membrane fusion assays, density measurements, and molecular modeling. The experimental results indicate that these LNP siRNA systems have an interior lipid core containing siRNA duplexes complexed to cationic lipid and that the interior core also contains phospholipid and cholesterol. Consistent with experimental observations, molecular modeling calculations indicate that the interior of LNP siRNA systems exhibits a periodic structure of aqueous compartments, where some compartments contain siRNA. It is concluded that LNP siRNA systems formulated by rapid mixing of an ethanol solution of lipid with an aqueous medium containing siRNA exhibit a nanostructured core. The results give insight into the mechanism whereby LNP siRNA systems are formed, providing an understanding of the high encapsulation efficiencies that can be achieved and information on methods of constructing more sophisticated LNP systems.

Membrane tethers are nanotubes formed by a lipid bilayer. They play important functional roles in cell biology and provide an experimental window on lipid properties. Tethers have been studied extensively in experiments and described by theoretical models, but their molecular structure remains unknown due to their small diameters and dynamic nature. We used molecular dynamics simulations to obtain molecular-level insight into tether formation. Tethers were pulled from single-component lipid bilayers by application of an external force to a lipid patch along the bilayer normal or by lateral compression of a confined bilayer. Tether development under external force proceeded by viscoelastic protrusion followed by viscous lipid flow. Weak forces below a threshold value produced only a protrusion. Larger forces led to a crossover to tether elongation, which was linear at a constant force. Under lateral compression, tethers formed from undulations of unrestrained bilayer area. We characterized in detail the tether structure and its formation process, and obtained the material properties of the membrane. To our knowledge, these results provide the first molecular view of membrane tethers.

The partitioning of amino acid sidechains into the membrane is a key aspect of membrane protein folding. However, lipid bilayers exhibit rapidly changing physicochemical properties over their nanometer-scale thickness, which complicates understanding the thermodynamics and microscopic details of membrane partitioning. Recent data from diverse approaches, including protein insertion by the Sec translocon, folding of a small beta-barrel membrane protein and computer simulations of the exact distribution of a variety of small molecules and peptides, have joined older hydrophobicity scales for membrane protein prediction. We examine the correlations among the scales and find that they are remarkably correlated even though there are large differences in magnitude. We discuss the implications of these scales for understanding membrane protein structure and function.

Computer simulations suggest that the translocation of arginine through the hydrocarbon core of a lipid membrane proceeds by the formation of a water-filled defect that keeps the arginine molecule hydrated even at the center of the bilayer. We show here that adding additional arginine molecules into one of these water defects causes only a small change in free energy. The barrier for transferring multiple arginines through the membrane is approximately the same as for a single arginine and may even be lower depending on the exact geometry of the system. We discuss these results in the context of arginine-rich peptides such as antimicrobial and cell-penetrating peptides.

We simulated spontaneous fusion of small unilamellar vesicles mediated by lung surfactant protein B (SP-B) using the MARTINI force field. An SP-B monomer triggers fusion events by anchoring two vesicles and facilitating the formation of a lipid bridge between the proximal leaflets. Once a lipid bridge is formed, fusion proceeds via a previously described stalk – hemifusion diaphragm – pore-opening pathway. In the absence of protein, fusion of vesicles was not observed in either unbiased simulations or upon application of a restraining potential to maintain the vesicles in close proximity. The shape of SP-B appears to enable it to bind to two vesicles at once, forcing their proximity, and to facilitate the initial transfer of lipids to form a high-energy hemifusion intermediate. Our results may provide insight into more general mechanisms of protein-mediated membrane fusion, and a possible role of SP-B in the secretory pathway and transfer of lung surfactant to the gas exchange interface.