Asking biological questions of physical systems: The device approach to emergent properties

H-Bonds in Crambin: Coherence in an α-Helix

We applied coherence analysis-used by engineers to identify linear interactions in stochastic systems-to molecular dynamics simulations of crambin, a thionin storage protein found in Abyssinian cabbage. A key advantage of coherence over other analyses is that it is robust, independent of the properties, or even the existence of probability distributions often relied on in statistical mechanics. For frequencies between 0.391 and 5.08 GHz (corresponding reciprocally to times of 2.56 and 0.197 ns), the displacements of oxygen and nitrogen atoms across α-helix H-bonds are strongly correlated, with a coherence greater than 0.9; the secondary structure causes the H-bonds to effectively act as a spring. Similar coherence behavior is observed for covalent bonds and other noncovalent interactions including H-bonds in β-sheets and salt bridges. In contrast, arbitrary pairs of atoms that are physically distant have uncorrelated motions and negligible coherence. These results suggest that coherence may be used to objectively identify atomic interactions without subjective thresholds such as H-bond lengths angles and angles. Strong coherence is also observed between the average position of adjacent leaves (groups of atoms) in an α-helix, suggesting that the harmonic analysis of classical molecular dynamics can successfully describe the propagation of allosteric interactions through the structure.

Setting Boundaries for Statistical Mechanics

Statistical mechanics has grown without bounds in space. Statistical mechanics of noninteracting point particles in an unbounded perfect gas is widely used to describe liquids like concentrated salt solutions of life and electrochemical technology, including batteries. Liquids are filled with interacting molecules. A perfect gas is a poor model of a liquid. Statistical mechanics without spatial bounds is impossible as well as imperfect, if molecules interact as charged particles, as nearly all atoms do. The behavior of charged particles is not defined until boundary structures and values are defined because charges are governed by Maxwell's partial differential equations. Partial differential equations require boundary structures and conditions. Boundary conditions cannot be defined uniquely 'at infinity' because the limiting process that defines 'infinity' includes such a wide variety of structures and behaviors, from elongated ellipses to circles, from light waves that never decay, to dipolar fields that decay steeply, to Coulomb fields that hardly decay at all. Boundaries and boundary conditions needed to describe matter are not prominent in classical statistical mechanics. Statistical mechanics of bounded systems is described in the EnVarA system of variational mechanics developed by Chun Liu, more than anyone else. EnVarA treatment does not yet include Maxwell equations.

Multiscale modeling shows that dielectric differences make NaV channels faster than KV channels

The generation of action potentials in excitable cells requires different activation kinetics of voltage-gated Na (Na_V) and K (K_V) channels. Na_V channels activate much faster and allow the initial Na⁺ influx that generates the depolarizing phase and propagates the signal. Recent experimental results suggest that the molecular basis for this kinetic difference is an amino acid side chain located in the gating pore of the voltage sensor domain, which is a highly conserved isoleucine in K_V channels but an equally highly conserved threonine in Na_V channels. Mutagenesis suggests that the hydrophobicity of this side chain in Shaker K_V channels regulates the energetic barrier that gating charges cross as they move through the gating pore and control the rate of channel opening. We use a multiscale modeling approach to test this hypothesis. We use high-resolution molecular dynamics to study the effect of the mutation on polarization charge within the gating pore. We then incorporate these results in a lower-resolution model of voltage gating to predict the effect of the mutation on the movement of gating charges. The predictions of our hierarchical model are fully consistent with the tested hypothesis, thus suggesting that the faster activation kinetics of Na_V channels comes from a stronger dielectric polarization by threonine (Na_V channel) produced as the first gating charge enters the gating pore compared with isoleucine (K_V channel).

Modeling the Device Behavior of Biological and Synthetic Nanopores with Reduced Models

Biological ion channels and synthetic nanopores are responsible for passive transport of ions through a membrane between two compartments. Modeling these ionic currents is especially amenable to reduced models because the device functions of these pores, the relation of input parameters (e.g., applied voltage, bath concentrations) and output parameters (e.g., current, rectification, selectivity), are well defined. Reduced models focus on the physics that produces the device functions (i.e., the physics of how inputs become outputs) rather than the atomic/molecular-scale physics inside the pore. Here, we propose four rules of thumb for constructing good reduced models of ion channels and nanopores. They are about (1) the importance of the axial concentration profiles, (2) the importance of the pore charges, (3) choosing the right explicit degrees of freedom, and (4) creating the proper response functions. We provide examples for how each rule of thumb helps in creating a reduced model of device behavior.

Application of a bipolar nanopore as a sensor: rectification as an additional device function

We model and simulate a nanopore sensor that selectively binds analyte ions. This binding leads to the modulation of the local concentrations of the ions of the background electrolyte (KCl), and, thus, to the modulation of the ionic current flowing through the pore. The nanopore's wall has a bipolar charge pattern with a larger positive buffer region determining the anions as the main charge carriers and a smaller negative binding region containing binding sites. This charge pattern proved to be an appropriate one as shown by a previous comparative study of varying charge patterns (Mádai et al. J. Mol. Liq., 2019, 283, 391-398.). Binding of the positive analyte ions attracts more anions in the pore thus increasing the current. The asymmetric nature of the pore results in an additional device function, rectification. Our model, therefore, is a dual response device. Using a reduced model of the nanopore studied by a hybrid computer simulation method (Local Equilibrium Monte Carlo coupled with the Nernst-Planck equation) we show that we can create a sensor whose underlying mechanisms are based on the changes in the local electric field as a response to changing thermodynamic conditions. The change in the electric field results in changes in the local ionic concentrations (depletion zones), and, thus, changes in ionic currents.

Multiscale analysis of the effect of surface charge pattern on a nanopore's rectification and selectivity properties: From all-atom model to Poisson-Nernst-Planck

We report a multiscale modeling study for charged cylindrical nanopores using three modeling levels that include (1) an all-atom explicit-water model studied with molecular dynamics, and reduced models with implicit water containing (2) hard-sphere ions studied with the Local Equilibrium Monte Carlo simulation method (computing ionic correlations accurately), and (3) point ions studied with Poisson-Nernst-Planck theory (mean-field approximation). We show that reduced models are able to reproduce device functions (rectification and selectivity) for a wide variety of charge patterns, that is, reduced models are useful in understanding the mesoscale physics of the device (i.e., how the current is produced). We also analyze the relationship of the reduced implicit-water models with the explicit-water model and show that diffusion coefficients in the reduced models can be used as adjustable parameters with which the results of the explicit- and implicit-water models can be related. We find that the values of the diffusion coefficients are sensitive to the net charge of the pore but are relatively transferable to different voltages and charge patterns with the same total charge.