Docking studies suggest ligand-specific δ-opioid receptor conformations

An automated docking procedure was used to study binding of a series of δ-selective ligands to three models of the δ-opioid receptor. These models are thought to represent the three ligand-specific receptor conformations. Docking results are in agreement with point mutation studies and suggest that different ligands—agonists and antagonists—may bind to the same binding site under different receptor conformations. Docking to different receptor models (conformations) also suggests that by changing to a receptor-specific conformation, the receptor may open or close different binding sites to other ligands. Figure  Ligands 5 (green) and 6 (orange) in bindingpocket BP1 of the R1 δ-opioid receptor model     

Interactions in pseudorotoxanes based on crown ether-secondary ammonium motifs. A theoretical study

A theoretical analysis of the nature of the interactions in dibenzo[24]crown-8 (DB24C8)-n-dibutylammonium (DBM)—pseudorotaxane complex at the MP2 and DFT levels shows that the main contribution to the binding energy is the electrostatic interaction with moderate (20–25%) correlation stabilization. The total binding energy in the DB24C8-DBM complex represents a sum of the binding energies of two NH–O and one CH–O hydrogen bonds and the latter constitutes about 25% of the total interaction energy, giving the total binding energy of −41.2 kcal mol⁻¹ at the BHandHLYP/6-311++G** level. Deprotonation of the DB24C8-DBM complex reduces the binding energy by some 50 kcal mol⁻¹, giving metastable complexes DB24C8-DBA-1 or DB24C8-DBA-2, which will dissociate to give free crown ether and n-dibutylamine because of the strong exchange repulsion that prevails in neutral complexes. Figure Formation of DB24C8-DBM pseudorotoxane complex     

ONIOM DFT/PM3 calculation on the interaction between STI-571 and abelson tyrosine kinase

The ONIOM2 (B3LYP/6–31G (d, p): PM3) and B3LYP/6–31G (d, p) methods were applied to investigate the interaction between STI-571 and abelson tyrosine kinase binding site. The complex of N-[4-methyl-3-(4-pyridin-3-yl-pyrimidin-2-ylamino)- phenyl]-benzamide (part of STI-571) and related 16 amino acid residues were found at B3LYP/6–31G (d, p) level to have hydrogen bonds and π....π stacking interaction, their binding energy via HAF optimization was −20.4 kcal mol⁻¹. The results derived from this study agreed well with the reported observation. Figure Optimized structure of STI-571 and Thr315 in abelson tyrosine kinase based on ONIOM2 method     

Pharmacophore mapping of a series of pyrrolopyrimidines,indolopyrimidines and their congeners as multidrug-resistance-associated protein (MRP1) modulators

Pharmacophore mapping studies were undertaken for a series of molecules belonging to pyrrolopyrimidines, indolopyrimidines and their congeners as multidrug resistance-associated protein (MRP1) modulators. A five-point pharmacophore with two hydrogen bond acceptors (A), one lipophilic/hydrophobic group (H), one positive ionic feature (P) and one aromatic ring (R) as pharmacophoric features was developed. The pharmacophore hypothesis yielded a statistically significant 3D-QSAR model, with a correlation coefficient of r ² = 0.799 for training set molecules. The model generated showed excellent predictive power, with a correlation coefficient Q ² = 0.679 for an external test set of 20 molecules. The pharmacophore was further validated using four structurally diverse compounds with MRP1 modulatory activity. These compounds mapped well onto four of the five features of the pharmacophore. The pharmacophore proposed here was then utilised for the successful retrieval of active molecules with diverse chemotypes from database search. The geometry and features of pharmacophore are expected to be useful for the design of selective MRP1 inhibitors. Figure Alignment of multidrug resistance-associated protein (MRP1) inhibitors with the developed pharmacophore. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The structure of 5a,6–anhydrotetracycline and its Mg2+ complexes in aqueous solution

