Compatibility of advanced water models with a united atom model of lipid in lipid bilayer simulation

Anomalous lateral diffusion of lipids during the fluid/gel phase transition of a lipid membrane

A lipid membrane undergoes a phase transition from fluid to gel phase upon changing external thermodynamic conditions, such as decreasing temperature and increasing pressure. Extremophilic organisms face the challenge of preventing this deleterious phase transition. The main focus of their adaptive strategy is to facilitate effective temperature sensing through sensor proteins, relying on the drastic changes in packing density and membrane fluidity during the phase transition. Although the changes in packing density parameters due to the fluid/gel phase transition are studied in detail, the impact on membrane fluidity is less explored in the literature. Understanding the lateral diffusive dynamics of lipids in response to temperature, particularly during the fluid/gel phase transition, is albeit crucial. Here we have simulated the phase transition of a single component lipid membrane composed of dipalmitoylphosphatidylcholine (DPPC) lipids using a coarse-grained (CG) model and studied the changes of the structural and dynamical properties. It is observed that near the phase transition point, both fluid and gel phase domains coexist together. The dynamics remains highly non-Gaussian for a long time even when the mean square displacement reaches the Fickian regime at a much earlier time. This Fickian yet non-Gaussian diffusion (FnGD) is a characteristic of a highly heterogeneous system, previously observed for the lateral diffusion of lipids in raft mimetic membranes having liquid-ordered and liquid-disordered phases co-existing together. We have analyzed the molecular trajectories and calculated the jump-diffusion of the lipids, stemming from sudden jump translations, using a translational jump-diffusion (TJD) approach. An overwhelming contribution of the jump-diffusion of the lipids is observed suggesting anomalous diffusion of lipids during fluid/gel phase transition of the membrane. These results are important in unravelling the intricate nature of lipid diffusion during the phase transition of the membrane and open up a new possibility of investigating the most significant change of membrane properties during phase transition, which can be effectively sensed by proteins.

Carbon Nanodots from an In Silico Perspective

Carbon nanodots (CNDs) are the latest and most shining rising stars among photoluminescent (PL) nanomaterials. These carbon-based surface-passivated nanostructures compete with other related PL materials, including traditional semiconductor quantum dots and organic dyes, with a long list of benefits and emerging applications. Advantages of CNDs include tunable inherent optical properties and high photostability, rich possibilities for surface functionalization and doping, dispersibility, low toxicity, and viable synthesis (top-down and bottom-up) from organic materials. CNDs can be applied to biomedicine including imaging and sensing, drug-delivery, photodynamic therapy, photocatalysis but also to energy harvesting in solar cells and as LEDs. More applications are reported continuously, making this already a research field of its own. Understanding of the properties of CNDs requires one to go to the levels of electrons, atoms, molecules, and nanostructures at different scales using modern molecular modeling and to correlate it tightly with experiments. This review highlights different in silico techniques and studies, from quantum chemistry to the mesoscale, with particular reference to carbon nanodots, carbonaceous nanoparticles whose structural and photophysical properties are not fully elucidated. The role of experimental investigation is also presented. Hereby, we hope to encourage the reader to investigate CNDs and to apply virtual chemistry to obtain further insights needed to customize these amazing systems for novel prospective applications.

Molecular dynamics study on the effects of charged amino acid distribution under low pH condition to the unfolding of hen egg white lysozyme and formation of beta strands

Aggregation of unfolded or misfolded proteins into amyloid fibrils can cause various diseases in humans. However, the fibrils synthesized in vitro can be developed toward useful biomaterials under some physicochemical conditions. In this study, atomistic molecular dynamics simulations were performed to address the mechanism of beta-sheet formation of the unfolded hen egg-white lysozyme (HEWL) under a high temperature and low pH. Simulations of the protonated HEWL at pH 2 and the non-protonated HEWL at pH 7 were performed at the highly elevated temperature of 450 K to accelerate the unfolding, followed by the 333 K temperature to emulate some previous in vitro studies. The simulations showed that HEWL unfolded faster, and higher beta-strand contents were observed at pH 2. In addition, one of the simulation replicas at pH 2 showed that the beta-strand forming sequence was consistent with the ‘K-peptide’, proposed as the core region for amyloidosis in previous experimental studies. Beta-strand formation mechanisms at the earlier stage of amyloidosis were explained in terms of the radial distribution of the amino acids. The separation between groups of positively charged sidechains from the hydrophobic core corresponded to the clustering of the hydrophobic residues and beta-strand formation.

