SEARCH · PubMed Health

Results for “Molecular Docking Simulation”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,207 records · Page 67Linked to original sources

PubMed2003

Effects of irbesartan on cloned potassium channels involved in human cardiac repolarization.

We studied the effects of irbesartan, a selective angiotensin II type 1 receptor antagonist, on human ether-a-go-go-related gene (HERG), KvLQT1+minK, hKv1.5, and Kv4.3 channels using the patch-clamp technique. Irbesartan exhibited a low affinity for HERG and KvLQT1+minK channels (IC(50) = 193.0 +/- 49.8 and 314.6 +/- 85.4 microM, respectively). In hKv1.5 channels, irbesartan produced two types of block, depending on the concentration tested. At 0.1 microM, irbesartan inhibited the current in a time-dependent manner (22 +/- 3.9% at +60 mV). The blockade increased steeply with channel activation increasing at more positive potentials. However, at 10 microM, irbesartan induced a time-independent blockade that occurred in the range of potentials of channel opening, reaching its maximum at approximately 0 mV, and remaining unchanged at more positive potentials (24.0 +/- 1.0% at +60 mV). In Kv4.3 currents, irbesartan produced a concentration-dependent block, which resulted in two IC(50) values (1.0 +/- 0.1 nM and 7.2 +/- 0.6 microM). At 1 microM, it inhibited the peak current and accelerated the time course of inactivation, decreasing the total charge crossing the membrane (36.6 +/- 7.8% at +50 mV). Irbesartan shifted the inactivation curve of Kv4.3 channels, the blockade increasing as the amount of inactivated channels increased. Molecular modeling was used to define energy-minimized dockings of irbesartan to hKv1.5 and HERG channels. In conclusion, irbesartan blocks Kv4.3 and hKv1.5 channels at therapeutic concentrations, whereas the blockade of HERG and KvLQT1+minK channels occurred only at supratherapeutic levels. In hKv1.5, a receptor site is apparent on each alpha-subunit of the channel, whereas in HERG channels a common binding site is present at the pore.

Animals↗

PubMed1999

Antiidiotypic antibody recognizes an amiloride binding domain within the alpha subunit of the epithelial Na+ channel.

We previously raised an antibody (RA6.3) by an antiidiotypic approach which was designed to be directed against an amiloride binding domain on the epithelial Na+ channel (ENaC). This antibody mimicked amiloride in that it inhibited transepithelial Na+ transport across A6 cell monolayers. RA6.3 recognized a 72-kDa polypeptide in A6 epithelia treated with tunicamycin, consistent with the size of nonglycosylated Xenopus laevis alphaENaC. RA6.3 specifically recognized an amiloride binding domain within the alpha-subunit of mouse and bovine ENaC. The deduced amino acid sequence of RA6.3 was used to generate a three-dimensional model structure of the antibody. The combining site of RA6.3 was epitope mapped using a novel computer-based strategy. Organic residues that potentially interact with the RA6.3 combining site were identified by data base screening using the program LUDI. Selected residues docked to the antibody in a manner corresponding to the ordered linear array of amino acid residues within an amiloride binding domain on the alpha-subunit of ENaC. A synthetic peptide spanning this domain inhibited the binding of RA6.3 to alphaENaC. This analysis provided a novel approach to develop models of antibody-antigen interaction as well as a molecular perspective of RA6.3 binding to an amiloride binding domain within alphaENaC.

Amiloride↗

PubMed2006

Analysis of interactions responsible for vasopressin binding to human neurohypophyseal hormone receptors-molecular dynamics study of the activated receptor-vasopressin-G(alpha) systems.

Vasopressin (CYFQNCPRG-NH(2), AVP) is a semicyclic endogenous peptide, which exerts a variety of biological effects in mammals. The main physiological roles of AVP are the regulation of water balance and the control of blood pressure and adrenocorticotropin hormone (ACTH) secretion, mediated via three different subtypes of vasopressin receptors: V1a, V1b and V2 receptors (V1aR, V1bR and V2R, respectively). They are the members of the class A, G-protein-coupled receptors (GPCRs). AVP also modulates several behavioral and social functions. In this study, the interactions responsible for AVP binding to vasopressin V1a and V2 receptors versus the closely related oxytocin ([I3,L8]AVP, OT) receptor (OTR) have been investigated. Three-dimensional models of the activated receptors were constructed using multiple sequence alignment, followed by homology modeling using the complex of activated rhodopsin with Gt(alpha) C-terminal peptide of transducin MII-Gt(338-350) prototype as a template. AVP was docked into the receptor-G(alpha) systems. The three lowest-energy pairs of receptor-AVP-G(alpha) (two complexes per each receptor) were selected. The 1-ns unconstrained molecular dynamics (MD) of complexes embedded into the fully hydrated 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC) lipid bilayer was conducted in the AMBER 7.0 force field. Six relaxed receptor-AVP-G(alpha) models were obtained. The residues responsible for AVP binding to vasopressin receptors have been identified and a different mechanism of AVP binding to V2R than to V1aR has been proposed.

