In silico screening for agonists and blockers of the β(2) adrenergic receptor: implications of inactive and activated state structures

Ten crystal structures of the β(2) adrenergic receptor have been published, reflecting different signaling states. Here, through controlled-docking experiments, we examined the implications of using inactive or activated structures on the in silico screening for agonists and blockers of the receptor. Specifically, we targeted the crystal structures solved in complex with carazolol (2RH1), the neutral antagonist alprenalol, the irreversible agonist FAUC50 (3PDS), and the full agonist BI-167017 (3P0G). Our results indicate that activated structures favor agonists over blockers, whereas inactive structures favor blockers over agonists. This tendency is more marked for activated than for inactive structures. Additionally, agonists tend to receive more favorable docking scores when docked at activated rather than inactive structures, while blockers do the opposite. Hence, the difference between the docking scores attained with an activated and an inactive structure is an excellent means for the classification of ligands into agonists and blockers as we determined through receiver operating characteristic curves and linear discriminant analysis. With respect to virtual screening, all structures prioritized well agonists and blockers over nonbinders. However, inactive structures worked better for blockers and activated structures worked better for agonists, respectively. Notably, the combination of individual docking experiments through receptor ensemble docking resulted in an excellent performance in the retrieval of both agonists and blockers. Finally, we demonstrated that the induced-fit docking of agonists is a viable way of modifying an inactive crystal structure and bias it toward the in silico recognition of agonists rather than blockers.     

Structure-based virtual screening of small-molecule antagonists of platelet integrin αIIbβ3 that do not prime the receptor to bind ligand

Integrin αIIbβ3 has emerged as an important therapeutic target for thrombotic vascular diseases owing to its pivotal role in mediating platelet aggregation through interaction with adhesive ligands. In the search for effective anti-thrombotic agents that can be administered orally without inducing the high-affinity ligand binding state, we recently discovered via high-throughput screening of 33,264 compounds a novel, αIIbβ3-selective inhibitor (RUC-1) of adenosine-5′-diphosphate (ADP) -induced platelet aggregation that exhibits a different chemical scaffold and mode of binding with respect to classical Arg-Gly-Asp (RGD)-mimicking αIIbβ3 antagonists. Most importantly, RUC-1 and its higher-affinity derivative, RUC-2, do not induce major conformational changes in the protein β3 subunit or prime the receptor to bind ligand. To identify additional αIIbβ3-selective chemotypes that inhibit platelet aggregation through similar mechanisms, we screened in silico over 2.5 million commercially available, ‘lead-like’ small molecules based on complementarity to the predicted binding mode of RUC-2 into the RUC-1-αIIbβ3 crystal structure. This first reported structure-based virtual screening application to the αIIbβ3 integrin led to the identification of 2 αIIbβ3-selective antagonists out of 4 tested, which compares favorably with the 0.003?% “hit rate” of our previous high-throughput chemical screening study. The newly identified compounds, like RUC-1 and RUC-2, showed specificity for αIIbβ3 compared to αVβ3 and did not prime the receptor to bind ligand. They thus may hold promise as αIIbβ3 antagonist therapeutic scaffolds.     

Structural basis for the β-adrenergic receptor subtype selectivity of the representative agonists and antagonists

The β(3)-adrenegic receptor (β(3)-AR) selectivity over β(1)- and β(2)-ARs has been the most important aspect for successful therapeutic agents for obesity and type-II diabetes, as the concomitant activation of β(1)- and β(2)-ARs would lead to undesirable side effects, such as increased heart rate. In order to explore the structural basis for the β-AR subtype selectivity of agonists and anatagonists, a three-dimensional structure of until date unresolved β(3)-AR has been modeled, compared with the resolved X-ray structures of β(1)- and β(2)-ARs, and used to study its stereoselective binding with until-date known diverse classes of representative agonists and antagonist. The obtained binding structures and calculated prime molecular mechanics-generalized Born surface area (MM-GBSA) binding free energies consistently reveal that while the subtype selectivity is strongly governed by the residues present in the extracellular ends of TM3, TM5, TM6, TM7 helices and of the ECL2 domain, the binding affinity is governed by the conserved residues present in the deep pocket limiting the degree of conformational and rotational freedoms to the bound ligand. The study demonstrates that the key structural requirements for the β(3)-selectivity are: (i) a negatively ionizable group (NIG) for direct interaction with β(3)-specific residue R315(6.58), (ii) a linker (9-10 ? length) between the protonated amine and NIG, and (iii) a substituted aryl ring directly attached to the β-hydroxyl carbon. The new computational insights acquired in this study are expected to be valuable in structure-based rational design of high-affinity agonists and antagonists with pronounced β(3)-selectivity for successful therapeutic agents for type-II diabetes and obesity.     

