Exploring the relationship between hub proteins and drug targets based on GO and intrinsic disorder

Protein–protein interactions (PPIs) play essential roles in many biological processes. In protein–protein interaction networks, hubs involve in numbers of PPIs and may constitute an important source of drug targets. The intrinsic disorder proteins (IDPs) with unstable structures can promote the promiscuity of hubs and also involve in many disease pathways, so they also could serve as potential drug targets. Moreover, proteins with similar functions measured by semantic similarity of gene ontology (GO) terms tend to interact with each other. Here, the relationship between hub proteins and drug targets based on GO terms and intrinsic disorder was explored. The semantic similarities of GO terms and genes between two proteins, and the rate of intrinsic disorder residues of each protein were extracted as features to characterize the functional similarity between two interacting proteins. Only using 8 feature variables, prediction models by support vector machine (SVM) were constructed to predict PPIs. The accuracy of the model on the PPI data from human hub proteins is as high as 83.72%, which is very promising compared with other PPI prediction models with hundreds or even thousands of features. Then, 118 of 142 PPIs between hubs are correctly predicted that the two interacting proteins are targets of the same drugs. The results indicate that only 8 functional features are fully efficient for representing PPIs. In order to identify new targets from IDP dataset, the PPIs between hubs and IDPs are predicted by the SVM model and the model yields a prediction accuracy of 75.84%. Further research proves that 3 of 5 PPIs between hubs and IDPs are correctly predicted that the two interacting proteins are targets of the same drugs. All results demonstrate that the model with only 8-dimensional features from GO terms and intrinsic disorder still gives a good performance in predicting PPIs and further identifying drug targets.     

The d'--d--d' vertical triad is less discriminating than the a'--a--a' vertical triad in the antiparallel coiled-coil dimer motif

Elucidating relationships between the amino-acid sequences of proteins and their three-dimensional structures, and uncovering non-covalent interactions that underlie polypeptide folding, are major goals in protein science. One approach toward these goals is to study interactions between selected residues, or among constellations of residues, in small folding motifs. The α-helical coiled coil has served as a platform for such studies because this folding unit is relatively simple in terms of both sequence and structure. Amino acid side chains at the helix-helix interface of a coiled coil participate in so-called "knobs-into-holes" (KIH) packing whereby a side chain (the knob) on one helix inserts into a space (the hole) generated by four side chains on a partner helix. The vast majority of sequence-stability studies on coiled-coil dimers have focused on lateral interactions within these KIH arrangements, for example, between an a position on one helix and an a' position of the partner in a parallel coiled-coil dimer, or between a--d' pairs in an antiparallel dimer. More recently, it has been shown that vertical triads (specifically, a'--a--a' triads) in antiparallel dimers exert a significant impact on pairing preferences. This observation provides impetus for analysis of other complex networks of side-chain interactions at the helix-helix interface. Here, we describe a combination of experimental and bioinformatics studies that show that d'--d--d' triads have much less impact on pairing preference than do a'--a--a' triads in a small, designed antiparallel coiled-coil dimer. However, the influence of the d'--d--d' triad depends on the lateral a'--d interaction. Taken together, these results strengthen the emerging understanding that simple pairwise interactions are not sufficient to describe side-chain interactions and overall stability in antiparallel coiled-coil dimers; higher-order interactions must be considered as well.     

Compensatory Adaptations of Structural Dynamics in an Intrinsically Disordered Protein Complex

Intrinsically disordered proteins (IDPs) play crucial roles in protein interaction networks and in this context frequently constitute important hubs and interfaces. Here we show by a combination of NMR and EPR spectroscopy that the binding of the cytokine osteopontin (OPN) to its natural ligand, heparin, is accompanied by thermodynamically compensating structural adaptations. The core segment of OPN expands upon binding. This “unfolding‐upon‐binding” is governed primarily through electrostatic interactions between heparin and charged patches along the protein backbone and compensates for entropic penalties due to heparin–OPN binding. It is shown how structural unfolding compensates for entropic losses through ligand binding in IDPs and elucidates the interplay between structure and thermodynamics of rapid substrate‐binding and ‐release events in IDP interaction networks.     

Mapping of two networks of residues that exhibit structural and dynamical changes upon binding in a PDZ domain protein

