Viral disorder or disordered viruses: do viral proteins possess unique features?

Many proteins or their regions are disordered in their native, biologically active states. Bioinformatics has revealed that these proteins/regions are highly abundant in different proteomes and carry out mostly regulatory functions related to molecular recognition, signal transduction, protein-protein, and protein-nucleic acid interactions. Viruses, these "organisms at the edge of life", have uniquely evolved to be highly adaptive for fast change in their biological and physical environment. To sustain these fast environmental changes, viral proteins elaborated multiple measures, from relatively low van der Waals contact densities, to inclusion of a large fraction of residues that are not arranged in well-defined secondary structural elements, to heavy use of short disordered regions, and to high resistance to mutations. On the other hand, viral proteins are rich in intrinsic disorder. Some of the intrinsically disordered regions are heavily used in the functioning of viral proteins. Others likely have evolved to help viruses accommodate to their hostile habitats. Still others evolved to help viruses in managing their economic usage of genetic material via alternative splicing, overlapping genes, and anti-sense transcription. In this review, we focus on structural peculiarities of viral proteins and on the role of intrinsic disorder in their functions.

Clusterin, a binding protein with a molten globule-like region.

Clusterin is a heterodimeric glycoprotein found in many tissues of the body and is the most abundant protein secreted by cultured rat Sertoli cells. The function of clusterin is unknown, but it has been associated with cellular injury, lipid transport, apoptosis, and it may be involved in the clearance of cellular debris caused by cell injury or death. Consistent with this last idea, clusterin has been shown to bind to a variety of molecules with high affinity including lipids, peptides, and proteins and the hydrophobic probe 1-anilino-8-naphthalenesulfonate (ANS). Given this variety of ligands, clusterin must have specific structural features that provide the protein with its promiscuous binding activity. Using sequence analyses, we show that clusterin likely contains three long regions of natively disordered or molten globule-like structures containing putative amphipathic alpha-helices. These disordered regions were highly sensitive to trypsin digestion, indicating a flexible nature. The effects of denaturation on the fluorescence of the clusterin-ANS complex were compared between proteins with structured binding pockets and molten globular forms of proteins. Clusterin bound ANS in a manner that was very similar to that of molten globular proteins. Furthermore, we found that, when bound to ANS, at least one cleavage site within the protease-sensitive disordered regions of clusterin was protected from trypsin digestion. In addition, we show that clusterin can function as a biological detergent that can solubilize bacteriorhodopsin. We propose that natively disordered regions with amphipathic helices form a dynamic, molten globule-like binding site and provide clusterin the ability to bind to a variety of molecules.

Beta-lactoglobulin molten globule induced by high pressure.

Beta-lactoglobulin (beta-LG) was treated with high hydrostatic pressure (HHP) at 600 MPa and 50 degrees C for selected times as long as 64 min. The intrinsic tryptophan fluorescence of beta-LG indicated that HHP treatment conditions induced a conformational change. HHP treatment conditions also promote a 3-fold increase in the extrinsic fluorescence of 1-anilinonaphthalene-8-sulfonate and a 2.6-fold decrease for cis-paraneric acid, suggesting an increase in accessible aromatic hydrophobicity and a decrease in aliphatic hydrophobicity. Far-ultraviolet circular dichroism (CD) spectra reveal that the secondary structure of beta-LG converts from native beta-sheets to non-native alpha-helices following HHP treatment, whereas near-ultraviolet CD spectra reveal that the native tertiary structure of beta-LG essentially disappears. Urea titrations reveal that native beta-LG unfolds cooperatively, but the pressure-treated molecule unfolds noncooperatively. The noncooperative state is stable for 3 months at 5 degrees C. The nonaccessible free thiol group of cysteine121 in native beta-LG became reactive to Ellman's reagent after adequate HHP treatment. Gel electrophoresis with and without beta-mercaptoethanol provided evidence that the exposed thiol group was lost concomitant with the formation of S-S-linked beta-LG dimers. Overall, these results suggest that HHP treatments induce beta-LG into hydrophobic molten globule structures that remain stable for at least 3 months.

Aberrant mobility phenomena of the DNA repair protein XPA.

