A coarse-grained (CG) model for peptides and proteins was developed as an extension of the Surface Property fItting Coarse grAined (SPICA) force field (FF). The model was designed to examine membrane proteins that are fully compatible with the lipid membranes of the SPICA FF. A preliminary version of this protein model was created using thermodynamic properties, including the surface tension and density in the SPICA (formerly called SDK) FF. In this study, we improved the CG protein model to facilitate molecular dynamics (MD) simulations with a reproduction of multiple properties from both experiments and all-atom (AA) simulations. An elastic network model was adopted to maintain the secondary structure within a single chain. The side-chain analogues reproduced the transfer free energy profiles across the lipid membrane and demonstrated reasonable association free energy (potential of mean force) in water compared to those from AA MD. A series of peptides/proteins adsorbed onto or penetrated into the membrane simulated by the CG MD correctly predicted the penetration depths and tilt angles of peripheral and transmembrane peptides/proteins as comparable to those in the orientations of proteins in membranes (OPM) database. In addition, the dimerization free energies of several transmembrane helices within a lipid bilayer were comparable to those from experimental estimation. Application studies on a series of membrane protein assemblies, scramblases, and poliovirus capsids demonstrated the good performance of the SPICA FF.

ABSTRACTA coarse-grained (CG) model for peptides and proteins was developed as an extension of the SPICA (Surface Property fItting Coarse grAined) force field (FF). The model was designed to examine membrane proteins that are fully compatible with the lipid membranes of the SPICA FF. A preliminary version of this protein model was created using thermodynamic properties, including the surface tension and density in the SPICA (formerly called SDK) FF. In this study, we improved the CG protein model to facilitate molecular dynamics (MD) simulation with a reproduction of multiple properties from both experiments and all-atom (AA) simulations. The side chain analogs reproduced the transfer free energy profiles across the lipid membrane and demonstrated reasonable dimerization free energies in water compared to those from AA-MD. A series of peptides/proteins adsorbed or penetrated into the membrane simulated by the CG-MD correctly predicted the penetration depths and tilt angles of peripheral and transmembrane peptides/proteins comparable to those in the orientation of protein in membrane (OPM) database. In addition, the dimerization free energies of several transmembrane helices within a lipid bilayer were comparable to those from experimental estimation. Application studies on a series of membrane protein assemblies, scramblases, and poliovirus capsids demonstrated a good performance of the SPICA FF.

The lipid matrix in the outer layer of mammalian skin, the stratum corneum, has been previously investigated by multiple biophysical techniques aimed at identifying hydrophilic and lipophilic pathways of permeation. Although consensus is developing over the microscopic structure of the lipid matrix, no molecular-resolution model describes the permeability of all chemical species simultaneously. Using molecular dynamics simulations of a model mixture of skin lipids, the self-assembly of the lipid matrix lamellae has been studied. At higher humidity, the resulting lamellar phase is maintained by partitioning excess water into isolated droplets of controlled size and spatial distribution. The droplets may fuse together to form intralamellar water channels, thereby providing a pathway for the permeation of hydrophilic species. These results reconcile competing data on the outer skin's structure and broaden the scope of molecular-based methods to improve the safety of topical products and to advance transdermal drug delivery.

New material development and commercial application is often quite complex due to the material properties and multiple transformations materials undergo in the supply chain, manufacturing process, and distribution of the finished product. In the fast-moving consumer goods industry of personal and household care products, these complexities are particularly acute due to the focus on and use of "commodity" materials that, at times, have significant variability in material properties. These materials are often formulated into complex liquids or assembled products, which undergo multiple transformations during making and can further undergo additional changes during distribution and use by the consumer (some desired, some not). At each stage of development, manufacturing, and distribution, materials models can be tremendously helpful in material and process selection and optimization. This chapter provides an overview of the current state-of-the-art in materials modeling as applied to the soft materials typically used in household and personal care products, with particular focus on modeling tools that span the length and time scales most relevant for modeling. We review the tools and methods in materials modeling and provide several examples where these tools have been used to guide the development of new materials. We conclude with commentary on additional advancements needed to drive practical application of these modeling tools more broadly for material development.