We describe the changes in structure and dynamics that occur in the second PDZ domain of human tyrosine phosphatase 1E upon binding the small peptide RA-GEF2 by an analysis of NMR data based on their use as ensemble-averaged restraints in molecular dynamics simulations. This approach reveals the presence of two interconnected networks of residues, the first exhibiting structural changes and the second dynamical changes upon binding, and it provides a detailed mapping of the regions of increased and decreased mobility upon binding. Analysis of the dynamical properties of the residues in these networks reveals that conformational changes are transmitted through pathways of coupled side-chain reorientations. These results illustrate how the strategy we described, in which NMR data are used in combination with molecular dynamics simulations, can be used to characterize in detail the complex organization of the changes in structure and dynamics that take place in proteins upon binding.     

Effect of side-chain disorder on the electronic structure of proteins

The influence of side-chain disorder on the electronic structure of proteins has been investigated in the case of polypeptides containing two or three different amino acid residues. It has been found that due to the different potentials of different side-chain groups, the original valence and conduction bands of the homopolypeptides are split into narrow bands. The comparison of the densities of electronic states in simple homopolypeptides and in composite polymers shows that new forbidden regions in the energy spectrum of proteins may develop. The consequences of these effects for the semiconductive properties of proteins are discussed.     

Exploring the Role of C–H….π Interactions on the Structural Stability of Single Chain “All-Alpha” Proteins

C–H….π interactions are known to be important contributors to protein stability. In this study, we have analyzed the influence of C–H….π interactions in single chain “all-alpha” proteins. In the data set, a total of 181 C–H….π interactions were observed. The most prominent representatives are the interactions between aromatic C–H donor groups and aromatic π acceptors. Eighty-one percent of the C–H….π interactions between side chain to side chain and remaining19% of the C–H….π interactions were observed between side-chain to side-chain five-member aromatic ring. The donor atom contribution to C–H….π interactions was mainly from Phe, Tyr, and Trp residues. The acceptor atom contribution to C–H….π interactions was mainly from Phe, Tyr, Trp, and His. The highest percentage of C–H….π interactions were observed form Phe residue. The secondary structure preference analysis of all C–H….π interacting residues showed that Phe, Tyr, Trp, and His preferred to be in helix. Long-range C–H….π interactions are the predominant type of interactions in single chain all-alpha proteins data set. All the C–H….π interactions forming residues in the data set preferred to be in the buried region. Seventy-three percent of the donor residues and 65% of the acceptor residues are highly conserved.     

A united-residue force field for off-lattice protein-structure simulations. I. Functional forms and parameters of long-range side-chain interaction potentials from protein crystal data

A two-stage procedure for the determination of a united-residue potential designed for protein simulations is outlined. In the first stage, the long-range and local-interaction energy terms of the total energy of a polypeptide chain are determined by analyzing protein-crystal data and averaging the all-atom energy surfaces. In the second stage (described in the accompanying article), the relative weights of the energy terms are optimized so as to locate the native structures of selected test proteins as the lowest energy structures. The goal of the work in the present study is to parameterize physically reasonable functional forms of the potentials of mean force for side-chain interactions. The potentials are of both radial and anisotropic type. Radial potentials include the Lennard-Jones and the shifted Lennard-Jones potential (with the shift parameter independent of orientation). To treat the angular dependence of side-chain interactions, three functional forms of the potential that were designed previously to describe anisotropic systems are evaluated: Berne-Pechukas (dilated Lennard-Jones); Gay-Berne (shifted Lennard-Jones with orientation-dependent shift parameters); and Gay-Berne-Vorobjev (the same as the preceding one, but with one more set of variable parameters). These functional forms were used to parameterize, within a short-distance range, the potentials of mean force for side-chain pair interactions that are related by the Boltzmann principle to the pair correlation functions determined from protein-crystal data. Parameter determination was formulated as a generalized nonlinear least-squares problem with the target function being the weighted sum of squares of the differences between calculated and “experimental” (i.e., estimated from protein-crystal data) angular, radial-angular, and radial pair correlation functions, as well as contact free energies. A set of 195 high-resolution nonhomologous structures from the Protein Data Bank was used to calculate the “experimental” values. The contact free energies were scaled by the slope of the correlation line between side-chain hydrophobicities, calculated from the contact free energies, and those determined by Fauchere and Pliška from the partition coefficients of amino acids between water and n-octanol. The methylene group served to define the reference contact free energy corresponding to that between the glycine methylene groups of backbone residues. Statistical analysis of the goodness of fit revealed that the Gay-Berne-Vorobjev anisotropic potential fits best to the experimental radial and angular correlation functions and contact free energies and therefore represents the free-energy surface of side-chain-side-chain interactions most accurately. Thus, its choice for simulations of protein structure is probably the most appropriate. However, the use of simpler functional forms is recommended, if the speed of computations is an issue. © 1997 by John Wiley & Sons, Inc. J Comput Chem 18: 849–873, 1997     