Amino Acid Sequence↗

PubMed2005

Intersurf: dynamic interface between proteins.

Protein docking is a fundamental biological process that links two proteins. This link is typically defined by an interaction between two large zones of the protein boundaries. Visualizing such an interface is useful to understand the process thanks to 3D protein structures, to estimate the quality of docking simulation results, and to classify interactions in order to predict docking affinity between classes of interacting zones. Since the interface may be defined by a surface that separates the two proteins, it is possible to create a map of interaction that allows comparisons to be performed in 2D. This paper presents a very fast algorithm that extracts an interface surface and creates a valid and low-distorted interaction map. Another benefit of our approach is that a pre-computed part of the algorithm enables the surface to be updated in real-time while residues are moved.

Algorithms↗

PubMed2005

CAPRI rounds 3-5 reveal promising successes and future challenges for RosettaDock.

CAPRI Rounds 3, 4, and 5 are the first public test of the published RosettaDock algorithm. The targets cover a wide range of sizes and shapes. For most targets, published biological information indicated the region of the binding site on at least one docking partner. The RosettaDock algorithm produced high accuracy predictions for three targets, medium-accuracy predictions for two targets, and an acceptable prediction for one target. RosettaDock predicted all five targets with less than 450 residues to high or medium accuracy, but it predicted only one of seven targets with above 450 residues to acceptable accuracy. RosettaDock's high-accuracy predictions for small to moderately large targets reveal the predictive power and fidelity of the algorithm, especially the high-resolution refinement and scoring protocol. In addition, RosettaDock can predict complexes from at least one homology-modeled docking partner with comparable accuracy to unbound cases of similar size. Larger targets present a more intensive sampling problem, and some large targets present repulsive barriers to entering the binding site. Ongoing improvements to RosettaDock's low-resolution search may alleviate this problem. This first public test suggests that RosettaDock can be useful in a significant range of applications in biochemistry and cell biology.

Algorithms↗

PubMed1993

Conservation of genetic information: a code for site-specific DNA recognition.

We present findings of genetic information conservation between the glucocorticoid response element (GRE) DNA and the cDNA encoding the glucocorticoid receptor (GR) DNA-binding domain (DBD). The regions of nucleotide sub-sequence similarity to the GRE in the GR DBD occur specifically at nucleotide sequences on the ends of exons 3,4, and 5 at their splice junction sites. These sequences encode the DNA recognition helix on exon 3, a beta-strand on exon 4, and a putative alpha-helix on exon 5, respectively. The nucleotide sequence of exon 5 that encodes the putative alpha-helix located on the carboxyl terminus of the GR DBD shares sequence similarity with the flanking nucleotide regions of the GRE. We generated a computer model of the GR DBD using atomic coordinates derived from nuclear magnetic resonance spectroscopy to which we attached the exon 5-encoded putative alpha-helix. We docked this GR DBD structure at the 39-base-pair nucleotide sequence containing the GRE binding site and flanking nucleotides, which contained conserved genetic information. We observed that amino acids of the DNA recognition helix, the beta-strand, and the putative alpha-helix are spatially aligned with trinucleotides identical to their cognate codons within the GRE and its flanking nucleotides.

Amino Acid Sequence↗

PubMed2006

Amber force field implementation, molecular modelling study, synthesis and MMP-1/MMP-2 inhibition profile of (R)- and (S)-N-hydroxy-2-(N-isopropoxybiphenyl-4-ylsulfonamido)-3-methylbutanamides.

