Codon usage is less optimized in eukaryotic gene segments encoding intrinsically disordered regions than in those encoding structural domains

Codon usage tends to be optimized in highly expressed genes. A plausible explanation for this phenomenon is that translational accuracy is increased in highly expressed genes with infrequent use of rare codons. Besides structural domains (SDs), eukaryotic proteins generally have intrinsically disordered regions (IDRs) that by themselves do not assume unique three-dimensional structures. As IDRs are free from structural constraint, they can probably accommodate more translational errors than SDs can. Thus, codon usage in IDRs is likely to be less optimized than that in SDs. Codon usage in all the genes of seven eukaryotes was examined in terms of both tRNA adaptation index and codon adaptation index. Different amino acid compositions in different protein regions were taken into account in calculating expected adaptation indices, to which observed indices were compared. Codon usage is less optimized in gene regions encoding IDRs than in those corresponding to SDs. The finding does not depend on whether IDRs are located at the N-terminus, in the middle, or at the C-terminus of proteins. Furthermore, the observation remains unchanged in two different algorithms used to predict IDRs in proteins. The result is consistent with the idea that IDRs tolerate more translational errors than SDs.     

Negatively Charged Disordered Regions are Prevalent and Functionally Important Across Proteomes

Intrinsically disordered regions (IDRs) of proteins are often characterized by a high fraction of charged residues, but differ in their overall net charge and in the organization of the charged residues. The function-encoding information stored via IDR charge composition and organization remains elusive. Here, we aim to decipher the sequence–function relationship in IDRs by presenting a comprehensive bioinformatic analysis of the charge properties of IDRs in the human, mouse, and yeast proteomes. About 50% of the proteins comprise at least a single IDR, which is either positively or negatively charged. Highly negatively charged IDRs are longer and possess greater net charge per residue compared with highly positively charged IDRs. A striking difference between positively and negatively charged IDRs is the characteristics of the repeated units, specifically, of consecutive Lys or Arg residues (K/R repeats) and Asp or Glu (D/E repeats) residues. D/E repeats are found to be about five times longer than K/R repeats, with the longest found containing 49 residues. Long stretches of consecutive D and E are found to be more prevalent in nucleic acid-related proteins. They are less common in prokaryotes, and in eukaryotes their abundance increases with genome size. The functional role of D/E repeats and the profound differences between them and K/R repeats are discussed.     

FG nucleoporins feature unique patterns that distinguish them from other IDPs

FG nucleoporins (FG Nups) are intrinsically disordered proteins and are the putative regulators of nucleocytoplasmic transport. They allow fast, yet selective, transport of molecules through the nuclear pore complex, but the underlying mechanism of nucleocytoplasmic transport is not yet fully discovered. As a result, FG Nups have been the subject of extensive research in the past two decades. Although most studies have been focused on analyzing the conformation and function of FG Nups from a biophysical standpoint, some recent studies have investigated the sequence-function relationship of FG Nups, with a few investigating amino acid sequences of a large number of FG Nups to understand common characteristics that might enable their function. Previously, we identified an evolutionarily conserved feature in FG Nup sequences, which are extended subsequences with low charge density, containing only positive charges, and located toward the N-terminus of FG Nups. We named these patterns longest positive like charge regions (lpLCRs). These patterns are specific to positively charged residues, and negatively charged residues do not demonstrate such a pattern. In this study, we compare FG Nups with other disordered proteins obtained from the DisProt and UniProt database in terms of presence of lpLCRs. Our results show that the lpLCRs are virtually exclusive to FG Nups and are not observed in other disordered proteins. Also, lpLCRs are what differentiate FG Nups from DisProt proteins in terms of charge distribution, meaning that excluding lpLCRs from the sequences of FG Nups make them similar to DisProt proteins in terms of charge distribution. We also previously showed the biophysical effect of lpLCRs in conformation of FG Nups. The results of this study are in line with our previous findings and imply that lpLCRs are virtually exclusive and functionally significant characteristics of FG Nups and nucleocytoplasmic transport.     

Exploring disorder in the human charged biased proteins

