Spectral gap optimization of order parameters for sampling complex molecular systems
Open accessAbstract
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.
Integrity signals
Registry
Cites
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.
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”
From Unsupervised machine learning in atomistic simulations, between predictions and understanding (2019) · cites this paper as [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.”
From Unsupervised machine learning in atomistic simulations, between predictions and understanding (2019) · cites this paper as [82]
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
A climbing image nudged elastic band method for finding saddle points and minimum energy paths
Henkelman, G., Uberuaga, B. P., & Jónsson, H. (2000) · The Journal of Chemical Physics · doi:10.1063/1.1329672 · Cited by 21,576 · see citing works
Canonical sampling through velocity rescaling
Bussi, G., Donadio, D., & Parrinello, M. (2007) · The Journal of Chemical Physics · doi:10.1063/1.2408420 · Cited by 18,034 · see citing works
GROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation
Hess, B., Kutzner, C., van der Spoel, D., & Lindahl, E. R. (2008) · Journal of Chemical Theory and Computation · doi:10.1021/ct700301q · Cited by 15,830 · see citing works
Information theory and statistical mechanics
Jaynes, E. T. (1957) · Physical Review · doi:10.1103/physrev.106.620 · Cited by 12,931 · see citing works
Unified Approach for Molecular Dynamics and Density-Functional Theory
Car, R., & Parrinello, M. (1985) · Physical Review Letters · doi:10.1103/physrevlett.55.2471 · Cited by 10,737 · see citing works
Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling
Torrie, G. M., & Valleau, J. P. (1977) · Journal of Computational Physics · doi:10.1016/0021-9991(77)90121-8 · Cited by 6,091 · see citing works
Laio, A., & Parrinello, M. (2002) · Proceedings of the National Academy of Sciences · doi:10.1073/pnas.202427399 · Cited by 5,626 · see citing works
Pecora, R. (1985) · Springer eBooks · doi:10.1007/978-1-4613-2389-1 · Cited by 4,863 · see citing works
Bond-orientational order in liquids and glasses
Steinhardt, P. J., Nelson, D. R., & Ronchetti, M. (1983) · Physical review. B, Condensed matter · doi:10.1103/physrevb.28.784 · Cited by 3,511 · see citing works
PLUMED 2: New feathers for an old bird
Tribello, G. A., Bonomi, M., Branduardi, D., Camilloni, C., & Bussi, G. (2013) · Computer Physics Communications · doi:10.1016/j.cpc.2013.09.018 · Cited by 3,427 · see citing works
Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method
Barducci, A., Bussi, G., & Parrinello, M. (2008) · Physical Review Letters · doi:10.1103/physrevlett.100.020603 · Cited by 3,184 · see citing works
Brownian dynamics with hydrodynamic interactions
Ermak, D., & McCammon, J. (1978) · The Journal of Chemical Physics · doi:10.1063/1.436761 · Cited by 2,238 · see citing works
Nonequilibrium Statistical Mechanics
Zwanzig, R. W. (2001) · doi:10.1093/oso/9780195140187.001.0001 · Cited by 2,075 · see citing works
Bolhuis, P. G., Chandler, D., Dellago, C., & Geissler, P. L. (2002) · Annual Review of Physical Chemistry · doi:10.1146/annurev.physchem.53.082301.113146 · Cited by 1,972 · see citing works
Coifman, R. R., Lafon, S., Lee, A. B., Maggioni, M., Nadler, B., Warner, F., & Zucker, S. W. (2005) · Proceedings of the National Academy of Sciences · doi:10.1073/pnas.0500334102 · Cited by 1,675 · see citing works
Drug-target residence time and its implications for lead optimization
Copeland, R. A., Pompliano, D. L., & Meek, T. D. (2006) · Nature Reviews Drug Discovery · doi:10.1038/nrd2082 · Cited by 1,464 · see citing works
Nonequilibrium Statistical Mechanics
Zwanzig, R. W. (1990) · doi:10.1142/1009 · Cited by 1,339 · see citing works
Constrained reaction coordinate dynamics for the simulation of rare events
Carter, E. A., Ciccotti, G., Hynes, J. T., & Kapral, R. (1989) · Chemical Physics Letters · doi:10.1016/s0009-2614(89)87314-2 · Cited by 1,023 · see citing works
Hyperdynamics: Accelerated Molecular Dynamics of Infrequent Events
Voter, A. F. (1997) · Physical Review Letters · doi:10.1103/physrevlett.78.3908 · Cited by 996 · see citing works
Identification of slow molecular order parameters for Markov model construction
Pérez-Hernández, G., Paul, F., Giorgino, T., De Fabritiis, G., & Noé, F. (2013) · The Journal of Chemical Physics · doi:10.1063/1.4811489 · Cited by 990 · cited 2× · see citing works
“The identification of such order parameters poses another complicated problem, but as routinely done in other methods aimed at optimizing CVs (15, 18, 24), we assume such order parameters are a priori known.”
BackgroundAdaptive biasing force method for scalar and vector free energy calculations.
Darve, E. F., Rodríguez-Gómez, D., & Pohorille, A. (2008) · Journal of Chemical Physics · doi:10.1063/1.2829861 · Cited by 793 · see citing works
Gear, C. W., Hyman, J. M., Kevrekidid, P. G., Kevrekidis, I. G., Runborg, O., & Theodoropoulos, C. (2003) · Communications in Mathematical Sciences · doi:10.4310/cms.2003.v1.n4.a5 · Cited by 778 · see citing works
Enhancing Important Fluctuations: Rare Events and Metadynamics from a Conceptual Viewpoint
Valsson, Ó., Tiwary, P., & Parrinello, M. (2016) · Annual Review of Physical Chemistry · doi:10.1146/annurev-physchem-040215-112229 · Cited by 619 · see citing works
Self-citationIdentification of slow molecular order parameters for Markov model constructionPreprint
Perez-hern, G., Paul, F., Toni Giorgino, A., Fabritiis, G. D., Noé, F., & Contribution, E. (2013) · Cited by 613 · see citing works
From A to B in free energy space.
Branduardi, D., Gervasio, F., & Parrinello, M. (2007) · Journal of Chemical Physics · doi:10.1063/1.2432340 · Cited by 596 · see citing works
Computing time scales from reaction coordinates by milestoning.
Faradjian, A., & Elber, R. (2004) · Journal of Chemical Physics · doi:10.1063/1.1738640 · Cited by 596 · see citing works
Hummer, G. (2005) · New Journal of Physics · doi:10.1088/1367-2630/7/1/034 · Cited by 530 · see citing works
Goldburg, W. I. (1999) · American Journal of Physics · doi:10.1119/1.19101 · Cited by 526 · see citing works
Ladbury, J. E. (1996) · Chemistry & Biology · doi:10.1016/s1074-5521(96)90164-7 · Cited by 521 · see citing works
A Time-Independent Free Energy Estimator for Metadynamics
Tiwary, P., & Parrinello, M. (2014) · The Journal of Physical Chemistry B · doi:10.1021/jp504920s · Cited by 516 · see citing works
Self-citationReaction coordinates and rates from transition paths
Best, R. B., & Hummer, G. (2005) · Proceedings of the National Academy of Sciences · doi:10.1073/pnas.0408098102 · Cited by 512 · see citing works
Tiwary, P., & Parrinello, M. (2013) · Physical Review Letters · doi:10.1103/physrevlett.111.230602 · Cited by 494 · see citing works
Self-citationDewetting and Hydrophobic Interaction in Physical and Biological Systems
Berne, B. J., Weeks, J. D., & Zhou, R. (2008) · Annual Review of Physical Chemistry · doi:10.1146/annurev.physchem.58.032806.104445 · Cited by 464 · see citing works
Self-citationThe Minimum Entropy Production Principle
Jaynes, E. (1980) · Annual Review of Physical Chemistry · doi:10.1146/annurev.pc.31.100180.003051 · Cited by 456 · see citing works
Obtaining reaction coordinates by likelihood maximization
Peters, B., & Trout, B. L. (2006) · The Journal of Chemical Physics · doi:10.1063/1.2234477 · Cited by 417 · see citing works
Rate processes with dynamical disorder
Zwanzig, R. W. (1990) · Accounts of Chemical Research · doi:10.1021/ar00173a005 · Cited by 409 · see citing works
Automatic Method for Identifying Reaction Coordinates in Complex Systems
Ma, A., & Dinner, A. R. (2005) · The Journal of Physical Chemistry B · doi:10.1021/jp045546c · Cited by 400 · see citing works
Principles of maximum entropy and maximum caliber in statistical physics
Pressé, S., Ghosh, K., Lee, J., & Dill, K. A. (2013) · Reviews of Modern Physics · doi:10.1103/revmodphys.85.1115 · Cited by 399 · see citing works
Hansmann, U. H. E., & Okamoto, Y. (1993) · Journal of Computational Chemistry · doi:10.1002/jcc.540141110 · Cited by 395 · see citing works
Kinetics of protein–ligand unbinding: Predicting pathways, rates, and rate-limiting steps
Tiwary, P., Limongelli, V., Salvalaglio, M., & Parrinello, M. (2015) · Proceedings of the National Academy of Sciences · doi:10.1073/pnas.1424461112 · Cited by 348 · see citing works
Self-citationSimplifying the representation of complex free-energy landscapes using sketch-map
Ceriotti, M., Tribello, G. A., & Parrinello, M. (2011) · Proceedings of the National Academy of Sciences · doi:10.1073/pnas.1108486108 · Cited by 340 · see citing works
Variational Approach to Enhanced Sampling and Free Energy Calculations
Valsson, O., & Parrinello, M. (2014) · Physical Review Letters · doi:10.1103/physrevlett.113.090601 · Cited by 336 · see citing works
Diffusion Maps, Reduction Coordinates, and Low Dimensional Representation of Stochastic Systems
Coifman, R., Kevrekidis, I., Lafon, S., Maggioni, M., & Nadler, B. (2008) · Multiscale Modeling & simulation · doi:10.1137/070696325 · Cited by 323 · see citing works
Metadynamics with Adaptive Gaussians
Branduardi, D., Bussi, G., & Parrinello, M. (2012) · Journal of Chemical Theory and Computation · doi:10.1021/ct3002464 · Cited by 306 · see citing works
Random walk in orthogonal space to achieve efficient free-energy simulation of complex systems
Zheng, L.-Q., Chen, M., & Yang, W. (2008) · Proceedings of the National Academy of Sciences of the United States of America · doi:10.1073/pnas.0810631106 · Cited by 301 · see citing works
Determination of reaction coordinates via locally scaled diffusion map
Rohrdanz, M. A., Zheng, W., Maggioni, M., & Clementi, C. (2011) · The Journal of Chemical Physics · doi:10.1063/1.3569857 · Cited by 292 · see citing works
A Kinetic Model of Trp-Cage Folding from Multiple Biased Molecular Dynamics Simulations
Marinelli, F., Pietrucci, F., Laio, A., & Piana, S. (2009) · PLoS Computational Biology · doi:10.1371/journal.pcbi.1000452 · Cited by 274 · see citing works
Parameter Space Compression Underlies Emergent Theories and Predictive Models
Machta, B. B., Chachra, R., Transtrum, M. K., & Sethna, J. P. (2013) · Science · doi:10.1126/science.1238723 · Cited by 273 · see citing works
Metainference: A Bayesian inference method for heterogeneous systems
Bonomi, M., Camilloni, C., Cavalli, A., & Vendruscolo, M. (2016) · Science Advances · doi:10.1126/sciadv.1501177 · Cited by 236 · see citing works
Large-scale conformational sampling of proteins using temperature-accelerated molecular dynamics
Abrams, C. F., & Vanden‐Eijnden, E. (2010) · Proceedings of the National Academy of Sciences · doi:10.1073/pnas.0914540107 · Cited by 209 · see citing works
Assessing the Reliability of the Dynamics Reconstructed from Metadynamics
Salvalaglio, M., Tiwary, P., & Parrinello, M. (2014) · Journal of Chemical Theory and Computation · doi:10.1021/ct500040r · Cited by 207 · see citing works
Self-citationFree Energies from Dynamic Weighted Histogram Analysis Using Unbiased Markov State Model
Rosta, E., & Hummer, G. (2015) · Journal of Chemical Theory and Computation · doi:10.1021/ct500719p · Cited by 144 · see citing works
Characterization of the free-energy landscapes of proteins by NMR-guided metadynamics
Granata, D., Camilloni, C., Vendruscolo, M., & Laio, A. (2013) · Proceedings of the National Academy of Sciences · doi:10.1073/pnas.1218350110 · Cited by 143 · see citing works
Electron transfer reaction dynamics in non-Debye solvents
Bicout, D., & Szabó, A. (1998) · The Journal of Chemical Physics · doi:10.1063/1.476800 · Cited by 143 · see citing works
Intramolecular rate process: Isomerization dynamics and the transition to chaos
Leon, N. D., & Berne, B. J. (1981) · The Journal of Chemical Physics · doi:10.1063/1.442459 · Cited by 129 · see citing works
Self-citationTime scale separation leads to position-dependent diffusion along a slow coordinate
Berezhkovskii, A. M., & Szabó, A. (2011) · The Journal of Chemical Physics · doi:10.1063/1.3626215 · Cited by 114 · see citing works
Numerical calculation of the rate of homogeneous gas–liquid nucleation in a Lennard-Jones system
ten Wolde, P. R., Ruiz-Montero, M. J., & Frenkel, D. (1999) · The Journal of Chemical Physics · doi:10.1063/1.477799 · Cited by 101 · see citing works
Molecular Reorientation in Liquids and Gases
Berne, B. J., Pechukas, P., & Harp, G. D. (1968) · The Journal of Chemical Physics · doi:10.1063/1.1670559 · Cited by 100 · see citing works
Self-citationLocating landmarks on high-dimensional free energy surfaces
Chen, M., Yu, T.-Q., & Tuckerman, M. E. (2015) · Proceedings of the National Academy of Sciences · doi:10.1073/pnas.1418241112 · Cited by 84 · see citing works
Role of water and steric constraints in the kinetics of cavity–ligand unbinding
Tiwary, P., Mondal, J., Morrone, J. A., et al. (2015) · Proceedings of the National Academy of Sciences · doi:10.1073/pnas.1516652112 · Cited by 82 · see citing works
Self-citationInferring Transition Rates of Networks from Populations in Continuous-Time Markov Processes
Dixit, P. D., Jain, A., Stock, G., & Dill, K. A. (2015) · Journal of Chemical Theory and Computation · doi:10.1021/acs.jctc.5b00537 · Cited by 67 · see citing works
Interplay between Hydrodynamics and the Free Energy Surface in the Assembly of Nanoscale Hydrophobes
Morrone, J. A., Li, J., & Berne, B. J. (2011) · The Journal of Physical Chemistry B · doi:10.1021/jp209568n · Cited by 54 · see citing works
Self-citationWell-Tempered Variational Approach to Enhanced Sampling
Valsson, Ó., & Parrinello, M. (2015) · Journal of Chemical Theory and Computation · doi:10.1021/acs.jctc.5b00076 · Cited by 49 · see citing works
Tiwary, P., & van de Walle, A. (2013) · Physical Review B · doi:10.1103/physrevb.87.094304 · Cited by 34 · see citing works
Self-citationMethod of the most probable path of evolution in the theory of stationary irreversible processes
Filyukov, A. A., & Karpov, V. Y. (1967) · Journal of Engineering Physics and Thermophysics · doi:10.1007/bf00828961 · Cited by 28 · see citing works
Random walk in orthogonal space to achieve efficient free-energy simulation of complex systems
Zhenga, L. (2016) · Cited by 9 · see citing works
Inferring transition rates on networks with incomplete knowledgePreprint
Dixit, P. D., Jain, A., Stock, G., & Dill, K. A. (2015) · arXiv (Cornell University) · doi:10.48550/arxiv.1504.01277 · Cited by 2 · see citing works
The Independence of the Junior Scientist's Mind: At What Price?
Scoles, G. (2016) · Annual Review of Physical Chemistry · doi:10.1146/annurev-physchem-040215-112236 · Cited by 2 · see citing works
Flagged first, then by citations
Research works printed but not linked: 11
- 9 C. F. Abrams and E. Vanden-Eijnden, Proc. Natl. Acad. Sci. 107 , 4961 (2010).
- 10 L. Zheng, M. Chen, and W. Yang, Proc. Natl. Acad. Sci. 105 , 20227 (2008).
- 16 R. B. Best and G. Hummer, Proc. Natl. Acad. Sci. 102 , 6732 (2005).
- 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).
- 19 A. Ma and A. R. Dinner, The Journal of Physical Chemistry B 109 , 6769 (2005).
- 21 M. Ceriotti, G. A. Tribello, and M. Parrinello, Proc. Natl. Acad. Sci. 108 , 13023 (2011).
- 22 M. Chen, T.-Q. Yu, and M. E. Tuckerman, Proc. Natl. Acad. Sci. 112 , 3235 (2015).
- 28 S. Press´ e, K. Ghosh, J. Lee, and K. A. Dill, Rev. Mod. Phys. 85 , 1115 (2013).
- 46 D. Granata, C. Camilloni, M. Vendruscolo, and A. Laio, Proc. Natl. Acad. Sci. 110 , 6817 (2013).
- 47 M. Bonomi, C. Camilloni, A. Cavalli, and M. Vendruscolo, arXiv preprint arXiv:1509.05684 (2015). · arXiv
- 49 P. Tiwary, V. Limongelli, M. Salvalaglio, and M. Parrinello, Proc. Natl. Acad. Sci. 112 , E386 (2015).
Cited by
Sorted by citations
Machine learning approaches for analyzing and enhancing molecular dynamics simulations
Wang, Y., Lamim Ribeiro, J. M., & Tiwary, P. (2020) · Current Opinion in Structural Biology · doi:10.1016/j.sbi.2019.12.016 · Cited by 303 · see citing works
Self-citationUnsupervised machine learning in atomistic simulations, between predictions and understanding
Ceriotti, M. (2019) · The Journal of Chemical Physics · doi:10.1063/1.5091842 · Cited by 173 · see citing works
Thermal Fluctuations of Singular Bar-Joint Mechanisms
Mannattil, M., Schwarz, J. M., & Santangelo, C. D. (2022) · Physical Review Letters · doi:10.1103/physrevlett.128.208005 · Cited by 7 · see citing works
- Authors
- Tiwary, P., & Berne, B. J.
- Affiliations
- Columbia University (US)
- ORCID
- Tiwary, Pratyush
- Venue
- Proceedings of the National Academy of Sciences
- Publisher
- National Academy of Sciences
- Published
- 2016-02-29
- Volume
- 113
- Issue
- 11
- Pages
- 2839-2844
- ISSN
- 0027-8424 (ISSN-L)
- 1091-6490
- PMC
- PMC4801247
- 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
- green
- www.ncbi.nlm.nih.gov
- 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