No abstract

A method is presented for the rigorous computation of the electric potential of molecules of arbitrary shape, under the assumption of continuous linear dielectric media. The computational technique involves finding the distribution of indued polarization charge on the molecular surface, and proceeds by an application of the method of boundary elements. The surface, which separates the molecular interior (of low dielectric constant) from the highly polar solvent, is given a piece-wise analytic representation as a collection of curvilinear elements. Given a set of internal fixed charges, the distribution of polarizationcharge is found as a continuous function over the surface elements, and the electric potential (including all polarization effects) is then easily computed at any point. The method is applied to a spherical interface, and to several small molecules of biological interest, including a hexapeptide. The resulting potentials show good convergence in all cases. The future application of the method to macromolecules is discussed.

A unifying principle of rational drug design is the use of either shape similarity or complementarity to identify compounds expected to be active against a given target. Shape similarity is the underlying foundation of ligand-based methods, which seek compounds with structure similar to known actives, while shape complementarity is the basis of most receptor-based design, where the goal is to identify compounds complementary in shape to a given receptor. These approaches can be extended to include molecular descriptors in addition to shape, such as lipophilicity or electrostatic potential. Here we introduce a new technique, which we call shape signatures, for describing the shape of ligand molecules and of receptor sites. The method uses a technique akin to ray-tracing to explore the volume enclosed by a ligand molecule, or the volume exterior to the active site of a protein. Probability distributions are derived from the ray-trace, and can be based solely on the geometry of the reflecting ray, or may include joint dependence on properties, such as the molecular electrostatic potential, computed over the surface. Our shape signatures are just these probability distributions, stored as histograms. They converge rapidly with the length of the ray-trace, are independent of molecular orientation, and can be compared quickly using simple metrics. Shape signatures can be used to test for both shape similarity between compounds and for shape complementarity between compounds and receptors and thus can be applied to problems in both ligand- and receptor-based molecular design. We present results for comparisons between small molecules of biological interest and the NCI Database using shape signatures under two different metrics. Our results show that the method can reliably extract compounds of shape (and polarity) similar to the query molecules. We also present initial results for a receptor-based strategy using shape signatures, with application to the design of new inhibitors predicted to be active against HIV protease.

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.