Formation mechanism of bound rubber in elastomer nanocomposites: a molecular dynamics simulation study

Bound rubber plays a key role in the mechanical reinforcement of elastomer nanocomposites. In the present work, we reveal the formation mechanism of bound rubber in elastomer nanocomposites, using the coarse-grained molecular dynamics simulations. For the polymer–nanoparticle system, the “chain bridge” connected with neighboring nanoparticles forms, once the gap between two neighboring nanoparticles is less than the polymer size. The polymer–nanoparticle–solvent systems, mimicking the oil-swollen rubber in the experiment, are simulated with three models. From the analysis of the potential energy, the static structure and dynamic diffusing processes, all the models indicate that the increase of the volume fraction of the nanoparticles and the polymer−nanoparticle interaction strength could promote the formation of the bound rubber. The existence of solvent disrupts the bound rubber, and eventually deteriorates the mechanical properties. These simulations could provide some theoretical guidance for a better understanding of the formation mechanism of the bound rubber, which is helpful for designing the elastomer materials with excellent mechanical properties.

The formation mechanism of the bound rubber in elastomer nanocomposites using the coarse-grained molecular-dynamics simulations.     

Exploring the pore charge dependence of K+ and Cl− permeation across a graphene monolayer: a molecular dynamics study

Selective permeation through graphene nanopores is attracting increasing interest as an efficient and cost-effective technique for water desalination and purification. In this work, using umbrella sampling and molecular dynamics simulations with constant electric field, we analyze the influence of pore charge on potassium and chloride ion permeation. As pore charge is increased, the barrier of the potential of mean force (PMF) gradually decreases until it turns into a well split in two subminima. While in the case of K⁺ this pattern can be explained as an increasing electrostatic compensation of the desolvation cost, in the case of Cl⁻ the pattern can be attributed to the accumulation of a concentration polarization layer of potassium ions screening pore charge. The analysis of potassium PMFs in terms of forces revealed a conflicting influence on permeation of van der Waals and electrostatic forces that both undergo an inversion of their direction as pore charge is increased. Even if the most important transition involves the interplay between the electrostatic forces exerted by graphene and water, the simulations also revealed an important role of the changing distribution of potassium and chloride ions. The influence of pore charge on the orientation of water molecules was also found to affect the van der Waals forces they exert on potassium.

Increase of graphene pore charge determines decrease of PMF barrier that turns into well: current increases, reaches plateau and declines.     

Molecular dynamics study of convective heat transfer mechanism in a nano heat exchanger

With the rapid development of micro/nano electro-mechanical systems, the convective heat transfer at the micro/nanoscale has been widely studied for the thermal management of micro/nano devices. Here we investigate the convective heat transfer mechanism of a nano heat exchanger by the employment of molecular dynamics simulation with a modified thermal pump method. First, the temperature jump and velocity slip are observed at the wall–fluid interfaces of the nano heat exchanger. Moreover, the larger Kapitza resistance in the entrance region weakens the convective heat transfer. Second, the heat transfer performance of the nano heat exchanger can be improved by increasing the surface wettability of the solid walls owing to more fluid atoms being involved in heat transport at the walls when the wall–fluid interaction is enhanced. Meanwhile, the strong surface wettability results in the appearance of the quasi-solid fluid layers, which improves the heat transfer between walls and fluids. Finally, we point out that when the surface wettability of the nano heat exchanger is weak, the heat transfer of the hot fluid side is better than that of the cold fluid side, while the convective heat transfer performances of the cold and hot fluid sides are reversed when the surface wettability is strong. This is because of the feebler temperature jump of the hot fluid side when wall–fluid interaction is small and the greater velocity slip of the cold fluid side for walls with large wall–fluid interaction.

The convective heat transfer mechanism in a nano heat exchanger is investigated using molecular dynamics simulation.     

Terpenes in ethanol: haloperidol permeation and partition through human skin and stratum corneum changes.

