The development of the CHARMM lipid force field (FF) can be traced back to the early 1990s with its current version denoted CHARMM36 (C36). The parametrization of C36 utilized high-level quantum ...mechanical data and free energy calculations of model compounds before parameters were manually adjusted to yield agreement with experimental properties of lipid bilayers. While such manual fine-tuning of FF parameters is based on intuition and trial-and-error, automated methods can identify beneficial modifications of the parameters via their sensitivities and thereby guide the optimization process. This work introduces a semi-automated approach to reparametrize the CHARMM lipid FF with consistent inclusion of long-range dispersion through the Lennard-Jones particle-mesh Ewald (LJ-PME) approach. The optimization method is based on thermodynamic reweighting with regularization with respect to the C36 set. Two independent optimizations with different topology restrictions are presented. Targets of the optimizations are primarily liquid crystalline phase properties of lipid bilayers and the compression isotherm of monolayers. Pair correlation functions between water and lipid functional groups in aqueous solution are also included to address headgroup hydration. While the physics of the reweighting strategy itself is well-understood, applying it to heterogeneous, complex anisotropic systems poses additional challenges. These were overcome through careful selection of target properties and reweighting settings allowing for the successful incorporation of the explicit treatment of long-range dispersion, and we denote the newly optimized lipid force field as C36/LJ-PME. The current implementation of the optimization protocol will facilitate the future development of the CHARMM and related lipid force fields.
The diversity of petal and leaf color can improve the ornamental value of rapeseed and promote the development of agriculture and tourism. The two copies of carotenoid isomerase gene (
) in
(
and
) ...was edited using the CRISPR/Cas9 system in the present study. The mutation phenotype of creamy white petals and yellowish leaves could be recovered only in targeted mutants of both
functional copies, indicating that the redundant roles of
and
are vital for the regulation of petal and leaf color. The carotenoid content in the petals and leaves of the
double mutant was significantly reduced. The chalcone content, a vital substance that makes up the yellow color, also decreased significantly in petals. Whereas, the contents of some carotenes (lycopene, α-carotene, γ-carotene) were increased significantly in petals. Further, transcriptome analysis showed that the targeted mutation of
resulted in the significant down-regulation of important genes
and
in the carotenoid and flavonoid synthesis pathways, respectively; however, the expression of other genes related to carotenes and xanthophylls synthesis, such as
,
,
and
, was up-regulated. This indicates that the molecular mechanism regulating petal color variation in
is more complicated than those reported in
and other
species. These results provide insight into the molecular mechanisms underlying flower color variation in rapeseed and provides valuable resources for rapeseed breeding.
Accurate lipid force field (FF) parameters used in molecular dynamics (MD) simulations are crucial for understanding the properties of lipid-containing systems and biological processes related to ...lipids. The last update of the CHARMM36 united atom chain model (C36UA) was in 2013 Lee, S. J. Phys. Chem. B 2014, 118, 547 556; it utilized CHARMM36 (C36) lipid FF parameters for headgroups and OPLS-UA Lennard-Jones (LJ) parameters for tails. Simulations with the FF were able to reproduce many experimental observables of lipid bilayers accurately, but to be more applicable for a wide range of lipids, additional FF parameter optimization was needed. In this work, we present an update of the model, named C36UAr. The parameterization included the LJ parameters for hydrocarbons and related dihedrals. Bulk liquid properties (density, heat of vaporization, isothermal compressibility, and diffusion constant) of model compounds were used as targets for the LJ parameter fitting, and dihedrals were fit to either quantum mechanical (QM) or potential of mean force (PMF) calculations using C36. Thermodynamic reweighting was used to further improve the parameters. Bilayer simulations of various lipid headgroups (phosphatidylcholine, phosphatidylethanolamine, and phosphatidylglycerol) and tails (saturated, monounsaturated, and polyunsaturated) were performed to validate the model, and significant improvements were seen in bilayer properties, including surface area, membrane thicknesses, NMR deuterium order parameters, and density profiles. C36UAr was also compared to the hydrogen mass repartitioning (HMR) method. The high accuracy and competitive efficiency shown in this study make C36UAr one of the best choices for studies of membrane structure and membrane-associated proteins.
The multilayered myelin sheath is a critical component of both central and peripheral nervous systems, forming a protective barrier against axonal damage and facilitating the movement of nervous ...impulses. It is primarily composed of cholesterol (CHL1), phosphatidylcholine (PC), phosphatidylethanolamine (PE), phosphatidylserine (PS), phosphatidylinositol (PI), sphingomyelin (SM), and galactosylceramide (GalCer) lipids. For rat sciatic nerve myelin (part of the peripheral nervous system, PNS), it has been found that cholesterol and unsaturated fatty acid contents are significantly lower in diabetic than in non-diabetic conditions. In this study, lipid compositions from experimental data are used to create four model rat sciatic nerve myelin lipid bilayers: PI-containing (non-diabetic and diabetic) and PS-containing (non-diabetic and diabetic), which were then simulated using the all-atom CHARMM36 force field. Simulation results of diabetic membranes indicate less rigid, more laterally expansive, and thinner bilayers as well as potentially reduced interactions between GalCer on opposing myelin leaflets, supporting a direct role of the cholesterol content decrease in instigating myelin deterioration and diabetic peripheral neuropathy. Compared to PI-lipids, PS-lipids were found to cause higher inter-lipid spacing and decreased order within membranes as a result of their smaller headgroup size and higher inter-lipid hydrogen bonding potential, which allow them to more frequently reside deeper in the membrane plane and produce pushing effects on other lipids. GalCer deuterium order parameters and non-diabetic headgroup-to-headgroup bilayer thicknesses were compared to experimental data, exhibiting close alignment, supporting the future usage of these models to study the PNS myelin sheath.
Arabidopsis thaliana is an important model organism, which has attracted many biologists. While most research efforts have been on studying the genetics and proteins of this organism, a systematic ...study of its lipidomics is lacking. Here, we present a novel, asymmetric model of its cell membrane with its lipid composition consisting of five glycerophospholipids, two glycolipids, and sitosterol determined from multiple independent experiments. A typical lipid type in plant membranes is glycosyl inositol phosphoryl ceramide (GIPC), which accounts for about 10% of the total lipids in the outer leaflet in our model. Two symmetric models representing the inner and outer leaflets of the membrane were built and simulated until equilibrium was reached and then combined to form the asymmetric model. Our results indicate that the outer leaflet is more rigid and tightly packed compared to the inner leaflet. Pressure profiles for the two leaflets are overall similar though the outer leaflet exhibits larger oscillations. A special focus on lipid organization is discussed and the interplay between glycolipids and sitosterols is found to be important. The current model provides a baseline for future modeling of similar membranes and can be used to study partitioning of small molecules in the membrane or further developed to study the interaction between plant membrane proteins and lipids.
Summary
Yellow‐seed is widely accepted as a good‐quality trait in Brassica crops. Previous studies have shown that the flavonoid biosynthesis pathway is essential for the development of seed colour, ...but its function in Brassica napus, an important oil crop, is poorly understood. To systematically explore the gene functions of the flavonoid biosynthesis pathway in rapeseed, several representative TRANSPARENT TESTA (TT) genes, including three structural genes (BnaTT7, BnaTT18, BnaTT10), two regulatory genes (BnaTT1, BnaTT2) and a transporter (BnaTT12), were selected for targeted mutation by CRISPR/Cas9 in the present study. Seed coat colour, lignin content, seed quality and yield‐related traits were investigated in these Bnatt mutants together with Bnatt8 generated previously. These Bnatt mutants produced seeds with an elevated seed oil content and decreased pigment and lignin accumulation in the seed coat without any serious defects in the yield‐related traits. In addition, the fatty acid (FA) composition was also altered to different degrees, i.e., decreased oleic acid and increased linoleic acid and α‐linolenic acid, in all Bnatt mutants except Bnatt18. Furthermore, gene expression analysis revealed that most of BnaTT mutations resulted in the down‐regulation of key genes related to flavonoid and lignin synthesis, and the up‐regulation of key genes related to lipid synthesis and oil body formation, which may contribute to the phenotype. Collectively, our study generated valuable resources for breeding programs, and more importantly demonstrated the functional divergence and overlap of flavonoid biosynthesis pathway genes in seed coat colour, oil content and FA composition of rapeseed.
Cell membranes are composed of a variety of lipids and proteins where they interact with each other to fulfill their roles. The first step in modeling these interactions in molecular simulations is ...to have reliable mimetics of the membrane’s lipid environment. This Feature Article presents our recent efforts to model complex cellular membranes using all-atom force fields. A short review of the CHARMM36 (C36) lipid force field and its recent update to incorporate the long-range dispersion is presented. Key examples of model membranes mimicking various species and organelles are given. These include single-celled organisms such as bacteria (E. coli., chlamydia, and P. aeruginosa) and yeast (plasma membrane, endoplasmic reticulum, and trans-Golgi network) and more advanced ones such as plants (soybean and Arabidopsis thaliana) and mammals (ocular lens, stratum corneum, and peripheral nerve myelin). Leaflet asymmetry in composition has also been applied to some of these models. With the increased lipid diversity in the C36 lipid FF, these complex models can better reflect the structural, mechanical, and dynamic properties of realistic membranes and open an opportunity to study biological processes involving other molecules.