The role of RNA sequence and structure in RNA--protein interactions

We investigate the sequence and structural properties of RNA-protein interaction sites in 211 RNA-protein chain pairs, the largest set of RNA-protein complexes analyzed to date. Statistical analysis confirms and extends earlier analyses made on smaller data sets. There are 24.6% of hydrogen bonds between RNA and protein that are nucleobase specific, indicating the importance of both nucleobase-specific and -nonspecific interactions. While there is no significant difference between RNA base frequencies in protein-binding and non-binding regions, distinct preferences for RNA bases, RNA structural states, protein residues, and protein secondary structure emerge when nucleobase-specific and -nonspecific interactions are considered separately. Guanine nucleobase and unpaired RNA structural states are significantly preferred in nucleobase-specific interactions; however, nonspecific interactions disfavor guanine, while still favoring unpaired RNA structural states. The opposite preferences of nucleobase-specific and -nonspecific interactions for guanine may explain discrepancies between earlier studies with regard to base preferences in RNA-protein interaction regions. Preferences for amino acid residues differ significantly between nucleobase-specific and -nonspecific interactions, with nonspecific interactions showing the expected bias towards positively charged residues. Irregular protein structures are strongly favored in interactions with the protein backbone, whereas there is little preference for specific protein secondary structure in either nucleobase-specific interaction or -nonspecific interaction. Overall, this study shows strong preferences for both RNA bases and RNA structural states in protein-RNA interactions, indicating their mutual importance in protein recognition.     

A new classification of membrane protein crystals

Although being much smaller than the number of soluble proteins in the Protein Data Bank, the number of membrane proteins therein now approaches 700, and a statistical analysis becomes meaningful. Such an analysis showed that the conventional subdivision into monotopic, β-barrel and α-helical membrane proteins is appropriate but should be amended by a classification according to the detergent micelle structure in the crystal, which can be derived from the packing of the membrane-immersed parts of the proteins. The crystal packing density is specific for the three conventional types of membrane proteins and soluble proteins. It is also specific for three observed detergent arrangements that are micelle pockets, micelle filaments and micelle sheets, demonstrating that the detergent structure affects crystallization. The packing density distribution of crystals from integral membrane proteins has approximately the same shape as that of soluble proteins but is by a factor of two broader and shifted to lower density. It seems unlikely that the differences can be explained by a mere solvent expansion due to the required detergent. The crystallized membrane proteins were further analyzed with respect to protein mass, oligomerization and crystallographic asymmetric unit, space group, crystal ordering and symmetry. The results provide a new view on membrane proteins.     

Data Deposition and Annotation at the Worldwide Protein Data Bank

The Protein Data Bank (PDB) is the repository for three-dimensional structures of biological macromolecules, determined by experimental methods. The data in the archive is free and easily available via the Internet from any of the worldwide centers managing this global archive. These data are used by scientists, researchers, bioinformatics specialists, educators, students, and general audiences to understand biological phenomenon at a molecular level. Analysis of this structural data also inspires and facilitates new discoveries in science. This chapter describes the tools and methods currently used for deposition, processing, and release of data in the PDB. References to future enhancements are also included. Shuchismita Dutta, Kyle Burkhardt, and Ganesh J. Swaminathan have contributed equally to this work.     

RCSB Protein Data bank: Tools for visualizing and understanding biological macromolecules in 3D

