Chromatin condensation is driven by the energetically favourable interaction between nucleosome core particles (NCPs). The close NCP-NCP contact, stacking, is a primary structural element of all ...condensed states of chromatin in vitro and in vivo. However, the molecular structure of stacked nucleosomes as well as the nature of the interactions involved in its formation have not yet been systematically studied. Here we undertake an investigation of both the structural and physico-chemical features of NCP structure and the NCP-NCP stacking. We introduce an "NCP-centred" set of parameters (NCP-NCP distance, shift, rise, tilt, and others) that allows numerical characterisation of the mutual positions of the NCPs in the stacking and in any other structures formed by the NCP. NCP stacking in more than 140 published NCP crystal structures were analysed. In addition, coarse grained (CG) MD simulations modelling NCP condensation was carried out. The CG model takes into account details of the nucleosome structure and adequately describes the long range electrostatic forces as well as excluded volume effects acting in chromatin. The CG simulations showed good agreement with experimental data and revealed the importance of the H2A and H4 N-terminal tail bridging and screening as well as tail-tail correlations in the stacked nucleosomes.
Piezoelectric semiconductors have emerged as redox catalysts, and challenges include effective conversion of mechanical energy to piezoelectric polarization and achieving high catalytic activity. The ...catalytic activity can be enhanced by simultaneous irradiation of ultrasound and light, but the existing piezoelectric semiconductors have trouble absorbing visible light. A piezoelectric catalyst is designed and tested for the generation of hydrogen peroxide (H2O2). It is based on Nb‐doped tetragonal BaTiO3 (BaTiO3:Nb) and is sensitized by carbon quantum dots (CDs). The photosensitizer injects electrons into the conduction band of the semiconductor, while the piezoelectric polarization directed electrons to the semiconductor surface, allowing for a high‐rate generation of H2O2. The piezoelectric polarization field restricts the recombination of photoinduced electron–hole pairs. A production rate of 1360 µmol gcatalyst−1 h−1 of H2O2 is achieved under visible light and ultrasound co‐irradiation. Individual piezo‐ and photocatalysis yielded lower production rates. Furthermore, the CDs enhance the piezocatalytic activity of the BaTiO3:Nb. It is noted that moderating the piezoelectricity of BaTiO3:Nb via microstructure modulation influences the piezophotocatalytic activity. This work shows a new methodology for synthesizing H2O2 by using visible light and mechanical energy.
Carbon quantum dots sensitized Nb‐doped tetragonal BaTiO3 nanopiezoelectrics are considerably effective in simultaneously utilizing visible light and vibration energy, thus achieving high hydrogen peroxide yields from ethanol and water suspensions. This is because the sensitized piezoelectrics‐mediated polarization field accelerates the migration of photoinduced charge carriers.
Rapid development of computer power during the last decade has made molecular simulations of lipid bilayers feasible for many research groups, which, together with the growing general interest in ...investigations of these very important biological systems has lead to tremendous increase of the number of research on the computational modeling of lipid bilayers. In this review, we give account of the recent progress in computer simulations of lipid bilayers covering mainly the period of the last 5 years, and covering several selected subjects: development of the force fields for lipid bilayer simulations, studies of the role of lipid unsaturation, the effect of cholesterol and other inclusions on properties of the bilayer, and use of coarse-grained models.
We review recent progress in computer simulations of lipid bilayers, covering questions of force field development, simulations of bilayers with unsaturated lipids, cholesterol and other inclusions, and use of coarse-grained models
Free energy calculations are vital for our understanding of biological processes on an atomistic scale and can offer insight to various mechanisms. However, in some cases, degrees of freedom (DOFs) ...orthogonal to the reaction coordinate have high energy barriers and/or long equilibration times, which prohibit proper sampling. Here we identify these orthogonal DOFs when studying the transfer of a solute from water to a model membrane. Important DOFs are identified in bulk liquids of different dielectric nature with metadynamics simulations and are used as reaction coordinates for the translocation process, resulting in two- and three-dimensional space of reaction coordinates. The results are in good agreement with experiments and elucidate the pitfalls of using one-dimensional reaction coordinates. The calculations performed here offer the most detailed free energy landscape of solutes embedded in lipid bilayers to date and show that free energy calculations can be used to study complex membrane translocation phenomena.
To provide insight into the molecular mechanisms of local anesthetic action, we have carried out an extensive investigation of two amide type local anesthetics, lidocaine and articaine in both ...charged and uncharged forms, interacting with DMPC lipid membrane. We have applied both standard molecular dynamics simulations and metadynamics simulations to provide a detailed description of the free energy landscape of anesthetics embedded in the lipid bilayer. The global minimum of the free energy surface (equilibrium position of anesthetics in the lipid membrane) occurred around 1nm of the bilayer center. The uncharged anesthetics show more affinity to bind to this region compared to the charged drugs. The binding free energy of uncharged lidocaine in the membrane (−30.3kJ/mol) is higher than uncharged articaine (−24.0kJ/mol), which is in good agreement with higher lipid solubility of lidocaine relative to the articaine. The octanol/water partition coefficient of uncharged drugs was also investigated using expanded ensemble simulations. In addition, complementary standard MD simulations were carried out to study the partitioning behavior of multiple anesthetics inside the lipid bilayer. The results obtained here are in line with previously reported simulations and suggest that the different forms of anesthetics induce different structural modifications in the lipid bilayer, which can provide new insights into their complex membrane translocation phenomena.
Display omitted
•Molecular dynamics and metadynamics simulations of local anesthetics lidocaine and articaine in a lipid bilayer have been carried out.•Localization and orientation of both charged and neutral forms of the drug molecules in the lipid membrane were investigated.•Binding free energies and partitioning properties of the drug molecules were determined.
Cardiolipin is a non-bilayer phospholipid with a unique dimeric structure. It localizes to negative curvature regions in bacteria and is believed to stabilize respiratory chain complexes in the ...highly curved mitochondrial membrane. Cardiolipin's localization mechanism remains unresolved, because important aspects such as the structural basis and strength for lipid curvature preferences are difficult to determine, partly due to the lack of efficient simulation methods. Here, we report a computational approach to study curvature preferences of cardiolipin by simulated membrane buckling and quantitative modeling. We combine coarse-grained molecular dynamics with simulated buckling to determine the curvature preferences in three-component bilayer membranes with varying concentrations of cardiolipin, and extract curvature-dependent concentrations and lipid acyl chain order parameter profiles. Cardiolipin shows a strong preference for negative curvatures, with a highly asymmetric chain order parameter profile. The concentration profiles are consistent with an elastic model for lipid curvature sensing that relates lipid segregation to local curvature via the material constants of the bilayers. These computations constitute new steps to unravel the molecular mechanism by which cardiolipin senses curvature in lipid membranes, and the method can be generalized to other lipids and membrane components as well.
The understanding of interactions between nanomaterials and biological molecules is of primary importance for biomedical applications of nanomaterials, as well as for the evaluation of their possible ...toxic effects. Here, we carried out extensive molecular dynamics simulations of the adsorption properties of about 30 small molecules representing biomolecular fragments at ZnS surfaces in aqueous media. We computed adsorption free energies and potentials of mean force of amino acid side chain analogs, lipids, and sugar fragments to ZnS (110) crystal surface and to a spherical ZnS nanoparticle. Furthermore, we investigated the effect of poly-methylmethacrylate (PMMA) coating on the adsorption preferences of biomolecules to ZnS. We found that only a few anionic molecules: aspartic and glutamic acids side chains, as well as the anionic form of cysteine show significant binding to pristine ZnS surface, while other molecules show weak or no binding. Spherical ZnS nanoparticles show stronger binding of these molecules due to binding at the edges between different surface facets. Coating of ZnS by PMMA changes binding preferences drastically: the molecules that adsorb to a pristine ZnS surface do not adsorb on PMMA-coated surfaces, while some others, particularly hydrophobic or aromatic amino-acids, show high binding affinity due to binding to the coating. We investigate further the hydration properties of the ZnS surface and relate them to the binding preferences of biomolecules.
Binary mixtures of ethanol and phospholipids DOPC and DOPE have been investigated in a composition range relevant for topical drug delivery applications. This was done using a combined computer ...simulation and experimental approach where molecular dynamics simulations of ethanol-lipid mixtures with different compositions were performed. Several key properties including diffusion coefficients, longitudinal relaxation times, and shear viscosity were computed. In addition, diffusion coefficients, viscosities and NMR longitudinal relaxation times were measured experimentally for comparison and in order to validate the results from simulation. Diffusion coefficients and relaxation times obtained from simulations are in good agreement with results from NMR and computed viscosities are in reasonable agreement with viscometry experiments indicating that the simulations provide a realistic description of the ethanol-phospholipid mixtures. Structural changes in the simulated systems were investigated using an analysis based on radial distribution functions. This showed that the structure of ethanol-DOPC mixtures remains essentially unchanged in the investigated concentration range while ethanol-DOPE mixtures undergo structural rearrangements with the tendency for forming small aggregates on the 100 ns time scale consisting of less than 10 lipids. Although our simulations and experiments indicate that no larger aggregates form, they also show that DOPE has stronger aggregation tendency than DOPC. This highlights the importance of the character of the lipid headgroup for lipid aggregation in ethanol and gives new insights into phase equilibrium, dynamics and rheology that could be valuable for the development of advanced topical drug delivery formulations.
Simulation insight and interpretation of results from NMR and viscometry experiments in lipid-ethanol mixtures.
A series of 19 hydroxylated polybrominated diphenyl ethers (OH-PBDEs) have been studied using density functional theory (DFT) and molecular dynamics simulations with the purpose of investigating ...eventual correlations between their physicochemical properties and toxic action. Dissociation constants (pK
), solvation free energies and octanol-water partition coefficients (log P) have been computed. Additionally, metadynamics simulations of OH-PBDEs passing through a lipid bilayer have been carried out for four OH-PBDE species. No correlations between computed pK
values and toxicity data have been found. Medium correlations were found between partition coefficients and the ability of OH-PBDEs to alter membrane potential in cell cultures, which is attributed to higher uptake of molecules with larger log P parameters. It was also demonstrated that in lipid bilayers, OH-PBDE molecules differ in their orientational distributions and can adopt different conformations which can affect the uptake of these molecules and influence the pathways of their toxic action.