Lipid nanoparticles (LNPs) containing ionizable aminolipids are among the leading platforms for the successful delivery of nucleic-acid-based therapeutics, including messenger RNA (mRNA). The two recently FDA-approved COVID-19 vaccines developed by Moderna and Pfizer/BioNTech belong to this category. Ionizable aminolipids, cholesterol, and DSPC lipids are among the key components of such formulations, crucially modulating physicochemical properties of these formulations and, consequently, the potency of these therapeutics. Despite the importance of these components, the distribution of these molecules in LNPs containing mRNA is not clear. In this study, we used all-atom molecular dynamics (MD) simulations to investigate the distribution and effects of the Lipid-5 (apparent pKa of the lipid nanoparticle = 6.56), a rationally designed and previously reported ionizable aminolipid by Moderna, on lipid bilayers [Mol. Ther. 2018, 26, 1509–1519]. The simulations were conducted with half of the aminolipids charged and half neutral approximately to the expected ionization in the microenvironment of the LNP surface. In all five simulated systems in this work, the cholesterol content was kept constant, whereas the DSPC and Lipid-5 concentrations were changed systematically. We found that at higher concentrations of the ionizable aminolipids, the neutral aminolipids form a disordered aggregate in the membrane interior that preferentially includes cholesterol. The rules underlying the lipid redistribution could be used to rationally choose lipids to optimize the LNP function.
The first structures of α-synuclein (αSyn) fibrils have recently been solved. Here, we use a unique combination of molecular dynamics simulation strategies to address the minimal nucleation size of the 11-amino acid NAC protofibril solved by X-ray and to interrogate the dynamic behavior of unexpected crystal waters in the steric zipper. We found that protofibrils of >8 chains are thermodynamically stabilized due to protection of the fibril core from solvent influx and ordering of the end strands by the fibril core. In these stable oligomers, water molecules resolved in the crystal structure freely exchange with bulk solvent but are, on average, stably coordinated along the β-sheet by inward-facing Thr72 and Thr75. We confirm the persistence of this water coordination via simulations of the full-length Greek-key structure solved by NMR and speculate that these Thr-water networks are important in the context of enhanced fibril nucleation in the familial A53T mutation.
Endogenous and exogenously administered S-nitrosothiols modulate the activities of central and peripheral systems that control breathing. We have unpublished data showing that the deleterious effects of morphine on arterial blood-gas chemistry (i.e., pH, pCO2, pO2, and sO2) and Alveolar-arterial gradient (i.e., index of gas exchange) were markedly diminished in anesthetized Sprague Dawley rats that received a continuous intravenous infusion of the endogenous S-nitrosothiol, S-nitroso-L-cysteine. The present study extends these findings by showing that unanesthetized adult male Sprague Dawley rats receiving an intravenous infusion of S-nitroso-L-cysteine (100 or 200 nmol/kg/min) markedly diminished the ability of intravenous injections of the potent synthetic opioid, fentanyl (10, 25, and 50 μg/kg), to depress the frequency of breathing, tidal volume, and minute ventilation. Our study also found that the ability of intravenously injected fentanyl (10, 25, and 50 μg/kg) to disturb eupneic breathing, which was measured as a marked increase of the non-eupneic breathing index, was substantially reduced in unanesthetized rats receiving intravenous infusions of S-nitroso-L-cysteine (100 or 200 nmol/kg/min). In contrast, the deleterious effects of fentanyl (10, 25, and 50 μg/kg) on frequency of breathing, tidal volume, minute ventilation and non-eupneic breathing index were fully expressed in rats receiving continuous infusions (200 nmol/kg/min) of the parent amino acid, L-cysteine, or the D-isomer, namely, S-nitroso-D-cysteine. In addition, the antinociceptive actions of the above doses of fentanyl as monitored by the tail-flick latency assay, were enhanced by S-nitroso-L-cysteine, but not L-cysteine or S-nitroso-D-cysteine. Taken together, these findings add to existing knowledge that S-nitroso-L-cysteine stereoselectively modulates the detrimental effects of opioids on breathing, and opens the door for mechanistic studies designed to establish whether the pharmacological actions of S-nitroso-L-cysteine involve signaling processes that include 1) the activation of plasma membrane ion channels and receptors, 2) selective intracellular entry of S-nitroso-L-cysteine, and/or 3) S-nitrosylation events. Whether alterations in the bioavailability and bioactivity of endogenous S-nitroso-L-cysteine is a key factor in determining the potency/efficacy of fentanyl on breathing is an intriguing question.
Recent isothiocyanate covalent labeling studies have suggested that a classical cannabinoid, (-)-7'-isothiocyanato-11-hydroxy-1',1'dimethylheptyl-hexahydrocannabinol (AM841), enters the cannabinoid CB2 receptor via the lipid bilayer (Pei, Y., Mercier, R. W., Anday, J. K., Thakur, G. A., Zvonok, A. M., Hurst, D., Reggio, P. H., Janero, D. R., and Makriyannis, A. (2008) Chem. Biol. 15, 1207-1219). However, the sequence of steps involved in such a lipid pathway entry has not yet been elucidated. Here, we test the hypothesis that the endogenous cannabinoid sn-2-arachidonoylglycerol (2-AG) attains access to the CB2 receptor via the lipid bilayer. To this end, we have employed microsecond time scale all-atom molecular dynamics (MD) simulations of the interaction of 2-AG with CB2 via a palmitoyl-oleoyl-phosphatidylcholine lipid bilayer. Results suggest the following: 1) 2-AG first partitions out of bulk lipid at the transmembrane alpha-helix (TMH) 6/7 interface; 2) 2-AG then enters the CB2 receptor binding pocket by passing between TMH6 and TMH7; 3) the entrance of the 2-AG headgroup into the CB2 binding pocket is sufficient to trigger breaking of the intracellular TMH3/6 ionic lock and the movement of the TMH6 intracellular end away from TMH3; and 4) subsequent to protonation at D3.49/D6.30, further 2-AG entry into the ligand binding pocket results in both a W6.48 toggle switch change and a large influx of water. To our knowledge, this is the first demonstration via unbiased molecular dynamics that a ligand can access the binding pocket of a class A G protein-coupled receptor via the lipid bilayer and the first demonstration via molecular dynamics of G protein-coupled receptor activation triggered by a ligand binding event.
1. Abstract Arrestins are important scaffolding proteins that are expressed in all vertebrate animals. They regulate cell signaling events upon binding to active G-protein coupled receptors ( GPCR ) and trigger endocytosis of active GPCRs. While many of the functional sites on arrestins have been characterized, the question of how these sites interact is unanswered. We used anisotropic network modelling ( ANM ) together with our covariance compliment techniques to survey all of the available structures of the non-visual arrestins to map how structural changes and protein-binding affect their structural dynamics. We found that activation and clathrin binding have a marked effect on arrestin dynamics, and that these dynamics changes are localized to a small number of distant functional sites. These sites include α-helix 1, the lariat loop, nuclear localization domain, and the C-domain β-sheets on the C-loop side. Our techniques suggest that clathrin binding and/or GPCR activation of arrestin perturb the dynamics of these sites independent of structural changes.