The DNA repair protein XPA recognizes a wide variety of bulky lesions and interacts with several other proteins during nucleotide excision repair. We recently identified regions of intrinsic order and disorder in full length Xenopus XPA (xXPA) protein using an experimental approach that combined time-resolved trypsin proteolysis and electrospray ionization interface coupled to a Fourier transform ion cyclotron resonance (ESI-FTICR) mass spectrometry (MS). MS data were consistent with the interpretation that xXPA contains no post-translational modifications. Here we characterize the discrepancy between the calculated molecular weight (31 kDa) for xXPA and its apparent molecular weight on SDS-PAGE (multiple bands from approximately 40-45 kDa) and gel filtration chromatography ( approximately 92 kDa), as well as the consequences of DNA binding on its anomalous mobility. Iodoacetamide treatment of xXPA prior to SDS-PAGE yielded a single 42-kDa band, showing that covalent modification of Cys did not correct aberrant mobility. Determination of sulfhydryl content in xXPA with Ellman's reagent revealed that all nine Cys in active protein are reduced. Unexpectedly, structural constraints induced by intramolecular glutaraldehyde crosslinks in xXPA produced a approximately 32-kDa monomer in closer agreement with its calculated molecular weight. To investigate whether binding to DNA alters xXPA's anomalous migration, we used gel filtration chromatography. For the first time, we purified stable complexes of xXPA and DNA +/- cisplatin +/- mismatches. xXPA showed at least 10-fold higher affinity for cisplatin DNA +/- mismatches compared to undamaged DNA +/- mismatches. In all cases, DNA binding did not correct xXPA's anomalous migration. To test predictions that a Glu-rich region (EEEEAEE) and/or disordered N- and C-terminal domains were responsible for xXPA's aberrant mobility, the molecular weights of partial proteolytic fragments from approximately 5 to 25 kDa separated by reverse-phase HPLC and precisely determined by ESI-FTICR MS were correlated with their migration on SDS-PAGE. Every partial tryptic fragment analyzed within this size range exhibited 10%-50% larger molecular weights than expected. Thus, both the disordered domains and the Glu-rich region in xXPA are primarily responsible for the aberrant mobility phenomena.

Identification of intrinsic order and disorder in the DNA repair protein XPA.

The DNA-repair protein XPA is required to recognize a wide variety of bulky lesions during nucleotide excision repair. Independent NMR solution structures of a human XPA fragment comprising approximately 40% of the full-length protein, the minimal DNA-binding domain, revealed that one-third of this molecule was disordered. To better characterize structural features of full-length XPA, we performed time-resolved trypsin proteolysis on active recombinant Xenopus XPA (xXPA). The resulting proteolytic fragments were analyzed by electrospray ionization interface coupled to a Fourier transform ion cyclotron resonance mass spectrometry and SDS-PAGE. The molecular weight of the full-length xXPA determined by mass spectrometry (30922.02 daltons) was consistent with that calculated from the sequence (30922.45 daltons). Moreover, the mass spectrometric data allowed the assignment of multiple xXPA fragments not resolvable by SDS-PAGE. The neural network program Predictor of Natural Disordered Regions (PONDR) applied to xXPA predicted extended disordered N- and C-terminal regions with an ordered internal core. This prediction agreed with our partial proteolysis results, thereby indicating that disorder in XPA shares sequence features with other well-characterized intrinsically unstructured proteins. Trypsin cleavages at 30 of the possible 48 sites were detected and no cleavage was observed in an internal region (Q85-I179) despite 14 possible cut sites. For the full-length xXPA, there was strong agreement among PONDR, partial proteolysis data, and the NMR structure for the corresponding XPA fragment.

Intrinsically disordered protein.

