Unveiling clinically significant PPARγ mutations for thiazolidinedione treatment responsiveness through atomistic simulations

In Type 2 diabetes, increased insulin sensitivity is induced by thiazolidinedione activation of the peroxisome proliferator-activated receptor gamma (PPARγ). Recent data indicate a relationship between SNPs in PPARγ and poor drug response. Therefore, understanding the pathogenic consequences of mutations in PPARγ-mediated protein-drug interactions will be prima-facie for establishing personalized medicine. The PPARG gene has 197 missense SNPs, 22 of which were determined to be both deleterious and destabilizing, employing in silico approaches. Molecular docking analysis suggested that the mutation influenced the binding energy of at least seven of the variants. The mutant R316H was identified as the most damaging and deleterious from the observed results. For a better understanding of the dynamic variation upon mutation at the atomic level, molecular dynamics simulations of the wild-type and R316H mutant PPARγ structure were performed. The analysis indicates that the mutation increased protein structural compactness while decreasing flexibility. The reduced dynamics in the mutant structure was further validated by principal component analysis. This mechanistic evaluation of the PPARγ protein variants provides insight into the relationship between genetic variation and interindividual variability of drug responsiveness and will facilitate the future studies for the development of tailored treatment regime for precision medicine.

Simulation toolkits at the molecular scale for trans-scale thermal signaling

Thermogenesis is a physiological activity of releasing heat that originates from intracellular biochemical reactions. Recent experimental studies discovered that externally applied heat changes intracellular signaling locally, resulting in global changes in cell morphology and signaling. Therefore, we hypothesize an inevitable contribution of thermogenesis in modulating biological system functions throughout the spatial scales from molecules to individual organisms. One key issue examining the hypothesis, namely, the "trans-scale thermal signaling," resides at the molecular scale on the amount of heat released via individual reactions and by which mechanism the heat is employed for cellular function operations. This review introduces atomistic simulation tool kits for studying the mechanisms of thermal signaling processes at the molecular scale that even state-of-the-art experimental methodologies of today are hardly accessible. We consider biological processes and biomolecules as potential heat sources in cells, such as ATP/GTP hydrolysis and multiple biopolymer complex formation and disassembly. Microscopic heat release could be related to mesoscopic processes via thermal conductivity and thermal conductance. Additionally, theoretical simulations to estimate these thermal properties in biological membranes and proteins are introduced. Finally, we envisage the future direction of this research field.

Design of short peptides and peptide amphiphiles as collagen mimics and an investigation of their interactions with collagen using molecular dynamics simulations and docking studies

Short peptide sequences and bolaamphiphiles derived from natural proteins are gaining importance due to their ability to form unique nanoscale architectures for a variety of biological applications. In this work, we have designed six short peptides (triplet or monomeric forms) and two peptide bolaamphiphiles that either incorporate the bioactive collagen motif (Gly-X-Y) or sequences where Gly, Pro, or hydroxyproline (Hyp) are replaced by Ala or His. For the bolaamphiphiles, a malate moiety was used as the aliphatic linker for connecting His with Hyp to create collagen mimics. Stability of the assemblies was assessed through molecular dynamics simulations and results indicated that (Pro-Ala-His)₃ and (Ala-His-Hyp)₃ formed the most stable structures, while the amphiphiles and the monomers showed some disintegration over the course of the 200 ns simulation, though most regained structural integrity and formed fibrillar structures, and micelles by the end of the simulation, likely due to the formation of more thermodynamically stable conformations. Multiple replica simulations (REMD) were also conducted where the sequences were simulated at different temperatures. Our results showed excellent convergence in most cases compared to constant temperature molecular dynamics simulation. Furthermore, molecular docking and MD simulations of the sequences bound to collagen triple helix structure revealed that several of the sequences had a high binding affinity and formed stable complexes, particularly (Pro-Ala-His)₃ and (Ala-His-Hyp)₃. Thus, we have designed new hybrid-peptide-based sequences which may be developed for potential applications as biomaterials for tissue engineering or drug delivery.

Using molecular simulation to understand the skin barrier