Ab initio calculations (B3LYP/Lanl2DZ level of theory) were performed in this study to determine all the structural and catalytic zinc parameters required in order to study MMPs and their complexes with hydroxamate inhibitors by means of the AMBER force field. The parameters thus obtained were used in order to study the docking of some known MMPi (Batimastat, CGS 27023A and Prinomastat) and our previously described inhibitor a which had shown an inhibitory activity for MMP-1, and -2, with the aim of explaining the different selectivity. On this basis the two enantiomers (R)-b and (S)-b were designed and synthesized, as more potent MMP-2 inhibitors than our previously described inhibitor a. Between these two enantiomers the eutomer (R)-b proved to be 24.7 times and 15.3 times more potent than CGS 27023A and the parent compound a on MMP-2, maintaining a higher index of MMP-2/MMP-1 selectivity compared with CGS 27023A and the more potent inhibitor Prinomastat. The hydroxamate (R)-b can be considered as a progenitor of a new class of biphenylsulfonamido-based inhibitors that differ from compound a in the presence of an alkyl side chain on the C alpha atom, and show different potency and selectivity profiles on the two MMPs considered.

Catalysis↗

PubMed2005

Cooperativity and pseudo-cooperativity in the glutathione S-transferase from Plasmodium falciparum.

Binding and catalytic properties of glutathione S-transferase from Plasmodium falciparum (PfGST) have been studied by means of fluorescence, steady state and pre-steady state kinetic experiments, and docking simulations. This enzyme displays a peculiar reversible low-high affinity transition, never observed in other GSTs, which involves the G-site and shifts the apparent K(D) for glutathione (GSH) from 200 to 0.18 mM. The transition toward the high affinity conformation is triggered by the simultaneous binding of two GSH molecules to the dimeric enzyme, and it is manifested as an uncorrected homotropic behavior, termed "pseudo-cooperativity." The high affinity enzyme is able to activate GSH, lowering its pK(a) value from 9.0 to 7.0, a behavior similar to that found in all known GSTs. Using 7-chloro-4-nitrobenzo-2-oxa-1,3-diazole, this enzyme reveals a potential optimized mechanism for the GSH conjugation but a low catalytic efficiency mainly due to a very low affinity for this co-substrate. Conversely, PfGST efficiently binds one molecule of hemin/monomer. The binding is highly cooperative (n(H) = 1.8) and occurs only when GSH is bound to the enzyme. The thiolate of GSH plays a crucial role in the intersubunit communication because no cooperativity is observed when S-methylglutathione replaces GSH. Docking simulations suggest that hemin binds to a pocket leaning into both the G-site and the H-site. The iron is coordinated by the amidic nitrogen of Asn-115, and the two carboxylate groups are in electrostatic interaction with the epsilon-amino group of Lys-15. Kinetic and structural data suggest that PfGST evolved by optimizing its binding property with the parasitotoxic hemin rather than its catalytic efficiency toward toxic electrophilic compounds.

4-Chloro-7-nitrobenzofurazan↗

PubMed2005

Withanolides, a new class of natural cholinesterase inhibitors with calcium antagonistic properties.

The withanolides 1-3 and 4-5 isolated from Ajuga bracteosa and Withania somnifera, respectively, inhibited acetylcholinesterase (AChE, EC 3.1.1.7) and butyrylcholinesterase (BChE, EC 3.1.1.8) enzymes in a concentration-dependent fashion with IC50 values ranging between 20.5 and 49,2 microm and 29.0 and 85.2 microm for AChE and BChE, respectively. Lineweaver-Burk as well as Dixon plots and their secondary replots indicated that compounds 1, 3, and 5 are the linear mixed-type inhibitors of AChE, while 2 and 4 are non-competitive inhibitors of AChE with K(i) values ranging between 20.0 and 45.0 microm. All compounds were found to be non-competitive inhibitors of BChE with K(i) values ranging between 27.7 and 90.6 microm. Molecular docking study revealed that all the ligands are completely buried inside the aromatic gorge of AChE, while compounds 1, 3, and 5 extend up to the catalytic triad. A comparison of the docking results showed that all ligands generally adopt the same binding mode and lie parallel to the surface of the gorge. The superposition of the docked structures demonstrated that the non-flexible skeleton of the ligands always penetrates the aromatic gorge through the six-membered ring A, allowing their simultaneous interaction with more than one subsite of the active center. The affinity of ligands with AChE was found to be the cumulative effects of number of hydrophobic contacts and hydrogen bonding. Furthermore, all compounds also displayed dose-dependent (0.005-1.0 mg/mL) spasmolytic and Ca2+ antagonistic potentials in isolated rabbit jejunum preparations, compound 4 being the most active with an ED50 value of 0.09 +/- 0.001 mg/mL and 0.22 +/- 0.01 microg/mL on spontaneous and K+ -induced contractions, respectively. The cholinesterase inhibitory potential along with calcium antagonistic ability and safe profile in human neutrophil viability assay could make compounds 1-5 possible drug candidates for further study to treat Alzheimer's disease and associated problems.

