Formulation and Skin Permeation of Active-Loaded Lipid Nanoparticles: Evaluation and Screening by Synergizing Molecular Dynamics Simulations and Experiments.

Encapsulation into nanoparticles (NPs) is a potential method to deliver pharmaceutical/cosmetic actives deep into the skin. However, understanding the NP formulations and underlying mechanism of active delivery to skin has scarcely been studied. We report a simulation platform that screens, evaluates, formulates, and provides atomic-resolution interpretation of NP-based formulations, and reveals the active permeation mechanism from NPs to skin. First, three actives, namely, ferulic acid (FA), clotrimazole (CZE), and tretinoin (TTN), and five lipid excipients' (Compritol, Precirol, Geleol, Gelot, Gelucire) combinations were screened by MD simulations for the best pairs. For each suggested pair, the actual active and lipid compositions for the synthesis of stable NP formulations were then obtained by experiments. MD simulations demonstrate that in NP formulations, FA and CZE actives are present at the surface of the NPs, whereas TTN actives are present at both the surface and interior of the NP core. The NP shapes obtained by simulation perfectly match with experiments. For each NP, separate MD simulations illustrate that active-loaded NPs approach the skin surface quickly, and then actives translocate from NP surface to skin surface followed by penetration of NPs through skin. The driving force for the translocation which initiates during the penetration process, is the stronger active-skin interaction compared to active-NP interaction. Permeation free energy indicates spontaneous transfer of actives from solution phase to the surface of the skin bilayer. The free energy barriers are increased in the order of FA < TTN < CZE. Significantly lower diffusions of actives are obtained in the main barrier region compared to bulk, and the average diffusion coefficients of actives are in the same order of magnitude (∼10-6 cm2/s). The estimated permeability coefficients (log P) of actives are mainly governed by free energy barriers. The study would facilitate the development of novel lipid-based NP formulations for personal-care/pharmaceutical applications.

Catalytic amino acid production from biomass-derived intermediates.

Amino acids are the building blocks for protein biosynthesis and find use in myriad industrial applications including in food for humans, in animal feed, and as precursors for bio-based plastics, among others. However, the development of efficient chemical methods to convert abundant and renewable feedstocks into amino acids has been largely unsuccessful to date. To that end, here we report a heterogeneous catalyst that directly transforms lignocellulosic biomass-derived α-hydroxyl acids into α-amino acids, including alanine, leucine, valine, aspartic acid, and phenylalanine in high yields. The reaction follows a dehydrogenation-reductive amination pathway, with dehydrogenation as the rate-determining step. Ruthenium nanoparticles supported on carbon nanotubes (Ru/CNT) exhibit exceptional efficiency compared with catalysts based on other metals, due to the unique, reversible enhancement effect of NH3 on Ru in dehydrogenation. Based on the catalytic system, a two-step chemical process was designed to convert glucose into alanine in 43% yield, comparable with the well-established microbial cultivation process, and therefore, the present strategy enables a route for the production of amino acids from renewable feedstocks. Moreover, a conceptual process design employing membrane distillation to facilitate product purification is proposed and validated. Overall, this study offers a rapid and potentially more efficient chemical method to produce amino acids from woody biomass components.

Water Desalination through Zeolitic Imidazolate Framework Membranes: Significant Role of Functional Groups.

A molecular simulation study is reported for water desalination through five zeolitic imidazolate framework (ZIF) membranes, namely ZIF-25, -71, -93, -96, and -97. The five ZIFs possess identical rho-topology but differ in functional groups. The rejection of salt (NaCl) is found to be around 97% in ZIF-25, and 100% in the other four ZIFs. The permeance ranges from 27 to 710 kg/(m(2)·h·bar), about one∼two orders of magnitude higher compared with commercial reverse osmosis membranes. Due to a larger aperture size da, ZIF-25, -71, and -96 exhibit a much higher water flux than ZIF-93 and -97; however, the flux in ZIF-25, -71, and -96 is governed by the polarity of functional group rather than da. With the hydrophobic CH3 group, ZIF-25 has the highest flux despite the smallest da among ZIF-25, -71, and -96. The lifetime of hydrogen bonding in ZIF-25 is shorter than that in ZIF-71 and -96. Furthermore, water molecules undergo a fast flushing motion in ZIF-25, but frequent jumping in ZIF-96 and particularly in ZIF-97. An Arrhenius-type relationship is found between water flux in ZIF-25 and temperature, and the activation energy is predicted to be 6.5 kJ/mol. This simulation study provides a microscopic insight into water desalination in a series of ZIFs, reveals the key factors (aperture size and polarity of functional group) governing water flux, and suggests that ZIF-25 might be an interesting reverse osmosis membrane for high-performance water desalination.

