On the interactions between RNA and titrateable lipid layers: implications for RNA delivery with lipid nanoparticles.

Characterising the interaction between cationic ionisable lipids (CIL) and nucleic acids (NAs) is key to understanding the process of RNA lipid nanoparticle (LNP) formation and release of NAs from LNPs. Here, we have used different surface techniques to reveal the effect of pH and NA type on the interaction with a model system of DOPC and the CIL DLin-MC3-DMA (MC3). At only 5% MC3, differences in the structure and dynamics of the lipid layer were observed. Both pH and %MC3 were shown to affect the absorption behaviour of erythropoietin mRNA, polyadenylic acid (polyA) and polyuridylic acid (polyU). The adsorbed amount of all studied NAs was found to increase with decreasing pH and increasing %MC3 but with different effects on the lipid layer, which could be linked to the NA secondary structure. For polyA at pH 6, adsorption to the surface of the layer was observed, whereas for other conditions and NAs, penetration of the NA into the layer resulted in the formation of a multilayer structure. By comparison to simulations excluding the secondary structure, differences in adsorption behaviours between polyA and polyU could be observed, indicating that the NA's secondary structure also affected the MC3-NA interactions.

New insights into the protein stabilizing effects of trehalose by comparing with sucrose.

Disaccharides are well known to be efficient stabilizers of proteins, for example in the case of lyophilization or cryopreservation. However, although all disaccharides seem to exhibit bioprotective and stabilizing properties, it is clear that trehalose is generally superior compared to other disaccharides. The aim of this study was to understand this by comparing how the structural and dynamical properties of aqueous trehalose and sucrose solutions influence the protein myoglobin (Mb). The structural studies were based on neutron and X-ray diffraction in combination with empirical potential structure refinement (EPSR) modeling, whereas the dynamical studies were based on quasielastic neutron scattering (QENS) and molecular dynamics (MD) simulations. The results show that the overall differences in the structure and dynamics of the two systems are small, but nevertheless there are some important differences which may explain the superior stabilizing effects of trehalose. It was found that in both systems the protein is preferentially hydrated by water, but that this effect is more pronounced for trehalose, i.e. trehalose forms less hydrogen bonds to the protein surface than sucrose. Furthermore, the rotational motion around dihedrals between the two glucose rings of trehalose is slower than in the case of the dihedrals between the glucose and fructose rings of sucrose. This leads to a less perturbed protein structure in the case of trehalose. The observations indicate that an aqueous environment closest to the protein molecules is beneficial for an efficient bioprotective solution.

Ionizable lipids penetrate phospholipid bilayers with high phase transition temperatures: perspectives from free energy calculations.

The efficacies of modern gene-therapies strongly depend on their contents. At the same time the most potent formulations might not contain the best compounds. In this work we investigated the effect of phospholipids and their saturation on the binding ability of (6Z,9Z,28Z,31Z)-heptatriacont-6,9,28,31-tetraene-19-yl 4-(dimethylamino) butanoate (DLin-MC3-DMA) to model membranes at the neutral pH. We discovered that DLin-MC3-DMA has affinity to the most saturated monocomponent lipid bilayer 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and an aversion to the unsaturated one 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC). The preference to a certain membrane was also well-correlated to the phase transition temperatures of phospholipid bilayers, and to their structural and dynamical properties. Additionally, in the case of the presence of DLin-MC3-DMA in the membrane with DOPC the ionizable lipid penetrated it, which indicates possible synergistic effects. Comparisons with other ionizable lipids were performed using a model lipid bilayer of 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC). Particularly, the lipids heptadecan-9-yl 8-[2-hydroxyethyl-(6-oxo-6-undecoxyhexyl)amino]octanoate (SM-102) and [(4-hydroxybutyl) azanediyl] di(hexane-6,1-diyl) bis(2-hexyldecanoate) (ALC-0315) from modern mRNA-vaccines against COVID-19 were investigated and force fields parameters were derived for those new lipids. It was discovered that ALC-0315 binds strongest to the membrane, while DLin-MC3-DMA is not able to reside in the bilayer center. The ability to penetrate the membrane POPC by SM-102 and ALC-0315 can be related to their saturation, comparing to DLin-MC3-DMA.

Two statins and cromolyn as possible drugs against the cytotoxicity of Aß(31-35) and Aß(25-35) peptides: a comparative study by advanced computer simulation methods.

