Simple but predictive protein models

The traditional approach to computational biophysics studies of molecular systems is brute force molecular dynamics simulations under the conditions of interest. The disadvantages of this approach are that the time and length scales that are accessible to computer simulations often do not reach biologically relevant scales. An alternative approach, which we call intuitive modeling, is hypothesis-driven and based on tailoring simplified protein models to the systems of interest. Using intuitive modeling, the length and time scales that can be achieved using simplified protein models exceed those of traditional molecular-dynamic simulations. Here, we describe several recent studies that signify the predictive power of simplified protein models within the intuitive-modeling approach.     

Extracting Curvature Preferences of Lipids Assembled in Flat Bilayers Shows Possible Kinetic Windows for Genesis of Bilayer Asymmetry and Domain Formation in Biological Membranes

Studies on the assembly of pure lipid components allow mechanistic insights toward understanding the structural and functional aspects of biological membranes. Molecular dynamic (MD) simulations on membrane systems provide molecular details on membrane dynamics that are difficult to obtain experimentally. A large number of MD studies have remained somewhat disconnected from a key concept of amphipathic assembly resulting in membrane structures—shape parameters of lipid molecules in those structures in aqueous environments. This is because most of the MD studies have been done on flat lipid membranes. With the above in view, we analyzed MD simulations of 26 pure lipid patches as a function of (1) lipid type(s) and (2) time of MD simulations along with 35–40 ns trajectories of five pure lipids. We report, for the first time, extraction of curvature preferences of lipids from MD simulations done on flat bilayers. Our results may lead to mechanistic insights into the possible origins of bilayer asymmetries and domain formation in biological membranes.     

Noninvasive cellular internalization of silver molecules by chitosan nanoneedles: a novel nanocarrier

We explore with molecular modeling, dynamics simulations, and a statistical model the ability of chitosan nanoneedles (CNNs) to be internalized into a model lipid bilayer as a function of their length, keeping in view of their applications in the field of biomedicine for advanced targeted drug delivery. In this study, we have computationally modeled and studied the structural geometry and the stability of CNNs formed by 4, 6, and 8 subunits. We reported the molecular surface analysis of the modeled CNNs along with molecular dynamic (MD) simulations studies toward revealing the noninvasive cellular internalization potential of these CNNs and a case study has been carried to study the ability of CNNs to translocate silver nanoparticles across membrane. The present results are strongly in support of further exploration of 8 subunits based CNNs for their application as target drug delivery vehicles. The hydrophilicity of the CNNs has been attributed as one of the key factors responsible for the internalization process. Moreover, our MD simulation studies marched the ability of CNNs to translocate silver nanoparticles through biological membrane in a similar manner that resembles cell-penetrating peptides.     

Multiscale simulations of heterogeneous model membranes

This review will focus on computer modeling aimed at providing insights into the existence, structure, size, and thermodynamic stability of localized domains in membranes of heterogeneous composition. Modeling the lateral organization within a membrane is problematic due to the relatively slow lateral diffusion rate for lipid molecules so that microsecond or longer time scales are needed to fully model the formation and stability of a raft in a membrane. Although atomistic simulations currently are not able to reach this scale, they can provide data on the intermolecular forces and correlations that are involved in lateral organization. These data can be used to define coarse grained models that are capable of predictions of lateral organization in membranes. In this paper, we review modeling efforts that use interaction data from MD simulations to construct coarse grained models for heterogeneous bilayers. In this review we will discuss MD simulations done with the aim of gaining the information needed to build accurate coarse-grained models. We will then review some of the coarse-graining work, emphasizing modeling that has resulted from or has a basis in atomistic simulations.     

Multiscale simulations of heterogeneous model membranes

This review will focus on computer modeling aimed at providing insights into the existence, structure, size, and thermodynamic stability of localized domains in membranes of heterogeneous composition. Modeling the lateral organization within a membrane is problematic due to the relatively slow lateral diffusion rate for lipid molecules so that microsecond or longer time scales are needed to fully model the formation and stability of a raft in a membrane. Although atomistic simulations currently are not able to reach this scale, they can provide data on the intermolecular forces and correlations that are involved in lateral organization. These data can be used to define coarse grained models that are capable of predictions of lateral organization in membranes. In this paper, we review modeling efforts that use interaction data from MD simulations to construct coarse grained models for heterogeneous bilayers. In this review we will discuss MD simulations done with the aim of gaining the information needed to build accurate coarse-grained models. We will then review some of the coarse-graining work, emphasizing modeling that has resulted from or has a basis in atomistic simulations.     