Evaluation of polymer-preservative interactions for preservation efficacy: molecular dynamics simulation and QSAR approaches.

Preservatives are critical ingredients in various pharmaceutical and consumer products. In particular, a high efficacy preservative system is essential in enhancing the shelf-life and safety of these products. However, the development of such a preservative system heavily relies on experimental approaches. In this study, molecular dynamics (MD) simulation was complemented with quantitative structure-activity relationship (QSAR) modelling to comprehensively evaluate polymer-preservative interactions between three different polymers (polyethylene terephthalate, PET; polypropylene, PP; and cellulose) and a series of preservatives from the classes of aliphatic, aromatic, and organic acids. First, adsorption of preservatives onto polymer surfaces was simulated in an aqueous environment. The preservatives did not adhere to hydrophilic cellulose, but most preservatives were adsorbed by PET and PP in distinct configurations. Interaction energies (IEs) between the preservatives and the polymers generally increase from cellulose to PP and PET. The diffusion coefficients of preservatives are dependent on polymer nature, preservative structure, and their resulting molecular interactions. Linear and low molecular weight preservatives exhibit higher diffusion coefficients in polymers. For a particular preservative, diffusion coefficients increased in the order of cellulose < PET < PP. Finally, using MD properties and molecular descriptors of preservatives, QSAR models were developed to identify key descriptors of preservatives and predict their IEs and diffusion coefficients in polymers. This study demonstrates a computational approach for identifying critical materials properties, and predicting polymer-preservative molecular interactions in water. Such an approach streamlines the rational selection and design of high efficacy preservative systems for various pharmaceutical, food and cosmetic products. Furthermore, the integrated computational strategy also reduces trial-and-error experimental efforts, thereby accelerating the development of high efficacy preservative systems.

Molecular Interactions between APIs and Enteric Polymeric Excipients in Solid Dispersion: Insights from Molecular Simulations and Experiments.

Solid dispersion of poorly soluble APIs is known to be a promising strategy to improve dissolution and oral bioavailability. To facilitate the development and commercialization of a successful solid dispersion formulation, understanding of intermolecular interactions between APIs and polymeric carriers is essential. In this work, first, we assessed the molecular interactions between various delayed-release APIs and polymeric excipients using molecular dynamics (MD) simulations, and then we formulated API solid dispersions using a hot melt extrusion (HME) technique. To assess the potential API-polymer pairs, three quantities were evaluated: (a) interaction energy between API and polymer [electrostatic (Ecoul), Lenard-Jones (ELJ), and total (Etotal)], (b) energy ratio (API-polymer/API-API), and (c) hydrogen bonding between API and polymer. The Etotal quantities corresponding to the best pairs: NPX-Eudragit L100, NaDLO-HPMC(P), DMF-HPMC(AS) and OPZ-HPMC(AS) were -143.38, -348.04, -110.42, and -269.43 kJ/mol, respectively. Using a HME experimental technique, few API-polymer pairs were successfully extruded. These extruded solid forms did not release APIs in a simulated gastric fluid (SGF) pH 1.2 environment but released them in a simulated intestinal fluid (SIF) pH 6.8 environment. The study demonstrates the compatibility between APIs and excipients, and finally suggests a potential polymeric excipient for each delayed-release API, which could facilitate the development of the solid dispersion of poorly soluble APIs for dissolution and bioavailability enhancement.

Metal-organic framework supported ionic liquid membranes for CO2 capture: anion effects.