Proteins can exist in a trinity of structures: the ordered state, the molten globule, and the random coil. The five following examples suggest that native protein structure can correspond to any of the three states (not just the ordered state) and that protein function can arise from any of the three states and their transitions. (1) In a process that likely mimics infection, fd phage converts from the ordered into the disordered molten globular state. (2) Nucleosome hyperacetylation is crucial to DNA replication and transcription; this chemical modification greatly increases the net negative charge of the nucleosome core particle. We propose that the increased charge imbalance promotes its conversion to a much less rigid form. (3) Clusterin contains an ordered domain and also a native molten globular region. The molten globular domain likely functions as a proteinaceous detergent for cell remodeling and removal of apoptotic debris. (4) In a critical signaling event, a helix in calcineurin becomes bound and surrounded by calmodulin, thereby turning on calcineurin's serine/threonine phosphatase activity. Locating the calcineurin helix within a region of disorder is essential for enabling calmodulin to surround its target upon binding. (5) Calsequestrin regulates calcium levels in the sarcoplasmic reticulum by binding approximately 50 ions/molecule. Disordered polyanion tails at the carboxy terminus bind many of these calcium ions, perhaps without adopting a unique structure. In addition to these examples, we will discuss 16 more proteins with native disorder. These disordered regions include molecular recognition domains, protein folding inhibitors, flexible linkers, entropic springs, entropic clocks, and entropic bristles. Motivated by such examples of intrinsic disorder, we are studying the relationships between amino acid sequence and order/disorder, and from this information we are predicting intrinsic order/disorder from amino acid sequence. The sequence-structure relationships indicate that disorder is an encoded property, and the predictions strongly suggest that proteins in nature are much richer in intrinsic disorder than are those in the Protein Data Bank. Recent predictions on 29 genomes indicate that proteins from eucaryotes apparently have more intrinsic disorder than those from either bacteria or archaea, with typically > 30% of eucaryotic proteins having disordered regions of length > or = 50 consecutive residues.

The protein non-folding problem: amino acid determinants of intrinsic order and disorder.

To investigate the determinants of protein order and disorder, three primary and one derivative database of intrinsically disordered proteins were compiled. The segments in each primary database were characterized by one of the following: X-ray crystallography, nuclear magnetic resonance (NMR), or circular dichroism (CD). The derivative database was based on homology. The three primary disordered databases have a combined total of 157 proteins or segments of length.30 with 18,010 residues, while the derivative database contains 572 proteins from 32 families with 52,688 putatively disordered residues. For the four disordered databases, the amino acid compositions were compared with those from a database of ordered structure. Relative to the ordered protein, the intrinsically disordered segments in all four databases were significantly depleted in W, C, F, I, Y, V, L and N, significantly enriched in A, R, G, Q, S, P, E and K, and inconsistently different in H, M, T, and D, suggesting that the first set be called order-promoting and the second set disorder-promoting. Also, 265 amino acid properties were ranked by their ability to discriminate order and disorder and then pruned to remove the most highly correlated pairs. The 10 highest-ranking properties after pruning consisted of 2 residue contact scales, 4 hydrophobicity scales, 3 scales associated with.-sheets and one polarity scale. Using these 10 properties for comparisons of the 3 primary databases suggests that disorder in all 3 databases is very similar, but with those characterized by NMR and CD being the most similar, those by CD and X-ray being next, and those by NMR and X-ray being the least similar.

Sequence complexity of disordered protein.

Intrinsic disorder refers to segments or to whole proteins that fail to self-fold into fixed 3D structure, with such disorder sometimes existing in the native state. Here we report data on the relationships among intrinsic disorder, sequence complexity as measured by Shannon's entropy, and amino acid composition. Intrinsic disorder identified in protein crystal structures, and by nuclear magnetic resonance, circular dichroism, and prediction from amino acid sequence, all exhibit similar complexity distributions that are shifted to lower values compared to, but significantly overlapping with, the distribution for ordered proteins. Compared to sequences from ordered proteins, these variously characterized intrinsically disordered segments and proteins, and also a collection of low-complexity sequences, typically have obviously higher levels of protein-specific subsets of the following amino acids: R, K, E, P, and S, and lower levels of subsets of the following: C, W, Y, I, and V. The Swiss Protein database of sequences exhibits significantly higher amounts of both low-complexity and predicted-to-be-disordered segments as compared to a non-redundant set of sequences from the Protein Data Bank, providing additional data that nature is richer in disordered and low-complexity segments compared to the commonness of these features in the set of structurally characterized proteins.

Isolation and characterization of cDNAs expressed in the early stages of flavonol-induced pollen germination in petunia.