Regulation of channel function due to physical energetic coupling with a lipid bilayer

Regulation of membrane protein functions due to hydrophobic coupling with a lipid bilayer has been investigated. An energy formula describing interactions between lipid bilayer and integral ion channels with different structures, which is based on the screened Coulomb interaction approximation, has been developed. Here the interaction energy is represented as being due to charge-based interactions between channel and lipid bilayer. The hydrophobic bilayer thickness channel length mismatch is found to induce channel destabilization exponentially while negative lipid curvature linearly. Experimental parameters related to channel dynamics are consistent with theoretical predictions. To measure comparable energy parameters directly in the system and to elucidate the mechanism at an atomistic level we performed molecular dynamics (MD) simulations of the ion channel forming peptide–lipid complexes. MD simulations indicate that peptides and lipids experience electrostatic and van der Waals interactions for short period of time when found within each other’s proximity. The energies from these two interactions are found to be similar to the energies derived theoretically using the screened Coulomb and the van der Waals interactions between peptides (in ion channel) and lipids (in lipid bilayer) due to mainly their charge properties. The results of in silico MD studies taken together with experimental observable parameters and theoretical energetic predictions suggest that the peptides induce ion channels inside lipid membranes due to peptide–lipid physical interactions. This study provides a new insight helping better understand of the underlying mechanisms of membrane protein functions in cell membrane leading to important biological implications.     

Constructing the suitable initial configuration of the membrane-protein system in molecular dynamics simulations

A method for constructing the suitable initial configuration of the membrane-protein system for molecular dynamics (MD) simulations is presented. This method could provide some hydrated initial configurations and help us to determine the best surface area of the system by contracting the surface area and comparing the optimized lowest energy of the system by energy minimization. The gramicidin A (GA) channel in;the fully hydrated dimyristoylphosphatidylcholine (DMPC) bilayer was used as our model. Three configurations with different surface areas were selected and applied for one 400 ps and two 300 ps MD simulations at constant pressure and temperature. All simulations were fairly stable without any constraints. Through analysis of the MD trajectories we found that the system with the best surface area was more stable than the other two systems, whose sizes were changed in the simulations. Further analysis of the bilayer normal length and the order parameters of the lipid alkyl tails indicates that the system with the best surface area shows some characteristics of the L_α phase, while both the smaller and the larger size systems have distinct deviations from the L_α phase that we expect. This illustrates that the correct surface area and the suitable initial configuration have an important influence on the phase of the membrane in the MD simulation. In addition, by comparing the root mean square differences of GA relative to the initial structure and interaction energy between different components of the system for all three systems, we find that the state of the DMPC bilayer has exerted a significant influence on the structure of GA. All these results demonstrate the validity of our method for constructing the initial configuration of the membrane-protein system for MD simulations. Received: 10 September 1998 / Revised version: 19 March 1999 / Accepted: 19 March 1999     

Structure and dynamics of K channel pore-lining helices: a comparative simulation study

Isolated pore-lining helices derived from three types of K-channel have been analyzed in terms of their structural and dynamic features in nanosecond molecular dynamics (MD) simulations while spanning a lipid bilayer. The helices were 1) M1 and M2 from the bacterial channel KcsA (Streptomyces lividans), 2) S5 and S6 from the voltage-gated (Kv) channel Shaker (Drosophila melanogaster), and 3) M1 and M2 from the inward rectifier channel Kir6.2 (human). In the case of the Kv and Kir channels, for which x-ray structures are not known, both short and long models of each helix were considered. Each helix was incorporated into a lipid bilayer containing 127 palmitoyloleoylphosphatidylcholine molecules, which was solvated with approximately 4000 water molecules, yielding approximately 20, 000 atoms in each system. Nanosecond MD simulations were used to aid the definition of optimal lengths for the helix models from Kv and Kir. Thus the study corresponds to a total simulation time of 10 ns. The inner pore-lining helices (M2 in KcsA and Kir, S6 in Shaker) appear to be slightly more flexible than the outer pore-lining helices. In particular, the Pro-Val-Pro motif of S6 results in flexibility about a molecular hinge, as was suggested by previous in vacuo simulations (, Biopolymers. 39:503-515). Such flexibility may be related to gating in the corresponding intact channel protein molecules. Analysis of H-bonds revealed interactions with both water and lipid molecules in the water/bilayer interfacial region. Such H-bonding interactions may lock the helices in place in the bilayer during the folding of the channel protein (as is implicit in the two-stage model of membrane protein folding). Aromatic residues at the extremities of the helices underwent complex motions on both short (<10 ps) and long (>100 ps) time scales.     