Acetylcholinesterase↗

PubMed2002

Can we separate active from inactive conformations?

Molecular modeling methodologies such as molecular docking, pharmacophore modeling, and 3D-QSAR, rely on conformational searches of small molecules as a starting point. All of these methodologies seek conformations of the small molecules as they bind to target proteins, i.e., their active conformations. Thus the question as to whether active conformations can be separated from inactive conformations is extremely relevant. In this paper, 3D-descriptors that separate random conformations from active conformations of small molecules are sought. To select appropriate descriptors, 65 protein-ligand complexes were taken from the protein data bank. For each ligand the active conformation was compared to randomly generated low energy conformations. Descriptors such as solvent accessible surface area, number of internal interactions and radius of gyration appear to be useful for separating the active conformations from the random conformations. The results with all these descriptors indicate that active conformations are less compact that random conformations, i.e., they have more solvent accessible surface area, fewer internal interactions and a larger radius of gyration than random conformations. Thus these descriptors could be useful as weights to bias conformational search procedures to conformations more likely to bind to proteins or as filters to eliminate conformations unlikely to bind to any protein.

Computer Simulation↗

PubMed2003

Ribosomal localization of translation initiation factor IF2.

Bacterial translation initiation factor IF2 is a GTP-binding protein that catalyzes binding of initiator fMet-tRNA in the ribosomal P site. The topographical localization of IF2 on the ribosomal subunits, a prerequisite for understanding the mechanism of initiation complex formation, has remained elusive. Here, we present a model for the positioning of IF2 in the 70S initiation complex as determined by cleavage of rRNA by the chemical nucleases Cu(II):1,10-orthophenanthroline and Fe(II):EDTA tethered to cysteine residues introduced into IF2. Two specific amino acids in the GII domain of IF2 are in proximity to helices H3, H4, H17, and H18 of 16S rRNA. Furthermore, the junction of the C-1 and C-2 domains is in proximity to H89 and the thiostrepton region of 23S rRNA. The docking is further constrained by the requisite proximity of the C-2 domain with P-site-bound tRNA and by the conserved GI domain of the IF2 with the large subunit's factor-binding center. Comparison of our present findings with previous data further suggests that the IF2 orientation on the 30S subunit changes during the transition from the 30S to 70S initiation complex.

Base Sequence↗

PubMed2004

A model of voltage gating developed using the KvAP channel crystal structure.

Having inspected the crystal structure of the complete KvAP channel protein, we suspect that the voltage-sensing domain is too distorted to provide reliable information about its native tertiary structure or its interactions with the central pore-forming domain. On the other hand, a second crystal structure of the isolated voltage-sensing domain may well correspond to a native open conformation. We also observe that the paddle model of gating developed from these two structures is inconsistent with many experimental results, and suspect it to be energetically unrealistic. Here we show that the isolated voltage-sensing domain crystal structure can be docked onto the pore domain portion of the full-length KvAP crystal structure in an energetically favorable way to create a model of the open conformation. Using this as a starting point, we have developed rather conventional models of resting and transition conformations based on the helical screw mechanism for the transition from the open to the resting conformation. Our models are consistent with both theoretical considerations and experimental results.

Amino Acid Sequence↗

PubMed2004

Human carbonyl reductase catalyzes reduction of 4-oxonon-2-enal.

