Langevin dynamics simulation of protein dynamics in nanopores at microsecond timescales

With rapid advancement in the fields of nanopore analysis of protein, it has become imperative to develop modeling framework for understanding the protein dynamics in nanopores. Such modeling framework should include the effects of electro-osmosis, as it plays significant role during protein translocation in confinement. Currently, the molecular dynamics simulations that include the hydrodynamic effects are limited to a timescale of few 100 ns. These simulations give insight about important events like protein unfolding which occurs in this timescale. But many electrophoresis experiments are limited by a detector resolution of approximately 2.5 microseconds. Analytical theory has been used to interpret protein dynamics at such large timescale. There is a need for molecular modeling of more complex environment and protein shapes which cannot be accounted for by analytical theory. We have developed a framework to study globular protein dynamics in nanopores by using langevin dynamics on a rigid body model of protein and the hydrodynamics is accounted by analytical theory for simple cylindrical nanopore geometry. This framework has been applied to study the dynamics of Ubiquitin translocation in SiNx nanopore by Nir et al. They have reported 7 times decrease in average dwell time of the protein inside the nanopore in response to a small change in pH from 7.0 to 7.2 and the modification of protein charge was attributed for such drastic change. Closer examination using our simulation revealed that the electro-osmotic effects originating due to very small change in the surface electrostatic potential of the nanopore could lead to such a drastic change in protein dynamics.