In this work, possible effective mechanisms of cromolyn, atorvastatin and lovastatin on the cytotoxicity of Aß(31-35) and Aß(25-35) peptides were investigated by classical molecular dynamics and well-tempered metadynamics simulations. The results demonstrate that all the drugs affect the behavior of the peptides, such as their ability to aggregate, and alter their secondary structures and their affinity to a particular drug. Our findings from the computed properties suggest that the best drug candidate is lovastatin. This medicine inhibits peptide aggregation, adsorbs the peptides on the surface of the drug clusters, changes the secondary structure and binds to MET35, which has been seen as the reason for the toxicity of the studied peptide sequences. Moreover, lovastatin is the drug which previously has demonstrated the strongest ability to penetrate the blood-brain barrier and makes lovastatin the most promising medicine among the three investigated drugs. Atorvastatin is also seen as a potential candidate if its penetration through the blood-brain barrier could be improved. Otherwise, its properties are even better than the ones demonstrated by lovastatin. Cromolyn appears to be less interesting as an anti-aggregant from the computational data, in comparison to the two statins.

Effect of encapsulated protein on the dynamics of lipid sponge phase: a neutron spin echo and molecular dynamics simulation study.

Lipid membranes are highly mobile systems with hierarchical, time and length scale dependent, collective motions including thickness fluctuations, undulations, and topological membrane changes, which play an important role in membrane interactions. In this work we have characterised the effect of encapsulating two industrially important enzymes, ß-galactosidase and aspartic protease, in lipid sponge phase nanoparticles on the dynamics of the lipid membrane using neutron spin echo (NSE) spectroscopy and molecular dynamics (MD) simulations. From NSE, reduced membrane dynamics were observed upon enzyme encapsulation, which were dependent on the enzyme concentration and type. By fitting the intermediate scattering functions (ISFs) with a modified Zilman and Granek model including nanoparticle diffusion, an increase in membrane bending rigidity was observed, with a larger effect for ß-galactosidase than aspartic protease at the same concentration. MD simulations for the system with and without aspartic protease showed that the lipids relax more slowly in the system with protein due to the replacement of the lipid carbonyl-water hydrogen bonds with lipid-protein hydrogen bonds. This indicates that the most likely cause of the increase in membrane rigidity observed in the NSE measurements was dehydration of the lipid head groups. The dynamics of the protein itself were also studied, which showed a stable secondary structure of protein over the simulation, indicating no unfolding events occurred.

Atomistic molecular dynamics simulations of tubulin heterodimers explain the motion of a microtubule.

Microtubules are essential parts of the cytoskeleton that are built by polymerization of tubulin heterodimers into a hollow tube. Regardless that their structures and functions have been comprehensively investigated in a modern soft matter, it is unclear how properties of tubulin heterodimer influence and promote the self-assembly. A detailed knowledge of such structural mechanisms would be helpful in drug design against neurodegenerative diseases, cancer, diabetes etc. In this work atomistic molecular dynamics simulations were used to investigate the fundamental dynamics of tubulin heterodimers in a sheet and a short microtubule utilizing well-equilibrated structures. The breathing motions of the tubulin heterodimers during assembly show that the movement at the lateral interface between heterodimers (wobbling) dominates in the lattice. The simulations of the protofilament curvature agrees well with recently published experimental data, showing curved protofilaments at polymerization of the microtubule plus end. The tubulin heterodimers exposed at the microtubule minus end were less curved and displayed altered interactions at the site of sheet closure around the outmost heterodimers, which may slow heterodimer binding and polymerization, providing a potential explanation for the limited dynamics observed at the minus end.

Component of Cannabis, Cannabidiol, as a Possible Drug against the Cytotoxicity of Aß(31-35) and Aß(25-35) Peptides: An Investigation by Molecular Dynamics and Well-Tempered Metadynamics Simulations.

In this work cannabidiol (CBD) was investigated as a possible drug against the cytotoxicity of Aß(31-35) and Aß(25-35) peptides with the help of atomistic molecular dynamics (MD) and well-tempered metadynamics simulations. Four interrelated mechanisms of possible actions of CBD are proposed from our computations. This implies that one mechanism can be a cause or/and a consequence of another. CBD is able to decrease the aggregation of peptides at certain concentrations of compounds in water. This particular action is more prominent for Aß(25-35), since originally Aß(31-35) did not exhibit aggregation properties in aqueous solutions. Interactions of CBD with the peptides affect secondary structures of the latter ones. Clusters of CBD are seen as possible adsorbents of Aß(31-35) and Aß(25-35) since peptides are tending to aggregate around them. And last but not least, CBD exhibits binding to MET35. All four mechanisms of actions can possibly inhibit the Aß-cytotoxicity as discussed in this paper. Moreover, the amount of water also played a role in peptide clustering: with a growing concentration of peptides in water without a drug, the aggregation of both Aß(31-35) and Aß(25-35) increased. The number of hydrogen bonds between peptides and water was significantly higher for simulations with Aß(25-35) at the higher concentration of peptides, while for Aß(31-35) that difference was rather insignificant. The presence of CBD did not substantially affect the number of hydrogen bonds in the simulated systems.

DOPC versus DOPE as a helper lipid for gene-therapies: molecular dynamics simulations with DLin-MC3-DMA.

