The transmembrane domain of Neu in a lipid bilayer: molecular dynamics simulations

The results of full-atom molecular dynamics simulations of the transmembrane domains (TMDs) of both native, and Glu⁶⁶⁴-mutant (either protonated or unprotonated) Neu in an explicit fully hydrated dimyristoylphosphatidylcholine (DMPC) lipid bilayer are presented. For the native TMD peptide, a 10.05 ns trajectory was collected, while for the mutant TMD peptides 5.05 ns trajectories were collected for each. The peptides in all three simulations display stable predominantly -helical hydrogen bonding throughout the trajectories. The only significant exception occurs near the C-terminal end of the native and unprotonated mutant TMDs just outside the level of the lipid headgroups, where -helical hydrogen bonding develops, introducing a kink in the backbone structure. However, there is no indication of the formation of a bulge within the hydrophobic region of either native or mutant peptides. Over the course of the simulation of the mutant peptide, it is found that a significant number of water molecules penetrate the hydrophobic region of the surrounding lipid molecules, effectively hydrating Glu⁶⁶⁴. If the energy cost of such water penetration is significant enough, this may be a factor in the enhanced dimerization affinity of Glu⁶⁶⁴-mutant Neu.     

Interaction between K+ channel gate modifier hanatoxin and lipid bilayer membranes analyzed by molecular dynamics simulation

Hanatoxin (HaTx) is an ellipsoidal-shaped peptide that binds to the voltage sensor of voltage-dependent channels. Of physicochemical interest, HaTx has a “ring” of charged residues around its periphery and a hydrophobic protrusion. It has previously been postulated that HaTx binds to and functions on the surface of membranes, but a recent fluorescent-quenching study has implied a fairly deep positioning of HaTx in the lipid bilayer membrane. We carried out numerous molecular dynamic simulations of HaTx1, a well-studied variant of HaTx, in fully hydrated phospholipid bilayers. The system reproduced the surface-binding mode of HaTx1, in which HaTx1 resided in the extracellular side (outer) of the water/membrane interface with the hydrophobic patch of HaTx1 facing the membrane interior. On the other hand, analyses with various parameter settings suggested that the surface-binding mode was unstable because of the substantial attractive electrostatic force between HaTx1 and the lipid head groups of the inner (opposite) leaflet. Compared with this electrostatic force, the energetic cost for membrane deformation involving meniscus formation appeared to be small. In an attempt to interpret the quenching data, we consider the possibility of dimpling (meniscus formation) that brings HaTx1 inward (only ~0.7–0.8 nm above the bilayer center), while accounting for the flexibility of both leaflets of the membrane and the long-range interaction between positively charged residues of the membrane-bound peptide and the polar head groups of the opposite leaflet of the membrane. It is suggested that molecular dynamics simulations taking into account the flexibility of the membrane surface is potentially useful in interpreting the fluorescence-quenching data.Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Hypoxia suppresses astrocyte glutamate transport independently of amyloid formation

Sustained hypoxia alters the expression of numerous proteins and predisposes individuals to Alzheimer’s disease (AD). We have previously shown that hypoxia in vitro alters Ca²⁺ homeostasis in astrocytes and promotes increased production of amyloid β peptides (Aβ) of AD. Indeed, alteration of Ca²⁺ homeostasis requires amyloid formation. Here, we show that electrogenic glutamate uptake by astrocytes is suppressed by hypoxia (1% O₂, 24 h) in a manner that is independent of amyloid β peptide formation. Thus, hypoxic suppression of glutamate uptake and expression levels of glutamate transporter proteins EAAT1 and EAAT2 were not mimicked by exogenous application of amyloid β peptide, or by prevention of endogenous amyloid peptide formation (using inhibitors of either β or γ secretase). Thus, dysfunction in glutamate homeostasis in hypoxic conditions is independent of Aβ production, but will likely contribute to neuronal damage and death associated with AD following hypoxic events.     

