Molecular dynamics simulations reveal structural insights into inhibitor binding modes and functionality in human Group IIA phospholipase A2

Human Group IIA phospholipase A₂ (hGIIA) promotes inflammation in immune-mediated pathologies by regulating the arachidonic acid pathway through both catalysis-dependent and -independent mechanisms. The hGIIA crystal structure, both alone and inhibitor-bound, together with structures of closely related snake-venom-derived secreted phospholipase enzymes has been well described. However, differentiation of biological and nonbiological contacts and the relevance of structures determined from snake venom enzymes to human enzymes are not clear. We employed molecular dynamics (MD) and docking approaches to understand the binding of inhibitors that selectively or nonselectively block the catalysis-independent mechanism of hGIIA. Our results indicate that hGIIA behaves as a monomer in the solution environment rather than a dimer arrangement that is in the asymmetric unit of some crystal structures. The binding mode of a nonselective inhibitor, KH064, was validated by a combination of the experimental electron density and MD simulations. The binding mode of the selective pentapeptide inhibitor FLSYK to hGIIA was stipulated to be different to that of the snake venom phospholipases A₂ of Daboia russelli pulchella (svPLA₂ ). Our data suggest that the application of MD approaches to crystal structure data is beneficial in evaluating the robustness of conclusions drawn based on crystal structure data alone. Proteins 2017; 85:827-842. © 2016 Wiley Periodicals, Inc.

Discovery of Rab1 binding sites using an ensemble of clustering methods

Targeting non-native-ligand binding sites for potential investigative and therapeutic applications is an attractive strategy in proteins that share common native ligands, as in Rab1 protein. Rab1 is a subfamily member of Rab proteins, which are members of Ras GTPase superfamily. All Ras GTPase superfamily members bind to native ligands GTP and GDP, that switch on and off the proteins, respectively. Rab1 is physiologically essential for autophagy and transport between endoplasmic reticulum and Golgi apparatus. Pathologically, Rab1 is implicated in human cancers, a neurodegenerative disease, cardiomyopathy, and bacteria-caused infectious diseases. We have performed structural analyses on Rab1 protein using a unique ensemble of clustering methods, including multi-step principal component analysis, non-negative matrix factorization, and independent component analysis, to better identify representative Rab1 proteins than the application of a single clustering method alone does. We then used the identified representative Rab1 structures, resolved in multiple ligand states, to map their known and novel binding sites. We report here at least a novel binding site on Rab1, involving Rab1-specific residues that could be further explored for the rational design and development of investigative probes and/or therapeutic small molecules against the Rab1 protein. Proteins 2017; 85:859-871. © 2016 Wiley Periodicals, Inc.

A Brief Review of RNA-Protein Interaction Database Resources

RNA–Protein interactions play critical roles in various biological processes. By collecting and analyzing the RNA–Protein interactions and binding sites from experiments and predictions, RNA–Protein interaction databases have become an essential resource for the exploration of the transcriptional and post-transcriptional regulatory network. Here, we briefly review several widely used RNA–Protein interaction database resources developed in recent years to provide a guide of these databases. The content and major functions in databases are presented. The brief description of database helps users to quickly choose the database containing information they interested. In short, these RNA–Protein interaction database resources are continually updated, but the current state shows the efforts to identify and analyze the large amount of RNA–Protein interactions.

A Machine Learning Approach to Predict Gene Regulatory Networks in Seed Development in Arabidopsis

Gene regulatory networks (GRNs) provide a representation of relationships between regulators and their target genes. Several methods for GRN inference, both unsupervised and supervised, have been developed to date. Because regulatory relationships consistently reprogram in diverse tissues or under different conditions, GRNs inferred without specific biological contexts are of limited applicability. In this report, a machine learning approach is presented to predict GRNs specific to developing Arabidopsis thaliana embryos. We developed the Beacon GRN inference tool to predict GRNs occurring during seed development in Arabidopsis based on a support vector machine (SVM) model. We developed both global and local inference models and compared their performance, demonstrating that local models are generally superior for our application. Using both the expression levels of the genes expressed in developing embryos and prior known regulatory relationships, GRNs were predicted for specific embryonic developmental stages. The targets that are strongly positively correlated with their regulators are mostly expressed at the beginning of seed development. Potential direct targets were identified based on a match between the promoter regions of these inferred targets and the cis elements recognized by specific regulators. Our analysis also provides evidence for previously unknown inhibitory effects of three positive regulators of gene expression. The Beacon GRN inference tool provides a valuable model system for context-specific GRN inference and is freely available at https://github.com/BeaconProjectAtVirginiaTech/beacon_network_inference.git.

Protein social behavior makes a stronger signal for partner identification than surface geometry

Cells are interactive living systems where proteins movements, interactions and regulation are substantially free from centralized management. How protein physico‐chemical and geometrical properties determine who interact with whom remains far from fully understood. We show that characterizing how a protein behaves with many potential interactors in a complete cross‐docking study leads to a sharp identification of its cellular/true/native partner(s). We define a sociability index, or S‐index, reflecting whether a protein likes or not to pair with other proteins. Formally, we propose a suitable normalization function that accounts for protein sociability and we combine it with a simple interface‐based (ranking) score to discriminate partners from non‐interactors. We show that sociability is an important factor and that the normalization permits to reach a much higher discriminative power than shape complementarity docking scores. The social effect is also observed with more sophisticated docking algorithms. Docking conformations are evaluated using experimental binding sites. These latter approximate in the best possible way binding sites predictions, which have reached high accuracy in recent years. This makes our analysis helpful for a global understanding of partner identification and for suggesting discriminating strategies. These results contradict previous findings claiming the partner identification problem being solvable solely with geometrical docking. Proteins 2016; 85:137–154. © 2016 Wiley Periodicals, Inc.

An Empirical Prior Improves Accuracy for Bayesian Estimation of Transcription Factor Binding Site Frequencies within Gene Promoters

A Bayesian method for sampling from the distribution of matches to a precompiled transcription factor binding site (TFBS) sequence pattern (conditioned on an observed nucleotide sequence and the sequence pattern) is described. The method takes a position frequency matrix as input for a set of representative binding sites for a transcription factor and two sets of noncoding, 5′ regulatory sequences for gene sets that are to be compared. An empirical prior on the frequency A (per base pair of gene-vicinal, noncoding DNA) of TFBSs is developed using data from the ENCODE project and incorporated into the method. In addition, a probabilistic model for binding site occurrences conditioned on λ is developed analytically, taking into account the finite-width effects of binding sites. The count of TFBS β (conditioned on the observed sequence) is sampled using Metropolis–Hastings with an information entropy-based move generator. The derivation of the method is presented in a step-by-step fashion, starting from specific conditional independence assumptions. Empirical results show that the newly proposed prior on β improves accuracy for estimating the number of TFBS within a set of promoter sequences.

Copper Ion Binding Site in β-Amyloid Peptide

β-Amyloid aggregates in the brain play critical roles in Alzheimer's disease, a chronic neurodegenerative condition. Amyloid-associated metal ions, particularly zinc and copper ions, have been implicated in disease pathogenesis. Despite the importance of such ions, the binding sites on the β-amyloid peptide remain poorly understood. In this study, we use scanning tunneling microscopy, circular dichroism, and surface-enhanced Raman spectroscopy to probe the interactions between Cu²⁺ ions and a key β-amyloid peptide fragment, consisting of the first 16 amino acids, and define the copper-peptide binding site. We observe that in the presence of Cu²⁺, this peptide fragment forms β-sheets, not seen without the metal ion. By imaging with scanning tunneling microscopy, we are able to identify the binding site, which involves two histidine residues, His13 and His14. We conclude that the binding of copper to these residues creates an interstrand histidine brace, which enables the formation of β-sheets.

