Protein-RNA interactions: structural analysis and functional classes

A data set of 89 protein-RNA complexes has been extracted from the Protein Data Bank, and the nucleic acid recognition sites characterized through direct contacts, accessible surface area, and secondary structure motifs. The differences between RNA recognition sites that bind to RNAs in functional classes has also been analyzed. Analysis of the complete data set revealed that van der Waals interactions are more numerous than hydrogen bonds and the contacts made to the nucleic acid backbone occur more frequently than specific contacts to nucleotide bases. Of the base-specific contacts that were observed, contacts to guanine and adenine occurred most frequently. The most favored amino acid-nucleotide pairings observed were lysine-phosphate, tyrosine-uracil, arginine-phosphate, phenylalanine-adenine and tryptophan-guanine. The amino acid propensities showed that positively charged and polar residues were favored as expected, but also so were tryptophan and glycine. The propensities calculated for the functional classes showed trends similar to those observed for the complete data set. However, the analysis of hydrogen bond and van der Waal contacts showed that in general proteins complexed with messenger RNA, transfer RNA and viral RNA have more base specific contacts and less backbone contacts than expected, while proteins complexed with ribosomal RNA have less base-specific contacts than the expected. Hence, whilst the types of amino acids involved in the interfaces are similar, the distribution of specific contacts is dependent upon the functional class of the RNA bound.     

Effective discrimination between biologically relevant contacts and crystal packing contacts using new determinants

In the structural models determined by X‐ray crystallography, contacts between molecules can be divided into two categories: biologically relevant contacts and crystal packing contacts. With the growth in the number and quality of available large crystal packing contacts structures, distinguishing crystal packing contacts from biologically relevant contacts remains a difficult task, which can lead to wrong interpretation of structural models. In this study, we performed a systematic analysis on the biologically relevant contacts and crystal packing contacts. The analysis results reveal that biologically contacts are more tightly packed than crystal packing contacts. This property of biologically contacts may contribute to the formation of their interfacial core region. Meanwhile, the differences between the core and surface region of biologically contacts in amino acid composition and evolutionary measure are more dramatic than crystal packing contacts and these differences appear to be useful in distinguishing these two categories of contacts. On the basis of the features derived from our analysis, we developed a random forest model to classify biological relevant contacts and crystal packing contacts. Our method can achieve a high receiver operating curve of 0.923 in the 5‐fold cross‐validation and accuracies of 91.4% and 91.7% for two different test sets. Moreover, in a comparison study, our model outperforms other existing methods, such as DiMoVo, Pita, Pisa, and Eppic. We believe that this study will provide useful help in the validation of oligomeric proteins and protein complexes. The model and all data used in this paper are freely available at http://cic.scu.edu.cn/bioinformatics/bio‐cry.zip . Proteins 2014; 82:3090–3100. © 2014 Wiley Periodicals, Inc.     

Non-randomness in side-chain packing: the distribution of interplanar angles

We analyze the distributions of interplanar angles between interacting side chains with well-defined planar regions, to see whether these distributions correspond to random packing or alternatively show orientational preferences. We use a non-homologous set of 79 high-resolution protein chain structures to show that the observed distributions are significantly different from the sinusoidal one expected for random packing. Overall, we see a relative excess of small angles and a paucity of large interplanar angles; the difference between the expected and observed distributions can be described as a shift of 5% of the interplanar angles from large (≥60°) to small (<30°) values. By grouping the residue pairs into categories based on chemical similarity, we find that some categories have very non-sinusoidal interplanar angle distributions, whereas other categories have distributions that are close to sinusoidal. For a few categories, observed deviations from a sinusoidal distribution can be explained by the electrostatic anisotropy of the isolated pair potential energy. In other cases, the observed distributions reflect the longer range effects of different possible interaction geometries. In particular, geometries that disrupt external hydrogen bonding are disfavored. Proteins 29:370–380, 1997. © 1997 Wiley-Liss, Inc.     

Hydration of protein-protein interfaces

