We investigated the translation of a protein through model nanopores using coarse-grained (CG) nonequilibrium molecular dynamics (NEMD) simulations and compared the mobilities with those obtained from previous coarse-grained equilibrium molecular dynamics model. We considered the effects of nanopore confinement and external force on the translation of streptavidin through nanopores of dimensions representative of experiments. As the nanopore radius approaches the protein hydrodynamic radius, r_{h}/r_{p}→1 (where r_{h} is the hydrodynamic radius of protein and r_{p} is the pore radius), the translation times are observed to increase by two orders of magnitude. The translation times are found to be in good agreement with the one-dimensional biased diffusion model. The results presented in this paper provide useful insights on nanopore designs intended to control the motion of biomolecules.

The rotational diffusion coefficient is an essential parameter in determining the mechanistic features of biomolecules in both crowded and confined environments. Understanding the influence of nanoconfinement on rotational diffusion is vital in conceptualizing dynamics of biomolecules (such as proteins) in nanopores. The control of the translational movement of biomolecules is practiced widely in nanopore experiments. However, the restrictions on the translational movement may affect other dynamic properties such as rotational diffusion. In this paper, we use a coarse-grained molecular dynamics approach to study the rotational dynamics of a sample protein under the influence of cylindrical nanopore confinement. Our simulation reveals a 2-fold reduction in magnitude from the bulk rotational diffusion coefficient value as the confinement radius reaches double the size of protein’s hydrodynamic radius. However, the changes in the rotational diffusion coefficient are relatively small compared to the changes in the translational diffusion coefficient. Interestingly, the rotational anisotropy also varies considerably when pore radii approach protein dimensions. Our simulations point out that the confinement effects cause the breakdown of small angular displacement theory when the pore radius is close to the protein hydrodynamic radius.

Combined kinetic theory-hydrodynamics treatment has been proven effective in the prediction of biomolecule dynamics, generally if a single biomolecule is present in the bulk solvent. But the validity of such a theory in many physiological conditions is controversial. In the present study, a sample protein surrounded by other large biomolecules is approximated as the protein in a cylindrical nanopore. The hydrodynamic radius of the protein is chosen as an indicator to check whether one of the widely used kinetic theoryhydrodynamics relation namely Stokes-Einstein-Debye relation, is genuine for confined conditions of the protein. It has been found that Stokes-Einstein-Debye relation cannot be satisfied by the protein if confinement dimensions are very close. The reason for the violation can be attributed to van der Waals interaction between pore and the protein.

No abstract

Viscosity variation of solvent in local regions near a solid surface, be it a biological surface of a protein or an engineered surface of a nanoconfinement, is a direct consequence of intermolecular interactions between the solid body and the solvent. The current coarse-grained molecular dynamics study takes advantage of this phenomenon to investigate the anomaly in a solvated protein’s rotational dynamics confined using a representative solid matrix. The concept of persistence time, the characteristic time of structural reordering in liquids, is used to compute the solvent’s local viscosity. With an increase in the degree of confinement, the confining matrix significantly influences the solvent molecule’s local viscosity present in the protein hydration layer through intermolecular interactions. This effect contributes to the enhanced drag force on protein motion, causing a reduction in the rotational diffusion coefficient. Simulation results suggest that the direct matrix-protein non-bonded interaction is responsible for the occasional jump and discontinuity in orientational motion when the protein is in very tight confinement.