IRMOF-1 supported ionic liquid (IL) membranes are investigated for CO(2) capture by atomistic simulation. The ILs consist of identical cation 1-n-butyl-3-methylimidazolium [BMIM](+), but four different anions, namely hexafluorophosphate [PF(6)](-), tetrafluoroborate [BF(4)](-), bis(trifluoromethylsulfonyl)imide [Tf(2)N](-), and thiocyanate [SCN](-). As compared with the cation, the anion has a stronger interaction with IRMOF-1 and a more ordered structure in IRMOF-1. The small anions [PF(6)](-), [BF(4)](-), and [SCN](-) prefer to locate near to the metal-cluster, particularly the quasi-spherical [PF(6)](-) and [BF(4)](-). In contrast, the bulky and chain-like [BMIM](+) and [Tf(2)N](-) reside near the phenyl ring. Among the four anions, [Tf(2)N](-) has the weakest interaction with IRMOF-1 and thus the strongest interaction with [BMIM](+). With increasing the weight ratio of IL to IRMOF-1 (W(IL/IRMOF-1)), the selectivity of CO(2)/N(2) at infinite dilution is enhanced. At a given W(IL/IRMOF-1), the selectivity increases as [Tf(2)N](-) < [PF(6)](-) < [BF(4)](-) < [SCN](-). This hierarchy is predicted by the COSMO-RS method, and largely follows the order of binding energy between CO(2) and anion estimated by ab initio calculation. In the [BMIM][SCN]/IRMOF-1 membrane with W(IL/IRMOF-1) = 1, [SCN](-) is identified to be the most favorable site for CO(2) adsorption. [BMIM][SCN]/IRMOF-1 outperforms polymer membranes and polymer-supported ILs in CO(2) permeability, and its performance surpasses Robeson's upper bound. This simulation study reveals that the anion has strong effects on the microscopic properties of ILs and suggests that MOF-supported ILs are potentially intriguing for CO(2) capture.

Molecular dynamics simulations to elucidate translocation and permeation of active from lipid nanoparticle to skin: complemented by experiments.

One of the most realistic approaches for delivering actives (pharmaceuticals/cosmetics) deep into skin layers is encapsulation into nanoparticles (NPs). Nonetheless, molecular-level mechanisms related to active delivery from NPs to the skin have scarcely been studied despite the large number of synthesis and characterization studies. We herein report the underlying mechanism of active translocation and permeation through the outermost layer of skin, the stratum corneum (SC), via molecular dynamics (MD) simulations complemented by experimental studies. A SC molecular model is constructed using current state-of-the-art methodology via incorporating the three most abundant skin lipids: ceramides, free fatty acids, and cholesterol. As a potent antioxidant, ferulic acid (FA) is used as the model active, and it is loaded into Gelucire 50/13 NP. MD simulations elucidate that, first, FA-loaded NP approaches the skin surface quickly, followed by slight penetration and adsorption onto the upper skin surface; FA then translocates from the NP surface to the skin surface due to stronger NP-skin interactions compared to the FA-NP interactions; then, once FA is released onto the skin surface, it slowly permeates deep into the skin bilayer. Both the free energy and resistance to permeation not only indicate the spontaneous transfer of FA from the bulk to the skin surface, but they also reveal that the main barrier against permeation exists in the middle of the lipid hydrophobic tails. Significantly lower diffusion of FA is obtained in the main barrier region compared to the bulk. The estimated permeability coefficient (log P) values are found to be higher than the experimental values. Importantly, the permeation process evaluated via MD simulations perfectly matches with experiments. The study suggests a molecular simulation platform that provides various crucial insights relating to active delivery from loaded NP to skin, and it could facilitate the design and development of novel NP-based formulations for transdermal delivery and the topical application of drugs/cosmetics.

Zeolitic Imidazolate Framework Membranes for Organic Solvent Nanofiltration: A Molecular Simulation Exploration.