Multiscale modeling of cellular actin filaments: from atomistic molecular to coarse-grained dynamics

In this article, we present a computational multiscale model for the characterization of subcellular proteins. The model is encoded inside a simulation tool that builds coarse-grained (CG) force fields from atomistic simulations. Equilibrium molecular dynamics simulations on an all-atom model of the actin filament are performed. Then, using the statistical distribution of the distances between pairs of selected groups of atoms at the output of the MD simulations, the force field is parameterized using the Boltzmann inversion approach. This CG force field is further used to characterize the dynamics of the protein via Brownian dynamics simulations. This combination of methods into a single computational tool flow enables the simulation of actin filaments with length up to 400 nm, extending the time and length scales compared to state-of-the-art approaches. Moreover, the proposed multiscale modeling approach allows to investigate the relationship between atomistic structure and changes on the overall dynamics and mechanics of the filament and can be easily (i) extended to the characterization of other subcellular structures and (ii) used to investigate the cellular effects of molecular alterations due to pathological conditions.     

Molecular Dynamics Simulations of Phospholipid Bilayers

Abstract

Molecular dynamics (MD) simulations at 37°C have been performed on three phospholipid bilayer systems composed of the lipids DLPE, DOPE, and DOPC. The model used included 24 explicit lipid molecules and explicit waters of solvation in the polar head group regions, together with constant-pressure periodic boundary conditions in three dimensions. Using this model, a MD simulation samples part of an infinite planar lipid bilayer. The lipid dynamics and packing behavior were characterized. Furthermore, using the results of the simulations, a number of diverse properties including bilayer structural parameters, hydrocarbon chain order parameters, dihedral conformations, electron density profile, hydration per lipid, and water distribution along the bilayer normal were calculated. Many of these properties are available for the three lipid systems chosen, making them well suited for evaluating the model and protocols used in these simulations by direct comparisons with experimental data. The calculated MD behavior, chain disorder, and lipid packing parameter, i.e. the ratio of the effective areas of hydrocarbon tails and head group per lipid (a_t/a_h), correctly predict the aggregation preferences of the three lipids observed experimentally at 37°C, namely: a gel bilayer for DLPE, a hexagonal tube for DOPE, and a liquid crystalline bilayer for DOPC. In addition, the model and conditions used in the MD simulations led to good agreement of the calculated properties of the bilayers with available experimental results, demonstrating the reliability of the simulations. The effects of the cis unsaturation in the hydrocarbon chains of DOPE and DOPC, compared to the fully saturated one in DLPE, as well as the effects of the different polar head groups of PC and PE with the same unsaturated chains on the lipid packing and bilayer structure have been investigated. The results of these studies indicate the ability of MD methods to provide molecular-level insights into the structure and dynamics of lipid assemblies.     

Coating of magnetic nanoparticles affects their interactions with model cell membranes