Carvacrol, linalool and alpha-terpineol (5% w/v) in 50% ethanol were used to enhance the permeation of haloperidol (HP) through human skin in vitro and their enhancement mechanism was investigated with HP-stratum corneum (SC) binding studies, fourier transform infrared spectroscopy (FT-IR) and differential scanning calorimetry (DSC). Carvacrol followed by terpineol and linalool enhanced flux and permeability coefficient but only carvacrol provided the required plasma concentration and the permeated daily doses. All terpenes increased the activity coefficient of HP in the skin. Carvacrol increased the lag time, which could be due to slow redistribution within SC. The thermogram of hydrated SC showed two lipid endotherms T1 and T2 at 65 and 78 degrees C and protein endotherm T3 at 97 degrees C. All endotherms were absent after SC treated for 48 h with 12 ml of terpene solutions and a decrease in melting points (m.p.) of lipids with a shift of protein endotherm were observed after 12 h treatment with 7 ml of terpene solutions. Linalool and terpineol decreased the m.p. of T1 to 33 degrees C. Carvacrol increased the T1 peak area, which was attributed to lateral lipid bilayer swelling. The IR spectra showed decreases in peak areas and heights of CH2 stretchings but did not show shift of these peaks, increase in their peak widths and shift in amide bands. All the three terpenes disrupted the lipid bilayer and extracted the lipids. Moreover, carvacrol increased the partition of HP whilst linalool and terpineol fluidized the lipids at skin temperature. There could be other possible protein-terpene interactions.     

Experimental and molecular dynamics study into the surfactant effect upon coal wettability

In order to improve the wettability and permeability of coal seams, the water injection efficiency of coal seams has to be boosted, the amount of dust generation has to be reduced, and coal and gas outburst must be prevented, and a surfactant is used to modulate the coal surface wettability. In this work, taking coal samples from Pingdingshan mine in Henan as the research object, their surface chemistry was initially scrutinized and then coal surface engineering via surfactants was inspected by a contact angle test. The coal wettability was ameliorated with surfactants, particularly using the 1 wt% non-ionic surfactant Triton X-100, which elicited a 47% lower contact angle than the raw coal. The surface free energy of the coal sample modified by 1.0 wt% Triton X-100 was increased from 44.51 mN m⁻¹ to 49.52 mN m⁻¹. The microstructural characteristics of coal samples allowed leveraging the Wiser model to construct three kinds of surfactant-coal adsorption models to dissect the adsorption configuration of the system. The results indicate that the addition of surfactants increases both the interaction of water with the coal and the diffusion coefficient of water molecules, resulting in the coal surface transformation from hydrophobicity to hydrophilicity. Our current work can provide salutary guidance and reference for coal water injection and dust suppression.

The experimental and molecular dynamics studies show that the effective wettability of Triton X-100 is controlled by the strong π–π adsorption between hydrophobic end and coal molecule, and the hydrogen bond between hydrophilic end and water.     

Molecular insights into the improved clinical performance of PEGylated interferon therapeutics: a molecular dynamics perspective

PEGylation is a widely adopted process to covalently attach a polyethylene glycol (PEG) polymer to a protein drug for the purpose of optimizing drug clinical performance. While the outcomes of PEGylation in imparting pharmacological advantages have been examined through experimental studies, the underlying molecular mechanisms remain poorly understood. Using interferon (IFN) as a representative model system, we carried out comparative molecular dynamics (MD) simulations of free PEGx, apo-IFN, and PEGx-IFN (x = 50, 100, 200, 300) to characterize the molecular-level changes in IFN introduced by PEGylation. The simulations yielded molecular evidence directly linked to the improved protein stability, bioavailability, retention time, as well as the decrease in protein bioactivity with PEG conjugates. Our results indicate that there is a tradeoff between the benefits and costs of PEGylation. The optimal PEG chain length used in PEGylation needs to strike a good balance among the competing factors and maximizes the overall therapeutic efficacy of the protein drug. We anticipate the study will have a broad implication for protein drug design and development, and provide a unique computational approach in the context of optimizing PEGylated protein drug conjugates.

We discovered molecular evidence that links PEGylation to improved clinical performance, yet at the expense of decreased bioactivity. Our computational approach will facilitate PEGylated protein drug design and optimize its overall therapeutic efficacy.     

