Counterbalance of ligand‐ and self‐coupled motions characterizes multispecificity of ubiquitin

Date hub proteins are a type of proteins that show multispecificity in a time‐dependent manner. To understand dynamic aspects of such multispecificity we studied Ubiquitin as a typical example of a date hub protein. Here we analyzed 9 biologically relevant Ubiquitin‐protein (ligand) heterodimer structures by using normal mode analysis based on an elastic network model. Our result showed that the self‐coupled motion of Ubiquitin in the complex, rather than its ligand‐coupled motion, is similar to the motion of Ubiquitin in the unbound condition. The ligand‐coupled motions are correlated to the conformational change between the unbound and bound conditions of Ubiquitin. Moreover, ligand‐coupled motions favor the formation of the bound states, due to its in‐phase movements of the contacting atoms at the interface. The self‐coupled motions at the interface indicated loss of conformational entropy due to binding. Therefore, such motions disfavor the formation of the bound state. We observed that the ligand‐coupled motions are embedded in the motions of unbound Ubiquitin. In conclusion, multispecificity of Ubiquitin can be characterized by an intricate balance of the ligand‐ and self‐coupled motions, both of which are embedded in the motions of the unbound form.     

Structural and functional characterization of a ubiquitin variant engineered for tight and specific binding to an alpha‐helical ubiquitin interacting motif

Ubiquitin interacting motifs (UIMs) are short α‐helices found in a number of eukaryotic proteins. UIMs interact weakly but specifically with ubiquitin conjugated to other proteins, and in so doing, mediate specific cellular signals. Here we used phage display to generate ubiquitin variants (UbVs) targeting the N‐terminal UIM of the yeast Vps27 protein. Selections yielded UbV.v27.1, which recognized the cognate UIM with high specificity relative to other yeast UIMs and bound with an affinity more than two orders of magnitude higher than that of ubiquitin. Structural and mutational studies of the UbV.v27.1‐UIM complex revealed the molecular details for the enhanced affinity and specificity of UbV.v27.1, and underscored the importance of changes at the binding interface as well as at positions that do not contact the UIM. Our study highlights the power of the phage display approach for selecting UbVs with unprecedented affinity and high selectivity for particular α‐helical UIM domains within proteomes, and it establishes a general approach for the development of inhibitors targeting interactions of this type.     

Exhaustive comparison and classification of ligand‐binding surfaces in proteins

Many proteins function by interacting with other small molecules (ligands). Identification of ligand‐binding sites (LBS) in proteins can therefore help to infer their molecular functions. A comprehensive comparison among local structures of LBSs was previously performed, in order to understand their relationships and to classify their structural motifs. However, similar exhaustive comparison among local surfaces of LBSs (patches) has never been performed, due to computational complexity. To enhance our understanding of LBSs, it is worth performing such comparisons among patches and classifying them based on similarities of their surface configurations and electrostatic potentials. In this study, we first developed a rapid method to compare two patches. We then clustered patches corresponding to the same PDB chemical component identifier for a ligand, and selected a representative patch from each cluster. We subsequently exhaustively as compared the representative patches and clustered them using similarity score, PatSim. Finally, the resultant PatSim scores were compared with similarities of atomic structures of the LBSs and those of the ligand‐binding protein sequences and functions. Consequently, we classified the patches into ～2000 well‐characterized clusters. We found that about 63% of these clusters are used in identical protein folds, although about 25% of the clusters are conserved in distantly related proteins and even in proteins with cross‐fold similarity. Furthermore, we showed that patches with higher PatSim score have potential to be involved in similar biological processes.     

Imidazo[1,2‐c]pyrimidin‐5(6H)‐one as a novel core of cyclin‐dependent kinase 2 inhibitors: Synthesis,activity measurement,docking, and quantum mechanical scoring