The use of fake ligands from computational solvent mapping in ligand and structure-based virtual screening

AIM

Virtual screening selects compounds that resemble a known modulator or compounds that fit into the binding site of a target protein. Computational solvent mapping defines important chemical features for binding to a target protein. Results/methodology: We have tested the ability to use solvent mapping for generating a 'fake' ligand that is a negative image of the binding site. We used this fake ligand as a query for the program ROCS and to define the search space of the docking programs FRED and HYBRID.

CONCLUSION

The fake ligands perform comparably to or better than the ligands from crystal structures across a set of ten targets. Thus, the approach is suitable for guiding virtual screening and hit-to-lead optimization.

Nucleotide Interdependency in Transcription Factor Binding Sites in the Drosophila Genome

A long-standing objective in modern biology is to characterize the molecular components that drive the development of an organism. At the heart of eukaryotic development lies gene regulation. On the molecular level, much of the research in this field has focused on the binding of transcription factors (TFs) to regulatory regions in the genome known as cis-regulatory modules (CRMs). However, relatively little is known about the sequence-specific binding preferences of many TFs, especially with respect to the possible interdependencies between the nucleotides that make up binding sites. A particular limitation of many existing algorithms that aim to predict binding site sequences is that they do not allow for dependencies between nonadjacent nucleotides. In this study, we use a recently developed computational algorithm, MARZ, to compare binding site sequences using 32 distinct models in a systematic and unbiased approach to explore nucleotide dependencies within binding sites for 15 distinct TFs known to be critical to Drosophila development. Our results indicate that many of these proteins have varying levels of nucleotide interdependencies within their DNA recognition sequences, and that, in some cases, models that account for these dependencies greatly outperform traditional models that are used to predict binding sites. We also directly compare the ability of different models to identify the known KRUPPEL TF binding sites in CRMs and demonstrate that a more complex model that accounts for nucleotide interdependencies performs better when compared with simple models. This ability to identify TFs with critical nucleotide interdependencies in their binding sites will lead to a deeper understanding of how these molecular characteristics contribute to the architecture of CRMs and the precise regulation of transcription during organismal development.

Tyrosine411 and Arginine410 of Human Serum Albumin Play an Important Role in the Binding of Sodium 4-Phenylbutyrate to Site II

Sodium 4-phenylbutyrate (PB) has many pharmacological activities; therefore extending its clinical use to the treatment of a wider variety of diseases would be desirable. However, our knowledge of the binding of PB to plasma proteins is not extensive. To address this issue in more detail, we characterized the protein binding of PB. Binding experiments showed that PB mainly binds to human serum albumin (HSA) in plasma. PB was also found to bind to a single site on HSA, which was identified as site II by fluorescent probe displacement experiment. Furthermore, an appropriate alkyl chain length and a carboxylic group in the PB structure were required for PB binding to HSA, suggesting that hydrophobic (and van der Waals) and electrostatic interactions are involved as binding modes. The contributions of hydrogen bonding and/or van der Waals interactions were also indicated by thermodynamic analyses. Tyrosine411 and arginine410 were identified as being involved in the binding of PB to site II, based on binding experiments using chemically modified- and mutant-HSA preparations. In conclusion, the available evidence indicates that PB binds to site II of HSA with assistance by multiple forces and that tyrosine411 and arginine410 both play important roles in this phenomenon.

Understanding the Robust Physisorption between Bovine Serum Albumin and Amphiphilic Polymer Coated Nanoparticles

The robust physisorption between nanoparticles (NPs) and proteins has attracted increasing attention due to the significance for both conjugation techniques and protein's corona formation at the bionano interface. In the present study, we first explored the possible binding sites of the bovine serum albumin (BSA) on amphiphilic polymer coated gold nanoparticles (AP-AuNPs). By using mass spectrometry, a 105-amino-acid peptide (12.2 kDa) is discovered as the possible "epitope" responsible for the robust physisorption between BSA and AP-AuNPs. Second, with the help of nanometal surface energy transfer (NSET) theory, we further found that the epitope peptide could insert at least 2.9 nm into the organic molecular layers of AP-AuNPs when the robust conjugates formed, which indicates how such a long epitope peptide can be accommodated by AP-AuNPs and resist protease's digestion. These findings might shed light on a new strategy for studying interactions between proteins and NPs, and further guide the rational design of NPs for safe and effective biomedical applications.

Identification of the Allosteric Site for Phenylalanine in Rat Phenylalanine Hydroxylase

Liver phenylalanine hydroxylase (PheH) is an allosteric enzyme that requires activation by phenylalanine for full activity. The location of the allosteric site for phenylalanine has not been established. NMR spectroscopy of the isolated regulatory domain (RDPheH(25-117) is the regulatory domain of PheH lacking residues 1-24) of the rat enzyme in the presence of phenylalanine is consistent with formation of a side-by-side ACT dimer. Six residues in RDPheH(25-117) were identified as being in the phenylalanine-binding site on the basis of intermolecular NOEs between unlabeled phenylalanine and isotopically labeled protein. The location of these residues is consistent with two allosteric sites per dimer, with each site containing residues from both monomers. Site-specific variants of five of the residues (E44Q, A47G, L48V, L62V, and H64N) decreased the affinity of RDPheH(25-117) for phenylalanine based on the ability to stabilize the dimer. Incorporation of the A47G, L48V, and H64N mutations into the intact protein increased the concentration of phenylalanine required for activation. The results identify the location of the allosteric site as the interface of the regulatory domain dimer formed in activated PheH.

[Web server for prediction of miRNAs and their precursors and binding sites]

A microRNA (miRNA) is a small noncoding RNA molecule about 22 nucleotides in length. The paper describes a web server for predicting miRNAs and their precursors and binding sites. The predictions are based on either sequence similarity to known miRNAs of 223 organisms or context-structural hidden Markov models. It has been shown that the proposed methods of prediction of human miRNAs and pre-miRNAs outperform the existing ones in accuracy. The average deviation of predicted 5'-ends of human miRNAs from actual positions is 3.13 nt in the case of predicting one pair of complementary miRNAs (miRNA-miRNA* duplex). A useful option for our application is the prediction of an additional miRNA pair. In this mode, the pairs closest to actual miRNA deviate by 1.61 nt on average. The proposed method also shows good performance in predicting mouse miRNAs. Binding sites for miRNAs are predicted by two known approaches based on complementarity and thermodynamic stability of the miRNA-mRNA duplex and on a new approach, which takes into account miRNAs competition for the site. The role of the secondary structure in miRNA processing is considered. The web server is available at http://wwwmgs.bionet.nsc.ru/mgs/programs/rnaanalys/.

The therapeutic effect of chlorogenic acid against Staphylococcus aureus infection through sortase A inhibition