Petunia (Petunia hybrida) pollen requires flavonols (Fl) to germinate. Adding kaempferol to Fl-deficient pollen causes rapid and synchronous germination and tube outgrowth. We exploited this system to identify genes responsive to Fls and to examine the changes in gene expression that occur during the first 0.5 h of pollen germination. We used a subtracted library and differential screening to identify 22 petunia germinating pollen clones. All but two were expressed exclusively in pollen and half of the clones were rare or low abundance cDNAs. RNA gel-blot analysis showed that the steady-state transcript levels of all the clones were increased in response to kaempferol. The sequences showing the greatest response to kaempferol encode proteins that have regulatory or signaling functions and include S/D4, a leucine-rich repeat protein, S/D1, a LIM-domain protein, and D14, a putative Zn finger protein with a heme-binding site. Eight of the clones were novel including S/D10, a cDNA only expressed very late in pollen development and highly up-regulated during the first 0.5 h of germination. The translation product of the S/D3 cDNA shares some features with a neuropeptide that regulates guidance and growth in the tips of extending axons. This study confirmed that the bulk of pollen mRNA accumulates well before germination, but that specific sequences are transcribed during the earliest moments of Fl-induced pollen germination.

Intrinsic protein disorder in complete genomes.

Intrinsic protein disorder refers to segments or to whole proteins that fail to fold completely on their own. Here we predicted disorder on protein sequences from 34 genomes, including 22 bacteria, 7 archaea, and 5 eucaryotes. Predicted disordered segments > or = 50, > or = 40, and > or = 30 in length were determined as well as proteins estimated to be wholly disordered. The five eucaryotes were separated from bacteria and archaea by having the highest percentages of sequences predicted to have disordered segments > or = 50 in length: from 25% for Plasmodium to 41% for Drosophila. Estimates of wholly disordered proteins in the bacteria ranged from 1% to 8%, averaging to 3 +/- 2%, estimates in various archaea ranged from 2 to 11%, plus an apparently anomalous 18%, averaging to 7 +/- 5% that drops to 5 +/- 3% if the high value is discarded. Estimates in the 5 eucarya ranged from 3 to 17%. The putative wholly disordered proteins were often ribosomal proteins, but in addition about equal numbers were of known and unknown function. Overall, intrinsic disorder appears to be a common, with eucaryotes perhaps having a higher percentage of native disorder than archaea or bacteria.

Comparing predictors of disordered protein.

More than 6,000 amino acid sequence attributes were ranked by their conditional probabilities for indicating ordered or disordered protein structure. The top 10 each from several different groups of attributes were merged with still other attributes and then subjected to selection by logistic regression. Evidently, the determination of order or disorder results from the interplay among several attributes, such as average Coordination Number, aromatic content and the numbers of non-polar amino acids, all of which favor the ordered state, and others like Net Charge, Flexibility Index, and the presence of certain polar amino acids, all of which favor disorder. The top 12 selected attributes were used as inputs for artificial neural network (ANN) predictors. Five predictors were developed, compared with each other, and with previous work. The best of these shows substantially improved generalization compared to our previously published predictor.

Isolation of chloroform-resistant mutants of filamentous phage: localization in models of phage structure.

Interaction of fd or M13 filamentous phage with a chloroform/water interface induces morphological change, contracting the filaments sequentially into shortened rods (I-forms), and then into spheroidal particles (S-forms). To further investigate this phage contraction, 34 and 26 chloroform-resistant isolates of fd and M13, respectively, were selected after chloroform treatment of wild-type phages at pH 8. 2 and 4 degrees C. DNA sequencing of gene VIII of the 34 fd isolates revealed five different mutants: these were D5H, M28L, V31L, I37T, and S50T. All 26 M13 isolates were I37T. These mutants exhibited variable sensitivity to chloroform, but all contracted much more slowly than wild-type phage during treatment at 4 degrees C. They all contracted like wild-type phage at 37 degrees C. Site-directed mutagenesis showed that the indicated single mutations carried the chloroform resistance. In structural models of the phage, the D5H locus is on the outside and the S50T locus is on the inside. The M28L and I37T loci are buried in a mostly hydrophobic region in the middle. Although these four mutants are spread out radially, they are localized in the axial direction into a thin disk in the model. The last mutant locus, V31L, is out of this disk, but this locus is proximal to the M28L and I37T loci and also in contact with the surface via a deep hydrophobic hole or depression. These five mutants, their locations, and their variable affects on contraction suggest that chloroform-induced contraction involves a specific mechanism rather than a generalized solvent-induced denaturation and that the critical structural changes occur in a localized level in the phage. These results add weight to suggestions that the sequential contraction of filaments-->I-forms-->S-forms mimic corresponding steps in phage penetration, and, in the reverse order, for phage assembly.

