In this article a new method is described to superimpose molecules using a crystallographic Fourier transform approach. Superimposed molecules, among Ž . other purposes, serve as a basis for three-dimensional 3D QSAR analyses in drug design and therefore an objective and reproducible method of molecule alignment is of major importance. Fourier data are generated for hypothetical crystals of cubic symmetry for the compounds under consideration. A Patterson-density-based similarity index is used to optimize rotational alignment of the molecules. After optimization of rotational orientation, an electron density derived similarity index is used to further optimize overlap of electron density as a function of translation of the molecules. Both similarity indices are maximized by a simple optimization routine, thus enabling automated superposition. The use of Fourier space offers several advantages. First, rotational and translational parameters can be optimized separately, thus providing a small parameter space. Second, a limited number of data already provide an Ž . adequate, continuous description of the electron or Patterson density distribution. Third, crystallography provides simple methods to calculate the Ž . Fourier transforms that are needed. The resolution of the Patterson electron * Author to whom all correspondence should be addressed at Utrecht University, Sorbonnelaan 16, 3584 Utrecht, The Netherlands. ( ) Journal of Computational Chemistry, Vol. 18, No. 5, 638]645 1997 Q 1997 by John Wiley & Sons CCC 0192-8651 / 97 / 050638-08 SUPERPOSITION OF MOLECULES density representation used for superposition can be varied in a straightforward manner. Results are shown for the superposition of two antiviral agents, 2rs1 and 2r04; the dihydrofolate reductase ligands, methotrexate and dihydrofolate; and a set of three «-thrombin inhibitors. Q 1997 by John Wiley & Sons, Inc.