Can Urea and Trimethylamine-N-oxide Prevent the Pressure-Induced Phase Transition of Lipid Membrane?

Organisms dwelling in ocean trenches are exposed to the high hydrostatic pressure of ocean water. Increasing pressure can alter the membrane packing density and fluidity and trigger the fluid-to-gel phase transition. To combat environmental stress, the organisms synthesize small polar solutes, which are known as osmolytes. Urea and trimethylamine-N-oxide (TMAO) are two such solutes found in deep-sea creatures. While TMAO stabilizes protein, urea induces protein denaturation. These solutes strongly influence the packing density and membrane fluidity of the lipid bilayer at different conditions. But can these solutes affect the pressure-induced phase transition of the lipid membrane? In the present work, we have studied the effect of these two solutes on pressure-induced fluid-to-gel phase transition based on the all-atom molecular dynamics (MD) simulation approach. A high-pressure-stimulated fluid-to-gel phase transition of the membrane is seen at 800 bar, which is consistent with previous experiments. We have also observed that in the low-pressure region (1-400 bar), urea slightly increases the membrane fluidity where TMAO decreases the same. However, the phase transition pressure remains almost unchanged on the addition of urea while TMAO shifts the phase transition toward a lower pressure. We have found that the hydrogen (H)-bond interaction between lipid and urea plays an important role in preserving the fluidity of the membrane in the low-pressure zone. However, at a higher pressure, both water and urea are excluded from the membrane surface. TMAO is also excluded from the interfacial region of the membrane at all pressures. Exclusion from the membrane surface further triggers the phase transition of the lipid membrane from the fluid to gel phase at a high pressure.

How does excess phenylalanine affect the packing density and fluidity of a lipid membrane?

Phenylketonuria (PKU) is an autosomal recessive error of phenylalanine (Phe) metabolism, where untreated Phe becomes cytotoxic. Previous experiments found that excess Phe decreases the packing density and increases the fluidity and permeability of a lipid membrane. It was proposed that Phe forms cytotoxic nanoscopic amyloid-like fibrils. In another study, the Phe fibrils were not visible near the lipid membrane. So, what leads to the deleterious effect of Phe on the lipid membrane? We put forward a molecular mechanism for the observed effect of excess Phe on the lipid membrane using all-atom molecular dynamics simulation. This study suggests that Phe monomers spontaneously intercalate into the membrane and form small hydrogen-bonded clusters, some of which locally perturb the membrane. These local effects result in an overall reduction in the membrane packing density, enhancement of membrane fluidity, and an increase of water permeability, observed in experiments. The present study does not observe any effect of the nanoscopic fibrillar structure of Phe on the membrane. This study, therefore, provides alternative insights into the excess Phe cytotoxicity in PKU disease.

How Do Urea and Trimethylamine N-Oxide Influence the Dehydration-Induced Phase Transition of a Lipid Membrane?

Living organisms are often exposed to extreme dehydration, which is detrimental to the structure and function of the cell membrane. The lipid membrane undergoes fluid-to-gel phase transition due to dehydration and thus loses fluidity and functionality. To protect the fluid phase of the bilayer these organisms adopt several strategies. Enhanced production of small polar organic solutes (also called osmolytes) is one such strategy. Urea and trimethylamine N-oxide (TMAO) are two osmolytes found in different organisms combating osmotic stress. Previous experiments have found that both these osmolytes have strong effects on lipid membrane under different hydration conditions. Urea prevents the dehydration-induced phase transition of the lipid membrane by directly interacting with the lipids, while TMAO does not inhibit the phase transition. To provide atomistic insights, we have carried out all-atom molecular dynamics (MD) simulation of a lipid membrane under varying hydration levels and studied the effect of these osmolytes on different structural and dynamic properties of the membrane. This study suggests that urea significantly inhibits the dehydration-induced fluid-to-gel phase transition by strongly interacting with the lipid membrane via hydrogen bonds, which balances the reduced lipid hydration due to the decreasing water content. In contrast, TMAO is excluded from the membrane surface due to unfavorable interaction with the lipids. This induces further dehydration of the lipids which reinforces the fluid-to-gel phase transition. We have also studied the counteractive role of TMAO on the effect of urea on lipid membrane when both the osmolytes are present. TMAO draws some urea molecules out of the membrane and thereby reduces the effect of urea on the lipid membrane at lower hydration levels. This is similar to the counteraction of urea's deleterious effects on protein by TMAO. All these observations are consistent with the experimental results and thus provide deep molecular insights into the role of these osmolytes in protecting the fluid phase of the membrane, the key survival strategy against osmotic-stress-induced dehydration.