Now in its 52nd year of continuous operations, the Protein Data Bank (PDB) is the premiere open‐access global archive housing three‐dimensional (3D) biomolecular structure data. It is jointly managed by the Worldwide Protein Data Bank (wwPDB) partnership. The Research Collaboratory for Structural Bioinformatics Protein Data Bank (RCSB PDB) is funded by the National Science Foundation, National Institutes of Health, and US Department of Energy and serves as the US data center for the wwPDB. RCSB PDB is also responsible for the security of PDB data in its role as wwPDB‐designated Archive Keeper. Every year, RCSB PDB serves tens of thousands of depositors of 3D macromolecular structure data (coming from macromolecular crystallography, nuclear magnetic resonance spectroscopy, electron microscopy, and micro‐electron diffraction). The RCSB PDB research‐focused web portal (RCSB.org) makes PDB data available at no charge and without usage restrictions to many millions of PDB data consumers around the world. The RCSB PDB training, outreach, and education web portal (PDB101.RCSB.org) serves nearly 700 K educators, students, and members of the public worldwide. This invited Tools Issue contribution describes how RCSB PDB (i) is organized; (ii) works with wwPDB partners to process new depositions; (iii) serves as the wwPDB‐designated Archive Keeper; (iv) enables exploration and 3D visualization of PDB data via RCSB.org; and (v) supports training, outreach, and education via PDB101.RCSB.org. New tools and features at RCSB.org are presented using examples drawn from high‐resolution structural studies of proteins relevant to treatment of human cancers by targeting immune checkpoints.     

Dynamics may significantly influence the estimation of interatomic distances in biomolecular X-ray structures

Atomic positions obtained by X-ray crystallography are time and space averages over many molecules in the crystal. Importantly, interatomic distances, calculated between such average positions and frequently used in structural and mechanistic analyses, can be substantially different from the more appropriate time-average and ensemble-average interatomic distances. Using crystallographic B-factors, one can deduce corrections, which have so far been applied exclusively to small molecules, to obtain correct average distances as a function of the type of atomic motion. Here, using 4774 high-quality protein X-ray structures, we study the significance of such corrections for different types of atomic motion. Importantly, we show that for distances shorter than 5 Å, corrections greater than 0.5 Å may apply, especially for noncorrelated or anticorrelated motion. For example, 14% of the studied structures have at least one pair of atoms with a correction of ≥ 0.5 Å in the case of noncorrelated motion. Using molecular dynamics simulations of villin headpiece, ubiquitin, and SH3 domain unit cells, we demonstrate that the majority of average interatomic distances in these proteins agree with noncorrelated corrections, suggesting that such deviations may be truly relevant. Importantly, we demonstrate that the corrections do not significantly affect stereochemistry and the overall quality of final refined X-ray structures, but can provide marked improvements in starting unrefined models obtained from low-resolution X-ray data. Finally, we illustrate the potential mechanistic and biological significance of the calculated corrections for KcsA ion channel and show that they provide indirect evidence that motions in its selectivity filter are highly correlated.     

Crystal structures of cobalamin-independent methionine synthase (MetE) from Streptococcus mutans: a dynamic zinc-inversion model

Cobalamin-independent methionine synthase (MetE) catalyzes the direct transfer of a methyl group from methyltetrahydrofolate to l-homocysteine to form methionine. Previous studies have shown that the MetE active site coordinates a zinc atom, which is thought to act as a Lewis acid and plays a role in the activation of thiol. Extended X-ray absorption fine structure studies and mutagenesis experiments identified the zinc-binding site in MetE from Escherichia coli. Further structural investigations of MetE from Thermotoga maritima lead to the proposition of two models: “induced fit” and “dynamic equilibrium”, to account for the catalytic mechanisms of MetE. Here, we present crystal structures of oxidized and zinc-replete MetE from Streptococcus mutans at the physiological pH. The structures reveal that zinc is mobile in the active center and has the possibility to invert even in the absence of homocysteine. These structures provide evidence for the dynamic equilibrium model.     

Topological classification and enumeration of RNA structures by genus

