Cation-π interactions in lipocalins: structural and functional implications

The cation-π interaction impacts protein folding, structural stability, specificity, and molecular recognition. Cation-π interactions have been overlooked in the lipocalin family. To fill this gap, these interactions were analyzed in the 113 crystal and solution structures from the lipocalin family. The cation-π interactions link previously identified structurally conserved regions and reveal new motifs, which are beyond the reach of a sequence alignment algorithm. Functional and structural significance of the interactions were tested experimentally in human tear lipocalin (TL). TL, a prominent and promiscuous lipocalin, has a key role in lipid binding at the ocular surface. Ligand binding modulation through the loop AB at the "open" end of the barrel has been erroneously attributed solely to electrostatic interactions. Data revealed that the interloop cation-π interaction in the pair Phe28-Lys108 contributes significantly to stabilize the holo-conformation of the loop AB. Numerous energetically significant and conserved cation-π interactions were uncovered in TL and throughout the lipocalin family. Cation-π interactions, such as the highly conserved Trp17-Arg118 pair in TL, were educed in low temperature experiments of mutants with Trp to Tyr substitutions.     

Composite structural motifs of binding sites for delineating biological functions of proteins

Most biological processes are described as a series of interactions between proteins and other molecules, and interactions are in turn described in terms of atomic structures. To annotate protein functions as sets of interaction states at atomic resolution, and thereby to better understand the relation between protein interactions and biological functions, we conducted exhaustive all-against-all atomic structure comparisons of all known binding sites for ligands including small molecules, proteins and nucleic acids, and identified recurring elementary motifs. By integrating the elementary motifs associated with each subunit, we defined composite motifs that represent context-dependent combinations of elementary motifs. It is demonstrated that function similarity can be better inferred from composite motif similarity compared to the similarity of protein sequences or of individual binding sites. By integrating the composite motifs associated with each protein function, we define meta-composite motifs each of which is regarded as a time-independent diagrammatic representation of a biological process. It is shown that meta-composite motifs provide richer annotations of biological processes than sequence clusters. The present results serve as a basis for bridging atomic structures to higher-order biological phenomena by classification and integration of binding site structures.     

Effect of water coordination on competition between π and non-π cation binding sites in aromatic amino acids: l-phenylalanine, l-tyrosine, and l-tryptophan Li+, Na+, and K+ complexes

Quantum chemistry methods have been applied to charged complexes of the alkali metals Li(+), Na(+), and K(+) with the aromatic amino acids (AAAs) phenylalanine (Phe), tyrosine (Tyr), and tryptophan (Trp). The geometries of 72 different complexes (Phe·M, Tyr·M, Trp·M, M is Li(+), Na(+), or K(+)) were completely optimized at the B3LYP/6-311+G(d,p) level of density functional theory. The solvent effect on the geometry and stability of individual complexes was studied by making use of a microsolvation model. The interaction enthalpies, entropies, and Gibbs energies of nine different complexes of the systems Phe·M, Tyr·M, and Trp·M (M is Li(+), Na(+), or K(+)) were also determined at the B3LYP density functional level of theory. The calculated Gibbs binding energies of the M(+)-AAA complexes follow the order Phe < Tyr < Trp for all three metal cations studied. Among the three AAAs studied, the indole ring of Trp is the best π donor for alkali metal cations. Our calculations demonstrated the existence of strong cation-π interactions between the alkali metals and the aromatic side chains of the three AAAs. These AAAs comprise about 8% of all known protein sequences. Thus, besides the potential for hydrogen-bond interaction, aromatic residues of Phe, Tyr, and Trp show great potential for π-donor interactions. The existence of cation-π interaction in proteins has also been demonstrated experimentally. However, more complex experimental studies of metal cation-π interaction in diverse biological systems will no doubt lead to more exact validation of these investigations.     

Very strong metal ligand aromatic cation-π interactions in transition metal complexes: intermolecular interaction in tetraphenylborate salts

Strong intermolecular interactions between ligands in cationic metal complexes and aromatic rings of tetraphenylborate anion, so-called metal ligand aromatic cation-π (MLACπ) interactions, were found by screening Cambridge Structural Database. Distances between phenyl ring and ligand are shorter in these structures than in previously reported MLACπ interactions by 0.2 Å.     

