On an exact starting expression for macromolecular hydrodynamic models

Intrinsic Viscosity of Proteins and Platonic Solids by Boundary Element Methods

The boundary element (BE) method is used to implement a very precise computation of the intrinsic viscosity for rigid molecules of arbitrary shape. The formulation, included in our program BEST, is tested against the analytical Simha formula for ellipsoids of revolution, and the results are essentially numerically exact. Previously unavailable, very precise results for a series of Platonic solids are also presented. The formulation includes the optional determination of the center of viscosity; however, for globular proteins, the difference compared to the computation based on the centroid is insignificant. The main application is to a series of 30 proteins ranging in molecular weight from 12 to 465 kD. The computation starts from the crystal structure as obtained from the Protein Data Bank, and a hydration thickness of 1.1 Å obtained in previous work with BEST was used. The results (extrapolated to an infinite number of triangular boundary elements) for the proteins are separated into two groups:  monomeric and multimeric proteins. The agreement with experimental measurements of the intrinsic viscosity in the case of monomeric proteins is excellent and within experimental error of 5%, demonstrating that the solution and crystal structure are hydrodynamically equivalent. However, for some multimeric proteins, we observe strong systematic deviations around -20%, which we interpret as a systematic deviation of the solution structure from the crystal structure. A possible description of the structural change is deduced by using simple ellipsoid model parameters. A method to obtain intrinsic viscosity values for proteins to 1-2% accuracy (better than experimental error) on the basis of a single BE computation (avoiding the need for an extrapolation on the number of surface triangles) is also presented.

Hydrodynamic effects on the relative rotational velocity of associating proteins

Hydrodynamic steering effects on the barnase-barstar association were studied through the analysis of the relative rotational velocity of the proteins. We considered the two proteins approaching each other in response to their electrostatic attraction and employed a method that accounts for the long-range and many-body character of the hydrodynamic interactions, as well as the complicated shapes of the proteins. Hydrodynamic steering effects were clearly seen when attractive forces were applied to the geometric centers of the proteins (resulting in zero torques) and the attraction acted along the line that connects centers of geometry of proteins in their crystallographic complex. When we rotated barstar relative to barnase around this line by an angle in the range from -90° to 60°, the rotational velocity arising solely from hydrodynamic interactions restored the orientation of the proteins in the crystal structure. However, because, in reality, both electrostatic forces and torques act on the proteins and these forces and torques depend on the protein-protein distance and the relative orientation of the binding partners, we also investigated more realistic situations employing continuum electrostatics calculations based on atomistic protein models. Overall, we conclude that hydrodynamic interactions aid barnase and barstar in assuming a proper relative orientation upon complex formation.

Recent advances in macromolecular hydrodynamic modeling

The modern implementation of the boundary element method [23] has ushered unprecedented accuracy and precision for the solution of the Stokes equations of hydrodynamics with stick boundary conditions. This article begins by reviewing computations with the program BEST of smooth surface objects such as ellipsoids, the dumbbell, and cylinders that demonstrate that the numerical solution of the integral equation formulation of hydrodynamics yields very high precision and accuracy. When BEST is used for macromolecular computations, the limiting factor becomes the definition of the molecular hydrodynamic surface and the implied effective solvation of the molecular surface. Studies on 49 different proteins, ranging in molecular weight from 9 to over 400kDa, have shown that a model using a 1.1Å thick hydration layer describes all protein transport properties very well for the overwhelming majority of them. In addition, this data implies that the crystal structure is an excellent representation of the average solution structure for most of them. In order to investigate the origin of a handful of significant discrepancies in some multimeric proteins (about -20% observed in the intrinsic viscosity), the technique of Molecular Dynamics simulation (MD) has been incorporated into the research program. A preliminary study of dimeric α-chymotrypsin using approximate implicit water MD is presented. In addition I describe the successful validation of modern protein force fields, ff03 and ff99SB, for the accurate computation of solution structure in explicit water simulation by comparison of trajectory ensemble average computed transport properties with experimental measurements. This work includes small proteins such as lysozyme, ribonuclease and ubiquitin using trajectories around 10ns duration. We have also studied a 150kDa flexible monoclonal IgG antibody, Trastuzumab, with multiple independent trajectories encompassing over 320ns of simulation. The close agreement within experimental error of the computed and measured properties allows us to conclude that MD does produce structures typical of those in solution, and that flexible molecules can be properly described using the method of ensemble averaging over a trajectory. We review similar work on the study of a transfer RNA molecule and DNA oligomers that demonstrate that within 3% a simple uniform hydration model 1.1Å thick provides agreement with experiment for these nucleic acids. In the case of linear oligomers, the precision can be improved close to 1% by a non-uniform hydration model that hydrates mainly in the DNA grooves, in agreement with high resolution X-ray diffraction. We conclude with a vista on planned improvements for the BEST program to decrease its memory requirements and increase its speed without sacrificing accuracy.

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.