Electrospray mass spectrometry for the direct accurate mass measurement of ligands in complex with the retinoid X receptor α ligand binding domain

Accurate mass measurements are often used in the structural determination of unknown compounds of low molecular mass (i.e., below approximately 500 Da). Recently, it has been shown that accurate mass measurements also can be made on small denatured proteins (i.e., M(r), approximately 17,000) to confirm their amino acid composition and identify the presence of isoforms. In the current report, we present nondenaturing electrospray (ES) mass spectrometry data on the direct accurate mass measurement of ligands in complex with the retinoid X receptor ligand binding domain (RXR LBD; M(r) 31,370.92). Average mass errors were below 0.198 Da, 6.3 ppm (standard deviation [SD], 0.146; n = 10) for low-affinity fatty acid agonists analyzed in complex with the RXR LBD. Protein consumption was less than 15 pmol, with fatty acid ligands present at concentrations corresponding to their median effective concentration value (low micromolar, determined in transfection assays). Although determination of fatty acid mass was only sufficiently accurate to give nominal mass values, measurements were of sufficient accuracy to assign fatty acid chain length, degree of unsaturation, or cyclization. Using 17beta-estradiol as a control, the ability to observe specific ligand binding is shown for both high- and low-affinity RXRalpha agonists. In addition, binding of a novel synthetic receptor agonist XCT0315908 to the RXRalpha LBD is reported. This compound showed a high degree of complex formation, and the receptor-ligand complex could be mass measured with an average mass error of -0.024 Da, 0.8 ppm (SD, 0.092; n = 9). Thus, specific binding of both nanomolar and micromolar affinity ligands to a nuclear receptor LBD can be directly observed using nondenaturing ES mass spectrometry and accurate mass measurements additionally can be made on intact complexes in the same experiment. This methodology also is applicable when ligands are present as components of mixtures.     

Intermolecular interactions and electrostatic properties of the β-hydroquinone apohost: implications for supramolecular chemistry

The crystal structure of the β-polymorph of hydroquinone (β-HQ), the apohost of a large family of clathrates, is reported with a specific focus on intermolecular interactions and the electrostatic nature of its cavity. Hirshfeld surface analysis reveals subtle close contacts between two interconnecting HQ networks, and the local packing and related close contacts were examined by breakdown of the fingerprint plot. An experimental multipole model containing anisotropic thermal parameters for hydrogen atoms has been successfully refined against 15(2) K single microcrystal synchrotron X-ray diffraction data. The experimental electron density model has been compared with a theoretical electron density calculated with the molecule embedded in its own crystal field. Hirshfeld charges, interaction energies and the electrostatic potential calculated for both models are qualitatively in good agreement, but small differences in the electrostatic potential persist due to charge transfer from all hydrogen atoms to the oxygen atoms in the theoretical model. The electrostatic potential in the center of the cavity is positive, very shallow and highly symmetric, suggesting that the inclusion of polar molecules in the void will involve a balance between opposing effects. The electric field is by symmetry zero in the center of the cavity, increasing to a value of 0.0185 e/?(2) (0.27 V/?) 1 ? along the 3-fold axis and 0.0105 e/?(2) (0.15 V/?) 1 ? along the perpendicular direction. While these values are substantial in a macroscopic context, they are quite small for a molecular cavity and are not expected to strongly polarize a guest molecule.     

Homology modelling of the human adenosine A2B receptor based on X-ray structures of bovine rhodopsin, the β2-adrenergic receptor and the human adenosine A2A receptor

A three-dimensional model of the human adenosine A_2B receptor was generated by means of homology modelling, using the crystal structures of bovine rhodopsin, the β₂-adrenergic receptor, and the human adenosine A_2A receptor as templates. In order to compare the three resulting models, the binding modes of the adenosine A_2B receptor antagonists theophylline, ZM241385, MRS1706, and PSB601 were investigated. The A_2A-based model was much better able to stabilize the ligands in the binding site than the other models reflecting the high degree of similarity between A_2A and A_2B receptors: while the A_2B receptor shares about 21% of the residues with rhodopsin, and 31% with the β₂-adrenergic receptor, it is 56% identical to the adenosine A_2A receptor. The A_2A-based model was used for further studies. The model included the transmembrane domains, the extracellular and the intracellular hydrophilic loops as well as the terminal domains. In order to validate the usefulness of this model, a docking analysis of several selective and nonselective agonists and antagonists was carried out including a study of binding affinities and selectivities of these ligands with respect to the adenosine A_2A and A_2B receptors. A common binding site is proposed for antagonists and agonists based on homology modelling combined with site-directed mutagenesis and a comparison between experimental and calculated affinity data. The new, validated A_2B receptor model may serve as a basis for developing more potent and selective drugs.     

