Bimolecular reaction simulation using Weighted Ensemble Brownian dynamics and the University of Houston Brownian Dynamics program.
We discuss here the implementation of the Weighted Ensemble Brownian (WEB) dynamics algorithm of Huber and Kim in the University of Houston Brownian Dynamics (UHBD) suite of programs and its application to bimolecular association problems. WEB dynamics is a biased Brownian dynamics (BD) algorithm that is more efficient than the standard Northrup-Allison-McCammon (NAM) method in cases where reaction events are infrequent because of intervening free energy barriers. Test cases reported here include the Smoluchowski rate for association of spheres, the association of the enzyme copper-zinc superoxide dismutase with superoxide anion, and the binding of the superpotent sweetener N-(p-cyanophenyl)-N'-(diphenylmethyl)-guanidinium acetic acid to a monoclonal antibody fragment, NC6.8. Our results show that the WEB dynamics algorithm is a superior simulation method for enzyme-substrate reaction encounters with large free energy barriers.