Folding minimal sequences: the lower bound for sequence complexity of globular proteins.

Alphabet size and informational entropy, two formal measures of sequence complexity, are herein applied to two prior studies on the folding of minimal proteins. These measures show a designed four-helix bundle to be unlike its natural counterparts but rather more like a coiled-coil dimer. Segments from a simplified sarc homology 3 domain and more than 2000000 segments from globular proteins both have lower bounds for alphabet size of 10 and for entropy near 2.9. These values are therefore suggested to be necessary and sufficient for folding into globular proteins having both rigid side chain packing and biological function.

Thousands of proteins likely to have long disordered regions.

Neural network predictors of protein disorder using primary sequence information were developed and applied to the Swiss Protein Database. More than 15,000 proteins were predicted to contain disordered regions of at least 40 consecutive amino acids, with more than 1,000 having especially high scores indicating disorder. These results support proposals that consideration of structure-activity relationships in proteins need to be broadened to include unfolded or disordered protein.

Protein disorder and the evolution of molecular recognition: theory, predictions and observations.

Observations going back more than 20 years show that regions in proteins with disordered backbones can play roles in their binding to other molecules; typically, the disordered regions become ordered upon complex formation. Thought-experiments with Schulz Diagrams, which are defined herein, suggest that disorder-to-order transitions are required for natural selection to operate separately on affinity and specificity. Separation of affinity and specificity may be essential for fine-tuning the molecular interaction networks that comprise the living state. For low affinity, high specificity interactions, our analysis suggests that natural selection would parse the amino acids conferring flexibility in the unbound state from those conferring specificity in the bound state. For high affinity, low specificity or for high affinity, multiple specificity interactions, our analysis suggests that the disorder-to-order transitions enable alternative packing interactions between side chains to accommodate the different binding targets. Disorder-to-order transitions upon binding also have significant kinetic implications as well, by having complex effects on both on- and off-rates. Current data are insufficient to decide on these proposals, but sequence and structure analysis on two examples support further investigations of the role of disorder-to-order transitions upon binding.

Effects of relative humidity and applied force on atomic force microscopy images of the filamentous phage fd.

The filamentous phage fd was studied by both contact- and tapping-mode atomic force microscopy under conditions of controlled variations in relative humidity and changes in the applied tip force. By spin-coating freshly cleaved mica with phage containing solutions having very low salt content followed by rapid humidity control, stable and reliable sample preparation was achieved. The apparent height of the phage varied by about 10-fold with a quadratic dependence on the stabilized relative humidity, extrapolating to 73% of the accepted X-ray diffraction-based height at 0% relative humidity. The variation in measured height with relative humidity largely reconciles previous widely varying atomic force microscopy estimates of this dimension for the filamentous phage. Our finding that contact-mode images of phage are more difficult to analyze than those acquired in tapping mode are consistent with previously published results on other biological specimens such as DNA.

Crystal structure of calsequestrin from rabbit skeletal muscle sarcoplasmic reticulum.

Calsequestrin, the major Ca2+ storage protein of muscle, coordinately binds and releases 40-50 Ca2+ ions per molecule for each contraction-relaxation cycle by an uncertain mechanism. We have determined the structure of rabbit skeletal muscle calsequestrin. Three very negative thioredoxin-like domains surround a hydrophilic center. Each monomer makes two extensive dimerization contacts, both of which involve the approach of many negative groups. This structure suggests a mechanism by which calsequestrin may achieve high capacity Ca2+ binding. The suggested mechanism involves Ca2+-induced collapse of the three domains and polymerization of calsequestrin monomers arising from three factors: N-terminal arm exchange, helix-helix contacts and Ca2+ cross bridges. This proposed structure-based mechanism accounts for the observed coupling of high capacity Ca2+ binding with protein precipitation.

Quaternion contact ribbons: a new tool for visualizing intra- and intermolecular interactions in proteins.

Protein side chain interactions between residues separated by at least one loop or turn or break in the amino acid sequence are called 'nonlocal contacts' in this manuscript, and contiguous sets of such interactions located between segments of secondary structure are called 'contact zones.' A new interactive program, the quaternion contact ribbon tool, has been developed to help protein chemists identify, straighten if twisted, and display contact zones between two neighboring segments of helix.