-
A Comparison of Classical Force-Fields for Molecular Dynamics Simulations of Lubricants
PubMed Central
Ewen, James P.; Gattinoni, Chiara; Thakkar, Foram M.; Morgan, Neal; Spikes, Hugh A.; Dini, Daniele
2016-01-01
For the successful development and application of lubricants, a full understanding of their complex nanoscale behavior under a wide range of external conditions is required, but this is difficult to obtain experimentally. Nonequilibrium molecular dynamics (NEMD) simulations can be used to yield unique insights into the atomic-scale structure and friction of lubricants and additives; however, the accuracy of the results depend on the chosen force-field. In this study, we demonstrate that the use of an accurate, all-atom force-field is critical in order to; (i) accurately predict important properties of long-chain, linear molecules; and (ii) reproduce experimental friction behavior of multi-component tribological systems. In particular, we focus on n-hexadecane, an important model lubricant with a wide range of industrial applications. Moreover, simulating conditions common in tribological systems, i.e., high temperatures and pressures (HTHP), allows the limits of the selected force-fields to be tested. In the first section, a large number of united-atom and all-atom force-fields are benchmarked in terms of their density and viscosity prediction accuracy of n-hexadecane using equilibrium molecular dynamics (EMD) simulations at ambient and HTHP conditions. Whilst united-atom force-fields accurately reproduce experimental density, the viscosity is significantly under-predicted compared to all-atom force-fields and experiments. Moreover, some all-atom force-fields yield elevated melting points, leading to significant overestimation of both the density and viscosity. In the second section, the most accurate united-atom and all-atom force-field are compared in confined NEMD simulations which probe the structure and friction of stearic acid adsorbed on iron oxide and separated by a thin layer of n-hexadecane. The united-atom force-field provides an accurate representation of the structure of the confined stearic acid film; however, friction coefficients are consistently under-predicted and the friction-coverage and friction-velocity behavior deviates from that observed using all-atom force-fields and experimentally. This has important implications regarding force-field selection for NEMD simulations of systems containing long-chain, linear molecules; specifically, it is recommended that accurate all-atom potentials, such as L-OPLS-AA, are employed. PMID:28773773
-
A Comparison of Classical Force-Fields for Molecular Dynamics Simulations of Lubricants.
PubMed
Ewen, James P; Gattinoni, Chiara; Thakkar, Foram M; Morgan, Neal; Spikes, Hugh A; Dini, Daniele
2016-08-02
For the successful development and application of lubricants, a full understanding of their complex nanoscale behavior under a wide range of external conditions is required, but this is difficult to obtain experimentally. Nonequilibrium molecular dynamics (NEMD) simulations can be used to yield unique insights into the atomic-scale structure and friction of lubricants and additives; however, the accuracy of the results depend on the chosen force-field. In this study, we demonstrate that the use of an accurate, all-atom force-field is critical in order to; (i) accurately predict important properties of long-chain, linear molecules; and (ii) reproduce experimental friction behavior of multi-component tribological systems. In particular, we focus on n -hexadecane, an important model lubricant with a wide range of industrial applications. Moreover, simulating conditions common in tribological systems, i.e., high temperatures and pressures (HTHP), allows the limits of the selected force-fields to be tested. In the first section, a large number of united-atom and all-atom force-fields are benchmarked in terms of their density and viscosity prediction accuracy of n -hexadecane using equilibrium molecular dynamics (EMD) simulations at ambient and HTHP conditions. Whilst united-atom force-fields accurately reproduce experimental density, the viscosity is significantly under-predicted compared to all-atom force-fields and experiments. Moreover, some all-atom force-fields yield elevated melting points, leading to significant overestimation of both the density and viscosity. In the second section, the most accurate united-atom and all-atom force-field are compared in confined NEMD simulations which probe the structure and friction of stearic acid adsorbed on iron oxide and separated by a thin layer of n -hexadecane. The united-atom force-field provides an accurate representation of the structure of the confined stearic acid film; however, friction coefficients are consistently under-predicted and the friction-coverage and friction-velocity behavior deviates from that observed using all-atom force-fields and experimentally. This has important implications regarding force-field selection for NEMD simulations of systems containing long-chain, linear molecules; specifically, it is recommended that accurate all-atom potentials, such as L-OPLS-AA, are employed.
-
Principal Component Analysis of Lipid Molecule Conformational Changes in Molecular Dynamics Simulations.
PubMed
Buslaev, Pavel; Gordeliy, Valentin; Grudinin, Sergei; Gushchin, Ivan
2016-03-08
Molecular dynamics simulations of lipid bilayers are ubiquitous nowadays. Usually, either global properties of the bilayer or some particular characteristics of each lipid molecule are evaluated in such simulations, but the structural properties of the molecules as a whole are rarely studied. Here, we show how a comprehensive quantitative description of conformational space and dynamics of a single lipid molecule can be achieved via the principal component analysis (PCA). We illustrate the approach by analyzing and comparing simulations of DOPC bilayers obtained using eight different force fields: all-atom generalized AMBER, CHARMM27, CHARMM36, Lipid14, and Slipids and united-atom Berger, GROMOS43A1-S3, and GROMOS54A7. Similarly to proteins, most of the structural variance of a lipid molecule can be described by only a few principal components. These major components are similar in different simulations, although there are notable distinctions between the older and newer force fields and between the all-atom and united-atom force fields. The DOPC molecules in the simulations generally equilibrate on the time scales of tens to hundreds of nanoseconds. The equilibration is the slowest in the GAFF simulation and the fastest in the Slipids simulation. Somewhat unexpectedly, the equilibration in the united-atom force fields is generally slower than in the all-atom force fields. Overall, there is a clear separation between the more variable previous generation force fields and significantly more similar new generation force fields (CHARMM36, Lipid14, Slipids). We expect that the presented approaches will be useful for quantitative analysis of conformations and dynamics of individual lipid molecules in other simulations of lipid bilayers.
-
Efficient molecular dynamics simulations with many-body potentials on graphics processing units
NASA Astrophysics Data System (ADS)
Fan, Zheyong; Chen, Wei; Vierimaa, Ville; Harju, Ari
2017-09-01
Graphics processing units have been extensively used to accelerate classical molecular dynamics simulations. However, there is much less progress on the acceleration of force evaluations for many-body potentials compared to pairwise ones. In the conventional force evaluation algorithm for many-body potentials, the force, virial stress, and heat current for a given atom are accumulated within different loops, which could result in write conflict between different threads in a CUDA kernel. In this work, we provide a new force evaluation algorithm, which is based on an explicit pairwise force expression for many-body potentials derived recently (Fan et al., 2015). In our algorithm, the force, virial stress, and heat current for a given atom can be accumulated within a single thread and is free of write conflicts. We discuss the formulations and algorithms and evaluate their performance. A new open-source code, GPUMD, is developed based on the proposed formulations. For the Tersoff many-body potential, the double precision performance of GPUMD using a Tesla K40 card is equivalent to that of the LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator) molecular dynamics code running with about 100 CPU cores (Intel Xeon CPU X5670 @ 2.93 GHz).
-
Coercion from the Air: The United States Use of Airpower to Influence End of Conflict Negotiations
DTIC Science & Technology
2017-05-25
shock of two atomic bombs drove the Japanese to surrender. In the months and days leading to the armistice in Korea, Far East Air Force (FEAF...mainland by the Twentieth Air Force and the shock of two atomic bombs drove the Japanese to surrender. In the months and days leading to the...increasing destruction brought upon the Japanese population and war-making capacity combined with the shock of two atomic bombs to drive the Japanese towards
-
Molecular dynamics simulations of polarizable DNA in crystal environment
NASA Astrophysics Data System (ADS)
Babin, Volodymyr; Baucom, Jason; Darden, Thomas A.; Sagui, Celeste
We have investigated the role of the electrostatic description and cell environment in molecular dynamics (MD) simulations of DNA. Multiple unrestrained MD simulations of the DNA duplex d(CCAACGTTGG)2 have been carried out using two different force fields: a traditional description based on atomic point charges and a polarizable force field. For the time scales probed, and given the ?right? distribution of divalent ions, the latter performs better than the nonpolarizable force field. In particular, by imposing the experimental unit cell environment, an initial configuration with ideal B-DNA duplexes in the unit cell acquires sequence-dependent features that very closely resemble the crystallographic ones. Simultaneously, the all-atom root-mean-square coordinates deviation (RMSD) with respect to the crystallographic structure is seen to decay. At later times, the polarizable force field is able to maintain this lower RMSD, while the nonpolarizable force field starts to drift away.
-
Insight into the Properties of Cardiolipin Containing Bilayers from Molecular Dynamics Simulations, Using a Hybrid All-Atom/United-Atom Force Field.
PubMed
Aguayo, Daniel; González-Nilo, Fernando D; Chipot, Christophe
2012-05-08
Simulation of three models of cardiolipin (CL) containing membranes using a new set of parameters for tetramyristoyl and tetraoleoyl CLs has been developed in the framework of the united-atom CHARMM27-UA and the all-atom CHARMM36 force fields with the aim of performing molecular dynamics (MD) simulations of cardiolipin-containing mixed-lipid membranes. The new parameters use a hybrid representation of all-atom head groups in conjunction with implicit-hydrogen united-atom (UA) to describe the oleoyl and myristoyl chains of the CLs, in lieu of the fully atomistic description, thereby allowing longer simulations to be undertaken. The physicochemical properties of the bilayers were determined and compared with previously reported data. Furthermore, using tetramyristoyl CL mixed with POPG and POPE lipids, a mitochondrial membrane was simulated. The results presented here show the different behavior of the bilayers as a result of the lipid composition, where the length of the acyl chain and the conformation of the headgroup can be associated with the mitochondrial membrane properties. The new hybrid CL parameters prove to be well suited for the simulation of the molecular structure of CL-containing bilayers and can be extended to other lipid bilayers composed of CLs with different acyl chains or alternate head groups.
-
Optimization of a Nucleic Acids united-RESidue 2-Point model (NARES-2P) with a maximum-likelihood approach
NASA Astrophysics Data System (ADS)
He, Yi; Liwo, Adam; Scheraga, Harold A.
2015-12-01
Coarse-grained models are useful tools to investigate the structural and thermodynamic properties of biomolecules. They are obtained by merging several atoms into one interaction site. Such simplified models try to capture as much as possible information of the original biomolecular system in all-atom representation but the resulting parameters of these coarse-grained force fields still need further optimization. In this paper, a force field optimization method, which is based on maximum-likelihood fitting of the simulated to the experimental conformational ensembles and least-squares fitting of the simulated to the experimental heat-capacity curves, is applied to optimize the Nucleic Acid united-RESidue 2-point (NARES-2P) model for coarse-grained simulations of nucleic acids recently developed in our laboratory. The optimized NARES-2P force field reproduces the structural and thermodynamic data of small DNA molecules much better than the original force field.
-
Quantifying the atomic-level mechanics of single long physisorbed molecular chains.
PubMed
Kawai, Shigeki; Koch, Matthias; Gnecco, Enrico; Sadeghi, Ali; Pawlak, Rémy; Glatzel, Thilo; Schwarz, Jutta; Goedecker, Stefan; Hecht, Stefan; Baratoff, Alexis; Grill, Leonhard; Meyer, Ernst
2014-03-18
Individual in situ polymerized fluorene chains 10-100 nm long linked by C-C bonds are pulled vertically from an Au(111) substrate by the tip of a low-temperature atomic force microscope. The conformation of the selected chains is imaged before and after manipulation using scanning tunneling microscopy. The measured force gradient shows strong and periodic variations that correspond to the step-by-step detachment of individual fluorene repeat units. These variations persist at constant intensity until the entire polymer is completely removed from the surface. Calculations based on an extended Frenkel-Kontorova model reproduce the periodicity and magnitude of these features and allow us to relate them to the detachment force and desorption energy of the repeat units. The adsorbed part of the polymer slides easily along the surface during the pulling process, leading to only small oscillations as a result of the high stiffness of the fluorenes and of their length mismatch with respect to the substrate surface structure. A significant lateral force also is caused by the sequential detachment of individual units. The gained insight into the molecule-surface interactions during sliding and pulling should aid the design of mechanoresponsive nanosystems and devices.
-
Introduction of steered molecular dynamics into UNRES coarse-grained simulations package.
PubMed
Sieradzan, Adam K; Jakubowski, Rafał
2017-03-30
In this article, an implementation of steered molecular dynamics (SMD) in coarse-grain UNited RESidue (UNRES) simulations package is presented. Two variants of SMD have been implemented: with a constant force and a constant velocity. The huge advantage of SMD implementation in the UNRES force field is that it allows to pull with the speed significantly lower than the accessible pulling speed in simulations with all-atom representation of a system, with respect to a reasonable computational time. Therefore, obtaining pulling speed closer to those which appear in the atomic force spectroscopy is possible. The newly implemented method has been tested for behavior in a microcanonical run to verify the influence of introduction of artificial constrains on keeping total energy of the system. Moreover, as time dependent artificial force was introduced, the thermostat behavior was tested. The new method was also tested via unfolding of the Fn3 domain of human contactin 1 protein and the I27 titin domain. Obtained results were compared with Gø-like force field, all-atom force field, and experimental results. © 2017 Wiley Periodicals, Inc. © 2017 Wiley Periodicals, Inc.
-
Probing atomic-scale friction on reconstructed surfaces of single-crystal semiconductors
NASA Astrophysics Data System (ADS)
Goryl, M.; Budzioch, J.; Krok, F.; Wojtaszek, M.; Kolmer, M.; Walczak, L.; Konior, J.; Gnecco, E.; Szymonski, M.
2012-02-01
Friction force microscopy (FFM) investigations have been performed on reconstructed (001) surfaces of InSb and Ge in an ultrahigh vacuum. On the c(8×2) reconstruction of InSb(001) atomic resolution is achieved under superlubric conditions, and the features observed in the lateral force images are precisely reproduced by numerical simulations, taking into account possible decorations of the probing tip. On the simultaneously acquired (1×3) reconstruction a significant disorder of the surface atoms is observed. If the loading force increases, friction becomes much larger on this reconstruction compared to the c(8×2) one. In FFM images acquired on the Ge(001)(2×1) characteristic substructures are resolved within the unit cells. In such a case, a strong dependence of the friction pattern on the scan direction is observed.
-
Computer simulation of concentrated solid solution strengthening
NASA Technical Reports Server (NTRS)
Kuo, C. T. K.; Arsenault, R. J.
1976-01-01
The interaction forces between a straight edge dislocation moving through a three-dimensional block containing a random array of solute atoms were determined. The yield stress at 0 K was obtained by determining the average maximum solute-dislocation interaction force that is encountered by edge dislocation, and an expression relating the yield stress to the length of the dislocation and the solute concentration is provided. The magnitude of the solid solution strengthening due to solute atoms can be determined directly from the numerical results, provided the dislocation line length that moves as a unit is specified.
-
Comparison of force fields on the basis of various model approaches--how to design the best model for the [CnMIM][NTf2] family of ionic liquids.
PubMed
Köddermann, Thorsten; Reith, Dirk; Ludwig, Ralf
2013-10-07
In this contribution, we present two new united-atom force fields (UA-FFs) for 1-alkyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide [C(n)MIM][NTf(2)] (n=1, 2, 4, 6, 8) ionic liquids (ILs). One is parametrized manually, and the other is developed with the gradient-based optimization workflow (GROW). By doing so, we wanted to perform a hard test to determine how researchers could benefit from semiautomated optimization procedures. As with our already published all-atom force field (AA-FF) for [C(n)MIM][NTf(2)] (T. Köddermann, D. Paschek, R. Ludwig, ChemPhysChem- 2007, 8, 2464), the new force fields were derived to fit experimental densities, self-diffusion coefficients, and NMR rotational correlation times for the IL cation and for water molecules dissolved in [C(2)MIM][NTf(2)]. In the manual force field, the alkyl chains of the cation and the CF3 groups of the anion were treated as united atoms. In the GROW force field, only the alkyl chains of the cation were united. All other parts of the structures of the ions remained unchanged to prevent any loss of physical information. Structural, dynamic, and thermodynamic properties such as viscosity, cation rotational correlation times, and heats of vaporization calculated with the new force fields were compared with values simulated with the previous AA-FF and the experimental data. All simulated properties were in excellent agreement with the experimental values. Altogether, the UA-FFs are slightly superior for speed-up reasons. The UA-FF speeds up the simulation by about 100 % and reduces the demanded disk space by about 78 %. More importantly, real time and efforts to generate force fields could be significantly reduced by utilizing GROW. The real time for the GROW parametrization in this work was 2 months. Manual parametrization, in contrast, may take up to 12 months, and this is, therefore, a significant increase in speed, though it is difficult to estimate the duration of manual parametrization. Copyright © 2013 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim.
-
Reconciling Structural and Thermodynamic Predictions Using All-Atom and Coarse-Grain Force Fields: The Case of Charged Oligo-Arginine Translocation into DMPC Bilayers
PubMed Central
2015-01-01
Using the translocation of short, charged cationic oligo-arginine peptides (mono-, di-, and triarginine) from bulk aqueous solution into model DMPC bilayers, we explore the question of the similarity of thermodynamic and structural predictions obtained from molecular dynamics simulations using all-atom and Martini coarse-grain force fields. Specifically, we estimate potentials of mean force associated with translocation using standard all-atom (CHARMM36 lipid) and polarizable and nonpolarizable Martini force fields, as well as a series of modified Martini-based parameter sets. We find that we are able to reproduce qualitative features of potentials of mean force of single amino acid side chain analogues into model bilayers. In particular, modifications of peptide–water and peptide–membrane interactions allow prediction of free energy minima at the bilayer–water interface as obtained with all-atom force fields. In the case of oligo-arginine peptides, the modified parameter sets predict interfacial free energy minima as well as free energy barriers in almost quantitative agreement with all-atom force field based simulations. Interfacial free energy minima predicted by a modified coarse-grained parameter set are −2.51, −4.28, and −5.42 for mono-, di-, and triarginine; corresponding values from all-atom simulations are −0.83, −3.33, and −3.29, respectively, all in units of kcal/mol. We found that a stronger interaction between oligo-arginine and the membrane components and a weaker interaction between oligo-arginine and water are crucial for producing such minima in PMFs using the polarizable CG model. The difference between bulk aqueous and bilayer center states predicted by the modified coarse-grain force field are 11.71, 14.14, and 16.53 kcal/mol, and those by the all-atom model are 6.94, 8.64, and 12.80 kcal/mol; those are of almost the same order of magnitude. Our simulations also demonstrate a remarkable similarity in the structural aspects of the ensemble of configurations generated using the all-atom and coarse-grain force fields. Both resolutions show that oligo-arginine peptides adopt preferential orientations as they translocate into the bilayer. The guiding theme centers on charged groups maintaining coordination with polar and charged bilayer components as well as local water. We also observe similar behaviors related with membrane deformations. PMID:25290376
-
Reconciling structural and thermodynamic predictions using all-atom and coarse-grain force fields: the case of charged oligo-arginine translocation into DMPC bilayers.
PubMed
Hu, Yuan; Sinha, Sudipta Kumar; Patel, Sandeep
2014-10-16
Using the translocation of short, charged cationic oligo-arginine peptides (mono-, di-, and triarginine) from bulk aqueous solution into model DMPC bilayers, we explore the question of the similarity of thermodynamic and structural predictions obtained from molecular dynamics simulations using all-atom and Martini coarse-grain force fields. Specifically, we estimate potentials of mean force associated with translocation using standard all-atom (CHARMM36 lipid) and polarizable and nonpolarizable Martini force fields, as well as a series of modified Martini-based parameter sets. We find that we are able to reproduce qualitative features of potentials of mean force of single amino acid side chain analogues into model bilayers. In particular, modifications of peptide-water and peptide-membrane interactions allow prediction of free energy minima at the bilayer-water interface as obtained with all-atom force fields. In the case of oligo-arginine peptides, the modified parameter sets predict interfacial free energy minima as well as free energy barriers in almost quantitative agreement with all-atom force field based simulations. Interfacial free energy minima predicted by a modified coarse-grained parameter set are -2.51, -4.28, and -5.42 for mono-, di-, and triarginine; corresponding values from all-atom simulations are -0.83, -3.33, and -3.29, respectively, all in units of kcal/mol. We found that a stronger interaction between oligo-arginine and the membrane components and a weaker interaction between oligo-arginine and water are crucial for producing such minima in PMFs using the polarizable CG model. The difference between bulk aqueous and bilayer center states predicted by the modified coarse-grain force field are 11.71, 14.14, and 16.53 kcal/mol, and those by the all-atom model are 6.94, 8.64, and 12.80 kcal/mol; those are of almost the same order of magnitude. Our simulations also demonstrate a remarkable similarity in the structural aspects of the ensemble of configurations generated using the all-atom and coarse-grain force fields. Both resolutions show that oligo-arginine peptides adopt preferential orientations as they translocate into the bilayer. The guiding theme centers on charged groups maintaining coordination with polar and charged bilayer components as well as local water. We also observe similar behaviors related with membrane deformations.
-
Further optimization of a hybrid united-atom and coarse-grained force field for folding simulations: Improved backbone hydration and interactions between charged side chains
PubMed Central
Han, Wei; Schulten, Klaus
2012-01-01
PACE, a hybrid force field which couples united-atom protein models with coarse-grained (CG) solvent, has been further optimized, aiming to improve itse ciency for folding simulations. Backbone hydration parameters have been re-optimized based on hydration free energies of polyalanyl peptides through atomistic simulations. Also, atomistic partial charges from all-atom force fields were combined with PACE in order to provide a more realistic description of interactions between charged groups. Using replica exchange molecular dynamics (REMD), ab initio folding using the new PACE has been achieved for seven small proteins (16 – 23 residues) with different structural motifs. Experimental data about folded states, such as their stability at room temperature, melting point and NMR NOE constraints, were also well reproduced. Moreover, a systematic comparison of folding kinetics at room temperature has been made with experiments, through standard MD simulations, showing that the new PACE may speed up the actual folding kinetics 5-10 times. Together with the computational speedup benefited from coarse-graining, the force field provides opportunities to study folding mechanisms. In particular, we used the new PACE to fold a 73-residue protein, 3D, in multiple 10 – 30 μs simulations, to its native states (Cα RMSD ~ 0.34 nm). Our results suggest the potential applicability of the new PACE for the study of folding and dynamics of proteins. PMID:23204949
-
High-speed atomic force microscopy combined with inverted optical microscopy for studying cellular events
PubMed Central
Suzuki, Yuki; Sakai, Nobuaki; Yoshida, Aiko; Uekusa, Yoshitsugu; Yagi, Akira; Imaoka, Yuka; Ito, Shuichi; Karaki, Koichi; Takeyasu, Kunio
2013-01-01
A hybrid atomic force microscopy (AFM)-optical fluorescence microscopy is a powerful tool for investigating cellular morphologies and events. However, the slow data acquisition rates of the conventional AFM unit of the hybrid system limit the visualization of structural changes during cellular events. Therefore, high-speed AFM units equipped with an optical/fluorescence detection device have been a long-standing wish. Here we describe the implementation of high-speed AFM coupled with an optical fluorescence microscope. This was accomplished by developing a tip-scanning system, instead of a sample-scanning system, which operates on an inverted optical microscope. This novel device enabled the acquisition of high-speed AFM images of morphological changes in individual cells. Using this instrument, we conducted structural studies of living HeLa and 3T3 fibroblast cell surfaces. The improved time resolution allowed us to image dynamic cellular events. PMID:23823461
-
High-speed atomic force microscopy combined with inverted optical microscopy for studying cellular events.
PubMed
Suzuki, Yuki; Sakai, Nobuaki; Yoshida, Aiko; Uekusa, Yoshitsugu; Yagi, Akira; Imaoka, Yuka; Ito, Shuichi; Karaki, Koichi; Takeyasu, Kunio
2013-01-01
A hybrid atomic force microscopy (AFM)-optical fluorescence microscopy is a powerful tool for investigating cellular morphologies and events. However, the slow data acquisition rates of the conventional AFM unit of the hybrid system limit the visualization of structural changes during cellular events. Therefore, high-speed AFM units equipped with an optical/fluorescence detection device have been a long-standing wish. Here we describe the implementation of high-speed AFM coupled with an optical fluorescence microscope. This was accomplished by developing a tip-scanning system, instead of a sample-scanning system, which operates on an inverted optical microscope. This novel device enabled the acquisition of high-speed AFM images of morphological changes in individual cells. Using this instrument, we conducted structural studies of living HeLa and 3T3 fibroblast cell surfaces. The improved time resolution allowed us to image dynamic cellular events.
-
Behavior of P85 and P188 Poloxamer Molecules: Computer Simulations Using United Atom Force Field.
DOE PAGES
Goliaei, Ardeshir; Lau, Edmond Y.; Adhikari, Upendra; ...
2016-05-27
To study the interaction between poloxamer molecules and lipid bilayers using molecular dynamics simulation technique with the united atom resolution, we augmented the GROMOS force field to include poloxamers. We validated the force field by calculating the radii of gyration of two poloxamers, P85 and P188, solvated in water and by considering the poloxamer density distributions at the air/water interface. The emphasis of our simulations was on the study of the interaction between poloxamers and lipid bilayer. At the water/lipid bilayer interface, we observed that both poloxamers studied, P85 and P188, behaved like surfactants: the hydrophilic blocks of poloxamers becamemore » adsorbed at the polar interface, while their hydrophobic block penetrated the interface into the aliphatic tail region of the lipid bilayer. We also observed that when P85 and P188 poloxamers interacted with damaged membranes that contained pores, the hydrophobic blocks of copolymers penetrated into the membrane in the vicinity of the pore and compressed the membrane. Lastly, due to this compression, water molecules were evacuated from the pore.« less
-
Precise control of surface electrostatic forces on polymer brush layers with opposite charges for resistance to protein adsorption.
PubMed
Sakata, Sho; Inoue, Yuuki; Ishihara, Kazuhiko
2016-10-01
Various molecular interaction forces are generated during protein adsorption process on material surfaces. Thus, it is necessary to control them to suppress protein adsorption and the subsequent cell and tissue responses. A series of binary copolymer brush layers were prepared via surface-initiated atom transfer radical polymerization, by mixing the cationic monomer unit and anionic monomer unit randomly in various ratios. Surface characterization revealed that the constructed copolymer brush layers exhibited an uniform super-hydrophilic nature and different surface potentials. The strength of the electrostatic interaction forces operating on these mixed-charge copolymer brush surfaces was evaluated quantitatively using force-versus-distance (f-d) curve measurements by atomic force microscopy (AFM) and probes modified by negatively charged carboxyl groups or positively charged amino groups. The electrostatic interaction forces were determined based on the charge ratios of the copolymer brush layers. Notably, the surface containing equivalent cationic/anionic monomer units hardly interacted with both the charged groups. Furthermore, the protein adsorption force and the protein adsorption mass on these surfaces were examined by AFM f-d curve measurement and surface plasmon resonance measurement, respectively. To clarify the influence of the electrostatic interaction on the protein adsorption behavior on the surface, three kinds of proteins having negative, positive, and relatively neutral net charges under physiological conditions were used in this study. We quantitatively demonstrated that the amount of adsorbed proteins on the surfaces would have a strong correlation with the strength of surface-protein interaction forces, and that the strength of surface-protein interaction forces would be determined from the combination between the properties of the electrostatic interaction forces on the surfaces and the charge properties of the proteins. Especially, the copolymer brush surface composed of equivalent cationic/anionic monomer units exhibited no significant interaction forces, and dramatically suppressed the adsorption of proteins regardless of their charge properties. We conclude that the established methodology could elucidate relationship between the protein adsorption behavior and molecular interaction, especially the electrostatic interaction forces, and demonstrated that the suppression of the electrostatic interactions with the ionic functional groups would be important for the development of new polymeric biomaterials with a high repellency of protein adsorption. Copyright © 2016 Elsevier Ltd. All rights reserved.