Refining near-native protein-protein docking decoys by local resampling and energy minimization

How to refine a near-native structure to make it closer to its native conformation is an unsolved problem in protein-structure and protein-protein complex-structure prediction. In this article, we first test several scoring functions for selecting locally resampled near-native protein-protein docking conformations and then propose a computationally efficient protocol for structure refinement via local resampling and energy minimization. The proposed method employs a statistical energy function based on a Distance-scaled Ideal-gas REference state (DFIRE) as an initial filter and an empirical energy function EMPIRE (EMpirical Protein-InteRaction Energy) for optimization and re-ranking. Significant improvement of final top-1 ranked structures over initial near-native structures is observed in the ZDOCK 2.3 decoy set for Benchmark 1.0 (74% whose global rmsd reduced by 0.5 angstrom or more and only 7% increased by 0.5 angstrom or more). Less significant improvement is observed for Benchmark 2.0 (38% versus 33%). Possible reasons are discussed.

SmC2P1, a C2 domain protein from Scophthalmus maximus that binds bacterial pathogen-infected lymphocytes and reduces bacterial survival

C2 domains are protein structural modules found in many eukaryotic proteins involved in signal transduction, membrane trafficking, and immune defense. Most of the studied C2 domain-containing proteins are multi-domained in structure, in which the C2 domain is an independently folded motif and plays an essential role in calcium-dependent membrane-targeting. Although C2 domains isolated from intact proteins have been studied for biological functions, no study on natural proteins containing C2 domain only has been documented. In this study, we identified a Scophthalmus maximus protein SmC2P1 that is comprised of a single C2 domain and lacks any other apparent domain structures. The deduced amino acid sequence of SmC2P1 contains 129 residues and shares 36-38% identities with the C2 domains of the perforins of several fish species. Like typical C2 domains, SmC2P1 is predicted to organize into eight beta-strands with a Ca2+-binding site located in inter-strand loops. SmC2P1 expression was detected, in deceasing order, in liver, spleen, blood, brain, muscle, kidney, gill, and heart. Experimental challenge of turbot with a bacterial pathogen significantly upregulated SmC2P1 expression in kidney in a time-dependent manner. Recombinant SmC2P1 purified from yeast exhibits no hemolytic activity but binds to pathogen-infected kidney lymphocytes in the presence of calcium. Furthermore, interaction of recombinant SmC2P1 with bacterium-infected lymphocytes reduced bacterial survival. These results indicate that SmC2P1 is a functional protein that is involved in host immune defense against bacterial infection. (C) 2010 Elsevier Ltd. All rights reserved.

Integration of relational and hierarchical network information for protein function prediction

BACKGROUND:In the current climate of high-throughput computational biology, the inference of a protein's function from related measurements, such as protein-protein interaction relations, has become a canonical task. Most existing technologies pursue this task as a classification problem, on a term-by-term basis, for each term in a database, such as the Gene Ontology (GO) database, a popular rigorous vocabulary for biological functions. However, ontology structures are essentially hierarchies, with certain top to bottom annotation rules which protein function predictions should in principle follow. Currently, the most common approach to imposing these hierarchical constraints on network-based classifiers is through the use of transitive closure to predictions.RESULTS:We propose a probabilistic framework to integrate information in relational data, in the form of a protein-protein interaction network, and a hierarchically structured database of terms, in the form of the GO database, for the purpose of protein function prediction. At the heart of our framework is a factorization of local neighborhood information in the protein-protein interaction network across successive ancestral terms in the GO hierarchy. We introduce a classifier within this framework, with computationally efficient implementation, that produces GO-term predictions that naturally obey a hierarchical 'true-path' consistency from root to leaves, without the need for further post-processing.CONCLUSION:A cross-validation study, using data from the yeast Saccharomyces cerevisiae, shows our method offers substantial improvements over both standard 'guilt-by-association' (i.e., Nearest-Neighbor) and more refined Markov random field methods, whether in their original form or when post-processed to artificially impose 'true-path' consistency. Further analysis of the results indicates that these improvements are associated with increased predictive capabilities (i.e., increased positive predictive value), and that this increase is consistent uniformly with GO-term depth. Additional in silico validation on a collection of new annotations recently added to GO confirms the advantages suggested by the cross-validation study. Taken as a whole, our results show that a hierarchical approach to network-based protein function prediction, that exploits the ontological structure of protein annotation databases in a principled manner, can offer substantial advantages over the successive application of 'flat' network-based methods.

