Volume 9, Issue 3, March 2015
Index of content:
- Theoretical Methods and Algorithms
142(2015); http://dx.doi.org/10.1063/1.4913399View Description Hide Description
A new theory and an exact computer algorithm for calculating kinetics and thermodynamic properties of a particle system are described. The algorithm avoids trapping in metastable states, which are typical challenges for Molecular Dynamics (MD) simulations on rough energy landscapes. It is based on the division of the full space into Voronoi cells. Prior knowledge or coarse sampling of space points provides the centers of the Voronoi cells. Short time trajectories are computed between the boundaries of the cells that we call milestones and are used to determine fluxes at the milestones. The flux function, an essential component of the new theory, provides a complete description of the statistical mechanics of the system at the resolution of the milestones. We illustrate the accuracy and efficiency of the exact Milestoning approach by comparing numerical results obtained on a model system using exact Milestoning with the results of long trajectories and with a solution of the corresponding Fokker-Planck equation. The theory uses an equation that resembles the approximate Milestoning method that was introduced in 2004 [A. K. Faradjian and R. Elber, J. Chem. Phys. 120(23), 10880-10889 (2004)]. However, the current formulation is exact and is still significantly more efficient than straightforward MD simulations on the system studied.
142(2015); http://dx.doi.org/10.1063/1.4908531View Description Hide Description
We recently derived a Markov model for macromolecular ligand binding dynamics from few physical assumptions and showed that its stationary distribution is the grand canonical ensemble [J. W. R. Martini, M. Habeck, and M. Schlather, J. Math. Chem. 52, 665 (2014)]. The transition probabilities of the proposed Markov process define a particular Glauber dynamics and have some similarity to the Metropolis-Hastings algorithm. Here, we illustrate that this model is the stochastic analog of (pseudo) rate equations and the corresponding system of differential equations. Moreover, it can be viewed as a limiting case of general stochastic simulations of chemical kinetics. Thus, the model links stochastic and deterministic approaches as well as kinetics and equilibrium described by the grand canonical ensemble. We demonstrate that the family of transition matrices of our model, parameterized by temperature and ligand activity, generates ligand binding kinetics that respond to changes in these parameters in a qualitatively similar way as experimentally observed kinetics. In contrast, neither the Metropolis-Hastings algorithm nor the Glauber heat bath reflects changes in the external conditions correctly. Both converge rapidly to the stationary distribution, which is advantageous when the major interest is in the equilibrium state, but fail to describe the kinetics of ligand binding realistically. To simulate cellular processes that involve the reversible stochastic binding of multiple factors, our pseudo rate equation model should therefore be preferred to the Metropolis-Hastings algorithm and the Glauber heat bath, if the stationary distribution is not of only interest.
Comments on the optical lineshape function: Application to transient hole-burned spectra of bacterial reaction centers142(2015); http://dx.doi.org/10.1063/1.4913685View Description Hide Description
The vibrational spectral density is an important physical parameter needed to describe both linear and non-linear spectra of multi-chromophore systems such as photosynthetic complexes. Low-temperature techniques such as hole burning (HB) and fluorescence line narrowing are commonly used to extract the spectral density for a given electronic transition from experimental data. We report here that the lineshape function formula reported by Hayes et al. [J. Phys. Chem. 98, 7337 (1994)] in the mean-phonon approximation and frequently applied to analyzing HB data contains inconsistencies in notation, leading to essentially incorrect expressions in cases of moderate and strong electron-phonon (el-ph) coupling strengths. A corrected lineshape function L(ω) is given that retains the computational and intuitive advantages of the expression of Hayes et al. [J. Phys. Chem. 98, 7337 (1994)]. Although the corrected lineshape function could be used in modeling studies of various optical spectra, we suggest that it is better to calculate the lineshape function numerically, without introducing the mean-phonon approximation. New theoretical fits of the P870 and P960 absorption bands and frequency-dependent resonant HB spectra of Rb. sphaeroides and Rps. viridis reaction centers are provided as examples to demonstrate the importance of correct lineshape expressions. Comparison with the previously determined el-ph coupling parameters [Johnson et al., J. Phys. Chem. 94, 5849 (1990); Lyle et al., ibid. 97, 6924 (1993); Reddy et al., ibid. 97, 6934 (1993)] is also provided. The new fits lead to modified el-ph coupling strengths and different frequencies of the special pair marker mode, ω sp, for Rb. sphaeroides that could be used in the future for more advanced calculations of absorption and HB spectra obtained for various bacterial reaction centers.