Geometry of interaction of the histidine ring with other planar and basic residues

Among the aromatic residues in protein structures, histidine (His) is unique, as it can exist in the neutral or positively charged form at the physiological pH. As such, it can interact with other aromatic residues as well as form hydrogen bonds with polar and charged (both negative and positive) residues. We have analyzed the geometry of interaction of His residues with nine other planar side chains containing aromatic (residues Phe, Tyr, Trp, and His), carboxylate (Asp and Glu), carboxamide (Asn and Gln) and guanidinium (Arg) groups in 432 polypeptide chains. With the exception of the aspartic (Asp) and glutamic (Glu) acid side-chains, all other residues prefer to interact in a face-to-face or offset-face-stacked orientation with the His ring. Such a geometry is different from the edge-to-face relative orientation normally associated with the aromatic-aromatic interaction. His-His pair prefers to interact in a face-to-face orientation; however, when both the residues bind the same metal ion, the interplanar angle is close to 90 degrees. The occurrence of different interactions (including the nonconventional N-H...pi and C-H...pi hydrogen bonds) have been correlated with the relative orientations between the interacting residues. Several structural motifs, mostly involved in binding metal ions, have been identified by considering the cases where His residues are in contact with four other planar moieties. About 10% of His residues used here are also found in sequence patterns in PROSITE database. There are examples of the amino end of the Lys side chain interacting with His residues in such a way that it is located on an arc around a ring nitrogen atom.     

Aromatic residues link binding and function of intrinsically disordered proteins

We have searched for intermolecular aromatic pairs in 77 protein-protein complexes of intrinsically disordered proteins (IDPs) to understand the role of π-π interactions in protein-protein interactions involving IDPs. We found that 40% of the complexes possess at least one intermolecular pair of aromatic residues. Analysis of composition, characteristics, location and the contribution to the free energy of binding showed that π-π interactions substantially contribute to binding by working as anchor residues, conformational locks, and ready-made recognition motifs required for specific binding. By using available experimental data we show that π-π interactions play a variety of roles that link binding of IDPs and their function in the cell. The results presented in this study pave the way to understand in atomic detail the inner workings of IDPs interaction networks.     

Interaction of DNA with clusters of amino acids in proteins

Protein-DNA interactions facilitate the fundamental functions of living cells and are universal in all living organisms. Several investigations have been carried out, essentially identifying pairs of interactions between the amino acid residues in proteins and the bases in DNA. In the present study, we have detected the recognition motifs that may constitute a cluster of spatially interacting residues in proteins, which interact with the bases of DNA. Graph spectral algorithm has been used to detect side chain clusters comprising Arg, Lys, Asn, Gln and aromatic residues from proteins interacting with DNA. We find that the interaction of proteins with DNA is through clusters in about half of the proteins in the dataset and through individual residues in the rest. Furthermore, inspection of the clusters has revealed additional interactions in a few cases, which have not been reported earlier. The geometry of the interaction between the DNA base and the protein residue is quantified by the distance d and the angle theta. These parameters have been identified for the cation-pi/H-bond stair motif that was reported earlier. Among the Arg, Lys, Asn and Gln residues, the range of (d, theta) values of the interacting Arg clearly falls into the cation-pi and the hydrogen bond interactions of the 'cation-pi/H-bond' stair motif. Analysis of the cluster composition reveals that the Arg residue is predominant than the Lys, Asn and Gln residues. The clusters are classified into Type I and Type II based on the presence or absence of aromatic residues (Phe, Tyr) in them. Residue conservation in these clusters has been examined. Apart from the conserved residues identified previously, a few more residues mainly Phe, Tyr and Arg have also been identified as conserved and interactive with the DNA. Interestingly, a few residues that are parts of interacting clusters and do not interact directly with the DNA have also been conserved. This emphasizes the importance of recognizing the protein side chain cluster motifs interacting with the DNA, which could serve as signatures of protein-DNA recognition in the families of DNA binding proteins.     

The fragment transformation method to detect the protein structural motifs