We present an analysis of the water molecules immobilized at the protein-protein interfaces of 115 homodimeric proteins and 46 protein-protein complexes, and compare them with 173 large crystal packing interfaces representing nonspecific interactions. With an average of 15 waters per 1000 A2 of interface area, the crystal packing interfaces are more hydrated than the specific interfaces of homodimers and complexes, which have 10-11 waters per 1000 A2, reflecting the more hydrophilic composition of crystal packing interfaces. Very different patterns of hydration are observed: Water molecules may form a ring around interfaces that remain "dry," or they may permeate "wet" interfaces. A majority of the specific interfaces are dry and most of the crystal packing interfaces are wet, but counterexamples exist in both categories. Water molecules at interfaces form hydrogen bonds with protein groups, with a preference for the main-chain carbonyl and the charged side-chains of Glu, Asp, and Arg. These interactions are essentially the same in specific and nonspecific interfaces, and very similar to those observed elsewhere on the protein surface. Water-mediated polar interactions are as abundant at the interfaces as direct protein-protein hydrogen bonds, and they may contribute to the stability of the assembly.     

Brownian dynamics of interactions between aldolase mutants and F-actin

Previous Brownian dynamics (BD) simulations (Ouporov IG, Knull HR and Thomasson KA 1999. Biophys. J. 76: 17-27) of complex formation between rabbit aldolase and F-actin have identified three lysine residues (K288, K293 and K341) on aldolase and acidic residues (DEDE) at the N-terminus of actin as important to binding. BD simulations of computer models of aldolase mutants with any of these lysine residues replaced by alanine show reduced binding energy; the greatest effect of a single substitution is for K341A, and replacement of all three lysines greatly reduces binding. BD simulations of wild-type rabbit aldolase vs altered F-actin show that binding is decreased if any one of the four N-terminal acidic residues is replaced by alanine and binding is greatly reduced if three or more of the N-terminal acidic residues are replaced; none of the four actin residues appear more critical for binding than the others.     

Brownian dynamics simulations of glycolytic enzyme subsets with F-actin

Previous Brownian dynamics (BD) simulations identified specific basic residues on fructose-1,6-bisphophate aldolase (aldolase) (I. V. Ouporov et al., Biophysical Journal, 1999, Vol. 76, pp. 17-27) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) (I. V. Ouporov et al., Journal of Molecular Recognition, 2001, Vol. 14, pp. 29-41) involved in binding F-actin, and suggested that the quaternary structure of the enzymes may be important. Herein, BD simulations of F-actin binding by enzyme dimers or peptides matching particular sequences of the enzyme and the intact enzyme triose phosphate isomerase (TIM) are compared. BD confirms the experimental observation that TIM has little affinity for F-actin. For aldolase, the critical residues identified by BD are found in surface grooves, formed by subunits A/D and B/C, where they face like residues of the neighboring subunit enhancing their electrostatic potentials. BD simulations between F-actin and aldolase A/D dimers give results similar to the native tetramer. Aldolase A/B dimers form complexes involving residues that are buried in the native structure and are energetically weaker; these results support the importance of quaternary structure for aldolase. GAPDH, however, placed the critical residues on the corners of the tetramer so there is no enhancement of the electrostatic potential between the subunits. Simulations using GAPDH dimers composed of either S/H or G/H subunits show reduced binding energetics compared to the tetramer, but for both dimers, the sets of residues involved in binding are similar to those found for the native tetramer. BD simulations using either aldolase or GAPDH peptides that bind F-actin experimentally show complex formation. The GAPDH peptide bound to the same F-actin domain as did the intact tetramer; however, unlike the tetramer, the aldolase peptide lacked specificity for binding a single F-actin domain.     

Conformation and packing properties of phosphatidic acid: The crystal structure of monosodium dimyristoylphosphatidate

