Browse
Submit Data
Databases
API
Help

Dataset Information

14 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Evaluation of Constrained and Restrained Molecular Dynamics Simulation Methods for Predicting Skin Lipid Permeability.

ABSTRACT: Recently, molecular dynamics (MD) simulations have been utilized to investigate the barrier properties of human skin stratum corneum (SC) lipid bilayers. Different MD methods and force fields have been utilized, with predicted permeabilities varying by few orders of magnitude. In this work, we compare constrained MD simulations with restrained MD simulations to obtain the potential of the mean force and the diffusion coefficient profile for the case of a water molecule permeating across an SC lipid bilayer. Corresponding permeabilities of the simulated lipid bilayer are calculated via the inhomogeneous solubility diffusion model. Results show that both methods perform similarly, but restrained MD simulations have proven to be the more robust approach for predicting the potential of the mean force profile. Critical to both methods are the sampling of the whole trans-bilayer axis and the following symmetrization process. Re-analysis of the previously reported free energy profiles showed that some of the discrepancies in the reported permeability values is due to misquotation of units, while some are due to the inaccurately obtained potential of the mean force. By using the existing microscopic geometrical models via the intercellular lipid pathway, the permeation through the whole SC is predicted from the MD simulation results, and the predicted barrier properties have been compared to experimental data from the literature with good agreement.

SUBMITTER: Piasentin N

PROVIDER: S-EPMC8717400 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Predicting chemically-induced skin reactions. Part II: QSAR models of skin permeability and the relationships between skin permeability and skin sensitization.

Project description:Skin permeability is widely considered to be mechanistically implicated in chemically-induced skin sensitization. Although many chemicals have been identified as skin sensitizers, there have been very few reports analyzing the relationships between molecular structure and skin permeability of sensitizers and non-sensitizers. The goals of this study were to: (i) compile, curate, and integrate the largest publicly available dataset of chemicals studied for their skin permeability; (ii) develop and rigorously validate QSAR models to predict skin permeability; and (iii) explore the complex relationships between skin sensitization and skin permeability. Based on the largest publicly available dataset compiled in this study, we found no overall correlation between skin permeability and skin sensitization. In addition, cross-species correlation coefficient between human and rodent permeability data was found to be as low as R(2)=0.44. Human skin permeability models based on the random forest method have been developed and validated using OECD-compliant QSAR modeling workflow. Their external accuracy was high (Q(2)ext=0.73 for 63% of external compounds inside the applicability domain). The extended analysis using both experimentally-measured and QSAR-imputed data still confirmed the absence of any overall concordance between skin permeability and skin sensitization. This observation suggests that chemical modifications that affect skin permeability should not be presumed a priori to modulate the sensitization potential of chemicals. The models reported herein as well as those developed in the companion paper on skin sensitization suggest that it may be possible to rationally design compounds with the desired high skin permeability but low sensitization potential.

| S-EPMC4408226 | biostudies-literature

Fragment contribution models for predicting skin permeability using HuskinDB.

Project description:Mathematical models to predict skin permeation tend to be based on animal derived experimental data as well as knowing physicochemical properties of the compound under investigation, such as molecular volume, polarity and lipophilicity. This paper presents a strikingly contrasting model to predict permeability, formed entirely from simple chemical fragment (functional group) data and a recently released, freely accessible human (i.e. non-animal) skin permeation database, known as the 'Human Skin Database - HuskinDB'. Data from within the database allowed development of several fragment-based models, each including a calculable effect for all of the most commonly encountered functional groups present in compounds within the database. The developed models can be applied to predict human skin permeability (logKp) for any compound containing one or more of the functional groups analysed from the dataset with no need to know any other physicochemical properties, solely the type and number of each functional group within the chemical structure itself. This approach simplifies mathematical prediction of permeability for compounds with similar properties to those used in this study.

| S-EPMC10667307 | biostudies-literature

Skin permeability prediction with MD simulation sampling spatial and alchemical reaction coordinates.

Project description:A molecular-level understanding of skin permeation may rationalize and streamline product development, and improve quality and control, of transdermal and topical drug delivery systems. It may also facilitate toxicity and safety assessment of cosmetics and skin care products. Here, we present new molecular dynamics simulation approaches that make it possible to efficiently sample the free energy and local diffusion coefficient across the skin's barrier structure to predict skin permeability and the effects of chemical penetration enhancers. In particular, we introduce a new approach to use two-dimensional reaction coordinates in the accelerated weight histogram method, where we combine sampling along spatial coordinates with an alchemical perturbation virtual coordinate. We present predicted properties for 20 permeants, and demonstrate how our approach improves correlation with ex vivo/in vitro skin permeation data. For the compounds included in this study, the obtained log KPexp-calc mean square difference was 0.9 cm2 h-2.