To identify functional structural motifs from protein structures of unknown function becomes increasingly important in recent years due to the progress of the structural genomics initiatives. Although certain structural patterns such as the Asp-His-Ser catalytic triad are easy to detect because of their conserved residues and stringently constrained geometry, it is usually more challenging to detect a general structural motifs like, for example, the betabetaalpha-metal binding motif, which has a much more variable conformation and sequence. At present, the identification of these motifs usually relies on manual procedures based on different structure and sequence analysis tools. In this study, we develop a structural alignment algorithm combining both structural and sequence information to identify the local structure motifs. We applied our method to the following examples: the betabetaalpha-metal binding motif and the treble clef motif. The betabetaalpha-metal binding motif plays an important role in nonspecific DNA interactions and cleavage in host defense and apoptosis. The treble clef motif is a zinc-binding motif adaptable to diverse functions such as the binding of nucleic acid and hydrolysis of phosphodiester bonds. Our results are encouraging, indicating that we can effectively identify these structural motifs in an automatic fashion. Our method may provide a useful means for automatic functional annotation through detecting structural motifs associated with particular functions.     

Structural basis for copper transfer by the metallochaperone for the Menkes/Wilson disease proteins

The Hah1 metallochaperone protein is implicated in copper delivery to the Menkes and Wilson disease proteins. Hah1 and the N-termini of its target proteins belong to a family of metal binding domains characterized by a conserved MT/HCXXC sequence motif. The crystal structure of Hah1 has been determined in the presence of Cu(I), Hg(II), and Cd(II). The 1.8 A resolution structure of CuHah1 reveals a copper ion coordinated by Cys residues from two adjacent Hah1 molecules. The CuHah1 crystal structure is the first of a copper chaperone bound to copper and provides structural support for direct metal ion exchange between conserved MT/HCXXC motifs in two domains. The structures of HgHah1 and CdHah1, determined to 1.75 A resolution, also reveal metal ion coordination by two MT/HCXXC motifs. An extended hydrogen bonding network, unique to the complex of two Hah1 molecules, stabilizes the metal binding sites and suggests specific roles for several conserved residues. Taken together, the structures provide models for intermediates in metal ion transfer and suggest a detailed molecular mechanism for protein recognition and metal ion exchange between MT/HCXXC containing domains.     

阳离子-π相互作用在两种典型蛋白质折叠型中的偏好性

阳离子-π相互作用是一种在阳离子和芳香性体系之间形成的一种作用力。在蛋白质中,带正电荷的氨基酸(Lys、Arg)和芳香族氨基酸(Phe、Tyr、Trp)之间可以形成阳离子-π相互作用。对α/β类蛋白中两种典型折叠类型――单绕和双绕的研究表明:(1)单绕结构中阳离子-π相互作用的分布密度是双绕结构的2.3倍。(2)Arg-Phe组合偏好在双绕中出现,Arg-Tyr组合偏好在单绕中出现。(3)在单绕中除Lys-Phe组合外,其余5种组合的阳离子-π相互作用能量高于双绕的对应组合,其中以Arg-Trp组合的能量最高。(4)在单绕结构中,样本所含氨基酸残基数量和样本中阳离子-π的数量有明显的相关性,在双绕结构中没有发现类似的相关性。(5)在单绕和双绕结构当中,把阳离子-π相互作用能量分解为静电能和范德华力能揭示出静电能与范德华力能之比接近2∶1,静电作用在阳离子-π相互作用中起主要作用。     

Using RNA secondary structures to guide sequence motif finding towards single-stranded regions

RNA binding proteins recognize RNA targets in a sequence specific manner. Apart from the sequence, the secondary structure context of the binding site also affects the binding affinity. Binding sites are often located in single-stranded RNA regions and it was shown that the sequestration of a binding motif in a double-strand abolishes protein binding. Thus, it is desirable to include knowledge about RNA secondary structures when searching for the binding motif of a protein. We present the approach MEMERIS for searching sequence motifs in a set of RNA sequences and simultaneously integrating information about secondary structures. To abstract from specific structural elements, we precompute position-specific values measuring the single-strandedness of all substrings of an RNA sequence. These values are used as prior knowledge about the motif starts to guide the motif search. Extensive tests with artificial and biological data demonstrate that MEMERIS is able to identify motifs in single-stranded regions even if a stronger motif located in double-strand parts exists. The discovered motif occurrences in biological datasets mostly coincide with known protein-binding sites. This algorithm can be used for finding the binding motif of single-stranded RNA-binding proteins in SELEX or other biological sequence data.     