The conformation and molecular packing of monosodium 1,2-dimyristoyl-sn-glycerophosphate (DMPA) has been determined by single crystal analysis (R = 0.107). The lipid crystallizes in the space group P2₁ with unit cell dimensions: $\text{a = 5.44, b = 7.95, c = 43.98}\text{A}\text{?}$ and β = 114.2°. The two molecules of the unit cell are related by a two-fold screw axis and pack tail-to-tail in a bilayer structure. The monosodium phosphate group packs with rather a small cross-section (24 Ų) relative to the two hydrocarbon chains. This unbalance in packing cross-section is overcome by an interdigitation of the phosphate head groups of adjacent bilayers and the formation of a single, common phosphate group layer at the bilayer interfaces. The phosphate groups are linked by hydrogen bonds to linear strands which laterally are separated by strands of sodium ions. The conformation of the molecules differs from that of other phospholipids. The glycerol chain is oriented parallel (instead of perpendicular) to the layer surface and the parallel stacking of the hydrocarbon chains is achieved by a bend of the γ-chain (instead of the β-chain). Otherwise the conformation of the glycerol dicarboxyl ester group displays the same preferred features as generally found in glycerophospholipids. The hydrocarbon chains pack according to the triclinic (T_∥) packing mode. The interaction and packing principles of the phosphate head group are discussed in relation to the structural behaviour of phosphatidic acid.     

Packing forces in nine crystal forms of cutinase

During the characterization of mutants and covalently inhibited complexes of Fusarium solani cutinase, nine different crystal forms have been obtained so far. Protein mutants with a different surface charge distribution form new intermolecular salt bridges or long-range electrostatic interactions that are accompanied by a change in the crystal packing. The whole protein surface is involved in the packing contacts and the hydrophobicities of the protein surfaces in mutual contact turned out to be noncorrelated, which indicates that the packing interactions are nonspecific. In the case of the hydrophobic variants, the packing contacts showed some specificity, as the protein in the crystal tends to form either crystallographic or noncrystallographic dimers, which shield the hydrophobic surface from the solvent. The likelihood of surface atoms to be involved in a crystal contact is the same for both polar and nonpolar atoms. However, when taking areas in the 200–600 Ų range, instead of individual atoms, the either highly hydrophobic or highly polar surface regions were found to have an increased probability of establishing crystal lattice contacts. The protein surface surrounding the active-site crevice of cutinase constitutes a large hydrophobic area that is involved in packing contacts in all the various crystalline contexts. Proteins 31:320–333, 1998. © 1998 Wiley-Liss, Inc.     

Generalized electrostatic model of the wrapping of DNA around oppositely charged proteins

Histonelike proteins in prokaryotes and histone octamers in eukaryotes carry large positive charges, which are responsible of strong electrostatic interactions with DNA. As a result, DNA wraps around proteins and genetic information is condensed. We describe a generalized model of these electrostatic interactions mediated by salt that explains the wrapping of DNA around the nucleosome octamer, around remodeling factors in eukaryotes and around histonelike proteins in prokaryotes. It comes out that small changes in protein dimension and charge produce large effects in the supramolecular DNA-protein architecture.     

Protein-nucleic acid recognition: statistical analysis of atomic interactions and influence of DNA structure

We analyzed structural features of 11,038 direct atomic contacts (either electrostatic, H-bonds, hydrophobic, or other van der Waals interactions) extracted from 139 protein-DNA and 49 protein-RNA nonhomologous complexes from the Protein Data Bank (PDB). Globally, H-bonds are the most frequent interactions (approximately 50%), followed by van der Waals, hydrophobic, and electrostatic interactions. From the protein viewpoint, hydrophilic amino acids are over-represented in the interaction databases: Positively charged amino acids mainly contact nucleic acid phosphate groups but can also interact with base edges. From the nucleotide point of view, DNA and RNA behave differently: Most protein-DNA interactions involve phosphate atoms, while protein-RNA interactions involve more frequently base edge and ribose atoms. The increased participation of DNA phosphate involves H-bonds rather than salt bridges. A statistical analysis was performed to find the occurrence of amino acid-nucleotide pairs most different from chance. These pairs were analyzed individually. Finally, we studied the conformation of DNA in the interaction sites. Despite the prevalence of B-DNA in the database, our results suggest that A-DNA is favored in the interaction sites.     

Interactions and space requirement of the phosphate head group of membrane lipids: The single crystal structures of a triclinic and a monoclinic form of hexadecyl-2-deoxyglycerophosphoric acid monohydrate