A study of comparative modelling,simulation and molecular dynamics of CXCR3 receptor with lipid bilayer

The G-coupled receptors seen on the cell surface are composites with a lipid bilayer. The chemokines are kind of G-coupled receptor which majorly involved in the activation and downstream signalling of the cell. In general, many G-coupled receptors lack their 3D structures which become a hurdle in the drug designing process. In this study, comparative modelling of the CXCR3 receptor was carried out, structure evaluation was done using various tools and softwares. Additionally, molecular dynamics and docking were performed to prove the structural quality and architecture. Interestingly, the studies like toggle switch mechanism, lipid dynamics, virtual screening were carried out to find the potent antagonist for the CXCR3 receptor. During virtual screening 14,303 similar molecules were retrieved among them only four compounds have an ability to interact with a crucial amino acid residue of an antagonist. Hence, these screened compounds can serve as a drug candidate for a CXCR3 receptor, but further in vitro and in vivo studies are ought to do to prove its same efficacy.     

Molecular dynamics simulations of local anesthetic articaine in a lipid bilayer

In order to investigate structural and dynamical properties of local anesthetic articaine in a model lipid bilayer, a series of molecular dynamics simulations have been performed. Simulations were carried out for neutral and charged (protonated) forms of articaine inserted in fully hydrated dimyristoylphosphatidylcholine (DMPC) lipid bilayer. For comparison purpose, a fully hydrated DMPC bilayer without articaine was also simulated. The length of each simulation was 200 ns. Various properties of the lipid bilayer systems in the presence of both charged and uncharged forms of articaine taken at two different concentrations have been examined: membrane area per lipid, mass density distributions, order parameters, radial distribution functions, head group tilt, diffusion coefficients, electrostatic potential, etc, and compared with results of previous simulations of DMPC bilayer in the presence of lidocaine. It was shown that addition of both charged and neutral forms of articaine causes increase of the dipole electrostatic potential in the membrane interior.     

Influence of a lipid bilayer on the conformational behavior of amphotericin B derivatives — A molecular dynamics study

Amphotericin B (AmB) is an effective but very toxic antifungal antibiotic. In our laboratory a series of AmB derivatives of improved selectivity of action was synthesized and tested. To understand molecular basis of this improvement, comparative conformational studies of amphotericin B and its two more selective derivatives were carried out in an aqueous solution and in a lipid membrane. These molecular simulation studies revealed that within a membrane environment the conformational behavior of the derivatives differs significantly from the one observed for the parent molecule. Possible reasons for such a difference are analyzed. Furthermore, we hypothesize that the observed conformational transition within the polar head of AmB derivatives may lead to destabilization of antibiotic-induced transmembrane channels. Consequently, the selective toxicity of the derivatives should increase as ergosterol-rich liquid-ordered domains are more rigid and conformationally ordered than their cholesterol-containing counterparts, and as such may better support less stable channel structure.     

Structure and dynamics of the full-length M1 muscarinic acetylcholine receptor studied by molecular dynamics simulations

The three-dimensional structure of full-length structure of the M₁ muscarinic receptor was obtained through the fragmental homology modeling procedure. A 10-ns molecular dynamics (MD) simulation of the protein imbedded in a lipid slab and surrounded by water molecules was further used to relax the model. It was found that the homology model corresponded to the conformation in the ground state, since no significant motions of the backbone of transmembrane domains were observed. Furthermore, the reliability of the model was validated by analyzing key inter-helical contacts, sidechain-sidechain interactions, the formation of stable aromatic microdomains (clusters) and the docking of acetylcholine to its binding site. Moreover, a few conserved interactions observed in the X-ray structure of rhodopsin, such as inter-helical sidechain-sidechain hydrogen bonds were accurately reproduced in the MD simulation. The coupling of ACh to its binding site was found to be dominated by π-cation and salt bridge interactions, while its conformational space was restrained through van der Waals and hydrogen bond interactions. In general, such features were in very good agreement with the available experimental as well as with theoretical data. Considering the above, the structural information obtained in this study can be used a starting point to investigate the activation mechanism of the receptor and the ability to develop selective agonists and allosteric modulators which could be used for the treatment of Alzheimer’s disease.     

