Topography of the active site of the noradrenaline neuronal membrane carrier based on the theoretical conformation analysis of inhibitors of neuronal catecholamine uptake

All equilibrium conformations for nine ligands, essentially different in structure, of the noradrenaline carrier through the synaptosomal membrane of rat hypothalamus were calculated by semi-empirical method. Among these compounds were amphetamine, methylphenidate, tricyclic antidepressants. The conformational energy minimization was performed in the space of torsional and bond angles. Geometrical characteristics of the conformers were determined in the cartesian coordinate system fixed relative to the benzene ring and nitrogen atom of the ligand. The selection of biologically active (productive) conformations was made according to the following criteria: 1) low conformational energy; 2) similarity of the nitrogen atoms and phenyl rings spatial disposition in all ligands; 3) accessibility for intermolecular interactions of the same sides of functional groupings in all ligands. The above criteria enabled the productive conformations for all ligands to be chosen unambiguously. The productive conformation of noradrenaline was found to have the Ph-C-C-N fragment in perpendicular trans-conformation. A topographic model for the carrier active site was suggested, its components being the nucleophilic and two arylophilic groups situated against the most accessible sides of the functional moieties of the productively bound ligands.     

Bioluminescence resonance energy transfer assays reveal ligand-specific conformational changes within preformed signaling complexes containing delta-opioid receptors and heterotrimeric G proteins

Heptahelical receptors communicate extracellular information to the cytosolic compartment by binding an extensive variety of ligands. They do so through conformational changes that propagate to intracellular signaling partners as the receptor switches from a resting to an active conformation. This active state has been classically considered unique and responsible for regulation of all signaling pathways controlled by a receptor. However, recent functional studies have challenged this notion and called for a paradigm where receptors would exist in more than one signaling conformation. This study used bioluminescence resonance energy transfer assays in combination with ligands of different functional profiles to provide in vivo physical evidence of conformational diversity of delta-opioid receptors (DORs). DORs and alpha(i1)beta(1)gamma(2) G protein subunits were tagged with Luc or green fluorescent protein to produce bioluminescence resonance energy transfer pairs that allowed monitoring DOR-G protein interactions from different vantage points. Results showed that DORs and heterotrimeric G proteins formed a constitutive complex that underwent structural reorganization upon ligand binding. Conformational rearrangements could not be explained by a two-state model, supporting the idea that DORs adopt ligand-specific conformations. In addition, conformational diversity encoded by the receptor was conveyed to the interaction among heterotrimeric subunits. The existence of multiple active receptor states has implications for the way we conceive specificity of signal transduction.     

Synthesis, alpha 1-adrenoceptor antagonist activity, and SAR study of novel arylpiperazine derivatives of phenytoin

In the search for new antiarrhythmic agents, some active 2-methoxyphenylpiperazine derivatives of phenytoin were obtained as a chemical modification of compound AZ-99 (3-ethyl-1-[2-hydroxy-3-(4-phenylpiperazin-1-yl)-propyl]-2,4-dioxo-5,5-diphenylimidazolidine). These compounds possessed structural properties similar to those of alpha(1)-adrenoceptor antagonists. In the present study, the affinities of the 2-methoxyphenylpiperazine derivatives (1a-3a) for alpha(1)- and alpha(2)-adrenoceptors were evaluated using radioligand ([(3)H]prazosin, [(3)H]clonidine) binding assays. In the next step, a new series of phenylpiperazine derivatives of phenytoin (4a-16a) containing 2-methoxyphenyl-, 2-ethoxyphenyl-, 2-pyridyl- or 2-furoylpiperazine moiety, as well as, various ester or alkyl substituents at 3-position of hydantoin ring were synthesized. The newly synthesized compounds were tested for their affinity to alpha(1)- and alpha(2)-adrenoceptors. They have shown affinities for alpha(1)-adrenoceptors at nanomolar to submicromolar range. Some compounds were moderately selective ligands of alpha(1)-adrenoceptors. Selected compounds (3a-5a, 7a, 13a, 14a) were also evaluated for their alpha(1)-adrenoceptor antagonistic properties in functional bioassays. A SAR study indicated that the most active compounds contain 2-alkoxyphenylpiperazine moieties and methyl or 2-methylpropionate substituent at 3-N position in hydantoin. The exchange of 2-alkoxyphenyl moiety into 2-furoyl or 2-pyridyl group significantly decreased affinities for alpha(1)-adrenoceptors. Molecular modelling results obtained using conformational analysis CONFLEX and PM5 method for geometry optimization, allowed for comparison of the spatial properties of tested compounds with pharmacophore model created by Barbaro et al. for the ideal alpha(1)-adrenoceptor antagonist.     

