Mode matches in hydrophobic free energy eigenfunctions predict peptide-protein interactions

The dominant statistical hydrophobic free energy inverse frequencies amino acid wavelengths as hydrophobic modes, of neurotensin (NT), cholescystokinin (CCK), the human dopamine D2 receptor [(DA)D2], and the human dopamine transporter (DAT) were determined using orthogonal decomposition of the autocovariance matrices of their amino acid sequences as hydrophobic free energy equivalents in kcal/mol. The leading eigenvalues-associated eigenvectors were convolved with the original series to construct eigenfunctions. Eigenfunctions were further analyzed using discrete trigonometric wavelet and all poles, maximum entropy power spectral transformations. This yielded clean representations of the dominant hydrophobic free energy modes, most of which are otherwise lost in the smoothing of hydropathy plots or contaminated by end effects and multimodality in conventional Fourier transformations. Mode matches were found between NT and (DA)D2 and between CCK and DAT, but not the converse. These mode matches successfully predicted the nonlinear kinetic interactions of NT-(DA)D2 in contrast with CCK-(DA) D2 on 3H-spiperone binding to (DA) D2, and by CCK-DAT but not NT-DAT on [N-methyl-3H]-WIN 35,428 binding to DAT in (DA)D2 and DAT cDNA stably transfected cell lines without known NT or CCK receptors. Computation of the dominant modes of hydrophobic free energy eigenfunctions may help predict functionally relevant peptide-membrane protein interactions, even across neurotransmitter families.