Enhancement of skin permeation of high molecular compounds by a combination of microneedle pretreatment and iontophoresis.

A combination of microneedle pretreatment and iontophoresis was evaluated for the potential to increase skin permeation of drugs. Two model compounds with low and high molecular D(2)O and fluorescein isothiocyanate (FITC)-dextrans (FD-4, FD-10, FD-40, FD-70 and FD-2000; average molecular weight of 3.8, 10.1, 39.0, 71.2 and 200.0 kDa), respectively, were used and the effect of microneedle pretreatment and iontophoresis on their in vitro permeability was evaluated using excised hairless rat skin with a 2-chamber diffusion cell. Convective solvent flow through the skin was measured using a set of calibrated capillaries attached to the diffusion cell. The following results were obtained: (1) convective solvent flow (electroosmosis) during iontophoresis through microneedle-pretreated skin, 2.62+/-0.32 microL/cm(2)/h, was almost the same as through intact skin, 2.71+/-0.25 microL/cm(2)/h, and (2) the combination of microneedle pretreatment and subsequent iontophoresis significantly enhanced FD flux compared with microneedle pretreatment alone or iontophoresis alone, whereas no synergistic effect was found on the flux of D(2)O. These results suggest that the combination of iontophoresis with microneedle pretreatment may be a useful means to increase skin permeation of high molecular compounds.     

Nanoscale insight on the durability of magnesium phosphate cement: a molecular dynamics study

The sustainable green building material magnesium phosphate cement (MPC) is widely used in the fields of solidifying heavy metals and nuclear waste and repair and reinforcement. Magnesium potassium phosphate hexahydrate (MKP) is the main hydration product of MPC. The transport of water and ions in MKP nanochannels determines the mechanical properties and durability of MPC materials. Herein, the interface models of MKP crystals with sodium chloride solution in the [001], [010] and [100] direction were established by molecular dynamics. The interaction of the MKP interface with water and ions was studied and the durability of MPC in sodium chloride solution was explained at the molecular level. The results show that a large number of water molecules are adsorbed on the MKP crystal surface through hydrogen bonds and Coulomb interactions; the surface water molecules have the bigger dipole moment and the dipole vector of most of the water molecules points to the solid matrix, when the crystal surfaces of the three models all show hydrophilicity. In addition, plenty of sodium ions are adsorbed at the MKP interface, and some potassium ions are desorbed from the matrix. In the MKP[001] model, the amount of potassium ions separated from the matrix and diffused into the solution is the highest and the interface crystal is the most disordered. Due to the attack of water and ions, the K–Os bond loses its chemical stability and the order of the MKP crystal is destroyed, which explains the decline of MPC performance after the erosion of sodium chloride solution at the molecular level. Besides, in the three models, the Na–Cl ion bond is more unstable than the K–Cl ion bond due to the smaller radius of the sodium atom. The stability of ionic bonds in the models is as follows: MKP[010] > MKP[100] > MKP[001].

The sustainable green building material magnesium phosphate cement (MPC) is widely used in the fields of solidifying heavy metals and nuclear waste and repair and reinforcement.     

The effect of surfactants on hydrate particle agglomeration in liquid hydrocarbon continuous systems: a molecular dynamics simulation study

Anti-agglomerants (AAs), both natural and commercial, are currently being considered for gas hydrate risk management of petroleum pipelines in offshore operations. However, the molecular mechanisms of the interaction between the AAs and gas hydrate surfaces and the prevention of hydrate agglomeration remain critical and complex questions that need to be addressed to advance this technology. Here, we use molecular dynamics (MD) simulations to investigate the effect of model surfactant molecules (polynuclear aromatic carboxylic acids) on the agglomeration behaviour of gas hydrate particles and disruption of the capillary liquid bridge between hydrate particles. The results show that the anti-agglomeration pathway can be divided into two processes: the spontaneous adsorption effect of surfactant molecules onto the hydrate surface and the weakening effect of the intensity of the liquid bridge between attracted hydrate particles. The MD simulation results also indicate that the anti-agglomeration effectiveness of surfactants is determined by the intrinsic nature of their molecular functional groups. Additionally, we find that surfactant molecules can affect hydrate growth, which decreases hydrate particle size and correspondingly lower the risk of hydrate agglomeration. This study provides molecular-level insights into the anti-agglomeration mechanism of surfactant molecules, which can aid in the ultimate application of natural or commercial AAs with optimal anti-agglomeration properties.