Structure-based secondary structure-independent approach to design protein ligands: Application to the design of Kv1.2 potassium channel blockers

We have developed a structure-based approach to the design of protein ligands. This approach is based on the transfer of a functional binding motif of amino acids, often referred as to the "hot spot", on a host protein able to reproduce the functional topology of these residues. The scaffolds were identified by a systematic in silico search in the Protein Data Bank for proteins possessing a group of residues in a topology similar to that adopted by the functional motif in a reference ligand of known 3D structure. In contrast to previously reported studies, this search is independent of the particular secondary structure supporting the functional motif. To take into account the global properties of the host protein, two additional criteria were taken into account in the selection process: (1) Only those scaffolds sterically compatible with the positioning of the functional motif as observed in a reference complex model were retained. (2) Host proteins displaying electrostatic potentials, in the region of the transferred functional motif, similar to that of the reference ligand were selected. This approach was applied to the development of protein ligands of the Kv1.2 channel using BgK, a small protein isolated from the sea anemone Bunodosoma granulifera, as the reference ligand. Four proteins obtained by this approach were produced for experimental evaluation. The X-ray structure of one of these proteins was determined to check for similarity of the transferred functional motif with the structure it adopts in the reference ligand. Three of these protein ligands bind the Kv1.2 channel with inhibition constants of 0.5, 1.5, and 1.6 microM. Several mutants of these designed protein ligands gave binding results consistent with the presumed binding mode. These results show that protein ligands can be designed by transferring a binding motif on a protein host selected to reproduce the functional topology of this motif, irrespective to the secondary structure supporting the functional motif, if the host protein possesses steric and electrostatic properties compatible with the binding to the target. This result opens the way to the design of protein ligands by taking advantage of the considerable structural repertoire of the Protein Data Bank.     

Differential ordering of the protein backbone and side chains during protein folding revealed by site-specific recombinant infrared probes

The time scale for ordering of the polypeptide backbone relative to the side chains is a critical issue in protein folding. The interplay between ordering of the backbone and ordering of the side chains is particularly important for the formation of β-sheet structures, as the polypeptide chain searches for the native stabilizing cross-strand interactions. We have studied these issues in the N-terminal domain of protein L9 (NTL9), a model protein with mixed α/β structure. We have developed a general approach for introducing site-specific IR probes for the side chains (azide) and backbone ((13)C═(18)O) using recombinant protein expression. Temperature-jump time-resolved IR spectroscopy combined with site-specific labeling enables independent measurement of the respective backbone and side-chain dynamics with single residue resolution. We have found that side-chain ordering in a key region of the β-sheet structure occurs on a slower time scale than ordering of the backbone during the folding of NTL9, likely as a result of the transient formation of non-native side-chain interactions.     

Protein folding using fragment assembly and physical energy function

We perform a systematic study of the effects of sequence-independent backbone interactions and sequence-dependent side-chain interactions on protein folding using fragment assembly and physical energy function. Structures for ten proteins belonging to various structural classes are predicted only with Lennard-Jones interaction between backbone atoms. We find nativelike structures for beta proteins, suggesting that for proteins in this class, the global tertiary structures can be determined mainly by sequence-independent backbone interactions. On the other hand, for alpha proteins, nonlocal hydrophobic side-chain interaction is also required to obtain nativelike structures.     

Residue-specific pKa determination of lysine and arginine side chains by indirect 15N and 13C NMR spectroscopy: application to apo calmodulin