Comparative molecular field analysis of some pyridazinone-containing alpha1-antagonists

Diverse series of piperazines linked at N1 to 4, 5, or 6 positions of 3-(2H)-pyridazinone ring and at N4, by a suitable alkyl spacer, to some classical alpha1-adrenoceptor pharmacophore moieties, were tested in vitro for their alpha1-adrenoceptor antagonist activity. The modeling of their biological activity (pKb) by comparative molecular field analysis led to the development of a statistically significant partial least squares (PLS) model able to detect at 3-D level the main physicochemical interactions responsible for alpha1-adrenoceptor antagonist activity.     

Thermodynamics of agonist and antagonist binding to the alpha 1-adrenoceptor studied using [125I]BE 2254

The thermodynamics of binding of [125I]BE 2254 to the alpha 1-adrenoceptor in guinea pig brain membranes have been investigated at four different temperatures between 0 and 37 degrees C. The affinity and binding capacity of the radioligand did not vary with temperature. Thus, the change in enthalpy upon binding was close to zero whereas the change in entropy was large and positive (delta S degrees approximately 45 cal/mol-deg). In addition, [125I]BE 2254 has been used as a reporter ligand to probe the thermodynamics of the interaction of a variety of alpha-adrenoceptor agonists and antagonists with the alpha 1-adrenoceptor. Binding of all ligands was associated with large positive changes in entropy (delta S degrees between 18 and 48 cal/mol-deg) and little, or no, change in enthalpy, a finding that provides no convincing evidence for conformational rearrangement of alpha 1-adrenoceptors upon ligand binding.     

Ligand design and synthesis of new imidazo[5,1-b]quinazoline derivatives as alpha1-adrenoceptor agonists and antagonists

A series of new imidazo[5,1-b]quinazoline derivatives (VII-IX) was designed, synthesized, and biologically evaluated for their in vivo hypotensive or hypertensive activities. The design of these compounds was based upon the molecular modeling simulation of the fitting values and conformational energy values of the best-fitted conformers to both the alpha(1)-adrenoceptor (alpha(1)-AR) agonist and alpha(1)-adrenoceptor (alpha(1)-AR) antagonist hypotheses. These hypotheses were generated from their corresponding lead compounds using CATALYST software. The simulation studies predicted that compounds IXa and IXe would have probable affinity for the alpha(1)-AR antagonist hypothesis, while compounds IXb, IXc, and IXg predicted a higher affinity for the alpha(1)-AR agonist hypothesis. In vivo biological evaluation of these compounds for their effects on the blood pressure of normotensive cats was consistent with the results of molecular modeling studies, where compounds IXa and IXe exhibited hypotensive activity, while compounds IXb, IXc, and IXg resulted in increasing the blood pressure of the experimental animals at different doses.     

Theoretical proton affinities of alpha1 adrenoceptor ligands

A systematic study has been performed of the proton affinity of a large family of agonists and antagonists of the alpha1-adrenoceptor at the B3LYP/6-31G* level of theory. After a conformational search, all the N atoms were considered as protonation sites and protonation energy values were determined. The inclusion of solvation by means of the Onsager model yielded stabilization in the proton affinity values obtained. In addition, a good correlation was found between the proton affinity values corresponding to the first protonation in gas phase of some of the compounds and their corresponding experimental affinity constants K(i) for the alpha1A adrenergic receptor.     

Conformation-activity relationships of cyclic dermorphin analogues

