Synthesis and structure-activity relationships of a new model of arylpiperazines. Study of the 5-HT(1a)/alpha(1)-adrenergic receptor affinity by classical hansch analysis, artificial neural networks, and computational simulation of ligand recognition.

A classical quantitative structure-activity relationship (Hansch) study and artificial neural networks (ANNs) have been applied to a training set of 32 substituted phenylpiperazines with affinity for 5-HT(1A) and alpha(1)-adrenergic receptors, to evaluate the structural requirements that are responsible for 5-HT(1A)/alpha(1) selectivity. The resulting models provide a significant correlation of electronic, steric, and hydrophobic parameters with the biological affinities. Although the derived linear Hansch correlations give good statistics and acceptable predictions, the introduction of nonlinear relationships in the analysis gives more solid models and more accurate predictions. In the ANN models on the basis of the obtained 3D plots, the 5-HT(1A) affinity has a nonlinear dependence on F, V(o), V(m), and pi(o), although the nonlinear relationship is not far from a planar one. The alpha(1)-adrenergic receptor affinity has a clear nonlinear dependence on F, V(o), V(m), pi(o), and pi(m). A comparison of both analyses gives an additional understanding for 5-HT(1A)/alpha(1) selectivity: (a) high F values increase the binding affinity for 5-HT(1A) receptors and decrease the affinity for alpha(1) sites; (b) the hydrophobicity at the meta-position has only influence for the alpha(1)-adrenergic receptor; (c) the meta-position seems to be implicated in the 5-HT(1A)/alpha(1) selectivity. While the 5-HT(1A) receptor is able to accommodate bulky substituents in the region of its active site, the steric requirements of the alpha(1)-adrenergic receptor at this position are more restricted. This information was used for the design of the new ligand EF-7412 (33) (5-HT(1A): K(i exptl) = 27 nM, alpha(1): K(i exptl) > 1000 nM; 5-HT(1A): K(i pred) (ANN) = 36 nM, alpha(1): K(i pred ANN) = 2745 nM) which was characterized as an antagonist in vivo in pre- and postsynaptic 5-HT(1A)R sites. Computational simulations of the complex between EF-7412 (33) and a 3D model of the transmembrane domain of the 5-HT(1A) receptor allowed us to define the molecular details of the ligand-receptor interaction that includes: (i) the ionic interaction between the protonated amine of the ligand and Asp 3.32; (ii) the hydrogen bonds between the m-NHSO(2)Et group of the ligand and Asn 7.39; and the hydrogen bonds between the hydantoin moiety of the ligand and (iii) Thr 3.37, (iv) Ser 5.42, and (v) Thr 5.43. These QSAR and ANN results in combination with computational simulations of ligand recognition will be useful for the design of potent selective 5-HT(1A) ligands.

[1]  Tomoo Aoyama,et al.  Reconstruction of Weight Matrices in Neural Networks : A Method of Correlating Outputs with Inputs , 1991 .

[2]  W. G. Richards,et al.  Application of Neural Networks: Quantitative Structure-Activity Relationships of the Derivatives of 2,4-Diamino-5-(substituted-benzyl) pyrimidines as DHFR Inhibitors. , 1993 .

[3]  C. Dourish,et al.  Silent 5-HT1A receptor antagonists: utility as research tools and therapeutic agents. , 1993, Trends in pharmacological sciences.

[4]  J. Ballesteros,et al.  [19] Integrated methods for the construction of three-dimensional models and computational probing of structure-function relations in G protein-coupled receptors , 1995 .

[5]  K. Schaper,et al.  FREE-WILSON-TYPE ANALYSIS OF NON-ADDITIVE SUBSTITUENT EFFECTS ON THPB DOPAMINE RECEPTOR AFFINITY USING ARTIFICIAL NEURAL NETWORKS , 1999 .

[6]  N. L. Franklin,et al.  Statistical Analysis in Chemistry and the Chemical Industry. , 1955 .

[7]  K. Palczewski,et al.  Crystal Structure of Rhodopsin: A G‐Protein‐Coupled Receptor , 2000, Science.

[8]  F. Fanelli,et al.  Amino acids of the α 1B‐adrenergic receptor involved in agonist binding: differences in docking catecholamines to receptor subtypes , 1996, FEBS letters.

[9]  James Devillers,et al.  Neural Networks in QSAR and Drug Design , 1996 .

[10]  Kenneth A Jacobson,et al.  Molecular architecture of G protein‐coupled receptors , 1996, Drug development research.

[11]  B. Kobilka,et al.  Amino acid substitutions at position 312 in the seventh hydrophobic segment of the beta 2-adrenergic receptor modify ligand-binding specificity. , 1993, Molecular pharmacology.

[12]  H. Lester,et al.  The role of conserved aspartate and serine residues in ligand binding and in function of the 5‐HT1A receptor: A site‐directed mutation study , 1992, FEBS letters.

[13]  E. Meng,et al.  Rhodopsin Sees the Light , 2000, Science.

[14]  H. Kubinyi.  QSAR : Hansch analysis and related approaches , 1993 .

[15]  Y. Cheng,et al.  Relationship between the inhibition constant (K1) and the concentration of inhibitor which causes 50 per cent inhibition (I50) of an enzymatic reaction. , 1973, Biochemical pharmacology.

[16]  Toshio Fujita,et al.  A Method for the Correlation of Biological Activity and Chemical Structure , 2001 .

[17]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[18]  R. Glennon,et al.  The binding of propranolol at 5-hydroxytryptamine1D beta T355N mutant receptors may involve formation of two hydrogen bonds to asparagine. , 1996, Molecular pharmacology.

[19]  Geoffrey E. Hinton,et al.  Learning representations by back-propagating errors , 1986, Nature.

[20]  Johann Gasteiger,et al.  Neural Networks for Chemists: An Introduction , 1993 .

[21]  O. H. Lowry,et al.  Protein measurement with the Folin phenol reagent. , 1951, The Journal of biological chemistry.

[22]  E. Ambrosio,et al.  [3H]prazosin binding to central nervous system regions of male and female rats , 1984, Neuroscience Letters.

[23]  T. A. Andrea,et al.  Applications of neural networks in quantitative structure-activity relationships of dihydrofolate reductase inhibitors. , 1991, Journal of medicinal chemistry.

[24]  P. Kraulis.  A program to produce both detailed and schematic plots of protein structures , 1991 .

[25]  B. Kobilka,et al.  A point mutation in the seventh hydrophobic domain of the alpha 2 adrenergic receptor increases its affinity for a family of beta receptor antagonists. , 1991, The Journal of biological chemistry.

[26]  B. Kobilka,et al.  Identification of a single amino acid residue responsible for the binding of a class of beta-adrenergic receptor antagonists to 5-hydroxytryptamine1A receptors. , 1992, Molecular pharmacology.

[27]  Peter A. Kollman,et al.  Application of the multimolecule and multiconformational RESP methodology to biopolymers: Charge derivation for DNA, RNA, and proteins , 1995, J. Comput. Chem..

[28]  C. B. Davis,et al.  Activity of aromatic substituted phenylpiperazines lacking affinity for dopamine binding sites in a preclinical test of antipsychotic efficacy. , 1989, Journal of medicinal chemistry.