The binding of resveratrol to monomer and fibril amyloid beta

As currently understood, Alzheimer’s disease (AD) is a chronic neurodegenerative disorder that is driven by the aggregation of amyloid beta (Aβ) protein. It has been shown that resveratrol (RES) may attenuate amyloid β peptide-induced toxicity, promote Aβ clearance and reduce senile plaques. However, it remains to be determined whether RES could interact directly with Aβ. The aim of the present study was to examine the direct binding of RES to monomer and fibril Aβ. Using surface plasmon resonance (SPR) and proton nuclear magnetic resonance (¹H NMR), our results identified the direct binding of RES to Aβ. The ability of RES to bind to both fibril and monomer Aβ(1–40 and 1–42) was further analyzed by SPR. The binding response of RES to fAβ(1–42) was higher than that to monomer Aβ(1–42), whereas the binding response of RES to fAβ(1–40) was lower than that to monomer Aβ(1–40). The K_D of RES for fibril Aβ(1–40 or 1–42) was higher than that for the corresponding monomer Aβ. Compared to the control compound Congo red (CR), the binding responses of RES to monomer Aβ(1–42) and Aβ(1–40) were stronger, but binding to fibril Aβ(1–42) was weaker, and the K_Ds of RES with both monomer and fibril Aβ(1–40) and Aβ(1–42) were higher than that of CR. When Aβ(1–40 or 1–42) was co-incubated with RES (50 μM), the thioflavin T fluorescence of the mixture was weakened, and the number and length of amyloid fibrils were decreased. Furthermore, the results of staining in consecutive brain slices from AD patients showed that RES (10⁻⁴ M) could stain senile plaques. These results indicated that RES could bind directly to Aβ in different states, which may provide new insight into the protective properties of RES against AD.     

Detection of Alzheimer's amyloid beta aggregation by capturing molecular trails of individual assemblies

Assembly of Amyloid beta (Aβ) peptides, in particular Aβ-42 is central to the formation of the amyloid plaques associated with neuro-pathologies such as Alzheimer’s disease (AD). Molecular assembly of individual Aβ-42 species was observed using a simple fluorescence microscope. From the molecular movements (aka Brownian motion) of the individual peptide assemblies, we calculated a temporal evolution of the hydrodynamic radius (R_H) of the peptide at physiological temperature and pH. The results clearly show a direct relationship between R_H of Aβ-42 and incubation period, corresponding to the previously reported peptide’s aggregation kinetics. The data correlates highly with in solution-based label-free electrochemical detection of the peptide’s aggregation, and Aβ-42 deposited on a solid surface and analysed using atomic force microscopy (AFM). To the best of our knowledge, this is the first analysis and characterisation of Aβ aggregation based on capturing molecular trails of individual assemblies. The technique enables both real-time observation and a semi-quantitative distribution profile of the various stages of Aβ assembly, at microM peptide concentration. Our method is a promising candidate for real-time observation and analysis of the effect of other pathologically-relevant molecules such as metal ions on pathways to Aβ oligomerisation and aggregation. The method is also a promising screening tool for AD therapeutics that target Aβ assembly.     

Conjugated double bonds in lipid bilayers: a molecular dynamics simulation study