A theoretical conformational analysis (molecular mechanics study) of nine cyclic tetrapeptides, structurally related to the highly mu-receptor-selective dermorphin analogue H-Tyr-D-Orn-Phe-Asp-NH2, was performed. These compounds display considerable diversity in their mu-receptor affinity and selectivity. A systematic search and subsequent energy minimization in absence of the exocyclic Tyr1 residue and Phe3 side chain revealed the constrained nature of the 11-13-membered ring structures contained in these analogues. No more than four low-energy conformers (within 2 kcal/mol of the lowest energy conformation) were found in each case. After attachment of the Tyr1 moiety and Phe3 side chain to the "bare" low-energy ring structures, a systematic search and energy minimization of these exocyclic moieties resulted in a limited number of low-energy conformational families for all compounds. Five analogues with high mu-receptor affinity--H-Tyr-D-Orn-Phe-Asp-NH2, H-Tyr-D-Orn-Phe-D-Asp-NH2, H-Tyr-D-Asp-Phe-Orn-NH2, H-Tyr-D-Asp-Phe-A2bu-NH2 (A2 bu: alpha, gamma-diaminobutyric acid) and H-Tyr-D-Cys-Phe-Cys-NH2--all showed a tilted stacking interaction between the Tyr1 and Phe3 aromatic rings in the lowest or second lowest energy conformation found. The same kind of stacking was not possible in low-energy conformers of the four analogues with poor affinity for the mu-receptor [H-Tyr-L-Orn-Phe-Asp-NH2, H-Tyr-D-Orn-D-Phe-Asp-NH2, H-Tyr-D-Orn-Phe(NMe)-Asp-NH2 [Phe(NMe): N alpha-methylphenylalanine], and H-Tyr-D-Orn-Phg-Asp-NH2 (Phg: phenylglycine)].(ABSTRACT TRUNCATED AT 250 WORDS)     

Analysis of efficacy of chiral adrenergic agonists

The origin of terms, affinity, intrinsic activity (or efficacy) and spare receptors has been reviewed. The Easson-Stedman theory (1933) in relation to the activation of adrenoceptors by agonists proved to be useful in the analysis of affinity and efficacy. Eudismic ratios of agonists provided critical information about the receptor-mediated activation. The evidence from circular dichroism spectroscopy with a fluorescent-tagged adrenoceptor agonist indicates a stereoselective interaction with the receptor. Thus, the simplest definition of efficacy may include the rate of change of the specific conformation of the receptor by the agonist, leading to the organized response. The functional groups of the potent enantiomer are postulated to interact in a "preferred" sequence with the receptor. The 7TM GPCR protein crystal structure of bovine rhodopsin was used as a model to construct the agonist interacting amino acid residues for alpha(1A)- and beta1-adrenoceptors. It was observed that both --+NH3 group and chiral --OH group of (-)-epinephrine interact with Asp106 TM III of alpha 1A-adrenoceptor. Similar interactions were observed for (+)-epinephrine but critical differences were observed. Enantiomers of epinephrine and oxymetazoline were also docked in the position at beta1-adrenoceptor to elucidate the conformational changes. Some unique information has emerged about the activation of adrenoceptors by agonists. The differences in the pharmacological efficacy of the enantiomers compare favorably with the dynamics of conformational changes by the agonist at alpha1A- and beta1 adrenoceptors.     

Label-free optical biosensor for ligand-directed functional selectivity acting on beta(2) adrenoceptor in living cells

Recent realization of ligand-directed functional selectivity demands high-resolution tools for studying receptor biology and ligand pharmacology. Here we use label-free optical biosensor to examine the dynamic mass redistribution of human epidermoid A431 cells in response to diverse beta(2)-adrenoceptor ligands. Multi-parameter analysis reveals distinct patterns in activation and signaling of the receptor induced by different agonists. Sequential and co-stimulation assays categorize various ligands for their ability to modulate signaling induced by catechol, a structural component of catecholamines. This study documents multiple ligand-specific states of the beta(2)-adrenoceptor and highlights the power of the biosensor assays for screening pathway-biased ligands.     

Examination of the conformational meaning of "delta-address" in the dermenkephalin sequence