Organic solvents are intensively used in chemical and pharmaceutical industries. Their separation and recovery account for a significant portion of energy consumption and capital cost in many industrial processes. In this study, three microporous crystalline zeolitic imidazolate frameworks (ZIF-25, ZIF-71, and ZIF-96) are investigated as organic solvent nanofiltration (OSN) membranes by molecular simulations. The fluxes of five solvents (methanol, ethanol, acetone, acetonitrile, and n-hexane) are predicted. Despite the smallest aperture size among the three ZIFs, ZIF-25 exhibits the highest flux for polar solvents (methanol, ethanol, acetone, and acetonitrile) because of its hydrophobic nature, whereas hydrophilic ZIF-96 shows the highest flux for nonpolar n-hexane. The analysis of structural information and interaction energy reveals that the solvent-framework interaction is crucial to determine solvent permeation. Good correlations between solvent permeances and a combination of solvent properties are found. In the presence of a model solute (paracetamol), solvent permeances are marginally affected; moreover, the rejection of paracetamol is 100% for the three ZIF membranes in all five solvents. This study highlights that the pore chemistry, in addition to pore size, plays an important role in solvent permeation, and it suggests that ZIFs are potential OSN membranes for the recovery of organic solvents.

Reversed thermo-switchable molecular sieving membranes composed of two-dimensional metal-organic nanosheets for gas separation.

It is highly desirable to reduce the membrane thickness in order to maximize the throughput and break the trade-off limitation for membrane-based gas separation. Two-dimensional membranes composed of atomic-thick graphene or graphene oxide nanosheets have gas transport pathways that are at least three orders of magnitude higher than the membrane thickness, leading to reduced gas permeation flux and impaired separation throughput. Here we present nm-thick molecular sieving membranes composed of porous two-dimensional metal-organic nanosheets. These membranes possess pore openings parallel to gas concentration gradient allowing high gas permeation flux and high selectivity, which are proven by both experiment and molecular dynamics simulation. Furthermore, the gas transport pathways of these membranes exhibit a reversed thermo-switchable feature, which is attributed to the molecular flexibility of the building metal-organic nanosheets.

Seawater Pervaporation through Zeolitic Imidazolate Framework Membranes: Atomistic Simulation Study.

An atomistic simulation study is reported for seawater pervaporation through five zeolitic imidazolate framework (ZIF) membranes including ZIF-8, -93, -95, -97, and -100. Salt rejection in the five ZIFs is predicted to be 100%. With the largest aperture, ZIF-100 possesses the highest water permeability of 5 × 10(-4) kg m/(m(2) h bar), which is substantially higher compared to commercial reverse osmosis membranes, as well as zeolite and graphene oxide pervaporation membranes. In ZIF-8, -93, -95, and -97 with similar aperture size, water flux is governed by framework hydrophobicity/hydrophilicity; in hydrophobic ZIF-8 and -95, water flux is higher than in hydrophilic ZIF-93 and -97. Furthermore, water molecules in ZIF-93 move slowly and remain in the membrane for a long time but undergo to-and-fro motion in ZIF-100. The lifetime of hydrogen bonds in ZIF-93 is found to be longer than in ZIF-100. This simulation study quantitatively elucidates the dynamic and structural properties of water in ZIF membranes, identifies the key governing factors (aperture size and framework hydrophobicity/hydrophilicity), and suggests that ZIF-100 is an intriguing membrane for seawater pervaporation.

Glucose recovery from aqueous solutions by adsorption in metal-organic framework MIL-101: a molecular simulation study.

A molecular simulation study is reported on glucose recovery from aqueous solutions by adsorption in metal-organic framework MIL-101. The F atom of MIL-101 is identified to be the most favorable adsorption site. Among three MIL-101-X (X = H, NH2 or CH3), the parent MIL-101 exhibits the highest adsorption capacity and recovery efficacy. Upon functionalization by -NH2 or -CH3 group, the steric hindrance in MIL-101 increases; consequently, the interactions between glucose and framework become less attractive, thus reducing the capacity and mobility of glucose. The presence of ionic liquid, 1-ethyl-3-methyl-imidazolium acetate, as an impurity reduces the strength of hydrogen-bonding between glucose and MIL-101, and leads to lower capacity and mobility. Upon adding anti-solvent (ethanol or acetone), a similar adverse effect is observed. The simulation study provides useful structural and dynamic properties of glucose in MIL-101, and it suggests that MIL-101 might be a potential candidate for glucose recovery.