Biswas Arumay, Eisert-Sasse Riley K, Okafor C Denise
Department of Chemistry, The Pennsylvania State University, University Park, PA, 16802, USA.
Department of Biochemistry and Molecular Biology, The Pennsylvania State University, University Park, PA, 16802, USA.
Chemphyschem. 2025 Feb 1;26(3):e202400783. doi: 10.1002/cphc.202400783. Epub 2024 Nov 20.
Molecular dynamics (MD) simulations are immensely valuable for studying protein structure, function and dynamics. Their ability to capture atomic-level behavior of molecules and describe their evolution over time makes it a powerful synergistic tool for biochemistry, structural biology and other life sciences. To advance research and knowledge on reasonable timescales, researchers must optimize the amount of useful information extracted from simulation data while often frugally managing computational resources. Often, this involves balancing the length of MD trajectories with the number of replicas of a given system, with the aim of maximizing sampling of the conformational landscape. However, identifying this balance is not always intuitive, and the lack of standards among researchers can produce large variability in results and predictions from MD measurements. Here, we investigate the variability in MD results when simulation length and replica numbers are varied. Using a 231-amino acid domain, we compare measurements from independent trajectories to a benchmark trajectory of 3, 1000-ns replicates. We perform these simulations on 27 protein-ligand complexes, allowing us to compare ligand-specific rankings of complexes across independent replicas. Our results reveal that some MD measurements are accurately ranked by single trajectories, while others are not. We uncover similar variability in the effects of trajectory lengths on measurements. Our findings suggest that a one-size-fits-all approach to MD simulations is not necessarily the best approach, and depending on the intended measurements and research question, it may be advantageous sometimes to prioritize longer trajectories over multiple replicas. This work provides important considerations for researchers while designing simulation studies.
分子动力学（MD）模拟对于研究蛋白质结构、功能和动力学极为重要。其捕捉分子原子水平行为并描述其随时间演变的能力，使其成为生物化学、结构生物学和其他生命科学中强大的协同工具。为了在合理的时间尺度上推进研究和知识，研究人员必须优化从模拟数据中提取的有用信息量，同时通常要节约地管理计算资源。通常，这涉及平衡MD轨迹的长度与给定系统的复制品数量，目的是最大化构象景观的采样。然而，确定这种平衡并不总是直观的，研究人员之间缺乏标准会导致MD测量结果和预测产生很大的变异性。在这里，我们研究了当模拟长度和复制品数量变化时MD结果的变异性。使用一个231个氨基酸的结构域，我们将独立轨迹的测量结果与3个1000纳秒复制品的基准轨迹进行比较。我们对27种蛋白质-配体复合物进行这些模拟，使我们能够比较不同独立复制品中复合物的配体特异性排名。我们的结果表明，一些MD测量可以通过单个轨迹准确排名，而其他测量则不能。我们还发现轨迹长度对测量的影响存在类似的变异性。我们的研究结果表明，一刀切的MD模拟方法不一定是最好的方法，根据预期的测量和研究问题，有时优先考虑更长的轨迹而不是多个复制品可能是有利的。这项工作为研究人员在设计模拟研究时提供了重要的考虑因素。