Water and side‐chain embedded π‐turns

Elucidating protein function from its structure is central to the understanding of cellular mechanisms. This involves deciphering the dependence of local structural motifs on sequence. These structural motifs may be stabilized by direct or water‐mediated hydrogen bonding among the constituent residues. π‐Turns, defined by interactions between (i) and (i + 5) positions, are large enough to contain a central space that can embed a water molecule (or a protein moiety) to form a stable structure. This work is an analysis of such embedded π‐turns using a nonredundant dataset of protein structures. A total of 2965 embedded π‐turns have been identified, as also 281 embedded Schellman motif, a type of π‐turn which occurs at the C‐termini of α‐helices. Embedded π‐turns and Schellman motifs have been classified on the basis of the protein atoms of the terminal turn residues that are linked by the embedded moiety, conformation, residue composition, and compared with the turns that have terminal residues connected by direct hydrogen bonds. Geometrically, the turns have been fitted to a circle and the position of the linker relative to its center analyzed. The hydroxyl group of Ser and Thr, located at (i + 3) position, is the most prominent linker for the side‐chain mediated π‐turns. Consideration of residue conservation among homologous sequences indicates the terminal and the linker positions to be the most conserved. The embedded π‐turn as a binding site (for the linker) is discussed in the context of “nest,” a concave depression that is formed in protein structures with adjacent residues having enantiomeric main‐chain conformations. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 441–453, 2014.     

Decomposing the space of protein quaternary structures with the interface fragment pair library

Background

The physical interactions between proteins constitute the basis of protein quaternary structures. They dominate many biological processes in living cells. Deciphering the structural features of interacting proteins is essential to understand their cellular functions. Similar to the space of protein tertiary structures in which discrete patterns are clearly observed on fold or sub-fold motif levels, it has been found that the space of protein quaternary structures is highly degenerate due to the packing of compact secondary structure elements at interfaces. Therefore, it is necessary to further decompose the protein quaternary structural space into a more local representation.

Results

Here we constructed an interface fragment pair library from the current structure database of protein complexes. After structural-based clustering, we found that more than 90% of these interface fragment pairs can be represented by a limited number of highly abundant motifs. These motifs were further used to guide complex assembly. A large-scale benchmark test shows that the native-like binding is highly likely in the structural ensemble of modeled protein complexes that were built through the library.

Conclusions

Our study therefore presents supportive evidences that the space of protein quaternary structures can be represented by the combination of a small set of secondary-structure-based packing at binding interfaces. Finally, after future improvements such as adding sequence profiles, we expect this new library will be useful to predict structures of unknown protein-protein interactions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0437-4) contains supplementary material, which is available to authorized users.     

Cation-π Interactions and their Functional Roles in Membrane Proteins

Cation-π interactions arise as a result of strong attractive forces between positively charged entities and the π-electron cloud of aromatic groups. The physicochemical characteristics of cation-π interactions are particularly well-suited to the dual hydrophobic/hydrophilic environment of membrane proteins. As high-resolution structural data of membrane proteins bring molecular features into increasingly sharper view, cation-π interactions are gaining traction as essential contributors to membrane protein chemistry, function, and pharmacology. Here we review the physicochemical properties of cation-π interactions and present several prominent examples which demonstrate significant roles for this specialized biological chemistry.     

Geometry of nonbonded interactions involving planar groups in proteins

Although hydrophobic interaction is the main contributing factor to the stability of the protein fold, the specificity of the folding process depends on many directional interactions. An analysis has been carried out on the geometry of interaction between planar moieties of ten side chains (Phe, Tyr, Trp, His, Arg, Pro, Asp, Glu, Asn and Gln), the aromatic residues and the sulfide planes (of Met and cystine), and the aromatic residues and the peptide planes within the protein tertiary structures available in the Protein Data Bank. The occurrence of hydrogen bonds and other nonconventional interactions such as C–H⋯π, C–H⋯O, electrophile–nucleophile interactions involving the planar moieties has been elucidated. The specific nature of the interactions constraints many of the residue pairs to occur with a fixed sequence difference, maintaining a sequential order, when located in secondary structural elements, such as α-helices and β-turns. The importance of many of these interactions (for example, aromatic residues interacting with Pro or cystine sulfur atom) is revealed by the higher degree of conservation observed for them in protein structures and binding regions. The planar residues are well represented in the active sites, and the geometry of their interactions does not deviate from the general distribution. The geometrical relationship between interacting residues provides valuable insights into the process of protein folding and would be useful for the design of protein molecules and modulation of their binding properties.     

