Simulation of Forces between Humid Amorphous Silica Surfaces: A Comparison of Empirical Atomistic Force Fields.
Ontology highlight
ABSTRACT:
SUBMITTER: Leroch S
PROVIDER: S-EPMC3558025 | biostudies-other | 2012 Dec
REPOSITORIES: biostudies-other
Similar Datasets
Project description:All-atom empirical molecular mechanics protein force fields, which have been developed to represent the energetics of peptide folding behavior in aqueous solution, have not been parameterized for protein interactions with solid material surfaces. As a result, their applicability for representing the adsorption behavior of proteins with functionalized material surfaces should not be assumed. To address this issue, we conducted replica-exchange molecular dynamics simulations of the adsorption behavior of structured peptides to functionalized surfaces using three protein force fields that are widely used for the simulation of peptide adsorption behavior: CHARMM22, AMBER94, and OPLS-AA. Simulation results for peptide structure both in solution and when adsorbed to the surfaces were compared to experimental results for similar peptide-surface systems to provide a means of evaluating and comparing the performance of these three force fields for this type of application. Substantial differences in both solution and adsorbed peptide conformations were found amongst these three force fields, with the CHARMM22 force field found to most closely match experimental results.
| S-EPMC4243546 | biostudies-literature
Project description:The dynamics of polymer chains in the polymer/solid interphase region have been a point of debate in recent years. Its understanding is the first step towards the description and the prediction of the properties of a wide family of commercially used polymeric-based nanostructured materials. Here, we present a detailed investigation of the conformational and dynamical features of unentangled and mildly entangled cis-1,4-polybutadiene melts in the vicinity of amorphous silica surface via atomistic simulations. Accounting for the roughness of the surface, we analyze the properties of the polymer chains as a function of their distance from the silica slab, their conformations and the chain molecular weight. Unlike the case of perfectly flat and homogeneous surfaces, the monomeric translational motion parallel to the surface was affected by the presence of the silica slab up to distances comparable with the extension of the density fluctuations. In addition, the intramolecular dynamical heterogeneities in adsorbed chains were revealed by linking the conformations and the structure of the adsorbed chains with their dynamical properties. Strong dynamical heterogeneities within the adsorbed layer are found, with the chains possessing longer sequences of adsorbed segments ("trains") exhibiting slower dynamics than the adsorbed chains with short ones. Our results suggest that, apart from the density-dynamics correlation, the configurational entropy plays an important role in the dynamical response of the polymers confined between the silica slabs.
| S-EPMC8401719 | biostudies-literature
Project description:Hydration forces are ubiquitous in nature and technology. Yet, the characterization of interfacial hydration structures and their dependence on the nature of the substrate and the presence of ions have remained challenging and controversial. We present a systematic study using dynamic Atomic Force Microscopy of hydration forces on mica surfaces and amorphous silica surfaces in aqueous electrolytes containing chloride salts of various alkali and earth alkaline cations of variable concentrations at pH values between 3 and 9. Our measurements with ultra-sharp AFM tips demonstrate the presence of both oscillatory and monotonically decaying hydration forces of very similar strength on both atomically smooth mica and amorphous silica surfaces with a roughness comparable to the size of a water molecule. The characteristic range of the forces is approximately 1 nm, independent of the fluid composition. Force oscillations are consistent with the size of water molecules for all conditions investigated. Weakly hydrated Cs+ ions are the only exception: they disrupt the oscillatory hydration structure and induce attractive monotonic hydration forces. On silica, force oscillations are also smeared out if the size of the AFM tip exceeds the characteristic lateral scale of the surface roughness. The observation of attractive monotonic hydration forces for asymmetric systems suggests opportunities to probe water polarization.
| S-EPMC10568262 | biostudies-literature
Project description:The balance between hydration and Derjaguin-Landau-Verwey-Overbeek (DLVO) forces at solid-liquid interfaces controls many processes, such as colloidal stability, wetting, electrochemistry, biomolecular self-assembly, and ion adsorption. Yet, the origin of molecular scale hydration forces and their relation to the surface charge density that controls the continuum scale electrostatic forces is poorly understood. We argue that these two types of forces are largely independent of each other. To support this hypothesis, we performed atomic force microscopy experiments using intermediate-sized tips that enable the simultaneous detection of DLVO and molecular scale oscillatory hydration forces at the interface between composite gibbsite:silica-aqueous electrolyte interfaces. We extract surface charge densities from forces measured at tip-sample separations of 1.5 nm and beyond using DLVO theory in combination with charge regulation boundary conditions for various pH values and salt concentrations. We simultaneously observe both colloidal scale DLVO forces and oscillatory hydration forces for an individual crystalline gibbsite particle and the underlying amorphous silica substrate for all fluid compositions investigated. While the diffuse layer charge varies with pH as expected, the oscillatory hydration forces are found to be largely independent of pH and salt concentration, supporting our hypothesis that both forces indeed have a very different origin. Oscillatory hydration forces are found to be distinctly more pronounced on gibbsite than on silica. We rationalize this observation based on the distribution of hydroxyl groups available for H bonding on the two distinct surfaces.
| S-EPMC8793142 | biostudies-literature
Project description:Computer simulation increasingly complements experimental efforts to describe nanoscale structure formation. Molecular mechanics simulations and related computational methods fundamentally rely on the accuracy of classical atomistic force fields for the evaluation of inter- and intramolecular energies. One indispensable component of such force fields, in particular for large organic molecules, is the accuracy of molecule-specific dihedral potentials which are the key determinants of molecular flexibility. We show in this work that non-local correlations of dihedral potentials play a decisive role in the description of the total molecular energy-an effect which is neglected in most state-of-the-art dihedral force fields. We furthermore present an efficient machine learning approach to compute intramolecular conformational energies. We demonstrate with the example of α-NPD, a molecule frequently used in organic electronics, that this approach outperforms traditional force fields by decreasing the mean absolute deviations by one order of magnitude to values smaller than 0.37 kcal/mol (16.0 meV) per dihedral angle.
| S-EPMC5803249 | biostudies-literature
Project description:Predicting the structure and dynamics of RNA molecules still proves challenging because of the relative scarcity of experimental RNA structures on which to train models and the very sensitive nature of RNA towards its environment. In the last decade, several atomistic force fields specifically designed for RNA have been proposed and are commonly used for simulations. However, it is not necessarily clear which force field is the most suitable for a given RNA molecule. In this contribution, we propose the use of the computational energy landscape framework to explore the energy landscape of RNA systems as it can bring complementary information to the more standard approaches of enhanced sampling simulations based on molecular dynamics. We apply the EL framework to the study of a small RNA pseudoknot, the Aquifex aeolicus tmRNA pseudoknot PK1, and we compare the results of five different RNA force fields currently available in the AMBER simulation software, in implicit solvent. With this computational approach, we can not only compare the predicted 'native' states for the different force fields, but the method enables us to study metastable states as well. As a result, our comparison not only looks at structural features of low energy folded structures, but provides insight into folding pathways and higher energy excited states, opening to the possibility of assessing the validity of force fields also based on kinetics and experiments providing information on metastable and unfolded states.Supplementary informationThe online version contains supplementary material available at 10.1007/s12551-024-01202-9.
| S-EPMC11297004 | biostudies-literature
Project description:Plasmonic nanoparticles in the quantum regime exhibit characteristic optical properties that cannot be described by classical theories. Time-dependent density functional theory (TDDFT) is rising as a versatile tool for study on such systems, but its application has been limited to very small clusters due to rapidly growing computational costs. We propose an atomistic dipole-interaction-model for quantum plasmon simulations as a practical alternative. Namely the atomic dipole approximation represents induced dipoles with atomic polarizabilities obtained from TDDFT without empirical parameters. It showed very good agreement with TDDFT for plasmonic spectra of small silver clusters at much lower computational cost, though it is not appropriate for molecular-like excitations. It could also reproduce the plasmonic band shift experimentally observed in sub-10 nm silver particles.
| S-EPMC5693991 | biostudies-literature
Project description:The high concentration of macromolecules in the crowded cellular interior influences different thermodynamic and kinetic properties of proteins, including their structural stabilities, intermolecular binding affinities and enzymatic rates. Moreover, various structural biology methods, such as NMR or different spectroscopies, typically involve samples with relatively high protein concentration. Due to large sampling requirements, however, the accuracy of classical molecular dynamics (MD) simulations in capturing protein behavior at high concentration still remains largely untested. Here, we use explicit-solvent MD simulations and a total of 6.4 µs of simulated time to study wild-type (folded) and oxidatively damaged (unfolded) forms of villin headpiece at 6 mM and 9.2 mM protein concentration. We first perform an exhaustive set of simulations with multiple protein molecules in the simulation box using GROMOS 45a3 and 54a7 force fields together with different types of electrostatics treatment and solution ionic strengths. Surprisingly, the two villin headpiece variants exhibit similar aggregation behavior, despite the fact that their estimated aggregation propensities markedly differ. Importantly, regardless of the simulation protocol applied, wild-type villin headpiece consistently aggregates even under conditions at which it is experimentally known to be soluble. We demonstrate that aggregation is accompanied by a large decrease in the total potential energy, with not only hydrophobic, but also polar residues and backbone contributing substantially. The same effect is directly observed for two other major atomistic force fields (AMBER99SB-ILDN and CHARMM22-CMAP) as well as indirectly shown for additional two (AMBER94, OPLS-AAL), and is possibly due to a general overestimation of the potential energy of protein-protein interactions at the expense of water-water and water-protein interactions. Overall, our results suggest that current MD force fields may distort the picture of protein behavior in biologically relevant crowded environments.
| S-EPMC4031056 | biostudies-literature
Project description:Silicon (Si) is suggested to be an important/essential nutrient for bone and connective tissue health. Silicon-substituted hydroxyapatite (Si-HA) has silicate ions incorporated into its lattice structure and was developed to improve attachment to bone and increase new bone formation. Here we investigated the direct adsorption of silicate species onto an HA coated surface as a cost effective method of incorporating silicon on to HA surfaces for improved implant osseointegration, and determined changes in surface characteristics and osteoblast cell adhesion. Plasma-sprayed HA-coated stainless steel discs were incubated in silica dispersions of different concentrations (0-42 mM Si), at neutral pH for 12 h. Adsorbed Si was confirmed by XPS analysis and quantified by ICP-OES analysis following release from the HA surface. Changes in surface characteristics were determined by AFM and measurement of surface wettability. Osteoblast cell adhesion was determined by vinculin plaque staining. Maximum Si adsorption to the HA coated disc occurred after incubation in the 6 mM silica dispersion and decreased progressively with higher silica concentrations, while no adsorption was observed with dispersions below 6 mM Si. Comparison of the Si dispersions that produced the highest and lowest Si adsorption to the HA surface, by TEM-based analysis, revealed an abundance of small amorphous nanosilica species (NSP) of ~1.5 nm in diameter in the 6 mM Si dispersion, with much fewer and larger NSP in the 42 mM Si dispersions. 29Si-NMR confirmed that the NSPs in the 6 mM silica dispersion were polymeric and similar in composition to the larger NSPs in the 42 mM Si dispersion, suggesting that the latter were aggregates of the former. Amorphous NSP adsorbed from the 6 mM dispersion on to a HA-coated disc surface increased the surface's water contact angle by 53°, whereas that adsorbed from the 42 mM dispersion decreased the contact angle by 18°, indicating increased and decreased hydrophobicity, respectively. AFM showed an increase in surface roughness of the 6 mM Si treated surface, which correlated well with an increase in number of vinculin plaques. These findings suggest that NSP of the right size (relative to charge) adsorb readily to the HA surface, changing the surface characteristics and, thus, improving osteoblast cell adhesion. This treatment provides a simple way to modify plasma-coated HA surfaces that may enable improved osseointegration of bone implants.
| S-EPMC4749379 | biostudies-literature
Project description:The behavior of biomolecules in crowded environments remains largely unknown due to the accuracy of simulation models and the limited experimental data for comparison. Here we chose a small crowder of tetraethylene glycol (PEG-4) to investigate the self-crowding of PEG-4 solutions and molecular crowding effects on the structure and diffusion of lysozyme at varied concentrations from dilute water to pure PEG-4 liquid. Two Amber-like force fields of Amber14SB and a99SB-disp were examined with TIP3P (fast diffusivity and low viscosity) and a99SB-disp (slow diffusivity and high viscosity) water models, respectively. Compared to the Amber14SB protein simulations, the a99SB-disp model yields more coordinated water and less PEG-4 molecules, less intramolecular hydrogen bonds (HBs), more protein-water HBs, and less protein-PEG HBs as well as stronger interactions and more hydrophilic and less hydrophobic contacts with solvent molecules. The a99SB-disp model offers comparable protein-solvent interactions in concentrated PEG-4 solutions to that in pure water. The PEG-4 crowding leads to a slow-down in the diffusivity of water, PEG-4, and protein, and the decline in the diffusion from atomistic simulations is close to or faster than the hard sphere model that neglects attractive interactions. Despite these differences, the overall structure of lysozyme appears to be maintained well at different PEG-4 concentrations for both force fields, except a slightly large deviation at 370 K at low concentrations with the a99SB-disp model. This is mainly attributed to the strong intramolecular interactions of the protein in the Amber14SB force field and to the large viscosity of the a99SB-disp water model. The results indicate that the protein force fields and the viscosity of crowder solutions affect the simulation of biomolecules under crowding conditions.
| S-EPMC9000840 | biostudies-literature