| S-EPMC9674988 | biostudies-literature

Molecular-dynamics simulation methods for macromolecular crystallography.

Project description:It is investigated whether molecular-dynamics (MD) simulations can be used to enhance macromolecular crystallography (MX) studies. Historically, protein crystal structures have been described using a single set of atomic coordinates. Because conformational variation is important for protein function, researchers now often build models that contain multiple structures. Methods for building such models can fail, however, in regions where the crystallographic density is difficult to interpret, for example at the protein-solvent interface. To address this limitation, a set of MD-MX methods that combine MD simulations of protein crystals with conventional modeling and refinement tools have been developed. In an application to a cyclic adenosine monophosphate-dependent protein kinase at room temperature, the procedure improved the interpretation of ambiguous density, yielding an alternative water model and a revised protein model including multiple conformations. The revised model provides mechanistic insights into the catalytic and regulatory interactions of the enzyme. The same methods may be used in other MX studies to seek mechanistic insights.

| S-EPMC9815100 | biostudies-literature

Permeability Barrier and Microstructure of Skin Lipid Membrane Models of Impaired Glucosylceramide Processing.

Project description:Ceramide (Cer) release from glucosylceramides (GlcCer) is critical for the formation of the skin permeability barrier. Changes in β-glucocerebrosidase (GlcCer'ase) activity lead to diminished Cer, GlcCer accumulation and structural defects in SC lipid lamellae; however, the molecular basis for this impairment is not clear. We investigated impaired GlcCer-to-Cer processing in human Cer membranes to determine the physicochemical properties responsible for the barrier defects. Minor impairment (5-25%) of the Cer generation from GlcCer decreased the permeability of the model membrane to four markers and altered the membrane microstructure (studied by X-ray powder diffraction and infrared spectroscopy), in agreement with the effects of topical GlcCer in human skin. At these concentrations, the accumulation of GlcCer was a stronger contributor to this disturbance than the lack of human Cer. However, replacement of 50-100% human Cer by GlcCer led to the formation of a new lamellar phase and the maintenance of a rather good barrier to the four studied permeability markers. These findings suggest that the major cause of the impaired water permeability barrier in complete GlcCer'ase deficiency is not the accumulation of free GlcCer but other factors, possibly the retention of GlcCer bound in the corneocyte lipid envelope.

| S-EPMC5527096 | biostudies-literature

Solution NMR refinement of a metal ion bound protein using metal ion inclusive restrained molecular dynamics methods.

Project description:Correctly calculating the structure of metal coordination sites in a protein during the process of nuclear magnetic resonance (NMR) structure determination and refinement continues to be a challenging task. In this study, we present an accurate and convenient means by which to include metal ions in the NMR structure determination process using molecular dynamics (MD) simulations constrained by NMR-derived data to obtain a realistic and physically viable description of the metal binding site(s). This method provides the framework to accurately portray the metal ions and its binding residues in a pseudo-bond or dummy-cation like approach, and is validated by quantum mechanical/molecular mechanical (QM/MM) MD calculations constrained by NMR-derived data. To illustrate this approach, we refine the zinc coordination complex structure of the zinc sensing transcriptional repressor protein Staphylococcus aureus CzrA, generating over 130 ns of MD and QM/MM MD NMR-data compliant sampling. In addition to refining the first coordination shell structure of the Zn(II) ion, this protocol benefits from being performed in a periodically replicated solvation environment including long-range electrostatics. We determine that unrestrained (not based on NMR data) MD simulations correlated to the NMR data in a time-averaged ensemble. The accurate solution structure ensemble of the metal-bound protein accurately describes the role of conformational sampling in allosteric regulation of DNA binding by zinc and serves to validate our previous unrestrained MD simulations of CzrA. This methodology has potentially broad applicability in the structure determination of metal ion bound proteins, protein folding and metal template protein-design studies.

| S-EPMC3773525 | biostudies-literature

Evaluation of regression methods when immunological measurements are constrained by detection limits.

Project description:BackgroundThe statistical analysis of immunological data may be complicated because precise quantitative levels cannot always be determined. Values below a given detection limit may not be observed (nondetects), and data with nondetects are called left-censored. Since nondetects cannot be considered as missing at random, a statistician faced with data containing these nondetects must decide how to combine nondetects with detects. Till now, the common practice is to impute each nondetect with a single value such as a half of the detection limit, and to conduct ordinary regression analysis. The first aim of this paper is to give an overview of methods to analyze, and to provide new methods handling censored data other than an (ordinary) linear regression. The second aim is to compare these methods by simulation studies based on real data.ResultsWe compared six new and existing methods: deletion of nondetects, single substitution, extrapolation by regression on order statistics, multiple imputation using maximum likelihood estimation, tobit regression, and logistic regression. The deletion and extrapolation by regression on order statistics methods gave biased parameter estimates. The single substitution method underestimated variances, and logistic regression suffered loss of power. Based on simulation studies, we found that tobit regression performed well when the proportion of nondetects was less than 30%, and that taken together the multiple imputation method performed best.ConclusionBased on simulation studies, the newly developed multiple imputation method performed consistently well under different scenarios of various proportion of nondetects, sample sizes and even in the presence of heteroscedastic errors.

| S-EPMC2592244 | biostudies-literature

Molecular dynamics simulation of lipid reorientation at bilayer edges.

Project description:Understanding cellular membrane processes is critical for the study of events such as viral entry, neurotransmitter exocytosis, and immune activation. Supported lipid bilayers are commonly used to model these membrane processes experimentally. Despite the relative simplicity of such a system, many important structural and dynamic parameters are not experimentally observable with current techniques. Computational approaches allow the development of a high-resolution model of bilayer processes. We have performed molecular dynamics simulations of dimyristoylphosphatidylcholine (DMPC) bilayers to model the creation of bilayer gaps-a common process in bilayer patterning-and to analyze their structure and dynamics. We propose a model for gap formation in which the bilayer edges form metastable micelle-like structures on a nanosecond timescale. Molecules near edges structurally resemble lipids in ungapped bilayers but undergo small-scale motions more rapidly. These data suggest that lipids may undergo rapid local rearrangements during membrane fusion, facilitating the formation of fusion intermediates thought key to the infection cycle of viruses such as influenza, Ebola, and HIV.

| S-EPMC1304275 | biostudies-literature

Molecular dynamic simulations of oxidized skin lipid bilayer and permeability of reactive oxygen species.

Project description:Lipid peroxidation by reactive oxygen species (ROS) during oxidative stress is non-enzymatic damage that affects the integrity of biological membrane, and alters the fluidity and permeability. We conducted molecular dynamic simulation studies to evaluate the structural properties of the bilayer after lipid peroxidation and to measure the permeability of distinct ROS. The oxidized membrane contains free fatty acid, ceramide, cholesterol, and 5α-hydroperoxycholesterol (5α-CH). The result of unconstrained molecular dynamic simulations revealed that lipid peroxidation causes area-per-lipid of the bilayer to increase and bilayer thickness to decrease. The simulations also revealed that the oxidized group of 5α-CH (-OOH) moves towards the aqueous layer and its backbone tilts causing lateral expansion of the bilayer membrane. These changes are detrimental to structural and functional properties of the membrane. The measured free energy profile for different ROS (H2O2, HO2, HO, and O2) across the peroxidized lipid bilayer showed that the increase in lipid peroxidation resulted in breaching barrier decrease for all species, allowing easy traversal of the membrane. Thus, lipid peroxidation perturbs the membrane barrier and imposes oxidative stress resulting into apoptosis. The collective insights increase the understanding of oxidation stress at the atomic level.

| S-EPMC6418262 | biostudies-literature

Predicting the Membrane Permeability of Fentanyl and Its Analogues by Molecular Dynamics Simulations.

Project description:The lipid membrane is considered a crucial component of opioid general anesthesia. The main drug used for the induction and maintenance of opioid anesthesia is fentanyl and its various analogues. However, these drugs have different clinical effects, and detailed atomic-level insight into the drug-membrane interactions could lead to a better understanding how these drugs exert their anesthetic properties. In this study, we have used extensive umbrella sampling molecular dynamics simulations to study the permeation process of fentanyl and three of its analogues into a variety of simple phospholipid membrane models. Our simulations show that we can accurately predict the permeability coefficients of these drug molecules, which is an important process in understanding how pharmaceuticals reach their molecular targets. We were also able to show that one phospholipid provides more accurate predictions than other lipids commonly used in these types of permeation studies, which will aid future studies of these types of processes.

| S-EPMC8389899 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data