Skin's effectiveness as a barrier to permeation of water and other chemicals rests almost entirely in the outermost layer of the epidermis, the stratum corneum (SC), which consists of layers of corneocytes surrounded by highly organized lipid lamellae. As the only continuous path through the SC, transdermal permeation necessarily involves diffusion through these lipid layers. The role of the SC as a protective barrier is supported by its exceptional lipid composition consisting of ceramides (CERs), cholesterol (CHOL), and free fatty acids (FFAs) and the complete absence of phospholipids, which are present in most biological membranes. Molecular simulation, which provides molecular level detail of lipid configurations that can be connected with barrier function, has become a popular tool for studying SC lipid systems. We review this ever-increasing body of literature with the goals of (1) enabling the experimental skin community to understand, interpret and use the information generated from the simulations, (2) providing simulation experts with a solid background in the chemistry of SC lipids including the composition, structure and organization, and barrier function, and (3) presenting a state of the art picture of the field of SC lipid simulations, highlighting the difficulties and best practices for studying these systems, to encourage the generation of robust reproducible studies in the future. This review describes molecular simulation methodology and then critically examines results derived from simulations using atomistic and then coarse-grained models.

Defects, diffusion and dopants in Li8SnO6

Octalithium tin (IV) oxide (Li₈SnO₆) is an important electrode material considered for lithium ion batteries (LIBs) because of its high lithium content. We employed atomistic simulations to examine the intrinsic defects, diffusion of Li-ions together with their migration energies and solution of potential dopants in Li₈SnO₆. The most thermodynamically favourable intrinsic defect is the Li Frenkel which increases the concentration of Li vacancies needed for the vacancy mediated diffusion of Li-ions in Li₈SnO₆. The calculated activation energy of migration of Li-ions (0.21eV) shows that the Li-ion conductivity in this material can be very fast. Promising isovalent dopants on the Li and Sn sites are Na and Ti, respectively. Doping of Ga on the Sn site can facilitate the formation of Li interstitials as well as oxygen vacancies in Li₈SnO₆. While the concentration of Li interstitials can enhance the capacity of this material, oxygen vacancies together with Li interstitials can lead to the loss of Li₂O in Li₈SnO₆.

Combining scattering analysis and atomistic simulation of wood-water interactions

Molecular-scale interactions between water and cellulose microfibril bundles in plant cell walls are not fully understood, despite their crucial role for many applications of plant biomass. Recent advances in X-ray and neutron scattering analysis allow more accurate interpretation of experimental data from wood cell walls. At the same time, microfibril bundles including hemicelluloses and water can be modelled at atomistic resolution. Computing scattering patterns from atomistic models enables a new, complementary approach to decipher some of the most fundamental questions at this level of the hierarchical cell wall structure. This article introduces studies related to moisture behavior of wood with small/wide-angle X-ray/neutron scattering and atomistic simulations, recent attempts to combine these two approaches, and perspectives and open questions for future research using this powerful combination. Finally, we discuss the opportunities of the combined method in relation to applications of lignocellulosic materials.

How to minimize dye-induced perturbations while studying biomembrane structure and dynamics: PEG linkers as a rational alternative

Organic dye-tagged lipid analogs are essential for many fluorescence-based investigations of complex membrane structures, especially when using advanced microscopy approaches. However, lipid analogs may interfere with membrane structure and dynamics, and it is not obvious that the properties of lipid analogs would match those of non-labeled host lipids. In this work, we bridged atomistic simulations with super-resolution imaging experiments and biomimetic membranes to assess the performance of commonly used sphingomyelin-based lipid analogs. The objective was to compare, on equal footing, the relative strengths and weaknesses of acyl chain labeling, headgroup labeling, and labeling based on poly-ethyl-glycol (PEG) linkers in determining biomembrane properties. We observed that the most appropriate strategy to minimize dye-induced membrane perturbations and to allow consideration of Brownian-like diffusion in liquid-ordered membrane environments is to decouple the dye from a membrane by a PEG linker attached to a lipid headgroup. Yet, while the use of PEG linkers may sound a rational and even an obvious approach to explore membrane dynamics, the results also suggest that the dyes exploiting PEG linkers interfere with molecular interactions and their dynamics. Overall, the results highlight the great care needed when using fluorescent lipid analogs, in particular accurate controls.

Simulation of kinematic Kikuchi diffraction patterns from atomistic structures

