Abstract P-2: Multiscale Molecular Dynamics Simulations of the Skin Membranes as a Tool for the Interpretation of the Transmission Electron Microscopy Images

Background: Human skin can inhibit chemical penetration which limits the clinical applications of transdermal drug delivery. The stratum corneum (SC) is the primary barrier and organized in lamellar membranes containing the lipids of ceramides (CER), free fatty acids (FFA), and cholesterol (CHOL). One of the most widely used ways to overcome the SC is the addition of chemical penetration enhancers (CPEs) to active ingredients. There are various methods, which have been employed to explore the mechanisms by which CPEs with drugs can change the morphology of SC including transmission electron microscopy. Here, we propose to use multiscale coarse-grained (CG) molecular dynamics (MD) simulations for the interpretation of the images of the SC from the electron microscopy experiments. Methods: We utilized the MARTINI force field for the CG simulations. We employed the mixed-lipid bilayer model of SC consisting of CER, CHOL, and FFA in a 1:1:1 molar ratio assembled with CHARMM-GUI web-service. The systems of the SC model membrane and various enhancers were simulated in the NPT ensemble with the polarizable water model and the reaction field approach for the long-range electrostatics with the usage of Gromacs 2019.4 software. Results: The membrane model was validated with standard characteristics: thickness, diffusion of the lipids, order parameters, and density profiles. After, we have added CPEs and active ingredients to the systems: menthol and osthole as control simulations, ethanol with linoleic acid and lidocaine as test simulations. We have observed the membrane desegregation in the case of menthol and osthole formulations similar to the published results while the permeation of lidocaine with ethanol and linoleic acid did not cause the disruption of the membranes but increased its fluidity and permeability properties. Conclusion: The method of multiscale coarse-grained molecular dynamics simulations can be utilized for the prediction and interpretation of morphology change of SC in addition to various substances

Stability and Radiation Damage of Protein Crystals as Studied by Means of Molecular Dynamics and Monte Carlo Simulation

Molecular Dynamics (MD) and Monte Carlo (MC) simulations of crystals can help in interpretation of experimental X-ray crystallography data. Particularly, they can be useful for understanding how various crystallization techniques affect protein conformational plasticity within the crystal lattice and the stability of biomolecular crystals. The latter has become especially important since the modern and extremely intense X-ray radiation sources (such as free electron lasers, FELs) appeared recently. In the present study we were able to show by means of computer simulations that the lysozyme crystals obtained using the Langmuir-Blodgett technique have an advantage over the classical ones (\u201cHanging Drop\u201d) in terms of their thermal stability as well as their stability against the radiation damage. We also demonstrate an important role of crystal water dynamics for stability of protein crystals

Effects of Mutations in the Receptor-Binding Domain of SARS-CoV-2 Spike on its Binding Affinity to ACE2 and Neutralizing Antibodies Revealed by Computational Analysis

SARS-CoV-2 causing coronavirus disease 2019 (COVID-19) is responsible for one of the most deleterious pandemics of our time. The interaction between the ACE2 receptors at the surface of human cells and the viral Spike (S) protein triggers the infection making the receptor-binding domain (RBD) of the SARS-CoV-2 S-protein a focal target for the neutralizing antibodies (Abs). Despite the recent progress in the development and deployment of vaccines, the emergence of novel variants of SARS-CoV-2 insensitive to Abs produced in response to the vaccine administration and/or monoclonal ones represents upcoming jeopardy. Here, we assessed the possible effects of single and multiple mutations in the RBD of SARS-CoV-2 S-protein on its binding energy to various antibodies and the human ACE2 receptor. The performed computational analysis indicates that while single amino acid replacements in RBD may only cause partial impairment of the Abs binding, moreover, limited to specific epitopes, some variants of SARS-CoV-2 (with as few as 8 mutations), which are already present in the population, may potentially result in a much broader antigenic escape. We also identified a number of point mutations, which, in contrast to the majority of replacements, reduce RBD affinity to various antibodies without affecting its binding to ACE2. Overall, the results provide guidelines for further experimental studies aiming at the identification of the high-risk RBD mutations allowing for an antigenic escape

Synergistic Effect of Chemical Penetration Enhancers on Lidocaine Permeability Revealed by Coarse-Grained Molecular Dynamics Simulations