A sequence-based computational model for the prediction of the solvent accessible surface area for α-helix and β-barrel transmembrane residues

Predicting the solvent accessible surface area (ASA) of transmembrane (TM) residues is of great importance for experimental researchers to elucidate diverse physiological processes. TM residues fall into two major structural classes (α-helix membrane protein and β-barrel membrane protein). The reported solvent ASA prediction models were developed for these two types of TM residues respectively. However, this prevents the general use of these methods because one cannot determine which model is suitable for a given TM residue without information of its type. To conquer this limitation, we developed a new computational model that can be used for predicting the ASA of both TM α-helix and β-barrel residues. The model was developed from 78 α-helix membrane protein chains and 24 β-barrel membrane protein. Its prediction ability was evaluated by cross validation method and its prediction result on an independent test set of 20 membrane protein chains. The results show that our model performs well for both types of TM residues and outperforms other prediction model which was developed for the specific type of TM residues. The prediction results also proved that the random forest model incorporating conservation score is an effective sequence-based computational approach for predicting the solvent ASA of TM residues.     

Liquid—liquid equilibria for ternary systems of water—phenol and solvents: new UNIFAC parameters for interactions between the groups ACOH/COOH and prediction

The UNIFAC group-contribution method is used to predict ternary liquid—liquid equilibrium data presented in a recent paper (Alvarez Gonzalez et al.) for the systems water/phenol/benzene, water/phenol/ethylbenzene, water/phenol/nonanoic acid, water/phenol/ethyl acetate, water/phenol/isopropyl acetate, water/phenol/n-butyl acetate, water/phenol/isoamyl acetate and water/phenol/cyclohexyl acetate at 25°C and water/phenol/n-hexyl acetate at 25, 35 and 45°C. New UNIFAC interaction parameters between the groups ACOH/COOH have been obtained.A comparison between the experimental and predicted values is presented.     

Exploring the conformational and biological versatility of β-turn-modified gramicidin S by using sugar amino acid homologues that vary in ring size

Monobenzylated sugar amino acids (SAAs) that differ in ether ring size (containing an oxetane, furanoid, and pyranoid ring) were synthesized and incorporated in one of the β-turn regions of the cyclo-decapeptide gramicidin S (GS). CD, NMR spectroscopy, modeling, and X-ray diffraction reveal that the ring size of the incorporated SAA moieties determines the spatial positioning of their cis-oriented carboxyl and aminomethyl substituents, thereby subtly influencing the amide linkages with the adjacent amino acids in the sequence. Unlike GS itself, the conformational behavior of the SAA-containing peptides is solvent dependent. The derivative containing the pyranoid SAA is slightly less hydrophobic and displays a diminished haemolytic activity, but has similar antimicrobial properties as GS.     

Bioactivity screening and mass spectrometric confirmation for the detection of PPARδ agonists that increase type 1 muscle fibres

Sensitive and robust bioassays able to detect nuclear receptor activation are very useful for veterinary and doping control, pharmaceutical industry and environmental scientists. Here, we used bioassays based on human leukemic monocyte lymphoma U937 and human liver hepatocellular carcinoma HepG2 cell lines to detect the ligand-induced activation of the peroxisome proliferator-activated receptor delta (PPARδ). Exposure of U937 cells to the PPARδ agonist GW501516 resulted in a marked increase in mRNA expression of the PPARδ target gene Angptl4 which was quantified by qRT-PCR analysis. Exposure of HepG2 cells transiently transfected with a PPARδ expression plasmid and a PPAR-response element-driven luciferase reporter plasmid to PPARδ agonists GW501516, GW610742 and L-165041 resulted in clear dose–response curves. Although the qRT-PCR resulted in higher fold inductions, the luciferase assay with transfected HepG2 cells is cheaper and quicker and about ten times more sensitive to GW501516 compared to analysis of Angptl4 mRNA expression in U937 cells by qRT-PCR. The HepG2-based luciferase assay was therefore used to screen GW501516-spiked supplements and feed and water samples. After liquid extraction and clean-up by solid phase extraction using a weak anion exchange column, extracts were screened in the HepG2 bioassay followed by confirmation with a newly developed UPLC-MS/MS method, using two transitions for each compound, i.e., for GW501516, 454.07>188.15 (collision energy (CE) 46 V) and 454.07>257.08 (CE 30 V); for GW610742, 472.07>206.2 (CE 48 V) and 472.07>275.08 (CE 30 V); and for L-165041, 401.2>193.15 (CE 26 V) and 401.2>343.2 (CE 20 V).     