Conjugated linoleic acids (CLA) are found naturally in dairy products. Two isomers of CLA, that differ only in the location of cis and trans double bonds, are found to have distinct and different biological effects. The cis 9 trans 11 (C9T11) isomer is attributed to have the anti-carcinogenic effects, while the trans 10 cis 12 (T10C12) isomer is believed to be responsible for the anti-obesity effects. Since dietary CLA are incorporated into membrane phospholipids, we have used Molecular Dynamics (MD) simulations to investigate the comparative effects of the two isomers on lipid bilayer structure. Specifically, simulations of phosphatidylcholine lipid bilayers in which the sn-2 chains contained one of the two isomers of CLA were performed. Force field parameters for the torsional potential of double bonds were obtained from ab initio calculations. From the MD trajectories we calculated and compared structural properties of the two lipid bilayers, including areas per molecule, density profiles, thickness of bilayers, tilt angle of tail chains, order parameters profiles, radial distribution function (RDF) and lateral pressure profiles. The main differences found between bilayers of the two CLA isomers, are (1) the order parameter profile for C9T11 has a dip in the middle of sn-2 chain while the profile for T10C12 has a deeper dip close to terminal of sn-2 chain, and (2) the lateral pressure profiles show differences between the two isomers. Our simulation results reveal localized physical structural differences between bilayers of the two CLA isomers that may contribute to different biological effects through differential interactions with membrane proteins or cholesterol.     

Molecular aspects of the interaction between amphotericin B and a phospholipid bilayer: molecular dynamics studies

Amphotericin B (AmB) is a polyene macrolide antibiotic used to treat systemic fungal infections. The molecular mechanism of AmB action is still only partly characterized. AmB interacts with cell-membrane components and forms membrane channels that eventually lead to cell death. The interaction between AmB and the membrane surface can be regarded as the first (presumably crucial) step on the way to channel formation. In this study molecular dynamics simulations were performed for an AmB–lipid bilayer model in order to characterize the molecular aspects of AmB–membrane interactions. The system studied contained a box of 200 dimyristoylphosphatidylcholine (DMPC) molecules, a single AmB molecule placed on the surface of the lipid bilayer and 8,065 water molecules. Two molecular dynamics simulations (NVT ensemble), each lasting 1 ns, were performed for the model studied. Two different programs, CHARMM and NAMD2, were used in order to test simulation conditions. The analysis of MD trajectories brought interesting information concerning interactions between polar groups of AmB and both DMPC and water molecules. Our studies show that AmB preferentially took a vertical position, perpendicular to the membrane surface, with no propensity to enter the membrane. Our finding may suggest that a single AmB molecule entering the membrane is very unlikely.Figure The figure presents the whole structure of the system simulated—starting point. AmB is presented as a space-filling model, DMPC molecules—green sticks, water molecules—red sticks     

A systematic molecular dynamics simulation study of temperature dependent bilayer structural properties

Although lipid force fields (FFs) used in molecular dynamics (MD) simulations have proved to be accurate, there has not been a systematic study on their accuracy over a range of temperatures. Motivated by the X-ray and neutron scattering measurements of common phosphatidylcholine (PC) bilayers (Ku?erka et al. BBA. 1808: 2761, 2011), the CHARMM36 (C36) FF accuracy is tested in this work with MD simulations of six common PC lipid bilayers over a wide range of temperatures. The calculated scattering form factors and deuterium order parameters from the C36 MD simulations agree well with the X-ray, neutron, and NMR experimental data. There is excellent agreement between MD simulations and experimental estimates for the surface area per lipid, bilayer thickness (D_B), hydrophobic thickness (D_C), and lipid volume (V_L). The only minor discrepancy between simulation and experiment is a measure of (D_B − D_HH) / 2 where D_HH is the distance between the maxima in the electron density profile along the bilayer normal. Additional MD simulations with pure water and heptane over a range of temperatures provide explanations of possible reasons causing the minor deviation. Overall, the C36 FF is accurate for use with liquid crystalline PC bilayers of varying chain types and over biologically relevant temperatures.     

Concentration effect of cimetidine with POPC bilayer: a molecular dynamics simulation study

