Nanopore sensing is attracting the attention of a large and varied scientific community. One of the main issues in nanopore sensing is how to associate the measured current signals to specific features of the molecule under investigation. This is particularly relevant when the translocating molecule is a protein and the pore is sufficiently narrow to necessarily involve unfolding of the translocating protein. Recent experimental results characterized the cotranslocational unfolding of Thioredoxin (Trx) passing through an α-hemolisin pore, providing evidence for the existence of a multistep process. In this study we report the results of all-atom molecular dynamics simulations of the same system. Our data indicate that Trx translocation involves two main barriers. The first one is an unfolding barrier associated with a translocation intermediate where the N-terminal region of Trx is stuck at the pore entrance in a conformation that strongly resembles the native one. After the abrupt unfolding of the N-terminal region, the Trx enters the α-hemolisin vestibule. During this stage, the constriction is occupied not only by the translocating residue but also by a hairpin-like structure forming a tangle in the constriction. The second barrier is associated with the disentangling of this region.

In nanopore sensing experiments, the properties of molecules are probed by the variation of ionic currents flowing through the nanopore. In this context, the electronic properties and the single-layer thickness of graphene constitute a major advantage for molecule characterization. Here we analyze the translocation pathway of the thioredoxin protein across a graphene nanopore, and the related ionic currents, by integrating two nonequilibrium molecular dynamics methods with a bioinformatic structural analysis. To obtain a qualitative picture of the translocation process and to identify salient features we performed unsupervised structural clustering on translocation conformations. This allowed us to identify some specific and robust translocation intermediates, characterized by significantly different ionic current flows. We found that the ion current strictly anticorrelates with the amount of pore occupancy by thioredoxin residues, providing a putative explanation of the multilevel current scenario observed in recently published translocation experiments.

No abstract

Molecular dynamics simulations allowed to characterize the co-translocational unfolding pathway of ubiquitin through a biological nanopore.