Search papers, labs, and topics across Lattice.
The paper introduces a computational method that uses short, unbiased molecular dynamics trajectories and Harmonic Linear Discriminant Analysis (HLDA) to predict the impact of point mutations on the free-energy surface (FES) of peptides. They demonstrate that the HLDA eigenvector derived from wild-type trajectories can predict residue-level propensities for mutations to stabilize or destabilize the folded state of the CLN025 beta-hairpin. Furthermore, they show that shifts in the leading HLDA eigenvalue correlate with changes in the free-energy difference between states, validated against Replica Exchange Molecular Dynamics simulations.
Predict which mutations will stabilize a protein's folded state using only short simulations of the wild-type sequence.
Engineering the free-energy surfaces (FES) of proteins and peptides is central to controlling conformational ensembles and their responses to perturbations. However, predicting how chemical modifications such as point mutations reshape the FES and shift conformational equilibria remains challenging, particularly in data-scarce settings. Building on the Collective Variables for Free Energy Surface Tailoring (CV-FEST) framework, we develop a computational approach that leverages short, unbiased molecular dynamics trajectories to guide mutation analysis. Using the ten-residue beta-hairpin CLN025 and a systematic library of its single-point mutants, we apply Harmonic Linear Discriminant Analysis (HLDA) to extract collective variables from the conformational data. We find that the HLDA eigenvector learned solely from short wild-type trajectories provides residue-level insight into the propensity of mutations at specific positions to thermodynamically stabilize or destabilize the folded state. Extending this analysis, we show that shifts in the leading HLDA eigenvalue across mutants, a measure of changes in separability between the conformational ensembles along the HLDA coordinate, correlate strongly with mutation-induced changes in the free-energy difference between states, as reflected in melting temperatures. Benchmarked against Replica Exchange Molecular Dynamics simulations, these findings suggest a promising and computationally affordable route toward guiding the engineering of biomolecular free-energy landscapes.