All-atom molecular dynamics simulations have been performed on cimetidine in the presence of a palmitoyloleoylphosphatidylcholine (POPC) bilayer. The free energy profile of a single cimetidine molecule passing across POPC bilayer displays a minimum at the interface of bilayer and water. Ten cimetidine molecules were inserted into POPC bilayer to obtain an 8 mol % drug model, and molecular dynamics results showed that cimetidine molecules reside at the polar region of POPC bilayer with sulphur atoms directing to the hydrophobic region. By comparing the one drug model with 8 mol % drug model, one can see that the central barrier to cross the membrane increases while the free energy in bulk water decreases, indicating that the ability of cimetidine passing across the POPC bilayer weakens at increased concentration. In addition, the free energy minimum shifts closer to the hydrophobic core. Our results indicate that with the increased drug concentration, it is more difficult for cimetidine to enter and pass across POPC bilayer.     

Triazole-containing BODIPY dyes as novel fluorescent probes for soluble oligomers of amyloid Aβ1-42 peptide

A straightforward functionalization of BODIPY dyes via incorporation of a triazole moiety produced fluorescent dyes that were capable of distinguishing between secondary structure conformations of soluble oligomeric species of Aβ1-42 peptide. Small concentrations of the dyes, relative to Aβ1-42, provided up to an 8-fold and 35-fold fluorescence increase in the presence of the unordered and ordered, β-sheet-rich conformations of soluble Aβ1-42 oligomers, respectively. These triazole-containing dyes could prove to be useful probes for monitoring amyloid conformational transitions in vitro.     

Pregnenolone protects the PC-12 cell line against amyloid beta peptide toxicity but its sulfate ester does not

Pregnenolone (P), the main precursor of the steroids, and its sulfate ester, pregnenolone sulfate (PS), are the major neurosteroids produced in the neural tissue. Many neuroendocrinological studies stressed the neuroprotective role of neurosteroids although it has been suggested that the inhibition of P and PS synthesis can delay neuronal cell death. The potential roles of P and PS in vital neuronal functions and in amyloid beta peptide (Aβ) toxicity are not clearly identified. This work aims to investigate the effects of P and PS on cell viability and Aβ peptide toxicity in a concentration and exposure time-dependent manner in rat PC-12 cells. The cells were treated with 20 μM Aβ peptide 25-35 and variable concentrations of P and PS ranging from 0.5 μM to 100 μM. To examine the effects of steroid treatment on Aβ peptide toxicity, 0.5 μM (low) and 50 μM (high) neurosteroids were used. The cell viability and lactate dehydrogenase release of cells were evaluated after 24, 48 and 72 h. Morphological changes of cells were also examined. The treatment with higher than 1 μM concentrations of P and PS significantly decreased the cell viability comparing to untreated cells. At lower concentrations, P and PS had no toxic actions until 72 h. The Aβ treatment resulted in a significant decrease in cell viability comparing to untreated cells. P showed a dose-dependent protective effect against Aβ peptide in PC-12 cells. But its sulfate ester did not have the same effect on Aβ peptide toxicity, even it significantly decreased cell viability in Aβ-treated cells. Consequently, the discrepant effects of P and PS on Aβ peptide toxicity may provide insight on the pathogenesis of Alzheimer’s disease.     

Interactions of amphotericin B derivatives with lipid membranes—A molecular dynamics study

Amphotericin B (AmB) is a well-known polyene macrolide antibiotic used to treat systemic fungal infections. AmB targets more efficiently fungal than animal membranes. However, there are only minor differences in the mode of action of AmB against both types of membranes, which is a source of AmB toxicity. In this work, we analyzed interactions of two low toxic derivatives of AmB (SAmE and PAmE), synthesized in our laboratory, with lipid membranes. Molecular dynamics simulations of the lipid bilayers containing ergosterol (fungal cells) or cholesterol (animal cells) and the studied antibiotic molecules were performed to compare the structural and dynamic properties of AmB derivatives and the parent drug inside the membrane. A number of differences was found for AmB and its derivatives' behavior in cholesterol- and ergosterol-containing membranes. We found that PAmE and SAmE can penetrate deeper into the hydrophobic region of the membrane compared to AmB. Modification of the amino and carboxyl group of AmB also resulted in the conformational transition within the antibiotic's polar head. Wobbling dynamics differentiation, depending on the sterol present, was discovered for the AmB derivatives. These differences may be interpreted as molecular factors responsible for the improved selectivity observed macroscopically for the studied AmB derivatives.     