The crystal structures of a triclinic form (HPA1) and a monoclinic form (HPA2) of hexadecyl-2-deoxyglycerophosphoric acid monohydrate were determined by single crystal analysis. The unit cell dimensions for HPA1 are $\text{a = 4.75, b = 5.72, c = 44.36}\text{A}\text{?}$ and α = 91.0, β = 101.5, γ = 100.5° (P$\text{1}$) and for HPA2, $\text{a = 4.75, b = 5.72, c = 88.72}\text{A}\text{?}$ and γ = 100.8° (P2₁). In both structures the molecules are fully extended and pack tail-to-tail in bilayers with tilting (47°) hydrocarbon chains. In HPA2, however, the chain tilt alternatingly changes direction in adjacent bilayers, giving rise to a doubled unit cell which spans two bilayers. The dihydrogen phosphate groups interact by hydrogen bonds and are arranged in rows. Laterally between these phosphate rows the water molecules are accommodated producing a compact two-dimensional network of hydrogen bonds. The packing cross-section in the layer plane of the dihydrogen phosphate monohydrate group is 26.7 Ų in both structures. The hydrocarbon chains pack according to the triclinic (T|) chain packing mode. In HPA2, however, the chain packing is somewhat less compact with accounts for a 2% increase in the molecular volume. In both structures the ether oxygen is accommodated into the hydrocarbon matrix without distortion of the chain packing.     

The crystal structure of cholesteryl dihydrogen phosphate