Structures of P-glycoprotein reveal its conformational flexibility and an epitope on the nucleotide-binding domain

P-glycoprotein (P-gp) is one of the best-known mediators of drug efflux-based multidrug resistance in many cancers. This validated therapeutic target is a prototypic, plasma membrane resident ATPBinding Cassette transporter that pumps xenobiotic compounds out of cells. The large, polyspecific drug-binding pocket of P-gp recognizes a variety of structurally unrelated compounds. The transport of these drugs across the membrane is coincident with changes in the size and shape of this pocket during the course of the transport cycle. Here, we present the crystal structures of three inward-facing conformations of mouse P-gp derived from two different crystal forms. One structure has a nanobody bound to the C-terminal side of the first nucleotide-binding domain. This nanobody strongly inhibits the ATP hydrolysis activity of mouse Pgp by hindering the formation of a dimeric complex between the ATP-binding domains, which is essential for nucleotide hydrolysis. Together, these inward-facing conformational snapshots of P-gp demonstrate a range of flexibility exhibited by this transporter, which is likely an essential feature for the binding and transport of large, diverse substrates. The nanobody-bound structure also reveals a unique epitope on P-gp.

A structural basis for different antifreeze protein roles

Antifreeze proteins (AFPs) are produced by a variety of organisms to either protect them from freezing or help them tolerate being frozen. Recent structural work has shown that AFPs bind to ice using ordered surface waters on a particular surface of the protein called the ice-binding site (IBS). These 'anchored clathrate' waters fuse to particular planes of an ice crystal and hence irreversibly bind the AFP to its ligand. An AFP isolated from the perennial ryegrass, Lolium perenne (LpAFP) was previously modelled as a right-handed beta helix with two proposed IBSs. Steric mutagenesis, where small side chains were replaced with larger ones, determined that only one of the putative IBSs was responsible for binding ice. The mutagenesis work also partly validated the fold of the computer-generated model of this AFP. In order to determine the structure of the protein, LpAFP was crystallized and solved to 1.4 Å resolution. The protein folds as an untwisted left-handed beta-helix, of opposite handedness to the model. The IBS identified by mutagenesis is remarkably flat, but less regular than the IBS of most other AFPs. Furthermore, several of the residues constituting the IBS are in multiple conformations. This irregularity may explain why LpAFP causes less thermal hysteresis than many other AFPs. Its imperfect IBS is also argued to be responsible for LpAFP's heightened ice-recrystallization inhibition activity. The structure of LpAFP is the first for a plant AFP and for a protein responsible for providing freeze tolerance rather than freeze resistance. To help understand what constitutes an IBS, a non-ice-binding homologue of type III AFP, sialic acid synthase (SAS), was engineered for ice binding. Point mutations were made to the germinal IBS of SAS to mimic key features seen in type III AFP. The crystal structures of some of the mutant proteins showed that the potential IBS became less charged and flatter as the mutations progressed, and ice affinity was gained. This proof-of-principle study highlights some of the difficulties in AFP engineering.

Population-based local search for protein folding simulation in the MJ energy model and cubic lattices

We present experimental results on benchmark problems in 3D cubic lattice structures with the Miyazawa-Jernigan energy function for two local search procedures that utilise the pull-move set: (i) population-based local search (PLS) that traverses the energy landscape with greedy steps towards (potential) local minima followed by upward steps up to a certain level of the objective function; (ii) simulated annealing with a logarithmic cooling schedule (LSA). The parameter settings for PLS are derived from short LSA-runs executed in pre-processing and the procedure utilises tabu lists generated for each member of the population. In terms of the total number of energy function evaluations both methods perform equally well, however. PLS has the potential of being parallelised with an expected speed-up in the region of the population size. Furthermore, both methods require a significant smaller number of function evaluations when compared to Monte Carlo simulations with kink-jump moves. (C) 2009 Elsevier Ltd. All rights reserved.

Antidepressant specificity of serotonin transporter suggested by three LeuT-SSRI structures.

