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

Title: Realistic Extension Algorithm via Covariance Hessian

Abstract

Coarse-graining of protein interactions provides a means of simulating large biological systems. Here, a coarse-graining method, REACH, is introduced, in which the force constants of a residue-scale elastic network model are calculated from the variance-covariance matrix obtained from atomistic molecular dynamics (MD) simulation. In test calculations, the C-atoms variance-covariance matrices are calculated from the ensembles of 1-ns atomistic MD trajectories in monomeric and dimeric myoglobin, and used to derive coarse-grained force constants for the local and nonbonded interactions. Construction of analytical model functions of the distance-dependence of the interresidue force constants allows rapid calculation of the REACH normal modes. The model force constants from monomeric and dimeric myoglobin are found to be similar in magnitude to each other. The MD intra- and intermolecular mean-square fluctuations and the vibrational density of states are well reproduced by the residue-scale REACH normal modes without requiring rescaling of the force constant parameters. The temperature-dependence of the myoglobin REACH force constants reveals that the dynamical transition in protein internal fluctuations arises principally from softening of the elasticity in the nonlocal interactions. The REACH method is found to be a reliable way of determining spatiotemporal protein motion without the need for expensive computations of long atomisticmore » MD simulations.« less

Authors:
 [1];  [2]
  1. University of Heidelberg
  2. ORNL
Publication Date:
Research Org.:
Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
Sponsoring Org.:
USDOE Laboratory Directed Research and Development (LDRD) Program
OSTI Identifier:
932204
DOE Contract Number:  
DE-AC05-00OR22725
Resource Type:
Journal Article
Resource Relation:
Journal Name: Biophysical Journal; Journal Volume: 93; Journal Issue: 10
Country of Publication:
United States
Language:
English
Subject:
59 BASIC BIOLOGICAL SCIENCES; ALGORITHMS; CONSTRUCTION; ELASTICITY; FLUCTUATIONS; MATRICES; MYOGLOBIN; PROTEINS; SIMULATION; TRAJECTORIES

Citation Formats

Moritsugu, K, and Smith, Jeremy C. Realistic Extension Algorithm via Covariance Hessian. United States: N. p., 2007. Web. doi:10.1529/biophysj.107.111898.
Moritsugu, K, & Smith, Jeremy C. Realistic Extension Algorithm via Covariance Hessian. United States. doi:10.1529/biophysj.107.111898.
Moritsugu, K, and Smith, Jeremy C. Mon . "Realistic Extension Algorithm via Covariance Hessian". United States. doi:10.1529/biophysj.107.111898.
@article{osti_932204,
title = {Realistic Extension Algorithm via Covariance Hessian},
author = {Moritsugu, K and Smith, Jeremy C},
abstractNote = {Coarse-graining of protein interactions provides a means of simulating large biological systems. Here, a coarse-graining method, REACH, is introduced, in which the force constants of a residue-scale elastic network model are calculated from the variance-covariance matrix obtained from atomistic molecular dynamics (MD) simulation. In test calculations, the C-atoms variance-covariance matrices are calculated from the ensembles of 1-ns atomistic MD trajectories in monomeric and dimeric myoglobin, and used to derive coarse-grained force constants for the local and nonbonded interactions. Construction of analytical model functions of the distance-dependence of the interresidue force constants allows rapid calculation of the REACH normal modes. The model force constants from monomeric and dimeric myoglobin are found to be similar in magnitude to each other. The MD intra- and intermolecular mean-square fluctuations and the vibrational density of states are well reproduced by the residue-scale REACH normal modes without requiring rescaling of the force constant parameters. The temperature-dependence of the myoglobin REACH force constants reveals that the dynamical transition in protein internal fluctuations arises principally from softening of the elasticity in the nonlocal interactions. The REACH method is found to be a reliable way of determining spatiotemporal protein motion without the need for expensive computations of long atomistic MD simulations.},
doi = {10.1529/biophysj.107.111898},
journal = {Biophysical Journal},
number = 10,
volume = 93,
place = {United States},
year = {Mon Jan 01 00:00:00 EST 2007},
month = {Mon Jan 01 00:00:00 EST 2007}
}