Chemical shifts in proteins: a shielding trajectory analysis of the fluorine nuclear magnetic resonance spectrum of the Escherichia coli galactose binding protein using a multipole shielding polarizability-local reaction field-molecular dynamics approach

Abstract

We report the first successful calculation of the fluorine nuclear magnetic resonance spectrum of a protein, the galactose binding protein from Escherichia coli, labeled with ( 5-'9F) tryptophan. Our results indicate that the experimental l9F chemical shifts are dominated by weak, or long-range, electrical interactions, which can be calculated by using the responses of the tensor elements to the uniform field components (&~d/dE,) and the nonuniform or gradient terms (t3uap/W,,), together with the average values of the fields (Ex) and field gradients (V,,) obtained from molecular dynamics (MD) trajectories. A series of shielding Au(Ex,Vxx,Vyy,Vzz)f(7), are obtained, and the mean values, A;, correlate well with the actual shift pattern and overall shift range observed experimentally (Luck, L. A.; Falke, J. J. Biochemistry 1991,30,42484256). The computed 19F NMR of the pentapeptide Gly- Gly-(5-F)Trp-Gly-Gly in HzO is close to, but somewhat more shielded than, that of the denatured protein. Almost all computed 19F chemical shifts are upfield of the field-free value, in accord with the results of ab initio calculations. The chemical shifts calculated are sensitive to the charge chosen for F in the LRF-MD trajectories, and best agreement with experiment is obtained with qF = -0.25. Calculations on the Salmonella typhimurium galactose binding protein yield extremely similar chemical shift spectra, consistent with the -7% difference in amino acid composition. The uniform field components make the largest contributions to the patterns observed, presumably because the gradient terms fall off more rapidly with distance. The exposed residue, Trp 284, has the largest amplitude of fluctuation associated with its trajectory, possibly due to the rapid and random movement of neighboring water molecules. Trp 284 is highly shielded due to interaction with water, although other buried residues (e.g. Trp 133) may also be highly shielded, due to electric field effects within the protein. Our results imply that van der Waals interactions do not play a major role for fluorine nonequivalence in proteins, since the experimental results can be reproduced by using solely the computed field and field gradient terms. The ability to predict protein NMR patterns and ranges offers promise for structural analysis and also provides a way of validating different methods of computing protein electrostatics. In this respect, it is instructive to note that the charge fields from ionized surface groups are found to be largely shielded by the dielectric response of the solvent and the protein.