BackgroundThe use of functionalized iron oxide nanoparticles of various chemical properties and architectures offers a new promising direction in theranostic applications. The increasing applications of nanoparticles in medicine require that these engineered nanomaterials will contact human cells without damaging essential tissues. Thus, efficient delivery must be achieved, while minimizing cytotoxicity during passage through cell membranes to reach intracellular target compartments.MethodsDifferential Scanning Calorimetry (DSC), molecular modeling, and atomistic Molecular Dynamics (MD) simulations were performed for two magnetite nanoparticles coated with polyvinyl alcohol (PVA) and polyarabic acid (ARA) in order to assess their interactions with model DPPC membranes.ResultsDSC experiments showed that both nanoparticles interact strongly with DPPC lipid head groups, albeit to a different degree, which was further confirmed and quantified by MD simulations. The two systems were simulated, and dynamical and structural properties were monitored. A bimodal diffusion was observed for both nanoparticles, representing the diffusion in the water phase and in the proximity of the lipid bilayer. Nanoparticles did not enter the bilayer, but caused ordering of the head groups and reduced the area per lipid compared to the pure bilayer, with MAG-PVA interacting more strongly and being closer to the lipid bilayer.ConclusionsResults of DSC experiments and MD simulations were in excellent agreement. Our findings demonstrate that the external coating is a key factor that affects nanoparticle-membrane interactions. Magnetite nanoparticles coated with PVA and ARA did not destabilize the model membrane and can be considered promising platforms for biomedical applications.General significanceUnderstanding the physico-chemical interactions of different nanoparticle coatings in contact with model cell membranes is the first step for assessing toxic response and could lead to predictive models for estimating toxicity. DSC in combination with MD simulations is an effective strategy to assess physico-chemical interactions of coated nanoparticles with lipid bilayers.     

Setting up and running molecular dynamics simulations of membrane proteins

Molecular dynamics simulations have become a popular and powerful technique to study lipids and membrane proteins. We present some general questions and issues that should be considered prior to embarking on molecular dynamics simulation studies of membrane proteins and review common simulation methods. We suggest a practical approach to setting up and running simulations of membrane proteins, and introduce two new (related) methods to embed a protein in a lipid bilayer. Both methods rely on placing lipids and the protein(s) on a widely spaced grid and then 'shrinking' the grid until the bilayer with the protein has the desired density, with lipids neatly packed around the protein. When starting from a grid based on a single lipid structure, or several potentially different lipid structures (method 1), the bilayer will start well-packed but requires more equilibration. When starting from a pre-equilibrated bilayer, either pure or mixed, most of the structure of the bilayer stays intact, reducing equilibration time (method 2). The main advantages of these methods are that they minimize equilibration time and can be almost completely automated, nearly eliminating one time consuming step in MD simulations of membrane proteins.     

Molecular dynamics simulation study of interaction between a class IIa bacteriocin and its immunity protein

Molecular dynamics (MD) simulations of carnobacteriocin B2 (CbnB2), a structurally well-characterized class IIa bacteriocin, and its immunity protein (ImB2) in lipid bilayer environment have been conducted to explore the interaction between them. Six 30-ns simulations were conducted in DPPC or POPG bilayer systems. In these simulations, ImB2 was placed in the aqueous layer with different orientations facing CbnB2 to sample all the faces of ImB2. The MD results indicate that (i) while CbnB2 remained embedded in the bilayer, it tends to move toward the interface, and (ii) the presence of CbnB2 in the DPPC bilayer attracts ImB2 toward the bilayer. In one of the orientations in DPPC bilayer system (simulation 1), ImB2 penetrates the bilayer and interacts with CbnB2 by ion-pair interaction. At several instances toward the later half of the simulation (15-30 ns), ImB2 and CbnB2 were found to form salt-bridge between Arg95 of ImB2 and Glu24 of CbnB2. Simulation in POPG bilayer displayed strong interaction between the positively charged ImB2 and the negatively charged polar head groups of the POPG molecules at the lipid-water interface. However, ImB2 was not able to penetrate the bilayer thereby preventing any interaction between ImB2 and CbnB2.     

Coarse-grained modeling of the actin filament derived from atomistic-scale simulations

A coarse-grained (CG) procedure that incorporates the information obtained from all-atom molecular dynamics (MD) simulations is presented and applied to actin filaments (F-actin). This procedure matches the averaged values and fluctuations of the effective internal coordinates that are used to define a CG model to the values extracted from atomistic MD simulations. The fluctuations of effective internal coordinates in a CG model are computed via normal-mode analysis (NMA), and the computed fluctuations are matched with the atomistic MD results in a self-consistent manner. Each actin monomer (G-actin) is coarse-grained into four sites, and each site corresponds to one of the subdomains of G-actin. The potential energy of a CG G-actin contains three bonds, two angles, and one dihedral angle; effective harmonic bonds are used to describe the intermonomer interactions in a CG F-actin. The persistence length of a CG F-actin was found to be sensitive to the cut-off distance of assigning intermonomer bonds. Effective harmonic bonds for a monomer with its third nearest neighboring monomers are found to be necessary to reproduce the values of persistence length obtained from all-atom MD simulations. Compared to the elastic network model, incorporating the information of internal coordinate fluctuations enhances the accuracy and robustness for a CG model to describe the shapes of low-frequency vibrational modes. Combining the fluctuation-matching CG procedure and NMA, the achievable time- and length scales of modeling actin filaments can be greatly enhanced. In particular, a method is described to compute the force-extension curve using the CG model developed in this work and NMA. It was found that F-actin is easily buckled under compressive deformation, and a writhing mode is developed as a result. In addition to the bending and twisting modes, this novel writhing mode of F-actin could also play important roles in the interactions of F-actin with actin-binding proteins and in the force-generation process via polymerization.     