We report on the synthesis, activity testing, docking, and quantum mechanical scoring of novel imidazo[1,2‐c]pyrimidin‐5(6H)‐one scaffold for cyclin‐dependent kinase 2 (CDK2) inhibition. A series of 26 compounds substituted with aromatic moieties at position 8 has been tested in in vitro enzyme assays and shown to inhibit CDK2. 2D structure‐activity relationships have ascertained that small substituents at position 8 (up to the size of naphtyl or methoxyphenyl) generally lead to single‐digit micromolar IC₅₀ values, whereas bigger substituents (substituted biphenyls) decreased the compounds' activities. The binding modes of the compounds obtained using Glide docking have exhibited up to 2 hinge‐region hydrogen bonds to CDK2 and differed in the orientation of the inhibitor core and the placement of the 8‐substituents. Semiempirical quantum mechanics‐based scoring identified probable favourable binding modes, which will serve for future structure‐based design and synthetic optimization of substituents of the heterocyclic core. In summary, we have identified a novel core for CDK2 inhibition and will explore it further to increase the potencies of the compounds and also monitor selectivities against other protein kinases.     

Ligand binding and homology modelling of insect odorant‐binding proteins

This review describes the main characteristics of odorant‐binding proteins (OBPs) for homology modelling and presents a summary of structure prediction studies on insect OBPs, along with the steps involved and some limitations and improvements. The technique involves a computing approach to model protein structures and is based on a comparison between a target (unknown structure) and one or more templates (experimentally determined structures). As targets for structure prediction, OBPs are considered to play a functional role for recognition, desorption, scavenging, protection and transportation of hydrophobic molecules (odourants) across an aqueous environment (lymph) to olfactory receptor neurones (ORNs) located in sensilla, the main olfactory units of insect antennae. Lepidopteran pheromone‐binding proteins, a subgroup of OBPs, are characterized by remarkable structural features, in which high sequence identities (approximately 30%) among these OBPs and a large number of available templates can facilitate the prediction of precise homology models. Approximately 30 studies have been performed on insect OBPs using homology modelling as a tool to predict their structures. Although some of the studies have assessed ligand‐binding affinity using structural information and biochemical measurements, few have performed docking and molecular dynamic (MD) simulations as a virtual method to predict best ligands. Docking and MD simulations are discussed in the context of discovery of novel semiochemicals (super‐ligands) using homology modelling to conceive further strategies in insect management.     

Exploring the binding diversity of intrinsically disordered proteins involved in one‐to‐many binding

Molecular recognition features (MoRFs) are intrinsically disordered protein regions that bind to partners via disorder‐to‐order transitions. In one‐to‐many binding, a single MoRF binds to two or more different partners individually. MoRF‐based one‐to‐many protein–protein interaction (PPI) examples were collected from the Protein Data Bank, yielding 23 MoRFs bound to 2–9 partners, with all pairs of same‐MoRF partners having less than 25% sequence identity. Of these, 8 MoRFs were bound to 2–9 partners having completely different folds, whereas 15 MoRFs were bound to 2–5 partners having the same folds but with low sequence identities. For both types of partner variation, backbone and side chain torsion angle rotations were used to bring about the conformational changes needed to enable close fits between a single MoRF and distinct partners. Alternative splicing events (ASEs) and posttranslational modifications (PTMs) were also found to contribute to distinct partner binding. Because ASEs and PTMs both commonly occur in disordered regions, and because both ASEs and PTMs are often tissue‐specific, these data suggest that MoRFs, ASEs, and PTMs may collaborate to alter PPI networks in different cell types. These data enlarge the set of carefully studied MoRFs that use inherent flexibility and that also use ASE‐based and/or PTM‐based surface modifications to enable the same disordered segment to selectively associate with two or more partners. The small number of residues involved in MoRFs and in their modifications by ASEs or PTMs may simplify the evolvability of signaling network diversity.     

Key sites residues involved in interacting with chemicals of pheromone‐binding proteins from Lymantria dispar