Electrostatic interactions in proteins can be probed experimentally through determination of residue-specific acidity constants. We describe here triple-resonance NMR techniques for direct determination of lysine and arginine side-chain protonation states in proteins. The experiments are based on detection of nonexchangeable protons over the full range of pH and temperature and therefore are well suited for pKa determination of individual amino acid side chains. The experiments follow the side-chain 15Nzeta (lysine) and 15Nepsilon or 13Czeta (arginine) chemical shift, which changes due to sizable changes in the heteronuclear electron distribution upon (de)protonation. Since heteronuclear chemical shifts are overwhelmed by the charge state of the amino acid side chain itself, these methods supersede 1H-based NMR in terms of accuracy, sensitivity, and selectivity. Moreover, the 15Nzeta and 15Nepsilon nuclei may be used to probe changes in the local electrostatic environment. Applications to three proteins are described: apo calmodulin, calbindin D9k, and FKBP12. For apo calmodulin, residue-specific pKa values of lysine side chains were determined to fall between 10.7 and 11.2 as a result of the high net negative charge on the protein surface. Ideal two-state titration behavior observed for all lysines indicates the absence of significant direct charge interactions between the basic residues. These results are compared with earlier studies based on chemical modification.     

Measurement of side-chain carboxyl pK(a) values of glutamate and aspartate residues in an unfolded protein by multinuclear NMR spectroscopy

A triple-resonance NMR pulse scheme is presented for measuring aspartic and glutamic acid side-chain pK(a) values in unfolded protein states where chemical shift overlap is limiting. The experiment correlates side-chain carboxyl carbon chemical shifts of these residues with the backbone amide proton chemical shift of the following residue. The methodology is applied to an (15)N, (13)C labeled sample of the N-terminal SH3 domain of the Drosophila protein drk, which exists in equilibrium between folded (F(exch)) and unfolded (U(exch)) states under nondenaturing conditions. Residue-specific pK(a) values of side-chain carboxyl groups are presented for the first time for an unfolded protein (drk U(exch) state), determined from a pH titration. Results indicate that deviations from pK(a) values measured for model compounds are likely due to local effects, while long-range electrostatic interactions appear to be of minor importance for this protein.     

A conserved helical capping hydrogen bond in PAS domains controls signaling kinetics in the superfamily prototype photoactive yellow protein

PAS domains form a divergent protein superfamily with more than 20?000 members that perform a wide array of sensing and regulatory functions in all three domains of life. Only nine residues are well-conserved in PAS domains, with an Asn residue at the start of α-helix 3 showing the strongest conservation. The molecular functions of these nine conserved residues are unknown. We use static and time-resolved visible and FTIR spectroscopy to investigate receptor activation in the photosensor photoactive yellow protein (PYP), a PAS domain prototype. The N43A and N43S mutants allow an investigation of the role of side-chain hydrogen bonding at this conserved position. The mutants exhibit a blue-shifted visible absorbance maximum and up-shifted chromophore pK(a). Disruption of the hydrogen bonds in N43A PYP causes both a reduction in protein stability and a 3400-fold increase in the lifetime of the signaling state of this photoreceptor. A significant part of this increase in lifetime can be attributed to the helical capping interaction of Asn43. This extends the known importance of helical capping for protein structure to regulating functional protein kinetics. A model for PYP activation has been proposed in which side-chain hydrogen bonding of Asn43 is critical for relaying light-induced conformational changes. However, FTIR spectroscopy shows that both Asn43 mutants retain full allosteric transmission of structural changes. Analysis of 30 available high-resolution structures of PAS domains reveals that the side-chain hydrogen bonding of residue 43 but not residue identity is highly conserved and suggests that its helical cap affects signaling kinetics in other PAS domains.     

CAMWI: Detecting protein complexes using weighted clustering coefficient and weighted density

Detection of protein complexes is very important to understand the principles of cellular organization and function. Recently, large protein–protein interactions (PPIs) networks have become available using high-throughput experimental techniques. These networks make it possible to develop computational methods for protein complex detection. Most of the current methods rely on the assumption that protein complex as a module has dense structure. However complexes have core-attachment structure and proteins in a complex core share a high degree of functional similarity, so it expects that a core has high weighted density. In this paper we present a Core-Attachment based method for protein complex detection from Weighted PPI Interactions using clustering coefficient and weighted density. Experimental results show that the proposed method, CAMWI improves the accuracy of protein complex detection.     

Arginine Side-Chain Hydrogen Exchange: Quantifying Arginine Side-Chain Interactions in Solution

