1
|
Blanco-González A, Wurl A, Mendes Ferreira T, Piñeiro Á, Garcia-Fandino R. Simulating Bacterial Membrane Models at the Atomistic Level: A Force Field Comparison. J Chem Theory Comput 2024. [PMID: 39226695 DOI: 10.1021/acs.jctc.4c00204] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/05/2024]
Abstract
Molecular dynamics (MD) simulations are currently an indispensable tool to understand both the dynamic and nanoscale organization of cell membrane models. A large number of quantitative parameters can be extracted from these simulations, but their reliability is determined by the quality of the employed force field and the simulation parameters. Much of the work on parametrizing and optimizing force fields for biomembrane modeling has been focused on homogeneous bilayers with a single phospholipid type. However, these may not perform effectively or could even be unsuitable for lipid mixtures commonly employed in membrane models. This work aims to fill this gap by comparing MD simulation results of several bacterial membrane models using different force fields and simulation parameters, namely, CHARMM36, Slipids, and GROMOS-CKP. Furthermore, the hydrogen isotope exchange (HIE) method, combined with GROMOS-CKP (GROMOS-H2Q), was also tested to check for the impact of this acceleration strategy on the performance of the force field. A common set of simulation parameters was employed for all of the force fields in addition to those corresponding to the original parametrization of each of them. Furthermore, new experimental order parameter values determined from NMR of several lipid mixtures are also reported to compare them with those determined from MD simulations. Our results reveal that most of the calculated physical properties of bacterial membrane models from MD simulations are substantially force field and lipid composition dependent. Some lipid mixtures exhibit nearly ideal behaviors, while the interaction of different lipid types in other mixtures is highly synergistic. None of the employed force fields seem to be clearly superior to the other three, each having its own strengths and weaknesses. Slipids are notably effective at replicating the order parameters for all acyl chains, including those in lipid mixtures, but they offer the least accurate results for headgroup parameters. Conversely, CHARMM provides almost perfect estimates for the order parameters of the headgroups but tends to overestimate those of the lipid tails. The GROMOS parametrizations deliver reasonable order parameters for entire lipid molecules, including multicomponent bilayers, although they do not reach the accuracy of Slipids for tails or CHARMM for headgroups. Importantly, GROMOS-H2Q stands out for its computational efficiency, being at least 3 times faster than GROMOS, which is already faster than both CHARMM and Slipids. In turn, GROMOS-H2Q yields much higher compressibilities compared to all other parametrizations.
Collapse
Affiliation(s)
- Alexandre Blanco-González
- Facultad de Física, University of Santiago de Compostela (USC), 15782 Santiago de Compostela, Spain
- Singular Research Centre in Chemical Biology and Molecular Materials, (CIQUS), Organic Chemistry Department, University of Santiago de Compostela (USC), 15782 Santiago de Compostela, Spain
- MD.USE Innovations S.L, Edificio Emprendia, 15782 Santiago de Compostela, Spain
| | - Anika Wurl
- Institute of Physics, Faculty of Natural Sciences II, Betty-Heimann-Str. 7, 06120 Halle/Saale, Germany
| | - Tiago Mendes Ferreira
- Institute of Physics, Faculty of Natural Sciences II, Betty-Heimann-Str. 7, 06120 Halle/Saale, Germany
| | - Ángel Piñeiro
- Facultad de Física, University of Santiago de Compostela (USC), 15782 Santiago de Compostela, Spain
| | - Rebeca Garcia-Fandino
- Singular Research Centre in Chemical Biology and Molecular Materials, (CIQUS), Organic Chemistry Department, University of Santiago de Compostela (USC), 15782 Santiago de Compostela, Spain
| |
Collapse
|
2
|
Richardson JD, Van Lehn RC. Free Energy Analysis of Peptide-Induced Pore Formation in Lipid Membranes by Bridging Atomistic and Coarse-Grained Simulations. J Phys Chem B 2024. [PMID: 39207202 DOI: 10.1021/acs.jpcb.4c03276] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/04/2024]
Abstract
Antimicrobial peptides (AMPs) are attractive materials for combating the antimicrobial resistance crisis because they can kill target microbes by directly disrupting cell membranes. Although thousands of AMPs have been discovered, their molecular mechanisms of action are still poorly understood. One broad mechanism for membrane disruption is the formation of membrane-spanning hydrophilic pores which can be stabilized by AMPs. In this study, we use molecular dynamics simulations to investigate the thermodynamics of pore formation in model single-component lipid membranes in the presence of one of three AMPs: aurein 1.2, melittin and magainin 2. To overcome the general challenge of modeling long time scale membrane-related behaviors, including AMP binding, clustering, and pore formation, we develop a generalizable methodology for sampling AMP-induced pore formation. This approach involves the long equilibration of peptides around a pore created with a nucleation collective variable by performing coarse-grained simulations, then backmapping equilibrated AMP-membrane configurations to all-atom resolution. We then perform all-atom simulations to resolve free energy profiles for pore formation while accurately modeling the interplay of lipid-peptide-solvent interactions that dictate pore formation free energies. Using this approach, we quantify free energy barriers for pore formation without direct biases on peptides or whole lipids, allowing us to investigate mechanisms of pore formation for these 3 AMPs that are a consequence of unbiased peptide diffusion and clustering. Further analysis of simulation trajectories then relates variations in pore lining by AMPs, AMP-induced lipid disruptions, and salt bridges between AMPs to the observed pore formation free energies and corresponding mechanisms. This methodology and mechanistic analysis have the potential to generalize beyond the AMPs in this study to improve our understanding of pore formation by AMPs and related antimicrobial materials.
Collapse
Affiliation(s)
- Joshua D Richardson
- Department of Chemical and Biological Engineering, University of Wisconsin─Madison, Madison, Wisconsin 53706, United States
| | - Reid C Van Lehn
- Department of Chemical and Biological Engineering, University of Wisconsin─Madison, Madison, Wisconsin 53706, United States
- Department of Chemistry, University of Wisconsin─Madison, Madison, Wisconsin 53706, United States
| |
Collapse
|
3
|
Rangubpit W, Sungted S, Wong-Ekkabut J, Distaffen HE, Nilsson BL, Dias CL. Pore Formation by Amyloid-like Peptides: Effects of the Nonpolar-Polar Sequence Pattern. ACS Chem Neurosci 2024. [PMID: 39172951 DOI: 10.1021/acschemneuro.4c00333] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/24/2024] Open
Abstract
One of the mechanisms accounting for the toxicity of amyloid peptides in diseases like Alzheimer's and Parkinson's is the formation of pores on the plasma membrane of neurons. Here, we perform unbiased all-atom simulations of the full membrane damaging pathway, which includes adsorption, aggregation, and perforation of the lipid bilayer accounting for pore-like structures. Simulations are performed using four peptides made with the same amino acids. Differences in the nonpolar-polar sequence pattern of these peptides prompt them to adsorb into the membrane with the extended conformations oriented either parallel [peptide labeled F1, Ac-(FKFE)2-NH2], perpendicular (F4, Ac-FFFFKKEE-NH2), or with an intermediate orientation (F2, Ac-FFKKFFEE-NH2, and F3, Ac-FFFKFEKE-NH2) in regard to the membrane surface. At the water-lipid interface, only F1 fully self-assembles into β-sheets, and F2 peptides partially fold into an α-helical structure. The β-sheets of F1 emerge as electrostatic interactions attract neighboring peptides to intermediate distances where nonpolar side chains can interact within the dry core of the bilayer. This complex interplay between electrostatic and nonpolar interactions is not observed for the other peptides. Although β-sheets of F1 peptides are mostly parallel to the membrane, some of their edges penetrate deep inside the bilayer, dragging water molecules with them. This precedes pore formation, which starts with the flow of two water layers through the membrane that expand into a stable cylindrical pore delimited by polar faces of β-sheets spanning both leaflets of the bilayer.
Collapse
Affiliation(s)
- Warin Rangubpit
- Department of Physics, New Jersey Institute of Technology, Newark, New Jersey 07102-1982, United States
| | - Siwaporn Sungted
- Department of Physics, Faculty of Science, Kasetsart University, Bangkok 10900, Thailand
- Computational Biomodelling Laboratory for Agricultural Science and Technology (CBLAST), Faculty of Science, Kasetsart University, Bangkok 10900, Thailand
| | - Jirasak Wong-Ekkabut
- Department of Physics, Faculty of Science, Kasetsart University, Bangkok 10900, Thailand
- Computational Biomodelling Laboratory for Agricultural Science and Technology (CBLAST), Faculty of Science, Kasetsart University, Bangkok 10900, Thailand
| | - Hannah E Distaffen
- Department of Chemistry, University of Rochester, Rochester, New York 14627-0216, United States
| | - Bradley L Nilsson
- Department of Chemistry, University of Rochester, Rochester, New York 14627-0216, United States
- Materials Science Program, University of Rochester, Rochester, New York 14627-0166, United States
| | - Cristiano L Dias
- Department of Physics, New Jersey Institute of Technology, Newark, New Jersey 07102-1982, United States
| |
Collapse
|
4
|
Shamaprasad P, Nădăban A, Iacovella CR, Gooris GS, Bunge AL, Bouwstra JA, McCabe C. The phase behavior of skin-barrier lipids: A combined approach of experiments and simulations. Biophys J 2024:S0006-3495(24)00477-6. [PMID: 39030908 DOI: 10.1016/j.bpj.2024.07.018] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2024] [Revised: 04/28/2024] [Accepted: 07/15/2024] [Indexed: 07/22/2024] Open
Abstract
Skin barrier function is localized in its outermost layer, the stratum corneum (SC), which is comprised of corneocyte cells embedded in an extracellular lipid matrix containing ceramides (CERs), cholesterol (CHOL), and free fatty acids (FFAs). The unique structure and composition of this lipid matrix are important for skin barrier function. In this study, experiments and molecular dynamics simulation were combined to investigate the structural properties and phase behavior of mixtures containing nonhydroxy sphingosine CER (CER NS), CHOL, and FFA. X-ray scattering for mixtures with varying CHOL levels revealed the presence of the 5.4 nm short periodicity phase in the presence of CHOL. Bilayers in coarse-grained multilayer simulations of the same compositions contained domains with thicknesses of approximately 5.3 and 5.8 nm that are associated with elevated levels, respectively, of CER sphingosine chains with CHOL, and CER acyl chains with FFA chains. The prevalence of the thicker domain increased with decreasing CHOL content. This might correspond to a phase with ∼5.8 nm spacing observed by x-rays (other details unknown) in mixtures with lower CHOL content. Scissoring and stretching frequencies from Fourier transform infrared spectroscopy (FTIR) also indicate interaction between FFA and CER acyl chains and little interaction between CER acyl and CER sphingosine chains, which requires CER molecules to adopt a predominantly extended conformation. In the simulated systems, neighbor preferences of extended CER chains align more closely with the FTIR observations than those of CERs with hairpin ceramide chains. Both FTIR and atomistic simulations of reverse mapped multilayer membranes detect a hexagonal to fluid phase transition between 65 and 80°C. These results demonstrate the utility of a collaborative experimental and simulation effort in gaining a more comprehensive understanding of SC lipid membranes.
Collapse
Affiliation(s)
- Parashara Shamaprasad
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, Tennessee
| | - Andreea Nădăban
- Division of BioTherapeutics, Leiden Academic Centre for Drug Research, Leiden University, Leiden, the Netherlands
| | - Christopher R Iacovella
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, Tennessee
| | - Gerrit S Gooris
- Division of BioTherapeutics, Leiden Academic Centre for Drug Research, Leiden University, Leiden, the Netherlands
| | - Annette L Bunge
- Department of Chemical and Biological Engineering, Colorado School of Mines, Golden, Colorado
| | - Joke A Bouwstra
- Division of BioTherapeutics, Leiden Academic Centre for Drug Research, Leiden University, Leiden, the Netherlands
| | - Clare McCabe
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, Tennessee; School of Engineering and Physical Science, Heriot-Watt University, Edinburgh, United Kingdom.
| |
Collapse
|
5
|
Mallick S, Agmon N. Lateral diffusion of ions near membrane surface. Phys Chem Chem Phys 2024; 26:19433-19449. [PMID: 38973628 DOI: 10.1039/d3cp04112c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/09/2024]
Abstract
Biological membranes isolate living cells from their environment, while allowing selective molecular transport between the inner and outer realms. For example, Na+ and K+ permeability through ionic channels contributes to neural conduction. Whether the ionic currents arise directly from cations in the bulk, or from the interface, is currently unclear. There are only scant results concerning lateral diffusion of ions on aquated membrane surfaces (and strong belief that this occurs through binding to a diffusing lipid). We performed classical molecular dynamics (MD) simulations of monovalent ions, Na+, K+, and Cl-, near the surface of the zwitterionic palmitoyl-oleoyl-phosphatidylcholine (POPC) membrane. Realistic force-fields for lipids (Amber's Lipid17 and Lipid21) and water (TIP4P-Ew) are tested for the mass and charge densities and the electrostatic potential across the membrane. These calculations reveal that the chloride can bind to the choline moiety through an intervening water molecule by forming a CH⋯OH hydrogen bond, while cations bind to both the phosphatic and carbonyl oxygens of phosphatidylcholine moieties. Upon transitioning from the bulk to the interface, a cation sheds some of its hydration water, which are replaced by headgroup atoms. Notably, an interfacial cation can bind 1-4 headgroup atoms, which is a key to understanding its surface hopping mechanism. We find that cation binding to three headgroup atoms immobilizes it, while binding to four energizes it. Consequently, the lateral cation diffusion rate is only 15-25 times slower than in the bulk, and 4-5 times faster than lipid self-diffusion. K+ diffusion is notably more anomalous than Na+, switching from sub- to super-diffusion after about 2 ns.
Collapse
Affiliation(s)
- Subhasish Mallick
- The Fritz Haber Research Center, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel.
| | - Noam Agmon
- The Fritz Haber Research Center, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel.
| |
Collapse
|
6
|
Arango AS, Park H, Tajkhorshid E. Topological Learning Approach to Characterizing Biological Membranes. J Chem Inf Model 2024; 64:5242-5252. [PMID: 38912752 DOI: 10.1021/acs.jcim.4c00552] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/25/2024]
Abstract
Biological membranes play key roles in cellular compartmentalization, structure, and its signaling pathways. At varying temperatures, individual membrane lipids sample from different configurations, a process that frequently leads to higher-order phase behavior and phenomena. Here, we present a persistent homology (PH)-based method for quantifying the structural features of individual and bulk lipids, providing local and contextual information on lipid tail organization. Our method leverages the mathematical machinery of algebraic topology and machine learning to infer temperature-dependent structural information on lipids from static coordinates. To train our model, we generated multiple molecular dynamics trajectories of dipalmitoyl-phosphatidylcholine membranes at varying temperatures. A fingerprint was then constructed for each set of lipid coordinates by PH filtration, in which interaction spheres were grown around the lipid atoms while tracking their intersections. The sphere filtration formed a simplicial complex that captures enduring key topological features of the configuration landscape using homology, yielding persistence data. Following fingerprint extraction for physiologically relevant temperatures, the persistence data were used to train an attention-based neural network for assignment of effective temperature values to selected membrane regions. Our persistence homology-based method captures the local structural effects, via effective temperature, of lipids adjacent to other membrane constituents, e.g., sterols and proteins. This topological learning approach can predict lipid effective temperatures from static coordinates across multiple spatial resolutions. The tool, called MembTDA, can be accessed at https://github.com/hyunp2/Memb-TDA.
Collapse
Affiliation(s)
- Andres S Arango
- Theoretical and Computational Biophysics Group, NIH Resource Center for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, Department of Biochemistry, and Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Hyun Park
- Theoretical and Computational Biophysics Group, NIH Resource Center for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, Department of Biochemistry, and Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Emad Tajkhorshid
- Theoretical and Computational Biophysics Group, NIH Resource Center for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, Department of Biochemistry, and Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
7
|
Oliveira AC, Filipe HAL, Geraldes CF, Voth GA, Moreno MJ, Loura LMS. Interaction of MRI Contrast Agent [Gd(DOTA)] - with Lipid Membranes: A Molecular Dynamics Study. Inorg Chem 2024; 63:10897-10914. [PMID: 38795015 PMCID: PMC11186012 DOI: 10.1021/acs.inorgchem.4c00972] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2024] [Revised: 04/16/2024] [Accepted: 04/30/2024] [Indexed: 05/27/2024]
Abstract
Contrast agents are important imaging probes in clinical MRI, allowing the identification of anatomic changes that otherwise would not be possible. Intensive research on the development of new contrast agents is being made to image specific pathological markers or sense local biochemical changes. The most widely used MRI contrast agents are based on gadolinium(III) complexes. Due to their very high charge density, they have low permeability through tight biological barriers such as the blood-brain barrier, hampering their application in the diagnosis of neurological disorders. In this study, we explore the interaction between the widely used contrast agent [Gd(DOTA)]- (Dotarem) and POPC lipid bilayers by means of molecular dynamics simulations. This metal complex is a standard reference where several chemical modifications have been introduced to improve key properties such as bioavailability and targeting. The simulations unveil detailed insights into the agent's interaction with the lipid bilayer, offering perspectives beyond experimental methods. Various properties, including the impact on global and local bilayer properties, were analyzed. As expected, the results indicate a low partition coefficient (KP) and high permeation barrier for this reference compound. Nevertheless, favorable interactions are established with the membrane leading to moderately long residence times. While coordination of one inner-sphere water molecule is maintained for the membrane-associated chelate, the physical-chemical attributes of [Gd(DOTA)]- as a MRI contrast agent are affected. Namely, increases in the rotational correlation times and in the residence time of the inner-sphere water are observed, with the former expected to significantly increase the water proton relaxivity. This work establishes a reference framework for the use of simulations to guide the rational design of new contrast agents with improved relaxivity and bioavailability and for the development of liposome-based formulations for use as imaging probes or theranostic agents.
Collapse
Affiliation(s)
- Alexandre C. Oliveira
- Coimbra
Chemistry Centre, Institute of Molecular
Sciences (CQC-IMS), 3004-535 Coimbra, Portugal
- Department
of Chemistry, University of Coimbra, 3004-535 Coimbra, Portugal
| | - Hugo A. L. Filipe
- Coimbra
Chemistry Centre, Institute of Molecular
Sciences (CQC-IMS), 3004-535 Coimbra, Portugal
- CPIRN-IPG—Center
of Potential and Innovation of Natural Resources, Polytechnic Institute of Guarda, 6300-559 Guarda, Portugal
| | - Carlos F.G.C. Geraldes
- Coimbra
Chemistry Centre, Institute of Molecular
Sciences (CQC-IMS), 3004-535 Coimbra, Portugal
- Department
of Life Sciences, University of Coimbra, Calçada Martim de Freitas, 3000-393 Coimbra, Portugal
- CIBIT/ICNAS
- Instituto de Ciências Nucleares Aplicadas à Saúde, Pólo das Ciências
da Saúde, Azinhaga de Santa Comba, 3000-548 Coimbra, Portugal
| | - Gregory A. Voth
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, James Franck
Institute, and Institute for Biophysical Dynamics, University of Chicago, Chicago, Illinois 60637, United States
| | - Maria João Moreno
- Coimbra
Chemistry Centre, Institute of Molecular
Sciences (CQC-IMS), 3004-535 Coimbra, Portugal
- Department
of Chemistry, University of Coimbra, 3004-535 Coimbra, Portugal
- CNC−Center
for Neuroscience and Cell Biology, University
of Coimbra, 3004-517 Coimbra, Portugal
| | - Luís M. S. Loura
- Coimbra
Chemistry Centre, Institute of Molecular
Sciences (CQC-IMS), 3004-535 Coimbra, Portugal
- Faculty
of Pharmacy, University of Coimbra, 3000-548 Coimbra, Portugal
- CNC−Center
for Neuroscience and Cell Biology, University
of Coimbra, 3004-517 Coimbra, Portugal
| |
Collapse
|
8
|
Stachowicz-Kuśnierz A, Rychlik P, Korchowiec J, Korchowiec B. Insights into Transfer of Supramolecular Doxorubicin/Congo Red Aggregates through Phospholipid Membranes. Molecules 2024; 29:2567. [PMID: 38893446 PMCID: PMC11173802 DOI: 10.3390/molecules29112567] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2024] [Revised: 05/26/2024] [Accepted: 05/27/2024] [Indexed: 06/21/2024] Open
Abstract
Doxorubicin (DOX) is a commonly used chemotherapeutic drug, from the anthracycline class, which is genotoxic to neoplastic cells via a DNA intercalation mechanism. It is effective and universal; however, it also causes numerous side effects. The most serious of them are cardiotoxicity and a decrease in the number of myeloid cells. For this reason, targeted DOX delivery systems are desirable, since they would allow lowering the drug dose and therefore limiting systemic side effects. Recently, synthetic dyes, in particular Congo red (CR), have been proposed as possible DOX carriers. CR is a planar molecule, built of a central biphenyl moiety and two substituted naphthalene rings, connected with diazo bonds. In water, it forms elongated ribbon-shaped supramolecular structures, which are able to selectively interact with immune complexes. In our previous studies, we have shown that CR aggregates can intercalate DOX molecules. In this way, they preclude DOX precipitation in water solutions and increase its uptake by MCF7 breast cancer cells. In the present work, we further explore the interactions between DOX, CR, and their aggregates (CR/DOX) with phospholipid membranes. In addition to neutral molecules, the protonated doxorubicin form, DXP, is also studied. Molecular dynamics simulations are employed to study the transfer of CR, DOX, DXP, and their aggregates through POPC bilayers. Interactions of CR, DOX, and CR/DOX with model monolayers are studied with Langmuir trough measurements. This study shows that CR may support the transfer of doxorubicin molecules into the bilayer. Both electrostatic and van der Waals interactions with lipids are important in this respect. The former promote the initial stages of the insertion process, the latter keep guest molecules inside the bilayer.
Collapse
Affiliation(s)
- Anna Stachowicz-Kuśnierz
- Faculty of Chemistry, Jagiellonian University, Gronostajowa 2, 30-387 Krakow, Poland; (P.R.); (J.K.)
| | | | | | - Beata Korchowiec
- Faculty of Chemistry, Jagiellonian University, Gronostajowa 2, 30-387 Krakow, Poland; (P.R.); (J.K.)
| |
Collapse
|
9
|
Marzouq A, Morgenstein L, Huang-Zhu CA, Yudovich S, Atkins A, Grupi A, Van Lehn RC, Weiss S. Long-Chain Lipids Facilitate Insertion of Large Nanoparticles into Membranes of Small Unilamellar Vesicles. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2024; 40:10477-10485. [PMID: 38710504 PMCID: PMC11272290 DOI: 10.1021/acs.langmuir.3c03471] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/08/2024]
Abstract
Insertion of hydrophobic nanoparticles into phospholipid bilayers is limited to small particles that can incorporate into a hydrophobic membrane core between two lipid leaflets. Incorporation of nanoparticles above this size limit requires the development of challenging surface engineering methodologies. In principle, increasing the long-chain lipid component in the lipid mixture should facilitate incorporation of larger nanoparticles. Here, we explore the effect of incorporating very long phospholipids (C24:1) into small unilamellar vesicles on the membrane insertion efficiency of hydrophobic nanoparticles that are 5-11 nm in diameter. To this end, we improve an existing vesicle preparation protocol and utilized cryogenic electron microscopy imaging to examine the mode of interaction and evaluate the insertion efficiency of membrane-inserted nanoparticles. We also perform classical coarse-grained molecular dynamics simulations to identify changes in lipid membrane structural properties that may increase insertion efficiency. Our results indicate that long-chain lipids increase the insertion efficiency by preferentially accumulating near membrane-inserted nanoparticles to reduce the thermodynamically unfavorable disruption of the membrane.
Collapse
Affiliation(s)
- Adan Marzouq
- Department of Chemistry, Bar-Ilan University, Ramat-Gan, 52900, Israel
- Institute for Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan, 52900, Israel
| | - Lion Morgenstein
- Institute for Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan, 52900, Israel
- Department of Physics, Bar-Ilan University, Ramat-Gan, 52900, Israel
- Azrieli Faculty of Medicine, Bar-Ilan University, Safed, 1311502, Israel
| | - Carlos A. Huang-Zhu
- Department of Chemical and Biological Engineering, University of Wisconsin – Madison, Madison, WI 53706, USA
| | - Shimon Yudovich
- Institute for Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan, 52900, Israel
- Department of Physics, Bar-Ilan University, Ramat-Gan, 52900, Israel
- Department of Molecular and Cell Biology, University of California, Berkeley, CA 94720, USA
| | - Ayelet Atkins
- Institute for Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan, 52900, Israel
| | - Asaf Grupi
- Institute for Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan, 52900, Israel
- Department of Physics, Bar-Ilan University, Ramat-Gan, 52900, Israel
| | - Reid C. Van Lehn
- Department of Chemical and Biological Engineering, University of Wisconsin – Madison, Madison, WI 53706, USA
- Department of Chemistry, University of Wisconsin – Madison, Madison, WI 53706, USA
| | - Shimon Weiss
- Institute for Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan, 52900, Israel
- Department of Physics, Bar-Ilan University, Ramat-Gan, 52900, Israel
- Department of Chemistry and Biochemistry, University of California Los Angeles, Los Angeles, CA 90095, USA
- California NanoSystems Institute, University of California Los Angeles, Los Angeles, CA 90095, USA
| |
Collapse
|
10
|
Tripathy M, Srivastava A. Non-affine deformation analysis and 3D packing defects: A new way to probe membrane heterogeneity in molecular simulations. Methods Enzymol 2024; 701:541-577. [PMID: 39025582 DOI: 10.1016/bs.mie.2024.03.011] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/20/2024]
Abstract
Here, we discuss a new framework developed over the last 5 years in our group to probe nanoscale membrane heterogeneity. The framework is based on the idea of characterizing lateral heterogeneity through non-affine deformation (NAD) measurements, transverse heterogeneity through three dimensional (3D) lipid packing defects, and using these approaches to formalize the seemingly trivial correlation between lateral organization and lipid packing in biological membranes. We find that measurements from NAD analysis, a prescription which is borrowed from Physics of glasses and granular material, can faithfully distinguish between liquid-ordered and disordered phases in membranes at molecular length scales and, can also be used to identify phase boundaries with high precision. Concomitantly, 3D-packing defects can not only distinguish between the two co-existing fluid phases based on their molecular scale packing (or membrane free volume), but also provide a route to connect the membrane domains to their functionality, such as exploring the molecular origins of inter-leaflet domain registration and peptide partitioning. The correlation between lateral membrane order and transverse packing presents novel molecular design-level features that can explain functions such as protein/peptide partitioning and small-molecule permeation dynamics in complex and heterogeneous membranes with high-fidelity. The framework allows us to explore the nature of lateral organization and molecular packing as a manifestation of intricate molecular interactions among a chemically rich variety of lipids and other molecules in a membrane with complex membrane composition and asymmetry across leaflets.
Collapse
Affiliation(s)
- Madhusmita Tripathy
- Department of Chemistry, Technical University of Darmstadt, Darmstadt, Germany.
| | - Anand Srivastava
- Molecular Biophysics Unit, Indian Institute of Science Bangalore, Karnataka, India.
| |
Collapse
|
11
|
Pavlyuk E, Yungerman I, Bliznyuk A, Moskovitz Y. Studying the Effects of Dissolved Noble Gases and High Hydrostatic Pressure on the Spherical DOPC Bilayer Using Molecular Dynamic Simulations. MEMBRANES 2024; 14:89. [PMID: 38668117 PMCID: PMC11052037 DOI: 10.3390/membranes14040089] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/15/2024] [Revised: 04/04/2024] [Accepted: 04/10/2024] [Indexed: 04/28/2024]
Abstract
Fine-grained molecular dynamics simulations have been conducted to depict lipid objects enclosed in water and interacting with a series of noble gases dissolved in the medium. The simple point-charge (SPC) water system, featuring a boundary composed of 1,2-Dioleoyl-sn-glycero-3-phosphocholine (DOPC) molecules, maintained stability throughout the simulation under standard conditions. This allowed for the accurate modeling of the effects of hydrostatic pressure at an ambient pressure of 25 bar. The chosen pressure references the 240 m depth of seawater: the horizon frequently used by commercial divers, who comprise the primary patient population of the neurological complication of inert gas narcosis and the consequences of high-pressure neurological syndrome. To quantify and validate the neurological effects of noble gases and discriminate them from high hydrostatic pressure, we reduced the dissolved gas molar concentration to 1.5%, three times smaller than what we previously tested for the planar bilayer (3.5%). The nucleation and growth of xenon, argon and neon nanobubbles proved consistent with the data from the planar bilayer simulations. On the other hand, hyperbaric helium induces only a residual distorting effect on the liposome, with no significant condensed gas fraction observed within the hydrophobic core. The bubbles were distributed over a large volume-both in the bulk solvent and in the lipid phase-thereby causing substantial membrane distortion. This finding serves as evidence of the validity of the multisite distortion hypothesis for the neurological effect of inert gases at high pressure.
Collapse
Affiliation(s)
- Eugeny Pavlyuk
- Laboratory of Multi-Scale Mathematical Modeling, Ural Federal University, Ekaterinburg 620002, Russia
| | - Irena Yungerman
- Department of Chemical Engineering, Technion—Israel Technological Institute, Technion City, Haifa 3200003, Israel
| | - Alice Bliznyuk
- Ilse Katz Institute for Nanoscale Science and Technology (IKI), Ben Gurion University of the Negev, Beer Sheva 8410501, Israel
| | - Yevgeny Moskovitz
- Laboratory of Multi-Scale Mathematical Modeling, Ural Federal University, Ekaterinburg 620002, Russia
- Department of Chemical Engineering, Technion—Israel Technological Institute, Technion City, Haifa 3200003, Israel
| |
Collapse
|
12
|
Nguyen ATP, Weigle AT, Shukla D. Functional regulation of aquaporin dynamics by lipid bilayer composition. Nat Commun 2024; 15:1848. [PMID: 38418487 PMCID: PMC10901782 DOI: 10.1038/s41467-024-46027-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2023] [Accepted: 02/12/2024] [Indexed: 03/01/2024] Open
Abstract
With the diversity of lipid-protein interactions, any observed membrane protein dynamics or functions directly depend on the lipid bilayer selection. However, the implications of lipid bilayer choice are seldom considered unless characteristic lipid-protein interactions have been previously reported. Using molecular dynamics simulation, we characterize the effects of membrane embedding on plant aquaporin SoPIP2;1, which has no reported high-affinity lipid interactions. The regulatory impacts of a realistic lipid bilayer, and nine different homogeneous bilayers, on varying SoPIP2;1 dynamics are examined. We demonstrate that SoPIP2;1's structure, thermodynamics, kinetics, and water transport are altered as a function of each membrane construct's ensemble properties. Notably, the realistic bilayer provides stabilization of non-functional SoPIP2;1 metastable states. Hydrophobic mismatch and lipid order parameter calculations further explain how lipid ensemble properties manipulate SoPIP2;1 behavior. Our results illustrate the importance of careful bilayer selection when studying membrane proteins. To this end, we advise cautionary measures when performing membrane protein molecular dynamics simulations.
Collapse
Affiliation(s)
- Anh T P Nguyen
- Department of Chemical and Biomolecular Engineering, University of Illinois Urbana-Champaign, Urbana, IL, 61801, USA
| | - Austin T Weigle
- Department of Chemistry, University of Illinois Urbana-Champaign, Urbana, IL, 61801, USA
| | - Diwakar Shukla
- Department of Chemical and Biomolecular Engineering, University of Illinois Urbana-Champaign, Urbana, IL, 61801, USA.
- Center for Biophysics and Quantitative Biology, University of Illinois Urbana-Champaign, Urbana, IL, 61801, USA.
- Department of Bioengineering, University of Illinois Urbana-Champaign, Urbana, IL, 61801, USA.
- Department of Plant Biology, University of Illinois Urbana-Champaign, Urbana, IL, 61801, USA.
| |
Collapse
|
13
|
Xie J, Pink DL, Jayne Lawrence M, Lorenz CD. Digestion of lipid micelles leads to increased membrane permeability. NANOSCALE 2024; 16:2642-2653. [PMID: 38229565 DOI: 10.1039/d3nr05083a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/18/2024]
Abstract
Lipid-based drug carriers are an attractive option to solubilise poorly water soluble therapeutics. Previously, we reported that the digestion of a short tail PC lipid (2C6PC) by the PLA2 enzyme has a significant effect on the structure and stability of the micelles it forms. Here, we studied the interactions of micelles of varying composition representing various degrees of digestion with a model ordered (70 mol% DPPC & 30 mol% cholesterol) and disordered (100% DOPC) lipid membrane. Micelles of all compositions disassociated when interacting with the two different membranes. As the percentage of digestion products (C6FA and C6LYSO) in the micelle increased, the disassociation occurred more rapidly. The C6FA inserts preferentially into both membranes. We find that all micelle components increase the area per lipid, increase the disorder and decrease the thickness of the membranes, and the 2C6PC lipid molecules have the most significant impact. Additionally, there is an increase in permeation of water into the membrane that accompanies the insertion of C6FA into the DOPC membranes. We show that the natural digestion of lipid micelles result in molecular species that can enhance the permeability of lipid membranes that in turn result in an enhanced delivery of drugs.
Collapse
Affiliation(s)
- Jun Xie
- Biological & Soft Matter Research Group, Department of Physics, Faculty of Natural, Mathematical & Engineering Sciences, King's College London, London, UK.
| | - Demi L Pink
- Biological & Soft Matter Research Group, Department of Physics, Faculty of Natural, Mathematical & Engineering Sciences, King's College London, London, UK.
| | - M Jayne Lawrence
- Division of Pharmacy and Optometry, School of Health Sciences, Faculty of Biology, Medicine and Health, University of Manchester, Stopford Building, Oxford Road, Manchester, UK
| | - Christian D Lorenz
- Biological & Soft Matter Research Group, Department of Physics, Faculty of Natural, Mathematical & Engineering Sciences, King's College London, London, UK.
| |
Collapse
|
14
|
Moral R, Paul S. Exploring Cyclic Peptide Nanotube Stability Across Diverse Lipid Bilayers and Unveiling Water Transport Dynamics. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2024; 40:882-895. [PMID: 38134046 DOI: 10.1021/acs.langmuir.3c03030] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2023]
Abstract
Cyclic Peptide Nanotubes (CPNTs) have emerged as compelling candidates for various applications, particularly as nanochannels within lipid bilayers. In this study, the stability of two CPNTs, namely 8 × [(Cys-Gly-Met-Gly)2] and 8 × [(Gly-Leu)4], are comprehensively investigated across different lipid bilayers, including 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), a mixed model membrane (POPE/POPG), and a realistic yeast model membrane. The results demonstrate that both CPNTs maintain their tubular structures in all lipid bilayers, with [(Cys-Gly-Met-Gly)2] showing increased stability over an extended period in these lipid membranes. The insertion of CPNTs shows negligible impact on lipid bilayer properties, including area per lipid, volume per lipid, and bilayer thickness. The study demonstrates that the CPNT preserves its two-line water movement pattern within all the lipid membranes, reaffirming their potential as water channels. The MSD curves further reveal that the dynamics of water molecules inside the nanotube are similar for all the bilayer systems with minor differences that arise due to different lipid environments.
Collapse
Affiliation(s)
- Rimjhim Moral
- Department of Chemistry, Indian Institute of Technology, Guwahati, Assam 781039, India
| | - Sandip Paul
- Department of Chemistry, Indian Institute of Technology, Guwahati, Assam 781039, India
| |
Collapse
|
15
|
Rout M, Mishra S, Panda S, Dehury B, Pati S. Lipid and cholesterols modulate the dynamics of SARS-CoV-2 viral ion channel ORF3a and its pathogenic variants. Int J Biol Macromol 2024; 254:127986. [PMID: 37944718 DOI: 10.1016/j.ijbiomac.2023.127986] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2023] [Revised: 11/06/2023] [Accepted: 11/07/2023] [Indexed: 11/12/2023]
Abstract
SARS-CoV-2 accessory protein, ORF3a is a putative ion channel which immensely contributes to viral pathogenicity by modulating host immune responses and virus-host interactions. Relatively high expression of ORF3a in diseased individuals and implication with inflammasome activation, apoptosis and autophagy inhibition, ratifies as an effective target for developing vaccines and therapeutics. Herein, we present the elusive dynamics of ORF3a-dimeric state using all-atoms molecular dynamics (MD) simulations at μ-seconds scale in a heterogeneous lipid-mimetic system in multiple replicates. Additionally, we also explore the effect of non-synonymous pathogenic mutations on ORF3a ion channel activity and viral pathogenicity in different SARS-CoV-2 variants using various structure-based protein stability (ΔΔG) tools and computational saturation mutagenesis. Our study ascertains the role of phosphatidylcholines and cholesterol in modulating the structure of ORF3a, which perturbs the size and flexibility of the polar cavity that allows permeation of large cations. Discrete trend in ion channel pore radius and area per lipid arises the premise that presence of lipids might also affect the overall conformation of ORF3a. MD structural-ensembles, in some replicates rationalize the crucial role of TM2 in maintaining the native structure of ORF3a. We also infer that loss of structural stability primarily grounds for pathogenicity in more than half of the pathogenic variants of ORF3a. Overall, the effect of mutation on alteration of ion permeability of ORF3a, proposed in this study brings mechanistic insights into variant consequences on viral membrane proteins of SARS-CoV-2, which can be utilized for the development of novel therapeutics to treat COVID-19 and other coronavirus diseases.
Collapse
Affiliation(s)
- Madhusmita Rout
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Chandrasekharpur, Bhubaneswar 751023, Odisha, India
| | - Sarbani Mishra
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Chandrasekharpur, Bhubaneswar 751023, Odisha, India
| | - Sunita Panda
- Mycology Division, ICMR-Regional Medical Research Centre, Chandrasekharpur, Bhubaneswar 751023, Odisha, India
| | - Budheswar Dehury
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Chandrasekharpur, Bhubaneswar 751023, Odisha, India.
| | - Sanghamitra Pati
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Chandrasekharpur, Bhubaneswar 751023, Odisha, India.
| |
Collapse
|
16
|
Islam MM, Nawagamuwage SU, Parshin IV, Richard MC, Burin AL, Rubtsov IV. Probing the Hydrophobic Region of a Lipid Bilayer at Specific Depths Using Vibrational Spectroscopy. J Am Chem Soc 2023; 145:26363-26373. [PMID: 37982703 DOI: 10.1021/jacs.3c10178] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2023]
Abstract
A novel spectroscopic approach for studying the flexibility and mobility in the hydrophobic interior of lipid bilayers at specific depths is proposed. A set of test compounds featuring an azido moiety and a cyano or carboxylic acid moiety, connected by an alkyl chain of different lengths, was synthesized. FTIR data and molecular dynamics calculations indicated that the test compounds in a bilayer are oriented so that the cyano or carboxylic acid moiety is located in the lipid head-group region, while the azido group stays inside the bilayer at the depth determined by its alkyl chain length. We found that the asymmetric stretching mode of the azido group (νN3) can serve as a reporter of the membrane interior dynamics. FTIR and two-dimensional infrared (2DIR) studies were performed at different temperatures, ranging from 22 to 45 °C, covering the Lβ-Lα phase transition temperature of dipalmitoylphosphatidylcholine (∼41 °C). The width of the νN3 peak was found to be very sensitive to the phase transition and to the temperature in general. We introduced an order parameter, SN3, which characterizes restrictions to motion inside the bilayer. 2DIR spectra of νN3 showed different extents of inhomogeneity at different depths in the bilayer, with the smallest inhomogeneity in the middle of the leaflet. The spectral diffusion dynamics of the N3 peak was found to be dependent on the depth of the N3 group location in the bilayer. The obtained results enhance our understanding of the bilayer dynamics and can be extended to investigate membranes with more complex compositions.
Collapse
Affiliation(s)
- Md Muhaiminul Islam
- Department of Chemistry, Tulane University, New Orleans, Louisiana 70118, United States
| | | | - Igor V Parshin
- Department of Chemistry, Tulane University, New Orleans, Louisiana 70118, United States
| | - Margaret C Richard
- Department of Chemistry, Tulane University, New Orleans, Louisiana 70118, United States
| | - Alexander L Burin
- Department of Chemistry, Tulane University, New Orleans, Louisiana 70118, United States
| | - Igor V Rubtsov
- Department of Chemistry, Tulane University, New Orleans, Louisiana 70118, United States
| |
Collapse
|
17
|
Arango AS, Park H, Tajkhorshid E. Topological Learning Approach to Characterizing Biological Membranes. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.11.28.569053. [PMID: 38076911 PMCID: PMC10705453 DOI: 10.1101/2023.11.28.569053] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/20/2023]
Abstract
Biological membranes play key roles in cellular compartmentalization, structure, and its signaling pathways. At varying temperatures, individual membrane lipids sample from different configurations, a process that frequently leads to higher-order phase behavior and phenomena. Here we present a persistent homology-based method for quantifying the structural features of individual and bulk lipids, providing local and contextual information on lipid tail organization. Our method leverages the mathematical machinery of algebraic topology and machine learning to infer temperature-dependent structural information of lipids from static coordinates. To train our model, we generated multiple molecular dynamics trajectories of DPPC membranes at varying temperatures. A fingerprint was then constructed for each set of lipid coordinates by a persistent homology filtration, in which interactions spheres were grown around the lipid atoms while tracking their intersections. The sphere filtration formed a simplicial complex that captures enduring key topological features of the configuration landscape, using homology, yielding persistence data. Following fingerprint extraction for physiologically relevant temperatures, the persistence data were used to train an attention-based neural network for assignment of effective temperature values to selected membrane regions. Our persistence homology-based method captures the local structural effects, via effective temperature, of lipids adjacent to other membrane constituents, e.g. sterols and proteins. This topological learning approach can predict lipid effective temperatures from static coordinates across multiple spatial resolutions. The tool, called MembTDA, can be accessed at https://github.com/hyunp2/Memb-TDA.
Collapse
Affiliation(s)
- Andres S Arango
- Theoretical and Computational Biophysics Group, NIH Resource Center for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, Department of Biochemistry, and Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Hyun Park
- Theoretical and Computational Biophysics Group, NIH Resource Center for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, Department of Biochemistry, and Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Emad Tajkhorshid
- Theoretical and Computational Biophysics Group, NIH Resource Center for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, Department of Biochemistry, and Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
18
|
Divine-Ayela C, Perez F, Striolo A. Hop, Skip, and Jump: Hydrogen Molecular Transport through Amorphous Polyethylene Matrices Studied via Molecular Dynamics Simulations. Ind Eng Chem Res 2023; 62:19893-19906. [PMID: 38037624 PMCID: PMC10682999 DOI: 10.1021/acs.iecr.3c02213] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2023] [Revised: 10/12/2023] [Accepted: 10/12/2023] [Indexed: 12/02/2023]
Abstract
In the pursuit of advancing and diversifying energy technologies for a more sustainable future, the possibilities of hydrogen (H2) usage will broaden, as will our understanding of its containment materials. Polyethylene (PE) has a vast assortment of uses and applications, which are growing with demands for alternative energy possibilities. One use of PE liner is as a prime candidate for nonmetallic piping and pressurized type IV storage devices. Such applications require PE to effectively prevent H2 transport through containment systems. To study the molecular transport mechanism of hydrogen through polymeric barriers, a system containing hydrogen molecules absorbed within amorphous PE is modeled here using molecular dynamics simulations. The simulations are conducted within a range of temperatures that span the glass transition temperature of amorphous PE. The simulated PE displays bulk density, radius of gyration, and self-diffusion coefficient that are consistent with experimental data. The simulated trajectories are interrogated to study the movement of the guest gas molecules. The results show that the diffusion coefficients increase with temperature, as expected. However, the mobility of the PE chains is found to affect the mobility of absorbed H2 molecules to a much lower extent than it affects that of CH4 molecules because of the much smaller size of the former than of the latter guest. From a molecular perspective, a "hopping" mechanism is observed, according to which H2 molecules hop between one vacant free volume space to another within the polymer matrix, in combination with longer, straight, undisturbed "jumps" or "skips" along directions aligned with regions of ordered PE chains. This suggests that the orientation of the crystallites within the semicrystalline PE matrix affects the H2 containment. Implications of these findings toward PE usage as containment material are discussed.
Collapse
Affiliation(s)
- Candice Divine-Ayela
- Department
of Chemical Engineering, University College London, London WC1E 7JE, United Kingdom
| | - Felipe Perez
- School
of Sustainable Chemical, Biological and Materials Engineering, The University of Oklahoma, Norman, Oklahoma 73019, United States
| | - Alberto Striolo
- Department
of Chemical Engineering, University College London, London WC1E 7JE, United Kingdom
- School
of Sustainable Chemical, Biological and Materials Engineering, The University of Oklahoma, Norman, Oklahoma 73019, United States
| |
Collapse
|
19
|
Xu J, Karra V, Large DE, Auguste DT, Hung FR. Understanding the Mechanical Properties of Ultradeformable Liposomes Using Molecular Dynamics Simulations. J Phys Chem B 2023; 127:9496-9512. [PMID: 37879075 PMCID: PMC10641833 DOI: 10.1021/acs.jpcb.3c04386] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2023] [Revised: 09/19/2023] [Accepted: 10/04/2023] [Indexed: 10/27/2023]
Abstract
Improving drug delivery efficiency to solid tumor sites is a central challenge in anticancer therapeutic research. Our previous experimental study (Guo et al., Nat. Commun. 2018, 9, 130) showed that soft, elastic liposomes had increased uptake and accumulation in cancer cells and tumors in vitro and in vivo respectively, relative to rigid particles. As a first step toward understanding how liposomes' molecular structure and composition modulates their elasticity, we performed all-atom and coarse-grained classical molecular dynamics (MD) simulations of lipid bilayers formed by mixing a long-tailed unsaturated phospholipid with a short-tailed saturated lipid with the same headgroup. The former types of phospholipids considered were 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) and 1,2-dipalmitoleoyl-sn-glycero-3-phosphocholine (termed here DPMPC). The shorter saturated lipids examined were 1,2-diheptanoyl-sn-glycero-3-phosphocholine (DHPC), 1,2-didecanoyl-sn-glycero-3-phosphocholine (DDPC), 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC), and 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC). Several lipid concentrations and surface tensions were considered. Our results show that DOPC or DPMPC systems having 25-35 mol % of the shortest lipids DHPC or DDPC are the least rigid, having area compressibility moduli KA that are ∼10% smaller than the values observed in pure DOPC or DPMPC bilayers. These results agree with experimental measurements of the stretching modulus and lysis tension in liposomes with the same compositions. These mixed systems also have lower areas per lipid and form more uneven x-y interfaces with water, the tails of both primary and secondary lipids are more disordered, and the terminal methyl groups in the tails of the long lipid DOPC or DPMPC wriggle more in the vertical direction, compared to pure DOPC or DPMPC bilayers or their mixtures with the longer saturated lipid DLPC or DMPC. These observations confirm our hypothesis that adding increasing concentrations of the short unsaturated lipid DHPC or DDPC to DOPC or DPMPC bilayers alters lipid packing and thus makes the resulting liposomes more elastic and less rigid. No formation of lipid nanodomains was noted in our simulations, and no clear trends were observed in the lateral diffusivities of the lipids as the concentration, type of secondary lipid, and surface tension were varied.
Collapse
Affiliation(s)
- Jiaming Xu
- Department
of Chemical Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| | - Vyshnavi Karra
- Department
of Chemical Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| | - Danielle E. Large
- Department
of Bioengineering, Northeastern University, Boston, Massachusetts 02115, United States
| | - Debra T. Auguste
- Department
of Chemical Engineering, Northeastern University, Boston, Massachusetts 02115, United States
- Department
of Bioengineering, Northeastern University, Boston, Massachusetts 02115, United States
| | - Francisco R. Hung
- Department
of Chemical Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| |
Collapse
|
20
|
Abdelmessih R, Xu J, Hung FR, Auguste DT. Integration of an LPAR1 Antagonist into Liposomes Enhances Their Internalization and Tumor Accumulation in an Animal Model of Human Metastatic Breast Cancer. Mol Pharm 2023; 20:5500-5514. [PMID: 37844135 PMCID: PMC10631474 DOI: 10.1021/acs.molpharmaceut.3c00348] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2023] [Revised: 09/29/2023] [Accepted: 10/02/2023] [Indexed: 10/18/2023]
Abstract
Lysophosphatidic acid receptor 1 (LPAR1) is elevated in breast cancer. The deregulation of LPAR1, including the function and level of expression, is linked to cancer initiation, progression, and metastasis. LPAR1 antagonists, AM095 or Ki16425, may be effective therapeutic molecules, yet their limited water solubility hinders in vivo delivery. In this study, we report on the synthesis of two liposomal formulations incorporating AM095 or Ki16425, embedded within the lipid bilayer, as targeted nanocarriers for metastatic breast cancer (MBC). The data show that the Ki16425 liposomal formulation exhibited a 50% increase in internalization by MBC mouse epithelial cells (4T1) and a 100% increase in tumor accumulation in a mouse model of MBC compared with that of a blank liposomal formulation (control). At the same time, normal mouse epithelial cells (EpH-4Ev) internalized the Ki16425 liposomal formulation 25% lesser than the control formulation. Molecular dynamics simulations show that the integration of AM095 or Ki16425 modified the physical and mechanical properties of the lipid bilayer, making it more flexible in these liposomal formulations compared with liposomes without drug. The incorporation of an LPAR1 antagonist within a liposomal drug delivery system represents a viable therapeutic approach for targeting the LPA-LPAR1 axis, which may hinder the progression of MBC.
Collapse
Affiliation(s)
- Rudolf
G. Abdelmessih
- Department of Chemical Engineering, Northeastern University, 360 Huntington Avenue, Boston, Massachusetts 02115, United States
| | - Jiaming Xu
- Department of Chemical Engineering, Northeastern University, 360 Huntington Avenue, Boston, Massachusetts 02115, United States
| | - Francisco R. Hung
- Department of Chemical Engineering, Northeastern University, 360 Huntington Avenue, Boston, Massachusetts 02115, United States
| | - Debra T. Auguste
- Department of Chemical Engineering, Northeastern University, 360 Huntington Avenue, Boston, Massachusetts 02115, United States
| |
Collapse
|
21
|
Kumar J, Chng CP, Huang C. Hydrophobic Matching Dictates over the Linear Rule of Mixtures in Binary Lipid Membranes. J Phys Chem B 2023; 127:7946-7954. [PMID: 37674349 DOI: 10.1021/acs.jpcb.3c04502] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/08/2023]
Abstract
Biological membranes feature heterogeneous mixtures of lipids with different head and tail characteristics. Their biophysical properties are dictated by the intimate interaction among different constituent lipids. Previous studies suggest that the membrane area-per-lipid (APL) deviates from the linear rule of mixtures (LRM) for binary lipid membranes, but the underlying mechanism remains elusive. Our molecular dynamics (MD) simulations of binary lipid membranes consisting of lipids with different tail characteristics reveal a competitive mechanism whereby lipids tend to deform each other to minimize the hydrophobic mismatch between their tails. Depending on the relative tail lengths and saturation levels, this may result in an either positive or negative deviation of APL from the LRM. As lipid packing plays an essential role in membrane fusion and peptide-membrane binding, our findings may help guide the selection of lipids for the effective rational design of nanoliposomes and membrane-targeting peptides.
Collapse
Affiliation(s)
- Jatin Kumar
- School of Mechanical and Aerospace Engineering, Nanyang Technological University, Singapore 639798, Republic of Singapore
| | - Choon-Peng Chng
- School of Mechanical and Aerospace Engineering, Nanyang Technological University, Singapore 639798, Republic of Singapore
| | - Changjin Huang
- School of Mechanical and Aerospace Engineering, Nanyang Technological University, Singapore 639798, Republic of Singapore
| |
Collapse
|
22
|
Campbell O, Monje-Galvan V. Lipid composition modulates interactions of p7 viroporin during membrane insertion. J Struct Biol 2023; 215:108013. [PMID: 37586469 DOI: 10.1016/j.jsb.2023.108013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2023] [Revised: 07/05/2023] [Accepted: 08/09/2023] [Indexed: 08/18/2023]
Abstract
Viral proteins interact with lipid membranes during various stages in the viral life cycle to propagate infection. p7 is an ion channel forming protein of Hepatitis C virus (HCV) that participates in viral assembly. Studies show that it has close ties to lipid metabolism in the cell and anionic phosphatidylserine (PS) lipids are suggested to be key for its permeabilizing function, but the mechanism of its interaction with the lipid environment is largely unknown. To begin unraveling the molecular processes of the protein, we evaluated the impact of lipid environment on the binding and insertion mechanism of p7 prior to channel formation and viral assembly using molecular dynamics simulations. It is seen that p7 is sensitive to its lipid environment and results in different remodeling patterns in membranes. Helix 1 (H1) is especially important for peptide insertion, with deeper entry taking place when the membrane contains phosphatidylserine (PS). Helix 2 (H2) and the adjacent loop connecting to Helix 3 (H3) prompts recruitment of phosphatidylethanolamine (PE) lipids to the protein binding site in membrane models with lower surface charge. This work provides perspectives on the interplay between protein-lipid dynamics and membrane composition, and insights on membrane reorganization in mechanisms of disease.
Collapse
Affiliation(s)
- Oluwatoyin Campbell
- Department of Chemical and Biological Engineering, School of Engineering and Applied Sciences, University at Buffalo, Buffalo, NY, USA
| | - Viviana Monje-Galvan
- Department of Chemical and Biological Engineering, School of Engineering and Applied Sciences, University at Buffalo, Buffalo, NY, USA.
| |
Collapse
|
23
|
Nguyen ATP, Weigle AT, Shukla D. Functional Regulation of Aquaporin Dynamics by Lipid Bilayer Composition. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.07.20.549977. [PMID: 37502896 PMCID: PMC10370204 DOI: 10.1101/2023.07.20.549977] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/29/2023]
Abstract
With the diversity of lipid-protein interactions, any observed membrane protein dynamics or functions directly depend on the lipid bilayer selection. However, the implications of lipid bilayer choice are seldom considered unless characteristic lipid-protein interactions have been previously reported. Using molecular dynamics simulation, we characterize the effects of membrane embedding on plant aquaporin SoPIP2;1, which has no reported high-affinity lipid interactions. The regulatory impacts of a realistic lipid bilayer, and nine different homogeneous bilayers, on varying SoPIP2;1 dynamics were examined. We demonstrate that SoPIP2;1s structure, thermodynamics, kinetics, and water transport are altered as a function of each membrane construct's ensemble properties. Notably, the realistic bilayer provides stabilization of non-functional SoPIP2;1 metastable states. Hydrophobic mismatch and lipid order parameter calculations further explain how lipid ensemble properties manipulate SoPIP2;1 behavior. Our results illustrate the importance of careful bilayer selection when studying membrane proteins. To this end, we advise cautionary measures when performing membrane protein molecular dynamics simulations.
Collapse
Affiliation(s)
- Anh T P Nguyen
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, IL 61801
| | - Austin T Weigle
- Department of Chemistry, University of Illinois at Urbana-Champaign, IL 61801
| | - Diwakar Shukla
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, IL 61801
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, IL 61801
- Department of Bioengineering, University of Illinois at Urbana-Champaign, IL 61801
- Department of Plant Biology, University of Illinois at Urbana-Champaign, IL 61801
| |
Collapse
|
24
|
Yadav A, Kelich P, Kallmyer N, Reuel NF, Vuković L. Characterizing the Interactions of Cell-Membrane-Disrupting Peptides with Lipid-Functionalized Single-Walled Carbon Nanotubes. ACS APPLIED MATERIALS & INTERFACES 2023; 15:24084-24096. [PMID: 37184257 DOI: 10.1021/acsami.3c01217] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/16/2023]
Abstract
Lipid-functionalized single-walled carbon nanotubes (SWNTs) have garnered significant interest for their potential use in a wide range of biomedical applications. In this work, we used molecular dynamics simulations to study the equilibrium properties of SWNTs surrounded by the phosphatidylcholine (POPC) corona phase and their interactions with three cell membrane disruptor peptides: colistin, TAT peptide, and crotamine-derived peptide. Our results show that SWNTs favor asymmetrical positioning within the POPC corona, so that one side of the SWNT, covered by the thinnest part of the corona, comes in contact with charged and polar functional groups of POPC and water. We also observed that colistin and TAT insert deeply into the POPC corona, while crotamine-derived peptide only adsorbs to the corona surface. In separate simulations, we show that three examined peptides exhibit similar insertion and adsorption behaviors when interacting with POPC bilayers, confirming that peptide-induced perturbations to POPC in conjugates and bilayers are similar in nature and magnitude. Furthermore, we observed correlations between the peptide-induced structural perturbations and the near-infrared emission of the lipid-functionalized SWNTs, which suggest that the optical signal of the conjugates transduces the morphological changes in the lipid corona. Overall, our findings indicate that lipid-functionalized SWNTs could serve as simplified cell membrane model systems for prescreening of new antimicrobial compounds that disrupt cell membranes.
Collapse
Affiliation(s)
- Anju Yadav
- Department of Chemistry and Biochemistry, University of Texas at El Paso, El Paso, Texas 79968, United States of America
| | - Payam Kelich
- Department of Chemistry and Biochemistry, University of Texas at El Paso, El Paso, Texas 79968, United States of America
| | | | - Nigel F Reuel
- Department of Chemical and Biological Engineering, Iowa State University, Ames, Iowa 50011, United States of America
| | - Lela Vuković
- Department of Chemistry and Biochemistry, University of Texas at El Paso, El Paso, Texas 79968, United States of America
| |
Collapse
|
25
|
Sepulveda RV, Sbarbaro C, Opazo MC, Duarte Y, González-Nilo F, Aguayo D. Insights into Early Steps of Decanoic Acid Self-Assemblies under Prebiotic Temperatures Using Molecular Dynamics Simulations. MEMBRANES 2023; 13:membranes13050469. [PMID: 37233530 DOI: 10.3390/membranes13050469] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/11/2023] [Revised: 04/18/2023] [Accepted: 04/25/2023] [Indexed: 05/27/2023]
Abstract
The origin of life possibly required processes in confined systems that facilitated simple chemical reactions and other more complex reactions impossible to achieve under the condition of infinite dilution. In this context, the self-assembly of micelles or vesicles derived from prebiotic amphiphilic molecules is a cornerstone in the chemical evolution pathway. A prime example of these building blocks is decanoic acid, a short-chain fatty acid capable of self-assembling under ambient conditions. This study explored a simplified system made of decanoic acids under temperatures ranging from 0 °C to 110 °C to replicate prebiotic conditions. The study revealed the first point of aggregation of decanoic acid into vesicles and examined the insertion of a prebiotic-like peptide in a primitive bilayer. The information gathered from this research provides critical insights into molecule interactions with primitive membranes, allowing us to understand the first nanometric compartments needed to trigger further reactions that were essential for the origin of life.
Collapse
Affiliation(s)
- Romina V Sepulveda
- Center for Bioinformatics and Integrative Biology, Facultad de Ciencias de la Vida, Universidad Andres Bello, Av. República 330, Santiago 8370146, Chile
| | - Christopher Sbarbaro
- Center for Bioinformatics and Integrative Biology, Facultad de Ciencias de la Vida, Universidad Andres Bello, Av. República 330, Santiago 8370146, Chile
| | - Ma Cecilia Opazo
- Instituto de Ciencias Naturales, Facultad de Medicina Veterinaria y Agronomía, Universidad de Las Américas, Manuel Montt 948, Providencia 7500000, Chile
| | - Yorley Duarte
- Center for Bioinformatics and Integrative Biology, Facultad de Ciencias de la Vida, Universidad Andres Bello, Av. República 330, Santiago 8370146, Chile
| | - Fernando González-Nilo
- Center for Bioinformatics and Integrative Biology, Facultad de Ciencias de la Vida, Universidad Andres Bello, Av. República 330, Santiago 8370146, Chile
| | - Daniel Aguayo
- Center for Bioinformatics and Integrative Biology, Facultad de Ciencias de la Vida, Universidad Andres Bello, Av. República 330, Santiago 8370146, Chile
- Agricultura Digital, Servicio Agrícola, Salinas y Fabres S.A., Ruta 5 Sur, Parcela 165, Hijuela Larga, Paine 9540000, Chile
| |
Collapse
|
26
|
Fluorescent Probes cis- and trans-Parinaric Acids in Fluid and Gel Lipid Bilayers: A Molecular Dynamics Study. Molecules 2023; 28:molecules28052241. [PMID: 36903487 PMCID: PMC10005308 DOI: 10.3390/molecules28052241] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2023] [Revised: 02/20/2023] [Accepted: 02/23/2023] [Indexed: 03/04/2023] Open
Abstract
Fluorescence probes are indispensable tools in biochemical and biophysical membrane studies. Most of them possess extrinsic fluorophores, which often constitute a source of uncertainty and potential perturbation to the host system. In this regard, the few available intrinsically fluorescent membrane probes acquire increased importance. Among them, cis- and trans-parinaric acids (c-PnA and t-PnA, respectively) stand out as probes of membrane order and dynamics. These two compounds are long-chained fatty acids, differing solely in the configurations of two double bonds of their conjugated tetraene fluorophore. In this work, we employed all-atom and coarse-grained molecular dynamics simulations to study the behavior of c-PnA and t-PnA in lipid bilayers of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), representative of the liquid disordered and solid ordered lipid phases, respectively. All-atom simulations indicate that the two probes show similar location and orientation in the simulated systems, with the carboxylate facing the water/lipid interface and the tail spanning the membrane leaflet. The two probes establish interactions with the solvent and lipids to a similar degree in POPC. However, the almost linear t-PnA molecules have tighter lipid packing around them, especially in DPPC, where they also interact more with positively charged lipid choline groups. Probably for these reasons, while both probes show similar partition (assessed from computed free energy profiles across bilayers) to POPC, t-PnA clearly partitions more extensively than c-PnA to the gel phase. t-PnA also displays more hindered fluorophore rotation, especially in DPPC. Our results agree very well with experimental fluorescence data from the literature and allow deeper understanding of the behavior of these two reporters of membrane organization.
Collapse
|
27
|
Yadav A, Kelich P, Kallmyer NE, Reuel NF, VukoviÄ L. Characterizing the Interactions of Cell Membrane-Disrupting Peptides with Lipid-Functionalized Single-Walled Carbon Nanotube Systems for Antimicrobial Screening. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.01.25.525557. [PMID: 36747775 PMCID: PMC9900920 DOI: 10.1101/2023.01.25.525557] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/25/2023]
Abstract
Lipid-functionalized single-walled carbon nanotubes (SWNTs) have garnered significant interest for their potential use in a wide range of biomedical applications. In this work, we used molecular dynamics simulations to study the equilibrium properties of SWNTs surrounded by the phosphatidylcholine (POPC) corona phase, and their interactions with three cell membrane disruptor peptides: colistin, TAT peptide, and crotamine-derived peptide. Our results show that SWNTs favor asymmetrical positioning within the POPC corona, so that one side of the SWNT, covered by the thinnest part of the corona, comes in contact with charged and polar functional groups of POPC and water. We also observed that colistin and TAT insert deeply into POPC corona, while crotamine-derived peptide only adsorbs to the corona surface. Compared to crotamine-derived peptide, colistin and TAT also induce larger perturbations in the thinnest region of the corona, by allowing more water molecules to directly contact the SWNT surface. In separate simulations, we show that three examined peptides exhibit similar insertion and adsorption behaviors when interacting with POPC bilayers, confirming that peptide-induced perturbations to POPC in conjugates and bilayers are similar in nature and magnitude. Furthermore, we observed correlations between the peptide-induced structural perturbations and the near-infrared emission of the lipid-functionalized SWNTs, which suggest that the optical signal of the conjugates transduces the morphological changes in the lipid corona. Overall, our findings indicate that lipid-functionalized SWNTs could serve as simplified cell membrane model systems for pre-screening of new antimicrobial compounds that disrupt cell membranes.
Collapse
|
28
|
Dong Y, Kampf N, Schilt Y, Cao W, Raviv U, Klein J. Dehydration does not affect lipid-based hydration lubrication. NANOSCALE 2022; 14:18241-18252. [PMID: 36468753 PMCID: PMC9753160 DOI: 10.1039/d2nr04799c] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/01/2022] [Accepted: 11/22/2022] [Indexed: 06/17/2023]
Abstract
Phosphatidylcholine (PC) lipid bilayers at surfaces massively reduce sliding friction, via the hydration lubrication mechanism acting at their highly-hydrated phosphocholine headgroups, a central paradigm of biological lubrication, particularly at articular cartilage surfaces where low friction is crucial for joint well-being. Nanotribological measurements probed the effect on such lubrication of dehydration by dimethyl sulfoxide (DMSO), known to strongly dehydrate the phosphocholine headgroups of such PC bilayers, i.e. reduce the thickness of the inter-bilayer water layer, and thus expected to substantially degrade the hydration lubrication. Remarkably, and unexpectedly, we found that the dehydration has little effect on the friction. We used several approaches, including atomic force microscopy, small- and wide-angle X-ray scattering and all-atom molecular dynamics simulations to elucidate this. Our results show that while DMSO clearly removes hydration water from the lipid head-groups, this is offset by both higher areal head-group density and by rigidity-enhancement of the lipid bilayers, both of which act to reduce frictional dissipation. This sheds strong light on the robustness of lipid-based hydration lubrication in biological systems, despite the ubiquitous presence of bio-osmolytes which compete for hydration water.
Collapse
Affiliation(s)
- Yihui Dong
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovot, 76100, Israel.
| | - Nir Kampf
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovot, 76100, Israel.
| | - Yaelle Schilt
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
| | - Wei Cao
- Department of Physical Chemistry, School of Chemistry, Tel Aviv University, Tel Aviv, 6997801, Israel
| | - Uri Raviv
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
| | - Jacob Klein
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovot, 76100, Israel.
| |
Collapse
|
29
|
Molecular dynamics simulations to study the role of biphenylalanine in promoting the antibacterial activity of ultrashort peptides. J Mol Graph Model 2022; 117:108282. [PMID: 35961218 DOI: 10.1016/j.jmgm.2022.108282] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2022] [Revised: 07/14/2022] [Accepted: 07/24/2022] [Indexed: 01/14/2023]
Abstract
The hydrophobic amino acid biphenylalanine (B) plays a key role in the antibacterial activity of ultrashort peptides. In this study, the interactions of tetrapeptide BRBR-NH2 (BRBR) and pentapeptide BRBRB-NH2 (BRBRB) with dioleoylphosphatidylcholine/dioleoylphosphatidylglycerol (DOPC/DOPG) mixed model membrane were studied by molecular dynamics simulation to assess the role of biphenylalanine in promoting the antibacterial activity of ultrashort peptides. At low peptide concentrations, both peptides presented amphiphilic conformations; residues B of the pentapeptide approached the membrane faster than those of the tetrapeptide and made more contacts with the membrane; BRBRB exhibited stronger membrane affinity than BRBR. However, due to the low peptide concentrations, the effects of these two peptides on the membrane were not significantly different. At high peptide concentrations, the strong affinity of BRBRB made it have more interaction with membrane than BRBR and most residues B of BRBRB inserted into the membrane; BRBRB was more prone to aggregation and caused the membrane more disordered and thinner than BRBR. Hydrophobic residues often act as anchors in the antibacterial activity of ultrashort antimicrobial peptides. Adding a hydrophobic residue B to the C-terminal of BRBR could improve the ability of the peptide to "grasp" the membrane. At high peptide concentrations, the addition of residue B might enhance the antibacterial activity of the peptide. Thus, our results will be helpful in designing efficient antibacterial drugs.
Collapse
|
30
|
Vermaas JV, Crowley MF, Beckham GT. Molecular simulation of lignin-related aromatic compound permeation through gram-negative bacterial outer membranes. J Biol Chem 2022; 298:102627. [PMID: 36273587 PMCID: PMC9720347 DOI: 10.1016/j.jbc.2022.102627] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2022] [Revised: 10/14/2022] [Accepted: 10/15/2022] [Indexed: 12/12/2022] Open
Abstract
Lignin, an abundant aromatic heteropolymer in secondary plant cell walls, is the single largest source of renewable aromatics in the biosphere. Leveraging this resource for renewable bioproducts through targeted microbial action depends on lignin fragment uptake by microbial hosts and subsequent enzymatic action to obtain the desired product. Recent computational work has emphasized that bacterial inner membranes are permeable to many aromatic compounds expected from lignin depolymerization processes. In this study, we expand on these findings through simulations for 42 lignin-related compounds across a gram-negative bacterial outer membrane model. Unbiased simulation trajectories indicate that spontaneous crossing for the full outer membrane is relatively rare at molecular simulation timescales, primarily due to preferential membrane partitioning and slow diffusion within the lipopolysaccharide layer within the outer membrane. Membrane partitioning and permeability coefficients were determined through replica exchange umbrella sampling simulations to overcome sampling limitations. We find that the glycosylated lipopolysaccharides found in the outer membrane increase the permeation barrier to many lignin-related compounds, particularly the most hydrophobic compounds. However, the effect is relatively modest; at industrially relevant concentrations, uncharged lignin-related compounds will readily diffuse across the outer membrane without the need for specific porins. Together, our results provide insight into the permeability of the bacterial outer membrane for assessing lignin fragment uptake and the future production of renewable bioproducts.
Collapse
Affiliation(s)
- Josh V. Vermaas
- Biosciences Center, National Renewable Energy Laboratory, Golden, Colorado, USA,National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee, USA,MSU-DOE Plant Research Laboratory, Michigan State University, East Lansing, Michigan, USA,Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan, USA,For correspondence: Josh V. Vermaas; Michael F. Crowley; Gregg T. Beckham
| | - Michael F. Crowley
- Renewable Resources and Enabling Sciences Center, National Renewable Energy, Laboratory, Golden, Colorado, USA,For correspondence: Josh V. Vermaas; Michael F. Crowley; Gregg T. Beckham
| | - Gregg T. Beckham
- Renewable Resources and Enabling Sciences Center, National Renewable Energy, Laboratory, Golden, Colorado, USA,For correspondence: Josh V. Vermaas; Michael F. Crowley; Gregg T. Beckham
| |
Collapse
|
31
|
Fehér B, Gascoigne L, Giezen SN, Voets IK. Impact of arginine modified SNARE peptides on interactions with phospholipid bilayers and coiled-coil formation: A molecular dynamics study. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.119972] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
32
|
Nguyen TQT, Lund FW, Zanjani AAH, Khandelia H. Magic mushroom extracts in lipid membranes. BIOCHIMICA ET BIOPHYSICA ACTA. BIOMEMBRANES 2022; 1864:183957. [PMID: 35561790 DOI: 10.1016/j.bbamem.2022.183957] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/08/2021] [Revised: 04/28/2022] [Accepted: 05/02/2022] [Indexed: 06/15/2023]
Abstract
The active hallucinogen of magic mushrooms, psilocin, is being repurposed to treat nicotine addiction and treatment-resistant depression. Psilocin belongs to the tryptamine class of psychedelic compounds which include the hormone serotonin. It is believed that psilocin exerts its effect by binding to the serotonin 5-HT2A receptor. However, recent in-vivo evidence suggests that psilocin may employ a different mechanism to exert its effects. Membrane-mediated receptor desensitization of neurotransmitter receptors is one such mechanism. We compare the impact of the neutral and charged versions of psilocin and serotonin on the properties of zwitterionic and anionic lipid membranes using molecular dynamics simulations and calorimetry. Both compounds partition to the lipid interface and induce membrane thinning. The tertiary amine in psilocin, as opposed to the primary amine in serotonin, limits psilocin's impact on the membrane although more psilocin partitions into the membrane than serotonin. Calorimetry corroborates that both compounds induce a classical melting point depression like anesthetics do. Our results also lend support to a membrane-mediated receptor-binding mechanism for both psilocin and serotonin and provide physical insights into subtle chemical changes that can alter the membrane-binding of psychedelic compounds.
Collapse
Affiliation(s)
- Teresa Quynh Tram Nguyen
- Physical Life Science, Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230 Odense M, Denmark
| | - Frederik Wendelboe Lund
- Physical Life Science, Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230 Odense M, Denmark
| | - Ali Asghar Hakami Zanjani
- Physical Life Science, Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230 Odense M, Denmark.
| | - Himanshu Khandelia
- Physical Life Science, Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230 Odense M, Denmark.
| |
Collapse
|
33
|
Shamaprasad P, Frame CO, Moore TC, Yang A, Iacovella CR, Bouwstra JA, Bunge AL, McCabe C. Using molecular simulation to understand the skin barrier. Prog Lipid Res 2022; 88:101184. [PMID: 35988796 PMCID: PMC10116345 DOI: 10.1016/j.plipres.2022.101184] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2022] [Revised: 08/11/2022] [Accepted: 08/12/2022] [Indexed: 11/15/2022]
Abstract
Skin's effectiveness as a barrier to permeation of water and other chemicals rests almost entirely in the outermost layer of the epidermis, the stratum corneum (SC), which consists of layers of corneocytes surrounded by highly organized lipid lamellae. As the only continuous path through the SC, transdermal permeation necessarily involves diffusion through these lipid layers. The role of the SC as a protective barrier is supported by its exceptional lipid composition consisting of ceramides (CERs), cholesterol (CHOL), and free fatty acids (FFAs) and the complete absence of phospholipids, which are present in most biological membranes. Molecular simulation, which provides molecular level detail of lipid configurations that can be connected with barrier function, has become a popular tool for studying SC lipid systems. We review this ever-increasing body of literature with the goals of (1) enabling the experimental skin community to understand, interpret and use the information generated from the simulations, (2) providing simulation experts with a solid background in the chemistry of SC lipids including the composition, structure and organization, and barrier function, and (3) presenting a state of the art picture of the field of SC lipid simulations, highlighting the difficulties and best practices for studying these systems, to encourage the generation of robust reproducible studies in the future. This review describes molecular simulation methodology and then critically examines results derived from simulations using atomistic and then coarse-grained models.
Collapse
Affiliation(s)
- Parashara Shamaprasad
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, TN 37235-1604, United States of America; Multiscale Modeling and Simulation (MuMS) Center, Vanderbilt University, Nashville, TN 37235-1604, United States of America
| | - Chloe O Frame
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, TN 37235-1604, United States of America; Multiscale Modeling and Simulation (MuMS) Center, Vanderbilt University, Nashville, TN 37235-1604, United States of America
| | - Timothy C Moore
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, TN 37235-1604, United States of America; Multiscale Modeling and Simulation (MuMS) Center, Vanderbilt University, Nashville, TN 37235-1604, United States of America
| | - Alexander Yang
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, TN 37235-1604, United States of America; Multiscale Modeling and Simulation (MuMS) Center, Vanderbilt University, Nashville, TN 37235-1604, United States of America
| | - Christopher R Iacovella
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, TN 37235-1604, United States of America; Multiscale Modeling and Simulation (MuMS) Center, Vanderbilt University, Nashville, TN 37235-1604, United States of America
| | - Joke A Bouwstra
- Division of BioTherapeutics, LACDR, Leiden University, 2333 CC Leiden, the Netherlands
| | - Annette L Bunge
- Department of Chemical and Biological Engineering, Colorado School of Mines, Golden, CO 80401, United States of America
| | - Clare McCabe
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, TN 37235-1604, United States of America; Multiscale Modeling and Simulation (MuMS) Center, Vanderbilt University, Nashville, TN 37235-1604, United States of America; School of Engineering and Physical Science, Heriot-Watt University, Edinburgh, United Kingdom.
| |
Collapse
|
34
|
Joodaki F, Martin LM, Greenfield ML. Generation and Computational Characterization of a Complex Staphylococcus aureus Lipid Bilayer. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2022; 38:9481-9499. [PMID: 35901279 DOI: 10.1021/acs.langmuir.2c00483] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Studies indicate a crucial cell membrane role in the antibiotic resistance of Staphylococcus aureus. To simulate its membrane structure and dynamics, a complex molecular-scale computational representation of the S. aureus lipid bilayer was developed. Phospholipid types and their amounts were optimized by reverse Monte Carlo to represent characterization data from the literature, leading to 19 different phospholipid types that combine three headgroups [phosphatidylglycerol, lysyl-phosphatidylglycerol (LPG), and cardiolipin] and 10 tails, including iso- and anteiso-branched saturated chains. The averaged lipid bilayer thickness was 36.7 Å, and area per headgroup was 67.8 Å2. Phosphorus and nitrogen density profiles showed that LPG headgroups tended to be bent and oriented more parallel to the bilayer plane. The water density profile showed that small amounts reached the membrane center. Carbon density profiles indicated hydrophobic interactions for all lipids in the middle of the bilayer. Bond vector order parameters along each tail demonstrated different C-H ordering even within distinct lipids of the same type; however, all tails followed similar trends in average order parameter. These complex simulations further revealed bilayer insights beyond those attainable with monodisperse, unbranched lipids. Longer tails often extended into the opposite leaflet. Carbon at and beyond a branch showed significantly decreased ordering compared to carbon in unbranched tails; this feature arose in every branched lipid. Diverse tail lengths distributed these disordered methyl groups throughout the middle third of the bilayer. Distributions in mobility and ordering reveal diverse properties that cannot be obtained with monodisperse lipids.
Collapse
Affiliation(s)
- Faramarz Joodaki
- Department of Chemical Engineering, University of Rhode Island, 360 Fascitelli Center for Advanced Engineering, Kingston, Rhode Island 02881, United States
| | - Lenore M Martin
- Department of Cell and Molecular Biology, University of Rhode Island, 120 Flagg Road, Kingston, Rhode Island 02881, United States
| | - Michael L Greenfield
- Department of Chemical Engineering, University of Rhode Island, 360 Fascitelli Center for Advanced Engineering, Kingston, Rhode Island 02881, United States
| |
Collapse
|
35
|
Zhang Y, Liu D, Zhang Y, Guo Y, Li W, THALMANN F. Coarse-grained description of monounsaturated peroxidized phospholipidbilayers. J Chem Phys 2022; 157:085103. [DOI: 10.1063/5.0096658] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Lipid peroxides result from a reaction between cis-unsaturated lipid chains and singlet oxygen molecules leading to the addition of a peroxide OOH side group next to the acyl-chain double bond. It is now established that HP-POPC (hydroperoxidized POPC) molecules form stable, thin and laterally expanded bilayers. The difference in structural organisation arises from the hydrophilic character of the OOH side group that has strong affinity with the water interface region, leading to significant reorganisation of the bilayer. In this Article we describe a coarse-grained (CG) model of POPC and DOPC lipid peroxides within the framework of the Martini CG force-field (v2.2), derived from experimental data. We then discuss extensively the predicted structure and the influence of hydration and show how shifting the position of the unsaturated bonds along the chain changes the structure. Finally, we provide electron and neutron scattering length density profiles of the simulated bilayers.
Collapse
Affiliation(s)
| | - DongHui Liu
- Nanjing University Department of Physics, China
| | | | - Yachong Guo
- Nanjing University Kuang Yaming Honors School, China
| | | | | |
Collapse
|
36
|
Castillo SR, Rickeard BW, DiPasquale M, Nguyen MHL, Lewis-Laurent A, Doktorova M, Kav B, Miettinen MS, Nagao M, Kelley EG, Marquardt D. Probing the Link between Pancratistatin and Mitochondrial Apoptosis through Changes in the Membrane Dynamics on the Nanoscale. Mol Pharm 2022; 19:1839-1852. [PMID: 35559658 DOI: 10.1021/acs.molpharmaceut.1c00926] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Pancratistatin (PST) is a natural antiviral alkaloid that has demonstrated specificity toward cancerous cells and explicitly targets the mitochondria. PST initiates apoptosis while leaving healthy, noncancerous cells unscathed. However, the manner by which PST induces apoptosis remains elusive and impedes the advancement of PST as a natural anticancer therapeutic agent. Herein, we use neutron spin-echo (NSE) spectroscopy, molecular dynamics (MD) simulations, and supporting small angle scattering techniques to study PST's effect on membrane dynamics using biologically representative model membranes. Our data suggests that PST stiffens the inner mitochondrial membrane (IMM) by being preferentially associated with cardiolipin, which would lead to the relocation and release of cytochrome c. Second, PST has an ordering effect on the lipids and disrupts their distribution within the IMM, which would interfere with the maintenance and functionality of the active forms of proteins in the electron transport chain. These previously unreported findings implicate PST's effect on mitochondrial apoptosis.
Collapse
Affiliation(s)
- Stuart R Castillo
- Department of Chemistry and Biochemistry, University of Windsor, Windsor, Ontario N9B 3P4, Canada
| | - Brett W Rickeard
- Department of Chemistry and Biochemistry, University of Windsor, Windsor, Ontario N9B 3P4, Canada
| | - Mitchell DiPasquale
- Department of Chemistry and Biochemistry, University of Windsor, Windsor, Ontario N9B 3P4, Canada
| | - Michael H L Nguyen
- Department of Chemistry and Biochemistry, University of Windsor, Windsor, Ontario N9B 3P4, Canada
| | - Aislyn Lewis-Laurent
- Department of Chemistry and Biochemistry, University of Windsor, Windsor, Ontario N9B 3P4, Canada
| | - Milka Doktorova
- Department of Molecular Physiology and Biological Physics, University of Virginia School of Medicine, Charlottesville, Virginia 22903, United States
| | - Batuhan Kav
- Max-Planck Institute of Colloids and Interfaces, Potsdam 14476, Germany.,Institute of Biological Information Processing: Structural Biochemistry (IBI-7), Forschungszentrum Julich, Julich 52428, Germany
| | | | - Michihiro Nagao
- National Institute of Standards and Technology, Center for Neutron Research, Gaithersburg, Maryland 20899, United States.,Department of Materials Science and Engineering, University of Maryland, College Park, Maryland 20742, United States.,Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States
| | - Elizabeth G Kelley
- National Institute of Standards and Technology, Center for Neutron Research, Gaithersburg, Maryland 20899, United States
| | - Drew Marquardt
- Department of Chemistry and Biochemistry, University of Windsor, Windsor, Ontario N9B 3P4, Canada.,Department of Physics, University of Windsor, Windsor, Ontario N9B 3P4, Canada
| |
Collapse
|
37
|
Zizzi EA, Cavaglià M, Tuszynski JA, Deriu MA. Alteration of lipid bilayer mechanics by volatile anesthetics: Insights from μs-long molecular dynamics simulations. iScience 2022; 25:103946. [PMID: 35265816 PMCID: PMC8898909 DOI: 10.1016/j.isci.2022.103946] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2021] [Revised: 01/25/2022] [Accepted: 02/15/2022] [Indexed: 11/24/2022] Open
Abstract
Very few drugs in clinical practice feature the chemical diversity, narrow therapeutic window, unique route of administration, and reversible cognitive effects of volatile anesthetics. The correlation between their hydrophobicity and their potency and the increasing amount of evidence suggesting that anesthetics exert their action on transmembrane proteins, justifies the investigation of their effects on phospholipid bilayers at the molecular level, given the strong functional and structural link between transmembrane proteins and the surrounding lipid matrix. Molecular dynamics simulations of a model lipid bilayer in the presence of ethylene, desflurane, methoxyflurane, and the nonimmobilizer 1,2-dichlorohexafluorocyclobutane (also called F6 or 2N) at different concentrations highlight the structural consequences of VA partitioning in the lipid phase, with a decrease of lipid order and bilayer thickness, an increase in overall lipid lateral mobility and area-per-lipid, and a marked reduction in the mechanical stiffness of the membrane, that strongly correlates with the compounds' hydrophobicity.
Collapse
Affiliation(s)
- Eric A. Zizzi
- PolitoMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, 10129 Turin, Italy
| | - Marco Cavaglià
- PolitoMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, 10129 Turin, Italy
| | - Jack A. Tuszynski
- PolitoMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, 10129 Turin, Italy
- Department of Physics, University of Alberta, Edmonton, AB, Canada
| | - Marco A. Deriu
- PolitoMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, 10129 Turin, Italy
| |
Collapse
|
38
|
Luís JP, Mata AI, Simões CJV, Brito RMM. Analysis of the Conformational Dynamics of the Soluble and Membrane-Bound Forms of Interleukin-1 Receptor Type-1: Insights into Linker Flexibility and Domain Rearrangements. Int J Mol Sci 2022; 23:ijms23052599. [PMID: 35269739 PMCID: PMC8910350 DOI: 10.3390/ijms23052599] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2021] [Revised: 02/21/2022] [Accepted: 02/25/2022] [Indexed: 12/24/2022] Open
Abstract
Interleukin-1 receptor type 1 (IL-1R1) is a key player in inflammation and immune responses. This receptor regulates IL-1 activity in two forms: as a membrane-bound form and as a soluble ectodomain. The details and differences between the conformational dynamics of the membrane-bound and the soluble IL-1R1 ectodomains (ECDs) remain largely elusive. Here, we study and compare the structural dynamics of the soluble and membrane-bound IL-1R1-ECDs using molecular dynamics (MD) simulations, focusing on the flexible interdomain linker of the ECD, as well as the spatial rearrangements between the Ig-like domains of the ECD. To explore the membrane-bound conformations, a full-length IL-1R1 structural model was developed and subjected to classical equilibrium MD. Comparative analysis of multiple MD trajectories of the soluble and the membrane-bound IL-1R1-ECDs reveals that (i) as somewhat expected, the extent of the visited “open-to-closed” transitional states differs significantly between the soluble and membrane-bound forms; (ii) the soluble form presents open-closed transitions, sampling a wider rotational motion between the Ig-like domains of the ECD, visiting closed and “twisted” conformations in higher extent, whereas the membrane-bound form is characterized by more conformationally restricted states; (iii) interestingly, the backbone dihedral angles of residues Glu202, Glu203 and Asn204, located in the flexible linker, display the highest variations during the transition between discrete conformational states detected in IL-1R1, thus appearing to work as the “central wheel of a clock’s movement”. The simulations and analyses presented in this contribution offer a deeper insight into the structure and dynamics of IL-1R1, which may be explored in a drug discovery setting.
Collapse
Affiliation(s)
- João P. Luís
- Coimbra Chemistry Center—Institute of Molecular Sciences (CQC-IMS), Department of Chemistry, University of Coimbra, 3004-535 Coimbra, Portugal; (A.I.M.); (C.J.V.S.)
- Correspondence: (J.P.L.); (R.M.M.B.)
| | - Ana I. Mata
- Coimbra Chemistry Center—Institute of Molecular Sciences (CQC-IMS), Department of Chemistry, University of Coimbra, 3004-535 Coimbra, Portugal; (A.I.M.); (C.J.V.S.)
| | - Carlos J. V. Simões
- Coimbra Chemistry Center—Institute of Molecular Sciences (CQC-IMS), Department of Chemistry, University of Coimbra, 3004-535 Coimbra, Portugal; (A.I.M.); (C.J.V.S.)
- BSIM Therapeutics, Instituto Pedro Nunes, 3030-199 Coimbra, Portugal
| | - Rui M. M. Brito
- Coimbra Chemistry Center—Institute of Molecular Sciences (CQC-IMS), Department of Chemistry, University of Coimbra, 3004-535 Coimbra, Portugal; (A.I.M.); (C.J.V.S.)
- BSIM Therapeutics, Instituto Pedro Nunes, 3030-199 Coimbra, Portugal
- Correspondence: (J.P.L.); (R.M.M.B.)
| |
Collapse
|
39
|
Järvenpää J, Perkkiö M, Laitinen R, Lahtela-Kakkonen M. PE and PET oligomers' interplay with membrane bilayers. Sci Rep 2022; 12:2234. [PMID: 35140293 PMCID: PMC8828856 DOI: 10.1038/s41598-022-06217-4] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2021] [Accepted: 01/25/2022] [Indexed: 11/08/2022] Open
Abstract
The prevalence of microplastic pollution in nature and foodstuffs is fairly well identified. However, studies of micro- or nanoplastics' cell membrane permeation and health effects in humans are lacking. Our study focuses on examining the interactions of polyethylene (PE) and polyethylene terephthalate (PET) with bilayer membranes. We have performed molecular dynamics simulations to study how plastic oligomers behave in bilayers. In addition, we have studied membrane permeation of PE and Bis(2-hydroxyethyl) terephthalate (BHET), a type of PET monomer, with Parallel Artificial Membrane Permeability Assay (PAMPA). As a result, in simulations the molecules exhibited different movements and preferred locations in membrane. PAMPA studies suggested similar preferences in membrane, especially for PE plastic. Our results suggest that passive diffusion could be an important transport mechanism into cells for some small plastic oligomers. Both molecular dynamics simulations and PAMPA have potential for micro- and nanoplastics research.
Collapse
Affiliation(s)
- Joni Järvenpää
- School of Pharmacy, Faculty of Health Sciences, University of Eastern Finland, 70210, Kuopio, Finland.
| | - Milla Perkkiö
- School of Pharmacy, Faculty of Health Sciences, University of Eastern Finland, 70210, Kuopio, Finland
| | - Riikka Laitinen
- School of Pharmacy, Faculty of Health Sciences, University of Eastern Finland, 70210, Kuopio, Finland
| | - Maija Lahtela-Kakkonen
- School of Pharmacy, Faculty of Health Sciences, University of Eastern Finland, 70210, Kuopio, Finland.
| |
Collapse
|
40
|
Hossain SI, Islam MZ, Saha SC, Deplazes E. Drug Meets Monolayer: Understanding the Interactions of Sterol Drugs with Models of the Lung Surfactant Monolayer Using Molecular Dynamics Simulations. Methods Mol Biol 2022; 2402:103-121. [PMID: 34854039 DOI: 10.1007/978-1-0716-1843-1_9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
The lung surfactant monolayer (LSM) is a thin layer of lipids and proteins that forms the air/water interface of the alveoli. The primary function of the LSM is to reduce the surface tension at the air/water interface during breathing. The LSM also forms the main biological barrier for any inhaled particles, including drugs, to treat lung diseases. Elucidating the mechanism by which these drugs bind to and absorb into the LSM requires a molecular-level understanding of any drug-induced changes to the morphology, structure, and phase changes of the LSM.Molecular dynamics simulations have been used extensively to study the structure and dynamics of the LSM. The monolayer is usually simulated in at least two states: the compressed state, mimicking exhalation, and the expanded state, mimicking inhalation. In this chapter, we provide detailed instructions on how to set up, run, and analyze coarse-grained MD simulations to study the concentration-dependent effect of a sterol drug on the LSM, both in the expanded and compressed state.
Collapse
Affiliation(s)
- Sheikh I Hossain
- School of Life Sciences, University of Technology Sydney, Ultimo, NSW, Australia
| | - Mohammad Z Islam
- School of Mechanical and Mechatronic Engineering, University of Technology Sydney, Ultimo, NSW, Australia
- Department of Mathematics, Jashore University of Science and Technology, Jashore, Bangladesh
| | - Suvash C Saha
- School of Mechanical and Mechatronic Engineering, University of Technology Sydney, Ultimo, NSW, Australia
| | - Evelyne Deplazes
- School of Life Sciences, University of Technology Sydney, Ultimo, NSW, Australia.
| |
Collapse
|
41
|
Canner SW, Feller SE, Wassall SR. Molecular Organization of a Raft-like Domain in a Polyunsaturated Phospholipid Bilayer: A Supervised Machine Learning Analysis of Molecular Dynamics Simulations. J Phys Chem B 2021; 125:13158-13167. [PMID: 34812629 DOI: 10.1021/acs.jpcb.1c06511] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Numerous health benefits are associated with omega-3 polyunsaturated fatty acids (n-3 PUFA) consumed in fish oils. An understanding of the mechanism remains elusive. The plasma membrane as a site of action is the focus in this study. With large-scale all-atom MD simulations run on a model membrane (1050 lipid molecules), we observed the evolution over time (6 μs) of a circular (raft-like) domain composed of N-palmitoylsphingomyelin (PSM) and cholesterol embedded into a surrounding (non-raft) patch composed of polyunsaturated 1-palmitoyl-2-docosahexaenoylphosphatylcholine (PDPC) (1:1:1 mol). A supervised machine learning algorithm was developed to characterize the migration of each lipid based on molecular conformation and the local environment. PDPC molecules were seen to infiltrate the ordered raft-like domain in a small amount, while a small concentration of PSM and cholesterol molecules was seen to migrate into the disordered non-raft region. Enclosing the raft-like domain, a narrow (∼2 nm in width) interfacial zone composed of PDPC, PSM, and cholesterol that buffers the substantial difference in order (ΔSCD ≈ 0.12) between raft-like and non-raft environments was seen to form. Our results suggest that n-3 PUFA regulate the architecture of lipid rafts enriched in sphingolipids and cholesterol with a minimal effect on order within their interior in membranes.
Collapse
Affiliation(s)
- Samuel W Canner
- Department of Physics, IUPUI, Indianapolis, Indiana 46202-3273, United States.,Department of Computer and Information Science, IUPUI, Indianapolis, Indiana 46202-5132, United States
| | - Scott E Feller
- Department of Chemistry, Wabash College, Crawfordsville, Indiana 47933, United States
| | - Stephen R Wassall
- Department of Physics, IUPUI, Indianapolis, Indiana 46202-3273, United States
| |
Collapse
|
42
|
Singhal A, Agur Sevink GJ. The role of size and nature in nanoparticle binding to a model lung membrane: an atomistic study. NANOSCALE ADVANCES 2021; 3:6635-6648. [PMID: 36132649 PMCID: PMC9417560 DOI: 10.1039/d1na00578b] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/24/2021] [Accepted: 09/16/2021] [Indexed: 05/09/2023]
Abstract
Understanding the uptake of nanoparticles (NPs) by different types of cellular membranes plays a pivotal role in the design of NPs for medical applications and in avoiding adverse effects that result in nanotoxicity. Yet, the role of key design parameters, such as the bare NP material, NP size and surface reactivity, and the nature of NP coatings, in membrane remodelling and uptake mechanisms is still very poorly understood, particularly towards the lower range of NP dimensions that are beyond the experimental imaging resolution. The same can be said about the role of a particular membrane composition. Here, we systematically employ biased and unbiased molecular dynamics simulations to calculate the binding energy for three bare materials (Ag/SiO2/TiO2) and three NP sizes (1/3/5 nm diameter) with a representative lung surfactant membrane, and to study their binding kinetics. The calculated binding energies show that irrespective of size, Ag nanoparticles bind very strongly to the bilayer, while the NPs made of SiO2 or TiO2 experience very low to no binding. The unbiased simulations provide insight into how the NPs and membrane affect each other in terms of the solvent-accessible surface area (SASA) of the NPs and the defect types and fluidity of the membrane. Using these systematic fine-grained results in coarsening procedures will pave the way for simulations considering NP sizes that are well beyond the membrane thickness, i.e. closer to experimental dimensions, for which different binding characteristics and more significant membrane remodelling are expected.
Collapse
Affiliation(s)
- Ankush Singhal
- Leiden Institute of Chemistry, Leiden University P.O. Box 9502 2300 RA Leiden The Netherlands
| | - G J Agur Sevink
- Leiden Institute of Chemistry, Leiden University P.O. Box 9502 2300 RA Leiden The Netherlands
| |
Collapse
|
43
|
Zhang L. Interaction of Human β Defensin Type 3 (hBD-3) with Different PIP2-Containing Membranes, a Molecular Dynamics Simulation Study. J Chem Inf Model 2021; 61:4670-4686. [PMID: 34473496 DOI: 10.1021/acs.jcim.1c00805] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
Human β defensin type 3 (hBD-3) is a cysteine-rich small antibacterial peptide. It belongs to the human innate immune system. hBD-3 has six cysteine residues, which form three pairs of disulfide bonds, and those bonds break in the reducing condition. It is known that hBD-3 can interact with bacterial membrane, and even eukaryotic cell membrane, which has a low concentration of phosphatidylinositol 4,5-bisphosphate (PIP2) lipids. PIP2 is a vital component in cell membranes and has been found to play important roles during antimicrobial peptide (AMP) interaction with membranes. To understand the functional mechanism of hBD-3 interacting with PIP2-containing membranes, the binding structures of hBD-3 on 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers mixed with 10% of PIP2 were predicted using two kinds of methods. The first one is by placing the hBD-3 monomer in different orientations above the POPC + 10%PIP2 membrane to set up five different initial simulation systems and performing long-term simulations on each to predict the most stable binding structure. It was found that hBD-3 analogue binds on the mixed lipid membrane on the two loop regions. The second method is by running long-term simulations on one or nine hBD-3 dimers binding on POPC mixed with 10%PIP2 lipid bilayer starting from the solid-state NMR (ssNMR)-suggested orientation. The dimer dissociated, and the most stable binding of hBD-3 in wild-type on the mixed membrane is also through the two loop regions, which agrees with the prediction from both the first method and the lipid self-assembly result. The PIP2 lipids can form long-lasting hydrogen bonds with positively charged residues such as Arg and Lys on hBD-3, thus forming clusters with hBD-3. As a comparison, hBD-3 dimers binding with a combined bilayer having 1,2-palmitoyl-oleoyl-sn-glycero-3-phosphoserine (POPS) on the upper and POPC on the lower leaflets and the combined POPS + POPC bilayer mixing with 10%PIP2 were also studied. The long-term simulation result shows that hBD-3 can bind with the heads of negatively charged POPS and PIP2 lipids and form hydrogen bonds. The stable binding sites of hBD-3 on PIP2 or POPS mixed bilayers are still on the two loop regions. On the combined POPS + POPC mixed with 10%PIP2 bilayer, the binding of hBD-3 with PIP2 lipids became less stable and fewer because of the competition of binding with the POPS lipids. Besides that, binding with hBD-3 can decrease the membrane thickness of the POPC + PIP2, POPS + POPC, and POPS + POPC + PIP2 bilayers and make POPS and PIP2 lipids more flexible based on the order parameter calculations. Our results supply molecular insight on AMP binding with different membranes and can help understand the functional mechanism of hBD-3 disrupting PIP2-containing membranes.
Collapse
Affiliation(s)
- Liqun Zhang
- Department of Chemical Engineering, Tennessee Technological University, Cookeville, Tennessee 38505, United States
| |
Collapse
|
44
|
Bouquiaux C, Castet F, Champagne B. Unravelling the Effects of Cholesterol on the Second-Order Nonlinear Optical Responses of Di-8-ANEPPS Dye Embedded in Phosphatidylcholine Lipid Bilayers. J Phys Chem B 2021; 125:10195-10212. [PMID: 34491062 DOI: 10.1021/acs.jpcb.1c05630] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
Cholesterol is known for its role in maintaining the correct fluidity and rigidity of the animals cell membranes and thus their functions. Assessing the content and the role of cholesterol in lipid bilayers is therefore of crucial importance for a deeper understanding and control of membrane functioning. In this computational work, we investigate bilayers built from three types of glycerophospholipid phosphatidylcholine (PC) lipids, namely dipalmitoylphosphatidylcholine (DPPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), and dioleoylphosphatidylcholine (DOPC), and containing different amounts of cholesterol by analyzing the second-harmonic generation (SHG) nonlinear optical (NLO) response of a probe molecule, di-8-ANEPPS, inserted into the membranes. This molecular property presents the advantage to be specific to interfacial regions such as lipid bilayers. To unravel these effects, Molecular Dynamics (MD) simulations have been performed on both DPPC and DOPC lipids by varying the cholesterol mole fraction (from 0 to 0.66), while POPC was only considered as a pure bilayer. In the case of the structural properties of the bilayers, all the analyses converge toward the same conclusion: as the mole fraction of cholesterol increases, the systems become more rigid, confirming the condensing effect of cholesterol. In addition, the chromophore is progressively more aligned with respect to the normal to the bilayer. On the contrary, addition of unsaturation disorders the lipid bilayers, with barely no impact on the alignment of the chromophore. Then, using the frames obtained from the MD simulations, the first hyperpolarizability β of the dye in its environment has been computed at the TDDFT level. On the one hand, the addition of cholesterol induces a progressive increase of the diagonal component the β tensor parallel to the bilayer normal. On the other hand, larger β values have been calculated for the unsaturated than for the saturated lipid systems. In summary, this study illustrates the relationship between the composition and structure of the bilayers and the NLO responses of the embedded dye.
Collapse
Affiliation(s)
- Charlotte Bouquiaux
- Theoretical Chemistry Lab, Unit of Theoretical and Structural Physical Chemistry, Namur Institute of Structured Matter, University of Namur, rue de Bruxelles, 61, B-5000 Namur, Belgium
| | - Frédéric Castet
- , Institut des Sciences Moléculaires, UMR 5255 CNRS, University of Bordeaux, cours de la Libération 351, F-33405 Talence Cedex, France
| | - Benoît Champagne
- Theoretical Chemistry Lab, Unit of Theoretical and Structural Physical Chemistry, Namur Institute of Structured Matter, University of Namur, rue de Bruxelles, 61, B-5000 Namur, Belgium
| |
Collapse
|
45
|
Jung I, Choi J, Nam J, No KT. Modeling lipid layers of atopic skin and observation of changes in lipid layer properties with changes in ceramide content. J Cosmet Dermatol 2021; 20:2924-2931. [PMID: 33238053 PMCID: PMC8451888 DOI: 10.1111/jocd.13861] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2020] [Revised: 07/09/2020] [Accepted: 11/18/2020] [Indexed: 01/11/2023]
Abstract
BACKGROUND Studies have shown that there is a high correlation between atopic dermatitis and decrease in ceramide content in the lipid bilayer of skin. Moreover, it has been shown that the reduction in ceramide content in the stratum corneum is unique to atopic dermatitis, indicating that there are particular structural differences between the lipid bilayers of normal and atopic skin. AIM This study aimed to compare the lipid bilayer of the atopic skin with that of the healthy skin and to establish a structural model of the lipid bilayer for atopy. METHODS Molecular dynamics simulations were performed using NAMD 2.8. Models of lipid bilayers of normal skin and atopic skin, and a model of lipid bilayer containing only ceramide were built with CHARMM-GUI. The thickness, area occupied per lipid, and alignment of lipids were compared among the three models. Potential mean force (PMF) of the sodium laureth sulfate (SLES) on lipid bilayers was calculated to predict the affinity between SLES and lipid bilayers. RESULTS Potential mean force calculations showed that the lipid bilayer of atopic skin was able to absorb the surfactant more easily than that of normal skin. CONCLUSIONS When the ceramide ratio is low, the thickness of lipid bilayer is reduced and its structure is weakened. Other structural differences between the lipid layers of normal and atopic skin included increased area per lipid and poor alignment of lipids. Further, the atopy lipid bilayer model was found to absorb more SLES than the normal skin lipid bilayer model.
Collapse
Affiliation(s)
- In‐Keun Jung
- Amore‐Pacific Research and Development CenterYonginKorea
- Yonsei UniversitySeoulKorea
| | - Joonho Choi
- Amore‐Pacific Research and Development CenterYonginKorea
| | - Jin Nam
- Amore‐Pacific Research and Development CenterYonginKorea
| | | |
Collapse
|
46
|
Kamel M, Löwe M, Schott-Verdugo S, Gohlke H, Kedrov A. Unsaturated fatty acids augment protein transport via the SecA:SecYEG translocon. FEBS J 2021; 289:140-162. [PMID: 34312977 DOI: 10.1111/febs.16140] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2021] [Revised: 07/01/2021] [Accepted: 07/23/2021] [Indexed: 12/22/2022]
Abstract
The translocon SecYEG and the associated ATPase SecA form the primary protein secretion system in the cytoplasmic membrane of bacteria. The secretion is essentially dependent on the surrounding lipids, but the mechanistic understanding of their role in SecA : SecYEG activity is sparse. Here, we reveal that the unsaturated fatty acids (UFAs) of the membrane phospholipids, including tetraoleoyl-cardiolipin, stimulate SecA : SecYEG-mediated protein translocation up to ten-fold. Biophysical analysis and molecular dynamics simulations show that UFAs increase the area per lipid and cause loose packing of lipid head groups, where the N-terminal amphipathic helix of SecA docks. While UFAs do not affect the translocon folding, they promote SecA binding to the membrane, and the effect is enhanced up to fivefold at elevated ionic strength. Tight SecA : lipid interactions convert into the augmented translocation. Our results identify the fatty acid structure as a notable factor in SecA : SecYEG activity, which may be crucial for protein secretion in bacteria, which actively change their membrane composition in response to their habitat.
Collapse
Affiliation(s)
- Michael Kamel
- Synthetic Membrane Systems, Institute for Biochemistry, Heinrich Heine University Düsseldorf, Germany
| | - Maryna Löwe
- Synthetic Membrane Systems, Institute for Biochemistry, Heinrich Heine University Düsseldorf, Germany
| | - Stephan Schott-Verdugo
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, Germany.,John von Neumann Institute for Computing (NIC), Jülich Supercomputing Centre (JSC), Institute of Biological Information Processing (IBI-7: Structural Bioinformatics), and Institute of Bio- and Geosciences (IBG-4: Bioinformatics), Forschungszentrum Jülich GmbH, Germany
| | - Holger Gohlke
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, Germany.,John von Neumann Institute for Computing (NIC), Jülich Supercomputing Centre (JSC), Institute of Biological Information Processing (IBI-7: Structural Bioinformatics), and Institute of Bio- and Geosciences (IBG-4: Bioinformatics), Forschungszentrum Jülich GmbH, Germany
| | - Alexej Kedrov
- Synthetic Membrane Systems, Institute for Biochemistry, Heinrich Heine University Düsseldorf, Germany
| |
Collapse
|
47
|
Borocci S, Bozzuto G, Bombelli C, Ceccacci F, Formisano G, Stringaro A, Molinari A, Mancini G. How stereochemistry of lipid components can affect lipid organization and the route of liposome internalization into cells. NANOSCALE 2021; 13:11976-11993. [PMID: 34212969 DOI: 10.1039/d1nr02175c] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Though liposome-based drugs are in clinical use, the mechanism of cell internalization of liposomes is yet an object of controversy. The present experimental investigation, carried out on human glioblastoma cells, indicated different internalization routes for two diastereomeric liposomes. Molecular dynamics simulations of the lipid bilayers of the two formulations indicated that the different stereochemistry of a lipid component controls some parameters such as area per lipid molecule and fluidity of lipid membranes, surface potential and water organization at the lipid/water interface, all of which affect the interaction with biomolecules and cell components.
Collapse
Affiliation(s)
- Stefano Borocci
- Dipartimento per la Innovazione nei sistemi Biologici, Agroalimentari e Forestali (DIBAF), Università degli Studi della Tuscia, L.go dell'Università, s.n.c., 01100 Viterbo, Italy.
| | | | | | | | | | | | | | | |
Collapse
|
48
|
Kjølbye LR, De Maria L, Wassenaar TA, Abdizadeh H, Marrink SJ, Ferkinghoff-Borg J, Schiøtt B. General Protocol for Constructing Molecular Models of Nanodiscs. J Chem Inf Model 2021; 61:2869-2883. [PMID: 34048229 DOI: 10.1021/acs.jcim.1c00157] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Nanodisc technology is increasingly being applied for structural and biophysical studies of membrane proteins. In this work, we present a general protocol for constructing molecular models of nanodiscs for molecular dynamics simulations. The protocol is written in python and based on geometric equations, making it fast and easy to modify, enabling automation and customization of nanodiscs in silico. The novelty being the ability to construct any membrane scaffold protein (MSP) variant fast and easy given only an input sequence. We validated and tested the protocol by simulating seven different nanodiscs of various sizes and with different membrane scaffold proteins, both circularized and noncircularized. The structural and biophysical properties were analyzed and shown to be in good agreement with previously reported experimental data and simulation studies.
Collapse
Affiliation(s)
- Lisbeth R Kjølbye
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark.,Novo Nordisk A/S, Novo Nordisk Park 1, 2760 Måløv, Denmark
| | | | - Tsjerk A Wassenaar
- Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Haleh Abdizadeh
- Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | | | - Birgit Schiøtt
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark.,Interdisciplinary Nanoscience Center, Aarhus University, Gustav Wieds vej 14, 8000 Aarhus C, Denmark
| |
Collapse
|
49
|
Wang H, Falcoz S, Berteau JP. Long-Chain Fatty Acids in Bones and Their Link to Submicroscopic Vascularization Network: NMR Assignment and Relaxation Studies under Magic Angle Spinning Conditions in Intramuscular Bones of Atlantic Herring Fish. J Phys Chem B 2021; 125:4585-4595. [PMID: 33914538 DOI: 10.1021/acs.jpcb.1c00186] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The long-lasting proton signals in bones are identified as long-chain fatty acids, including saturated, mono-, and di-unsaturated fatty acids, with direct nuclear magnetic resonance evidence. We used intramuscular bones from Atlantic Herring fish to avoid interference from lipid-rich marrows. The key is to recognize that these signals are from mobile phase materials and study them with J-coupled correlation spectroscopies under magic angle spinning conditions. We kept extensive 1H-spin-echo records that allowed us to examine the effect of magic angle spinning on the transverse relaxation time of water and lipids over time. While it is impossible to distinguish based on chemical shifts, the relaxation data suggest that the signals are more consistent with the interpretation of phospholipid membranes than triglycerides in lipid droplets. In particular, the simultaneous T2 changes in water and lipids suggest that the centrifugal impact of magic angle spinning alters the lipid's structure in very tight spaces. Additional evidence of phospholipid membranes came from the choline-γ resonance at 3.2 ppm in fresh samples, which disappears with magic angle spinning. Thus, the fatty acid signals are at least partially from membrane bilayer structures, and we propose that they are linked to the submicroscopic vascularization channels similar to the dense canaliculi network in mammalian bones. Our detection of phospholipids from bones depended critically on two factors: (1) the elimination of the overwhelming triglyceride signals from marrows and (2) the preservation of water that biomembranes require. The relaxation data reveal aspects of lipid fluidity that have not been elucidated by previous order parameter studies on model membranes. Relaxation times have long been considered difficult to interpret. A robust and renewed understanding may be beneficial.
Collapse
Affiliation(s)
- Hsin Wang
- Department of Chemistry and Biochemistry, The City College of New York and CUNY Institute for Macromolecular Assemblies, 85 St. Nicholas Terrace, New York, New York 10031, United States
| | - Steve Falcoz
- Department of Physical Therapy, The College of Staten Island, 2800 Victory Blvd, Staten Island, New York 10314, United States
| | - Jean-Philippe Berteau
- Department of Physical Therapy, The College of Staten Island, 2800 Victory Blvd, Staten Island, New York 10314, United States.,New York Centre for Biomedical Engineering, City University of New York - City College of New York, New York, New York 10031, United States.,Nanosciences Initiative, City University of New York - Advance Science Research Center, New York, New York 10031, United States
| |
Collapse
|
50
|
Martins IBS, Viegas TG, Dos Santos Alvares D, de Souza BM, Palma MS, Ruggiero Neto J, de Araujo AS. The effect of acidic pH on the adsorption and lytic activity of the peptides Polybia-MP1 and its histidine-containing analog in anionic lipid membrane: a biophysical study by molecular dynamics and spectroscopy. Amino Acids 2021; 53:753-767. [PMID: 33890127 DOI: 10.1007/s00726-021-02982-0] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2021] [Accepted: 04/12/2021] [Indexed: 01/09/2023]
Abstract
Antimicrobial peptides (AMPs) are part of the innate immune system of many species. AMPs are short sequences rich in charged and non-polar residues. They act on the lipid phase of the plasma membrane without requiring membrane receptors. Polybia-MP1 (MP1), extracted from a native wasp, is a broad-spectrum bactericide, an inhibitor of cancer cell proliferation being non-hemolytic and non-cytotoxic. MP1 mechanism of action and its adsorption mode is not yet completely known. Its adsorption to lipid bilayer and lytic activity is most likely dependent on the ionization state of its two acidic and three basic residues and consequently on the bulk pH. Here we investigated the effect of bulk acidic (pH 5.5) and neutral pH (7.4) solution on the adsorption, insertion, and lytic activity of MP1 and its analog H-MP1 to anionic (7POPC:3POPG) model membrane. H-MP1 is a synthetic analog of MP1 with lysines replaced by histidines. Bulk pH changes could modulate this peptide efficiency. The combination of different experimental techniques and molecular dynamics (MD) simulations showed that the adsorption, insertion, and lytic activity of H-MP1 are highly sensitive to bulk pH in opposition to MP1. The atomistic details, provided by MD simulations, showed peptides contact their N-termini to the bilayer before the insertion and then lay parallel to the bilayer. Their hydrophobic faces inserted into the acyl chain phase disturb the lipid-packing.
Collapse
Affiliation(s)
- Ingrid Bernardes Santana Martins
- Department of Physics, IBILCE, UNESP-São Paulo State University, Cristóvão Colombo, 2265-Jardim Nazareth, São José do Rio Preto, SP, 15054-000, Brazil
| | - Taisa Giordano Viegas
- Department of Physics, IBILCE, UNESP-São Paulo State University, Cristóvão Colombo, 2265-Jardim Nazareth, São José do Rio Preto, SP, 15054-000, Brazil
| | - Dayane Dos Santos Alvares
- Department of Physics, IBILCE, UNESP-São Paulo State University, Cristóvão Colombo, 2265-Jardim Nazareth, São José do Rio Preto, SP, 15054-000, Brazil
| | - Bibiana Monson de Souza
- Department of Basic and Applied Biology, Institute of Biosciences, UNESP-São Paulo State University, Rio Claro, SP, Brazil
| | - Mário Sérgio Palma
- Department of Basic and Applied Biology, Institute of Biosciences, UNESP-São Paulo State University, Rio Claro, SP, Brazil
| | - João Ruggiero Neto
- Department of Physics, IBILCE, UNESP-São Paulo State University, Cristóvão Colombo, 2265-Jardim Nazareth, São José do Rio Preto, SP, 15054-000, Brazil.
| | - Alexandre Suman de Araujo
- Department of Physics, IBILCE, UNESP-São Paulo State University, Cristóvão Colombo, 2265-Jardim Nazareth, São José do Rio Preto, SP, 15054-000, Brazil.
| |
Collapse
|