The architecture of a biological membrane hinges upon the fundamental fact that its properties are determined by more than the sum of its individual components. Studies on model membranes have shown the need to characterize in molecular detail how properties such as thickness, fluidity, and macroscopic bending rigidity are regulated by the interactions between individual molecules in a non-trivial fashion. Simulation-based approaches are invaluable to this purpose but are typically limited to short sampling times and model systems that are often smaller than the required properties. To alleviate both limitations, the use of coarse-grained (CG) models is nowadays an established computational strategy. We here present a new CG force field for cholesterol, which was developed by using measured properties of small molecules, and can be used in combination with our previously developed force field for phospholipids. The new model performs with precision comparable to atomistic force fields in predicting the properties of cholesterol-rich phospholipid bilayers, including area per lipid, bilayer thickness, tail order parameter, increase in bending rigidity, and propensity to form liquid-ordered domains in ternary mixtures. We suggest the use of this model to quantify the impact of cholesterol on macroscopic properties and on microscopic phenomena involving localization and trafficking of lipids and proteins on cellular membranes.

The first line of defense of the human body against harmful agents such as toxic chemicals, viruses and bacteria is the stratum corneum (SC), the skin's outer layer. A simple "brick and mortar" analogy is often used to describe how the lipid matrix (the "mortar") of the stratum corneum holds together keratin-rich and largely impenetrable corneocytes (the "bricks"). Primarily due to uncertainty in the molecular structure of the lipid matrix, the pathways of permeation of toxic chemicals or drugs administered transdermally are not completely understood. We've used molecular dynamics to simulate a multi-lamellae 30-nm cross-section of the lipid matrix, totaling nearly 1 million atoms, under varying conditions of acidity and salt concentration for more than two microseconds. During our simulations, the skin's lipids form spontaneously a multilamellar structure where small pockets of water are separated by regions of partial fusion between the lamellae. The simulations suggest a potential assembly mechanism of the SC having a delicate balance between the lipid, water and salinity contents, which we've interrogated further using our recently developed coarse grain parameters for ceramides and cholesterol. In addition to visualizing a key event in the assembly of the largest human organ, our results further clarify the criteria to develop new delivery vehicles for pharmaceutical treatment.

Ceramides are indispensable constituents of the stratum corneum (SC), the uppermost impermeable layer of human skin. Ceramides with shorter (four- to eight-carbon acyl chains) fatty acid chains increase skin and model membrane permeability, while further shortening of the chain leads to increased resistance to penetration almost as good as that of ceramides from healthy skin (24 carbons long on average). Here we address the extent to which the atomistic CHARMM36 and coarse-grain MARTINI molecular dynamics (MD) simulations reflect the skin permeability data. As a result, we observed the same bell-shaped permeability trend for water that was observed in the skin and multilayer membrane experiments for model compounds. We showed that the enhanced permeability of the short ceramides is mainly caused by the disturbance of their headgroup conformation because of their inability to accommodate the shorter lipid acyl chain into a typical hairpin conformation, which further led to their destabilization and phase separation. As MD simulations described well delicate structural features of SC membranes, they seem to be suitable for further studies of the SC superstructure, including the development of skin penetration enhancers for transdermal drug delivery and skin toxicity risk assessment studies.

