Parameter Optimization Method in Multidimensional Umbrella Sampling

Umbrella sampling (US) is an effective method for calculating free-energy landscapes (FELs). However, the complexity of controlling the sampling positions complicates multidimensional FEL calculations. In this study, we proposed a method for controlling sampling by optimizing the US parameters. This method comprises the introduction of a target point and the optimization of the parameters to sample a window around this point. We approximated each window to normal distributions using an umbrella integration method and calculated the divergences between the window distributions and the state distributed at the target position by a variationally enhanced sampling method. Thus, the minimization of the divergence facilitated sampling around the target point, after which the parameters could be optimized on the fly while performing equilibration simulation. In practice, our method employs bias potentials with off-diagonal terms, ensuring a more efficient calculation of multidimensional FEL. Additionally, we developed an algorithm for determining the target point for automated FEL search; the algorithm samples in a specified direction while controlling the overlap of distributions. We performed three different FEL calculations as examples: (1) the calculation of the permeation of a water molecule through a lipid bilayer (one-dimensional FEL), (2) the calculation of the internal structural changes in alanine dipeptide in water (two-dimensional FEL), and (3) the calculation of the internal structural changes from a β-strand structure to an α-helix structure in alanine decapeptide (Ala10, 16-dimensional FEL). These results confirmed that our method could control the number of US windows and calculate the high-dimensional FELs that could not be evaluated by the conventional US method.