Comprehensive energy calculations were applied to four opioid-related peptides with different receptor selectivities, namely the delta-selective dermenkephalin (Tyr-D-Met-Phe-His-Leu-Met-Asp-NH2, DRE), the mu-selective dermorphin (Tyr-D-Ala-Phe-Gly-Tyr-Pro-Ser-NH2, DRM) and their "hybrid" peptides DRM/DRE (Tyr-D-Ala-Phe-Gly-Leu-Met-Asp-NH2) and DRE/DRM (Tyr-D-Met-Phe-His-Tyr-Pro-Ser-NH2). It was shown that the N-terminal tripeptide "mu-messages" in the delta-selective ligands DRE and DRM/DRE can possess similar low energy space arrangements of their functionally important elements (the N-terminal alpha-amino group and the aromatic moieties of Tyr and Phe), but that these are different from the space arrangement of these moieties in mu-selective DRM and DRE/DRM. These results suggest that the C-terminal tripeptide "delta-address" in DRE may influence the conformation of the "mu-message" in DRM. A refined model for the delta-receptor-bound conformation of DRE is proposed based on these calculations which is similar to that previously suggested for the cyclic delta-selective peptide [D-Pen2, D-Pen5]enkephalin (DPDPE). This model also has partial correspondence with the structure of the delta-selective alkaloid naltrindole.     

Conformationally restricted macrocyclic analogues of combretastatins

New analogues of combretastatins have been evaluated as inhibitors of tubulin polymerization. These compounds present a macrocyclic structure, in which the para positions of the aromatic moieties have been linked by a 5- or 6-atoms chain, in order to produce a conformational restriction. This could contribute to determine the active conformation for these ligands. Such a conformational restriction and/or the steric hindrance makes them less potent inhibitors than the model compound CA-4.     

Alpha(1)-adrenoceptor activation: a comparison of 4-(anilinomethyl)imidazoles and 4-(phenoxymethyl)imidazoles to related 2-imidazolines

Literature reports suggest that disruption of an interhelical salt bridge is critical for alpha(1)-adrenoceptor activation, and the basic amine found in adrenergic receptor ligands is responsible for the disruption. Novel 4-(anilinomethyl)imidazoles and 4-(phenoxymethyl)imidazoles are agonists of the cloned human alpha(1)-adrenoceptors in vitro, and potent, selective alpha(1A)-adrenoceptor agonists have been identified in this series. These imidazoles demonstrate similar potencies and alpha(1)-subtype selectivities as the corresponding 2-substituted imidazolines. The extremely close SAR suggests that, in spite of the large difference in basicity, these imidazoles and imidazolines may establish the same interactions to activate alpha(1)-adrenoceptors.     

Synthesis and biological activity of new 1,4-benzodioxan-arylpiperazine derivatives. Further validation of a pharmacophore model for alpha(1)-adrenoceptor antagonists.

A series of WB4101 (1)-related benzodioxanes (2-17) have been synthesized by replacing the phenoxyethyl moiety of 1 with a N-alkyl piperazine bearing a cyclic substituent (a substituted or unsubstituted phenyl group, a pyridine or pyridazinone ring, a furoyl moiety) at the second nitrogen atom. The binding profile of these compounds has been assessed by radioligand receptor binding assay at alpha(1)- and alpha(2)-adrenoceptors, in comparison to prazosin and rauwolscine, respectively. Moreover, structure-activity relationships have been derived for compounds 2-17 based on their fitting to a pharmacophore model for alpha(1)-adrenoceptor antagonists recently proposed by our research group. In a parallel way, the same compounds have been used to further test the predictive power and statistical significance of the model itself. The accuracy of the results obtained also in this case revealed the robustness of the calculated pharmacophore model and led to the identification of the molecular structural moieties which are thought to contribute to the biological activity.     

Ligand-stabilized conformational states of human beta(2) adrenergic receptor: insight into G-protein-coupled receptor activation

