Fluid Phase Coexistence in Biological Membrane: Insights from Local Nonaffine Deformation of Lipids

Membrane Catalyzed Formation of Nucleotide Clusters and Their Role in the Origins of Life: Insights from Molecular Simulations and Lattice Modeling

One of the mysteries in studying the molecular "Origin of Life" is the emergence of RNA and RNA-based life forms, where nonenzymatic polymerization of nucleotides is a crucial hypothesis in formation of large RNA chains. The nonenzymatic polymerization can be mediated by various environmental settings, such as cycles of hydration and dehydration, temperature variations, and proximity to a variety of organizing matrices, such as clay, salt, fatty acids, lipid membrane, and mineral surface. In this work, we explore the influence of different phases of the lipid membrane toward nucleotide organization and polymerization in a simulated prebiotic setting. Our molecular simulations quantify the localization propensity of a mononucleotide, uridine monophosphate (UMP), in distinct membrane settings. We perform all-atom molecular dynamics (MD) simulations to estimate the role of the monophasic and biphasic membranes in modifying the behavior of UMPs localization and their clustering mechanism. Based on the interaction energy of mononucleotides with the membrane and their diffusion profile from our MD calculations, we developed a lattice-based model to explore the thermodynamic limits of the observations made from the MD simulations. The mathematical model substantiates our hypothesis that the lipid layers can act as unique substrates for "catalyzing" polymerization of mononucleotides due to the inherent spatiotemporal heterogeneity and phase change behavior.

Enthalpic and entropic contributions to interleaflet coupling drive domain registration and antiregistration in biological membrane

Biological membrane is a complex self-assembly of lipids, sterols, and proteins organized as a fluid bilayer of two closely stacked lipid leaflets. Differential molecular interactions among its diverse constituents give rise to heterogeneities in the membrane lateral organization. Under certain conditions, heterogeneities in the two leaflets can be spatially synchronized and exist as registered domains across the bilayer. Several contrasting theories behind mechanisms that induce registration of nanoscale domains have been suggested. Following a recent study showing the effect of position of lipid tail unsaturation on domain registration behavior, we decided to develop an analytical theory to elucidate the driving forces that create and maintain domain registry across leaflets. Towards this, we formulated a Hamiltonian for a stacked lattice system where site variables capture the lipid molecular properties such as the position of unsaturation and various other interactions that could drive phase separation and interleaflet coupling. We solve the Hamiltonian using Monte Carlo simulations and create a complete phase diagram that reports the presence or absence of registered domains as a function of various Hamiltonian parameters. We find that the interleaflet coupling should be described as a competing enthalpic contribution due to interaction of lipid tail termini, primarily due to saturated-saturated interactions, and an interleaflet entropic contribution from overlap of unsaturated tail termini. A higher position of unsaturation is seen to provide weaker interleaflet coupling. Thermodynamically stable nanodomains could also be observed for certain points in the parameter space in our bilayer model, which were further verified by carrying out extended Monte Carlo simulations. These persistent noncoalescing registered nanodomains close to the lower end of the accepted nanodomain size range also point towards a possible "nanoscale" emulsion description of lateral heterogeneities in biological membrane leaflets.

Modulation of Phospholipid Bilayer Properties by Simvastatin

Simvastatin (Zocor) is one of the most prescribed drugs for reducing high cholesterol. Although simvastatin is ingested in its inactive lactone form, it is converted to its active dihydroxyheptanoate form by carboxylesterases in the liver. The dihydroxyheptanoate form can also be converted back to its original lactone form. Unfortunately, some of the side effects associated with the intake of simvastatin and other lipophilic statins at higher doses include statin-associated myopathy (SAM) and, in more severe cases, kidney failure. While the cause of SAM is unknown, it is hypothesized that these side effects are dependent on the localization of statins in lipid bilayers and their impact on bilayer properties. In this work, we carry out all-atom molecular dynamics simulations on both the lactone and dihydroxyheptanoate forms of simvastatin (termed "SN" and "SA", respectively) with a pure 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid bilayer and a POPC/cholesterol (30 mol %) binary mixture as membrane models. Additional simulations were carried out with multiple simvastatin molecules to mimic in vitro conditions that produced pleiotropic effects. Both SN and SA spontaneously diffused into the lipid bilayer, and a longer simulation time of 4 μs was needed for the complete incorporation of multiple SAs into the bilayer. By constructing potential mean force and electron density profiles, we find that SN localizes deeper within the hydrophobic interior of the bilayer and that SA has a greater tendency to form hydrogen-bonding interactions with neighboring water molecules and lipid headgroups. For the pure POPC bilayer, both SN and SA increase membrane order, while membrane fluidity increases for the POPC/cholesterol bilayer.