The emergence and wide spread of multi-drug resistant Staphylococcus aureus (S. aureus) requires the development of new therapeutic agents with alternative modes of action. Anti-virulence strategies are hoped to meet that need. Sortase A (SrtA) has attracted great interest as a potential drug target to treat infections caused by S. aureus, as many of the surface proteins displayed by SrtA function as virulence factors by mediating bacterial adhesion to specific organ tissues, invasion of host cells, and evasion of the host-immune responses. It has been suggested that inhibitors of SrtA might be promising candidates for the treatment and/or prevention of S. aureus infections. In this study, we report that chlorogenic acid (CHA), a natural compound that lacks significant anti-S. aureus activity, inhibit the activity of SrtA in vitro (IC50 = 33.86 ± 5.55 μg/ml) and the binding of S. aureus to fibrinogen (Fg). Using molecular dynamics simulations and mutagenesis assays, we further demonstrate that CHA binds to the binding sites of C184 and G192 in the SrtA. In vivo studies demonstrated that CHA prevent mice from S. aureus-induced renal abscess, resulting in a significant survival advantage. These findings indicate that CHA is a promising therapeutic compound against SrtA during S. aureus infections.

Identification of several potential chromatin binding sites of HOXB7 and its downstream target genes in breast cancer

HOXB7 encodes a transcription factor that is overexpressed in a number of cancers and encompasses many oncogenic functions. Previous results have shown it to promote cell proliferation, angiogenesis, epithelial–mesenchymal transition, DNA repair and cell survival. Because of its role in many cancers and tumorigenic processes, HOXB7 has been suggested to be a potential drug target. However, HOXB7 binding sites on chromatin and its targets are poorly known. The aim of our study was to identify HOXB7 binding sites on breast cancer cell chromatin and to delineate direct target genes located nearby these binding sites. We found 1,504 HOXB7 chromatin binding sites in BT‐474 breast cancer cell line that overexpresses HOXB7. Seventeen selected binding sites were validated by ChIP‐qPCR in several breast cancer cell lines. Furthermore, we analyzed expression of a large number of genes located nearby HOXB7 binding sites and found several new direct targets, such as CTNND2 and SCGB1D2. Identification of HOXB7 chromatin binding sites and target genes is essential to understand better the role of HOXB7 in breast cancer and mechanisms by which it regulates tumorigenic processes.

What's new?

The transcription factor HOXB7 is overexpressed in various cancers, but it's not yet known just which genes HOXB7 activates. How does it influence cancer on a molecular level? This study found 1500 sequences where HOXB7 binds the chromatin in breast cancer cells. They went on to identify several potential target genes near the HOXB7 binding sites. Not only will finding these genes help explain how HOXB7 overexpression promotes tumor growth, it will help understand what side effects might result from hindering HOXB7 expression.

Probing interactions between human lung adenocarcinoma A549 cell and its aptamers at single-molecule resolution

Because cell-specific aptamers have high potential for biomedical applications, investigation of the interaction between cell and its aptamers may be of key importance for an improved understanding of biochemical processes. Herein, the interaction between human lung adenocarcinoma A549 cell and its four aptamers was explored using single-molecule force spectroscopy (SMFS). The values of the unbinding force varied from 117.1 to 171.0 pN at the loading rate of 1.8 × 10(5)  pN/s. Based on the dependence of singe molecule force on the atomic force microscopy loading rate, the corresponding kinetic parameters were obtained. The results revealed two activation barriers and two transient states in the unbinding process of aptamer/cell interaction. More importantly, the binding sites on A549 cells with its four aptamers were defined to be different using SMFS and flow cytometry. This work demonstrated that SMFS can be used as a powerful tool for exploring the aptamer/cell binding behavior at the single-molecule level, and may provide valuable information for the design and application of aptamer probes.

Enhancement of the binding affinity of methylene blue to site I in human serum albumin by cupric and ferric ions

In this work, the binding characteristics of methylene blue (MB) to human serum albumin (HSA) and the influence of Cu(2+) and Fe(3+) on the binding affinity of MB to HSA were investigated using fluorescence, absorption, circular dichroism (CD) spectroscopy and molecular modelling. The results of competitive binding experiments using the site probes ketoprofen and ibuprofen as specific markers suggested that MB was located in site I within sub-domain IIA of HSA. The molecular modelling results agreed with the results of competitive site marker experiments and the results of CD spectra indicated that the interaction between MB and HSA caused the conformational changes in HSA. The binding affinity of MB to HSA was enhanced but to a different extent in the presence of Cu(2+) and Fe(3+), respectively, which indicated that the influence of different metal ions varied. Enhancement of the binding affinity of MB to HSA in the presence of Cu(2+) is due to the formation of Cu(2+)-HSA complex leading to the conformational changes in HSA, whereas in the presence of Fe(3+), enhancement of the binding affinity is due to the greater stability of the Fe(3+)-HSA-MB complex compared with the MB-HSA complex.

Competitive mode and site of interaction of ticagrelor at the human platelet P2Y12 -receptor

BACKGROUND

The G-protein-coupled P2Y12 -receptor plays a crucial role in platelet aggregation. Recently, ticagrelor was licensed as the first perorally active and reversible P2Y12 -receptor antagonist.

OBJECTIVE

The present study investigated the site and the antagonistic mode of action of ticagrelor at wild-type or mutant human P2Y12 -receptors.

METHODS

Recombinant wild-type or mutant human P2Y12 -receptors were stably expressed in Chinese hamster ovary Flp-In cells. Receptor function was assessed by quantification of ADP- and 2-methylthio-ADP-mediated inhibition of forskolin-induced cellular cAMP production either using a [(3) H]cAMP-radioaffinity assay or a cAMP response element-driven luciferase reporter gene assay.

RESULTS

The natural agonist ADP inhibited forskolin-induced cAMP formation at the wild-type P2Y12 -receptor with a lower potency (EC50 209 nm) than the synthetic agonist 2-methylthio-ADP (EC50 1.0 nm). Ticagrelor shifted the concentration-response curves of both agonists in a parallel and surmountable manner to the right. Increasing concentrations of ticagrelor caused increasing shifts. Schild-plot analysis revealed pA2 values of 8.85 for ticagrelor against ADP, and 8.69 against 2-methylthio-ADP, and slopes of the regression lines not different from unity. In cells expressing a recombinant C194A(5.43) -mutant P2Y12 -receptor construct, ticagrelor lost antagonistic potency when tested against ADP or 2-methylthio-ADP.

CONCLUSIONS

The experiments reveal a surmountable and competitive mode of antagonism of ticagrelor at P2Y12 -receptors activated by either the natural agonist ADP or the synthetic agonist 2-methylthio-ADP. Cys194(5.43) is likely to be involved in the interaction of ticagrelor with ADP and 2-methylthio-ADP. The data give new insights into the site and mode of action of ticagrelor at the human P2Y12 -receptor.

Assessment of ligand binding site predictions in CASP10

The identification of amino acid residues in proteins involved in binding small molecule ligands is an important step for their functional characterization, as the function of a protein often depends on specific interactions with other molecules. The accuracy of computational methods aiming to predict such binding residues was evaluated within the "function prediction (prediction of binding sites, FN)" category of the critical assessment of protein structure prediction (CASP) experiment. In the last edition of the experiment (CASP10), 17 research groups participated in this category, and their predictions were evaluated on 13 prediction targets containing biologically relevant ligands. The results of this experiment indicate that several methods achieved an overall good performance, showing the usefulness of such methods in predicting ligand binding residues. As in previous years, methods based on a homology transfer approach were dominating. In comparison to CASP9, a larger fraction of the top predictors are automated servers. However, due to the small number of targets and the characteristics of the prediction format, the differences observed among the first ten methods were not statistically significant and it was also not possible to analyze differences in accuracy for different ligand types or overall structure, difficulty. To overcome these limitations and to allow for a more detailed evaluation, in future editions of CASP, methods in the FN category will no longer be evaluated on the "normal" CASP targets, but assessed continuously by CAMEO (continuous automated model evaluation) based on weekly prereleased sequences from the PDB.

