Telescoping-Solvation-Box Protocol-Based Umbrella Sampling Method for Potential Mean Force Estimation

Abstract

Previously, it was shown that the telescoping box scheme, in combination with adaptive steered molecule dynamics (ASMD), can be used to estimate the potential of mean force (PMF) with a decrease in computational cost associated with large solvation boxes. Since ASMD reduces to umbrella sampling (US) in the limit of infinitely slow pulling velocity, a hypothesis was made that the telescoping box scheme can be extended to include the US method. The hypothesis was tested using the unfolding pathway of a polyalanine peptide in a water box and translocation of α-tocopherol through the human membrane. Two different approaches were tried: telescoping US (TELUS), in which the number of solvent molecules was linearly coupled to the reaction coordinate, and Block-TELUS, which was a compromise between the fixed solvation box of the US and the window-wise variable solvation box of TELUS. In the case of polyalanine peptide in a water box, both approaches gave overlapping potential of mean force (PMF) concerning the benchmark US-PMF. Window-wise comparison of properties like root-mean-square inner product, Ramachandran plot, α-helix content, and hydrogen bond formation was used to verify that both approaches captured the same dynamics as the US method. Applying the Block-TELUS protocol in the system with diffusing α-tocopherol through the bilayer resulted in overlapping PMF to its US benchmark. A comparison between the window-wise orientation of the chromanol headgroup also produced almost identical results. This study concluded that with the careful application of telescoping solvation boxes, a less computationally expensive US could be performed for systems where the effect of distant solvent molecules on the configurational space sampled in the window depends on the value of the reaction coordinate.

Uncovering the Structural and Binding Insights of Dual Inhibitors Simultaneously Targeting Two Distinct Sites on EGFR Kinase

“…A total of 6 ns simulation was performed in each window followed by short energy minimization and equilibration to avoid the unnecessary overlapping of the atom coordinates. The WHAM − calculation was performed for each umbrella window along the RC, − and the Monte Carlo bootstrap , was performed to estimate the statistical uncertainty of the PMF profile, and the tolerance for iteration was set to be 0.000001 to obtain a convergent free energy plot. To keep the simulation consistent with the prior one, the temperature was set to be 300 K, a commonly used temperature for modeling biomolecular systems.…”