Schematic of anti-agglomeration effect of surfactants promoting gas hydrate particle dispersion.     

Self-assembly of glycerol monooleate with the antimicrobial peptide LL-37: a molecular dynamics study

Over the past decade, the rapid increase in the incidence of antibiotic-resistant bacteria has promoted research towards alternative therapeutics such as antimicrobial peptides (AMPs), but their biodegradability limits their application. Encapsulation into nanocarriers based on the self-assembly of surfactant-like lipids is emerging as a promising strategy for the improvement of AMPs'' stability and their protection against degradation when in biological media. An in-depth understanding of the interactions between the structure-forming lipids and AMPs is required for the design of nanocarriers. This in silico study, demonstrates the self-assembly of the amphiphilic lipid glycerol monooleate (GMO) with the antimicrobial peptide LL-37 into nanocarriers on the molecular scale. Molecular dynamics (MD) simulations show the formation of direct micelles, with either one or two interacting LL-37, and vesicles in this two-component system in agreement with experimental results from small-angle X-ray scattering studies. The hydrophobic contacts between LL-37 and GMOs in water appear responsible for the formation of these nanoparticles. The results also suggest that the enhanced antimicrobial efficiency of LL-37 in these nanocarriers that was previously observed experimentally can be explained by the availability of its side chains with charged amino acids, an increase of the electrostatic interaction and a decrease of the peptide''s conformational entropy upon interacting with GMO. The results of this study contribute to the fundamental understanding of lipid–AMP interactions and may guide the comprehensive design of lipid-based self-assembled nanocarriers for antimicrobial peptides.

Molecular dynamics simulations of glycerol-monooleate (GMO)/LL-37 nanocarriers show that hydrophobic interactions among the molecules drive the formation of GMO/LL-37 micelles.     

Tuning the mechanical properties of silicene nanosheet by auxiliary cracks: a molecular dynamics study

Silicene has become a topic of interest nowadays due to its potential application in various electro-mechanical nanodevices. In our previous work on silicene, fracture stresses of single crystal and polycrystalline silicene have been investigated. Existence of defects in the form of cracks reduces the fracture strength of silicene nanosheets to a great extent. In this study, an engineering way has been proposed for improving the fracture stress of silicene nanosheets with a pre-existing crack by incorporating auxiliary cracks symmetrically in a direction perpendicular to the main crack. We call this mechanism the “Failure shielding mechanism”. An extensive molecular dynamics simulation based analysis has been performed to capture the atomic level auxiliary crack-main crack interactions. It is found that the main crack tip stress distribution is significantly changed with the presence of auxiliary cracks for loading along both armchair and zigzag directions. The effects of temperature and the crack propagation speed of silicene have also been studied. Interestingly, in the case of loading along the zigzag direction, SW defect formation is observed at the tip of main crack. This leads to a reduction of the tip stress resulting in a more prominent failure shielding in case of zigzag loading than in armchair loading. Moreover, the position and length of the cracks as well as the loading directions have significant impacts on the tip stress distribution. Finally, this study opens the possibilities of strain engineering for silicene by proposing an engineering way to tailor the fracture strength of silicene.

Inclusion of auxiliary cracks increases the fracture stress of silicene nanosheets with a pre existing crack.     

Conformational dynamics of superoxide dismutase (SOD1) in osmolytes: a molecular dynamics simulation study