4-Oxonon-2-enal (4ONE) was demonstrated to be a product of lipid peroxidation, and previous studies found that it was highly reactive toward DNA and protein. The present study sought to determine whether carbonyl reductase (CR) catalyzes reduction of 4ONE, representing a potential pathway for metabolism of the lipid peroxidation product. Recombinant CR was cloned from a human liver cDNA library, expressed in Escherichia coli, and purified by metal chelate chromatography. Both 4ONE and its glutathione conjugate were found to be substrates for CR, and kinetic parameters were calculated. TLC analysis of reaction products revealed the presence of three compounds, two of which were identified as 4-hydroxynon-2-enal (4HNE) and 1-hydroxynon-2-en-4-one (1HNO). GC/MS analysis confirmed 4HNE and 1HNO and identified the unknown reaction product as 4-oxononanal (4ONA). Analysis of oxime derivatives of the reaction products via LC/MS confirmed the unknown as 4ONA. The time course for CR-mediated, NADPH-dependent 4ONE reduction and appearance of 4HNE and 1HNO was determined using HPLC, demonstrating 4HNE to be a major product and 1HNO and 4ONA to be minor products. Simulated structures of 4ONE in the active site of CR/NADPH calculated via docking experiments predict the ketone positioned as primary hydride acceptor. Results of the present study demonstrate that 4ONE is a substrate for CR/NADPH and the enzyme may represent a pathway for biotransformation of the lipid. Furthermore, these findings reveal that CR catalyzes hydride transfer selectively to the ketone but also to the aldehyde and C=C of 4ONE, resulting in 4HNE, 1HNO, and 4ONA, respectively.

Alcohol Oxidoreductases↗

PubMedDate not supplied

Structure-based QSAR analysis of a set of 4-hydroxy-5,6-dihydropyrones as inhibitors of HIV-1 protease: an application of the receptor-dependent (RD) 4D-QSAR formalism.

Receptor-dependent (RD) 4D-QSAR models were constructed for a set of 39 4-hydroxy-5,6-dihydropyrone analogue HIV-1 protease inhibitors. The receptor model used in this QSAR analysis was derived from the HIV-1 protease (PDB ID ) crystal structure. The bound ligand in the active site of the enzyme, also a 4-hydroxy-5,6-dihydropyrone analogue, was used as the reference ligand for docking the data set compounds. The optimized RD 4D-QSAR models are not only statistically significant (r(2) = 0.86, q(2) = 0.80 for four- and greater-term models) but also possess reasonable predictivity based on test set predictions. The proposed "active" conformations of the docked analogues in the active site of the enzyme are consistent in overall molecular shape with those suggested from crystallographic studies. Moreover, the RD 4D-QSAR models also "capture" the existence of specific induced-fit interactions between the enzyme active site and each specific inhibitor. Hydrophobic interactions, steric shape requirements, and hydrogen bonding of the 4-hydroxy-5,6-dihydropyrone analogues with the HIV-1 protease binding site model dominate the RD 4D-QSAR models in a manner again consistent with experimental conclusions. Some possible hypotheses for the development of new lead HIV-1 protease inhibitors can be inferred from the RD 4D-QSAR models.

Computer Simulation↗

PubMed2006

Interaction between the HIV-1 protein Vpr and the adenine nucleotide translocator.

The HIV-1 protein Vpr circulates in the serum of seropositive individuals and in the cerebrospinal fluid of AIDS patients with neurological disorders. Vpr triggers apoptosis of numerous cell types after extracellular addition, vpr gene transfer or in the context of viral infection. Moreover, in vivo, transgenic mice over-expressing Vpr have enhanced T lymphocytes apoptosis. In previous studies, we suggested that the Vpr apoptotic activities were because of its binding to the adenine nucleotide translocator (ANT), a mitochondrial ATP/ADP antiporter. To specify this interaction, fragments of both proteins were synthesized and used in biochemical and biophysical experiments. We demonstrate here that in vitro, the (27-51) and (71-82) Vpr peptides bind to a region encompassing the first ANT intermembrane space loop and part of its second and third transmembrane helices. Computational analysis using a docking program associated to dynamic simulations enabled us to construct a three-dimensional model of the Vpr-ANT complex. In this model, the N-terminus of Vpr plunges in the ANT cavity whereas the Vpr C-terminal extremity is located at the surface of the ANT allowing possible interactions with a third partner. These results could be used to design molecules acting as pro-apoptotic Vpr analogs or as apoptosis inhibitors preventing the Vpr-ANT interaction.

Animals↗

PubMed2005

Molecular design and characterization of recombinant long half-life mutants of human tissue factor pathway inhibitor.

