Derivation of 13C chemical shift surfaces for the anomeric carbons of oligosaccharides and glycopeptides using ab initio methodology.

The dependence between the anomeric carbon chemical shift and the glycosidic bond (phi, psi) dihedral angles in oligosaccharide and glycopeptide model compounds was studied by Gauge-Including Atomic Orbital (GIAO) ab initio calculations. Complete chemical shift surfaces versus phi and psi for D-Glcp-D-Glcp disaccharides with (1-->1), (1-->2), (1-->3), and (1-->4) linkages in both alpha- and beta-configurations were computed using a 3-21G basis set, and scaled to reference results from calculations at the 6-311G** level of theory. Similar surfaces were obtained for GlcNAcThr and GlcNAcSer model glycopeptides in alpha- and beta-configurations, using in this case different conformations for the peptide moiety. The results obtained for both families of model compounds are discussed. We also present the determination of empirical formulas of the form 13Cdelta = f (phi,psi) obtained by fitting the raw ab initio data to trigonometric series expansions suitable for use in molecular mechanics and dynamics simulations. Our investigations are consistent with experimental observations and earlier calculations performed on smaller glycosidic bond models, and show the applicability of chemical shift surfaces in the study of the conformational behavior of oligosaccharides and glycopeptides.

Computational studies on HIV-1 protease inhibitors: influence of calculated inhibitor-enzyme binding affinities on the statistical quality of 3D-QSAR CoMFA models.

A theoretical study was performed on a set of 38 human immunodeficiency type 1 (HIV-1) protease inhibitors that are structurally similar to the AIDS drug Indinavir. Comparison between the computed binding energies and experimental activity data (pIC(50)) found a high degree of correlation (r(2)() = 0.82). Three-dimensional quantitative structure-activity relationship (3D-QSAR) models using comparative molecular field analysis (CoMFA) yielded predicted activities that were in excellent agreement with the corresponding experimentally determined values. Inclusion of the calculated enzyme-inhibitor binding energy as an additional descriptor in the CoMFA model yielded a significant improvement in the internal predictive ability of our model (q(2)() = 0.45 to q(2)() = 0.69). Separate CoMFA models were constructed to evaluate the influence of different alignment schemes (Atom Fit and Field Fit) and different partial atomic charge assignment schemes (Discover CVFF, Gasteiger-Marsili, and AM1-ESP) on the statistical quality of the models.

Evidence for a strong sulfur-aromatic interaction derived from crystallographic data.

We have uncovered new evidence for a significant interaction between divalent sulfur atoms and aromatic rings. Our study involves a statistical analysis of interatomic distances and other geometric descriptors derived from entries in the Cambridge Crystallographic Database (F. H. Allen and O. Kennard, Chem. Design Auto. News, 1993, Vol. 8, pp. 1 and 31-37). A set of descriptors was defined sufficient in number and type so as to elucidate completely the preferred geometry of interaction between six-membered aromatic carbon rings and divalent sulfurs for all crystal structures of nonmetal-bearing organic compounds present in the database. In order to test statistical significance, analogous probability distributions for the interaction of the moiety X-CH(2)-X with aromatic rings were computed, and taken a priori to correspond to the null hypothesis of no significant interaction. Tests of significance were carried our pairwise between probability distributions of sulfur-aromatic interaction descriptors and their CH(2)-aromatic analogues using the Smirnov-Kolmogorov nonparametric test (W. W. Daniel, Applied Nonparametric Statistics, Houghton-Mifflin: Boston, New York, 1978, pp. 276-286), and in all cases significance at the 99% confidence level or better was observed. Local maxima of the probability distributions were used to define a preferred geometry of interaction between the divalent sulfur moiety and the aromatic ring. Molecular mechanics studies were performed in an effort to better understand the physical basis of the interaction. This study confirms observations based on statistics of interaction of amino acids in protein crystal structures (R. S. Morgan, C. E. Tatsch, R. H. Gushard, J. M. McAdon, and P. K. Warme, International Journal of Peptide Protein Research, 1978, Vol. 11, pp. 209-217; R. S. Morgan and J. M. McAdon, International Journal of Peptide Protein Research, 1980, Vol. 15, pp. 177-180; K. S. C. Reid, P. F. Lindley, and J. M. Thornton, FEBS Letters, 1985, Vol. 190, pp. 209-213), as well as studies involving molecular mechanics (G. Nemethy and H. A. Scheraga, Biochemistry and Biophysics Research Communications, 1981, Vol. 98, pp. 482-487) and quantum chemical calculations (B. V. Cheney, M. W. Schulz, and J. Cheney, Biochimica Biophysica Acta, 1989, Vol. 996, pp.116-124; J. Pranata, Bioorganic Chemistry, 1997, Vol. 25, pp. 213-219)-all of which point to the possible importance of the sulfur-aromatic interaction. However, the preferred geometry of the interaction, as determined from our analysis of the small-molecule crystal data, differs significantly from that found by other approaches.