One of the limitations of atomistic simulations is that many of the computational tools used to extract structural information from atomic trajectories provide metrics that are not directly compatible with experiments for validation. In this work, to bridge between simulation and experiment, a method is presented to produce simulated Kikuchi diffraction patterns using data from atomistic simulations, without requiring a priori specification of the crystal structure or defect periodicity. The Kikuchi pattern simulation is based on the kinematic theory of diffraction, with Kikuchi line intensities computed via a discrete structure factor calculation. Reciprocal lattice points are mapped to Kikuchi lines using a geometric projection of the reciprocal space data. This method is validated using single crystal atomistic models, and the novelty of this approach is emphasized by simulating kinematic Kikuchi diffraction patterns from an atomistic model containing a nanoscale dislocation loop. Deviations in kinematic Kikuchi line intensities are explained considering the displacement field of the dislocation loop, as is done in diffraction contrast theory.

Influence of Elastic Stiffness and Surface Adhesion on Bouncing of Nanoparticles

Granular collisions are characterized by a threshold velocity, separating the low-velocity regime of grain sticking from the high-velocity regime of grain bouncing: the bouncing velocity, v _b . This parameter is particularly important for nanograins and has applications for instance in astrophysics where it enters the description of collisional dust aggregation. Analytic estimates are based on the macroscopic Johnson-Kendall-Roberts (JKR) theory, which predicts the dependence of v _b on the radius, elastic stiffness, and surface adhesion of grains. Here, we perform atomistic simulations with model potentials that allow us to test these dependencies for nanograin collisions. Our results not only show that JKR describes the dependence on materials parameters qualitatively well, but also point at considerable quantitative deviations. These are the most pronounced for small adhesion, where elastic stiffness does not influence the value of the bouncing velocity.

Mechanical Deformation Mechanisms and Properties of Prion Fibrils Probed by Atomistic Simulations

Prion fibrils, which are a hallmark for neurodegenerative diseases, have recently been found to exhibit the structural diversity that governs disease pathology. Despite our recent finding concerning the role of the disease-specific structure of prion fibrils in determining their elastic properties, the mechanical deformation mechanisms and fracture properties of prion fibrils depending on their structures have not been fully characterized. In this work, we have studied the tensile deformation mechanisms of prion and non-prion amyloid fibrils by using steered molecular dynamics simulations. Our simulation results show that the elastic modulus of prion fibril, which is formed based on left-handed β-helical structure, is larger than that of non-prion fibril constructed based on right-handed β-helix. However, the mechanical toughness of prion fibril is found to be less than that of non-prion fibril, which indicates that infectious prion fibril is more fragile than non-infectious (non-prion) fibril. Our study sheds light on the role of the helical structure of amyloid fibrils, which is related to prion infectivity, in determining their mechanical deformation mechanisms and properties.

Membrane pore formation in atomistic and coarse-grained simulations

Biological cells and their organelles are protected by ultra thin membranes. These membranes accomplish a broad variety of important tasks like separating the cell content from the outer environment, they are the site for cell-cell interactions and many enzymatic reactions, and control the in- and efflux of metabolites. For certain physiological functions e.g. in the fusion of membranes and also in a number of biotechnological applications like gene transfection the membrane integrity needs to be compromised to allow for instance for the exchange of polar molecules across the membrane barrier. Mechanisms enabling the transport of molecules across the membrane involve membrane proteins that form specific pores or act as transporters, but also so-called lipid pores induced by external fields, stress, or peptides. Recent progress in the simulation field enabled to closely mimic pore formation as supposed to occur in vivo or in vitro. Here, we review different simulation-based approaches in the study of membrane pores with a focus on lipid pore properties such as their size and energetics, poration mechanisms based on the application of external fields, charge imbalances, or surface tension, and on pores that are induced by small molecules, peptides, and lipids. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Recent progress in atomistic simulation of electrical current DNA sequencing

We review recent advances in the DNA sequencing method based on measurements of transverse electrical currents. Device configurations proposed in the literature are classified according to whether the molecular fingerprints appear as the major (Mode I) or perturbing (Mode II) current signals. Scanning tunneling microscope and tunneling electrode gap configurations belong to the former category, while the nanochannels with or without an embedded nanopore belong to the latter. The molecular sensing mechanisms of Modes I and II roughly correspond to the electron tunneling and electrochemical gating, respectively. Special emphasis will be given on the computer simulation studies, which have been playing a critical role in the initiation and development of the field. We also highlight low-dimensional nanomaterials such as carbon nanotubes, graphene, and graphene nanoribbons that allow the novel Mode II approach. Finally, several issues in previous computational studies are discussed, which points to future research directions toward more reliable simulation of electrical current DNA sequencing devices.