Steered molecular dynamics simulations reveal important mechanisms in reversible monoamine oxidase B inhibition

The monotopic membrane protein monoamine oxidase B (MAO B) is an important drug target for Parkinson's disease. In order to design more specific, and thereby more effective, inhibitors for this enzyme, it is necessary to determine what factors govern inhibitor specificity and the inhibitor binding process, including the roles of the lipid bilayer, the active site loop, and several key residues within the binding pocket. Atomistic molecular dynamics simulations of MAO B either embedded in a lipid bilayer or free in solution have been performed. The simulations suggest that the bilayer controls the availability of the active site cavity by regulating the degree of fluctuation in two key loops that form the greater part of the active site entrance (residues 85-110 and 155-165). In turn, the enzyme itself causes local thinning and a decrease in area per lipid of the surrounding bilayer environment. Additional MD simulations of MAO B in complex with seven different reversible inhibitors followed by nonequilibrium steered MD simulations of the inhibitor unbinding have also been performed. The simulations demonstrate that the average energy of interaction between inhibitor and MAO B residues during inhibitor egress is an effective indicator of inhibitor strength and is also useful for identifying key residues that govern inhibitor specificity. These data provide researchers with valuable tools for designing effective MAO B inhibitors as well as outline a method that can be translated to the study of other enzyme-inhibitor complexes.     

Numerical studies of the membrane fluorescent dyes dynamics in ground and excited states

Fluorescence methods are widely used in studies of biological and model membranes. The dynamics of membrane fluorescent markers in their ground and excited electronic states and correlations with their molecular surrounding within the fully hydrated phospholipid bilayer are still not well understood. In the present work, Quantum Mechanical (QM) calculations and Molecular Dynamics (MD) simulations are used to characterize location and interactions of two membrane polarity probes (Prodan; 6-propionyl-2-dimethylaminonaphthalene and its derivative Laurdan; 2-dimethylamino-6-lauroylnaphthalene) with the dioleoylphosphatidylcholine (DOPC) lipid bilayer model. MD simulations with fluorophores in ground and excited states are found to be a useful tool to analyze the fluorescent dye dynamics and their immediate vicinity. The results of QM calculations and MD simulations are in excellent agreement with available experimental data. The calculation shows that the two amphiphilic dyes initially placed in bulk water diffuse within 10 ns towards their final location in the lipid bilayer. Analysis of solvent relaxation process in the aqueous phase occurs on the picoseconds timescale whereas it takes nanoseconds at the lipid/water interface. Four different relaxation time constants, corresponding to different relaxation processes, where observed when the dyes were embedded into the membrane.     

Application of biasing‐potential replica‐exchange simulations for loop modeling and refinement of proteins in explicit solvent

Comparative or homology modeling of a target protein based on sequence similarity to a protein with known structure is widely used to provide structural models of proteins. Depending on the target‐template similarity these model structures may contain regions of limited structural accuracy. In principle, molecular dynamics (MD) simulations can be used to refine protein model structures and also to model loop regions that connect structurally conserved regions but it is limited by the currently accessible simulation time scales. A recently developed biasing potential replica exchange (BP‐REMD) method was used to refine loops and complete decoy protein structures at atomic resolution including explicit solvent. In standard REMD simulations several replicas of a system are run in parallel at different temperatures allowing exchanges at preset time intervals. In a BP‐REMD simulation replicas are controlled by various levels of a biasing potential to reduce the energy barriers associated with peptide backbone dihedral transitions. The method requires much fewer replicas for efficient sampling compared with T‐REMD. Application of the approach to several protein loops indicated improved conformational sampling of backbone dihedral angle of loop residues compared to conventional MD simulations. BP‐REMD refinement simulations on several test cases starting from decoy structures deviating significantly from the native structure resulted in final structures in much closer agreement with experiment compared to conventional MD simulations. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Multiscale modeling of lipids and lipid bilayers