Pheromone‐binding proteins (PBPs) are distributed widely on the antennae of insects, and they are believed to be involved in the process of chemical signal transduction, but their interaction with chemicals is largely unknown. Here, we present our findings on the key amino acid residues of PBPs in the gypsy moth, Lymantria dispar. Potential key residues were screened with the Calculate Mutation Energy program and molecular docking methods. Mutated proteins were obtained by mutating residues to alanine via site‐directed mutagenesis. Circular dichroism (CD) spectroscopy showed that the mutated proteins formed α‐helix, and the stability of protein structure was influenced due to mutations. Fluorescence binding assays were further conducted with the mutated proteins, sex pheromones and analogues. Results showed that to PBP 1, tyrosine at position 41 and phenylalanine at position 76 could be the key amino acid residues influencing the stability of structure; in addition, phenylalanine at 36 and lysine at position 94 could be key amino acid residues interacting with chemicals. To PBP 2, glycine at position 49, phenylalanine at position 76 and lysine at position 121 could be the key amino acid residues in the structural stability. These results shed light on the relationship between the specific amino acids and functions of PBPs in transmitting the chemical signals.     

Dynamic features of homodimer interfaces calculated by normal‐mode analysis

Knowledge of the dynamic features of protein interfaces is necessary for a deeper understanding of protein–protein interactions. We performed normal‐mode analysis (NMA) of 517 nonredundant homodimers and their protomers to characterize dimer interfaces from a dynamic perspective. The motion vector calculated by NMA for each atom of a dimer was decomposed into internal and external motion vectors in individual component subunits, followed by the averaging of time‐averaged correlations between these vectors over atom pairs in the interface. This averaged correlation coefficient (ACC) was defined for various combinations of vectors and investigated in detail. ACCs decrease exponentially with an increasing interface area and r‐value, that is, interface area divided by the entire subunit surface area. As the r‐value reflects the nature of dimer formation, the result suggests that both the interface area and the nature of dimer formation are responsible for the dynamic properties of dimer interfaces. For interfaces with small or medium r‐values and without intersubunit entanglements, ACCs are found to increase on dimer formation when compared with those in the protomer state. In contrast, ACCs do not increase on dimer formation for interfaces with large r‐values and intersubunit entanglements such as in interwinding dimers. Furthermore, relationships between ACCs for intrasubunit atom pairs and for intersubunit atom pairs are found to significantly differ between interwinding and noninterwinding dimers for external motions. External motions are considered as an important factor for characterizing dimer interfaces.     

Ionic strength‐dependent conformations of a ubiquitin‐like small archaeal modifier protein (SAMP1) from Haloferax volcanii

Eukaryotic ubiquitin and ubiquitin‐like systems play crucial roles in various cellular biological processes. In this work, we determined the solution structure of SAMP1 from Haloferax volcanii by NMR spectroscopy. Under low ionic conditions, SAMP1 presented two distinct conformations, one folded β‐grasp and the other disordered. Interestingly, SAMP1 underwent a conformational conversion from disorder to order with ion concentration increasing, indicating that the ordered conformation is the functional form of SAMP1 under the physiological condition of H. volcanii. Furthermore, SAMP1 could interact with proteasome‐activating nucleotidase B, supposing a potential role of SAMP1 in the protein degradation pathway mediated by proteasome.     

Structure‐based in silico identification of ubiquitin‐binding domains provides insights into the ALIX‐V:ubiquitin complex and retrovirus budding

The ubiquitylation signal promotes trafficking of endogenous and retroviral transmembrane proteins. The signal is decoded by a large set of ubiquitin (Ub) receptors that tether Ub‐binding domains (UBDs) to the trafficking machinery. We developed a structure‐based procedure to scan the protein data bank for hidden UBDs. The screen retrieved many of the known UBDs. Intriguingly, new potential UBDs were identified, including the ALIX‐V domain. Pull‐down, cross‐linking and E3‐independent ubiquitylation assays biochemically corroborated the in silico findings. Guided by the output model, we designed mutations at the postulated ALIX‐V:Ub interface. Biophysical affinity measurements using microscale‐thermophoresis of wild‐type and mutant proteins revealed some of the interacting residues of the complex. ALIX‐V binds mono‐Ub with a K_d of 119 μM. We show that ALIX‐V oligomerizes with a Hill coefficient of 5.4 and IC₅₀ of 27.6 μM and that mono‐Ub induces ALIX‐V oligomerization. Moreover, we show that ALIX‐V preferentially binds K63 di‐Ub compared with mono‐Ub and K48 di‐Ub. Finally, an in vivo functionality assay demonstrates the significance of ALIX‐V:Ub interaction in equine infectious anaemia virus budding. These results not only validate the new procedure, but also demonstrate that ALIX‐V directly interacts with Ub in vivo and that this interaction can influence retroviral budding.     