Tissue factor pathway inhibitor (TFPI) is a physiological inhibitor of extrinsic pathway of coagulation and has biological functions of anticoagulation and anti-inflammation. Although TFPI has been proved to be a good therapeutic agent of sepsis, inflammatory shock, and DIC, the clinical application and therapeutic effects of TFPI are impeded because of its short half-life in vivo. In order to prolong the half-life of TFPI, homology modeling and molecule docking were performed on a computer workstation principally in protein structural biology and binding characteristics between TFPI and its receptor LRP (low-density lipoprotein receptor related protein). Two recombinant long half-life human TFPI mutants coined TFPI-Mut1 and TFPI-Mut4 were designed and expressed in E.coli. In comparison with the wild-type TFPI, TFPI-Mut1 and TFPI-Mut4 presented a few of changes in spatial configuration and a decrease in relative Gibbs free energy of docking complex by 17.3% and 21.5%, respectively, as indicated by a computer simulation. After refolding and purification, anticoagulant activities, anti-TF/FVIIa and anti-FXa activities of the mutants were found to be the same as those of wide-type TFPI. The pharmacokinetics research indicated that alpha phase half-life (t1/2 alpha) of TFPI-Mut1 and TFPI-Mut4 were prolonged 1.33-fold and 1.96-fold respectively, beta phase half-life (t1/2 beta) of TFPI-Mut1 and TFPI-Mut4 were prolonged 1.62-fold and 4.22-fold respectively. These results suggested that TFPI-Mut1 and TFPI-Mut4 maintained the bioactivities of wild-type TFPI, prolonged half-life in vivo simultaneously and were expected for better clinical value and therapeutic effect.

Animals↗

PubMed1996

Small structural changes account for the high thermostability of 1[4Fe-4S] ferredoxin from the hyperthermophilic bacterium Thermotoga maritima.

BACKGROUND: The characterization of the structural features that account for the high thermostability of some proteins is of great scientific and biotechnological interest. Proteins from hyperthermophilic organisms with optimum growth temperatures of 80 degrees C and higher generally show high intrinsic stabilities. The comparison of high resolution X-ray structures of these proteins with their counterparts from mesophilic organisms has therefore helped to identify potentially stabilizing forces in a number of cases. Small monomeric proteins which comprise only a single domain, such as ferredoxins, are especially suitable for such comparisons since the search for determinants of protein stability is considerably simplified. RESULTS: The 1.75 A crystal structure of the extremely thermostable 1[4Fe-4S] ferredoxin from Thermotoga maritima (FdTm) was determined and compared with other monocluster-containing ferredoxins with different degrees of thermostability. CONCLUSIONS: A comparison of the three-dimensional structure of FdTm with that of ferredoxins from mesophilic organisms suggests that the very high thermostability of FdTm is unexpectedly achieved without large changes of the overall protein structure. Instead, an increased number of potentially stabilizing features is observed in FdTm, compared with mesophilic ferredoxins. These include stabilization of alpha helices, replacement of residues in strained conformation by glycines, strong docking of the N-terminal methionine and an overall increase in the number of hydrogen bonds. Most of these features stabilize several secondary structure elements and improve the overall rigidity of the polypeptide backbone. The decreased flexibility will certainly play a relevant role in shielding the iron-sulfur cluster against physiologically high temperatures and further improve the functional integrity of FdTm.

Amino Acid Sequence↗

PubMed2005

Model structures of alpha-2 adrenoceptors in complex with automatically docked antagonist ligands raise the possibility of interactions dissimilar from agonist ligands.

Antagonist binding to alpha-2 adrenoceptors (alpha2-ARs) is not well understood. Structural models were constructed for the three human alpha2-AR subtypes based on the bovine rhodopsin X-ray structure. Twelve antagonist ligands (including covalently binding phenoxybenzamine) were automatically docked to the models. A hallmark of agonist binding is the electrostatic interaction between a positive charge on the agonist and the negatively charged side chain of D3.32. For antagonist binding, ion-pair formation would require deviations of the models from the rhodopsin structural template, e.g., a rotation of TM3 to relocate D3.32 more centrally within the binding cavity, and/or creation of new space near TM2/TM7 such that antagonists would be shifted away from TM5. Thus, except for the quinazolines, antagonist ligands automatically docked to the model structures did not form ion-pairs with D3.32. This binding mode represents a valid alternative, whereby the positive charge on the antagonists is stabilized by cation-pi interactions with aromatic residues (e.g., F6.51) and antagonists interact with D3.32 via carboxylate-aromatic interactions. This binding mode is in good agreement with maps derived from a molecular interaction library that predicts favorable atomic contacts; similar interaction environments are seen for unrelated proteins in complex with ligands sharing similarities with the alpha2-AR antagonists.

Adrenergic alpha-2 Receptor Agonists↗