To an RNA pseudoknot structure is naturally associated a topological surface, which has its associated genus, and structures can thus be classified by the genus. Based on earlier work of Harer–Zagier, we compute the generating function $\mathbf{D}_{g,\sigma }(z)=\sum _{n}\mathbf{d}_{g,\sigma }(n)z^n$ for the number $\mathbf{d}_{g,\sigma }(n)$ of those structures of fixed genus $g$ and minimum stack size $\sigma $ with $n$ nucleotides so that no two consecutive nucleotides are basepaired and show that $\mathbf{D}_{g,\sigma }(z)$ is algebraic. In particular, we prove that $\mathbf{d}_{g,2}(n)\sim k_g\,n^{3(g-\frac{1}{2})} \gamma _2^n$ , where $\gamma _2\approx 1.9685$ . Thus, for stack size at least two, the genus only enters through the sub-exponential factor, and the slow growth rate compared to the number of RNA molecules implies the existence of neutral networks of distinct molecules with the same structure of any genus. Certain RNA structures called shapes are shown to be in natural one-to-one correspondence with the cells in the Penner–Strebel decomposition of Riemann’s moduli space of a surface of genus $g$ with one boundary component, thus providing a link between RNA enumerative problems and the geometry of Riemann’s moduli space.     

The capsid of the small RNA phage PRR1 is stabilized by metal ions

Many nonenveloped virus particles are stabilized by calcium ions bound in the interfaces between the protein subunits. These ions may have a role in the disassembly process. The small RNA phages of the Leviviridae family have T = 3 quasi-symmetry and are unique among simple viruses in that they have a coat protein with a translational repressor activity and a fold that has not been observed in other viruses. The crystal structure of phage PRR1 has been determined to 3.5 Å resolution. The structure shows a tentative binding site for a calcium ion close to the quasi-3-fold axis. The RNA-binding surface used for repressor activity is mostly conserved. The structure does not show any significant differences between quasi-equivalent subunits, which suggests that the assembly is not controlled by conformational switches as in many other simple viruses.     

Analysis of Four-Way Junctions in RNA Structures

RNA secondary structures can be divided into helical regions composed of canonical Watson-Crick and related base pairs, as well as single-stranded regions such as hairpin loops, internal loops, and junctions. These elements function as building blocks in the design of diverse RNA molecules with various fundamental functions in the cell. To better understand the intricate architecture of three-dimensional (3D) RNAs, we analyze existing RNA four-way junctions in terms of base-pair interactions and 3D configurations. Specifically, we identify nine broad junction families according to coaxial stacking patterns and helical configurations. We find that helices within junctions tend to arrange in roughly parallel and perpendicular patterns and stabilize their conformations using common tertiary motifs such as coaxial stacking, loop-helix interaction, and helix packing interaction. Our analysis also reveals a number of highly conserved base-pair interaction patterns and novel tertiary motifs such as A-minor-coaxial stacking combinations and sarcin/ricin motif variants. Such analyses of RNA building blocks can ultimately help in the difficult task of RNA 3D structure prediction.     

Tertiary Motifs Revealed in Analyses of Higher-Order RNA Junctions

RNA junctions are secondary-structure elements formed when three or more helices come together. They are present in diverse RNA molecules with various fundamental functions in the cell. To better understand the intricate architecture of three-dimensional (3D) RNAs, we analyze currently solved 3D RNA junctions in terms of base-pair interactions and 3D configurations. First, we study base-pair interaction diagrams for solved RNA junctions with 5 to 10 helices and discuss common features. Second, we compare these higher-order junctions to those containing 3 or 4 helices and identify global motif patterns such as coaxial stacking and parallel and perpendicular helical configurations. These analyses show that higher-order junctions organize their helical components in parallel and helical configurations similar to lower-order junctions. Their sub-junctions also resemble local helical configurations found in three- and four-way junctions and are stabilized by similar long-range interaction preferences such as A-minor interactions. Furthermore, loop regions within junctions are high in adenine but low in cytosine, and in agreement with previous studies, we suggest that coaxial stacking between helices likely forms when the common single-stranded loop is small in size; however, other factors such as stacking interactions involving noncanonical base pairs and proteins can greatly determine or disrupt coaxial stacking. Finally, we introduce the ribo-base interactions: when combined with the along-groove packing motif, these ribo-base interactions form novel motifs involved in perpendicular helix-helix interactions. Overall, these analyses suggest recurrent tertiary motifs that stabilize junction architecture, pack helices, and help form helical configurations that occur as sub-elements of larger junction networks. The frequent occurrence of similar helical motifs suggest nature's finite and perhaps limited repertoire of RNA helical conformation preferences. More generally, studies of RNA junctions and tertiary building blocks can ultimately help in the difficult task of RNA 3D structure prediction.     

