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.

Distinguishing between smooth and rough free energy barriers in protein folding.

Chevron Behavior and Isostable Enthalpic Barriers in Protein Folding: Successes and Limitations of Simple Gō-like Modeling.

Desolvation Barrier Effects Are a Likely Contributor to the Remarkable Diversity in the Folding Rates of Small Proteins.

Solvation and Desolvation Effects in Protein Folding: Native Flexibility, Kinetic Cooperativity and Enthalpic Barriers under Isostability Conditions

The Folding Transition State of Protein L is Extensive with Nonnative Interactions (and Not Small and Polarized)

Excluded volume, local structural cooperativity,and the polymer physics of protein folding rates

Reliable protein folding on non-funneled energy landscapes: the free energy reaction path

Probing Possible Downhill Folding: Native Contact Topology Likely Places a Significant Constraint on the Folding Cooperativity of Proteins with ∼40 Residues

Transition states for folding of circular-permuted proteins.

Measuring the energy landscape roughness and the transition state location of biomolecules using single molecule mechanical unfolding experiments

Topological Quantum Computation and Error Correction by Biological Cells

Kinetic and thermodynamic analysis of proteinlike heteropolymers: Monte Carlo histogram technique

Force Dependence of Proteins' Transition State Position and the Bell-Evans Model

Single-molecule chemo-mechanical unfolding reveals multiple transition state barriers in a small single-domain protein

All-atom Replica Exchange Molecular Simulation of Protein BBL

Statics, metastable states and barriers in protein folding: A replica variational approach

Watching proteins fold one molecule at a time

Disentangling folding from energetic traps in simulations of disordered proteins

Statistical mechanics of a correlated energy landscape model for protein folding funnels

Funnels, Pathways and the Energy Landscape of Protein Folding: A Synthesis

The Clustering of AGN and Galaxies at Intermediate Redshift