Semiempirical molecular orbital theory has been used for a systematic scan of the binding positions for a Mg²⁺ ion with 5a,6–anhydrotetracycline taking both conformational flexibility and possible different tautomeric forms into account. The magnesium ion has been calculated alone and with four or five complexed water molecules in order to simulate the experimental situation more closely. The results are analyzed by comparing the behavior of the title compound with that of tetracycline itself and possible causes for the stronger induction of the Tetracycline Receptor (TetR) by 5a,6–anhydrotetracycline than by tetracycline are considered. Energetically favored 3D -structure of the zwitteranionic 5a,6-anhydrotetracycline magnesium complex in solution Electronic Supplementary Material Supplementary material is available for this article at     

Alpha particle chemistry. On the formation of stable complexes between He2+ and other simple species: implications for atmospheric and interstellar chemistry

The possibility that stable complexes may be formed between alpha particles (He²⁺) and small molecules is investigated using QCISD quantum mechanical calculations. Implications for their presence in the terrestrial atmosphere and/or in interstellar space are discussed. Figure Optimized structure of a stable H₂OHe²⁺ complex     

Computer-based design of novel HIV-1 entry inhibitors: neomycin conjugated to arginine peptides at two specific sites

Aminoglycoside–arginine conjugates (AAC and APAC) are multi-target inhibitors of human immunodeficiency virus type-1 (HIV-1). Here, we predict new conjugates of neomycin with two arginine peptide chains binding at specific sites on neomycin [poly-arginine-neomycin-poly-arginine (PA-Neo-PA)]. The rationale for the design of such compounds is to separate two short arginine peptides with neomycin, which may extend the binding region of the CXC chemokine receptor type 4 (CXCR4). We used homology models of CXCR4 and unliganded envelope glycoprotein 120 (HIV-1_IIIB gp120) and docked PA-Neo-PAs and APACs to these using a multistep docking procedure. The results indicate that PA-Neo-PAs spread over two negatively charged patches of CXCR4. PA-Neo-PA–CXCR4 complexes are energetically more favorable than AACs/APAC–CXCR4 complexes. Notably, our CXCR4 model and docking procedure can be applied to predict new compounds that are either inhibitors of gp120–CXCR4 binding without affecting stromal cell-derived factor 1α (SDF-1α) chemotaxis activity, or inhibitors of SDF-1α–CXCR4 binding resulting in an anti-metastasis effect. We also predict that PA-Neo-PAs and APACs can interfere with CD4–gp120 binding in unliganded conformation. Figure The r5-Neo-r5-CXCR4 complex. CXCR4 is shown in CPK representation. The negatively charged residues are shown in red and positively charged residues in blue. The r5-Neo-r5 is shown in stick representation, neomycin core is colored yellow and arginine moieties are colored magenta. Two negatively charged patches separated by neutral and positively charged residues are visible.     

Conformational stability and three-dimensional model of the δ-opioid pharmacophore for the extended antiparallel dimer structure of Met-enkephalin in water

The conformational stability of the extended antiparallel dimer structure of Met-enkephalin in water was analyzed by examining the hydration structure of enkephalin using molecular dynamics simulations. The result shows that, despite of the hydrophicility of the terminal atoms in the pentapeptide, the main contributor for the stability of the dimer in water is the four intermolecular hydrogen bonds between the Gly² and Phe⁴ groups. The three-dimensional model of the δ-opioid pharmacophore for this dimer structure was also established. Such a model was demonstrated to match the δ-opioid pharmacophore query derived from the non-peptides SIOM, TAN-67, and OMI perfectly. This result thus strongly supports the assumption that the dimer structure of Met-enkephalin is a possible δ-receptor binding conformation. Figure Schematic model of the extended antiparallel dimer structure of Met-enkephalin     

The platinum-olefin binding energy in series of (PH<Subscript>3</Subscript>)<Subscript>2</Subscript>Pt(olefin) complexes - a theoretical study

