skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Determining Atomistic SAXS Models of Tri-Ubiquitin Chains from Bayesian Analysis of Accelerated Molecular Dynamics Simulations

Authors:
 [1];  [2];  [1];  [3]; ORCiD logo [4]; ORCiD logo [1]
  1. Department of Physics and Center for the Molecular Study of Condensed Soft Matter, Illinois Institute of Technology, Chicago, Illinois 60616, United States
  2. Department of Chemistry, University of Massachusetts−Amherst, Amherst, Massachusetts 01003, United States; Department of Chemistry, University of Wisconsin−Madison, Madison, Wisconsin 53706, United States
  3. Department of Chemistry, University of Wisconsin−Madison, Madison, Wisconsin 53706, United States
  4. Department of Chemistry, University of Massachusetts−Amherst, Amherst, Massachusetts 01003, United States; Department of Biochemistry and Molecular Biology, University of Massachusetts−Amherst, Amherst, Massachusetts 01003, United States
Publication Date:
Research Org.:
Argonne National Lab. (ANL), Argonne, IL (United States). Advanced Photon Source (APS)
Sponsoring Org.:
National Institutes of Health (NIH)
OSTI Identifier:
1395364
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Chemical Theory and Computation; Journal Volume: 13; Journal Issue: 6
Country of Publication:
United States
Language:
ENGLISH

Citation Formats

