Comparing Enhanced Sampling Methods in Exploring the Conformational Space of β-Catenin 17–48
Open Access
- 24 June 2026
- journal article
- Published by American Chemical Society (ACS) in The Journal of Physical Chemistry B
- Vol. 130 (29) , 7323-7333
Abstract
Intrinsically disordered proteins (IDPs) play critical roles in cellular signaling and regulation, yet their dynamic conformational landscapes make them difficult to characterize experimentally and computationally. Phosphorylation, one of the most common post-translational modifications, frequently occurs within intrinsically disordered regions and can modulate protein structure and function. Enhanced sampling molecular dynamics methods offer a potential route to more efficiently explore the diverse conformations accessible to IDPs, but systematic comparisons of their performance sampling the conformational landscape of such systems remain limited. Here, we evaluate the conformational sampling of the intrinsically disordered β-catenin17-48 peptide in both its nonphosphorylated and phosphorylated states using three enhanced sampling approaches: Gaussian-accelerated molecular dynamics (GaMD), metadynamics (METAD), and weighted ensemble simulations (WESTPA), in comparison with conventional molecular dynamics simulations. Two collective variables (CVs) were explored to guide sampling: the ϕ dihedral angles of the phosphorylation sites Ser33 and Ser37 and the end-to-end distance of the peptide. We found that different enhanced sampling methods explored distinct regions of conformational space rather than converging to a single ensemble, with GaMD largely overlapping with unbiased simulations. Notably, METAD and WESTPA more readily accessed conformational regions not observed in unbiased simulations. Analysis of the combined conformational ensembles identified intermediate conformations connecting the nonphosphorylated and phosphorylated states, which are preferentially sampled in simulations employing adaptive strategies. Additionally, the Ser33/Ser37 ϕ angle CV more effectively captures phosphorylation-dependent conformational shifts than the end-to-end distance metric. Together, these results highlight how both the choice of enhanced sampling strategy and the selection of collective variables influence the exploration of IDP conformational landscapes.Keywords
Funding Information
- National Institute of Food and Agriculture (VA-160211)
- National Institute of General Medical Sciences (R35GM133754)
This publication has 58 references indexed in Scilit:
- Weighted Ensemble Simulations With the Drude Polarizable ModelJournal of Computational Chemistry, 2025
- Base pair dynamics, electrostatics, and thermodynamics at the LTR-III quadruplex:duplex junctionBiophysical Journal, 2024
- MDTraj: A Modern Open Library for the Analysis of Molecular Dynamics TrajectoriesBiophysical Journal, 2015
- MDAnalysis: A toolkit for the analysis of molecular dynamics simulationsJournal of Computational Chemistry, 2011
- SPARTA+: a modest improvement in empirical NMR chemical shift prediction by means of an artificial neural networkJournal of Biomolecular NMR, 2010
- TauopathiesCellular and Molecular Life Sciences, 2007
- Assessing protein disorder and induced foldingProteins: Structure, Function, and Bioinformatics, 2005
- Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water modelsJournal of Computational Chemistry, 1992
- Dictionary of protein secondary structure: Pattern recognition of hydrogen‐bonded and geometrical featuresBiopolymers, 1983
- Numerical integration of the cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanesJournal of Computational Physics, 1977