A considerable interest has been put in the identification of biased regions in proteins. These regions are frequently associated with a structural role in the cell and particularly with protein disorder. Here, we have investigated the intrinsically disordered regions (IDRs) in the human charged biased proteins identified in our earlier work. We found that 65% of charged biased proteins contained significant IDRs involved particularly in DNA and RNA binding. Also, we have observed that these proteins are well conserved in metazoans and more particularly in mammalian. In addition, the IDRs are located largely in N-terminal, C-terminal sequence flanking the functional domains (FD) and slightly less in (FD) itself. Our work also supports the association between protein disorder and protein–protein/DNA interaction. An example will be described.     

Intrinsically disordered regions in autophagy proteins

Autophagy is an essential eukaryotic pathway required for cellular homeostasis. Numerous key autophagy effectors and regulators have been identified, but the mechanism by which they carry out their function in autophagy is not fully understood. Our rigorous bioinformatic analysis shows that the majority of key human autophagy proteins include intrinsically disordered regions (IDRs), which are sequences lacking stable secondary and tertiary structure; suggesting that IDRs play an important, yet hitherto uninvestigated, role in autophagy. Available crystal structures corroborate the absence of structure in some of these predicted IDRs. Regions of orthologs equivalent to the IDRs predicted in the human autophagy proteins are poorly conserved, indicating that these regions may have diverse functions in different homologs. We also show that IDRs predicted in human proteins contain several regions predicted to facilitate protein–protein interactions, and delineate the network of proteins that interact with each predicted IDR‐containing autophagy protein, suggesting that many of these interactions may involve IDRs. Lastly, we experimentally show that a BCL2 homology 3 domain (BH3D), within the key autophagy effector BECN1 is an IDR. This BH3D undergoes a dramatic conformational change from coil to α‐helix upon binding to BCL2s, with the C‐terminal half of this BH3D constituting a binding motif, which serves to anchor the interaction of the BH3D to BCL2s. The information presented here will help inform future in‐depth investigations of the biological role and mechanism of IDRs in autophagy proteins. Proteins 2014; 82:565–578. © 2013 Wiley Periodicals, Inc.     

An archaebacterial homolog of pelota, a meiotic cell division protein in eukaryotes

Abstract An open reading frame ( pelA ) specifying a homolog of pelota and DOM34, proteins required for meiotic cell division in Drosophila melanogaster and Saccharomyces cerevisiae , respectively, has been cloned, sequenced and identified from the archaebacterium Sulfolobus solfataricus . The S. solfataricus PelA protein is about 20% identical with pelota, DOM34 and the hypothetical protein R74.6 of Caenorhabditis elegans . The presence of a pelota homolog in archaebacteria implies that the meiotic functions of the eukaryotic protein were co-opted from, or added to, other functions existing before the emergence of eukaryotes. The nuclear localization signal and negatively charged carboxy-terminus characteristic of eukaryotic pelota-like proteins are absent from the S. solfataricus homolog, and hence may be indicative of the acquired eukaryotic function(s).     

Charge Segregation and Low Hydrophobicity Are Key Features of Ribosomal Proteins from Different Organisms

Ribosomes are large and highly charged macromolecular complexes consisting of RNA and proteins. Here, we address the electrostatic and nonpolar properties of ribosomal proteins that are important for ribosome assembly and interaction with other cellular components and may influence protein folding on the ribosome. We examined 50 S ribosomal subunits from 10 species and found a clear distinction between the net charge of ribosomal proteins from halophilic and non-halophilic organisms. We found that ∼67% ribosomal proteins from halophiles are negatively charged, whereas only up to ∼15% of ribosomal proteins from non-halophiles share this property. Conversely, hydrophobicity tends to be lower for ribosomal proteins from halophiles than for the corresponding proteins from non-halophiles. Importantly, the surface electrostatic potential of ribosomal proteins from all organisms, especially halophiles, has distinct positive and negative regions across all the examined species. Positively and negatively charged residues of ribosomal proteins tend to be clustered in buried and solvent-exposed regions, respectively. Hence, the majority of ribosomal proteins is characterized by a significant degree of intramolecular charge segregation, regardless of the organism of origin. This key property enables the ribosome to accommodate proteins within its complex scaffold regardless of their overall net charge.     

Sequence-to-Conformation Relationships of Disordered Regions Tethered to Folded Domains of Proteins