A multiscale modeling approach is applied for simulations of lipids and lipid assemblies on mesoscale. First, molecular dynamics simulation of initially disordered system of lipid molecules in water within all-atomic model was carried out. On the next stage, structural data obtained from the molecular dynamics (MD) simulation were used to build a coarse-grained (ten sites) lipid model, with effective interaction potentials computed by the inverse Monte Carlo method. Finally, several simulations of the coarse-grained model on longer length- and time-scale were performed, both within Monte Carlo and molecular dynamics simulations: a periodical sample of lipid molecules ordered in bilayer, a free sheet of such bilayer without periodic boundary conditions, formation of vesicle from a plain membrane, process of self-assembly of lipids randomly dispersed in volume. It was shown that the coarse-grained model, developed exclusively from all-atomic simulation data, reproduces well all the basic features of lipids in water solution.     

Implementation of residue-level coarse-grained models in GENESIS for large-scale molecular dynamics simulations

Residue-level coarse-grained (CG) models have become one of the most popular tools in biomolecular simulations in the trade-off between modeling accuracy and computational efficiency. To investigate large-scale biological phenomena in molecular dynamics (MD) simulations with CG models, unified treatments of proteins and nucleic acids, as well as efficient parallel computations, are indispensable. In the GENESIS MD software, we implement several residue-level CG models, covering structure-based and context-based potentials for both well-folded biomolecules and intrinsically disordered regions. An amino acid residue in protein is represented as a single CG particle centered at the Cα atom position, while a nucleotide in RNA or DNA is modeled with three beads. Then, a single CG particle represents around ten heavy atoms in both proteins and nucleic acids. The input data in CG MD simulations are treated as GROMACS-style input files generated from a newly developed toolbox, GENESIS-CG-tool. To optimize the performance in CG MD simulations, we utilize multiple neighbor lists, each of which is attached to a different nonbonded interaction potential in the cell-linked list method. We found that random number generations for Gaussian distributions in the Langevin thermostat are one of the bottlenecks in CG MD simulations. Therefore, we parallelize the computations with message-passing-interface (MPI) to improve the performance on PC clusters or supercomputers. We simulate Herpes simplex virus (HSV) type 2 B-capsid and chromatin models containing more than 1,000 nucleosomes in GENESIS as examples of large-scale biomolecular simulations with residue-level CG models. This framework extends accessible spatial and temporal scales by multi-scale simulations to study biologically relevant phenomena, such as genome-scale chromatin folding or phase-separated membrane-less condensations.     

Dimerization of Aβ40 inside dipalmitoylphosphatidylcholine bilayer and its effect on bilayer integrity: Atomistic simulation at three temperatures

Amyloid-beta (Aβ) protein is related to Alzheimer disease (AD), and various experiments have shown that oligomers as small as dimers are cytotoxic. Recent studies have concluded that interactions of Aβ with neuronal cell membranes lead to disruption of membrane integrity and toxicity and they play a key role in the development of AD. Molecular dynamics (MD) simulations have been used to investigate Aβ in aqueous solution and membranes. We have previously studied monomeric Aβ40 embedded in dipalmitoylphosphatidylcholine (DPPC) membrane using MD simulations. Here, we explore interactions of two Aβ40 peptides in DPPC bilayer and its consequences on dimer distribution in a lipid bilayer and on the secondary structure of the peptides. We explored that N-terminals played an important role in dimeric Aβ peptide aggregations and Aβ-bilayer interactions, while C-terminals bound peptides to bilayer like anchors. We did not observe exiting of peptides in our simulations although we observed insertion of peptides into the core of bilayer in some of our simulations. So it seems that the presence of Aβ on membrane surface increases its aggregation rate, and as diffusion occurs in two dimensions, it can increase the probability of interpeptide interactions. We found that dimeric Aβ, like monomeric one, had the ability to cause structural destabilization of DPPC membrane, which in turn might ultimately lead to cell death in an in vivo system. This information could have important implications for understanding the affinity of Aβ oligomers (here dimer) for membranes and the mechanism of Aβ oligomer toxicity in AD.