Ionizable lipids are important compounds of modern therapeutic lipid nano-particles (LNPs). One of the most promising ionizable lipids (or amine lipids) is DLin-MC3-DMA. Depending on their pharmaceutical application these LNPs can also contain various helper lipids, such as phospho- and pegylated lipids, cholesterol and nucleic acids as a cargo. Due to their complex compositions the structures of these therapeutics have not been refined properly. Therefore, the role of each lipid in the pharmacological properties of LNPs has not been determined. In this work an atomistic model for the neutral form of DLin-MC3-DMA was derived and all-atom molecular dynamics (MD) simulations were carried out in order to investigate the effect of the phospholipid headgroup on the possible properties of the shell-membranes of LNPs. Bilayers containing either DOPC or DOPE lipids at two different ratios of DLin-MC3-DMA (5 mol% and 15 mol%) were constructed and simulated at neutral pH 7.4. The results from the analysis of MD trajectories revealed that DOPE lipid headgroups associated strongly with lipid tails and carbonyl oxygens of DLin-MC3-DMA, while for DOPC lipid headgroups no significant associations were observed. Furthermore, the strong associations between DOPE and DLin-MC3-DMA result in the positioning of DLin-MC3-DMA at the surface of the membrane. Such an interplay between the lipids slows down the lateral diffusion of all simulated bilayers, where a more dramatic decrease of the diffusion rate is observed in membranes with DOPE. This can explain the low water penetration of lipid bilayers with phosphatidylethanolamines and, probably, can relate to the bad transfection properties of LNPs with DOPE and DLin-MC3-DMA.

Modelling of interactions between Aß(25-35) peptide and phospholipid bilayers: effects of cholesterol and lipid saturation.

Aggregation of amyloid beta (Aß) peptides in neuronal membranes is a known promoter of Alzheimer's disease. To gain insight into the molecular details of Aß peptide aggregation and its effect on model neuronal membranes, we carried out molecular dynamics simulations of the Aß(25-35) fragment of the amyloid precursor protein in phospholipid bilayers composed of either fully saturated or highly unsaturated lipids, in the presence or absence of cholesterol. It was found that the peptide does not penetrate through any of the considered membranes, but can reside in the headgroup region and upper part of the lipid tails showing a clear preference to a polyunsaturated cholesterol-free membrane. Due to the ordering and condensing effect upon addition of cholesterol, membranes become more rigid facilitating peptide aggregation on the surface. Except for the case of the cholesterol-free saturated lipid bilayer, the peptides have a small effect on the membrane structure and ordering. It was also found that the most "active" amino-acid for peptide-lipid and peptide-cholesterol interaction is methionine-35, followed by asparagine-27 and serine-26, which form hydrogen bonds between peptides and polar atoms of lipid headgroups. These amino acids are also primarily responsible for peptide aggregation. This work will be relevant for designing strategies to develop drugs to combat Alzheimer's disease.

Effect of lipid saturation on amyloid-beta peptide partitioning and aggregation in neuronal membranes: molecular dynamics simulations.

Aggregation of amyloid-[Formula: see text] (Aß) peptides, cleaved from the amyloid precursor protein, is known as a precursor of the Alzheimer's disease (AD). It is also known that Alzheimer's disease is characterized by a substantial decrease of the amount of polyunsaturated lipids in the neuronal membranes of the frontal gray matter. To get insight into possible interconnection of these phenomena, we have carried out molecular dynamics simulations of two fragments of A[Formula: see text] peptide, A[Formula: see text][Formula: see text] and A[Formula: see text][Formula: see text], in four different lipid bilayers: two monocomponent ones (14:0-14:0 PC, 18:0-22:6 PC), and two bilayers containing mixtures of 18:0-18:0 PE, 22:6-22:6 PE, 16:0-16:0 PC and 18:1-18:1 PC lipids of composition mimicking neuronal membranes in a "healthy" and "AD" brain. The simulations showed that the presence of lipids with highly unsaturated 22:6cis fatty acids chains strongly affects the interaction of amyloid-[Formula: see text] peptides with lipid membranes. The polyunsaturated lipids cause stronger adsorption of A[Formula: see text]-peptides by the membrane and lead to weaker binding between peptides when the latter form aggregates. This difference in the behaviour observed in monocomponent bilayers is propagated in a similar fashion to the mixed membranes mimicking composition of neuronal membranes in "healthy" and "AD" brains, with "healthy" membrane having higher fraction of polyunsaturated lipids. Our simulations give strong indication that it can be physical-chemical background of the interconnection between amyloid fibrillization causing Alzheimer's disease, and content of polyunsaturated lipids in the neuronal membranes.

Cholesterol in phospholipid bilayers: positions and orientations inside membranes with different unsaturation degrees.