Assessing side-chain perturbations of the protein backbone: a knowledge-based classification of residue Ramachandran space

Grouping the 20 residues is a classic strategy to discover ordered patterns and insights about the fundamental nature of proteins, their structure, and how they fold. Usually, this categorization is based on the biophysical and/or structural properties of a residue's side-chain group. We extend this approach to understand the effects of side chains on backbone conformation and to perform a knowledge-based classification of amino acids by comparing their backbone phi, psi distributions in different types of secondary structure. At this finer, more specific resolution, torsion angle data are often sparse and discontinuous (especially for nonhelical classes) even though a comprehensive set of protein structures is used. To ensure the precision of Ramachandran plot comparisons, we applied a rigorous Bayesian density estimation method that produces continuous estimates of the backbone phi, psi distributions. Based on this statistical modeling, a robust hierarchical clustering was performed using a divergence score to measure the similarity between plots. There were seven general groups based on the clusters from the complete Ramachandran data: nonpolar/beta-branched (Ile and Val), AsX (Asn and Asp), long (Met, Gln, Arg, Glu, Lys, and Leu), aromatic (Phe, Tyr, His, and Cys), small (Ala and Ser), bulky (Thr and Trp), and, lastly, the singletons of Gly and Pro. At the level of secondary structure (helix, sheet, turn, and coil), these groups remain somewhat consistent, although there are a few significant variations. Besides the expected uniqueness of the Gly and Pro distributions, the nonpolar/beta-branched and AsX clusters were very consistent across all types of secondary structure. Effectively, this consistency across the secondary structure classes implies that side-chain steric effects strongly influence a residue's backbone torsion angle conformation. These results help to explain the plasticity of amino acid substitutions on protein structure and should help in protein design and structure evaluation.     

A hydrogen bond regulates slow motions in ubiquitin by modulating a β-turn flip

Proteins exist as conformational ensembles composed of multiple interchanging substates separated by kinetic barriers. Interconverting conformations are often difficult to probe, owing to their sparse population and transient nature. Here, we report the identification and characterization of a subset of conformations in ubiquitin that participate in microsecond-to-millisecond motions in the amides of Ile23, Asn25, and Thr55. A novel side chain to the backbone hydrogen bond that regulates these motions has also been identified. Combining our NMR studies with the available X-ray data, we have unearthed the physical process underlying slow motions—the interconversion of a type I into a type II β-turn flip at residues Glu51 through Arg54. Interestingly, the dominant conformer of wild-type ubiquitin observed in solution near neutral pH is only represented by about 22% of the crystal structures. The conformers generated as a result of the dynamics of the hydrogen bond appear to be correlated to ligand recognition by ubiquitin.     

Structural Analysis of Peroxide-Soaked MnSOD Crystals Reveals Side-On Binding of Peroxide to Active-Site Manganese

The superoxide dismutase (SOD) enzymes are important antioxidant agents that protect cells from reactive oxygen species. The SOD family is responsible for catalyzing the disproportionation of superoxide radical to oxygen and hydrogen peroxide. Manganese- and iron-containing SOD exhibit product inhibition whereas Cu/ZnSOD does not. Here, we report the crystal structure of Escherichia coli MnSOD with hydrogen peroxide cryotrapped in the active site. Crystallographic refinement to 1.55 Å and close inspection revealed electron density for hydrogen peroxide in three of the four active sites in the asymmetric unit. The hydrogen peroxide molecules are in the position opposite His26 that is normally assumed by water in the trigonal bipyramidal resting state of the enzyme. Hydrogen peroxide is present in active sites B, C, and D and is side-on coordinated to the active-site manganese. In chains B and D, the peroxide is oriented in the plane formed by manganese and ligands Asp167 and His26. In chain C, the peroxide is bound, making a 70° angle to the plane. Comparison of the peroxide-bound active site with the hydroxide-bound octahedral form shows a shifting of residue Tyr34 towards the active site when peroxide is bound. Comparison with peroxide-soaked Cu/ZnSOD indicates end-on binding of peroxide when the SOD does not exhibit inhibition by peroxide and side-on binding of peroxide in the product-inhibited state of MnSOD.     

