Stagnation of flow in protein cavities by boundary element microhydrodynamics.
In this work, we apply the boundary element method to describe the fluid velocity profiles in pockets in protein surfaces that are crucial to their function as enzymes. First, we study a simplified model, that of a dimpled sphere, in order to properly interpret the behavior in more complex surfaces such as proteins. In that case, we are able to observe the difference between an unphysical sharp edge for the dimple and a smooth edge. The sharp edge produces extra dissipation in the fluid, accounting for much more friction for all types of body motions. We were able to observe the direct correlation of the stagnation depth with the depth of the dimple in this simple case, allowing us to interpret this feature in a similar fashion for proteins. We have found that the fluid in the protein pockets translates with the body, irrespective of the direction body motion, for a distance comparable to the size of the pocket, and that such stagnation volumes are larger for motions parallel to the pocket axis. Outside of these pockets, the fluid velocity profile decays to that of the surrounding fluid far away from the protein (taken to be zero in our case, for convenience), as the Oseen tensor requires. We have also found that there is weak local motion of fluid inside of the pockets, with velocities about 1% of those of the body. This study suggests that there may be a role for the hydrodynamics of solvent inside of pockets for the transport of substrates to protein active sites. If solvent is effectively stagnant inside of a pocket, then transport must occur by diffusion near the pocket surface even if the fluid around the protein is stirred. The weak local motions inside of the pocket may also be relevant in this transport process, but these may be easily overwhelmed by any electrostatic interactions that are likely present at active sites.
Aragon
Sergio R
SR
Department of Chemistry and Biochemistry, San Francisco State University, San Francisco, California, USA. aragons@sfsu.edu
Hahn
David K
DK
D000418:Albumins / Q000737:chemistry; D000465:Algorithms; D000818:Animals; D001703:Biophysics / Q000379:methods*; D020134:Catalytic Domain; D038101:Insulator Elements; D015233:Models, Statistical; D008968:Molecular Conformation; D009113:Muramidase / Q000737:chemistry; D009211:Myoglobin / Q000737:chemistry; D011487:Protein Conformation; D011506:Proteins / Q000737:chemistry*; D040901:Proteomics / Q000379:methods; D012997:Solvents / Q000737:chemistry; D013314:Stress, Mechanical
D000418:Albumins; D009211:Myoglobin; D011506:Proteins; D012997:Solvents; D009113:Muramidase
J Chem Theory Comput. 2006 Sep;2(5):1416-28
26626849
J Comput Chem. 2004 Jul 15;25(9):1191-205
15116362
Biopolymers. 1986 Apr;25(4):627-37
3708107
Proc Natl Acad Sci U S A. 1993 May 1;90(9):3835-9
8483901
J Mol Graph. 1993 Jun;11(2):139-41
8347567
Science. 1983 Aug 19;221(4612):709-13
6879170
Biophys J. 2006 Sep 1;91(5):1591-603
16714342
Biophys J. 1995 Jun;68(6):2261-70
7647233