Comparison of ring current methods for use in molecular modeling refinement of NMR derived three-dimensional structures.

A comparison between three different methods commonly used to estimate ring current effects on chemical shifts is presented. Haigh-Mallion, Johnson-Bovey, and classical point-dipole approximations were used to estimate the ring current contribution to chemical shifts for protons in several proteins for which both detailed X-ray crystal structures and chemical shift assignments were available. For the classical point-dipole model, new proportionality constants were calculated by fitting to ring current estimations from both the quantum-mechanical Haigh-Mallion and semiclassical Johnson-Bovey methods and compared with the previously used point-dipole constant of Perkins and Dwek. Statistical analysis of the predictions obtained by all methods indicates that the point-dipole approximation parametrized against quantum-mechanical data is superior to the previously used classical model, comparable to Johnson-Bovey calculations, and slightly poorer than predictions from the Haigh-Mallion theory. The implementation of a pseudoenergy penalty term for use in structure refinement from chemical shift data based on the classical point-dipole model is described, and its usefulness in cases where other NMR information is limited is discussed with a specific example.

SMART: a solvent-accessible triangulated surface generator for molecular graphics and boundary element applications.

An algorithm is presented for generating a representation of the solvent-accessible molecular surface as a smooth triangulated manifold. The algorithm, called SMART (SMooth moleculAR surface Triangulator), divides the contact and reentrant portions of the solvent-accessible molecular surface into curvilinear three-sided elements. In contrast to the author's earlier implementation of this general approach [Zauhar, R.J. and Morgan, R.S., J. Comput. Chem., 11 (1990) 603], the SMART algorithm defines elements directly on the appropriate geometric surface types (rather than using interpolation over cubic elements), and has special features to handle highly distorted regions which often appear in deep crevices and internal cavities. While the method is designed for use with boundary element techniques in continuum electrostatics, it can also be applied to the accurate computation of molecular surface areas and volumes, and the generation of shaded surfaces for display with interactive computer graphics.

Europium(III) luminescence and tyrosine to terbium(III) energy-transfer studies of invertebrate (octopus) calmodulin.

The effects of minor differences in the amino acid sequences between a vertebrate (bovine testes) and an invertebrate (octopus) calmodulin on metal ion binding were investigated via laser-induced Eu3+ and Tb3+ luminescence. Amino acid substitutions at residues which are coordinated to the metal ion do not produce any detectable changes in the 7F0----5D0 excitation spectrum of the Eu3+ ion bound to octopus calmodulin relative to bovine testes calmodulin; only minor differences in the excited-state lifetime values in D2O solution are observed. The dissociation constants for Eu3+ (1.0 +/- 0.2 microM) and Tb3+ (5 +/- 1 microM) from the weak lanthanide binding sites (III and IV, numbered from the amino terminus) of octopus calmodulin were measured using luminescence techniques. Both values agree well with those reported previously for bovine testes calmodulin [Mulqueen, P. M., Tingey, J. M., & Horrocks, W. D., Jr. (1985) Biochemistry 24, 6639-6645]. The measured dissociation constant of Eu3+ bound in the tight lanthanide binding sites (I and II) is 6 +/- 2 nM for octopus calmodulin and 12 +/- 2 nM for bovine testes calmodulin. The distances between sites I and II (12.4 +/- 0.5 A) and sites III and IV (11.7 +/- 0.8 A) were determined from Förster-type energy transfer in D2O solutions of octopus calmodulin containing bound Eu3+ donor and Nd3+ acceptor ions. Förster theory parameters for nonradiative energy transfer between Tyr138 and Tb3+ ions bound at sites III and IV of octopus calmodulin were comprehensively evaluated, including a dynamics simulation of the orientation factor kappa 2. This theory is found to account quantitatively for the observed energy-transfer efficiency as evaluated from the observed sensitized Tb3+ emission.

A new method for computing the macromolecular electric potential.

A general methodology is developed for the rigorous computation of the electrostatic potential for a protein of arbitrary shape, assuming the presence of linear dielectric media. The theory proceeds by considering the distribution of induced polarization charge at the dielectric interface, rather than by attempting a direct solution of Poisson's equation (as in the finite-difference approach of Warwicker & Watson). The method is applied to a study of two-dimensional model proteins, where it is shown that the presence of a cleft is associated with a region of relatively high potential in the solvent medium. The results of a preliminary calculation in three dimensions for the protein lysozyme are also discussed; again, a region of enhanced potential is observed near the cleft at the active site. Our computational evidence supports the suggestion of Warwicker & Watson that clefts are associated with important electrostatic effects.