, 1009 (2010); 328 Science et al. Virgil Percec and Other Complex Architectures Self-Assembly of Janus Dendrimers into Uniform Dendrimersomes This copy is for your personal, non-commercial use only. clicking here. colleagues, clients, or customers by , you can order high-quality copies for your If you wish to distribute this article to others here. following the guidelines can be obtained by Permission to republish or repurpose articles or portions of articles ): August 21, 2014 www.sciencemag.org (this information is current as of The following resources related to this article are available online at http://www.sciencemag.org/content/328/5981/1009.full.html version of this article at: including high-resolution figures, can be found in the online Updated information and services, http://www.sciencemag.org/content/suppl/2010/05/18/328.5981.1009.DC1.html can be found at: Supporting Online Material http://www.sciencemag.org/content/328/5981/1009.full.html#ref-list-1 , 8 of which can be accessed free: cites 34 articles This article http://www.sciencemag.org/content/328/5981/1009.full.html#related-urls 5 articles hosted by HighWire Press; see: cited by This article has been http://www.sciencemag.org/cgi/collection/chemistry Chemistry subject collections: This article appears in the following

The partitioning behavior of drug-like molecules into biomembranes has a crucial impact on the design and efficacy of therapeutic drugs. Thermodynamic properties connected with the interaction of molecules with membranes can be evaluated by calculating free-energy profiles normal to the membrane surface. We calculated the free-energy profiles of 25 drug-like molecules in a 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) membrane and free energies of solvation in water and heptane using two methods, molecular dynamics (MD) simulations with the Berger lipid force field and COSMOmic, based on a continuum conductor-like screening model for realistic solvation (COSMO-RS). The biased MD simulations (in total ∼22 μs long) were relatively computationally expensive, whereas the COSMOmic approach offered a significantly less expensive alternative. Both methods provided similar results and showed that the studied amphiphilic drug-like molecules accumulate in the membrane, with the majority localized below the head group region. The MD simulations were more lipophilic and gave free-energy profiles that were systematically deeper than those calculated by COSMOmic. To investigate the physical nature of the increased lipophilicity, we analyzed a water/heptane system and identified that it is most likely caused by overestimation of the attractive term of the Lennard-Jones potential in lipid tails. We concluded that COSMOmic can be successfully used for high-throughput computations of global thermodynamic properties, for example, partition coefficients and energy barrier heights, in phosphocholine membranes. In contrast, MD is better for investigating local properties like molecular positioning and orientation in the membrane because they more accurately reflect the complex structure of lipid bilayers. MD is also useful for studies of highly complex systems, for example, drug-membrane-protein interactions.

Since computing resources have advanced enough to allow routine molecular simulation studies of drug molecules interacting with biologically relevant membranes, a considerable amount of work has been carried out with fluid phospholipid systems. However, there is very little work in the literature on drug interactions with gel phase lipids. This poses a significant limitation for understanding permeation through the stratum corneum where the primary pathway is expected to be through a highly ordered lipid matrix. To address this point, we analyzed the interactions of p-aminobenzoic acid (PABA) and its ethyl (benzocaine) and butyl (butamben) esters with two membrane bilayers, which differ in their fluidity at ambient conditions. We considered a dioleoylphosphatidylcholine (DOPC) bilayer in a fluid state and a ceramide 2 (CER2, ceramide NS) bilayer in a gel phase. We carried out unbiased (100 ns long) and biased z-constraint molecular dynamics simulations and calculated the free energy profiles of all molecules along the bilayer normal. The free energy profiles converged significantly slower for the gel phase. While the compounds have comparable affinities for both membranes, they exhibit penetration barriers almost 3 times higher in the gel phase CER2 bilayer. This elevated barrier and slower diffusion in the CER2 bilayer, which are caused by the high ordering of CER2 lipid chains, explain the low permeability of the gel phase membranes. We also compared the free energy profiles from MD simulations with those obtained from COSMOmic. This method provided the same trends in behavior for the guest molecules in both bilayers; however, the penetration barriers calculated by COSMOmic did not differ between membranes. In conclusion, we show how membrane fluid properties affect the interaction of drug-like molecules with membranes.