Modeling and active site refinement for G protein-coupled receptors: application to the β-2 adrenergic receptor

It is well known that G protein-coupled receptors are prime targets for drug discovery. At the present time there is only one protein from this class that has an X-ray crystal structure, bovine rhodopsin. Crystal structures of rhodopsin have become invaluable templates for the modeling of class-A G protein-coupled receptors as they likely represent the overall topology of this family of proteins. However, because of low sequence homology within the class and the inherent mobility of integral membrane proteins, it is unlikely that this single structural template reflects the ensemble of conformations accessible for any given receptor. We have devised a procedure based upon comparative modeling that uses induced fit modeling coupled with binding site expansion. The modeling protocol enables an ensemble approach to binding mode prediction. The utility of models for β-2 adrenergic receptor will be discussed.     

New functionalization of graphene oxide with N2O2 ligand for efficient loading of Cu nanostructures as a heterogeneous nanocatalyst for the synthesis of β-hydroxy-1,2,3-triazoles

A new and heterogeneous copper complex immobilized on graphene oxide (GO) was prepared. This was achieved through organic functionalization of GO using 1,8-diamino-3,6-dioxaoctane (DADO) and then inorganic coordination of copper on the edges and basal plane of the functionalized GO (GO-DADO-Cu), which was reduced to Cu(0). The chemical structure of the prepared nanocatalyst was analyzed using various techniques. Most of the analyses confirmed the successful anchoring of copper and organic ligand on the GO surface. Moreover, the synthesized nanocatalyst has shown high catalytic activity in the synthesis of β-hydroxy-1,2,3-triazole derivatives under mild reaction conditions (water and room temperature) resulting in good to excellent yields.     

Integration of genetic virtual screening patterns and latent multivariate modeling techniques for QSAR optimization based on combinations and/or interactions between peptides and proteins

Both the concept and the model of snug quantitative structure-activity relationship (QSAR) were pro-posed and developed for molecular design through constructing QSAR based on some known mode of receptor/ligand interactions. Many disadvantages of traditional models can be avoided by using the proposed method because the traditional models only determined upon molecular structural features in sample sets themselves. A genetic virtual screening of peptide/protein combinations (GVSPPC) is proposed for the first time by utilizing this idea to examine peptide/protein affinity activities. A genetic algorithm (GA) was developed for screening combinative targets with an interaction mode for virtual receptors. GVSPPC succeeds in disposing difficulties in rational QSAR,in order to search for the ligand/receptor interactions on conditions of unknown structures. Some bioactive oligo-/poly-peptide systems covering 58 angiotensin converting enzyme (ACE) inhibitors and 18 double site mutation residues in camel antibody protein cAb-Lys3 were investigated by GVSPPC with satisfactory results (R 2 cu>0.91,Q 2 cv > 0.86,ERMS=0.19-0.95),respectively,which demonstrates that GVSPPC is more inter-pretable in the ligand-receptor interaction than the traditional QSAR method.     

A label-free MALDI TOF MS-based method for studying the kinetics and inhibitor screening of the Alzheimer’s disease drug target β-secretase

Analytical and Bioanalytical Chemistry - Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI TOF MS) is a well-established method with a unique set of qualities...     

Xphos ligand and platinum catalysts: A versatile catalyst for the synthesis of functionalized β-(E)-vinylsilanes from terminal alkynes

Hydrosilylation of functionalized terminal arylalkynes with a variety of silanes catalyzed by PtCl₂ or PtO₂ in the presence of the air-stable and bulky Xphos ligand was investigated. Regardless of the electronic nature (electron withdrawing or donating group) and the position (o, m, p) of the substituents on the aromatic ring, a single β-(E)-styrylsilanes was obtained in good to excellent yields. The regioselectivity of the H-Si bond addition was found to be governed by steric effects induced by the bulky Xphos ligand. A dramatic regioselectivity was also observed when functionalized terminal aliphatic alkynes were employed as a substrate and in these cases regioisomeric β-(E)-vinylsilanes were generated with excellent selectivity.     

Preparation,characterization and catalytic activity of MoO3/CeO2–ZrO2 solid heterogeneous catalyst for the synthesis of β-enaminones

