1
|
Sahu S, Garg A, Saini R, Debnath A. Interface Water Assists in Dimethyl Sulfoxide Crossing and Poration in Model Bilayer. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2024; 40:5764-5775. [PMID: 38445595 DOI: 10.1021/acs.langmuir.3c03498] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/07/2024]
Abstract
Understanding the mechanism of transport and pore formation by a commonly used cryoprotectant, dimethyl sulfoxide (DMSO), across cell membranes is fundamentally crucial for drug delivery and cryopreservation. To shed light on the mechanism and thermodynamics of pore formation and crossing behavior of DMSO, extensive all-atom molecular dynamics simulations of 1,2-dimyristoyl-rac-glycero-3-phosphocholine (DMPC) bilayers are performed at various concentrations of DMSO at a temperature above the physiological temperature. Our results unveil that DMSO partially depletes water from the interface and positions itself between lipid heads without full dehydration. This induces a larger area per headgroup, increased disorder, and enhanced fluidity without any disintegration even at the highest DMSO concentration studied. The enhanced disorder fosters local fluctuations at the interface that nucleate dynamic and transient pores. The potential of mean force (PMF) of DMSO crossing is derived from two types of biased simulations: a single DMSO pulling using the umbrella sampling technique and a cylindrical pore formation using the recently developed chain reaction coordinate method. In both cases, DMSO crossing encounters a barrier attributed to unfavorable polar nonpolar interactions between DMSO and lipid tails. As the DMSO concentration increases, the barrier height reduces along with the faster lateral and perpendicular diffusion of DMSO suggesting favorable permeation. Our findings suggest that the energy required for pore formation decreases when water assists in the formation of DMSO pores. Although DMSO displaces water from the interface toward the far interface region without complete dehydration, the presence of interface water diminishes pore formation free energy. The existence of interface water leads to the formation of a two-dimensional percolated water-DMSO structure at the interface, which is absent otherwise. Overall, these insights into the mechanism of DMSO crossing and pore formation in the bilayer will contribute to understanding cryoprotectant behavior under supercooled conditions in the future.
Collapse
Affiliation(s)
- Samapika Sahu
- Department of Chemistry, Indian Institute of Technology Jodhpur, Jodhpur 342037, India
| | - Avinash Garg
- Department of Chemistry, Indian Institute of Technology Jodhpur, Jodhpur 342037, India
| | - Renu Saini
- Department of Chemistry, Indian Institute of Technology Jodhpur, Jodhpur 342037, India
| | - Ananya Debnath
- Department of Chemistry, Indian Institute of Technology Jodhpur, Jodhpur 342037, India
| |
Collapse
|
2
|
Sahu S, Srinivasan H, Jadhav SE, Sharma VK, Debnath A. Aspirin-Induced Ordering and Faster Dynamics of a Cationic Bilayer for Drug Encapsulation. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2023; 39:16432-16443. [PMID: 37948158 DOI: 10.1021/acs.langmuir.3c02241] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2023]
Abstract
The lipid dynamics and phase play decisive roles in drug encapsulation and delivery to the intracellular target. Thus, understanding the dynamic and structural alterations of membranes induced by drugs is essential for targeted delivery. To this end, united-atom molecular dynamics simulations of a model bilayer, dioctadecyldimethylammonium bromide (DODAB), are performed in the absence and presence of the usual nonsteroidal anti-inflammatory drug (NSAID), aspirin, at 298, 310, and 345 K. At 298 and 310 K, the bilayers are in the interdigitated two-dimensional square phases, which become rugged in the presence of aspirin, as evident from height fluctuations. At 345 K, the bilayer is in the fluid phase in both the absence and presence of aspirin. Aspirin is preferentially located near the oppositely charged headgroup and creates void space, which leads to an increase in the interdigitation and order parameters. Although the center of mass of lipids experiences structural arrest, they reach the diffusive regime faster and have higher lateral diffusion constants in the presence of aspirin. Results are found to be consistent with recent quasi-elastic neutron scattering studies that reveal that aspirin acts as a plasticizer and enhances lateral diffusion of lipids in both ordered and fluid phases. Different relaxation time scales of the bonds along the alkyl tails of DODAB due to the multitude of lipid motions become faster upon the addition of aspirin. Our results show that aspirin insertion is most favorable at physiological temperature. Thus, the ordered, more stable, and faster DODAB bilayer can be a potential drug carrier for the protected encapsulation of aspirin, followed by targeted and controlled drug release with antibacterial activity in the future.
Collapse
Affiliation(s)
- Samapika Sahu
- Department of Chemistry, Indian Institute of Technology, Jodhpur 342037, India
| | - Harish Srinivasan
- Solid State Physics Division, Bhabha Atomic Research Center, Mumbai 400085, India
- Homi Bhabha National Institute, Mumbai 400094, India
| | - Sankalp E Jadhav
- Department of Chemistry, Indian Institute of Technology, Jodhpur 342037, India
| | - Veerendra K Sharma
- Solid State Physics Division, Bhabha Atomic Research Center, Mumbai 400085, India
- Homi Bhabha National Institute, Mumbai 400094, India
| | - Ananya Debnath
- Department of Chemistry, Indian Institute of Technology, Jodhpur 342037, India
| |
Collapse
|
3
|
Saini R, Debnath A. Thylakoid Composition Facilitates Chlorophyll a Dimerization through Stronger Interlipid Interactions. J Phys Chem B 2023; 127:9082-9094. [PMID: 37819861 DOI: 10.1021/acs.jpcb.3c04942] [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: 10/13/2023]
Abstract
Plant thylakoid membrane serves as a crucial matrix for the aggregation of chlororophyll a (CLA) pigments, essential for light harvesting. To understand the role of lipid compositions in the stability of CLA aggregates, dimerization of chlorophyll a molecules (CLA) is studied in the presence of the thylakoid and the bilayers comprising either the least or the highest unsaturated lipids by using coarse-grained molecular dynamics simulations. The thylakoid membrane enhances the stability of the CLA dimer compared with other membranes due to very strong lipid-lipid interactions. The thylakoid exhibits a distinct distribution of lipids around the CLA dimer. Less unsaturated lipids reside in close proximity to the dimer, promoting increased order and efficient packing. Conversely, higher unsaturated lipids are depleted from the dimer, imparting flexibility to the membrane. The combination of tight packing near the dimer and membrane flexibility away from the dimer enhances the stability of the dimer in the thylakoid membrane. Our results suggest that lipid mixing, rather than lipid unsaturation, plays a critical role in facilitating CLA dimerization by modulating the membrane microenvironment through stronger lipid-lipid interactions. These insights will be useful in understanding how lipid compositions affect efficient light absorption and energy transfer during photosynthesis in the future.
Collapse
Affiliation(s)
- Renu Saini
- Department of Chemistry, Indian Institute of Technology, Jodhpur, Rajasthan 342030, India
| | - Ananya Debnath
- Department of Chemistry, Indian Institute of Technology, Jodhpur, Rajasthan 342030, India
| |
Collapse
|
4
|
Malik S, Karmakar S, Debnath A. Quantifying dynamical heterogeneity length scales of interface water across model membrane phase transitions. J Chem Phys 2023; 158:091103. [PMID: 36889951 DOI: 10.1063/5.0137727] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2023] Open
Abstract
All-atom molecular dynamics simulations of 1,2-dimyristoyl-sn-glycero-3-phosphocholine lipid membranes reveal a drastic growth in the heterogeneity length scales of interface water (IW) across fluid to ripple to gel phase transitions. It acts as an alternate probe to capture the ripple size of the membrane and follows an activated dynamical scaling with the relaxation time scale solely within the gel phase. The results quantify the mostly unknown correlations between the spatiotemporal scales of the IW and membranes at various phases under physiological and supercooled conditions.
Collapse
Affiliation(s)
- Sheeba Malik
- Department of Chemistry, IIT Jodhpur, Jodhpur, Rajasthan, India
| | - Smarajit Karmakar
- Centre for Interdisciplinary Sciences, Tata Institute of Fundamental Research, Hyderabad, India
| | - Ananya Debnath
- Department of Chemistry, IIT Jodhpur, Jodhpur, Rajasthan, India
| |
Collapse
|
5
|
Phan A, Stamatakis M, Koh CA, Striolo A. Wetting Properties of Clathrate Hydrates in the Presence of Polycyclic Aromatic Compounds: Evidence of Ion-Specific Effects. J Phys Chem Lett 2022; 13:8200-8206. [PMID: 36006399 PMCID: PMC9442800 DOI: 10.1021/acs.jpclett.2c01846] [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] [Indexed: 05/14/2023]
Abstract
Polycyclic aromatic hydrocarbons (PAHs) have attracted remarkable multidisciplinary attention due to their intriguing π-π stacking configurations, showing enormous opportunity for their use in a variety of advanced applications. To secure progress, detailed knowledge on PAHs' interfacial properties is required. Employing molecular dynamics, we probe the wetting properties of brine droplets (KCl, NaCl, and CaCl2) on sII methane-ethane hydrate surfaces immersed in various oil solvents. Our simulations show synergistic effects due to the presence of PAHs compounded by ion-specific effects. Our analysis reveals phenomenological correlations between the wetting properties and a combination of the binding free-energy difference and entropy changes upon oil solvation for PAHs at oil/brine and oil/hydrate interfaces. The detailed thermodynamic analysis conducted upon the interactions between PAHs and various interfaces identifies molecular-level mechanisms responsible for wettability alterations, which could be applicable for advancing applications in optics, microfluidics, biotechnology, medicine, as well as hydrate management.
Collapse
Affiliation(s)
- Anh Phan
- Department
of Chemical and Process Engineering, Faculty of Engineering and Physical
Sciences, University of Surrey, Guildford, Surrey GU2 7XH, United
Kingdom
| | - Michail Stamatakis
- Department
of Chemical Engineering, University College
London, London WC1E 7JE, United Kingdom
| | - Carolyn A. Koh
- Center
for Hydrate Research, Chemical & Biological Engineering Department, Colorado School of Mines, Golden, Colorado 80401, United States
| | - Alberto Striolo
- Department
of Chemical Engineering, University College
London, London WC1E 7JE, United Kingdom
- School
of Chemical, Biological and Materials Engineering, University of Oklahoma, Norman, Oklahoma 73019, United States
| |
Collapse
|
6
|
Gupta S, Mandal T. Simulation study of domain formation in a model bacterial membrane. Phys Chem Chem Phys 2022; 24:18133-18143. [PMID: 35856570 DOI: 10.1039/d2cp01873j] [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/2022]
Abstract
Recent experimental studies revealed that functional membrane microdomains (FMMs) are formed in prokaryotic cells which are structurally and functionally similar to the lipid rafts formed in eukaryotic cells. In this study, we employ coarse-grained molecular dynamics simulations to investigate the mechanism of domain formation and its physiochemical properties in a model methicillin-resistant staphylococcus aureus (MRSA) cell membrane. We find that domains are formed through lateral segregation of staphyloxanthin (STX), a carotenoid which shields the bacteria from the host's immune because of its antioxidant nature. Simulation results suggest that membrane integrity increases with the size of the domain, which is assessed by computing bond order parameter of the lipid tails, membrane expansion modulus and water permeability across the membrane. Various membrane domain proteins such as flotillin-like protein floA and penicillin binding protein (PBP2a) preferentially bind with the STX and accumulate in the membrane domain which is consistent with the recent experimental results.
Collapse
Affiliation(s)
- Shivam Gupta
- Department of Physics, Indian Institute of Technology Kanpur, Kanpur-208016, India.
| | - Taraknath Mandal
- Department of Physics, Indian Institute of Technology Kanpur, Kanpur-208016, India.
| |
Collapse
|
7
|
Bhendale M, Srivastava A, Singh JK. Insights into the Phase Diagram of Pluronic L64 Using Coarse-Grained Molecular Dynamics Simulations. J Phys Chem B 2022; 126:4731-4744. [PMID: 35708274 DOI: 10.1021/acs.jpcb.2c02429] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We investigate the concentration-dependent phase diagram of pluronic L64 in aqueous media at 300 and 320 K using coarse-grained (CG) molecular dynamics (MD) simulations. The CG model is derived by adapting the Martini model for nonbonded interactions along with the Boltzmann inversion (BI) of bonded interactions from all-atom (AA) simulations. Our derived CG model successfully captures the experimentally observed micellar-, hexagonal-, lamellar-, and polymer-rich solution phase. The end-to-end distance reveals the conformational change from an open-chain structure in the micellar phase to a folded-chain structure in the lamellar phase, increasing the orientational order. An increase in temperature leads to expulsion of water molecules from the L64 moiety, suggesting an increase in L64 hydrophobicity. Thermodynamic analysis using the two-phase thermodynamics (2PT) method suggests the entropy of the system decreases with increasing L64 concentration and the decrease in free energy (F) with temperature is mainly driven by the entropic factor (-TS). Further, the increase in aggregation number at lower concentrations and self-assembly at very high concentrations is energetically driven, whereas the change from the micellar phase to the lamellar phase with increasing L64 concentration is entropically driven. Our model provides molecular insights into L64 phases which can be further explored to design functionality-based suprastructures for drug delivery and tissue engineering applications.
Collapse
Affiliation(s)
- Mangesh Bhendale
- Department of Chemical Engineering, Indian Institute of Technology Kanpur, Kanpur, Uttar Pradesh 208016, India
| | - Arpita Srivastava
- Department of Chemical Engineering, Indian Institute of Technology Kanpur, Kanpur, Uttar Pradesh 208016, India
| | - Jayant K Singh
- Department of Chemical Engineering, Indian Institute of Technology Kanpur, Kanpur, Uttar Pradesh 208016, India.,Prescience Insilico Private Limited, Fifth Floor, Novel MSR Building, Marathahalli, Bengaluru, Karnataka 560037, India
| |
Collapse
|
8
|
Gupta C, Sarkar D, Tieleman DP, Singharoy A. The ugly, bad, and good stories of large-scale biomolecular simulations. Curr Opin Struct Biol 2022; 73:102338. [PMID: 35245737 DOI: 10.1016/j.sbi.2022.102338] [Citation(s) in RCA: 27] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2021] [Revised: 11/29/2021] [Accepted: 01/24/2022] [Indexed: 12/20/2022]
Abstract
Molecular modeling of large biomolecular assemblies exemplifies a disruptive area holding both promises and contentions. Propelled by peta and exascale computing, several simulation methodologies have now matured into user-friendly tools that are successfully employed for modeling viruses, membranous nano-constructs, and key pieces of the genetic machinery. We present three unifying biophysical themes that emanate from some of the most recent multi-million atom simulation endeavors. Despite connecting molecular changes with phenotypic outcomes, the quality measures of these simulations remain questionable. We discuss the existing and upcoming strategies for constructing representative ensembles of large systems, how new computing technologies will boost this area, and make a point that integrative modeling guided by experimental data is the future of biomolecular computations.
Collapse
Affiliation(s)
- Chitrak Gupta
- School of Molecular Sciences, Center for Applied Structural Discovery, Arizona State University at Tempe, Tempe, AZ, 85282, USA; Biodesign Institute, Tempe, AZ, 85281, USA. https://twitter.com/ChitrakGupta2
| | - Daipayan Sarkar
- School of Molecular Sciences, Center for Applied Structural Discovery, Arizona State University at Tempe, Tempe, AZ, 85282, USA; MSU-DOE Plant Research Laboratory, Michigan State University, East Lansing, MI, 48824-1319, USA. https://twitter.com/17Dsarkar
| | - D Peter Tieleman
- Centre for Molecular Simulation and Department of Biological Sciences, University of Calgary, Calgary, AB T2N 1N4, Canada.
| | - Abhishek Singharoy
- School of Molecular Sciences, Center for Applied Structural Discovery, Arizona State University at Tempe, Tempe, AZ, 85282, USA; Biodesign Institute, Tempe, AZ, 85281, USA.
| |
Collapse
|
9
|
Phan A, Stamatakis M, Koh CA, Striolo A. Correlating Antiagglomerant Performance with Gas Hydrate Cohesion. ACS APPLIED MATERIALS & INTERFACES 2021; 13:40002-40012. [PMID: 34382786 DOI: 10.1021/acsami.1c06309] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
Although inhibiting hydrate formation in hydrocarbon-water systems is paramount in preventing pipe blockage in hydrocarbon transport systems, the molecular mechanisms responsible for antiagglomerant (AA) performance are not completely understood. To better understand why macroscopic performance is affected by apparently small changes in the AA molecular structure, we perform molecular dynamics simulations. We quantify the cohesion energy between two gas hydrate nanoparticles dispersed in liquid hydrocarbons in the presence of different AAs, and we achieve excellent agreement against experimental data obtained at high pressure using the micromechanical force apparatus. This suggests that the proposed simulation approach could provide a screening method for predicting, in silico, the performance of new molecules designed to manage hydrates in flow assurance. Our results suggest that entropy and free energy of solvation of AAs, combined in some cases with the molecular orientation at hydrate-oil interfaces, are descriptors that could be used to predict performance, should the results presented here be reproduced for other systems as well. These insights could help speed up the design of new AAs and guide future experiments.
Collapse
Affiliation(s)
- Anh Phan
- Department of Chemical Engineering, University College London, London WC1E 7JE, U.K
| | - Michail Stamatakis
- Department of Chemical Engineering, University College London, London WC1E 7JE, U.K
| | - Carolyn A Koh
- Center for Hydrate Research, Chemical & Biological Engineering Department, Colorado School of Mines, Golden, Colorado 80401, United States
| | - Alberto Striolo
- Department of Chemical Engineering, University College London, London WC1E 7JE, U.K
| |
Collapse
|
10
|
Sharma P, Desikan R, Ayappa KG. Evaluating Coarse-Grained MARTINI Force-Fields for Capturing the Ripple Phase of Lipid Membranes. J Phys Chem B 2021; 125:6587-6599. [PMID: 34081861 DOI: 10.1021/acs.jpcb.1c03277] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Phospholipids, which are an integral component of cell membranes, exhibit a rich variety of lamellar phases modulated by temperature and composition. Molecular dynamics (MD) simulations have greatly enhanced our understanding of phospholipid membranes by capturing experimentally observed phases and phase transitions at molecular resolution. However, the ripple (Pβ') membrane phase, observed as an intermediate phase below the main gel-to-liquid crystalline transition with some lipids, has been challenging to capture with MD simulations, both at all-atom and coarse-grained (CG) resolutions. Here, with an aggregate ∼2.5 μs all-atom and ∼122 μs CGMD simulations, we systematically assess the ability of six CG MARTINI 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) lipid and water force-field (FF) variants, parametrized to capture the DPPC gel and fluid phases, for their ability to capture the Pβ' phase, and compared observations with those from an all-atom FF. Upon cooling from the fluid phase to below the phase transition temperature with smaller (380-lipid) and larger (>2200-lipid) MARTINI and all-atom (CHARMM36 FF) DPPC lipid bilayers, we observed that smaller bilayers with both all-atom and MARTINI FFs sampled interdigitated Pβ' and ripple-like states, respectively. However, while all-atom simulations of the larger DPPC membranes exhibited the formation of the Pβ' phase, MARTINI membranes did not sample interdigitated ripple-like states at larger system sizes. We then demonstrated that the ripple-like states in smaller MARTINI membranes were kinetically trapped structures caused by finite size effects rather than being representative of true Pβ' phases. We showed that a MARTINI FF variant that could capture the tilted Lβ' gel phase, a prerequisite for stabilizing the Pβ' phase, was unable to capture the rippled phase upon cooling. Our study reveals that the current MARTINI FFs (including MARTINI3) may require specific reparametrization of the interaction potentials to stabilize lipid interdigitation, a characteristic of the ripple phase.
Collapse
Affiliation(s)
- Pradyumn Sharma
- Department of Chemical Engineering, Indian Institute of Science, Bengaluru 560012, India
| | - Rajat Desikan
- Department of Chemical Engineering, Indian Institute of Science, Bengaluru 560012, India
| | - K Ganapathy Ayappa
- Department of Chemical Engineering, Indian Institute of Science, Bengaluru 560012, India.,Centre for BioSystems Science and Engineering, Indian Institute of Science, Bengaluru 560012, India
| |
Collapse
|