The influence of cholesterol on interactions and dynamics of ibuprofen in a lipid bilayer

In this work, molecular dynamics (MD) simulations with atomistic details were performed to examine the influence of the cholesterol on the interactions and the partitioning of the hydrophobic drug ibuprofen in a fully hydrated 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) bilayer. Analysis of MD simulations indicated that ibuprofen molecules prefer to be located in the hydrophobic acyl chain region of DMPC/cholesterol bilayers. This distribution decreases the lateral motion of lipid molecules. The presence of ibuprofen molecules in the bilayers with 0 and 25 mol% cholesterol increases the ordering of hydrocarbon tails of lipids whereas for the bilayers with 50 mol% cholesterol, ibuprofen molecules perturb the flexible chains of DMPC lipids which leads to the reduction of the acyl chain order parameter. The potential of the mean force (PMF) method was used to calculate the free energy profile for the transferring of an ibuprofen molecule from the bulk water into the DMPC/cholesterol membranes. The PMF studies indicated that the presence of 50 mol% cholesterol in the bilayers increases the free energy barrier and slows down the permeation of the ibuprofen drug across the DMPC bilayer. This can be due to the condensing and ordering effects of the cholesterol on the bilayer.     

Phosphorylation of amyloid precursor protein (APP) at Tyr687 regulates APP processing by alpha- and gamma-secretase

Abnormal proteolytic processing of amyloid precursor protein (APP) is a pathologic feature of Alzheimer’s disease. Recent studies have demonstrated that serine/threonine phosphorylation specifically at amino-acid residue Thr668 (APP695 numbering) regulates APP processing. In this study, we investigated the possibility that tyrosine phosphorylation of APP regulates APP processing. A tyrosine kinase inhibitor decreased expression of the C83 fragment which is a cleaved product of APP by α-secretase. By overexpressing APP mutant proteins, Tyr687 was found to be the major tyrosine kinase phosphorylation site. Expression of the C83 fragment was decreased in APPY687A-expressing cells relative to APP wild-type (APPWT)-expressing cells, which likely reflects the different cellular localization patterns of these two proteins. Expression of APP intracellular domain (AICD) which is a cleaved product of the C83 fragment by γ-secretase was decreased in C83Y687A-expressing cells. These results suggest that phosphorylation of APP at Tyr687 regulates APP processing by α- and γ-secretases, determining the expression level of AICD.     

Cold-active enzymes studied by comparative molecular dynamics simulation

Enzymes from cold-adapted species are significantly more active at low temperatures, even those close to zero Celsius, but the rationale of this adaptation is complex and relatively poorly understood. It is commonly stated that there is a relationship between the flexibility of an enzyme and its catalytic activity at low temperature. This paper gives the results of a study using molecular dynamics simulations performed for five pairs of enzymes, each pair comprising a cold-active enzyme plus its mesophilic or thermophilic counterpart. The enzyme pairs included α-amylase, citrate synthase, malate dehydrogenase, alkaline protease and xylanase. Numerous sites with elevated flexibility were observed in all enzymes; however, differences in flexibilities were not striking. Nevertheless, amino acid residues common in both enzymes of a pair (not present in insertions of a structure alignment) are generally more flexible in the cold-active enzymes. The further application of principle component analysis to the protein dynamics revealed that there are differences in the rate and/or extent of opening and closing of the active sites. The results indicate that protein dynamics play an important role in catalytic processes where structural rearrangements, such as those required for active site access by substrate, are involved. They also support the notion that cold adaptation may have evolved by selective changes in regions of enzyme structure rather than in global change to the whole protein. Figure Collective motions in C^α atoms of the active site of cold-active xylanase Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.