Intrinsically disordered proteins and regions (IDPs/IDRs) are characterized by well-defined sequence-to-conformation relationships (SCRs). These relationships refer to the sequence-specific preferences for average sizes, shapes, residue-specific secondary structure propensities, and amplitudes of multiscale conformational fluctuations. SCRs are discerned from the sequence-specific conformational ensembles of IDPs. A vast majority of IDPs are actually tethered to folded domains (FDs). This raises the question of whether or not SCRs inferred for IDPs are applicable to IDRs tethered to FDs. Here, we use atomistic simulations based on a well-established forcefield paradigm and an enhanced sampling method to obtain comparative assessments of SCRs for 13 archetypal IDRs modeled as autonomous units, as C-terminal tails connected to FDs, and as linkers between pairs of FDs. Our studies uncover a set of general observations regarding context-independent versus context-dependent SCRs of IDRs. SCRs are minimally perturbed upon tethering to FDs if the IDRs are deficient in charged residues and for polyampholytic IDRs where the oppositely charged residues within the sequence of the IDR are separated into distinct blocks. In contrast, the interplay between IDRs and tethered FDs has a significant modulatory effect on SCRs if the IDRs have intermediate fractions of charged residues or if they have sequence-intrinsic conformational preferences for canonical random coils. Our findings suggest that IDRs with context-independent SCRs might be independent evolutionary modules, whereas IDRs with context-dependent SCRs might co-evolve with the FDs to which they are tethered.     

Benchmarks for flexible and rigid transcription factor-DNA docking

Background

Although structural domains in proteins (SDs) are important, half of the regions in the human proteome are currently left with no SD assignments. These unassigned regions consist not only of novel SDs, but also of intrinsically disordered (ID) regions since proteins, especially those in eukaryotes, generally contain a significant fraction of ID regions. As ID regions can be inferred from amino acid sequences, a method that combines SD and ID region assignments can determine the fractions of SDs and ID regions in any proteome.

Results

In contrast to other available ID prediction programs that merely identify likely ID regions, the DICHOT system we previously developed classifies the entire protein sequence into SDs and ID regions. Application of DICHOT to the human proteome revealed that residue-wise ID regions constitute 35%, SDs with similarity to PDB structures comprise 52%, while SDs with no similarity to PDB structures account for the remaining 13%. The last group consists of novel structural domains, termed cryptic domains, which serve as good targets of structural genomics. The DICHOT method applied to the proteomes of other model organisms indicated that eukaryotes generally have high ID contents, while prokaryotes do not. In human proteins, ID contents differ among subcellular localizations: nuclear proteins had the highest residue-wise ID fraction (47%), while mitochondrial proteins exhibited the lowest (13%). Phosphorylation and O-linked glycosylation sites were found to be located preferentially in ID regions. As O-linked glycans are attached to residues in the extracellular regions of proteins, the modification is likely to protect the ID regions from proteolytic cleavage in the extracellular environment. Alternative splicing events tend to occur more frequently in ID regions. We interpret this as evidence that natural selection is operating at the protein level in alternative splicing.

Conclusions

We classified entire regions of proteins into the two categories, SDs and ID regions and thereby obtained various kinds of complete genome-wide statistics. The results of the present study are important basic information for understanding protein structural architectures and have been made publicly available at http://spock.genes.nig.ac.jp/~genome/DICHOT.     

A hydrophilic microenvironment in the substrate-translocating groove of the YidC membrane insertase is essential for enzyme function

The YidC family of proteins are membrane insertases that catalyze the translocation of the periplasmic domain of membrane proteins via a hydrophilic groove located within the inner leaflet of the membrane. All homologs have a strictly conserved, positively charged residue in the center of this groove. In Bacillus subtilis, the positively charged residue has been proposed to be essential for interacting with negatively charged residues of the substrate, supporting a hypothesis that YidC catalyzes insertion via an early-step electrostatic attraction mechanism. Here, we provide data suggesting that the positively charged residue is important not for its charge but for increasing the hydrophilicity of the groove. We found that the positively charged residue is dispensable for Escherichia coli YidC function when an adjacent residue at position 517 was hydrophilic or aromatic, but was essential when the adjacent residue was apolar. Additionally, solvent accessibility studies support the idea that the conserved positively charged residue functions to keep the top and middle of the groove sufficiently hydrated. Moreover, we demonstrate that both the E. coli and Streptococcus mutans YidC homologs are functional when the strictly conserved arginine is replaced with a negatively charged residue, provided proper stabilization from neighboring residues. These combined results show that the positively charged residue functions to maintain a hydrophilic microenvironment in the groove necessary for the insertase activity, rather than to form electrostatic interactions with the substrates.     