The rate with which labile backbone hydrogen atoms in proteins exchange with the solvent has long been used to probe protein interactions in aqueous solutions. Arginine, an essential amino acid found in many interaction interfaces, is capable of an impressive range of interactions via its guanidinium group. The hydrogen exchange rate of the guanidinium hydrogens therefore becomes an important measure to quantify side-chain interactions. Herein we present an NMR method to quantify the hydrogen exchange rates of arginine side-chain ¹H^ϵ protons and thus present a method to gauge the strength of arginine side-chain interactions. The method employs ¹³C-detection and the one-bond deuterium isotope shift observed for ¹⁵N^ϵ to generate two exchanging species in ¹H₂O/²H₂O mixtures. An application to the protein T4 Lysozyme is shown, where protection factors calculated from the obtained exchange rates correlate well with the interactions observed in the crystal structure. The methodology presented provides an important step towards characterising interactions of arginine side-chains in enzymes, in phase separation, and in protein interaction interfaces in general.     

Noncovalent cross-links in context with other structural and functional elements of proteins

Proteins are heteropolymers with evolutionary selected native sequences of residues. These native sequences code for unique and stable 3D structures indispensable for biochemical activity and for proteolysis resistance, the latter which guarantees an appropriate lifetime for the protein in the protease rich cellular environment. Cross-links between residues close in space but far in the primary structure are required to maintain the folded structure of proteins. Some of these cross-links are covalent, most frequently disulfide bonds, but the majority of the cross-links are sets of cooperative noncovalent long-range interactions. In this paper we focus on special clusters of noncovalent long-range interactions: the Stabilization Centers (SCs). The relation between the SCs and secondary structural elements as well as the relation between SCs and functionally important regions of proteins are presented to show a detailed picture of these clusters, which are believed to be primarily responsible for major aspects of protein stability.     

Structural and Dynamic Characterization of the Molecular Hub Early Region 1A (E1A) from Human Adenovirus

The small‐DNA human adenovirus encodes one of the most versatile molecular hubs, the E1A protein. This protein is essential for productive viral infection in human cells and a vast amount of biologically relevant data are available on its interactions with host proteins. Up to now, however, no high‐resolution structural and dynamic information on E1A is available despite its important biological role. Among the different spliced variants of E1A, two are expressed at high level in the early stage of infection. These are 243 and 289 residues isoforms. Herein, we present their NMR characterization, showing that they are both highly disordered, but also demonstrate a certain heterogeneous behavior in terms of structural and dynamic properties. Furthermore, we present the characterization of the isolated domain of the longer variant, known as CR3. This study opens the way to understanding at the molecular level how E1A functions.     

A new method for side-chain conformation prediction using a Hopfield network and reproduced rotamers

We present a new side-chain prediction method based on energy minimization using a Hopfield network, focusing on the buried residues of proteins. In this method, the network is composed of automata assigned to each rotamer to restrict side-chain conformational space. We reproduced a rotamer library that enabled us to more widely cover the space for side-chain conformations than those previously produced. The accuracy of the side-chain modeling was estimated by three standards: root mean square deviations (rmsds) between the modeled and the crystal structures, the percentages of correctly predicted side-chain torsion angles, and the percentages of correctly predicted hydrogen bonds. The average rmsd for buried side chains of 21 proteins was 1.10 Å. The value was almost always improved relative to the previous works. The percentage of side-chain X₁ angles for buried residues was 87.3%. By considering the hydrogen bond energy, the average percentage of correctly predicted hydrogen bonds rose from 33% without hydrogen bond energy to 52% with the bond energy. We applied this method to homology modeling, where the protein backbone used to predict side-chain conformations deviates from the correct conformation, and could predict side-chain conformations as correctly as those using the correct backbones. © 1996 by John Wiley & Sons, Inc.     

Photocaged Quinone Methide Crosslinkers for Light‐Controlled Chemical Crosslinking of Protein–Protein and Protein–DNA Complexes

Small‐molecule crosslinkers are invaluable for probing biomolecular interactions and for crosslinking mass spectrometry. Existing chemical crosslinkers target only a small selection of amino acids, while conventional photo‐crosslinkers target almost all residues non‐specifically, complicating data analysis. Herein, we report photocaged quinone methide (PQM)‐based crosslinkers that target nine nucleophilic residues through Michael addition, including Gln, Arg, and Asn, which are inaccessible to existing chemical crosslinkers. PQM crosslinkers were used in vitro, in Escherichia coli, and in mammalian cells to crosslink dimeric proteins and endogenous membrane receptors. The heterobifunctional crosslinker NHQM could crosslink proteins to DNA, for which few crosslinkers exist. The photoactivatable reactivity of these crosslinkers and their ability to target multiple amino acids will enhance the use of chemical crosslinking for studies of protein–protein and protein–DNA networks and for structural biology.