G-protein-coupled receptors (GPCRs) are known to exist in dynamic equilibrium between inactive- and several active-state conformations, even in the absence of a ligand. Recent experimental studies on the β₂ adrenergic receptor (β₂AR) indicate that structurally different ligands with varying efficacies trigger distinct conformational changes and stabilize different receptor conformations. We have developed a computational method to study the ligand-induced rotational orientation changes in the transmembrane helices of GPCRs. This method involves a systematic spanning of the rotational orientation of the transmembrane helices (TMs) that are in the vicinity of the ligand for predicting the helical rotations that occur on ligand binding. The predicted ligand-stabilized receptor conformations are characterized by a simultaneous lowering of the ligand binding energy and a significant gain in interhelical and receptor-ligand hydrogen bonds. Using the β₂AR as a model, we show that the receptor conformational state depends on the structure and efficacy of the ligand for a given signaling pathway. We have studied the ligand-stabilized receptor conformations of five different ligands, a full agonist, norepinephrine; a partial agonist, salbutamol; a weak partial agonist, dopamine; a very weak agonist, catechol; and an inverse agonist, ICI-115881. The predicted ligand-stabilized receptor models correlate well with the experimentally observed conformational switches in β₂AR, namely, the breaking of the ionic lock between R131^3.50 at the intracellular end of TM3 (part of the DRY motif) and E268^6.30 on TM6, and the rotamer toggle switch on W286^6.48 on TM6. In agreement with trp-bimane quenching experiments, we found that norepinephrine and dopamine break the ionic lock and engage the rotamer toggle switch, whereas salbutamol, a noncatechol partial agonist only breaks the ionic lock, and the weak agonist catechol only engages the rotamer toggle switch. Norepinephrine and dopamine occupy the same binding region, between TM3, TM5, and TM6, whereas the binding site of salbutamol is shifted toward TM4. Catechol binds deeper into the protein cavity compared to the other ligands, making contact with TM5 and TM6. A part of the catechol binding site overlaps with those of dopamine and norepinephrine but not with that of salbutamol. Virtual ligand screening on 10,060 ligands on the norepinephrine-stabilized receptor conformation shows an enrichment of 38% compared to ligand unbound receptor conformation. These results show that ligand-induced conformational changes are important for developing functionally specific drugs that will stabilize a particular receptor conformation. These studies represent the first step toward a more universally applicable computational method for studying ligand efficacy and GPCR activation.     

Model of three-dimensional structure of vitamin D receptor and its binding mechanism with 1alpha,25-dihydroxyvitamin D(3).

Comparative modeling of the vitamin D receptor three-dimensional structure and computational docking of 1alpha,25-dihydroxyvitamin D(3) into the putative binding pocket of the two deletion mutant receptors: (207-423) and (120-422, Delta [164-207]) are reported and evaluated in the context of extensive mutagenic analysis and crystal structure of holo hVDR deletion protein published recently. The obtained molecular model agrees well with the experimentally determined structure. Six different conformers of 1alpha,25-dihydroxyvitamin D(3) were used to study flexible docking to the receptor. On the basis of values of conformational energy of various complexes and their consistency with functional activity, it appears that 1alpha,25-dihydroxyvitamin D(3) binds the receptor in its 6-s-trans form. The two lowest energy complexes obtained from docking the hormone into the deletion protein (207-423) differ in conformation of ring A and orientation of the ligand molecule in the VDR pocket. 1alpha,25-Dihydroxyvitamin D(3) possessing the A-ring conformation with axially oriented 1alpha-hydroxy group binds receptor with its 25-hydroxy substituent oriented toward the center of the receptor cavity, whereas ligand possessing equatorial conformation of 1alpha-hydroxy enters the pocket with A ring directed inward. The latter conformation and orientation of the ligand is consistent with the crystal structure of hVDR deletion mutant (118-425, Delta [165-215]). The lattice model of rVDR (120-422, Delta [164-207]) shows excellent agreement with the crystal structure of the hVDR mutant. The complex obtained from docking the hormone into the receptor has lower energy than complexes for which homology modeling was used. Thus, a simple model of vitamin D receptor with the first two helices deleted can be potentially useful for designing a general structure of ligand, whereas the advanced lattice model is suitable for examining binding sites in the pocket.     

Cell surface delivery and structural re-organization by pharmacological chaperones of an oligomerization-defective alpha(1b)-adrenoceptor mutant demonstrates membrane targeting of GPCR oligomers