Cholesterol is an essential component of all animal cell membranes and plays an important role in maintaining the membrane structure and physical-chemical properties necessary for correct cell functioning. The presence of cholesterol is believed to be responsible for domain formation (lipid rafts) due to different interactions of cholesterol with saturated and unsaturated lipids. In order to get detailed atomistic insight into the behaviour of cholesterol in bilayers composed of lipids with varying degrees of unsaturation, we have carried out a series of molecular dynamics simulations of saturated and polyunsaturated lipid bilayers with different contents of cholesterol, as well as well-tempered metadynamics simulations with a single cholesterol molecule in these bilayers. From these simulations we have determined distributions of cholesterol across the bilayer, its orientational properties, free energy profiles, and specific interactions of molecular groups able to form hydrogen bonds. Both molecular dynamics and metadynamics simulations showed that the most unsaturated bilayer with 22:6 fatty acid chains shows behaviour which is most different from other lipids. In this bilayer, cholesterol is relatively often found in a "flipped" configuration with the hydroxyl group oriented towards the membrane middle plane. This bilayer has also the highest (least negative) binding free energy among liquid phase bilayers, and the lowest reorientation barrier. Furthermore, cholesterol molecules in this bilayer are often found to form head-to-tail contacts which may lead to specific clustering behaviour. Overall, our simulations support ideas that there can be a subtle interconnection between the contents of highly unsaturated fatty acids and cholesterol, deficiency or excess of each of them is related to many human afflictions and diseases.

Molecular Dynamics Simulations of Furfural and 5-Hydroxymethylfurfural at Ambient and Hydrothermal Conditions.

In this work, we present results from molecular dynamics simulations of aqueous solutions of furfural and 5-hydroxymethylfurfural, which are important intermediates in the hydrothermal carbonization processes of biomass conversion. The computations were performed both at ambient and hydrothermal conditions using a two-level factorial design varying concentration, temperature, and pressure. A number of equilibrium and dynamic properties have been computed including enthalpies and free energies of vaporization, free energies of solvation, diffusion coefficients, and rotational/reorientational correlation times. Structural properties of solutions were analyzed using radial and spatial distribution functions. It was shown that the formation of hydrogen bonds among 5-hydroxymethylfurfural molecules is preferred compared to hydrogen bonding between 5-hydroxymethylfurfural and water. In addition, our results suggest that the oxygen atoms in the furan rings of furfural and 5-hydroxymethylfurfural do not participate in hydrogen bonding to the same extent as the oxygen atoms in the hydroxyl and carbonyl groups. It is also observed that furfural molecules aggregate under certain conditions, and we show how this is affected by changes in temperature, pressure, and concentration in agreement with experimental solubility data. The analysis of the computational results provides useful insight into the structure and dynamics of the considered molecules at conditions of hydrothermal carbonization, as well as at ambient conditions.

Quantum chemical and molecular dynamics modelling of hydroxylated polybrominated diphenyl ethers.

A series of 19 hydroxylated polybrominated diphenyl ethers (OH-PBDEs) have been studied using density functional theory (DFT) and molecular dynamics simulations with the purpose of investigating eventual correlations between their physicochemical properties and toxic action. Dissociation constants (pKa), solvation free energies and octanol-water partition coefficients (log P) have been computed. Additionally, metadynamics simulations of OH-PBDEs passing through a lipid bilayer have been carried out for four OH-PBDE species. No correlations between computed pKa values and toxicity data have been found. Medium correlations were found between partition coefficients and the ability of OH-PBDEs to alter membrane potential in cell cultures, which is attributed to higher uptake of molecules with larger log P parameters. It was also demonstrated that in lipid bilayers, OH-PBDE molecules differ in their orientational distributions and can adopt different conformations which can affect the uptake of these molecules and influence the pathways of their toxic action.

Extension of the Slipids Force Field to Polyunsaturated Lipids.

The all-atomic force field Slipids (Stockholm Lipids) for lipid bilayers simulations has been extended to polyunsaturated lipids. Following the strategy adopted in the development of previous versions of the Slipids force field, the parametrization was essentially based on high-level ab initio calculations. Atomic charges and torsion angles related to polyunsaturated lipid tails were parametrized using structures of dienes molecules. The new parameters of the force field were validated in simulations of bilayers composed of seven polyunsaturated lipids. An overall good agreement was found with available experimental data on the areas per lipids, volumetric properties of bilayers, deuterium order parameters, and scattering form factors. Furthermore, simulations of bilayers consisting of highly polyunsaturated lipids and cholesterol molecules have been carried out. The majority of cholesterol molecules were found in a position parallel to bilayer normal with the hydroxyl group directed to the membrane surface, while a small fraction of cholesterol was found in the bilayer center parallel to the membrane plane. Furthermore, a tendency of cholesterol molecules to form chain-like clusters in polyunsaturated bilayers was qualitatively observed.