Brownian dynamics simulations of biomolecular diffusional association processes

Abstract

Brownian dynamics (BD) is a computational method to simulate molecular diffusion processes. Although the BD method has been developed over several decades and is well established, new methodological developments are improving its accuracy, widening its scope, and increasing its application. In biological applications, BD is used to investigate the diffusive behavior of molecules subject to forces due to intermolecular interactions or interactions with material surfaces. BD can be used to compute rate constants for diffusional association, generate structures of encounter complexes for molecular binding partners, and examine the transport properties of geometrically complex molecules. Often, a series of simulations is performed, for example, for different protein mutants or environmental conditions, so that the effects of the changes on diffusional properties can be estimated. While biomolecules are commonly described at atomic resolution and internal molecular motions are typically neglected, coarse‐graining and the treatment of conformational flexibility are increasingly employed. Software packages for BD simulations of biomolecules are growing in capabilities, with several new packages providing novel features that expand the range of questions that can be addressed. These advances, when used in concert with experiment or other simulation methods, such as molecular dynamics, open new opportunities for application to biochemical and biological systems. Here, we review some of the latest developments in the theory, methods, software, and applications of BD simulations to study biomolecular diffusional association processes and provide a perspective on their future use and application to outstanding challenges in biology, bioengineering, and biomedicine. This article is categorized under: Structure and Mechanism > Computational Biochemistry and Biophysics Molecular and Statistical Mechanics > Molecular Dynamics and Monte‐Carlo Methods Software > Simulation Methods

Transient and elusive intermediate states in self‐assembly processes: An overview

“…[44] The methodology rests upon a detailed treatment of the stochastic motion of particles within a solvent, where this motion encompasses both Brownian motion and hydrodynamic interactions. [45] This treatment is achieved via the Langevin equations, which capture the translational and rotational behavior of the particles and provide a means to bridge between thermal motion induced by solvent molecule collisions and the translational motion of the particles. [46] The applicability of this methodology is a key aspect of BD simulations, particularly for complex particle types, such as colloids that exhibit surface charges or anisotropic shapes.…”

“…With the development of BD simulation, researchers have been able to model the self‐assembly of nanoparticle and colloidal entities with tremendous precision [44] . The methodology rests upon a detailed treatment of the stochastic motion of particles within a solvent, where this motion encompasses both Brownian motion and hydrodynamic interactions [45] . This treatment is achieved via the Langevin equations, which capture the translational and rotational behavior of the particles and provide a means to bridge between thermal motion induced by solvent molecule collisions and the translational motion of the particles [46] .…”

An Introductory Tutorial to the SEEKR2 (Simulation Enabled Estimation of Kinetic Rates V. 2) Multiscale Milestoning Software [Article v1.0]

“…Brownian dynamics describes motion at a mesoscopic level, in which models of large molecules are usually simplified to smaller collections of rigid bodies, and the solvent is treated using continuum theories [19] instead of explicitly with individual molecules. A recent review summarizes theoretical details, common methods and software, and some of the most recent applications of BD [20]. Considering Newton's equations of motion and assuming a separation between the time scales of the macromolecules and the solvent molecules, one can derive the following stochastic differential equation of motion [21]:…”

Editorial: Mechanisms, thermodynamics and kinetics of ligand binding revealed from molecular simulations and machine learning

“…These simulations can then be analyzed to compute both thermodynamic free energies and kinetic rates of ligand binding (Pang and Zhou, 2017;Tang et al, 2017;Nunes-Alves et al, 2020;Wang et al, 2022). In addition, Brownian dynamics simulations have been very efficient in generating a large number of ligand binding trajectories and estimating the binding kinetic rates (Huber and McCammon, 2019;Muñiz-Chicharro et al, 2022). Finally, emerging machine learning techniques have greatly enhanced molecular simulations and facilitated analysis of the simulation trajectories (Glielmo et al, 2021).This Research Topic is focused on studies of the pathways, mechanisms, free energies and kinetics of ligand binding to target receptors.…”