Amyotrophic lateral sclerosis (ALS) is a progressive neurodegenerative disease caused by the misfolding of Cu, Zn superoxide dismutase (SOD1). Several earlier studies have shown that monomeric apo SOD1 undergoes significant local unfolding dynamics and is the predecessor for aggregation. Here, we have employed atomistic molecular dynamics (MD) simulations to study the structure and dynamics of monomeric apo and holo SOD1 in water, aqueous urea and aqueous urea–TMAO (trimethylamine oxide) solutions. Loop IV (zinc-binding loop) and loop VII (electrostatic loop) of holo SOD1 are considered as functionally important loops as they are responsible for the structural stability of holo SOD1. We found larger local unfolding of loop IV and VII of apo SOD1 as compared to holo SOD1 in water. Urea induced more unfolding in holo SOD1 than apo SOD1, whereas the stabilization of both the form of SOD1 was observed in ternary solution (i.e. water/urea/TMAO solution) but the extent of stabilization was higher in holo SOD1 than apo SOD1. The partially unfolded structures of apo SOD1 in water, urea and holo SOD1 in urea were identified by the exposure of the hydrophobic cores, which are highly dynamic and these may be the initial events of aggregation in SOD1. Our simulation studies support the formation of aggregates by means of the local unfolding of monomeric apo SOD1 as compared to holo SOD1 in water.

Change in conformations of apo and holo SOD1 in water and in osmolytes in terms of configurational entropy (S).     

Insights into the molecular mechanism underlying CD4-dependency and neutralization sensitivity of HIV-1: a comparative molecular dynamics study on gp120s from isolates with different phenotypes

The envelope (Env) of HIV-1 plays critical roles in viral infection and immune evasion. Although structures of prefusion Env have been determined and phenotypes relevant to the CD4 dependency and the neutralization sensitivity for various HIV-1 isolates have been identified, the detailed structural dynamics and energetics underlying these two phenotypes have remained elusive. In this study, two unliganded structural models of gp120, one from the CD4-dependent, neutralization-resistant isolate H061.14 and the other from the CD4-independent, neutralization-sensitive R2 strain, were constructed, and subsequently were subjected to multiple-replica molecular dynamics (MD) simulations followed by free energy landscape (FEL) construction. Comparative analyses of MD trajectories reveal that during simulations R2-gp120 demonstrated larger structural fluctuations/deviations and higher global conformational flexibility than H061.14-gp120. Close comparison of local conformational flexibility shows that some of the structural regions involving direct interactions with gp41 and adjacent gp120 subunits in the context of the closed trimeric Env exhibit significantly higher flexibility in R2-gp120 than in H061.14-gp120, thus likely increasing the probability for R2-Env to open the trimer crown and prime gp41 fusogenic properties without induction by CD4. Collective motions derived from principal component analysis (PCA) reveal that R2-gp120 is prone to spontaneous transition to the neutralization-sensitive CD4-bound state while H061.14-gp120 tends to maintain the neutralization-resistant unliganded state. Finally, comparison between FELs reveals that R2-gp120 has larger conformational entropy, richer conformational diversity, and lower thermostability than H061.14-gp120, thus explaining why R2-gp120 is more structurally unstable and conformationally flexible, and has a higher propensity to transition to the CD4-bound state than H061.14-gp120. The present results reveal that the differences in dynamics and energetics between R2-gp120 and H061.14-gp120 impart Env trimers with distinct capacities to sample different states (i.e., R2-Env samples more readily the open state while H061.14-Env is more inclined to maintain the closed state), thus shedding light on the molecular mechanism underlying the HIV-1 phenotype associated with CD4 dependency/neutralization sensitivity.

The envelope (Env) of HIV-1 plays critical roles in viral infection and immune evasion.     

Tunable ductility of a nano-network from few-layered graphene bonded with benzene: a molecular dynamics study

Developing novel graphene-based materials with unique mechanical properties is of significance to meet the requirements in new applications. The pristine graphene shows a brittle fracture when the stretching strain on it exceeds the critical value. Further, it fails to bear the external load. Herein, to enhance the ductility of the pristine graphene, we proposed a corrugated sandwich carbon network based on few-layered graphene, in which the two surface layers are bonded with several corrugated core layers via benzene molecules. The effects of factors such as the geometry, temperature, and strain rate on the ductility of the carbon network were evaluated using the uniaxial tension tests by molecular dynamics simulations. Results show that the new carbon material has more than one peak fracture strain in stretching. The second peak fracture strain is proportional to the length difference between the surface layers and core layers. Hence, the carbon network has a tunable ductility, which suggests a flexible design of such novel materials in a nanostructure/nanodevice with large deformation.