Many G-protein-coupled receptors, including the alpha(1b)-adrenoceptor, form homo-dimers or oligomers. Mutation of hydrophobic residues in transmembrane domains I and IV alters the organization of the alpha(1b)-adrenoceptor oligomer, with transmembrane domain IV playing a critical role. These mutations also result in endoplasmic reticulum trapping of the receptor. Following stable expression of this alpha(1b)-adrenoceptor mutant, cell surface delivery, receptor function and structural organization were recovered by treatment with a range of alpha(1b)-adrenoceptor antagonists that acted at the level of the endoplasmic reticulum. This was accompanied by maturation of the mutant receptor to a terminally N-glycosylated form, and only this mature form was trafficked to the cell surface. Co-expression of the mutant receptor with an otherwise wild-type form of the alpha(1b)-adrenoceptor that is unable to bind ligands resulted in this wild-type variant also being retained in the endoplasmic reticulum. Ligand-induced cell surface delivery of the mutant alpha(1b)-adrenoceptor now allowed co-recovery to the plasma membrane of the ligand-binding-deficient mutant. These results demonstrate that interactions between alpha(1b)-adrenoceptor monomers occur at an early stage in protein synthesis, that ligands of the alpha(1b)-adrenoceptor can act as pharmacological chaperones to allow a structurally compromised form of the receptor to pass cellular quality control, that the mutated receptor is not inherently deficient in function and that an oligomeric assembly of ligand-binding-competent and -incompetent forms of the alpha(1b)-adrenoceptor can be trafficked to the cell surface by binding of a ligand to only one component of the receptor oligomer.     

Modification of human hemoglobin by glutathione. III. Perturbations of hemoglobin conformation analyzed by computer modeling

The perturbations of the conformation of human deoxyhemoglobin induced by the covalent attachment of glutathione at cysteine beta 93 have been investigated by computer simulation in conjunction with molecular graphics. In the first phase of the analysis, a systematic search was carried out of the conformational space of glutathione attached to deoxyhemoglobin. In this search, the conformation of the hemoglobin molecule was held constant, while the relative energies of a series of 186,624 glutathione conformations involving systematic variation of six dihedral angels were calculated. From this search, the most favorable conformation was selected as the starting conformation for energy minimization of the glutathionyl hemoglobin molecule as a function of all Cartesian coordinates. In order to provide a reference state, an independent minimization by the same procedures was carried out for deoxyhemoglobin in the absence of glutathione. Comparison of the minimized structures with and without glutathione attached revealed a number of significant differences. The most conspicuous difference in the protein moiety concerned the salt bridge between aspartate beta 94 and histidine beta 146 which is destabilized upon minimization of the glutathionyl-hemoglobin complex due to interactions of the aspartate residue with the glycyl NH group of glutathione. Other observed differences in the minimized structures are located at the alpha 1-beta 2 interface and include displacement of the carboxyl group of aspartate beta 99. In the minimized complex, the glutathione portion assumes a quasi-cyclic conformation stabilized through interactions between the free (gamma-glutamyl) amino and (glycyl) carboxyl ends of the tripeptide and between this carboxyl end and the epsilon amino group of lysine alpha 40. In a parallel conformational study of glutathione alone, a similar structure was found as the lowest energy form. These quasi-cyclic conformations contrast with the extended structures reported by Wright (Wright, W.B. (1955) Acta Crystallogr. 11, 632-642) for crystals of glutathione where interactions between molecules play a major role. The conclusions of our analysis are in agreement with the experimental investigations reported in the two preceding papers and permit, moreover, a coherent interpretation of the observed functional and structural changes in deoxyhemoglobin induced by glutathione.     

Trans-dominant inhibition of integrin function.

Occupancy of integrin adhesion receptors can alter the functions of other integrins and cause partition of the ligand-occupied integrin into focal adhesions. Ligand binding also changes the conformation of integrin extracellular domains. To explore the relationship between ligand-induced conformational change and integrin signaling, we examined the effect of ligands specific for integrin alpha IIb beta 3 on the functions of target integrins alpha 5 beta 1 and alpha 2 beta 1. We report that binding of integrin-specific ligands to a suppressive integrin can inhibit the function of other target integrins (trans-dominant inhibition). Trans-dominant inhibition is due to a blockade of integrin signaling. Furthermore, this inhibition involves both a conformational change in the extracellular domain and the presence of the beta cytoplasmic tail in the suppressive integrin. Similarly, ligand-induced recruitment of alpha IIb beta 3 to focal adhesions also involves a conformational rearrangement of its extracellular domain. These findings imply that the ligand-induced conformational changes can propagate from an integrin's extracellular to its intracellular face. Trans-dominant inhibition by integrin ligands may coordinate integrin signaling and can lead to unexpected biological effects of integrin-specific inhibitors.