Abstract:Protein conformational heterogeneity plays an essential role in a myriad of different biological processes. Extensive conformational heterogeneity is especially characteristic of intrinsically disordered proteins and protein regions (collectively IDRs), which lack a well-defined three-dimensional structure and instead rapidly exchange between a diverse ensemble of configurations. An emerging paradigm recognizes that the conformational biases encoded in IDR ensembles can play a central role in their biological function, necessitating understanding these sequence-ensemble relations. All-atom simulations have provided critical insight into our modern understanding of the solution behavior of IDRs. However, effectively exploring the accessible conformational space associated with large, heterogeneous ensembles is challenging. In particular, identifying poorly sampled or energetically trapped regions of disordered proteins in simulations often relies on qualitative assessment based on visual inspection of simulations and/or analysis data. These approaches, while convenient, run the risk of masking poorly-sampled simulations. In this work, we present an algorithm for quantifying per-residue local conformational heterogeneity in protein simulations. Our work builds on prior work and compares the similarity between backbone dihedral angle distributions generated from molecular simulations in a limiting polymer model and across independent all-atom simulations. In this regime, the polymer model serves as a statistical reference model for extensive conformational heterogeneity in a real chain. Quantitative comparisons of probability vectors generated from these simulations reveal the extent of conformational sampling in a simulation, enabling us to distinguish between situations in which protein regions are well-sampled, poorly-sampled, or folded. To demonstrate the effectiveness of this approach, we apply our algorithm to several toy, synthetic, and biological systems. Accurately assessing local conformational sampling in simulations of IDRs will help better quantify new enhanced sampling methods, ensure force field comparisons are equivalent, and provide confidence that conclusions drawn from simulations are robust.

A coarse‐grained model for disordered and multi‐domain proteins

A coarse-grained model for disordered and multi-domain proteins

Learning Effective Molecular Models from Experimental Observables.

Combining molecular dynamics simulations with small-angle X-ray and neutron scattering data to study multi-domain proteins in solution.

A Discontinuous Potential Model for Protein-Protein Interactions.

A coarse-grained model of disordered RNA for simulations of biomolecular condensates

COCOMO2: A coarse-grained model for interacting folded and disordered proteins

SOP-MULTI: A Self-Organized Polymer-Based Coarse-Grained Model for Multidomain and Intrinsically Disordered Proteins with Conformation Ensemble Consistent with Experimental Scattering Data

Evaluation of Essential Dynamics and Fixed-Length Coarse Graining for Multidomain Proteins

SOP-MULTI: A self-organized polymer based coarse-grained model for multi-domain and intrinsically disordered proteins with conformation ensemble consistent with experimental scattering data

Cooperativity, Local-Nonlocal Coupling, and Nonnative Interactions: Principles of Protein Folding from Coarse-Grained Models

Monte Carlo Sampling of Protein Folding by Combining an All-Atom Physics-Based Model with a Native State Bias

Cooperativity and Folding Kinetics in a Multi-Domain Protein with Interwoven Chain Topology

Consistent View of Protein Fluctuations from All-Atom Molecular Dynamics and Coarse-Grained Dynamics with Knowledge-Based Force-Field.

Unbiased simulation of structural transitions in calmodulin

Disentangling folding from energetic traps in simulations of disordered proteins

Contact-based molecular dynamics of structured and disordered proteins in a coarse-grained model: fixed contacts, switchable contacts and those described by pseudo-improper-dihedral angles

Minimal models for proteins and RNA: From folding to function

Transferable coarse-grained potential for $\textit{de novo}$ protein folding and design

Double-Well Ultra-Coarse-Grained Model to Describe Protein Conformational Transitions

Coarse-graining protein structures with local multivariate features from molecular dynamics.