The ductility of the corrugated sandwich carbon nano-network is tunable and higher than that of the pristine graphene.     

SMGA gels for the skin permeation of haloperidol.

Small molecule gelling agent (SMGA) gels were developed using the gelator GP-1 in the solvents, namely, isostearyl alcohol (ISA) and propylene glycol (PG), to deliver haloperidol through the skin. The concentrations of the drug, haloperidol, the enhancer, farnesol and the gelator, GP-1 are 3 mg/ml, 5% (w/v) and 5% (w/v), respectively. The study employed a three-factor full factorial statistical design to investigate the influence of factor level changes on the permeability coefficient and permeation lag-time of haloperidol. Gels were prepared by raising temperature to 120 degrees C, followed by natural cooling under room temperature of 22+/-1 degrees C. The rheological properties of the gels were examined with a strain-controlled dynamic mechanical method. The in vitro permeation study was conducted with automated flow-through type cells. The gels successfully incorporated the drug and enhancer without losing their aesthetic properties. The in vitro human skin permeation study showed the permeation of the drug in ISA-based gels reached the pseudo steady state faster than PG-based gels and were less affected by gelator. PG-based gels delivered the drug at a faster rate with the incorporation of the enhancer. GP-1 did not influence the drug permeation rate but it increased permeation lag-time. The co-existence of gelator or enhancer increased the lag-time to a larger extent than when used separately. The novel SMGA gels are suitable for topical or transdermal delivery.     

Molecular dynamics simulation of droplet nucleation and growth on a rough surface: revealing the microscopic mechanism of the flooding mode

Droplet nucleation and growth have a significant influence on dropwise condensation heat transfer. Controlling the droplet nucleation and growth with a high-precision surface to realize the dropwise condensation heat transfer enhancement is a promising method. Molecular dynamics simulation is employed to investigate the effects of heat flux, surface wettability and adjacent droplet size on the new droplet''s nucleation and growth. Simulation results indicate that the high heat flux can lead to droplet nucleation and growth inside the rough structure and finally a Wenzel droplet will form due to the coalescence between the inside droplet and the initial existing droplet. However, for a surface with a larger contact angle, the droplet in the Wenzel state will transfer to the Cassie state due to droplet coalescence. In addition, it is also related to the size of the existing droplet whether or not the nucleation process occurs. For the first time, the droplet nucleation radius is introduced to quantitatively determine the droplet nucleation state (inside or outside the nanostructures) and whether the droplet could achieve the state transition from Wenzel to Cassie or not in the growth process.

Droplet nucleation and growth have a significant influence on dropwise condensation heat transfer.     

Enhanced skin permeation of a lipophilic drug using supersaturated formulations.

Supersaturation was used to enhance the permeation of a lipophilic model compound (a lavendustin derivative, LAP) through excised pig skin in vitro. The drug was dissolved in a series of liquid and semisolid vehicles (in which it had different solubilities) and which were prepared using either (i) the method of mixed cosolvents, (ii) the method of solvent evaporation, or (iii) the method of dissolving the drug with heating. Saturated formulations showed comparable permeation rates through the skin, independent of the absolute concentration of the drug in the vehicle. Supersaturated solutions at a degree of saturation of two resulted in a doubling of the drug permeation rate. These experiments show, therefore, that the percutaneous absorption of LAP may be consistently increased using supersaturated formulations, independent of the type and composition of the vehicles and independent of their method of preparation.     

自制脂质体超声纳米微泡的特性和显影效果