Pathway of binding of the intrinsically disordered mitochondrial inhibitor protein to F1-ATPase

The hydrolysis of ATP by the ATP synthase in mitochondria is inhibited by a protein called IF1. Bovine IF1 has 84 amino acids, and its N-terminal inhibitory region is intrinsically disordered. In a known structure of bovine F1-ATPase inhibited with residues 1-60 of IF1, the inhibitory region from residues 1-50 is mainly α-helical and buried deeply at the α(DP)β(DP)-catalytic interface, where it forms extensive interactions with five of the nine subunits of F1-ATPase but mainly with the β(DP)-subunit. As described here, on the basis of two structures of inhibited complexes formed in the presence of large molar excesses of residues 1-60 of IF1 and of a version of IF1 with the mutation K39A, it appears that the intrinsically disordered inhibitory region interacts first with the αEβE-catalytic interface, the most open of the three catalytic interfaces, where the available interactions with the enzyme allow it to form an α-helix from residues 31-49. Then, in response to the hydrolysis of an ATP molecule and the associated partial closure of the interface to the αTPβTP state, the extent of the folded α-helical region of IF1 increases to residues 23-50 as more interactions with the enzyme become possible. Finally, in response to the hydrolysis of a second ATP molecule and a concomitant 120° rotation of the γ-subunit, the interface closes further to the α(DP)β(DP)-state, allowing more interactions to form between the enzyme and IF1. The structure of IF1 now extends to its maximally folded state found in the previously observed inhibited complex.

Binding site multiplicity with fatty acid ligands: implications for the regulation of PKR kinase autophosphorylation with palmitate

Saturated long chain-free fatty acids (FFAs), especially palmitate, have been implicated in apoptosis by inhibiting the activity of PKR (double-stranded RNA-dependent protein kinase). We recently found evidence that palmitate interacts directly with the kinase domain of PKR, subsequently inhibiting the autophosphorylation of PKR. To investigate the interactions of palmitate with PKR and its effects on PKR autophosphorylation, we performed extensive unbiased MD simulations combined with biochemical and biophysical experiments. The simulations predict multiple putative binding sites of palmitate on both the phosphorylated and unphosphorylated PKR with similar binding affinities. Ligand-protein interactions involving a large variety of different binding modes challenge the conventional view of highly specific, single binding sites. Key interactions of palmitate involve the αC-helix of PKR, especially near residue R307. Experimental mutation of R307 was found to affect palmitate binding and reduce its inhibitory effect. Based on this study a new allosteric mechanism is proposed where palmitate binding to the αC-helix prevents the inactive-to-active transition of PKR and subsequently reduces its ability to autophosphorylate.

Related genes and potential biomarkers for early diagnosis of Alzheimer's disease: a preliminary study based on DNA microarray

AIM

The aim of this study is to extend our understanding of the molecular mechanism of Alzheimer's disease (AD).

METHODS