Structural insights into fusion mechanisms of small extracellular vesicles with model plasma membranes

Extracellular vesicles (EVs) are a potent intercellular communication system. Such small vesicles transport biomolecules between cells and throughout the body, strongly influencing the fate of recipient cells. Due to their specific biological functions they have been proposed as biomarkers for various diseases and as optimal candidates for therapeutic applications. Despite their extreme biological relevance, their mechanisms of interaction with the membranes of recipient cells are still hotly debated. Here, we propose a multiscale investigation based on atomic force microscopy, small angle X-ray scattering, small angle neutron scattering and neutron reflectometry to reveal structure-function correlations of purified EVs in interaction with model membrane systems of variable complex compositions and to spot the role of different membrane phases on the vesicle internalization routes. Our analysis reveals strong interactions of EVs with the model membranes and preferentially with the borders of protruding phase domains. Moreover, we found that upon vesicle breaking on the model membrane surface, the biomolecules carried by/on EVs diffuse with different kinetics rates, in a process distinct from simple fusion. The biophysical platform proposed here has clear implications on the modulation of EV internalization routes by targeting specific domains at the plasma cell membrane and, as a consequence, on EV-based therapies.

Three-Dimensional Packing Defects in Lipid Membrane as a Function of Membrane Order

Lipid membrane packing defects are considered to be an essential parameter that regulates specific membrane binding of several peripheral proteins. In the absence of direct experimental characterization, lipid packing defects and their role in the binding of peripheral proteins are generally investigated through computational studies, which have been immensely successful in unraveling the key steps of the membrane-binding process. However, packing defects are calculated using two-dimensional (2D) projections and the crucial information on their depths is generally overlooked. Here, we present a simple yet computationally efficient algorithm, which identifies these defects in three dimensions. We validate the algorithm on a number of model membrane systems that are previously studied using 2D defect calculations and find that the defect size and the defect depth may not always be directly correlated. We employ the algorithm to understand the nature of packing defects in flat bilayer membranes exhibiting liquid-ordered (L_o), liquid-disordered (L_d), and co-existing (L_o/L_d) phases. Our results indicate the presence of shallower, smaller, and spatially localized defects in the L_o phase membranes as compared to the defects in L_d and mixed L_o/L_d phase membranes. Such analyses can elucidate the molecular-scale mechanisms that drive the preferential localization of certain proteins to either of the liquid phases or their interface. We also analyze the membrane sensing and anchoring process of a peptide in terms of the three-dimensional defects, which provides additional insights into the process with respect to depth distributions across the bilayer leaflets.

Degeneracy in molecular scale organization of biological membranes

The scale-rich spatiotemporal organization in biological membranes has its origin in the differential inter- and intra-molecular interactions among their constituents. In this work, we explore the molecular-origin behind that variety and possible degeneracy in lateral organization in membranes. For our study, we post-process microsecond long all-atom molecular dynamics trajectories for three systems that exhibit fluid phase coexistence: (i) PSM/POPC/Chol (0.47/0.32/0.21), (ii) PSM/DOPC/Chol (0.43/0.38/0.19) and (iii) DPPC/DOPC/Chol (0.37/0.36/0.27). To distinguish the liquid ordered and disordered regions at molecular scales, we calculate the degree of non-affineness of individual lipids in their neighbourhood and track their topological rearrangements. Disconnectivity graph analysis with respect to membrane organization shows that the DPPC/DOPC/Chol and PSM/DOPC/Chol systems exhibit funnel-like energy landscapes as opposed to a highly frustrated energy landscape for the more biomimetic PSM/POPC/Chol system. We use these measurements to develop a continuous lattice Hamiltonian and evolve that using Monte Carlo simulated annealing to explore the possibility of structural degeneracy in membrane organization. Our data show that model membranes with lipid constituents that are biomimetic (PSM/POPC/Chol) have the ability to access a large range of membrane sub-structure space (higher degeneracy) as compared to the other two systems, which form only one kind of substructure even with changing composition. Since the spatiotemporal organization in biological membranes dictates the "molecular encounters" and in turn larger scale biological processes such as molecular transport, trafficking and cellular signalling, we posit that this structural degeneracy could enable access to a larger repository to functionally important molecular organization in systems with physiologically relevant compositions.

Non-affinity in multi-material mechanical metamaterials