Studies of drug-membrane interactions witness an ever-growing interest, as penetration, accumulation, and positioning of drugs play a crucial role in drug delivery and metabolism in human body. Molecular dynamics simulations complement nicely experimental measurements and provide us with new insight into drug-membrane interactions; however, the quality of the theoretical data dramatically depends on the quality of the force field used. We calculated the free energy profiles of 11 molecules through a model dimyristoylphosphatidylcholine (DMPC) membrane bilayer using five force fields, namely Berger, Slipids, CHARMM36, GAFFlipids, and GROMOS 43A1-S3. For the sake of comparison, we also employed the semicontinuous tool COSMOmic. High correlation was observed between theoretical and experimental partition coefficients (log K). Partition coefficients calculated by all-atomic force fields (Slipids, CHARMM36, and GAFFlipids) and COSMOmic differed by less than 0.75 log units from the experiment and Slipids emerged as the best performing force field. This work provides the following recommendations (i) for a global, systematic and high throughput thermodynamic evaluations (e.g., log K) of drugs COSMOmic is a tool of choice due to low computational costs; (ii) for studies of the hydrophilic molecules CHARMM36 should be considered; and (iii) for studies of more complex systems, taking into account all pros and cons, Slipids is the force field of choice.

The level of hydration controls the cohesion between apposed lamellae of saturated free fatty acids found in the lipid matrix of stratum corneum, the outermost layer of mammalian skin. This multilamellar lipid matrix is highly impermeable to water and ions, so that the local hydration shell of its fatty acids may not always be in equilibrium with the acidity and relative humidity, which significantly change over a course of days during skin growth. The homeostasis of the stratum corneum at each moment of its growth likely requires a balance between two factors, which affect in opposite ways the diffusion of hydrophilic species through the stratum corneum: (i) an increase in water order as the lipid lamellae come in closer contact, and (ii) a decrease in water order as the fraction of charged fatty acids is lowered by pH. Herein molecular dynamics simulations are employed to estimate the impact of both effects on water molecules confined between lamellae of fatty acids. Under conditions where membrane undulations are energetically favorable, the charged fatty acids are able to sequester cations around points of contact between lamellae that are fully dehydrated, while essentially maintaining a multilamellar structure for the entire system. This observation suggests that the undulations of the fatty acid lamellae control the diffusion of hydrophilic species through the water phase by altering the positional and rotational order of water molecules in the embedded/occluded “droplets.”

The effect of carbon-based nanoparticles (CNPs) on biological systems is currently of great concern. Yet, few experimental techniques are capable of directly imaging and probing the energetics of such nano-bio systems. Here, we use coarse grain molecular dynamics simulations to study spherical fullerene molecules interacting with dipalmitoyl phosphatidylcholine (DPPC) lipid membranes. Using free energy calculations we show that all the tested fullerene molecules can spontaneously diffuse into both a lipid bilayer and a lipid monolayer. In addition, we establish that large fullerene molecules tend to partition preferentially into bilayers, which affects the lipid monolayer-to-bilayer transition during the respiration cycle. Our results identify a possible CNP perturbation to the function of the pulmonary monolayer membrane and suggest a potential pathway for CNP entry into the body through lung inhalation.

Coarse-grained (CG) molecular models are now widely used to understand the structure and functionality of macromolecular self-assembling systems. In the last few years, significant efforts have been devoted to construct quantitative CG models based on data from molecular dynamics (MD) simulations with more detailed all-atom (AA) intermolecular force fields as well as experimental thermodynamic data. We review some of the recent progress pertaining to the MD simulation of self-assembling macromolecular systems, using as illustrations the application of CG models to probe surfactant and lipid self-assembly including liposome and dendrimersome formation as well as the interaction of biomembranes with nanoparticles.

Functionalized single-walled carbon nanotubes (SWNTs) are widely applied in biomedical science. To understand the interaction between SWNTs and biological systems, various studies have attempted to use coarse-grained molecular dynamics (CGMD). However, there is limited validation of the existing CG models of SWNTs. Here, we present CG models for both pristine and carboxylated SWNTs which are validated against experimental dispersion data. In addition, we present the first ever DLVO analysis of the colloidal stability of parallel SWNTs and establish that the solvent-induced repulsion between fullerenes, which is not considered in DLVO theory, is crucial to obtain a correct physical picture of SWNT dispersibility. The results presented here provide physical insight into the colloidal stability of SWNTs and can be applied to large-scale MD studies of biological systems.