A series of CeO₂–ZrO₂ with different molybdenum (8–20 wt% MoO₃) loaded materials were prepared by homogeneous co-precipitation followed by impregnation method. The prepared materials were tested for their catalytic activity performance in the synthesis of β-enaminones by condensation of various anilines with dimedone under solvent-free conditions in microwave providing excellent yields within short reaction time. An obtained result reveals that, catalytic activity increases with increase in Mo wt% loading. The particle size of prepared materials was estimated using Debye–Scherrer equation. Particle size increases with increase in Mo wt% loading providing nanosized particle ranging from 7.11 to 42.09 nm. The synthesized materials were characterized by means of X-ray powder diffraction (XRD), Fourier-transform infrared spectroscopy (FT-IR), scanning electron microscopy (SEM) and energy dispersive spectroscopy (EDS) techniques.     

3D-QSAR methods on the basis of ligand–receptor complexes. Application of COMBINE and GRID/GOLPE methodologies to a series of CYP1A2 ligands

Many heterocyclic amines (HCA) present in cooked food exert a genotoxic activity when they are metabolised (N-oxidated) by the human cytochrome P450 1A2 (CYP1A2h). In order to rationalize the observed differences in activity of this enzyme on a series of 12 HCA, 3D-QSAR methods were applied on the basis of models of HCA–CYP1A2h complexes. The CYP1A2h enzyme model has been previously reported and was built by homology modeling based on cytochrome P450 BM3. The complexes were automatically generated applying the AUTODOCK software and refined using AMBER. A COMBINE analysis on the complexes identified the most important enzyme–ligand interactions that account for the differences in activity within the series. A GRID/GOLPE analysis was then performed on just the ligands, in the conformations and orientations found in the modeled complexes. The results from both methods were concordant and confirmed the advantages of incorporating structural information from series of ligand–receptor complexes into 3D-QSAR methodologies.     

Comparing the peroxo/superoxo nature of the interaction between molecular O2 and β-diketiminato-copper and nickel complexes

Adducts resulting from the interaction between molecular oxygen and β-diketiminato-copper and nickel complexes have been recently described in the literature as peroxo and superoxo complexes, respectively. The nature of the interaction is analyzed by means of DDCI calculations and an orthogonal valence bond reading of the ground state wavefunction for each system. Our results reveal that there is not any substantial difference between these systems, both presenting a marked leading superoxo nature, which is in line with the fact that LCu-O(2) and LNi-O(2) adducts present similar O-O distances and quite close O-O stretching vibration modes.     

Combination of β-elimination and liquid chromatography/quadrupole time-of-flight mass spectrometry for the determination of O-glycosylation sites

Determination of O-glycosylation sites in glycopeptides was developed by using two model compounds designed from mucin2 tandem repeat motif and erythropoietin. β-Elimination/addition reaction using dimethylamine on glycosylated site through a Michael-type condensation produced efficient deglycosylation with appropriate chemical modification. The use of dimethylamine was efficient to release the O-linked glycan in a reaction time period of 2-6 h at 55 °C. Peptide sequencing was then performed using the liquid chromatography/quadrupole time-of-flight mass spectrometry and MS-MS experiments. Interpretation of fragmentation pathways of the β-elimination/addition products enabled straightforward recognition of glycosylation site. Compared to the fragmentation of corresponding native peptides, mass shift of −18 Da or +27 Da was clearly observed for the two kinds of β-elimination/addition products of the glycosylated threonine. Dimethylamine was found to provide higher efficiency of β-elimination/addition than methylamine and ammonia.     

Effect of the β-diketonate ligand on the spin states of [Ni(β-dkt)2(NH2-quin)] complexes

Three new complexes [Ni(β-dkt)₂(NH₂-quin)] {β-dkt = 2,2,6,6-tetramethylheptane-3,3-dionate (tmhd) 1, hexafluoroacetylacetonate (hfac) 2, 1,3-diphenylpropanedionate (dbm) 3} have been prepared by reacting [Ni(β-dkt)₂(H₂O)₂] with 8-aminoquinoline (NH₂-quin). [Ni(tmhd)₂(NH₂-quin)] is found to be solvatochromic exhibiting a square planar geometry in CH₂Cl₂, acetone, and acetonitrile with the tmhd ligand acting unusually as a counteranion while in THF and DMSO the complex is octahedral. In contrast, 2 and 3 are octahedral in all solvents. Single crystal X-ray diffraction studies reveal octahedral nickel centres with a cis arrangement of the β-diketonates. The molecular packing consists of hydrogen bonded dimers which in the case of 2 and 3 are connected to one another via π···π and C-H?π interactions, respectively. Cyclic voltammetry shows 1 and 3 oxidise irreversibly at 0.46 and 1.17 V, respectively.