Non-affine deformations enable mechanical metamaterials to achieve their unusual properties while imposing implications for their structural integrity. The presence of multiple phases with different mechanical properties results in additional non-affinity of the deformations, a phenomenon that has never been studied before in the area of extremal mechanical metamaterials. Here, we studied the degree of non-affinity, [Formula: see text], resulting from the random substitution of a fraction of the struts,[Formula: see text], that make up a lattice structure and are printed using a soft material (elastic modulus = [Formula: see text]) by those printed using a hard material ([Formula: see text]). Depending on the unit cell angle (i.e., [Formula: see text] = 60°, 90°, or 120°), the lattice structures exhibited negative, near-zero, or positive values of the Poisson's ratio, respectively. We found that the auxetic structures exhibit the highest levels of non-affinity, followed by the structures with positive and near-zero values of the Poisson's ratio. We also observed an increase in [Formula: see text] with [Formula: see text] and [Formula: see text] until [Formula: see text] =10⁴ and [Formula: see text]= 75%-90% after which [Formula: see text] saturated. The dependency of [Formula: see text] upon [Formula: see text] was therefore found to be highly asymmetric. The positive and negative values of the Poisson's ratio were strongly correlated with [Formula: see text]. Interestingly, achieving extremely high or extremely low values of the Poisson's ratio required highly affine deformations. In conclusion, our results clearly show the importance of considering non-affinity when trying to achieve a specific set of mechanical properties and underscore the structural integrity implications in multi-material mechanical metamaterials.

Untangling Direct and Domain-Mediated Interactions Between Nicotinic Acetylcholine Receptors in DHA-Rich Membranes

At the neuromuscular junction (NMJ), the nicotinic acetylcholine receptor (nAChR) self-associates to give rise to rapid muscle movement. While lipid domains have maintained nAChR aggregates in vitro, their specific roles in nAChR clustering are currently unknown. In the present study, we carried out coarse-grained molecular dynamics simulations (CG-MD) of 1-4 nAChR molecules in two membrane environments: one mixture containing domain-forming, homoacidic lipids, and a second mixture consisting of heteroacidic lipids. Spontaneous dimerization of nAChRs was up to ten times more likely in domain-forming membranes; however, the effect was not significant in four-protein systems, suggesting that lipid domains are less critical to nAChR oligomerization when protein concentration is higher. With regard to lipid preferences, nAChRs consistently partitioned into liquid-disordered domains occupied by the omega-3 ([Formula: see text]-3) fatty acid, docosahexaenoic acid (DHA); enrichment of DHA boundary lipids increased with protein concentration, particularly in homoacidic membranes. This result suggests dimer formation blocks access of saturated chains and cholesterol, but not polyunsaturated chains, to boundary lipid sites.

Boundary lipids of the nicotinic acetylcholine receptor: Spontaneous partitioning via coarse-grained molecular dynamics simulation

Reconstituted nicotinic acetylcholine receptors (nAChRs) exhibit significant gain-of-function upon addition of cholesterol to reconstitution mixtures, and cholesterol affects the organization of nAChRs within domain-forming membranes, but whether nAChR partitions to cholesterol-rich liquid-ordered ("raft" or l_o) domains or cholesterol-poor liquid-disordered (l_do) domains is unknown. We use coarse-grained molecular dynamics simulations to observe spontaneous interactions of cholesterol, saturated lipids, and polyunsaturated (PUFA) lipids with nAChRs. In binary Dipalmitoylphosphatidylcholine:Cholesterol (DPPC:CHOL) mixtures, both CHOL and DPPC acyl chains were observed spontaneously entering deep "non-annular" cavities in the nAChR TMD, particularly at the subunit interface and the β subunit center, facilitated by the low amino acid density in the cryo-EM structure of nAChR in a native membrane. Cholesterol was highly enriched in the annulus around the TMD, but this effect extended over (at most) 5-10 Å. In domain-forming ternary mixtures containing PUFAs, the presence of a single receptor did not significantly affect the likelihood of domain formation. nAChR partitioned to any cholesterol-poor l_do domain that was present, regardless of whether the l_do or l_o domain lipids had PC or PE headgroups. Enrichment of PUFAs among boundary lipids was positively correlated with their propensity for demixing from cholesterol-rich phases. Long n-3 chains (tested here with Docosahexaenoic Acid, DHA) were highly enriched in annular and non-annular embedded sites, partially displacing cholesterol and completely displacing DPPC, and occupying sites even deeper within the bundle. Shorter n-6 chains were far less effective at displacing cholesterol from non-annular sites.