Sertraline and fluoxetine are selective serotonin re-uptake inhibitors (SSRIs) that are widely prescribed to treat depression. They exert their effects by inhibiting the presynaptic plasma membrane serotonin transporter (SERT). All SSRIs possess halogen atoms at specific positions, which are key determinants for the drugs' specificity for SERT. For the SERT protein, however, the structural basis of its specificity for SSRIs is poorly understood. Here we report the crystal structures of LeuT, a bacterial SERT homolog, in complex with sertraline, R-fluoxetine or S-fluoxetine. The SSRI halogens all bind to exactly the same pocket within LeuT. Mutation at this halogen-binding pocket (HBP) in SERT markedly reduces the transporter's affinity for SSRIs but not for tricyclic antidepressants. Conversely, when the only nonconserved HBP residue in both norepinephrine and dopamine transporters is mutated into that found in SERT, their affinities for all the three SSRIs increase uniformly. Thus, the specificity of SERT for SSRIs is dependent largely on interaction of the drug halogens with the protein's HBP.

Rhodopsin and the Others: A Historical Perspective on Structural Studies of G Protein-Coupled Receptors

The role of rhodopsin as a structural prototype for the study of the whole superfamily of G protein-coupled receptors (GPCRs) is reviewed in an historical perspective. Discovered at the end of the nineteenth century, fully sequenced since the early 1980s, and with direct three-dimensional information available since the 1990s, rhodopsin has served as a platform to gather indirect information on the structure of the other superfamily members. Recent breakthroughs have elicited the solution of the structures of additional receptors, namely the beta 1- and beta 2-adrenergic receptors and the A(2A) adenosine receptor, now providing an opportunity to gauge the accuracy of homology modeling and molecular docking techniques and to perfect the computational protocol. Notably, in coordination with the solution of the structure of the A(2A) adenosine receptor, the first "critical assessment of GPCR structural modeling and docking" has been organized, the results of which highlighted that the construction of accurate models, although challenging, is certainly achievable. The docking of the ligands and the scoring of the poses clearly emerged as the most difficult components. A further goal in the field is certainly to derive the structure of receptors in their signaling state, possibly in complex with agonists. These advances, coupled with the introduction of more sophisticated modeling algorithms and the increase in computer power, raise the expectation for a substantial boost of the robustness and accuracy of computer-aided drug discovery techniques in the coming years.

Mechanism of Activation of a G Protein-coupled Receptor, the Human Cholecystokinin-2 Receptor

G protein-coupled receptors (GPCRs) represent a major focus in functional genomics programs and drug development research, but their important potential as drug targets contrasts with the still limited data available concerning their activation mechanism. Here, we investigated the activation mechanism of the cholecystokinin-2 receptor (CCK2R). The three-dimensional structure of inactive CCK2R was homology-modeled on the basis of crystal coordinates of inactive rhodopsin. Starting from the inactive CCK2R modeled structure, active CCK2R (namely cholecystokinin-occupied CCK2R) was modeled by means of steered molecular dynamics in a lipid bilayer and by using available data from other GPCRs, including rhodopsin. By comparing the modeled structures of the inactive and active CCK2R, we identified changes in the relative position of helices and networks of interacting residues, which were expected to stabilize either the active or inactive states of CCK2R. Using targeted molecular dynamics simulations capable of converting CCK2R from the inactive to the active state, we delineated structural changes at the atomic level. The activation mechanism involved significant movements of helices VI and V, a slight movement of helices IV and VII, and changes in the position of critical residues within or near the binding site. The mutation of key amino acids yielded inactive or constitutively active CCK2R mutants, supporting this proposed mechanism. Such progress in the refinement of the CCK2R binding site structure and in knowledge of CCK2R activation mechanisms will enable target-based optimization of nonpeptide ligands.

The respiratory syncytial virus small hydrophobic protein is phosphorylated via a mitogen-activated protein kinase p38-dependent tyrosine kinase activity during virus infection