Comparison analysis of primary ligand‐binding sites in seven‐helix membrane proteins

Seven‐helix transmembrane proteins, including the G‐protein‐coupled receptors (GPCRs), mediate a broad range of fundamental cellular activities through binding to a wide range of ligands. Understanding the structural basis for the ligand‐binding selectivity of these proteins is of significance to their structure‐based drug design. Comparison analysis of proteins' ligand‐binding sites provides a useful way to study their structure‐activity relationships. Various computational methods have been developed for the binding‐site comparison of soluble proteins. In this work, we applied this approach to the analysis of the primary ligand‐binding sites of 92 seven‐helix transmembrane proteins. Results of the studies confirmed that the binding site of bacterial rhodopsins is indeed different from all GPCRs. In the latter group, further comparison of the binding sites indicated a group of residues that could be responsible for ligand‐binding selectivity and important for structure‐based drug design. Furthermore, unexpected binding‐site dissimilarities were observed among adrenergic and adenosine receptors, suggesting that the percentage of the overall sequence identity between a target protein and a template protein alone is not sufficient for selecting the best template for homology modeling of seven‐helix membrane proteins. These results provided novel insight into the structural basis of ligand‐binding selectivity of seven‐helix membrane proteins and are of practical use to the computational modeling of these proteins. © 2010 Wiley Periodicals, Inc. Biopolymers 95: 31–38, 2011.     

Iterative cluster-NMA: A tool for generating conformational transitions in proteins

Computational models provide insight into the structure-function relationship in proteins. These approaches, especially those based on normal mode analysis, can identify the accessible motion space around a given equilibrium structure. The large magnitude, collective motions identified by these methods are often well aligned with the general direction of the expected conformational transitions. However, these motions cannot realistically be extrapolated beyond the local neighborhood of the starting conformation. In this article, the iterative cluster-NMA (icNMA) method is presented for traversing the energy landscape from a starting conformation to a desired goal conformation. This is accomplished by allowing the evolving geometry of the intermediate structures to define the local accessible motion space, and thus produce an appropriate displacement. Following the derivation of the icNMA method, a set of sample simulations are performed to probe the robustness of the model. A detailed analysis of beta1,4-galactosyltransferase-T1 is also given, to highlight many of the capabilities of icNMA. Remarkably, during the transition, a helix is seen to be extended by an additional turn, emphasizing a new unknown role for secondary structures to absorb slack during transitions. The transition pathway for adenylate kinase, which has been frequently studied in the literature, is also discussed.     

General RNA‐binding proteins have a function in poly(A)‐binding protein‐dependent translation

The interaction between the poly(A)‐binding protein (PABP) and eukaryotic translational initiation factor 4G (eIF4G), which brings about circularization of the mRNA, stimulates translation. General RNA‐binding proteins affect translation, but their role in mRNA circularization has not been studied before. Here, we demonstrate that the major mRNA ribonucleoprotein YB‐1 has a pivotal function in the regulation of eIF4F activity by PABP. In cell extracts, the addition of YB‐1 exacerbated the inhibition of 80S ribosome initiation complex formation by PABP depletion. Rabbit reticulocyte lysate in which PABP weakly stimulates translation is rendered PABP‐dependent after the addition of YB‐1. In this system, eIF4E binding to the cap structure is inhibited by YB‐1 and stimulated by a nonspecific RNA. Significantly, adding PABP back to the depleted lysate stimulated eIF4E binding to the cap structure more potently if this binding had been downregulated by YB‐1. Conversely, adding nonspecific RNA abrogated PABP stimulation of eIF4E binding. These data strongly suggest that competition between YB‐1 and eIF4G for mRNA binding is required for efficient stimulation of eIF4F activity by PABP.     

