State-dependent Lipid Interactions with the A2a Receptor Revealed by MD Simulations Using In Vivo-Mimetic Membranes.
Ontology highlight
ABSTRACT:
SUBMITTER: Song W
PROVIDER: S-EPMC7031699 | biostudies-literature | 2019 Feb
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Ceramides play a key modulatory role in many cellular processes, which results from their effect on the structure and dynamics of biological membranes. In this study, we investigate the influence of C16-ceramide (C16) on the biophysical properties of DMPC lipid bilayers using solid-state NMR and atomistic molecular dynamics (MD) simulations. MD simulations and NMR measurements were carried out for a pure DMPC bilayer and for a 20% DMPC-C16 mixture. Calculated key structural properties, namely area per lipid, chain order parameters, and mass density profiles, indicate that C16 has an ordering effect on the DMPC bilayer. Furthermore, the simulations predict that specific hydrogen-bonds form between DMPC and C16 molecules. Multi-nuclear solid-state NMR was used to verify these theoretical predictions. Chain order parameters extracted from (13)C(1)H dipole couplings were measured for both lipid and ceramide and follow the trend suggested by the MD simulations. Furthermore, (1)H-MAS NMR experiments showed a direct contact between ceramide and lipids.
| S-EPMC4137758 | biostudies-literature
Project description:Understanding the interactions between carbon nanoparticles (CNPs) and biological membranes is critically important for applications of CNPs in biomedicine and toxicology. Due to the complexity and diversity of the systems, most molecular simulation studies have focused on the interactions of CNPs and single component bilayers. In this work, we performed coarse-grained molecular dynamic (CGMD) simulations to investigate the behaviors of fullerenes in the presence of multiple lipid components in the plasma membranes with varying fullerene concentrations. Our results reveal that fullerenes can spontaneously penetrate the plasma membrane. Interestingly, fullerenes prefer to locate themselves in the region of the highly unsaturated lipids that are enriched in the inner leaflet of the plasma membrane. This causes fullerene aggregation even at low concentrations. When increasing fullerene concentrations, the fullerene clusters grow, and budding may emerge at the inner leaflet of the plasma membrane. Our findings suggest by tuning the lipid composition, fullerenes can be loaded deeply inside the plasma membrane, which can be useful for designing drug carrier liposomes. Moreover, the mechanisms of how fullerenes perturb multicomponent cell membranes and how they directly enter the cell are proposed. These insights can help to determine fullerene toxicity in living cells.
| S-EPMC9138838 | biostudies-literature
Project description:Anionic lipids and cholesterols (CHOL) are critical to the function of nicotinic acetylcholine receptors (nAChR). We investigated their interactions with an open- and closed-channel alpha4beta2 nAChR by over 10 ns molecular dynamics simulations in a ternary lipid mixture of 1-palmitoyl-2-oleoyl phosphatidylcholine (POPC), 1-palmitoyl-2-oleoyl phosphatidic acid (POPA), and CHOL with a ratio of 3:1:1 (Haddadian et al., J. Phys. Chem. B 2008, 112, 13981). On average there were 65 and 74 interfacial lipids around the closed- and open-channel alpha4beta2 nAChR, respectively, in the equilibrated simulation systems. In the open-channel system, 42% of the interfacial POPA had acyl chains partially inserted into intra- or intersubunit cavities, as compared to only 7% in the closed-channel alpha4beta2. No CHOL was found in cavities within single subunits, though some CHOL infiltrated into the gaps between subunits. Because of its smaller headgroup, POPA could access some nonannular sites where POPC could not easily reach due to steric exclusion. Furthermore, POPA acted not only as an acceptor for hydrogen bonding (H bonding) as POPC did, but also as a donor through its hydroxyl group for H bonding with the backbone of the protein. The charged headgroup of POPA allowed the lipid to form stable salt bridges with conserved Arg and Lys residues at the interfaces of the transmembrane (TM) and extracellular (EC) or intracellular (IC) domains of the alpha4beta2. A higher number of salt bridges and hydrogen bonds (H bonds) between POPA and the alpha4beta2 nAChR were found in the open system than in the closed system, suggesting a potential role of POPA in the equilibrium between different channel states. Most interfacial POPA molecules showed lower order parameters than the bulk POPA due to the mixed effect of gauche defects, hydrophobic mismatch, and the lipid orientations near the magic angle. These unique properties enable the interfacial POPA to achieve what POPC cannot with regard to specific interactions with the protein, thereby making POPA essential for the function of nAChR.
| S-EPMC2698126 | biostudies-literature
Project description:This work reports the effects of the bioflavinoids genistein and daidzein on lipid bilayers as determined by volume measurements, X-ray scattering, and molecular dynamics simulations. The experimental and simulated total molecular volumes were found to be in outstanding agreement with each other before the addition of genistein and daidzein and also after their addition. Both bioflavinoids inserted into the hydrocarbon region of both DOPC and diphytanoylPC near the carbonyls of the lipids and both decreased the bilayer thicknesses. The long axes of both bioflavinoids were oriented nearly parallel to the plane of the bilayer with their carbonyl groups preferentially pointed toward the proximal surface. A difference is that daidzein had a solubility limit of ∼0.14 mol fraction in DOPC (∼0.12 mol fraction in diphytanoylPC), whereas genistein was soluble at least to 0.20 mol fraction in both lipid membranes. Measurements of bending modulus K(C) and simulation results for area compressibility modulus K(A) indicate that both bioflavinoids soften bilayers.
| S-EPMC3320743 | biostudies-literature
Project description:Protein kinase C delta (PKC-δ) is an important signaling molecule in human cells that has both proapoptotic as well as antiapoptotic functions. These conflicting activities can be modulated by two classes of ligands, phorbol esters and bryostatins. Phorbol esters are known tumor promoters, while bryostatins have anti-cancer properties. This is despite both ligands binding to the C1b domain of PKC-δ (δC1b) with a similar affinity. The molecular mechanism behind this discrepancy in cellular effects remains unknown. Here, we have used molecular dynamics simulations to investigate the structure and intermolecular interactions of these ligands bound to δC1b with heterogeneous membranes. We observed clear interactions between the δC1b-phorbol complex and membrane cholesterol, primarily through the backbone amide of L250 and through the K256 side-chain amine. In contrast, the δC1b-bryostatin complex did not exhibit interactions with cholesterol. Topological maps of the membrane insertion depth of the δC1b-ligand complexes suggest that insertion depth can modulate δC1b interactions with cholesterol. The lack of cholesterol interactions suggests that bryostatin-bound δC1b may not readily translocate to cholesterol-rich domains within the plasma membrane, which could significantly alter the substrate specificity of PKC-δ compared to δC1b-phorbol complexes.
| S-EPMC10002858 | biostudies-literature
Project description:The interaction of ?-helical peptides with lipid bilayers is central to our understanding of the physicochemical principles of biological membrane organization and stability. Mutations that alter the position or orientation of an ?-helix within a membrane, or that change the probability that the ?-helix will insert into the membrane, can alter a range of membrane protein functions. We describe a comparative coarse-grained molecular dynamics simulation methodology, based on self-assembly of a lipid bilayer in the presence of an ?-helical peptide, which allows us to model membrane transmembrane helix insertion. We validate this methodology against available experimental data for synthetic model peptides (WALP23 and LS3). Simulation-based estimates of apparent free energies of insertion into a bilayer of cystic fibrosis transmembrane regulator-derived helices correlate well with published data for translocon-mediated insertion. Comparison of values of the apparent free energy of insertion from self-assembly simulations with those from coarse-grained molecular dynamics potentials of mean force for model peptides, and with translocon-mediated insertion of cystic fibrosis transmembrane regulator-derived peptides suggests a nonequilibrium model of helix insertion into bilayers.
| S-EPMC3077683 | biostudies-literature
Project description:Independent experimental and computational approaches show agreement concerning arginine/membrane interactions when a single arginine is introduced at selected positions within the membrane-spanning region of acetyl-GGALW(5)LALALAL(12)AL(14)ALALW(19)LAGA-ethanolamide, designated GWALP23. Peptide sequence isomers having Arg in position 12 or position 14 display markedly different behaviors, as deduced by both solid-state NMR experiments and coarse-grained molecular dynamics (CG-MD) simulations. With respect to the membrane normal of DOPC or DPPC lipid bilayer membranes, GWALP23-R14 shows one major state whose apparent average tilt is approximately 10 degrees greater than that of GWALP23. The presence of R14 furthermore induces bilayer thinning and peptide displacement to "lift" the charged guanidinium toward the bilayer surface. By contrast, GWALP23-R12 exhibits multiple states that are in slow exchange on the NMR time scale, with CG-MD simulations indicating two distinct positions with different screw rotation angles in the membrane, along with an increased tendency to exit the lipid bilayer.
| S-EPMC2888137 | biostudies-literature
Project description:Herein, we introduce a new strategy to estimate binding free energies using end-state molecular dynamics simulation trajectories. The method is adopted from linear interaction energy (LIE) and ANI-2x neural network potentials (machine learning) for the atomic simulation environment (ASE). It predicts the single-point interaction energies between ligand-protein and ligand-solvent pairs at the accuracy of the wb97x/6-31G* level for the conformational space that is sampled by molecular dynamics (MD) simulations. Our results on 54 protein-ligand complexes show that the method can be accurate and have a correlation of R = 0.87-0.88 to the experimental binding free energies, outperforming current end-state methods with reduced computational cost. The method also allows us to compare BFEs of ligands with different scaffolds. The code is available free of charge (documentation and test files) at https://github.com/otayfuroglu/deepQM.
| S-EPMC9472276 | biostudies-literature
Project description:Interactions of charged molecules with biomembranes regulate many of their biological activities, but their binding affinities to lipid bilayers are difficult to measure experimentally and model theoretically. Classical molecular dynamics (MD) simulations have the potential to capture the complex interactions determining how charged biomolecules interact with membranes, but systematic overbinding of sodium and calcium cations in standard MD simulations raises the question of how accurately force fields capture the interactions between lipid membranes and charged biomolecules. Here, we evaluate the binding of positively charged small molecules, etidocaine, and tetraphenylphosphonium to a phosphatidylcholine (POPC) lipid bilayer using the changes in lipid head-group order parameters. We observed that these molecules behave oppositely to calcium and sodium ions when binding to membranes: (i) their binding affinities are not overestimated by standard force field parameters, (ii) implicit inclusion of electronic polarizability increases their binding affinity, and (iii) they penetrate into the hydrophobic membrane core. Our results can be explained by distinct binding mechanisms of charged small molecules with hydrophobic moieties and monoatomic ions. The binding of the former is driven by hydrophobic effects, while the latter has direct electrostatic interactions with lipids. In addition to elucidating how different kinds of charged biomolecules bind to membranes, we deliver tools for further development of MD simulation parameters and methodology.
| S-EPMC9483918 | biostudies-literature
Project description:The dynamic and structural properties of membrane proteins are intimately affected by the lipid bilayer. One property of membrane proteins is uniaxial rotational diffusion, which depends on the membrane viscosity and thickness. This rotational diffusion is readily manifested in solid-state NMR spectra as characteristic line shapes and temperature-dependent line narrowing or broadening. We show here that this whole-body uniaxial diffusion is suppressed in lipid bilayers mimicking the composition of eukaryotic cell membranes, which are rich in cholesterol and sphingomyelin. We demonstrate this membrane-induced immobilization on the transmembrane peptide of the influenza A M2 (AM2-TM) proton channel protein. At physiological temperature, AM2-TM undergoes uniaxial diffusion faster than approximately 10(5) s(-1) in DLPC, DMPC, and POPC bilayers, but the motion is slowed by 2 orders of magnitude, to <10(3) s(-1), in a cholesterol-rich virus envelope-mimetic membrane ("viral membrane"). The immobilization is manifested as near rigid-limit (2)H quadrupolar couplings and (13)C-(1)H, (15)N-(1)H, and (13)C-(15)N dipolar couplings for all labeled residues. The immobilization suppresses intermediate time scale broadening of the NMR spectra, thus allowing high-sensitivity and high-resolution spectra to be measured at physiological temperature. The conformation of the protein in the viral membrane is more homogeneous than in model PC membranes, as evidenced by the narrow (15)N lines. The immobilization of the M2 helical bundle by the membrane composition change indicates the importance of studying membrane proteins in environments as native as possible. It also suggests that eukaryote-mimetic lipid membranes may greatly facilitate structure determination of membrane proteins by solid-state NMR.
| S-EPMC4082982 | biostudies-literature