Dataset Information

A comparison of multiscale methods for the analysis of molecular dynamics simulations.

ABSTRACT: Molecular dynamics (MD) is the only technique available for obtaining dynamic protein data at atomic spatial resolution and picosecond or finer temporal resolution. In recent years, the cost of computational resources has decreased exponentially while the number of known protein structures, many of which are not characterized biochemically, has increased rapidly. These events have led to an increase in the use of MD in biological research, both to examine phenomena that cannot be resolved experimentally and to generate hypotheses that direct further experimental research. In fact, several databases of MD simulations have arisen in recent years. MD simulations, and especially MD simulation databases, contain massive amounts of data, yet interesting phenomena often occur over very short time periods and on the scale of only a few atoms. Analysis of such data must balance these fine-detail events with the global picture they create. Here, we address the multiscale nature of the problem by comparing several MD analysis methods to show their strengths and weaknesses at various scales using the wild-type and R282W mutant forms of the DNA-binding domain of protein p53. By leveraging these techniques together, we are able to pinpoint fine-detail and big picture differences between the protein's variants. Our analyses indicate that the R282W mutation of p53 destabilizes the L1 loop and loosens the H2 helix conformation, but the loosened L1 loop can be rescued by residue H115, preventing the R282W mutation from completely destabilizing the protein or abolishing activity.

SUBMITTER: Benson NC

PROVIDER: S-EPMC3406285 | biostudies-literature | 2012 Jul

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

A comparison of multiscale methods for the analysis of molecular dynamics simulations.

Benson Noah C NC Daggett Valerie V

The journal of physical chemistry. B 20120424 29

Molecular dynamics (MD) is the only technique available for obtaining dynamic protein data at atomic spatial resolution and picosecond or finer temporal resolution. In recent years, the cost of computational resources has decreased exponentially while the number of known protein structures, many of which are not characterized biochemically, has increased rapidly. These events have led to an increase in the use of MD in biological research, both to examine phenomena that cannot be resolved experi ...[more]

PMID: 22494262

Similar Datasets

Project description:When photoactive molecules interact strongly with confined light modes in optical cavities, new hybrid light-matter states form. They are known as polaritons and correspond to coherent superpositions of excitations of the molecules and of the cavity photon. The polariton energies and thus potential energy surfaces are changed with respect to the bare molecules, such that polariton formation is considered a promising paradigm for controlling photochemical reactions. To effectively manipulate photochemistry with confined light, the molecules need to remain in the polaritonic state long enough for the reaction on the modified potential energy surface to take place. To understand what determines this lifetime, we have performed atomistic molecular dynamics simulations of room-temperature ensembles of rhodamine chromophores strongly coupled to a single confined light mode with a 15 fs lifetime. We investigated three popular experimental scenarios and followed the relaxation after optically pumping (i) the lower polariton, (ii) the upper polariton, or (iii) uncoupled molecular states. The results of the simulations suggest that the lifetimes of the optically accessible lower and upper polaritons are limited by (i) ultrafast photoemission due to the low cavity lifetime and (ii) reversible population transfer into the "dark" state manifold. Dark states are superpositions of molecular excitations but with much smaller contributions from the cavity photon, decreasing their emission rates and hence increasing their lifetimes. We find that population transfer between polaritonic modes and dark states is determined by the overlap between the polaritonic and molecular absorption spectra. Importantly, excitation can also be transferred "upward" from the lower polariton into the dark-state reservoir due to the broad absorption spectra of the chromophores, contrary to the common conception of these processes as a "one-way" relaxation from the dark states down to the lower polariton. Our results thus suggest that polaritonic chemistry relying on modified dynamics taking place within the lower polariton manifold requires cavities with sufficiently long lifetimes and, at the same time, strong light-matter coupling strengths to prevent the back-transfer of excitation into the dark states.

Project description: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.

Dataset Information

A comparison of multiscale methods for the analysis of molecular dynamics simulations.

Publications

A comparison of multiscale methods for the analysis of molecular dynamics simulations.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets