Surprisal Metrics for Quantifying Perturbed Conformational Dynamics in Markov State Models
Vincent A Voelz1, Brandon Elman1, Asghar M Razavi1
1Department of Chemistry, Temple University , Philadelphia, Pennsylvania 19122, United States.
Journal of Chemical Theory and Computation
|November 20, 2015
Summary
We developed a new surprisal metric to quantify changes in biomolecular dynamics transitions between states. This method helps identify key conformational changes due to mutations or environmental shifts.
Area of Science:
- Computational Biology
- Biophysics
- Molecular Dynamics
Background:
- Markov state models (MSMs) are widely used to analyze biomolecular conformational dynamics.
- Assessing changes in transition rates due to perturbations (e.g., mutations) is crucial but challenging.
- Existing methods may not adequately distinguish changes in transition rates from shifts in state definitions.
Purpose of the Study:
- To introduce a novel surprisal metric for quantifying differences in transition rates between two related Markov state models.
- To provide a method for assessing the impact of perturbations on biomolecular dynamics independent of state definition changes.
- To demonstrate the utility of this metric in identifying affected conformational states.
Main Methods:
- Development of a surprisal metric based on relative entropy and Jensen-Shannon divergence.
- Calculation of transition counts and statistical uncertainties for MSMs.
- Application of the metric to diverse biomolecular systems, including lattice models and all-atom simulations.
- Exploration of surprisal-based adaptive sampling for uncertainty reduction.
Main Results:
- The surprisal metric effectively quantifies differences in transition rates between closely related MSMs.
- The metric successfully identifies conformational states most sensitive to perturbations in protein hairpin and peptide models.
- Surprisal-based adaptive sampling significantly reduces statistical uncertainty in divergence calculations.
Conclusions:
- The surprisal metric offers a robust way to analyze changes in biomolecular dynamics under perturbation.
- This approach provides insights into how mutations and environmental factors alter molecular kinetics.
- Surprisal-based adaptive sampling can enhance the efficiency of molecular simulation and aid in ab initio design.
Related Concept Videos
¹H NMR of Conformationally Flexible Molecules: Temporal Resolution
1.4K
At room temperature, the chair conformer of cyclohexane undergoes rapid ring flipping between two equivalent chair conformers at a rate of approximately 105 times per second. These two chair conformers are in equilibrium. The rapid ring flipping results in the interconversion of the axial proton to an equatorial proton and an equatorial to the axial proton. Such interconversions are too rapid and cannot be detected on the NMR timescale. Hence, the NMR spectrometer cannot distinguish between the...
1.4K
¹H NMR of Conformationally Flexible Molecules: Variable-Temperature NMR
1.8K
The axial and equatorial protons in cyclohexane can be distinguished by performing a variable-temperature NMR experiment. In this process, except for one proton, the remaining eleven protons are replaced by deuterium. The deuterium substitution avoids the possible peak splitting caused by the spin-spin coupling between the adjacent protons. The remaining proton flips between the axial and equatorial positions.
1.8K
Entropy Change in Reversible Processes
3.4K
In the Carnot engine, which achieves the maximum efficiency between two reservoirs of fixed temperatures, the total change in entropy is zero. The observation can be generalized by considering any reversible cyclic process consisting of many Carnot cycles. Thus, it can be stated that the total entropy change of any ideal reversible cycle is zero.
The statement can be further generalized to prove that entropy is a state function. Take a cyclic process between any two points on a p-V diagram.
The statement can be further generalized to prove that entropy is a state function. Take a cyclic process between any two points on a p-V diagram.
3.4K
Standard Entropy Change for a Reaction
25.8K
Entropy is a state function, so the standard entropy change for a chemical reaction (ΔS°rxn) can be calculated from the difference in standard entropy between the products and the reactants.
25.8K
Cooperative Allosteric Transitions
2.8K
2.8K
Cooperative Allosteric Transitions
9.4K
Cooperative allosteric transitions can occur in multimeric proteins, where each subunit of the protein has its own ligand-binding site. When a ligand binds to any of these subunits, it triggers a conformational change that affects the binding sites in the other subunits; this can change the affinity of the other sites for their respective ligands. The ability of the protein to change the shape of its binding site is attributed to the presence of a mix of flexible and stable segments in the...
9.4K