History of Protein Data Bank Japan: standing at the beginning of the age of structural genomics

Prof. Haruki Nakamura, who is the former head of Protein Data Bank Japan (PDBj) and an expert in computational biology, retired from Osaka University at the end of March 2018. He founded PDBj at the Institute for Protein Research, together with other faculty members, researchers, engineers, and annotators in 2000, and subsequently established the worldwide Protein Data Bank (wwPDB) in 2003 to manage the core archive of the Protein Data Bank (PDB), collaborating with RCSB-PDB in the USA and PDBe in Europe. As the former head of PDBj and also an expert in structural bioinformatics, he has grown PDBj to become a well-known data center within the structural biology community and developed several related databases, tools and integrated with new technologies, such as the semantic web, as primary services offered by PDBj.     

Mapping of protein-protein interaction sites by the 'absence of interference' approach

Protein-protein interactions are critical to most biological processes, and locating protein-protein interfaces on protein structures is an important task in molecular biology. We developed a new experimental strategy called the ‘absence of interference’ approach to determine surface residues involved in protein-protein interaction of established yeast two-hybrid pairs of interacting proteins. One of the proteins is subjected to high-level randomization by error-prone PCR. The resulting library is selected by yeast two-hybrid system for interacting clones that are isolated and sequenced. The interaction region can be identified by an absence or depletion of mutations. For data analysis and presentation, we developed a Web interface that analyzes the mutational spectrum and displays the mutational frequency on the surface of the structure (or a structural model) of the randomized protein†. Additionally, this interface might be of use for the display of mutational distributions determined by other types of random mutagenesis experiments. We applied the approach to map the interface of the catalytic domain of the DNA methyltransferase Dnmt3a with its regulatory factor Dnmt3L. Dnmt3a was randomized with high mutational load. A total of 76 interacting clones were isolated and sequenced, and 648 mutations were identified. The mutational pattern allowed to identify a unique interaction region on the surface of Dnmt3a, which comprises about 500-600 Å². The results were confirmed by site-directed mutagenesis and structural analysis. The absence-of-interference approach will allow high-throughput mapping of protein interaction sites suitable for functional studies and protein docking.     

On the Structure of the Proton-Binding Site in the Fo Rotor of Chloroplast ATP Synthases

The recently reported crystal structures of the membrane-embedded proton-dependent c-ring rotors of a cyanobacterial F₁F_o ATP synthase and a chloroplast F₁F_o ATP synthase have provided new insights into the mechanism of this essential enzyme. While the overall features of these c-rings are similar, a discrepancy in the structure and hydrogen-bonding interaction network of the H⁺ sites suggests two distinct binding modes, potentially reflecting a mechanistic differentiation. Importantly, the conformation of the key glutamate side chain to which the proton binds is also altered. To investigate the nature of these differences, we use molecular dynamics simulations of both c-rings embedded in a phospholipid membrane. We observe that the structure of the c₁₅ ring from Spirulina platensis is unequivocally stable within the simulation time. By contrast, the proposed structure of the H⁺ site in the chloroplast c₁₄ ring changes rapidly and consistently into that reported for the c₁₅ ring, indicating that the latter represents a common binding mode. To assess this hypothesis, we have remodeled the c₁₄ ring by molecular replacement using the published structure factors. The resulting structure provides clear evidence in support of a common binding site conformation and is also considerably improved statistically. These findings, taken together with a sequence analysis of c-subunits in the ATP synthase family, indicate that the so-called proton-locked conformation observed in the c₁₅ ring may be a common characteristic not only of light-driven systems such as chloroplasts and cyanobacteria but also of a selection of other bacterial species.     

