Cited by 259Advanced Physical and Chemical Molecular InteractionsProtein Structure and DynamicsMarkov Chains and Monte Carlo Methods

Spectral gap optimization of order parameters for sampling complex molecular systems

Open access
Tiwary, Pratyush, Berne, B. J.·2016·Proceedings of the National Academy of Sciences

Abstract

In modern-day simulations of many-body systems, much of the computational complexity is shifted to the identification of slowly changing molecular order parameters called collective variables (CVs) or reaction coordinates. A vast array of enhanced-sampling methods are based on the identification and biasing of these low-dimensional order parameters, whose fluctuations are important in driving rare events of interest. Here, we describe a new algorithm for finding optimal low-dimensional CVs for use in enhanced-sampling biasing methods like umbrella sampling, metadynamics, and related methods, when limited prior static and dynamic information is known about the system, and a much larger set of candidate CVs is specified. The algorithm involves estimating the best combination of these candidate CVs, as quantified by a maximum path entropy estimate of the spectral gap for dynamics viewed as a function of that CV. The algorithm is called spectral gap optimization of order parameters (SGOOP). Through multiple practical examples, we show how this postprocessing procedure can lead to optimization of CV and several orders of magnitude improvement in the convergence of the free energy calculated through metadynamics, essentially giving the ability to extract useful information even from unsuccessful metadynamics runs.

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1073/pnas.1600917113

Cites

68 of 68 references verified0 citation errors · 0 not checkable

Statistics

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

  • SGOOP singles out the linear combination of input features that correspond to the direction that is most relevant for the long-time dynamics of the system being studied.

    “the spectral gap optimization of order parameters (SGOOP) method to single out the linear combination of input features that correspond to the direction that is most relevant for the long-time dynamics of the system being studied 82”
  • Dimensionality reduction algorithms such as TICA or SGOOP, that incorporate information on the system's dynamics, implicitly ensure that the order parameter(s) perform well as descriptors of the slowest molecular motion.

    “Dimensionality reduction algorithms as TICA 80,81 or SGOOP 82 , that incorporate information on the system's dynamics, implicitly ensure that the order parameter(s) perform well as descriptors of the slowest molecular motion.”

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Cites

Flagged first, then by citations

Research works printed but not linked: 11

  1. 9 C. F. Abrams and E. Vanden-Eijnden, Proc. Natl. Acad. Sci. 107 , 4961 (2010).
  2. 10 L. Zheng, M. Chen, and W. Yang, Proc. Natl. Acad. Sci. 105 , 20227 (2008).
  3. 16 R. B. Best and G. Hummer, Proc. Natl. Acad. Sci. 102 , 6732 (2005).
  4. 17 R. R. Coifman, S. Lafon, A. B. Lee, M. Maggioni, B. Nadler, F. Warner, and S. W. Zucker, Proc. Natl. Acad. Sci. 102 , 7426 (2005).
  5. 19 A. Ma and A. R. Dinner, The Journal of Physical Chemistry B 109 , 6769 (2005).
  6. 21 M. Ceriotti, G. A. Tribello, and M. Parrinello, Proc. Natl. Acad. Sci. 108 , 13023 (2011).
  7. 22 M. Chen, T.-Q. Yu, and M. E. Tuckerman, Proc. Natl. Acad. Sci. 112 , 3235 (2015).
  8. 28 S. Press´ e, K. Ghosh, J. Lee, and K. A. Dill, Rev. Mod. Phys. 85 , 1115 (2013).
  9. 46 D. Granata, C. Camilloni, M. Vendruscolo, and A. Laio, Proc. Natl. Acad. Sci. 110 , 6817 (2013).
  10. 47 M. Bonomi, C. Camilloni, A. Cavalli, and M. Vendruscolo, arXiv preprint arXiv:1509.05684 (2015). · arXiv
  11. 49 P. Tiwary, V. Limongelli, M. Salvalaglio, and M. Parrinello, Proc. Natl. Acad. Sci. 112 , E386 (2015).

Cited by

Sorted by citations

Authors
Tiwary, P., & Berne, B. J.
Affiliations
Columbia University (US)
Venue
Proceedings of the National Academy of Sciences
Published
2016-02-29
Volume
113
Issue
11
Pages
2839-2844
ISSN
  • 0027-8424 (ISSN-L)
  • 1091-6490
OpenAlex
W2246932694
Cited by
259
References
68
Citations by year
2026: 14 · 2025: 20 · 2024: 27 · 2023: 29 · 2022: 32 · 2021: 27 · 2020: 24 · 2019: 28 · 2018: 29 · 2017: 19 · 2016: 8 · 2015: 1 · 2012: 1
FWCI
19.11
Type
journal article
Language
en
Open access
Locations
  • Proceedings of the National Academy of Sciences (publishedVersion)
  • Europe PMC (PubMed Central) (publishedVersion)
  • PubMed Central (publishedVersion)
  • PubMed (publishedVersion)
  • arXiv (Cornell University) (submittedVersion)
  • arXiv (Cornell University)
Funders
  • HHS | NIH | National Institute of General Medical Sciences (NIH-GM4330)
  • NSF | CISE | Advanced Cyberinfrastructure (TG-MCA08X002)
Topics
  • Physical Sciences › Chemistry › Physical and Theoretical Chemistry › Advanced Physical and Chemical Molecular Interactions
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Protein Structure and Dynamics
  • Physical Sciences › Mathematics › Statistics and Probability › Markov Chains and Monte Carlo Methods
Keywords
collective variables, reaction coordinate, enhanced sampling, metadynamics, umbrella sampling, path entropy, free energy landscape, rare events, molecular dynamics
MeSH terms
Algorithms, Models, Molecular, Models, Theoretical, Peptides, Entropy, Molecular Dynamics Simulation

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review