Polarizability and Orientation Dynamics of Small Proteins

Abstract: Proteins often carry an intrinsic electric dipole moment, which can interact with external electric fields and cause protein motion. Previous research has found that the orientation of small proteins in gas phase can be controlled in a static electric field. This effect is hoped to benefit applications such as single-particle imaging, and possibly other techniques involving proteins in electric fields. With the purpose of improving our understanding and modeling of protein orientation, this project investigated the scarcely explored quantum mechanical aspects of the process, namely the polarizability. Ground-state electronic structure simulations of three small model proteins, ubiquitin, Trp-cage and lysozyme, under the influence of electric fields were performed in vacuum. The electric dipole moments of the proteins were extracted from simulations with an applied electric field of strength 1 V/nm for varying angles, with respect to a body fixed reference frame. A Python program was written to analyze and visualize the results. The results point to a connection between the polarizability and the structure of the proteins, as well as size. Next a 3D rigid rotor model was developed using Mathematica in order to study the orientation dynamics classically in a simplified and time efficient way, with the possibility of including the previous quantum results. A comparison between a simulation of ubiquitin with and without polarizability concluded that the polarizability seems to have a damping effect on the orientation dynamics, at least for the initial conditions tested in this study. Further research is necessary to validate the model and perform statistical analysis of many simulations with varying initial conditions.