Protein body formation in leaves of Nicotiana benthamiana: a concentration‐dependent mechanism influenced by the presence of fusion tags

Protein bodies (PBs) are endoplasmic reticulum (ER) derived organelles originally found in seeds whose function is to accumulate seed storage proteins. It has been shown that PB formation is not limited to seeds and green fluorescent protein (GFP) fused to either elastin‐like polypeptide (ELP) or hydrophobin (HFBI) fusion tags induce the formation of PBs in leaves of N. benthamiana. In this study, we compared the ELP‐ and HFBI‐induced PBs and showed that ELP‐induced PBs are larger than HFBI‐induced PBs. The size of ELP‐ and HFBI‐induced PBs increased over time along with the accumulation levels of their fused protein. Our results show that PB formation is a concentration‐dependent mechanism in which proteins accumulating at levels higher than 0.2% of total soluble protein are capable of inducing PBs in vivo. Our results show that the presence of fusion tags is not necessary for the formation of PBs, but affects the distribution pattern and size of PBs. This was confirmed by PBs induced by fluorescent proteins as well as fungal xylanases. We noticed that in the process of PB formation, secretory and ER‐resident molecules are passively sequestered into the lumen of PBs. We propose to use this property of PBs as a tool to increase the accumulation levels of erythropoietin and human interleukin‐10 by co‐expression with PB‐inducing proteins.     

Role of Ubiquitination in Endocytic Trafficking of G‐Protein‐Coupled Receptors

Lysyl ubiquitination has long been known to target cytoplasmic proteins for proteasomal degradation, and there is now extensive evidence that ubiquitination functions in vacuolar/lysosomal targeting of membrane proteins from both the biosynthetic and endocytic pathways. G‐protein‐coupled receptors (GPCRs) represent the largest and most diverse family of membrane proteins, whose function is of fundamental importance both physiologically and therapeutically. In this review, we discuss the role of ubiquitination in the vacuolar/lysosomal downregulation of GPCRs through the endocytic pathway, with a primary focus on lysosomal trafficking in mammalian cells. We will summarize evidence indicating that mammalian GPCRs are regulated by ubiquitin‐dependent mechanisms conserved in budding yeast, and then consider evidence for additional ubiquitin‐dependent and ‐independent regulation that may be specific to animal cells.     

Perinatal nutrition interacts with genetic background to alter behavior in a parent‐of‐origin‐dependent manner in adult Collaborative Cross mice

Previous studies in animal models and humans have shown that exposure to nutritional deficiencies in the perinatal period increases the risk of psychiatric disease. Less well understood is how such effects are modulated by the combination of genetic background and parent‐of‐origin (PO). To explore this, we exposed female mice from 20 Collaborative Cross (CC) strains to protein deficient, vitamin D deficient, methyl donor enriched or standard diet during the perinatal period. These CC females were then crossed to a male from a different CC strain to produce reciprocal F1 hybrid females comprising 10 distinct genetic backgrounds. The adult F1 females were then tested in the open field, light/dark, stress‐induced hyperthermia, forced swim and restraint stress assays. Our experimental design allowed us to estimate effects of genetic background, perinatal diet, PO and their interactions on behavior. Genetic background significantly affected all assessed phenotypes. Perinatal diet exposure interacted with genetic background to affect body weight, basal body temperature, anxiety‐like behavior and stress response. In 8 of 9 genetic backgrounds, PO effects were observed on multiple phenotypes. Additionally, we identified a small number of diet‐by‐PO effects on body weight, stress response, anxiety‐ and depressive‐like behavior. Our data show that rodent behaviors that model psychiatric disorders are affected by genetic background, PO and perinatal diet, as well as interactions among these factors.