The phosphorylation status of the small hydrophobic (SH) protein of respiratory syncytial virus (RSV) was examined in virus-infected Vero cells. The SH protein v.,as isolated from [S-35]methionine- and [P-33]orthophosphate-labelled IRSV-infected cells and analysed by SDS-PAGE. In each case, a protein product of the expected size for the SH protein was observed. Phosphoamino acid analysis and reactivity with the phosphotyrosine specific antibody PY20 showed that the SH protein was modified by tyrosine phosphorylation. The role or tyrosine kinase activity in SH protein phosphorylation was confirmed by the use of genistein, a broad-spectrum tyrosine kinase inhibitor, to inhibit SH protein phosphorylation. Further analysis showed that the different glycosylated forms of the SH protein were phosphorylated, as was the oligomeric form of the protein. Phosphorylation of the SH protein was specifically inhibited by the mitogen-activated protein kinase (MAPK) p38 inhibitor SB203580, suggesting that SH protein phosphorylation occurs via a MAPK p38-dependent pathway. Analysis of virus-infected cells using fluorescence microscopy showed that, although the SH protein was distributed throughout the cytoplasm, it appeared to accumulate, at low levels, in the endoplasmic reticulum/Golgi complex, confirming recent observations. However, in the presence of SB203580. an increased accumulation of the SH protein in the Golgi complex was observed, although other virus structures, such as virus filaments and inclusion bodies, remained largely unaffected. These results showed that during RSV infection, the SH protein is modified by an MAPK p38-dependant tyrosine kinase activity and that this modification influences its cellular distribution.

Fasciola hepatica: Histological changes in the reproductive structures of triclabendazole (TCBZ)-sensitive and TCBZ-resistant flukes after treatment in vivo with TCBZ and the related benzimidazole derivative, Compound Alpha

Twenty-four shed-reared lambs were each infected orally with 250 metacercariae of Fasciola hepatica, using either the triclabendazole (TCBZ)-sensitive Cullompton isolate or the TCBZ-resistant Sligo isolate. Twelve weeks after infection the lambs were treated with TCBZ (10 mg/kg) or with the experimental fasciolicide, Compound Alpha (Cpd alpha), a benzimidazole derivative of TCBZ (15 mg/kg). The lambs were euthanised 48,72 and 96 h after TCBZ treatment, or 24, 48 and 72 h after Cpd a treatment, and flukes were collected from the liver and/or gall bladder of each animal. Untreated animals harbouring 12-week infections were euthanised 24 h after administration of anthelmintic to the treatment groups, and the untreated flukes provided control material. A semi-quantitative assessment of the degree of histological change induced by the two drugs after different times of exposure was achieved by scoring the intensity of three well-defined lesions that developed in the testes and uteri of a representative sample of flukes from each lamb. In general, it was found that in those tissues where active meiosis and/or mitosis occurred (testis, ovary, and vitelline follicles), there was progressive loss of cell content due to apparent failure of cell division to keep pace with expulsion of the mature or effete products. Further, actively dividing cell types tended to become individualised, rounded and condensed, characteristic of apoptotic cell death. Protein synthetic activity was apparently inhibited in the Mehlis' secretory cells. In the uterus, where successful formation of shelled eggs represents the culmination of a complex sequence of cytokinetic, cytological and synthetic activity involving the vitelline follicles, the ovary and the Mehlis' gland, histological evidence indicating failure of ovigenesis was evident from 24 h post-treatment onwards. The development of these lesions may be related to the known antitubulin activity of the benzimidazole class of anthelmintics, to the induction of apoptosis in cells where mitosis or meiosis has aborted due to failure of spindle formation, and to drug-induced inhibition of protein synthesis. The semi-quantitative findings indicated that Cpd a is slightly less efficacious than TCBZ itself in causing histological damage to the reproductive structures of TCBZ-sensitive flukes, and that, like TCBZ, it caused no histological damage in flukes of the TCBZ-resistant isolate. This study illustrates the potential utility of histological techniques for conveniently screening representative samples of flukes in field trials designed to validate instances of drug resistance or to test the efficacy of new products against known drug-resistant and drug-susceptible fluke isolates. It also provides reference criteria for drug-induced histopathological changes in fluke reproductive structures which may aid interpretation of TEM findings. (C) 2009 Elsevier B.V. All rights reserved.

Immunohistochemical localization of a Shaker-related voltage-gated potassium channel protein in Schistosoma mansoni (Trematoda: Digenea)