Theoretical investigation of Pt(0)-olefin organometallic complexes containing tertiary phosphine ligands was focused on the strength of platinum-olefin electronic interaction. DFT theoretical study of electronic effects in a substantial number of ethylene derivatives was evaluated in terms of the Pt-olefin binding energy using MP2 correlation theory. Organometallics bearing coordinated olefins with general formula (R¹R²C = CR³R⁴)Pt(PH₃)₂ [R = various substituents] had been selected, including olefins containing both electron-donor substituents as well as electron-withdrawing groups. The stability of the corresponding complexes increases with a strengthening electron-withdrawal ability of the olefin substituents. Figure Representation of (CH₂ = CHR)Pt(PPh₃)₂ and the stability chart     

Modeling study of the influences of the aromatic transitions and the local environment on the far-UV rotational strengths in TEM-1 β-lactamase

Rotational strengths in the far-UV of TEM-1 β-lactamase have been investigated with two theoretical models based on the matrix method. The first model excludes, and a second includes, effects of the local electrostatic interactions on the chromophore energies. Special attention is given to the contributions of the aromatic side-chain chromophores, and the mechanisms of generation of rotational strengths are analyzed. The sensitivity of the computational models with respect to the structural changes of the protein are discussed. Figure Structure of TEM-1 β-lactamase. Both domains—α and αβ, the secondary structural elements and the aromatic and disulfide chromophores are shown     

G2, G3, and complete basis set calculations of the thermodynamic properties of cis- and trans-triazene

Following our recent study on triazane, we present a follow-up study on the thermodynamic properties of triazane’s unsaturated analog, triazene. We predict optimized structural parameters, vibrational frequencies, enthalpies of formation, enthalpies of combustion, specific enthalpies of combustion, and proton affinities. Our results indicate that the cis form of triazene has a specific enthalpy of combustion of −15.2 kJ g⁻¹ and the trans form has a specific enthalpy of combustion of −14.7 kJ g⁻¹. Figure Structures of cis- and trans-triazane, N₃H₃     

Computational studies of the binding site of α<Subscript>1A</Subscript>-adrenoceptor antagonists

Aimed at achieving a good understanding of the 3-dimensional structures of human α_1A-adrenoceptor (α_1A-AR), we have successfully developed its homology model based on the crystal structure of β₂-AR. Subsequent structural refinements were performed to mimic the receptor’s natural membrane environment by using molecular mechanics (MM) and molecular dynamics (MD) simulations in the GBSW implicit membrane model. Through molecular docking and further simulations, possible binding modes of subtype-selective α_1A-AR antagonists, Silodosin, RWJ-69736 and (+)SNAP-7915, were examined. Results of the modeling and docking studies are qualitatively consistent with available experimental data from mutagenesis studies. The homology model built should be very useful for designing more potent subtype-selective α_1A-AR antagonists and for guiding further mutagenesis studies. Figure The superposition of β₂-AR crystal structure (gold ribbons) and α_1A-AR homology model (blue ribbons)     

Design and synthesis of thiourea based receptor containing naphthalene as oxalate selective sensor

The thiourea based receptor containing naphthalene groups (1), has been successfully designed and synthesized for application as an oxalate receptor. A density functional theory at B3LYP/6-31G(d,p) level of theory has been applied to predict the binding ability between 1 and selected anions, i.e., oxalate, malonate, succinate, glutarate, dihydrogen phosphate, and hydrogen sulphate. Calculation results point out that receptor 1 shows the strongest interaction to oxalate ion with the binding free energy of 172.48 kcal mol⁻¹. The recognition ability of 1 to the selected anions has been also investigated by means of the absorption and emission techniques. Experimental results are in excellent agreement with the calculation data in which receptor 1 shows highly selective for oxalate ion over the other anions with logβ of 3.82 (0.02) M⁻¹ by means of the size of binding cavity.      