目的　测定自制纳米级脂质体超声微泡造影剂(nano-liposomal bubble,NB)的基本特性.方法　采用逆向蒸发法制备脂质体,用声震法制备NB后观察其形态,检测其平均粒径、表面电位及浓度等物理特性;在脱气水中分别注入生理盐水、NB及SonoVue (200μl),观察显影效果;分别经兔耳缘静脉团注生理盐水、NB及SonoVue (2.0 ml/kg),动态观察兔心脏及肝的增强情况.结果　镜下观察,NB粒径范围133.1～199.5 nm,平均粒径为(171.60±30.82)nm,平均电位为- (1.92±0.65)mV,分布均匀,浓度为(3.8～5.6)×108/ml.常温下造影剂放置1周、1月后观察,基本特性无明显改变.体外显影结果显示:NB与SonoVue一样具有显著的显影效果.体内造影后观察:NB与SonoVue均能明显增强兔心脏及肝的二维灰阶显像.结论　NB性质稳定,形态好,粒径均一,具有良好的显影效果.由于其粒径小且基于脂质体基础上易被修饰,因此在超声分子影像及基因/药物传递方面具有广阔的应用前景.     

Shedding light on the structural properties of lipid bilayers using molecular dynamics simulation: a review study

In recent years, a massive increase has been observed in the number of published articles describing accurate and reliable molecular dynamics simulations of lipid bilayers. This is due to several reasons, including the development of fast and efficient methods for treating long-range electrostatic interactions, significant progress in computer hardware, progress in atomistic simulation algorithms and the development of well-validated empirical molecular mechanical force fields. Although molecular dynamics is an effective approach for investigating different aspects of lipid bilayers, to the best of our knowledge, there is no review in the literature that explains the different analyses that can be carried out with membrane simulation. This review gives an overview about the some of the most important possible analyses, technical challenges, and existing protocols that can be performed on the biological membrane by molecular dynamics simulation. The reviewed analyses include the degree of membrane disruption, average area per lipid, probability distributions for the area per lipid molecule, membrane thickness, membrane area compressibility, lateral diffusion, rotational diffusion, order parameters, head group tilt, electron density profile, mass density profile, electrostatic potential profile, ordering of vicinity waters, number of hydrogen bonds, and radial distribution function.

This review gives an overview about the some of the most important possible analyzes, technical challenges, and existing protocols that can be performed on the biological membrane by the molecular dynamics simulation.     

Comparison of the effect of fatty alcohols on the permeation of melatonin between porcine and human skin.

Melatonin (MT) is a hormone secreted by the pineal gland that plays an important role in the regulation of the circadian sleep-wake cycle. It would be advantageous to administer MT using a transdermal delivery system for the treatment of sleep disorders such as delayed sleep syndrome, jet lag in travelers, cosmonauts and shift workers. The porcine skin has been found to have similar morphological and functional characteristics as human skin. The elastic fibres in the dermis, enzyme pattern of the epidermis, epidermal tissue turnover time, keratinous proteins and thickness of epidermis of porcine skin are similar to human skin. However, the fat deposition and vascularisation of the cutaneous glands of porcine skin are different from human skin. In addition, porcine skin has been found to have a close permeability character to human skin. However, the comparative effect of chemical penetration enhancers on the permeation of drugs between porcine and human skin has not been reported. The purpose of this study was to compare the effect of fatty alcohols on the permeability of porcine and human skin using MT as a model compound. The effect of saturated fatty alcohols (octanol, nonanol, decanol, undecanol, lauryl alcohol, tridecanol, myristyl alcohol) and unsaturated fatty alcohols (oleyl alcohol, linoleyl alcohol, linolenyl alcohol) at 5% concentration was tested across dermatomed porcine and human skin. Our studies showed a parabolic relationship between the carbon chain length of saturated fatty alcohols and permeation enhancement of MT with both porcine and human skin. Maximum permeation of MT was observed when fatty alcohol carbon chain length was 10. In general, as the level of unsaturation increased from one to two double bonds, there was an increase in the permeation of MT both in porcine and human skin. However, a decrease in the permeation was observed with three double bonds. Regression analysis using the steady state flux data showed a significant positive correlation between porcine and human skin for saturated fatty alcohols (r(2)=0.8868, P=0.0005). However, though a positive correlation was observed between the porcine and human skin (r(2)=0.8638), the correlation was statistically insignificant (P=0.0706). The static diffusion cell system employed in this study has major artifact compared to a flow through system. In conclusion, the permeability of porcine skin to MT in the presence of saturated and unsaturated fatty alcohols was qualitatively similar to human skin but quantitatively different with some fatty alcohols.