Dependence of prePhoA-phospholipid interaction in vivo and in vitro on charge of signal peptide N-terminus and content of anionic phospholipids in membranes

Replacement of the positively charged signal peptide with neutral or negatively charged peptides due to substitution of Lys(–20) in the N-terminal region of the signal peptide leads to decreases in the rate of prePhoA membrane translocation in vivo and in the efficiency of prePhoA insertion into liposomes in vitro. The effect of anionic phospholipids on prePhoA insertion into model membranes is determined by the signal peptide N-terminus charge, while the dependence of prePhoA translocation across the cytoplasmic membrane in vivo is not, under the studied variations in the content of anionic phospholipids. This is evidence of the possibility of direct electrostatic interaction between the signal peptide N-terminus and anionic phospholipids, which in vivo, however, seems to involve some proteins of the Sec machinery.     

Anionic lipids stimulate Sec-independent insertion of a membrane protein lacking charged amino acid side chains

We have investigated the influence of the different lipid classes of Escherichia coli on Sec-independent membrane protein insertion, using an assay in which a mutant of the single-spanning Pf3 coat protein is biosynthetically inserted into liposomes. It was found that phosphatidylethanolamine and other non-bilayer lipids do not have a significant effect on insertion. Surprisingly, the anionic lipids phosphatidylglycerol and cardiolipin stimulate N-terminal translocation of the protein, even though it has no charged amino acid side chains. This novel effect is general for anionic lipids and depends on the amount of charge on the lipid headgroup. Since the N-terminus of the protein is at least partially positively charged due to a helix dipole moment, apparently negatively charged lipids can stimulate translocation of slightly positively charged protein segments in a direction opposite to the positive-inside rule. A mechanism is proposed to explain these results.     

Specific sequence of motifs of mitochondrial uncoupling proteins

We have searched for the exclusivity of common sequence motifs of the mitochondrial uncoupling proteins (UCP1, UCP2, UCP3, UCP4, BMCP1, and plant UCP [PUMP]) within the gene family of mitochondrial anion carrier proteins. The UCP-specific sequences, "UCP signatures", were found in the first, second, and fourth alpha-helices. First: Ala/Ser-Cys/Thr/n-n/Phe-Ala/Gly-[negatively charged residue]-n/Phe-n/Cys-Thr-Phe/n; second: Gly/Ala-Ile/Leu-Gln/X-[positively charged residue]-NH-n/Cys-Ser/nphi/X-n/Ser-OH/Gly-n-[positively charged residue]-Ile/Met-Gly/Val-n/Thr; fourth: Pro-Asn/ Thr-n-X-[positively charged residue]-Asn/Ser/Ala-n-n-Ile/Leu-n-Asn/Val-Cys/n-n/Thr-[negatively charged residue]-n-n/Thr/Pro-OH/Val (n, nonpolar; phi, aromatic; (positively charged residue/negatively charged residue, charged residue). The second and part of the third signature are also present in the yeast dicarboxylate transporter. The UCP signature excluding BMCP1 was also found in the second matrix segment: [positively charged residue]-(Pro/ del-Leu/del)-[positively charged residue]-phi-X-Gly/Ser-Thr/n-X-NH/[negatively charged residue]-Ala-phi. These UCP signatures are thought to be involved in fatty acid anion binding and translocation.     

Expression, purification and the 1.8 angstroms resolution crystal structure of human neuron specific enolase

Human neuron-specific enolase (NSE) or isozyme gamma has been expressed with a C-terminal His-tag in Escherichia coli. The enzyme has been purified, crystallized and its crystal structure determined. In the crystals the enzyme forms the asymmetric complex NSE x Mg2 x SO4/NSE x Mg x Cl, where "/" separates the dimer subunits. The subunit that contains the sulfate (or phosphate) ion and two magnesium ions is in the closed conformation observed in enolase complexes with the substrate or its analogues; the other subunit is in the open conformation observed in enolase subunits without bound substrate or analogues. This indicates negative cooperativity for ligand binding between subunits. Electrostatic charge differences between isozymes alpha and gamma, -19 at physiological pH, are concentrated in the regions of the molecular surface that are negatively charged in alpha, i.e. surface areas negatively charged in alpha are more negatively charged in gamma, while areas that are neutral or positively charged tend to be charge-conserved.