Free-Energy and Molecular Dynamics Analysis of Ser195 and Thr195 Rotameric States in Thrombin
- 1University of Arkansas
Abstract
Ser195 is highly conserved as the catalytic nucleophile in trypsin like serine proteases, yet the molecular basis for the strong preference for serine over the chemically similar threonine residue remains unclear. Here, we used well tempered metadynamics and microsecond long molecular dynamics (MD) simulations to characterize the conformational energetics of S195 and T195 in thrombin. Well tempered metadynamics revealed distinct rotameric free energy landscapes in both the apo and the ligand bound systems. In the apo state, S195 preferentially occupied a lower angle rotameric state, while T195 exhibited two low energy states with similar free energies. Unbiased simulations further showed greater conformational heterogeneity and increased higher angle sampling for T195. Evaluation of the His57 to residue 195 hydrogen bond further revealed that S195 maintained close hydrogen bond distances across both rotameric regions, whereas T195 exhibited longer and more variable His57 to hydroxyl distances with reduced distance based occupancy. Analysis of the local environment showed that the additional T195 methyl group was associated with redistributed protein contacts. Additional rotamer conditioned analyses demonstrated differences in local packing within corresponding conformational states. Furthermore, solvent accessibility and water contact analyses revealed that the effects of the substitution on the catalytic side chain environment depend on rotameric conformation. Together, these findings demonstrate that two complementary simulation approaches converge on similar conformational preferences and suggest that the S195T substitution alters the rotameric free energy landscape in association with changes in catalytic geometry, local packing, and solvation. These results provide insight into the evolutionary preference for serine over threonine at the catalytic position of thrombin.
Related
- Structure-Guided Optimization of a LipL32-Targeting DNA Aptamer through Molecular Dynamics and Free-Energy Analysis
- Gaussian accelerated Molecular Dynamics - Thermodynamic Integration (GaMD-TI): Improved alchemical free energy calculations with enhanced sampling
- Structural insights into an aptainhibitor of the SARS-CoV-2 E protein using molecular docking, molecular dynamics, and biophysical analyses
- Allosteric Communication in the Tubulin Dimer Revealed through Molecular Dynamics and Network Analysis
- Edge-Aware Graph Attention Networks for Interpreting Biophysical Mechanisms from Molecular Dynamics Simulations
BibTeX
@misc{sauve2026freeenergy,
title = {{Free-Energy and Molecular Dynamics Analysis of Ser195 and Thr195 Rotameric States in Thrombin}},
author = {S. Sauve and P. Chaturvedi and M. Moradi},
howpublished = {bioRxiv},
year = {2026},
month = oct,
doi = {10.64898/2026.10.13.758922},
url = {https://www.biorxiv.org/content/10.64898/2026.10.13.758922v1},
note = {CC BY},
}
RIS
TY - UNPB
TI - Free-Energy and Molecular Dynamics Analysis of Ser195 and Thr195 Rotameric States in Thrombin
AU - S. Sauve
AU - P. Chaturvedi
AU - M. Moradi
PY - 2026
DA - 2026/10/13
PB - bioRxiv
DO - 10.64898/2026.10.13.758922
UR - https://www.biorxiv.org/content/10.64898/2026.10.13.758922v1
AB - Ser195 is highly conserved as the catalytic nucleophile in trypsin like serine proteases, yet the molecular basis for the strong preference for serine over the chemically similar threonine residue remains unclear. Here, we used well tempered metadynamics and microsecond long molecular dynamics (MD) simulations to characterize the conformational energetics of S195 and T195 in thrombin. Well tempered metadynamics revealed distinct rotameric free energy landscapes in both the apo and the ligand bound systems. In the apo state, S195 preferentially occupied a lower angle rotameric state, while T195 exhibited two low energy states with similar free energies. Unbiased simulations further showed greater conformational heterogeneity and increased higher angle sampling for T195. Evaluation of the His57 to residue 195 hydrogen bond further revealed that S195 maintained close hydrogen bond distances across both rotameric regions, whereas T195 exhibited longer and more variable His57 to hydroxyl distances with reduced distance based occupancy. Analysis of the local environment showed that the additional T195 methyl group was associated with redistributed protein contacts. Additional rotamer conditioned analyses demonstrated differences in local packing within corresponding conformational states. Furthermore, solvent accessibility and water contact analyses revealed that the effects of the substitution on the catalytic side chain environment depend on rotameric conformation. Together, these findings demonstrate that two complementary simulation approaches converge on similar conformational preferences and suggest that the S195T substitution alters the rotameric free energy landscape in association with changes in catalytic geometry, local packing, and solvation. These results provide insight into the evolutionary preference for serine over threonine at the catalytic position of thrombin.
ER -