Effect of TMAO on the Structure and Phase Transition of Lipid Membranes: Potential Role of TMAO in Stabilizing Cell Membranes under Osmotic Stress

Extremophiles adopt strategies to deal with different environmental stresses, some of which are severely damaging to their cell membrane. To combat high osmotic stress, deep-sea organisms synthesize osmolytes, small polar organic molecules, like trimethylamine-N-oxide (TMAO), and incorporate them in the cell. TMAO is known to protect cells from high osmotic or hydrostatic pressure. Several experimental and simulation studies have revealed the roles of such osmolytes on stabilizing proteins. In contrast, the effect of osmolytes on the lipid membrane is poorly understood and broadly debated. A recent experiment has found strong evidence of the possible role of TMAO in stabilizing lipid membranes. Using the molecular dynamics (MD) simulation technique, we have demonstrated the effect of TMAO on two saturated fully hydrated lipid membranes in their fluid and gel phases. We have captured the impact of TMAO's concentration on the membrane's structural properties along with the fluid/gel phase transition temperatures. On increasing the concentration of TMAO, we see a substantial increase in the packing density of the membrane (estimated by area, thickness, and volume) and enhancement in the orientational order of lipid molecules. Having repulsive interaction with the lipid head group, the TMAO molecules are expelled away from the membrane surface, which induces dehydration of the lipid head groups, increasing the packing density. The addition of TMAO also increases the fluid/gel phase transition temperature of the membrane. All of these results are in close agreement with the experimental observations. This study, therefore, provides a molecular-level understanding of how TMAO can influence the cell membrane of deep-sea organisms and help in combating the osmotic stress condition.

A New Lipid Force Field (FUJI)

To explore inhomogeneous and anisotropic systems such as lipid bilayers, the Lennard-Jones particle mesh Ewald (LJ-PME) method has been applied without a conventional isotropic dispersion correction. As the popular AMBER and CHARMM lipid force fields were developed using a cutoff scheme, their lipid bilayers unacceptably shrink when using the LJ-PME method. In this study, a new all-atom lipid force field (FUJI) was developed on the basis of the AMBER force-field scheme including the Lipid14 van der Waals parameters. Point charges were calculated using the restrained electrostatic potentials of many lipid conformers. Further, torsion energy profiles were calculated using high-level ab initio molecular orbitals (LCCSD(T)/Aug-cc-pVTZ//LMP2/Aug-cc-pVTZ), following which the molecular mechanical dihedral parameters were derived through a fast Fourier transform. By incorporation of these parameters into a new lipid force field without fitting experimental data, the desired lipid characteristics such as the area per lipid and lateral diffusion coefficients were obtained through GROMACS molecular dynamics simulations using the LJ-PME method and virtual hydrogen sites. The calculated area per lipid and lateral diffusion coefficients showed satisfactory agreement with experimental data. Furthermore, the electron-density profiles along the membrane normal were calculated for pure lipid bilayers, and the resulting membrane thicknesses agreed well with the experimental values. As the new lipid force field is compatible with FUJI for protein and small molecules, the new FUJI force field will offer accurate modeling for complex systems consisting of various membrane proteins and lipids.

Translocation of a hydroxyl functionalized carbon dot across a lipid bilayer: an all-atom molecular dynamics simulation study

Passive permeation of CD across lipid bilayer is almost impossible. Forced permeation results membrane rupture.