Precise boundary element computation of protein transport properties: Diffusion tensors, specific volume, and hydration

A precise boundary element method for the computation of hydrodynamic properties has been applied to the study of a large suite of 41 soluble proteins ranging from 6.5 to 377 kDa in molecular mass. A hydrodynamic model consisting of a rigid protein excluded volume, obtained from crystallographic coordinates, surrounded by a uniform hydration thickness has been found to yield properties in excellent agreement with experiment. The hydration thickness was determined to be delta = 1.1 +/- 0.1 A. Using this value, standard deviations from experimental measurements are: 2% for the specific volume; 2% for the translational diffusion coefficient, and 6% for the rotational diffusion coefficient. These deviations are comparable to experimental errors in these properties. The precision of the boundary element method allows the unified description of all of these properties with a single hydration parameter, thus far not achieved with other methods. An approximate method for computing transport properties with a statistical precision of 1% or better (compared to 0.1-0.2% for the full computation) is also presented. We have also estimated the total amount of hydration water with a typical -9% deviation from experiment in the case of monomeric proteins. Both the water of hydration and the more precise translational diffusion data hint that some multimeric proteins may not have the same solution structure as that in the crystal because the deviations are systematic and larger than in the monomeric case. On the other hand, the data for monomeric proteins conclusively show that there is no difference in the protein structure going from the crystal into solution.

Transient electric birefringence of wormlike macromolecules in electric fields of arbitrary strength: A computer simulation study

We have studied the birefringence decay of linear models of macromolecules for two different types of flexibility, the broken-rod chain and the wormlike chain, using a computer simulation of a transient electric birefringence experiment. We have paid particular attention to the influence of the intensity of the orienting field, including two orienting mechanisms, the induced dipole, and the permanent dipole. We have compared wormlike and broken-rod models of the same radius of gyration, finding that they present a different decay curve under the influence of the same intensity of the field. We have seen that these differences are due to the faster relaxation times (smaller in the wormlike chain model) and amplitudes, because, regardless of the type of flexibility, the overall size of a molecule (measured by the radius of gyration) essentially determines the longest relaxation time. We have also analyzed how the relaxation process is affected by the degree of flexibility, the orientation mechanisms, and the intensity of the field. Studying a different aspect, we have paid attention to the deformation of a molecule in a transient electric birefringence experiment as a source of information. In this work we have developed equations to characterize this deformation in terms of one of the components of the gyration tensor, if a dynamic light scattering experiment under the influence of an electric field could be performed. To develop this work we have simulated the Brownian dynamics of the different models, relaxing after the removal of an orienting external electric field of arbitrary strength. A comparison with other methods such a the rigid body treatment or the correlation analysis of Brownian trajectories has also been included. We have seen that differences between the two Brownian dynamics methods are small and that the rigid-body treatment is only an acceptable approximation to obtain the longest relaxation time.

A precise boundary element method for macromolecular transport properties

A very precise boundary element numerical solution of the exact formulation of the hydrodynamic resistance problem with stick boundary conditions is presented. BEST, the Fortran 77 program developed for this purpose, computes the full transport tensors in the center of resistance or the center of diffusion for an arbitrarily shaped rigid body, including rotation-translation coupling. The input for this program is a triangulation of the solvent-defined surface of the molecule of interest, given by Connolly's MSROLL or other suitable triangulator. The triangulation is prepared for BEST by COALESCE, a program that allows user control over the quality and number of triangles to describe the surface. High numerical precision is assured by effectively exact integration of the Oseen tensor over triangular surface elements, and by scaling the hydrodynamic computation to the precise surface area of the molecule. Efficiency of computation is achieved by the use of public domain LAPACK routines that call BLAS Level 3 hardware-optimized subroutines available for most processors. A protein computation can be done in less than 10 min of CPU time in a modern Pentium IV processor. The present work includes a complete analysis of the sources of error in the numerical work and techniques to eliminate these errors. The operation of BEST is illustrated with applications to ellipsoids of revolution, and Lysozyme, a small protein. The typical numerical accuracy achieved is 0.05% compared to analytical theory. The numerical precision for a protein is better than 1%, much better than experimental errors in these quantities, and more than 10 times better than traditional bead-based methods.

Predicting protein diffusion coefficients

Diffusion coefficients for proteins in water are predicted. The numerical method developed is general enough to be applied to a wide range of protein surface shapes, from rodlike to globular. Results are presented for lysozyme and tobacco mosaic virus, and they are compared with actual data and with predictions made by less general methods.