Computer simulation methods spanning several temporal and spatial scales are reviewed, focusing on their applications on the neuromuscular synapse. Quantum mechanics treats the enzymatic catalysis of neurotransmitters on the picometer scale. Molecular dynamics reveals conformational changes of the enzyme acetylcholinesterase for nanoseconds. Brownian dynamics follow the substrate molecule in its diffusion on the microsecond level. Methods such as finite elements describe the diffusion of neurotransmitters as a changing concentration continuum in the synapse. Promising directions for future research include integration of methods on several scales, and applying these methods to the acetylcholine receptor.
OVERVIEWWith the development of both theory and computer technology, simulation has become a viable way of investigating the behavior of molecules. A series of simulation methods are now available to address chemical questions at several time and length scales [1]. The success of these methods in "pure", theoretical chemistry has been matched by their applications in materials science.Biochemistry, the molecule-based study of life, also requires understanding on several scales, some of which are not readily reached by experiment. Encouragingly, it is possible to apply the simulation methods in biochemistry (Fig. 1, p. 298). Here, some of these endeavors-as applied on an example, the synapse in the nervous system-are discussed, after a short presentation of several simulation methods for different temporal and spatial scales.Quantum mechanics (QM) calculations have been firmly established as a rigorous methodology; their pioneers were awarded a Nobel Prize in chemistry in 1998 [2,3]. Combined with larger-grain methods such as molecular dynamics (see below), treatment of enzymatic catalysis to a degree comparable to that for small-molecule reactions is feasible.First applied to proteins more than 25 years ago by McCammon and coworkers [4], the method of Newtonian molecular dynamics (MD) has now entered the mainstream of biochemistry [5]. MD aptly reveals the macromolecular conformational changes on the nanosecond-level, filling a time-resolution gap (until recently) long left by experimental methods.The methods MD, Brownian dynamics (BD) [6], Monte Carlo (MC) simulation, and the finite elements (FE) [7] together offer a spectrum of tools, describing diffusion either in a discrete fashionmolecule-by-molecule-or as a continuum of concentration. Recently, other mesoscale and multiscale methods have emerged, with the promise of bringing up the scale of simulations to the cellular level. Ambitious projects aim at yet larger scales, but they lie beyond the scope of this review [8,9].