Enhanced validation of small-molecule ligands and carbohydrates in the Protein Data Bank

1. Download : Download high-res image (176KB)
2. Download : Download full-size image

     

Crystal structure of scytalidoglutamic peptidase with its first potent inhibitor provides insights into substrate specificity and catalysis

Scytalidoglutamic peptidase (SGP) from Scytalidium lignicolum is the founding member of the newly discovered\ family of peptidases, G1, so far found exclusively in fungi. The crystal structure of SGP revealed a previously undescribed fold for peptidases and a unique catalytic dyad of residues Gln53 and Glu136. Surprisingly, the beta-sandwich structure of SGP is strikingly similar to members of the carbohydrate-binding concanavalin A-like lectins/glucanases superfamily. By analogy with the active sites of aspartic peptidases, a mechanism employing nucleophillic attack by a water molecule activated by the general base functionality of Glu136 has been proposed. Here, we report the first crystal structures of SGP in complex with two transition state peptide analogs designed to mimic the tetrahedral intermediate of the proteolytic reaction. Of these two analogs, the one containing a central S-hydroxyl group is a potent sub-nanomolar inhibitor of SGP. The inhibitor binds non-covalently to the concave surface of the upper beta-sheet and enables delineation of the S4 to S3' substrate specificity pockets of the enzyme. Structural differences in these pockets account for the unique substrate preferences of SGP among peptidases having an acidic pH optimum. Inhibitor binding is accompanied by a structuring of the region comprising residues Tyr71-Gly80 from being mostly disordered in the apoenzyme and leading to positioning of crucial active site residues for establishing enzyme-inhibitor contacts. In addition, conformational rearrangements are seen in a disulfide bridged surface loop (Cys141-Cys148), which moves inwards, partially closing the open substrate binding cleft of the native enzyme. The non-hydrolysable scissile bond analog of the inhibitor is located in the active site forming close contacts with Gln53 and Glu136. The nucleophilic water molecule is displaced and a unique mode of binding is observed with the S-OH of the inhibitor occupying the oxyanion binding site of the proposed tetrahedral intermediate. Details of the enzyme-inhibitor interactions and mechanistic interpretations are discussed.     

High-resolution crystal structures of trimeric and rod phycocyanin

The phycobilisome light-harvesting antenna in cyanobacteria and red algae is assembled from two substructures: a central core composed of allophycocyanin surrounded by rods that always contain phycocyanin (PC). Unpigmented proteins called linkers are also found within the rods and core. We present here two new structures of PC from the thermophilic cyanobacterium Thermosynechococcus vulcanus. We have determined the structure of trimeric PC to 1.35 Å, the highest resolution reported to date for this protein. We also present a structure of PC isolated in its intact and functional rod form at 1.5 Å. Analysis of rod crystals showed that in addition to the α and β PC subunit, there were three linker proteins: the capping rod linker (L_R^8.7), the rod linker (L_R), and only one of three rod-core linkers (L_RC, CpcG4) with a stoichiometry of 12:12:1:1:1. This ratio indicates that the crystals contained rods composed of two hexamers. The crystallographic parameters of the rod crystals are nearly identical with that of the trimeric form, indicating that the linkers do not affect crystal packing and are completely embedded within the rod cavities. Absorption and fluorescence emission spectra were red-shifted, as expected for assembled rods, and this could be shown for the rod in solution as well as in crystal using confocal fluorescence microscopy. The crystal packing imparts superimposition of the three rod linkers, canceling out their electron density. However, analysis of B-factors and the conformations of residues facing the rod channel indicate the presence of linkers. Based on the experimental evidence presented here and a homology-based model of the L_R protein, we suggest that the linkers do not in fact link between rod hexamers but stabilize the hexameric assembly and modify rod energy absorption and transfer capabilities.