Evaluating Charge Equilibration Methods To Generate Electrostatic Fields in Nanoporous Materials.
Ontology highlight
ABSTRACT:
SUBMITTER: Ongari D
PROVIDER: S-EPMC6328974 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:We present results of molecular dynamics simulations of a model DPPC-water monolayer using charge equilibration (CHEQ) force fields, which explicitly account for electronic polarization in a classical treatment of intermolecular interactions. The surface pressure, determined as the difference between the monolayer and pure water surface tensions at 323 K, is predicted to be 22.92 ±1.29 dyne/cm, just slightly below the broad range of experimental values reported for this system. The surface tension for the DPPC-water monolayer is predicted to be 42.35 ±1.16 dyne/cm, in close agreement with the experimentally determined value of 40.9 dyne/cm. This surface tension is also consistent with the value obtained from DPPC monolayer simulations using state-of-the-art nonpolarizable force fields. The current results of simulations predict a monolayer-water potential difference relative to the pure water-air interface of 0.64 ±0.02 Volts, an improved prediction compared to the fixed-charge CHARMM27 force field, yet still overestimating the experimental range of 0.3 to 0.45 Volts. As the charge equilibration model is a purely charge-based model for polarization, the current results suggest that explicitly modeled polarization effects can offer improvements in describing interfacial electrostatics in such systems.
| S-EPMC3488352 | biostudies-literature
Project description:All-atom molecular dynamics simulations have been applied in the recent past to explore the free energetics underlying ion transport processes in biological ion channels. Roux and co-workers, Kuyucak and co-workers, Busath and co-workers, and others have performed rather elegant and extended time scale molecular dynamics simulations using current state-of-the-art fixed-charge (nonpolarizable) force fields to calculate the potential of mean force defining the equilibrium flux of ions through prototypical channels such as gramicidin A. An inescapable conclusion of such studies has been the gross overestimation of the equilibrium free energy barrier, generally predicted to be from 10 to 20 kcal/mol depending on the force field and simulation protocol used in the calculation; this translates to an underestimation of experimentally measurable single channel conductances by several orders of magnitude. Next-generation polarizable force fields have been suggested as possible alternatives for more quantitative predictions of the underlying free energy surface in such systems. (1) Presently, we consider ion permeation energetics in the gramicidin A channel using a novel polarizable force field. Our results predict a peak barrier height of 6 kcal/mol relative to the channel entrance; this is significantly lower than the uncorrected value of 12 kcal/mol for nonpolarizable force fields such as GROMOS and CHARMM27 which do not account for electronic polarization. These results provide promising initial indications substantiating the long-conjectured importance of polarization effects in describing ion-protein interactions in narrow biological channels.
| S-EPMC2818143 | biostudies-literature
Project description:Polarizable force fields for lipid and solvent environments are used for molecular dynamics simulations of a fully hydrated dipalmitoylphosphatidylcholine (DPPC) bilayer and gramicidin A (gA) dimer embedded in a dimyristoylphosphatidylcholine (DMPC) bilayer. The lipid bilayer is modeled using the CHARMM charge equilibration (CHEQ) polarizable force field for lipids and the TIP4P-FQ force field to represent solvent. For the DPPC bilayer system, results are compared to the same system simulated using the nonpolarizable CHARMM27r (C27r) force field and TIP3P water. Calculated atomic and electron density profiles, head group orientations as measured by the phosphorus-nitrogen vector orientation, and deuterium order parameters are found to be consistent with previous simulations and with experiment. The CHEQ model exhibits greater water penetration into the bilayer interior, as demonstrated by the potential of mean force calculated from the water density profile. This is a result of the variation of the water molecular dipole from 2.55 D in the bulk to 1.88 D in the interior. We discuss this finding in the context of previous studies (both simulation and experiment) that have investigated the extent of penetration of water into DPPC bilayers. We also discuss the effects of including explicit polarization on the water dipole moment variation as a function of distance from the bilayer. We show distributions of atomic charges over the course of the simulation since the CHEQ model allows the charges to fluctuate. We have calculated the interfacial dipole potential, which the CHEQ model predicts to be 0.95 V compared to 0.86 V as predicted by the C27r model. We also discuss dielectric permittivity profiles and the differences arising between the two models. We obtain bulk values of 72.77 for the CHEQ model (TIP4P-FQ water) and 91.22 for C27r (TIP3P), and values approaching unity in the membrane interior. Finally, we present results of simulations of gA embedded in a DMPC bilayer using the CHEQ model and discuss structural properties.
| S-EPMC2746983 | biostudies-literature
Project description:In this work, we have studied the pH-dependent surface charge nature of nanoporous graphene. This has been investigated by membrane potential and by streaming current measurements, both with varying pH. We observed a lowering of the membrane potential with decreasing pH for a fixed concentration gradient of potassium chloride (KCl) in the Donnan dominated regime. Interestingly, the potential reverses its sign close to pH 4. The fitted value of effective fixed ion concentration (C̅ R) in the membrane also follows the same trend. The streaming current measurements show a similar trend with sign reversal around pH 4.2. The zeta potential data from the streaming current measurement is further analyzed using a 1-pK model. The model is used to determine a representative pK (acid-base equilibrium constant) of 4.2 for the surface of these perforated graphene membranes. In addition, we have also theoretically investigated the effect of the PET support in our membrane potential measurement using numerical simulations. Our results indicate that the concentration drop inside the PET support can be a major contributor (up to 85%) for a significant deviation of the membrane potential from the ideal Nernst potential.
| S-EPMC7859941 | biostudies-literature
Project description:Over the last few years it has been understood that the interface between living cells and the underlying materials can be a powerful tool to manipulate cell functions. In this study, we explore the hypothesis that the electrical cell/material interface can regulate the differentiation of cancer stem-like cells (CSCs). Electrospun polymer fibres, either polyamide 66 or poly(lactic acid), with embedded graphene nanoplatelets (GnPs), have been fabricated as CSC scaffolds, providing both the 3D microenvironment and a suitable electrical environment favorable for CSCs adhesion, growth and differentiation. We have investigated the impact of these scaffolds on the morphological, immunostaining and electrophysiological properties of CSCs extracted from human glioblastoma multiform (GBM) tumor cell line. Our data provide evidence in favor of the ability of GnP-incorporating scaffolds to promote CSC differentiation to the glial phenotype. Numerical simulations support the hypothesis that the electrical interface promotes the hyperpolarization of the cell membrane potential, thus triggering the CSC differentiation. We propose that the electrical cell/material interface can regulate endogenous bioelectrical cues, through the membrane potential manipulation, resulting in the differentiation of CSCs. Material-induced differentiation of stem cells and particularly of CSCs, can open new horizons in tissue engineering and new approaches to cancer treatment, especially GBM.
| S-EPMC9765404 | biostudies-literature
Project description:A procedure to derive the electrostatic potential (ESP) for dynamic charge densities obtained from structure models or maximum-entropy densities is introduced. The ESP essentially is obtained by inverse Fourier transform of the dynamic structure factors of the total charge density corresponding to the independent atom model, the multipole model or maximum-entropy densities, employing dedicated software that will be part of the BayMEM software package. Our approach is also discussed with respect to the Ewald summation method. It is argued that a meaningful ESP can only be obtained if identical thermal smearing is applied to the nuclear (positive) and electronic (negative) parts of the dynamic charge densities. The method is applied to structure models of dl-serine at three different temperatures of 20, 100 and 298 K. The ESP at locations near the atomic nuclei exhibits a drastic reduction with increasing temperature, the largest difference between the ESP from the static charge density and the ESP of the dynamic charge density being at T = 20 K. These features demonstrate that zero-point vibrations are sufficient for changing the spiky nature of the ESP at the nuclei into finite values. On 0.5 e Å-3 isosurfaces of the electron densities (taken as the molecular surface relevant to intermolecular interactions), the dynamic ESP is surprisingly similar at all temperatures, while the static ESP of a single molecule has a slightly larger range and is shifted towards positive potential values.
| S-EPMC5713142 | biostudies-literature
Project description:In most applications of nanoporous materials the pore structure is as important as the chemical composition as a determinant of performance. For example, one can alter performance in applications like carbon capture or methane storage by orders of magnitude by only modifying the pore structure. For these applications it is therefore important to identify the optimal pore geometry and use this information to find similar materials. However, the mathematical language and tools to identify materials with similar pore structures, but different composition, has been lacking. We develop a pore recognition approach to quantify similarity of pore structures and classify them using topological data analysis. This allows us to identify materials with similar pore geometries, and to screen for materials that are similar to given top-performing structures. Using methane storage as a case study, we also show that materials can be divided into topologically distinct classes requiring different optimization strategies.
| S-EPMC5457500 | biostudies-literature
Project description:Electrospray ionization mass spectrometry is increasingly applied to study the structures and interactions of membrane protein complexes. However, the charging mechanism is complicated by the presence of detergent micelles during ionization. Here, we show that the final charge of membrane proteins can be predicted by their molecular weight when released from the non-charge reducing saccharide detergents. Our data indicate that PEG detergents lower the charge depending on the number of detergent molecules in the surrounding micelle, whereas fos-choline detergents may additionally participate in ion-ion reactions after desolvation. The supercharging reagent sulfolane, on the other hand, has no discernible effect on the charge of detergent-free membrane proteins. Taking our observations into the context of protein-detergent interactions in the gas phase, we propose a charge equilibration model for the generation of native-like membrane protein ions. During ionization of the protein-detergent complex, the ESI charges are distributed between detergent and protein according to proton affinity of the detergent, number of detergent molecules, and surface area of the protein. Charge equilibration influenced by detergents determines the final charge state of membrane proteins. This process likely contributes to maintaining a native-like fold after detergent release and can be harnessed to stabilize particularly labile membrane protein complexes in the gas phase.
| S-EPMC8972943 | biostudies-literature
Project description:Machine learning potentials (MLPs) have revolutionized the field of atomistic simulations by describing atomic interactions with the accuracy of electronic structure methods at a small fraction of the cost. Most current MLPs construct the energy of a system as a sum of atomic energies, which depend on information about the atomic environments provided in the form of predefined or learnable feature vectors. If, in addition, nonlocal phenomena like long-range charge transfer are important, fourth-generation MLPs need to be used, which include a charge equilibration (Qeq) step to take the global structure of the system into account. This Qeq can significantly increase the computational cost and thus can become a computational bottleneck for large systems. In this Article, we present a highly efficient formulation of Qeq that does not require the explicit computation of the Coulomb matrix elements, resulting in a quasi-linear scaling method. Moreover, our approach also allows for the efficient calculation of energy derivatives, which explicitly consider the global structure-dependence of the atomic charges as obtained from Qeq. Due to its generality, the method is not restricted to MLPs and can also be applied within a variety of other force fields.
| S-EPMC11360134 | biostudies-literature
Project description:Pluripotency can be induced from different cell types using different methods. Several studies have reported that iPSC retain an epigenetic memory from the tissue from which they were derived. However, the influence of the method used for reprogramming in the final pluripotent state has not been investigated yet. Therefore, we have taken advantage of germline stem cells (GSCs) that can be converted into germline pluripotent stem cells (gPSCs) using specific culture conditions or into induced pluripotent stem cells (iPSCs) after the forced expression of four transcription factors. Here, we compare the gene expression profile, imprinting methylation pattern and differentiation potential between GSC-derived iPSCs and gPSCs. Our analysis demonstrates that both reprogramming methods induce equivalent levels of pluripotency and that a donor GSC transcriptional memory is retained only at early steps after reprogramming. In summary, our results show that different reprogramming methods generate undistinguishable pluripotent states and exclude a method-dependent epigenetic memory. Cell Culture. All the cell lines used in this study were derived from OG2 mice. Unless otherwise noted, all cell culture reagents were purchased from Gibco or PAA. Neonatal GSCs derived from 10 days postpartum (d.p.p) murine testis and adult GSCs derived from 35 d.p.p murine testis were derived as previously described and cultured on irradiated C3H feeder cells in a medium composed of StemPro-34 SFM (including supplement), N2 supplement (1:100), 1 % heat-inactivated fetal calf serum (FCS), 5 mg ml-1 bovine serum albumin (BSA) fraction V, 100 U ml-1 penicillin, 100 µg ml-1 streptomycin, 2 mM Lâglutamine, 0.1 mM non-essential amino acids (NEAA), 0.1 mM βâmercaptoethanol, MEM vitamins (1:100), 6 mg ml-1 Dâ(+)âglucose (Sigma), 30 mg ml-1 pyruvic acid, 0.1 % DLâlactic acid (Sigma), 60 ng ml-1 progesterone (Sigma), 30 ng ml-1 βâestradiol (Sigma), 15 ng ml-1 mouse EGF (Peprotech), 10 ng ml-1 human bFGF (Peprotech), 10 ng ml-1 human GDNF (Peprotech), and 1000 U ml-1 mouse LIF (inâhouse preparation). For maintenance culture, GSCs were passaged at a 1:2âratio once a week after digestion with 0.05 % trypsin/EDTA. In addition, ESCs, iPSC and gPSC were cultured on irradiated C3H feeders in the following medium: Knockout DMEM, 20% Knockout Serum Replacement, penicillin, streptomycin, Lâglutamine, NEAA, βâmercaptoethanol, and LIF. Generation of pluripotent cell lines. iPSC were generated using tetracycline-inducible lentiviruses. For the lentiviral production, 293T cells were plated at 50 % confluency in 10-cm dishes one day prior to transfection. The next day, the lentiviral vector (3 µg) together with the packaging plasmids psPAX2 (2 µg) and pMD2.G (1 µg) (Addgene plasmids #12260 and #12259, respectively) were transfected using FuGENE 6 reagent (Roche) according to the manufacturerâs protocol. After 48 h, viral-containing supernatants were filtered (0.45 µm), concentrated around 30-fold by centrifugation (26,000 à g, 2 h, 4 °C) and stored at -80 °C. For the generation of iPSCs, a single-cell suspension of 50,000 GSCs was infected overnight with five tetracycline-inducible lentiviruses (5â20 µl concentrated supernatant each) encoding for Oct4, Sox2, Klf4, c-Myc, and the reverse tetracycline-transactivator M2-rtTA (Addgene plasmids #20323, #20326, #20322, #20324, and #20342) in the presence of 4 µg ml-1 polybrene (Sigma). One week after the infection, the expression of the viral transgenes was induced by adding 4 µg ml-1 doxycycline (Sigma) to the medium. MEFs-derived iPSCs were generated using the same protocol. gPSC were generated following a previous reported protocol. Briefly, GSCs were seeded at a density of approx. 1000 cells cm-2 on feeder cells and were not splitted until gPSC colonies were observed. Both iPSC and gPSC colonies were identified by Oct4-GFP expression, manually picked and expanded under ESC conditions. In vitro differentiation of pluripotent cells. Embryoid body (EB) formation was induced in suspension culture adopting the hanging drop method (600 cells in 20 µl) for 5 days in ESC medium without LIF. EBs were seeded on gelatine-coated plates and further cultured for 14 days in the same medium. Chimera formation. Aggregations of putative pluripotent cells with denuded 8-cell stage embryos were performed as previously reported. After overnight incubation at 37 °C and 5 % CO2, 11â14 aggregates were transferred into the uterine horns of 2.5 d.p.c pseudo-pregnant surrogate mothers (mated with vasectomized males). The chimeric pups were obtained and dissected at 14.5 d.p.c. Microarray gene expression analysis. Purified cRNA samples were prepared with the linear TotalPrep RNA Amplification Kit (Ambion) from 500 ng original RNA, involving synthesis of double-stranded cDNA by T7-promoter-linked oligo(dT) primers and 14 hours of in-vitro transcription incorporating biotin-labelled nucleotides. Hybridizations onto MouseRef 8 v2.0 Expression BeadChips (Illumina) was performed according to the manufacturerâs recommendations. Two replicates were hybridized from each sample. Chips were stained with streptavidin-Cy3 (GE Healthcare) and scanned using an iScan reader (Illumina). The bead intensities were mapped to gene information using BeadStudio 3.2 (Illumina). Background correction was performed using the Affymetrix robust multiarray analysis (RMA) background correction model. Variance stabilization was performed by log2 scaling, and gene expression normalization was calculated with the method implemented in the lumi package of R-Bioconductor. Samples were hybridised as biological replicates. 26 samples were analyzed. iPSC[GSCneonatal]#1: iPSC day 10 from GSC, clone 1, 2 biological rep iPSC[GSCneonatal]#2: iPSC day 10 from GSC, clone 2, 2 biological rep iPSC[GSCadult]#1: iPSC day 40 from GSC, clone 1, 2 biological rep iPSC[GSCadult]#2: iPSC day 40 from GSC, clone 2, 2 biological rep gPSC[GSCneonatal]#1: gPSC day 10 from GSC, clone 1, 2 biological rep gPSC[GSCneonatal]#2: gPSC day 10 from GSC, clone 2, 2 biological rep gPSC[GSCadult]#1: gPSC day 40 from GSC, clone 1, 2 biological rep gPSC[GSCadult]#2: gPSC day 40 from GSC, clone 2, 2 biological rep ESC: OG2 ESCs, 2 biological rep iPSC[MEF]: iPSC (Tet-4F) from OG2 MEFs, 2 biological rep GSC[neonatal]: GSC day 10, 2 biological rep GSC[adult]: GSC day 40, 2 biological rep MEF: Mouse embryonic fibroblasts, 2 biological rep
2015-12-17 | E-GEOD-42993 | biostudies-arrayexpress