Bowerman, Samuel, Rana, Ambar S. J. B., Rice, Amy, Pham, Grace H., Strieter, Eric R., and Wereszczynski, Jeff. Determining Atomistic SAXS Models of Tri-Ubiquitin Chains from Bayesian Analysis of Accelerated Molecular Dynamics Simulations. United States: N. p., 2017. Web. doi:10.1021/acs.jctc.7b00059.
Bowerman, Samuel, Rana, Ambar S. J. B., Rice, Amy, Pham, Grace H., Strieter, Eric R., & Wereszczynski, Jeff. Determining Atomistic SAXS Models of Tri-Ubiquitin Chains from Bayesian Analysis of Accelerated Molecular Dynamics Simulations. United States. doi:10.1021/acs.jctc.7b00059.
Bowerman, Samuel, Rana, Ambar S. J. B., Rice, Amy, Pham, Grace H., Strieter, Eric R., and Wereszczynski, Jeff. 2017. "Determining Atomistic SAXS Models of Tri-Ubiquitin Chains from Bayesian Analysis of Accelerated Molecular Dynamics Simulations". United States. doi:10.1021/acs.jctc.7b00059.
@article{osti_1395364,
title = {Determining Atomistic SAXS Models of Tri-Ubiquitin Chains from Bayesian Analysis of Accelerated Molecular Dynamics Simulations},
author = {Bowerman, Samuel and Rana, Ambar S. J. B. and Rice, Amy and Pham, Grace H. and Strieter, Eric R. and Wereszczynski, Jeff},
abstractNote = {},
doi = {10.1021/acs.jctc.7b00059},
journal = {Journal of Chemical Theory and Computation},
number = 6,
volume = 13,
place = {United States},
year = 2017,
month = 5
}
  • Molecular dynamics (MD) simulations were performed on dense liquids of polyethylene chains of 24 and 66 united atom CH{sub 2} units. A series of models was studied ranging in atomistic detail from coarse-grained, freely-jointed, tangent site chains to realistic, overlapping site models subjected to bond angle restrictions and torsional potentials. These same models were also treated with the self-consistent, polymer reference interaction site model (PRISM) theory. The intramolecular and total structure factors, as well as, the intermolecular radial distribution functions g(r) and direct correlation functions C(r) were obtained from theory and simulation. Angular correlation functions were also simulation obtained frommore » the MD simulations. Comparisons between theory and reveal that PRISM theory works well for computing the intermolecular structure of coarse-grained chain models, but systematically underpredicts the extent of intermolecular packing as more atomistic details are introduced into the model. A consequence of g(r) having insufficient structure is that the theory yields an isothermal compressibility that progressively becomes larger, relative to the simulations, as overlapping the PRISM sites and angular restrictions are introduced into the model. We found that theory could be considerably improved by adding a tail function to C(r) beyond the effective hard core diameter. The range of this tail function was determined by requiring the theory to yield the correct compressibility.« less
  • The typical ionic liquid, 1-ethyl-3-methyl imidazolium nitrate (EMIM-NO{sub 3}), was examined by molecular dynamics simulations of an all-atomistic model to show the characteristics of networks of cages and/or bonds in the course of vitrification of this fragile glass-former. The system shows changes of dynamics at two characteristic temperatures, T{sub B} (or T{sub c}) and the glass transition temperature T{sub g}, found in other fragile glass forming liquids [K. L. Ngai and J. Habasaki, J. Chem. Phys. 141, 114502 (2014)]. On decreasing temperature, the number of neighboring cation-anion pairs, N{sub B}, within the first minimum of the pair correlation function, g(r){submore » min}, increases. On crossing T{sub B} (>T{sub g}), the system volume and diffusion coefficient both show changes in temperature dependence, and as usual at T{sub g}. The glass transition temperature, T{sub g}, is characterized by the saturation of the total number of “bonds,” N{sub B} and the corresponding decrease in degree of freedom, F = [(3N − 6) − N{sub B}], of the system consisting of N particles. Similar behavior holds for the other ion-ion pairs. Therefore, as an alternative, the dynamics of glass transition can be interpreted conceptually by rigidity percolation. Before saturation occurring at T{sub g}, the number of bonds shows a remarkable change at around T{sub B}. This temperature is associated with the disappearance of the loosely packed coordination polyhedra of anions around cation (or vice versa), related to the loss of geometrical freedom of the polyhedra, f{sub g}, of each coordination polyhedron, which can be defined by f{sub g} = [(3N{sub V} − 6) − N{sub b}]. Here, 3N{sub v} is the degree of freedom of N{sub V} vertices of the polyhedron, and N{sub b} is number of fictive bonds. The packing of polyhedra is characterized by the soft percolation of cages, which allows further changes with decreasing temperature. The power spectrum of displacement of the central ion in the cage is found to be correlated with the fluctuation of N{sub b} of cation-cation (or anion-anion) pairs in the polyhedron, although the effect from the coordination shells beyond the neighboring ions is not negligible.« less
  • We have determined the interfacial properties of tetrahydrofuran (THF) from direct simulation of the vapor-liquid interface. The molecules are modeled using six different molecular models, three of them based on the united-atom approach and the other three based on a coarse-grained (CG) approach. In the first case, THF is modeled using the transferable parameters potential functions approach proposed by Chandrasekhar and Jorgensen [J. Chem. Phys. 77, 5073 (1982)] and a new parametrization of the TraPPE force fields for cyclic alkanes and ethers [S. J. Keasler et al., J. Phys. Chem. B 115, 11234 (2012)]. In both cases, dispersive and coulombicmore » intermolecular interactions are explicitly taken into account. In the second case, THF is modeled as a single sphere, a diatomic molecule, and a ring formed from three Mie monomers according to the SAFT-γ Mie top-down approach [V. Papaioannou et al., J. Chem. Phys. 140, 054107 (2014)]. Simulations were performed in the molecular dynamics canonical ensemble and the vapor-liquid surface tension is evaluated from the normal and tangential components of the pressure tensor along the simulation box. In addition to the surface tension, we have also obtained density profiles, coexistence densities, critical temperature, density, and pressure, and interfacial thickness as functions of temperature, paying special attention to the comparison between the estimations obtained from different models and literature experimental data. The simulation results obtained from the three CG models as described by the SAFT-γ Mie approach are able to predict accurately the vapor-liquid phase envelope of THF, in excellent agreement with estimations obtained from TraPPE model and experimental data in the whole range of coexistence. However, Chandrasekhar and Jorgensen model presents significant deviations from experimental results. We also compare the predictions for surface tension as obtained from simulation results for all the models with experimental data. The three CG models predict reasonably well (but only qualitatively) the surface tension of THF, as a function of temperature, from the triple point to the critical temperature. On the other hand, only the TraPPE united-atoms models are able to predict accurately the experimental surface tension of the system in the whole temperature range.« less
  • A self-adaptive accelerated molecular dynamics method is developed to model infrequent atomic- scale events, especially those events that occur on a rugged free-energy surface. Key in the new development is the use of the total displacement of the system at a given temperature to construct a boost-potential, which is slowly increased to accelerate the dynamics. The temperature is slowly increased to accelerate the dynamics. By allowing the system to evolve from one steady-state con guration to another by overcoming the transition state, this self-evolving approach makes it possible to explore the coupled motion of species that migrate on vastly differentmore » time scales. The migrations of single vacancy (V) and small He-V clusters, and the growth of nano-sized He-V clusters in Fe for times in the order of seconds are studied by this new method. An interstitial- assisted mechanism is rst explored for the migration of a helium-rich He-V cluster, while a new two-component Ostwald ripening mechanism is suggested for He-V cluster growth.« less