Multiple Scales in the Simulation of Ion Channels and Proteins

Bob Eisenberg

DOI: https://doi.org/10.1021/jp106760t

2019-03-05

Abstract:Computation of biological processes creates great promise for everyday life and great challenges for physical scientists. Simulations of molecular dynamics appeal to biologists as a natural extension of structural biology. Once biologists see a structure, they want to see it move. Molecular biology shows that a few atoms, often messenger ions like Ca$^{2+}$, control biological function on the scale of cells, sometimes organisms. Enormously concentrated ions (~20 M) in protein channels and enzymes can control important characteristics of living systems, just as highly concentrated ions near electrodes control important characteristics of electrochemical systems. The scale differences needed to simulate all the atoms of biological cells are $10^7$ in linear dimension, $10^{21}$ in three dimensions, $10^{9}$ in resolution, $10^{11}$ in time, and $10^{13}$ in particle number (to deal with concentrations of Ca$^{2+}$). $\mathbf{These}$ $\mathbf{scales}$ $\mathbf{must}$ $\mathbf{be}$ $\mathbf{dealt}$ $\mathbf{with}$ $\mathbf{simultaneously}$ if the simulation is to deal with most biological functions. We suggest a computational approach using explicit multiscale analysis instead of implicit simulation of all scales. Successful computation of ions concentrated in special places may be a significant step to understanding the defining characteristics of biological and electrochemical systems.

Biomolecules

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the multi - scale challenges faced by molecular dynamics simulations when dealing with ion channels and proteins. Specifically, the author discusses the huge differences on multiple scales such as space, time and concentration when simulating biological systems, and these differences make traditional all - atom simulations almost impossible to achieve. For example, from the scale of a single atom to the scale of cells, tissues and even the whole organism, the spatial scale spans multiple orders of magnitude; the time scale also has a huge gap from atomic vibrations at the femtosecond (fs) level to biological functions at the second level. In addition, the range of variation in ion concentration is also very wide, from the micromolar level to the hundreds of millimolar level.

To solve these problems, the author proposes an explicit multi - scale analysis method. This method is based on the energy variational principle (EnVarA) proposed by Chun Liu and can automatically handle multi - component and multi - scale interactions without additional unknown parameters. Through this method, highly concentrated ion solutions can be simulated more effectively, which is crucial for understanding the characteristics of biological and electrochemical systems. This method not only helps to simulate the functions of ion channels, but also promotes in - depth understanding of biomolecules and their functions.