We downloaded the gene expression profile GSE18309 from Gene Expression Omnibus database, which includes 3 genechips from patients with mild cognitive impairment (MCI), 3 genechips from patients with AD, and 3 genechips from normal controls (NC). Linear Models for Microarray Data package was used to identify differentially expressed genes (DEGs) in MCI versus NC group and AD versus NC group. Then, we extracted the overlapping DEGs of 2 groups for functional and pathway enrichment analysis using FuncAssociate software accompanied by gene ontology and expressing analysis systematic explorer, respectively. Further, AutoDock4 (http://autodock.scripps.edu/) was used to predict the docking site between small molecule ligands and proteins of a key DEG.

RESULTS

A total of 60 DEGs were identified. Biological processes associated with nutrient response and muscle development were significantly dysregulated in AD and MCI. In addition, we identified 2 active binding sites (A5 and L30) on protein structure of cholecystokinin A receptor (CCKAR) for drug design.

CONCLUSION

The DEGs including CCKAR might be used as biomarkers for early diagnosis of AD. However, further experimental studies are needed to confirm our results.

Exploring the ATP-binding site of P2X receptors

P2X receptors are ATP-gated non-selective cation channels involved in many different physiological processes, such as synaptic transmission, inflammation, and neuropathic pain. They form homo- or heterotrimeric complexes and contain three ATP-binding sites in their extracellular domain. The recent determination of X-ray structures of a P2X receptor solved in two states, a resting closed state and an ATP-bound, open-channel state, has provided unprecedented information not only regarding the three-dimensional shape of the receptor, but also on putative conformational changes that couple ATP binding to channel opening. These data provide a structural template for interpreting the huge amount of functional, mutagenesis, and biochemical data collected during more than fifteen years. In particular, the interfacial location of the ATP binding site and ATP orientation have been successfully confirmed by these structural studies. It appears that ATP binds to inter-subunit cavities shaped like open jaws, whose tightening induces the opening of the ion channel. These structural data thus represent a firm basis for understanding the activation mechanism of P2X receptors.

Analysis of the binding sites of porcine sialoadhesin receptor with PRRSV

Porcine reproductive and respiratory syndrome virus (PRRSV) can infect pigs and cause enormous economic losses to the pig industry worldwide. Porcine sialoadhesin (pSN) and CD163 have been identified as key viral receptors on porcine alveolar macrophages (PAM), a main target cell infected by PRRSV. In this study, the protein structures of amino acids 1-119 from the pSN and cSN (cattle sialoadhesin) N-termini (excluding the 19-amino acid signal peptide) were modeled via homology modeling based on mSN (mouse sialoadhesin) template structures using bioinformatics tools. Subsequently, pSN and cSN homology structures were superposed onto the mSN protein structure to predict the binding sites of pSN. As a validation experiment, the SN N-terminus (including the wild-type and site-directed-mutant-types of pSN and cSN) was cloned and expressed as a SN-GFP chimera protein. The binding activity between SN and PRRSV was confirmed by WB (Western blotting), FAR-WB (far Western blotting), ELISA (enzyme-linked immunosorbent assay) and immunofluorescence assay. We found that the S107 amino acid residue in the pSN N-terminal played a crucial role in forming a special cavity, as well as a hydrogen bond for enhancing PRRSV binding during PRRSV infection. S107 may be glycosylated during PRRSV infection and may also be involved in forming the cavity for binding PRRSV along with other sites, including W2, Y44, S45, R97, R105, W106 and V109. Additionally, S107 might also be important for pSN binding with PRRSV. However, the function of these binding sites must be confirmed by further studies.

Protein function annotation by local binding site surface similarity

Hundreds of protein crystal structures exist for proteins whose function cannot be confidently determined from sequence similarity. Surflex-PSIM, a previously reported surface-based protein similarity algorithm, provides an alternative method for hypothesizing function for such proteins. The method now supports fully automatic binding site detection and is fast enough to screen comprehensive databases of protein binding sites. The binding site detection methodology was validated on apo/holo cognate protein pairs, correctly identifying 91% of ligand binding sites in holo structures and 88% in apo structures where corresponding sites existed. For correctly detected apo binding sites, the cognate holo site was the most similar binding site 87% of the time. PSIM was used to screen a set of proteins that had poorly characterized functions at the time of crystallization, but were later biochemically annotated. Using a fully automated protocol, this set of 8 proteins was screened against ∼60,000 ligand binding sites from the PDB. PSIM correctly identified functional matches that predated query protein biochemical annotation for five out of the eight query proteins. A panel of 12 currently unannotated proteins was also screened, resulting in a large number of statistically significant binding site matches, some of which suggest likely functions for the poorly characterized proteins.

The androgen receptor transcriptional program in castration-resistant prostate cancer: cell lines vs. tissue samples

The androgen receptor (AR) is central to the initiation and progression of prostate cancer, even after castration. Its transcriptional activity has previously been studied in cell lines. A group at the University of Cambridge recently outlined the AR transcriptional program in tissue samples, with an emphasis on castration-resistant tumors. AR binding sites, gene-expression changes (in xenografts), and potential transcription factor interactions were notably different from those observed in cultured cells. These discrepancies suggest a distinct signaling network for the AR in vivo and serve as a reminder that results from in vitro models should be checked against clinical realities.

GABA(A) receptor transmembrane amino acids are critical for alcohol action: disulfide cross-linking and alkyl methanethiosulfonate labeling reveal relative location of binding sites

Alcohols and inhaled anesthetics modulate GABA(A) receptor (GABA(A)R) function via putative binding sites within the transmembrane regions. The relative position of the amino acids lining these sites could be either inter- or intra-subunit. We introduced cysteines in relevant TM locations and tested the proximity of cysteine pairs using oxidizing and reducing agents to induce or break disulfide bridges between cysteines, and thus change GABA-mediated currents in wild-type and mutant α1β2γ2 GABA(A)Rs expressed in Xenopus laevis oocytes. We tested for: (i) inter-subunit cross-linking: a cysteine located in α1TM1 [either α1(Q229C) or α1(L232C)] was paired with a cysteine in different positions of β2TM2 and TM3; (ii) intra-subunit cross-linking: a cysteine located either in β2TM1 [β2(T225C)] or in TM2 [β2(N265C)] was paired with a cysteine in different locations along β2TM3. Three inter-subunit cysteine pairs and four intra-subunits cross-linked. In three intra-subunit cysteine combinations, the alcohol effect was reduced by oxidizing agents, suggesting intra-subunit alcohol binding. We conclude that the structure of the alcohol binding site changes during activation and that potentiation or inhibition by binding at inter- or intra-subunit sites is determined by the specific receptor and ligand.

Spectroscopic and molecular modeling studies of the interaction between morin and polyamidoamine dendrimer

Interactions between the polyamidoamine (PAMAM) dendrimer and drug molecules are of interest for their potential biomedical applications. The goal of this work is to examine the interaction of PAMAM-C12 25% dendrimer with morin. The ultraviolet-visible, fluorescence spectroscopic methods as well as molecular modeling were used to analyze drug-binding mode, binding constants and binding sites, etc. The experimental data showed that the binding constant of morin-PAMAM-C12 25% is about 10(5) L/mol. The interaction of morin with PAMAM-C12 25% is mainly driven by the hydrophobic, electrostatic, hydrogen bonds and van der Waals forces. There are mainly three classes of binding site of morin at the interface of PAMAM-C12 25%. These results provided some useful information for self-assembling and disassembling the PAMAM dendrimer as well as efficient drug delivery and therapeutic applications.

Promiscuity of the euonymus carbohydrate-binding domain

Plants synthesize small amounts of carbohydrate-binding proteins on exposure to stress. For example, on exposure to drought, high salt, wounding and by treatment with some plant hormones or by pathogen attack. In contrast to the 'classical' plant lectins that are mostly located in the vacuolar compartment, this new class of inducible lectins is present in the cytoplasm and in the nucleus. Taking into account that any physiological role of plant lectins most likely relies on their specific carbohydrate-binding activity and specificity, the discovery of these stress-related lectins provides strong evidence for the importance of protein-carbohydrate-interactions in plant cells. Hitherto, six families of such nucleocytoplasmic lectins have been identified in plants. This review will focus on the nucleocytoplasmic lectins with one or more Euonymus lectin (EUL) domain(s). The carbohydrate-binding specificity of EUL proteins from a monocot, a dicot and a lower plant has been compared. Furthermore, modeling of the different EUL domains revealed a similar ß-trefoil fold consisting of three bundles of ß-sheet organized around a pseudo three-fold symmetry axis. Despite the sequence similarity and the conserved amino acids in the binding site, glycan array analyses showed that the EUL domain has a promiscuous carbohydrate-binding site capable of accommodating high mannose N-glycans, blood group B related structures and galactosylated epitopes.

Insights to the evolution of Nucleobase-Ascorbate Transporters (NAT/NCS2 family) from the Cys-scanning analysis of xanthine permease XanQ

The nucleobase-ascorbate transporter or nucleobase-cation symporter-2 (NAT/NCS2) family is one of the five known families of transporters that use nucleobases as their principal substrates and the only one that is evolutionarily conserved and widespread in all major taxa of organisms. The family is a typical paradigm of a group of related transporters for which conservation in sequence and overall structure correlates with high functional variations between homologs. Strikingly, the human homologs fail to recognize nucleobases or related cytotoxic compounds. This fact allows important biomedical perspectives for translation of structure-function knowledge on this family to the rational design of targeted antimicrobial purine-related drugs. To date, very few homologs have been characterized experimentally in detail and only two, the xanthine permease XanQ and the uric acid/xanthine permease UapA, have been studied extensively with site-directed mutagenesis. Recently, the high-resolution structure of a related homolog, the uracil permease UraA, has been solved for the first time with crystallography. In this review, I summarize current knowledge and emphasize how the systematic Cys-scanning mutagenesis of XanQ, in conjunction with existing biochemical and genetic evidence for UapA and the x-ray structure of UraA, allow insight on the structure-function and evolutionary relationships of this important group of transporters. The review is organized in three parts referring to (I) the theory of use of Cys-scanning approaches in the study of membrane transporter families, (II) the state of the art with experimental knowledge and current research on the NAT/NCS2 family, (III) the perspectives derived from the Cys-scanning analysis of XanQ.

Neurotoxins and their binding areas on voltage-gated sodium channels

Voltage-gated sodium channels (VGSCs) are large transmembrane proteins that conduct sodium ions across the membrane and by doing so they generate signals of communication between many kinds of tissues. They are responsible for the generation and propagation of action potentials in excitable cells, in close collaboration with other channels like potassium channels. Therefore, genetic defects in sodium channel genes can cause a wide variety of diseases, generally called “channelopathies.” The first insights into the mechanism of action potentials and the involvement of sodium channels originated from Hodgkin and Huxley for which they were awarded the Nobel Prize in 1963. These concepts still form the basis for understanding the function of VGSCs. When VGSCs sense a sufficient change in membrane potential, they are activated and consequently generate a massive influx of sodium ions. Immediately after, channels will start to inactivate and currents decrease. In the inactivated state, channels stay refractory for new stimuli and they must return to the closed state before being susceptible to a new depolarization. On the other hand, studies with neurotoxins like tetrodotoxin (TTX) and saxitoxin (STX) also contributed largely to our today’s understanding of the structure and function of ion channels and of VGSCs specifically. Moreover, neurotoxins acting on ion channels turned out to be valuable lead compounds in the development of new drugs for the enormous range of diseases in which ion channels are involved. A recent example of a synthetic neurotoxin that made it to the market is ziconotide (Prialt^®, Elan). The original peptide, ω-MVIIA, is derived from the cone snail Conus magus and now FDA/EMA-approved for the management of severe chronic pain by blocking the N-type voltage-gated calcium channels in pain fibers. This review focuses on the current status of research on neurotoxins acting on VGSC, their contribution to further unravel the structure and function of VGSC and their potential as novel lead compounds in drug development.

Binding of the inhibitor protein IF(1) to bovine F(1)-ATPase

In the structure of bovine F(1)-ATPase inhibited with residues 1-60 of the bovine inhibitor protein IF(1), the α-helical inhibitor interacts with five of the nine subunits of F(1)-ATPase. In order to understand the contributions of individual amino acid residues to this complex binding mode, N-terminal deletions and point mutations have been introduced, and the binding properties of each mutant inhibitor protein have been examined. The N-terminal region of IF(1) destabilizes the interaction of the inhibitor with F(1)-ATPase and may assist in removing the inhibitor from its binding site when F(1)F(o)-ATPase is making ATP. Binding energy is provided by hydrophobic interactions between residues in the long α-helix of IF(1) and the C-terminal domains of the β(DP)-subunit and β(TP)-subunit and a salt bridge between residue E30 in the inhibitor and residue R408 in the C-terminal domain of the β(DP)-subunit. Several conserved charged amino acids in the long α-helix of IF(1) are also required for establishing inhibitory activity, but in the final inhibited state, they are not in contact with F(1)-ATPase and occupy aqueous cavities in F(1)-ATPase. They probably participate in the pathway from the initial interaction of the inhibitor and the enzyme to the final inhibited complex observed in the structure, in which two molecules of ATP are hydrolysed and the rotor of the enzyme turns through two 120° steps. These findings contribute to the fundamental understanding of how the inhibitor functions and to the design of new inhibitors for the systematic analysis of the catalytic cycle of the enzyme.

xFITOM: a generic GUI tool to search for transcription factor binding sites

Locating transcription factor binding sites in genomic sequences is a key step in deciphering transcription networks. Currently available software for site search is mostly server-based, limiting the range and flexibility of this type of analysis. xFITOM is a fully customizable program for locating binding sites in genomic sequences written in C++. Through an easy-to-use interface, xFITOM that allows users an unprecedented degree of flexibility in site search. Among other features,it enables users to define motifs by mixing real sites and IUPAC consensus sequences,to search the annotated sequences of unfinished genomes and to choose among 11 different search algorithms.

AVAILABILITY

XFITOM IS AVAILABLE FOR DOWNLOAD AT: http://research.umbc.edu/˜erill.

Salt-triggered expression of the ANAC092-dependent senescence regulon in Arabidopsis thaliana

The NAC domain transcription factor ANAC092 plays a central role in leaf senescence in Arabidopsis thaliana. We recently identified 170 genes whose expression increases upon activation of ANAC092 in a chemically (estradiol) controlled experimental set-up, 78 of which are known senescence-associated genes (SAGs). In accordance with the well-known phenomenon that salt stress promotes early leaf senescence in many plant species, we previously observed salt stress-enhanced expression of many SAGs of the ANAC092 regulon. Global expression profiling now revealed that 36 genes, representing 46% of all ANAC092 downstream SAGs, are induced by long-term (4 days) salt stress in shoots of Arabidopsis, whereas short-term stress (6 hours) only slightly affects gene expression. Expression analysis also showed that 14 of the 36 genes are induced by hydrogen peroxide (H2O2) treatment. Additionally, 15 senescence-associated NAC genes (senNACs), including ANAC092, respond to H2O2 exposure. Our data support the model that salt-triggered senescence is at least partly mediated through the ANAC092 gene regulatory network. Other senNACs most likely contribute to the coordination of this process, potentially in concert with H2O2-mediated signaling.

Binding site on human immunoglobulin G for the affinity ligand HWRGWV

Affinity ligand HWRGWV has demonstrated the ability to isolate human immunoglobulin G (hIgG) from mammalian cell culture media. The ligand specifically binds hIgG through its Fc portion. This work shows that deglycosylation of hIgG has no influence on its binding to the HWRGWV ligand and the ligand does not compete with Protein A or Protein G in binding hIgG. It is suggested by the mass spectrometry (MS) data and docking simulation that HWRGWV binds to the pFc portion of hIgG and interacts with the amino acids in the loop Ser383-Asn389 (SNGQPEN) located in the C(H)3 domain. Subsequent modeling has suggested a possible three-dimensional minimized solution structure for the interaction of hIgG and the HWRGWV ligand. The results support the fact that a peptide as small as a hexamer can have specific interactions with large proteins such as hIgG.

An Inductive Logic Programming Approach to Validate Hexose Binding Biochemical Knowledge

Hexoses are simple sugars that play a key role in many cellular pathways, and in the regulation of development and disease mechanisms. Current protein-sugar computational models are based, at least partially, on prior biochemical findings and knowledge. They incorporate different parts of these findings in predictive black-box models. We investigate the empirical support for biochemical findings by comparing Inductive Logic Programming (ILP) induced rules to actual biochemical results. We mine the Protein Data Bank for a representative data set of hexose binding sites, non-hexose binding sites and surface grooves. We build an ILP model of hexose-binding sites and evaluate our results against several baseline machine learning classifiers. Our method achieves an accuracy similar to that of other black-box classifiers while providing insight into the discriminating process. In addition, it confirms wet-lab findings and reveals a previously unreported Trp-Glu amino acids dependency.

An in silico approach to map the binding site of doxorubicin on hemoglobin

Binding modalities of doxorubicin (DOX), a widely used antineoplastic anthracyline antibiotic with hemoglobin (Hb) have been studied. The protein and the ligand were prepared using CORINA and protonated with insight II. The best conformation was sought by employing GOLDV. Molecular modeling calculations showed that DOX binds Hb to a non-classical drug binding site. The alpha subunit of Hb has been assigned to posses the binding site for DOX with a binding affinity (Ka) = 16.98 x10(3) mol(-1). The interaction was found to be thermodynamically favorable (DeltaG degrees = -66.23 KJmol(-1)). The analysis of DOX binding site to Hb suggested that the types of interactions that contribute in this binding are hydrophobic contacts, hydrogen bonding and electrostatic interactions.

Diversity of aminopeptidases, derived from four lepidopteran gene duplications, and polycalins expressed in the midgut of Helicoverpa armigera: identification of proteins binding the delta-endotoxin, Cry1Ac of Bacillus thuringiensis

Helicoverpa armigera midgut proteins that bind the Bacillus thuringiensis (Bt) delta-endotoxin Cry1Ac were purified by affinity chromatography. SDS-PAGE showed that several proteins were eluted with N-acetylgalactosamine and no further proteins were detected after elution with urea. Tandem mass spectral data for tryptic peptides initially indicated that the proteins resembled aminopeptidases (APNs) from other lepidopterans and cDNA sequences for seven APNs were isolated from H. armigera through a combination of cloning with primers derived from predicted peptide sequences and established EST libraries. Phylogenetic analysis showed lepidopteran APN genes in nine clades of which five were part of a lepidopteran-specific radiation. The Cry1Ac-binding proteins were then identified with four of the seven HaAPN genes. Three of those four APNs are likely orthologs of APNs characterised as Cry1Ac-binding proteins in other lepidopterans. The fourth Cry1Ac-binding APN has orthologs not previously identified as Cry1Ac-binding partners. The HaAPN genes were expressed predominantly in the midgut through larval development. Each showed consistent expression along the length of the midgut but five of the genes were expressed at levels about two orders of magnitude greater than the remaining two. The remaining mass spectral data identified sequences encoding polycalin proteins with multiple lipocalin-like domains. A polycalin has only been previously reported in another lepidopteran, Bombyx mori, but polycalins in both species are now linked with binding of Bt Cry toxins. This is the first report of hybrid, lipocalin-like domains in shorter polycalin sequences that are not present in the longest sequence. We propose that these hybrid domains are generated by alternative splicing of the mRNA.

Cysteine-390 is the binding site of luminous substance with symplectin, a photoprotein from Okinawan squid, Symplectoteuthis oualaniensis

Symplectin is a photoprotein from a luminous squid, Symplectoteuthis oualaniensis. It has a luminous substrate, dehydrocoelenterazine (DCZ), linked through a thioether bond with a cysteine residue. We have proven the binding site of luminous substrate in symplectin by using an artificial analogue of DCZ, ortho-fluoro-DCZ (F-DCZ). F-DCZ-symplectin emitting strong blue light was reconstituted from apo-symplectin and F-DCZ. Proteolytic digestion of the reconstituted F-DCZ-symplectin afforded peptides including C(390)GLK-F-DCZ (amide), which was detected with a house assembled nano-LC-ESI-Q-TOF-MS. The chromo-peptide derived from the F-DCZ-symplectin after luminescence showed the lower molecular mass than that before the luminescence by 12 mass units, corresponding to the loss of one carbon atom upon emitting light. Thus, we have concluded that F-DCZ analogue binds to Cys390 in symplectin so as to emit light.

Cysteine-390 is the binding site of luminous substance with symplectin, a photoprotein from Okinawan squid, Symplectoteuthis oualaniensis

Symplectin is a photoprotein from a luminous squid, Symplectoteuthis oualaniensis. It has a luminous substrate, dehydrocoelenterazine (DCZ), linked through a thioether bond with a cysteine residue. We have proven the binding site of luminous substrate in symplectin by using an artificial analogue of DCZ, ortho-fluoro-DCZ (F-DCZ). F-DCZ-symplectin emitting strong blue light was reconstituted from apo-symplectin and F-DCZ. Proteolytic digestion of the reconstituted F-DCZ-symplectin afforded peptides including C(390)GLK-F-DCZ (amide), which was detected with a house assembled nano-LC-ESI-Q-TOF-MS. The chromo-peptide derived from the F-DCZ-symplectin after luminescence showed the lower molecular mass than that before the luminescence by 12 mass units, corresponding to the loss of one carbon atom upon emitting light. Thus, we have concluded that F-DCZ analogue binds to Cys390 in symplectin so as to emit light.

Being a binding site: characterizing residue composition of binding sites on proteins

The Protein Data Bank contains the description of more than 45,000 three-dimensional protein and nucleic-acid structures today. Started to exist as the computer-readable depository of crystallographic data complementing printed articles, the proper interpretation of the content of the individual files in the PDB still frequently needs the detailed information found in the citing publication. This fact implies that the fully automatic processing of the whole PDB is a very hard task. We first cleaned and re-structured the PDB data, then analyzed the residue composition of the binding sites in the whole PDB for frequency and for hidden association rules. Main results of the paper: (i) the cleaning and repairing algorithm (ii) redundancy elimination from the data (iii) application of association rule mining to the cleaned non-redundant data set. We have found numerous significant relations of the residue-composition of the ligand binding sites on protein surfaces, summarized in two figures. One of the classical data-mining methods for exploring implication-rules, the association-rule mining, is capable to find previously unknown residue-set preferences of bind ligands on protein surfaces. Since protein-ligand binding is a key step in enzymatic mechanisms and in drug discovery, these uncovered preferences in the study of more than 19,500 binding sites may help in identifying new binding protein-ligand pairs.

Three dimensional modeling of N-terminal region of galanin and its interaction with the galanin receptor

The neuropeptide galanin comes under the powerful and versatile modulators of classical neurotransmitters and is present in brain tissues, which are intimately involved in epileptogenesis. It acts as appealing targets for studying basic mechanisms of seizure initiation and arrest, and for the development of novel approaches for various neurodegenerative diseases. Galanin is widely distributed in the mammalian brain which controls various processes such as sensation of pain, learning, feeding, sexual behaviour, carcinogenesis, pathophysiology of neuroendocrine tumors and others. The function of galanin can be exploited through its interaction with three G-protein coupled receptors subtypes such as GalR1, GalR2 and GalR3. The N-terminal region of galanin comprises about highly conserved 15 amino acid residues, which act as the crucial region for agonist-receptor binding. We have constructed a theoretical structural model for the N-terminal region of galanin from Homo sapiens by homology modeling. The stereochemistry of the model was checked using PROCHECK. The functionally conserved regions were identified by surface mapping of phylogenetic information generated by online web algorithm ConSurf. The docking studies on the pharmacologically important galanin receptors with the theoretical model of N-terminal region of galanin predicted crucial residues for binding which would be useful in the development of novel leads for neurodegenerative disorders.

Biological activity and binding site characteristics of the PA1b Entomotoxin on insects from different orders

The aim of this work was to investigate both the biological activity of an entomotoxin, the pea albumin 1b (PA1b), and the presence or absence of its binding site within an array of insect species. The data obtained showed that insect sensitivity was not related to its taxonomic position. Moreover, PA1b was not toxic to several tested microorganisms. However, the binding site was found to be conserved among very different insects, displaying similar thermodynamic constants regardless of the in vivo species sensitivity. The binding site alone was, therefore, not sufficient for toxicity. One exception was the pea weevil, Bruchus pisorum, which was the only tested species without any detectable binding activity. These findings indicate that the binding site probably has an important endogenous function in insects and that adaptation to pea seeds resulted in the elimination of the toxin binding activity in two independent insect lineages. Other mechanisms are likely to interact with the toxin effects, although they are still largely unknown, but there is no evidence of any specific degradation of PA1b in the midgut of insects insensitive to the toxin, such as Drosophila melanogaster or Mamestra brassicae.

Docking studies of agonists and antagonists suggest an activation pathway of the A3 adenosine receptor

Structural determinants of ligand efficacy in the human A(3) adenosine receptor (AR) were studied using pharmacophore and docking analyses of various categories of A(3) selective ligands: inverse agonist, neutral antagonist (nonnucleoside and nucleoside), and agonist (partial and full). The homology modeling of GPCRs was adapted to provide two templates: the rhodopsin-based resting state for antagonist binding and a putative Meta I state, conformationally altered at a key residue (W6.48), for agonist binding. The preferential binding domains and/or local conformational changes associated with docking of three high affinity A(3)AR ligands were compared: inverse agonist PSB-11 1 ((R)-8-ethyl-4-methyl-2-phenyl-imidazo[2,1-i]purin-5-one); neutral antagonist MRE-3008F20 7 (5-[[(4-methoxyphenyl)amino]carbonyl]amino-8-methyl-2-(2-furyl)pyra-zolo[4,3-e]1,2,4-triazolo[1,5-c]pyrimidine), and full agonist Cl-IB-MECA 21 (2-chloro-N(6)-(3-iodobenzyl)-5'-N-methylcarboxamidoadenosine) to define a distinct recognition mode for each. Ribose-containing agonists were more hydrophilic than nonnucleoside antagonists, and H-bonding ability at the ribose 3'- and 5'-positions was required for agonism. From the receptor perspective, common requirements for activation included the destabilization of H-bond networks at W6.48 and H7.43, the specific interactions of the ribose moiety in its putative hydrophilic pocket at T3.36, S7.42, and H7.43, the stabilization of the complex by inward movement of F5.43, and the characteristic rotation of W6.48. By analogy, outward rotation of the W6.48 side-chain upon activation of an internally-crosslinking mutant M(3) muscarinic receptor was indicated by constrained molecular dynamics (MD). Our results are consistent with an anti-clockwise rotation (from the extracellular view) of transmembrane domains 3, 5, 6, and 7, as proposed for other Family A GPCRs. Thus, the putative conformational changes associated with A(3)AR activation indicate a shared mechanism of GPCR activation similar to rhodopsin.

Molecular transport through channels and pores: effects of in-channel interactions and blocking

Facilitated translocation of molecules through channels and pores is of fundamental importance for transmembrane transport in biological systems. Several such systems have specific binding sites inside the channel, but a clear understanding of how the interaction between channel and molecules affects the flow is still missing. We present a generic analytical treatment of the problem that relates molecular flow to the first passage time across and the number of particles inside the channel. Both quantities depend in different ways on the channel properties. For the idealized case of noninteracting molecules, we find an increased flow whenever there is a binding site in the channel, despite an increased first passage time. In the more realistic case that molecules may block the channel, we find an increase of flow only up to a certain threshold value of the binding strength and a dependence on the sign of the concentration gradient, i.e., asymmetric transport. The optimal binding strength in that case is analyzed. In all cases the reason for transport facilitation is an increased occupation probability of a particle inside the channel that overcomes any increase in the first passage time because of binding.

Glutamate receptor trafficking: endoplasmic reticulum quality control involves ligand binding and receptor function

The glutamate receptor (GluR) agonist-binding site consists of amino acid residues in the extracellular S1 and S2 domains in the N-terminal and M3-M4 loop regions, respectively. In the present study, we sought to confirm that the conserved ligand-binding residues identified in the AMPA receptor S1S2 domains also participate in ligand binding of GluR6 kainate receptors. Amino acid substitutions were made in the GluR6 parent at R523, T690, and E738 to alter their potential interactions with ligand. Mutant receptors were expressed in human embryonic kidney 293 cells, confirmed by Western blot analysis, and tested by [3H]kainate binding and patch-clamp recording. Each of the binding site mutations was sufficient to reduce [3H]kainate binding to undetectable levels and eliminate functional responses to glutamate or kainate. As with our studies of other nonfunctional mutants (Fleck et al., 2003), immunocytochemical staining and cell-surface biotinylation studies showed that the mutant receptors were retained intracellularly and did not traffic to the cell surface. Endoglycosidase-H digests and colocalization with endoplasmic reticulum (ER) markers demonstrated that the mutant receptors are immaturely glycosylated and retained in the ER. Immunoprecipitation, native PAGE, and functional studies confirmed that the GluR6-binding site mutants are capable of multimeric assembly, indicating their retention in the ER does not result from a gross protein folding error. Together, these results confirm the role of R523, T690, and E738 directly in ligand binding to GluR6 and further support our previous report that nonfunctional GluRs are retained intracellularly by a functional checkpoint in ER quality control.

What's new in the IGF-binding proteins?

Since their initial discovery over 25 years ago as IGF carrier proteins, the insulin-like growth factor binding protein (IGFBP) family has grown to six members, ranging in size from 216 to 289 amino acids. The assumption over the years has been that this family of proteins, having higher affinities for IGF-I and IGF-II than does the IGF-IR, serves to block access of these ligands to the receptor. Although the need for such regulatory proteins is consistent with the constitutive secretion of IGFs from many cell types, it is not surprising that additional functions have begun to be uncovered for these proteins. This review will examine new and old actions of the IGFBPs from a biochemical and cell biological perspective.

Influence of pore residues on permeation properties in the Kv2.1 potassium channel. Evidence for a selective functional interaction of K+ with the outer vestibule

The Kv2.1 potassium channel contains a lysine in the outer vestibule (position 356) that markedly reduces open channel sensitivity to changes in external [K(+)]. To investigate the mechanism underlying this effect, we examined the influence of this outer vestibule lysine on three measures of K(+) and Na(+) permeation. Permeability ratio measurements, measurements of the lowest [K(+)] required for interaction with the selectivity filter, and measurements of macroscopic K(+) and Na(+) conductance, were all consistent with the same conclusion: that the outer vestibule lysine in Kv2.1 interferes with the ability of K(+) to enter or exit the extracellular side of the selectivity filter. In contrast to its influence on K(+) permeation properties, Lys 356 appeared to be without effect on Na(+) permeation. This suggests that Lys 356 limited K(+) flux by interfering with a selective K(+) binding site. Combined with permeation studies, results from additional mutagenesis near the external entrance to the selectivity filter indicated that this site was located external to, and independent from, the selectivity filter. Protonation of a naturally occurring histidine in the same outer vestibule location in the Kv1.5 potassium channel produced similar effects on K(+) permeation properties. Together, these results indicate that a selective, functional K(+) binding site (e.g., local energy minimum) exists in the outer vestibule of voltage-gated K(+) channels. We suggest that this site is the location of K(+) hydration/dehydration postulated to exist based on the structural studies of KcsA. Finally, neutralization of position 356 enhanced outward K(+) current magnitude, but did not influence the ability of internal K(+) to enter the pore. These data indicate that in Kv2.1, exit of K(+) from the selectivity filter, rather than entry of internal K(+) into the channel, limits outward current magnitude. We discuss the implications of these findings in relation to the structural basis of channel conductance in different K(+) channels.

Mechanism of lonidamine inhibition of the CFTR chloride channel

1. The cystic fibrosis transmembrane conductance regulator (CFTR) Cl(-) channel is blocked by a broad range of organic anionic compounds. Here we investigate the effects of the indazole compound lonidamine on CFTR channels expressed in mammalian cell lines using patch clamp recording. 2. Application of lonidamine to the intracellular face of excised membrane patches caused a voltage-dependent block of CFTR currents, with an apparent K(d) of 58 micro M at -100 mV. 3. Block by lonidamine was apparently independent of channel gating but weakly sensitive to the extracellular Cl(-) concentration. 4. Intracellular lonidamine led to the introduction of brief interruptions in the single channel current at hyperpolarized voltages, leading to a reduction in channel mean open time. Lonidamine also introduced a new component of macroscopic current variance. Spectral analysis of this variance suggested a blocker on rate of 1.79 micro M(-1) s(-1) and an off-rate of 143 s(-1). 5. Several point mutations within the sixth transmembrane region of CFTR (R334C, F337S, T338A and S341A) significantly weakened block of macroscopic CFTR current, suggesting that lonidamine enters deeply into the channel pore from its intracellular end. 6. These results identify and characterize lonidamine as a novel CFTR open channel blocker and provide important information concerning its molecular mechanism of action.

Efficient docking of peptides to proteins without prior knowledge of the binding site

Reliability in docking of ligand molecules to proteins or other targets is an important challenge for molecular modeling. Applications of the docking technique include not only prediction of the binding mode of novel drugs, but also other problems like the study of protein-protein interactions. Here we present a study on the reliability of the results obtained with the popular AutoDock program. We have performed systematical studies to test the ability of AutoDock to reproduce eight different protein/ligand complexes for which the structure was known, without prior knowledge of the binding site. More specifically, we look at factors influencing the accuracy of the final structure, such as the number of torsional degrees of freedom in the ligand. We conclude that the Autodock program package is able to select the correct complexes based on the energy without prior knowledge of the binding site. We named this application blind docking, as the docking algorithm is not able to "see" the binding site but can still find it. The success of blind docking represents an important finding in the era of structural genomics.