The search for new formulations for transdermal drug delivery (TDD) is an important field in medicine and cosmetology. Molecules with specific physicochemical properties which can increase the permeability of active ingredients across the stratum corneum (SC) are called chemical penetration enhancers (CPEs), and it was shown that some CPEs can act synergistically. In this study, we performed coarse-grained (CG) molecular dynamics (MD) simulations of the lidocaine delivery facilitated by two CPEs—linoleic acid (LA) and ethanol—through the SC model membrane containing cholesterol, N-Stearoylsphingosine (DCPE), and behenic acid. In our simulations, we probed the effects of individual CPEs as well as their combination on various properties of the SC membrane and the lidocaine penetration across it. We demonstrated that the addition of both CPEs decreases the membrane thickness and the order parameters of the DPCE hydrocarbon chains. Moreover, LA also enhances diffusion of the SC membrane components, especially cholesterol. The estimated potential of mean force (PMF) profiles for the lidocaine translocation across SC in the presence/absence of two individual CPEs and their combination demonstrated that while ethanol lowers the free energy barrier for lidocaine to enter SC, LA decreases the depth of the free energy minima for lidocaine inside SC. These two effects supposedly result in synergistic penetration enhancement of drugs. Altogether, the present simulations provide a detailed molecular picture of CPEs’ action and their synergistic effect on the penetration of small molecular weight therapeutics that can be beneficial for the design of novel drug and cosmetics formulations

Computational Analysis of Mutations in the Receptor-Binding Domain of SARS-CoV-2 Spike and Their Effects on Antibody Binding

Currently, SARS-CoV-2 causing coronavirus disease 2019 (COVID-19) is responsible for one of the most deleterious pandemics of our time. The interaction between the ACE2 receptors at the surface of human cells and the viral Spike (S) protein triggers the infection, making the receptor-binding domain (RBD) of the SARS-CoV-2 S-protein a focal target for the neutralizing antibodies (Abs). Despite the recent progress in the development and deployment of vaccines, the emergence of novel variants of SARS-CoV-2 insensitive to Abs produced in response to the vaccine administration and/or monoclonal ones represent a potential danger. Here, we analyzed the diversity of neutralizing Ab epitopes and assessed the possible effects of single and multiple mutations in the RBD of SARS-CoV-2 S-protein on its binding affinity to various antibodies and the human ACE2 receptor using bioinformatics approaches. The RBD-Ab complexes with experimentally resolved structures were grouped into four clusters with distinct features at sequence and structure level. The performed computational analysis indicates that while single amino acid replacements in RBD may only cause partial impairment of the Abs binding, moreover, limited to specific epitopes, the variants of SARS-CoV-2 with multiple mutations, including some which were already detected in the population, may potentially result in a much broader antigenic escape. Further analysis of the existing RBD variants pointed to the trade-off between ACE2 binding and antigenic escape as a key limiting factor for the emergence of novel SAR-CoV-2 strains, as the naturally occurring mutations in RBD tend to reduce its binding affinity to Abs but not to ACE2. The results provide guidelines for further experimental studies aiming to identify high-risk RBD mutations that allow for an antigenic escape

Comparison of available computational studies of the interaction of C₆₀ with lipid bilayer.

<p>Comparison of available computational studies of the interaction of C₆₀ with lipid bilayer.</p

Properties of the DPPC membrane in simulations with different amount of fullerenes.

<p>Properties of the DPPC membrane in simulations with different amount of fullerenes.</p

D₃ stereoisomer of tris-malonic fullereneB

<p>A. Free energy profile of the process of the D₃ penetration into the model eukaryotic membrane. B. Orientation (corresponding to the global energy minimum of the free energy profile) of the D₃ molecule adsorbed to the membrane.</p

Comparative Computational Study of Interaction of C₆₀-Fullerene and Tris-Malonyl-C₆₀-Fullerene Isomers with Lipid Bilayer: Relation to Their Antioxidant Effect

<div><p>Oxidative stress induced by excessive production of reactive oxygen species (ROS) has been implicated in the etiology of many human diseases. It has been reported that fullerenes and some of their derivatives–carboxyfullerenes–exhibits a strong free radical scavenging capacity. The permeation of C₆₀-fullerene and its amphiphilic derivatives–C₃-tris-malonic-C₆₀-fullerene (C₃) and D₃-tris-malonyl-C₆₀-fullerene (D₃)–through a lipid bilayer mimicking the eukaryotic cell membrane was studied using molecular dynamics (MD) simulations. The free energy profiles along the normal to the bilayer composed of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) for C₆₀, C₃ and D₃ were calculated. We found that C₆₀ molecules alone or in clusters spontaneously translocate to the hydrophobic core of the membrane and stay inside the bilayer during the whole period of simulation time. The incorporation of cluster of fullerenes inside the bilayer changes properties of the bilayer and leads to its deformation. In simulations of the tris-malonic fullerenes we discovered that both isomers, C₃ and D₃, adsorb at the surface of the bilayer but only C₃ tends to be buried in the area of the lipid headgroups forming hydrophobic contacts with the lipid tails. We hypothesize that such position has implications for ROS scavenging mechanism in the specific cell compartments.</p></div

C₃ stereoisomer of tris-malonic fullerene.

<p>A. Free energy profile of the process of the C₃ penetration into the model eukaryotic membrane. B. Intermediate orientation of the C₃ molecule adsorbed to the membrane with its solvent shell retained. C. Stable conformation (corresponding to the global energy minimum of the free energy profile) of C₃ adsorbed to the membrane and established hydrophobic contact with the lipid tails region.</p