Model of the whole rat AT1 receptor and the ligand-binding site

We present a three-dimensional model of the rat type 1 receptor (AT₁) for the hormone angiotensin II (Ang II). Ang II and the AT₁ receptor play a critical role in the cell-signaling process responsible for the actions of renin–angiotensin system in the regulation of blood pressure, water-electrolyte homeostasis and cell growth. Development of improved therapeutics would be significantly enhanced with the availability of a 3D-structure model for the AT₁ receptor and of the binding site for agonists and antagonists. This model was constructed using a combination of computation and homology-modeling techniques starting with the experimentally determined three-dimensional structure of bovine rhodopsin (PDB#1F88) as a template. All 359 residues and two disulfide bonds in the rat AT₁ receptor have been accounted for in this model. Ramachandran-map analysis and a 1 nanosecond molecular dynamics simulation of the solvated receptor with and without the bound ligand, Ang II, lend credence to the validity of the model. Docking calculations were performed with the agonist, Ang II and the antihypertensive antagonist, losartan.      

Structures of inclusion complexes of halogenbenzoic acids and α-cyclodextrin based on AM1 calculations

Semi-empirical quantum mechanics calculations using AM1 (Austin Method 1) were carried out for various host-guest combinations of α-cyclodextrin and mono-halogen benzoic acids. The energetically favorable inclusion structures were identified. The AM1 results show that α-cyclodextrin complexes with mono-halogen benzoic acid acids (where the halogen is chlorine, bromide, iodine) as guest compounds are more stable in the “head first” position than in the “tail-first” position for meta and para isomers while ortho mono-halogen benzoic acids complexes with α-cyclodextrin are more stable in “tail-first” position. The calculated structures were found to be in good agreement with those obtained from crystalographic databases.      

Homology modeling and examination of the effect of the D92E mutation on the H5N1 nonstructural protein NS1 effector domain

Virulent H5N1 strains of influenza virus often harbor a D92E point mutation in the nonstructural protein NS1. This crucial mutation has been correlated with increased virulence and/or cytokine resistance, but the structural implications of such a change are still unclear. Furthermore, NS1 protein could also be a potential target for the development of novel antiviral agents against H5N1 strains. Therefore, a reasonable 3D model of H5N1 NS1 is important for the understanding of the molecular basis of increased virulence and the design of novel antiviral agents. Based on the crystal structure of a non-H5N1 NS1 protein, a model of H5N1 NS1 was developed by homology modeling, molecular mechanics and molecular dynamics simulations. It was found that the D92E mutation could result in weakened interactions of the carboxylate side chain with other phosphorylated residues, thereby activating phosphorylation of NS1. Figure Superposition of snapshots picked from the two molecular dynamic (MD) trajectories: a H5N1 NS1 homology model and b non-H5N1 NS1 crystal structure after 0 (green ribbon), 5 (blue ribbon) and 10 ns (pink ribbon) MD simulation     

Acetylcholine nicotinic receptors: finding the putative binding site of allosteric modulators using the “blind docking” approach

Allosteric potentiation of acetylcholine nicotinic receptors is considered to be one of the most promising approaches for the treatment of Alzheimer’s disease. However, the exact localization of the allosteric binding site and the potentiation mechanism at the molecular level are presently unknown. We have performed the “blind docking” of three known allosteric modulators (galanthamine, codeine and eserine) with the Acetylcholine Binding Protein and models of human α7, α3β4 and α4β2 nicotinic receptors, created by homology modeling. Three putative binding sites were identified in the channel pore, each one showing different affinities for the ligands. One of these sites is localized opposite to the agonist binding site and is probably implicated in the potentiation process. On the basis of these results, a possible mechanism for nicotinic acetylcholine receptor (nAChRs) activation is proposed. The present findings may represent an important advance for understanding the allosteric modulation mechanism of nAChRs. Electronic supplementary material Supplementary material is available for this article at