Contextual specificity in peptide-mediated protein interactions

Most biological processes are regulated through complex networks of transient protein interactions where a globular domain in one protein recognizes a linear peptide from another, creating a relatively small contact interface. Although sufficient to ensure binding, these linear motifs alone are usually too short to achieve the high specificity observed, and additional contacts are often encoded in the residues surrounding the motif (i.e. the context). Here, we systematically identified all instances of peptide-mediated protein interactions of known three-dimensional structure and used them to investigate the individual contribution of motif and context to the global binding energy. We found that, on average, the context is responsible for roughly 20% of the binding and plays a crucial role in determining interaction specificity, by either improving the affinity with the native partner or impeding non-native interactions. We also studied and quantified the topological and energetic variability of interaction interfaces, finding a much higher heterogeneity in the context residues than in the consensus binding motifs. Our analysis partially reveals the molecular mechanisms responsible for the dynamic nature of peptide-mediated interactions, and suggests a global evolutionary mechanism to maximise the binding specificity. Finally, we investigated the viability of non-native interactions and highlight cases of potential cross-reaction that might compensate for individual protein failure and establish backup circuits to increase the robustness of cell networks.     

A phage display-based method for determination of relative affinities of mutants. Application of the actin-binding motifs in thymosin beta 4 and the villin headpiece

We propose phage display combined with enzyme-linked immunosorbent assay as a tool for the systematic analysis of protein-protein interactions by investigating the binding behavior of variants to a partner protein. Via enzyme-linked immunosorbent assay we determine both the amount of fusion protein presented at the phage surface and the amount of complex formed, the ratio of which is proportional to the affinity. Hence this method enables us to calculate the relative affinities of a large number of mutants. As model systems, we investigated actin-binding motifs conserved in a number of proteins binding monomeric or filamentous actin. The hexapeptide motifs LKKTET, present in thymosin beta4, and LKKEKG, present in the villin headpiece, were mutated, and the variants were analyzed. Study of the positional tolerance allows postulating that the motifs, although similar in primary structures adopt different conformations when bound to actin. In addition, our data show that the second and the fourth amino acid of the thymosin beta4 motif and the first three residues of the villin headpiece motif are most important for actin binding. The latter result challenges the charged crown hypothesis for the villin headpiece filamentous actin interaction.     

Mutagenesis and biochemical studies on AuaA confirmed the importance of the two conserved aspartate-rich motifs and suggested difference in the amino acids for substrate binding in membrane-bound prenyltransferases

AuaA is a membrane-bound farnesyltransferase from the myxobacterium Stigmatella aurantiaca involved in the biosynthesis of aurachins. Like other known membrane-bound aromatic prenyltransferases, AuaA contains two conserved aspartate-rich motifs. Several amino acids in the first motif NXxxDxxxD were proposed to be responsible for prenyl diphosphate binding via metal ions like Mg(2+). Site-directed mutagenesis experiments demonstrated in this study that asparagine, but not the arginine residue in NRxxDxxxD, is important for the enzyme activity of AuaA, differing from the importance of NQ or ND residues in the NQxxDxxxD or NDxxDxxxD motifs observed in some membrane-bound prenyltransferases. The second motif of known membrane-bound prenyltransferases was proposed to be involved in the binding of their aromatic substrates. KDIxDxEGD, also found in AuaA, had been previously speculated to be characteristic for binding of flavonoids or homogenisate. Site-directed mutagenesis experiments with AuaA showed that KDIxDxEGD was critical for the enzyme activity. However, this motif is very likely not specific for flavonoid or homogenisate prenyltransferases, because none of the tested flavonoids was accepted by AuaA or its mutant R53A in the presence of farnesyl, geranyl or dimethylallyl diphosphate.