We have recently isolated a cDNA (SKV1.1) encoding a Shakei-related K+ channel from the human parasitic trematode Schistosoma mansoni. In order to better understand the functions of SKv1.1 protein, the distribution of SKv1.1 protein in adult S. mansoni was analyzed by immunohistochemistry using a region-specific antibody. SKV1.1 proteins were widely expressed in the nervous and muscular systems. The strongest immunoreactivity (IR) was observed in the nervous system of both male and female. In the nervous system, IR for SKv1.1 proteins was localized in cell bodies and nerve fibers of the anterior ganglia, the central commissure, and the main nerve cords. IR was also observed in the dorsal and the ventral peripheral nerve nets, fine nerve fibers entering into a variety of structures such as the dorsal tubercles, longitudinal and ventral muscle fibers, and oral and ventral suckers. In the muscular system, SKv1.1 proteins were localized to the longitudinal, circular, and ventral muscle fibers of male as well as in isolated muscle fibers where native A-type K+ currents were measured. Moderate IR was also seen in a large number of cell bodies in the parenchyma. These results indicate that SKv1.1 protein may play an important role in the regulation of the excitability of neurons and muscle cells of S. mansoni. (C) 1995 Academic Press, Inc.

Functional analysis of the Campylobacter jejuni N-linked protein glycosylation pathway

We describe in this report the characterization of the recently discovered N-linked glycosylation locus of the human bacterial pathogen Campylobacter jejuni, the first such system found in a species from the domain Bacteria. We exploited the ability of this locus to function in Escherichia coli to demonstrate through mutational and structural analyses that variant glycan structures can be transferred onto protein indicating the relaxed specificity of the putative oligosaccharyltransferase PglB. Structural data derived from these variant glycans allowed us to infer the role of five individual glycosyltransferases in the biosynthesis of the N-linked heptasaccharide. Furthermore, we show that C. jejuni- and E. coli-derived pathways can interact in the biosynthesis of N-linked glycoproteins. In particular, the E. coli encoded WecA protein, a UDP-GlcNAc: undecaprenylphosphate GlcNAc-1-phosphate transferase involved in glycolipid biosynthesis, provides for an alternative N-linked heptasaccharide biosynthetic pathway bypassing the requirement for the C. jejuni-derived glycosyltransferase PglC. This is the first experimental evidence that biosynthesis of the N-linked glycan occurs on a lipid-linked precursor prior to transfer onto protein. These findings provide a framework for understanding the process of N-linked protein glycosylation in Bacteria and for devising strategies to exploit this system for glycoengineering.

Unraveling the Structure and Function of G Protein-Coupled Receptors Through NMR Spectroscopy

G protein-coupled receptors (GPCRs) are a large superfamily of signaling proteins expressed on the plasma membrane. They are involved in a wide range of physiological processes and, therefore, are exploited as drug targets in a multitude of therapeutic areas. In this extent, knowledge of structural and functional properties of GPCRs may greatly facilitate rational design of modulator compounds. Solution and solid-state nuclear magnetic resonance (NMR) spectroscopy represents a powerful method to gather atomistic insights into protein structure and dynamics. In spite of the difficulties inherent the solution of the structure of membrane proteins through NMR, these methods have been successfully applied, sometimes in combination with molecular modeling, to the determination of the structure of GPCR fragments, the mapping of receptor-ligand interactions, and the study of the conformational changes associated with the activation of the receptors. In this review, we provide a summary of the NMR contributions to the study of the structure and function of GPCRs, also in light of the published crystal structures.

Organization of mammary epithelial cells into 3D acinar structures requires glucocorticoid and JNK signaling

Mammary epithelial cells cultured on a concentrated laminin-rich extracellular matrix formed 3D acinar structures that matured to polarized monolayers surrounding a lumen. In the absence of glucocorticoids mature acinus formation failed and the expression of an acinus-associated, activator protein 1 (AP1) and nuclear factor kappaB transcription factor DNA-binding profile was lost. Treatment with the JNK inhibitor, SP600125, caused similar effects, whereas normal organization of the mammary epithelial cells as acini caused JNK activation in a glucocorticoid-dependent manner. The forming acini expressed BRCA1, GADD45beta, MEKK4, and the JNK activating complex GADD 45beta-MEKK4 in a glucocorticoid-dependent fashion. JNK catalyzed phosphorylation of c-Jun was also detected in the acini. In addition, expression of beta4 integrin and in situ occupation of its promoter by AP1 components, c-Jun and Fos, was glucocorticoid dependent. These results suggest that glucocortocoid signaling regulates acinar integrity through a pathway involving JNK regulation of AP1 transcription factors and beta4 integrin expression.