The computational design of advanced materials based on surfactant self-assembly without ever stepping foot in the laboratory is an important goal, but there are significant barriers to this approach, because of the limited spatial and temporal scales accessible by computer simulations. In this paper, we report our work to bridge the gap between laboratory and computational time scales by implementing the coarse-grained (CG) force field previously reported by Shinoda et al. [Shinoda, W.; DeVane, R.; Klein, M. L. Mol. Simul. 2007, 33, 27-36] into the HOOMD-Blue graphical processing unit (GPU)-accelerated molecular dynamics (MD) software package previously reported by Anderson et al. [Anderson, J. A.; Lorenz, C. D.; Travesset, A. J. Comput. Phys. 2008, 227, 5342-5359]. For a system of 25 750 particles, this implementation provides performance on a single GPU, which is superior to that of a widely used parallel MD simulation code running on an optimally sized CPU-based cluster. Using our GPU setup, we have collected 0.6 ms of MD trajectory data for aqueous solutions of 7 different nonionic polyethylene glycol (PEG) surfactants, with most of the systems studied representing ∼1 000 000 atoms. From this data, we calculated various properties as a function of the length of the hydrophobic tails and PEG head groups. Specifically, we determined critical micelle concentrations (CMCs), which are in good agreement with experimental data, and characterized the size and shape of micelles. However, even with the microsecond trajectories employed in this study, we observed that the micelles composed of relatively hydrophobic surfactants are continuing to grow at the end of our simulations. This suggests that the final micelle size distributions of these systems are strongly dependent on initial conditions and that either longer simulations or advanced sampling techniques are needed to properly sample their equilibrium distributions. Nonetheless, the combination of coarse-grained modeling and GPU acceleration marks a significant step toward the computational prediction of the thermodynamic properties of slowly evolving surfactant systems.

A coarse-grained (CG) model has been developed for the anionic surfactant class, linear alkylbenzene sulfonates (LAS), which are the most widely used synthetic surfactants. The development work started from a systematic examination of tens of CG water models with different resolutions, interaction potentials (Lennard-Jones and Morse), and cut-off distances. The relationships between the parameters under specific choices of the above options and the thermodynamic properties, such as density, surface tension, and compressibility, were found to fit simple mathematical equations. The limits of applicability of these CG water models were explored by checking the melting temperature. Considering both efficiency and accuracy, a CG water model which includes three water molecules in one CG site was chosen. Correspondingly, the LAS molecules were mapped into CG sites each contains approximately three heavy atoms and connected hydrogens. Structural data obtained from atomistic simulations and thermodynamic data from experiments were used as targets to parameterize standard potential forms for bonded and non-bonded interactions. An extensive evaluation of the CG model for a series of different alkane molecules (aliphatic or aromatic, linear or branched) shows that the present model is not only reliable, but also transferable. This point is crucial to assure that the model is capable of representing different isomers and homologues in the LAS family. The resulting model is easily implemented into standard MD codes. The added computational efficiency permits the simulation of the self-assembly of LAS solutions starting from a random configuration. The model is shown to accurately reproduce the phase behavior of solutions of pure isomers of sodium dodecylbenzene sulfonate, despite the fact that phase behavior was not directly taken into account in the parameterization.

Many unusual phenomena in surfactant solutions below the critical micelle concentration have been attributed to the existence of premicelles. Here we investigate the properties of premicelles in aqueous solutions of non-ionic surfactants surrounding the critical micelle concentration (CMC) using advanced computer simulations. The concentration of premicelles increases with total concentration up to the CMC, above which the concentrations of free monomers and premicelles remain nearly constant. Analysis shows premicellar structures ranging from spherical to chain-like. Monomer exchange kinetics are observed to depend strongly on the aggregate size, with premicelles exhibiting extremely short (1.3ns) relaxation times relative to full sized micelles (72ns).