Vibrational Analysis of Brucite Surfaces and the Development of an Improved Force Field for Molecular Simulation of Interfaces.
Ontology highlight
ABSTRACT:
SUBMITTER: Zeitler TR
PROVIDER: S-EPMC3993912 | biostudies-other | 2014 Apr
REPOSITORIES: biostudies-other
Similar Datasets
Project description:Polylactide is a biodegradable polymer that is widely used for biomedical applications, and it is a replacement for some petroleum based polymers in applications that range from packaging to carpeting. Efforts to characterize and further enhance polylactide based systems using molecular simulations have to this point been hindered by the lack of accurate atomistic models for the polymer. Thus, we present force field parameters specifically suited for molecular modeling of PLA. The model, which we refer to as PLAFF3, is based on a combination of the OPLS and CHARMM force fields, with modifications to bonded and nonbonded parameters. Dihedral angle parameters were adjusted to reproduce DFT data using newly developed CMAP dihedral cross terms, and the model was further adjusted to reproduce experimentally resolved crystal structure conformations, melt density, volume expansivity, and the glass transition temperature of PLA. We recommend the use of PLAFF3 in modeling PLA in its crystalline or amorphous states and have provided the necessary input files required for the publicly available molecular dynamics code GROMACS.
| S-EPMC3237685 | biostudies-literature
Project description:Propane dehydrogenation (PDH) is an on-purpose catalytic technology to produce propylene from propane that operates at high temperatures, 773-973 K. Several key industry players have been active in developing new catalysts and processes with improved carbon footprint and economics, where Pt-based catalysts have played a central role. The optimization of these catalytic systems through computational and atomistic simulations requires large-scale models that account for their reactivity and dynamic properties. To address this challenge, we developed a new reactive ReaxFF force field (2023-Pt/C/H) that enables large-scale simulations of PDH reactions catalyzed on Pt surfaces. The optimization of force-field parameters relies on a large training set of density functional theory (DFT) calculations of Pt-catalyzed PDH mechanism, including geometries, adsorption and relative energies of reaction intermediates, and key C-H and C-C bond-breaking/forming reaction steps on the Pt(111) surface. The internal validation supports the accuracy of the developed 2023-Pt/C/H force-field parameters, resulting in mean absolute errors (MAE) against DFT data of 14 and 12 kJ mol-1 for relative energies of intermediates and energy barriers, respectively. We demonstrated the applicability of the 2023-Pt/C/H force field with reactive molecular dynamics simulations of propane on different Pt surface topologies and temperatures. The simulations successfully model the formation of propene in the gas phase as well as competitive, unproductive reactions such as deep dehydrogenation and C-C bond cleavage that produce H, C1 and C2 adsorbed species responsible of catalytic deactivation of Pt surface. Results show the following reactivity order: Pt(111) < Pt(100) < Pt(211), and that for the stepped Pt(211) surface, propane activation occurs on low-coordinated Pt atoms at the steps. The measured selectivity as a function of surface topology follows the same trend as activity, the Pt(211) facet being the most selective. The 2023-Pt/C/H reactive force field can also describe the increase of reactivity with the temperature. From these simulations, we were able to estimate the Arrhenius activation energy, 73 kJ mol-1, whose value is close to those reported experimentally for PDH catalyzed by large, supported Pt nanoparticles . The newly developed 2023-Pt/C/H reactive force field can be used in subsequent investigations of different Pt topologies and of collective effects such as temperature, propane pressure, or H surface coverage.
| S-EPMC10895921 | biostudies-literature
Project description:The physics-based molecular force field (PMFF) was developed by integrating a set of potential energy functions in which each term in an intermolecular potential energy function is derived based on experimental values, such as the dipole moments, lattice energy, proton transfer energy, and X-ray crystal structures. The term "physics-based" is used to emphasize the idea that the experimental observables that are considered to be the most relevant to each term are used for the parameterization rather than parameterizing all observables together against the target value. PMFF uses MM3 intramolecular potential energy terms to describe intramolecular interactions and includes an implicit solvation model specifically developed for the PMFF. We evaluated the PMFF in three ways. We concluded that the PMFF provides reliable information based on the structure in a biological system and interprets the biological phenomena accurately by providing more accurate evidence of the biological phenomena.
| S-EPMC7217328 | biostudies-literature
Project description:BACKGROUND:After the discovery of membrane-reversible electroporation decades ago, the procedure has been used extensively in biology, biotechnology and medicine. The research on the basic mechanism has increasingly attracted attention. Although most research has focused on models that consider all atomic and molecular interactions and much atomic-level information can be obtained, the huge computational demand limits the models to simulations of only a few nanometers on the spatial scale and a few nanoseconds on the time scale. In order to more comprehensively study the reversible electroporation mechanism of phospholipid membrane on the nanoscale and at longer time intervals of up to 100 ns, we developed a dipalmitoylphosphatidylcholine (DPPC) phospholipid membrane model with the coarse-grained Martini force field. The model was tested by separately examining the morphology of the phospholipid membrane, the hydrophilic channel size, the distribution of the voltage potential on both sides of the membrane, and the movement of water molecules and ions during electroporation. RESULTS:The results showed that the process went through several stages: (1) the formation of the pore with defects originating on the surface. (2) The maintenance of the pore. The defects expanded to large pores and the size remains unchanged for several nanoseconds. (3) Pore healing stage due to self-assembly. Phospholipid membrane shrunk and the pore size decreased until completely closed. The pores were not circular in cross-section for most of the time and the potential difference across the membrane decreased dramatically after the pores formed, with almost no restoration of membrane integrity even when the pores started to close. CONCLUSIONS:The mechanism of the reversible electroporation process on the nanoscale level, including defects, expansion, stability, and pore closing stages on a longer time scale of up to 100 ns was demonstrated more comprehensively with the coarse-grained Martini force field, which took both the necessary molecular information and the calculation efficiency into account.
| S-EPMC6933919 | biostudies-literature
Project description:Molecular dynamics simulation of the oxygenase component of nitrobenzene dioxygenase (NBDO) system, a member of the naphthalene family of Rieske nonheme iron dioxygenases, has been carried out using the AMBER force field combined with a new set of parameters for the description of the mononuclear nonheme iron center and iron-sulfur Rieske cluster. Simulation results provide information on the structure and dynamics of nitrobenzene dioxygenase in an aqueous environment and shed light on specific interactions that occur in its catalytic center. The results suggest that the architecture of the active site is stabilized by key hydrogen bonds, and Asn258 positions the substrate for oxidation. Analysis of protein-water interactions reveal the presence of a network of solvent molecules at the entrance to the active site, which could be of potential catalytic importance.
| S-EPMC4059247 | biostudies-literature
Project description:Coarse-grained (CG) force fields have become promising tools for studies of protein behavior, but the balance of speed and accuracy is still a challenge in the research of protein coarse graining methodology. In this work, 20 CG beads have been designed based on the structures of amino acid residues, with which an amino acid can be represented by one or two beads, and a CG solvent model with five water molecules was adopted to ensure the consistence with the protein CG beads. The internal interactions in protein were classified according to the types of the interacting CG beads, and adequate potential functions were chosen and systematically parameterized to fit the energy distributions. The proposed CG force field has been tested on eight proteins, and each protein was simulated for 1000 ns. Even without any extra structure knowledge of the simulated proteins, the Cα root mean square deviations (RMSDs) with respect to their experimental structures are close to those of relatively short time all atom molecular dynamics simulations. However, our coarse grained force field will require further refinement to improve agreement with and persistence of native-like structures. In addition, the root mean square fluctuations (RMSFs) relative to the average structures derived from the simulations show that the conformational fluctuations of the proteins can be sampled.
| S-EPMC3509591 | biostudies-literature
Project description:Molecular simulation is a valuable and complementary tool that may assist with the interpretation of single-molecule Förster resonance energy transfer (FRET) experiments, if the energy function is of sufficiently high quality. Here we present force-field parameters for one of the most common pairs of chromophores used in experiments, AlexaFluor 488 and 594. From microsecond molecular-dynamics simulations, we are able to recover both experimentally determined equilibrium constants and association/dissociation rates of the chromophores with free tryptophan, as well as the decay of fluorescence anisotropy of a labeled protein. We find that it is particularly important to obtain a correct balance of solute-water interactions in the simulations in order to faithfully capture the experimental anisotropy decays, which provide a sensitive benchmark for fluorophore mobility. Lastly, by a combination of experiment and simulation, we address a potential complication in the interpretation of experiments on polyproline, used as a molecular ruler for FRET experiments, namely the potential association of one of the chromophores with the polyproline helix. Under conditions where simulations accurately capture the fluorescence anisotropy decay, we find at most a modest, transient population of conformations in which the chromophores associate with the polyproline. Explicit calculation of FRET transfer efficiencies for short polyprolines yields results in good agreement with experiment. These results illustrate the potential power of a combination of molecular simulation and experiment in quantifying biomolecular dynamics.
| S-EPMC4457477 | biostudies-literature
Project description:Pantetheine is ubiquitous in nature in various forms of pantetheine-containing ligands (PCLs), including coenzyme A and phosphopantetheine. Lack of scalable force field libraries for PCLs has hampered the computational studies of biological macromolecules containing PCLs. We describe here the development of the first generation Pantetheine Force Field (PFF) library that is compatible with Amber force fields; parameterized using Gasteiger, AM1-BCC, or RESP charging methods combined with gaff2 and ff14SB parameter sets. In addition, a "plug-and-play" strategy was employed to enable the systematic charging of computationally expensive molecules sharing common substructural motifs. The validation studies performed on the PFF library showed promising performance where molecular dynamics (MD) simulations results were compared with experimental data of three representative systems. The PFF library represents the first force field library capable of modeling systems containing PCLs in silico and will aid in various applications including protein engineering and drug discovery.
| S-EPMC8266206 | biostudies-literature
Project description:The hydration free energy, structure, and dynamics of the zinc divalent cation are studied using a polarizable force field in molecular dynamics simulations. Parameters for the Zn(2+) are derived from gas-phase ab initio calculation of Zn(2+)-water dimer. The Thole-based dipole polarization is adjusted based on the Constrained Space Orbital Variations (CSOV) calculation while the Symmetry Adapted Perturbation Theory (SAPT) approach is also discussed. The vdW parameters of Zn(2+) have been obtained by comparing the AMOEBA Zn(2+)-water dimerization energy with results from several theory levels and basis sets over a range of distances. Molecular dynamics simulations of Zn(2+) solvation in bulk water are subsequently performed with the polarizable force field. The calculated first-shell water coordination number, water residence time and free energy of hydration are consistent with experimental and previous theoretical values. The study is supplemented with extensive Reduced Variational Space (RVS) and Electron Localization Function (ELF) computations in order to unravel the nature of the bonding in Zn(2+)(H(2)O)(n) (n=1,6) complexes and to analyze the charge transfer contribution to the complexes. Results show that the importance of charge transfer decreases as the size of Zn-water cluster grows due to anticooperativity and to changes in the nature of the metal-ligand bonds. Induction could be dominated by polarization when the system approaches condensed-phase and the covelant effects are eliminated from the Zn(II)-water interaction. To construct an "effective" classical polarizable potential for Zn(2+) in bulk water, one should therefore avoid over-fitting to the ab initio charge transfer energy of Zn(2+)-water dimer. Indeed, in order to avoid overestimation of condensed-phase many-body effects, which is crucial to the transferability of polarizable molecular dynamics, charge transfer should not be included within the classical polarization contribution and should preferably be either incorporated in to the pairwise van der Waals contribution or treated explicitly.
| S-EPMC2992432 | biostudies-literature
Project description:Many photoinduced excited states' relaxation processes and chemical reactions occur at interfaces and surfaces, including charge transfer, energy transfer, proton transfer, proton-coupled electron transfer, configurational dynamics, conical intersections, etc. Of them, interactions of electronic and vibrational motions, namely, vibronic couplings, are the main determining factors for the relaxation processes or reaction pathways. However, time-resolved electronic-vibrational spectroscopy for interfaces and surfaces is lacking. Here we develop interface/surface-specific two-dimensional electronic-vibrational sum frequency generation spectroscopy (2D-EVSFG) for time-dependent vibronic coupling of excited states at interfaces and surfaces. We further demonstrate the fourth-order technique by investigating vibronic coupling, solvent correlation, and time evolution of the coupling for photoexcited interface-active molecules, crystal violet (CV), at the air/water interface as an example. The two vibronic absorption peaks for CV molecules at the interface from the 2D-EVSFG experiments were found to be more prominent than their counterparts in bulk from 2D-EV. Quantitative analysis of the vibronic peaks in 2D-EVSFG suggested that a non-Condon process participates in the photoexcitation of CV at the interface. We further reveal vibrational solvent coupling for the zeroth level on the electronic state with respect to that on the ground state, which is directly related to the magnitude of its change in solvent reorganization energy. The change in the solvent reorganization energy at the interface is much smaller than that in bulk methanol. Time-dependent center line slopes (CLSs) of 2D-EVSFG also showed that kinetic behaviors of CV at the air/water interface are significantly different from those in bulk methanol. Our ultrafast 2D-EVSFG experiments not only offer vibrational information on both excited states and the ground state as compared with the traditional doubly resonant sum frequency generation and electronic-vibrational coupling but also provide vibronic coupling, dynamical solvent effects, and time evolution of vibronic coupling at interfaces.
| S-EPMC10375875 | biostudies-literature