Crystals of cholesteryl dihydrogen phosphate grown from 1,4-dioxane solution are monoclinic, space group C2 with $\text{a = 24.40, b = 6.27, c = 40.86}\text{A}\text{?}\text{and}\text{β = 102.7°}$. The asymmetric unit contains two molecules of cholesteryl phosphate CP and one dioxane molecule of the solvent. The CP molecules pack tail to tail in a bilayer structure. Within the layer they are arranged in double rows with their phosphate groups linked to ribbons by hydrogen bonds. Laterally the double strands of phosphate groups are separated by rows of dioxane molecules. The dioxane serves as hydrogen bond acceptor and as a spacer molecule that compensates the differences in cross-sectional area of the cholesteryl residue (38.4 Ų and the phosphate group (24 Ų). In the cholesterol matrix the CP molecules joined to double rows have packing contact with the smooth side of their skeleta and interdigitate with their annular methyl groups with those of molecules of the adjacent double rows. The branched cholesteryl side chains facing the bilayer center are loosely packed and show considerable disorder and/or thermal motion.     

Modification of protein crystal packing by systematic mutations of surface residues: Implications on biotemplating and crystal porosity

Bioinspired nano‐scale biotemplating for the development of novel composite materials has recently culminated in several demonstrations of nano‐structured hybrid materials. Protein crystals, routinely prepared for the elucidation of protein 3D structures by X‐ray crystallography, present an ordered and highly accurate 3D array of protein molecules. Inherent to the 3D arrangement of the protein “building blocks” in the crystal, a complementary 3D array of interconnected cavities—voids array, exhibiting highly ordered porosity is formed. The porous arrays of protein crystal may serve as a nano‐structured, accurate biotemplate by a “filling” process. These cavities arrays are shaped by the mode of protein packing throughout the crystallization process. Here we propose and demonstrate feasibility of targeting site specific mutations to modify protein's surface to affect protein crystal packing, enabling the generation of a series of protein crystal “biotemplates” all originating from same parent protein. The selection of these modification sites was based on in silico analysis of protein–protein interface contact areas in the parent crystal. The model protein selected for this study was the N‐terminal type II cohesin from the cellulosomal scaffold in ScaB subunit of Acetivibrio cellulolyticus and mutations were focused on lysine residues involved in protein packing as prime target. The impact of systematically mutating these lysine residues on protein packing and its resulting interconnected cavities array were found to be most significant when surface lysine residues were substituted to tryptophan residues. Our results demonstrate the feasibility of using pre‐designed site directed mutations for the generation of a series of protein crystal biotemplates from a “parent” protein. Biotechnol. Bioeng. 2009; 104: 444–457 © 2009 Wiley Periodicals, Inc.     

Charged residues at protein interaction interfaces: unexpected conservation and orchestrated divergence

Protein-protein interactions play an essential role in the functioning of cell. The importance of charged residues and their diverse role in protein-protein interactions have been well studied using experimental and computational methods. Often, charged residues located in protein interaction interfaces are conserved across the families of homologous proteins and protein complexes. However, on a large scale, it has been recently shown that charged residues are significantly less conserved than other residue types in protein interaction interfaces. The goal of this work is to understand the role of charged residues in the protein interaction interfaces through their conservation patterns. Here, we propose a simple approach where the structural conservation of the charged residue pairs is analyzed among the pairs of homologous binary complexes. Specifically, we determine a large set of homologous interactions using an interaction interface similarity measure and catalog the basic types of conservation patterns among the charged residue pairs. We find an unexpected conservation pattern, which we call the correlated reappearance, occurring among the pairs of homologous interfaces more frequently than the fully conserved pairs of charged residues. Furthermore, the analysis of the conservation patterns across different superkingdoms as well as structural classes of proteins has revealed that the correlated reappearance of charged residues is by far the most prevalent conservation pattern, often occurring more frequently than the unconserved charged residues. We discuss a possible role that the new conservation pattern may play in the long-range electrostatic steering effect.     

A protein-RNA docking benchmark (I): nonredundant cases

We have developed a nonredundant protein-RNA docking benchmark dataset, which is derived from the available bound and unbound structures in the Protein Data Bank involving polypeptide and nucleic acid chains. It consists of nine unbound-unbound cases where both the protein and the RNA are available in the free form. The other 36 cases are of unbound-bound type where only the protein is available in the free form. The conformational change upon complex formation is calculated by a distance matrix alignment method, and based on that, complexes are classified into rigid, semi-flexible, and full flexible. Although in the rigid body category, no significant conformational change accompanies complex formation, the fully flexible test cases show large domain movements, RNA base flips, etc. The benchmark covers four major groups of RNA, namely, t-RNA, ribosomal RNA, duplex RNA, and single-stranded RNA. We find that RNA is generally more flexible than the protein in the complexes, and the interface region is as flexible as the molecule as a whole. The structural diversity of the complexes in the benchmark set should provide a common ground for the development and comparison of the protein-RNA docking methods. The benchmark can be freely downloaded from the internet.     

Molecular packing and intermolecular interactions in two structural polymorphs of N-palmitoylethanolamine, a type 2 cannabinoid receptor agonist

The molecular structure, packing properties, and intermolecular interactions of two structural polymorphs of N-palmitoylethanolamine (NPEA) have been determined by single-crystal X-ray diffraction. Polymorphs alpha and beta crystallized in monoclinic space group P2(1)/c and orthorhombic space group Pbca, respectively. In both polymorphs, NPEA molecules are organized in a tail-to-tail manner, resembling a bilayer membrane. Although the molecular packing in polymorph alpha is similar to that in N-myristoylethanolamine and N-stearoylethanolamine, polymorph beta is a new form. The acyl chains in both polymorphs are tilted by approximately 35 degrees with respect to the bilayer normal, with their hydrocarbon moieties packed in an orthorhombic subcell. In both structures, the hydroxy group of NPEA forms two hydrogen bonds with the hydroxy groups of molecules in the opposite leaflet, resulting in extended, zig-zag type H-bonded networks along the b-axis in polymorph alpha and along the a-axis in polymorph beta. Additionally, the amide N-H and carbonyl groups of adjacent molecules are involved in N-H...O hydrogen bonds that connect adjacent molecules along the b-axis and a-axis, respectively, in alpha and beta. Whereas in polymorph alpha the L-shaped NPEA molecules in opposite layers are arranged to yield a Z-like organization, in polymorph beta one of the two NPEA molecules is rotated 180 degrees , leading to a W-like arrangement. Lattice energy calculations indicate that polymorph alpha is more stable than polymorph beta by approximately 2.65 kcal/mol.     

RBPBind: Quantitative Prediction of Protein-RNA Interactions

There are hundreds of RNA binding proteins in the human genome alone and their interactions with messenger and other RNAs in a cell regulate every step in an RNA's life cycle. To understand this interplay of proteins and RNA it is important to be able to know which protein binds which RNA how strongly and where. Here, we introduce RBPBind, a web-based tool for the quantitative prediction of the interaction of single-stranded RNA binding proteins with target RNAs that fully takes into account the effect of RNA secondary structure on binding affinity. Given a user-specified RNA and a protein selected from a set of several RNA-binding proteins, RBPBind computes their binding curve and effective binding constant. The server also computes the probability that, at a given protein concentration, a protein molecule will bind to any particular nucleotide along the RNA. The sequence specificity of the protein-RNA interaction is parameterized from public RNAcompete experiments and integrated into the recursions of the Vienna RNA package to simultaneously take into account protein binding and RNA secondary structure. We validate our approach by comparison to experimentally determined binding affinities of the HuR protein for several RNAs of different sequence contexts from the literature, showing that integration of raw sequence affinities into RNA secondary structure prediction significantly improves the agreement between computationally predicted and experimentally measured binding affinities. Our resource thus provides a quick and easy way to obtain reliable predicted binding affinities and locations for single-stranded RNA binding proteins based on RNA sequence alone.     

Kinase inhibitors and the case for CH...O hydrogen bonds in protein-ligand binding

Although the hydrogen bond is known to be an important mediator of intermolecular interactions, there has yet to be an analysis of the role of CH...O hydrogen bonds in protein-ligand complexes. In this work, we present evidence for such nonstandard hydrogen bonds from a survey of aromatic ligands in 184 kinase crystal structures and 358 high-resolution structures from the Protein Data Bank. CH groups adjacent to the positively charged nitrogen of nicotinamide exhibit geometric preferences strongly suggestive of hydrogen bonding interactions, as do heterocyclic CH groups in kinase ligands, while other aromatic CH groups do not exhibit these characteristics. Ab initio calculations reveal a considerable range of CH...O hydrogen bonding potentials among different aromatic ring systems, with nicotinamide and heterocycles preferred in kinase inhibitors showing particularly favorable interactions. These results provide compelling evidence for the existence of CH...O hydrogen bonds in protein-ligand interactions, as well as information on the relative strength of various aromatic CH donors. Such knowledge will be of considerable value in protein modeling, ligand design, and structure-activity analysis.     

Engineered variants of InlB with an additional leucine‐rich repeat discriminate between physiologically relevant and packing contacts in crystal structures of the InlB:MET complex

The physiological relevance of contacts in crystal lattices often remains elusive. This was also the case for the complex between the invasion protein internalin B (InlB) from Listeria monocytogenes and its host cell receptor, the human receptor tyrosine kinase (RTK) MET. InlB is a MET agonist and induces bacterial host cell invasion. Activation of RTKs generally involves ligand‐induced dimerization of the receptor ectodomain. The two currently available crystal structures of the InlB:MET complex show the same arrangement of InlB and MET in a 1:1 complex, but different dimeric 2:2 assemblies. Only one of these 2:2 assemblies is predicted to be stable by a computational procedure. This assembly is mainly stabilized by a contact between the Cap domain of InlB from one and the Sema domain of MET from another 1:1 complex. Here, we probe the physiological relevance of this interaction. We generated variants of the leucine‐rich repeat (LRR) protein InlB by inserting an additional repeat between the first and the second LRR. This should allow formation of the 1:1 complex but disrupt the potential 2:2 complex involving the Cap‐Sema contact due to steric distortions. A crystal structure of one of the engineered proteins showed that it folded properly. Binding affinity to MET was comparable to that of wild‐type InlB. The InlB variant induced MET phosphorylation and cell scatter like wild‐type InlB. These results suggest that the Cap‐Sema interaction is not physiologically relevant and support the previously proposed assembly, in which a 2:2 InlB:MET complex is built around a ligand dimer.