Interior and surface of monomeric proteins

Deciphering the shape and deformation of secondary structures through local conformation analysis

Background

Protein deformation has been extensively analysed through global methods based on RMSD, torsion angles and Principal Components Analysis calculations. Here we use a local approach, able to distinguish among the different backbone conformations within loops, α-helices and β-strands, to address the question of secondary structures' shape variation within proteins and deformation at interface upon complexation.

Results

Using a structural alphabet, we translated the 3 D structures of large sets of protein-protein complexes into sequences of structural letters. The shape of the secondary structures can be assessed by the structural letters that modeled them in the structural sequences. The distribution analysis of the structural letters in the three protein compartments (surface, core and interface) reveals that secondary structures tend to adopt preferential conformations that differ among the compartments. The local description of secondary structures highlights that curved conformations are preferred on the surface while straight ones are preferred in the core. Interfaces display a mixture of local conformations either preferred in core or surface. The analysis of the structural letters transition occurring between protein-bound and unbound conformations shows that the deformation of secondary structure is tightly linked to the compartment preference of the local conformations.

Conclusion

The conformation of secondary structures can be further analysed and detailed thanks to a structural alphabet which allows a better description of protein surface, core and interface in terms of secondary structures' shape and deformation. Induced-fit modification tendencies described here should be valuable information to identify and characterize regions under strong structural constraints for functional reasons.

A Contribution from the Hematology at Biomaterial Interfaces Research Group, Kao P, Parhi P, Krishnan A, Noh H, Haider W, Tadigadapa S, Allara DL, Vogler EA. Volumetric interpretation of protein adsorption: interfacial packing of protein adsorbed to hydrophobic surfaces from surface-saturating solution concentrations

The maximum capacity of a hydrophobic adsorbent is interpreted in terms of square or hexagonal (cubic and face-centered-cubic, FCC) interfacial packing models of adsorbed blood proteins in a way that accommodates experimental measurements by the solution-depletion method and quartz-crystal-microbalance (QCM) for the human proteins serum albumin (HSA, 66 kDa), immunoglobulin G (IgG, 160 kDa), fibrinogen (Fib, 341 kDa), and immunoglobulin M (IgM, 1000 kDa). A simple analysis shows that adsorbent capacity is capped by a fixed mass/volume (e.g. mg/mL) surface-region (interphase) concentration and not molar concentration. Nearly analytical agreement between the packing models and experiment suggests that, at surface saturation, above-mentioned proteins assemble within the interphase in a manner that approximates a well-ordered array. HSA saturates a hydrophobic adsorbent with the equivalent of a single square or hexagonally-packed layer of hydrated molecules whereas the larger proteins occupy two-or-more layers, depending on the specific protein under consideration and analytical method used to measure adsorbate mass (solution depletion or QCM). Square or hexagonal (cubic and FCC) packing models cannot be clearly distinguished by comparison to experimental data. QCM measurement of adsorbent capacity is shown to be significantly different than that measured by solution depletion for similar hydrophobic adsorbents. The underlying reason is traced to the fact that QCM measures contribution of both core protein, water of hydration, and interphase water whereas solution depletion measures only the contribution of core protein. It is further shown that thickness of the interphase directly measured by QCM systematically exceeds that inferred from solution-depletion measurements, presumably because the static model used to interpret solution depletion does not accurately capture the complexities of the viscoelastic interfacial environment probed by QCM.

Protein homomers in point-group assembly: symmetry making and breaking are specific and distinctive in their codes of chemical alphabet in side chains

Oligomerizing to point-group symmetry, protein oligomers need to have the symmetry broken for biologically crucial functions, such as, allosteric regulation, enzyme catalysis, and so forth. In the making of symmetry, based on self assembly, and the breaking of symmetry, based on intermolecular interactions, proteins may manifest, like their other functions, specific scripts over the coding alphabet in side chains. To address the possibility, we analyzed 82 protein homodimers in their C(2)-symmetry-related side chains across noncrystallographic interfaces, to know if they may be identical or distinct in conformation, and thus conserved or broken in symmetry. We find the propensity to conformational mismatch across interfaces correlated with side-chain chemical structure, low to very low in aromatic Trp, Tyr, His, Phe, and Arg, and high to very high in aliphatic Val, Pro, Met, Glu, Ser, Lys, Gln, Asn, and Asp, related not to polarity but, interestingly, to aromaticity of the structure. The organizational plan having aromatics embedded in a hub of aliphatic-nonpolar groups and a surrounding rim of aliphatic-polar groups, called "hotspot," has been known to direct protein-protein interaction. Finding conformational-mismatch propensities of side chains congruous with their specific chemical roles in protein-protein interaction, we propose that aromatic side chains will drive protein homomers to high symmetry, while polar- and nonpolar aliphatic side chains will drive them to the functionally-necessitated breaks of symmetry. Side chains are in their roles as protein-coding alphabet illuminated in the physics, which is discussed.

A generic approach to evaluate how B-cell epitopes are surface-exposed on protein structures

Methods that predict antibody epitopes could help to promote the development of diagnostic tools, vaccines or immunotherapies by affecting the epitope binding of antibodies during an immunological response to antigens.

It is generally assumed that there is a direct relationship between antibody accessibility to antigens and accessible surface of proteins. Based on this assumption, prediction systems often includes solvent accessibility values calculated from the primary sequence of proteins or from their three dimensional structures as a predictive criterion. However, the current prediction systems seem weakly efficient in view of benchmark tests.

We were interested in evaluating how amino acids that have been experimentally identified as epitopic elements could differ from the rest of the antigenic molecule at the level of surface exposure, hence we assessed the average accessibility of epitopes. The approach used here utilises published epitopes deduced from numerous identification techniques, including sequence scanning and structure visualisation after crystallography, and it involves many types of antigens from toxins to allergens. Our results show that epitopic residues are not distributed among any specific Relative Surface Accessibility and Protrusion Index values and that, in some cases, epitopes cover the entire antigenic sequence.

These results led to the conclusion that the classification of known epitopes with respect to the experimental conditions used to identify them should be introduced before attempting to characterise epitopic areas in a generic way.

Effect of bioparticle size on dispersion and retention in monolithic and perfusive beds

Single-component pulse response studies were used to compare the retention and transport behavior of small molecules, proteins, and a virus on commercially available monolithic and perfusive ion-exchangers. Temporal distortion and extra-column effects were corrected for using a simple algorithm based on the method of moments. It was found that temporal distortion is inversely related to the number of theoretical plates. With increasing bioparticle size, retention increased and the transition from a non-eluting to a non-adsorbing state with increasing ionic strength became more abrupt. Both of these observations are qualitatively explained by calculations of particle-surface electrostatic attractive energy. Calculations also suggest that, for sufficiently large bioparticles, such as viruses or cells, hydrodynamic drag can promote elution. Under non-adsorbing conditions, plate height increased only weakly with flow rate and the skew remained unchanged. With increasing retention, plate height increased dramatically for proteins. Plate height was scaled by permeability rather than bead diameter to enable comparison among different stationary phases.

Phospho.ELM: a database of phosphorylation sites--update 2011

The Phospho.ELM resource (http://phospho.elm.eu.org) is a relational database designed to store in vivo and in vitro phosphorylation data extracted from the scientific literature and phosphoproteomic analyses. The resource has been actively developed for more than 7 years and currently comprises 42 574 serine, threonine and tyrosine non-redundant phosphorylation sites. Several new features have been implemented, such as structural disorder/order and accessibility information and a conservation score. Additionally, the conservation of the phosphosites can now be visualized directly on the multiple sequence alignment used for the score calculation. Finally, special emphasis has been put on linking to external resources such as interaction networks and other databases.

Phospho3D 2.0: an enhanced database of three-dimensional structures of phosphorylation sites

Phospho3D is a database of three-dimensional (3D) structures of phosphorylation sites (P-sites) derived from the Phospho.ELM database, which also collects information on the residues surrounding the P-site in space (3D zones). The database also provides the results of a large-scale structural comparison of the 3D zones versus a representative dataset of structures, thus associating to each P-site a number of structurally similar sites. The new version of Phospho3D presents an 11-fold increase in the number of 3D sites and incorporates several additional features, including new structural descriptors, the possibility of selecting non-redundant sets of 3D structures and the availability for download of non-redundant sets of structurally annotated P-sites. Moreover, it features P3Dscan, a new functionality that allows the user to submit a protein structure and scan it against the 3D zones collected in the Phospho3D database. Phospho3D version 2.0 is available at: http://www.phospho3d.org/.

Underexposed polar residues and protein stabilization

Increasing protein stability is interesting for practical reasons and because it tests our understanding of protein energetics. We explore here the feasibility of stabilizing proteins by replacing underexposed polar residues by apolar ones of similar size and shape. We have compared the stability of wild-type apoflavodoxin with that of a few carefully selected mutants carrying Y → F, Q → L, T → V or K → M replacements. Although a clear inverse correlation between native solvent exposures of replaced polar residues and stability of mutants is observed, most mutations fail to stabilize the protein. The promising exceptions are the two Q → L mutations tested, which characteristically combine the greatest reduction in polar burial with the greatest increase in apolar burial relative to wild type. Analysis of published stability data corresponding to a variety of mutant proteins confirms that, unlike Y → F or T → V replacements, Q → L mutations tend to be stabilizing, and it suggests that N → L mutations might be stabilizing as well. On the other hand, we show that the stability changes associated to the apoflavodoxin mutations can be rationalized in terms of differential polar and apolar burials upon folding plus a generic destabilizing penalty term. Simple equations combining these contributions predict stability changes in a large data set of 113 mutants (Y → F, Q → L or T → V) similarly well as more complex algorithms available on the Internet.

A solution NMR investigation into the murine amelogenin splice-variant LRAP (Leucine-Rich Amelogenin Protein)

Amelogenins are the dominant proteins present in ameloblasts during the early stages of enamel biomineralization, making up >90% of the matrix protein. Along with the full-length protein there are several splice-variant isoforms of amelogenin present including LRAP (Leucine-Rich Amelogenin Protein), a protein that consists of the first 33 and the last 26 residues of full-length amelogenin. Using solution-state NMR spectroscopy we have assigned the (1)H-(15)N HSQC spectrum of murine LRAP (rp(H)LRAP) in 2% acetic acid at pH 3.0 by making extensive use of previous chemical shift assignments for full-length murine amelogenin (rp(H)M180). This correlation was possible because LRAP, like the full-length protein, is intrinsically disordered under these solution conditions. The major difference between the (1)H-(15)N HSQC spectra of rp(H)M180 and rp(H)LRAP was an additional set of amide resonances for each of the seven non-proline residues between S12 and Y12 near the N-terminus of rp(H)LRAP indicating that the N-terminal region of LRAP exists in two different conformations. Analysis of the proline carbon chemical shifts suggests that the molecular basis for the two states is not a cis-trans isomerization of one or more of the proline residues in the N-terminal region. Starting from 2% acetic acid, where rp(H)LRAP was monomeric in solution, NaCl addition effected residue specific changes in molecular dynamics manifested by the reduction in intensity and disappearance of (1)H-(15)N HSQC cross peaks. As observed for the full-length protein, these perturbations may signal early events governing supramolecular self-assembly of rp(H)LRAP into nanospheres. However, the different patterns of (1)H-(15)N HSQC cross peak perturbation between rp(H)LRAP and rp(H)M180 in high salt suggest that the termini may behave differently in their respective nanospheres, and perhaps, these differences contribute to the cell signaling properties attributable to LRAP but not to the full-length protein.

Adhesion through single peptide aptamers

Aptamer and antibody mediated adhesion is central to biological function and is valuable in the engineering of "lab on a chip" devices. Single molecule force spectroscopy using optical tweezers enables direct nonequilibrium measurement of these noncovalent interactions for three peptide aptamers selected for glass, polystyrene, and carbon nanotubes. A comprehensive examination of the strong attachment between antifluorescein 4-4-20 and fluorescein was also carried out using the same assay. Bond lifetime, barrier width, and free energy of activation are extracted from unbinding histogram data using three single molecule pulling models. The evaluated aptamers appear to adhere stronger than the fluorescein antibody under no- and low-load conditions, yet weaker than antibodies at loads above ∼25 pN. Comparison to force spectroscopy data of other biological linkages shows the diversity of load dependent binding and provides insight into linkages used in biological processes and those designed for engineered systems.

Systematic analysis of short internal indels and their impact on protein folding

Background

Protein sequence insertions/deletions (indels) can be introduced during evolution or through alternative splicing (AS). Alternative splicing is an important biological phenomenon and is considered as the major means of expanding structural and functional diversity in eukaryotes. Knowledge of the structural changes due to indels is critical to our understanding of the evolution of protein structure and function. In addition, it can help us probe the evolution of alternative splicing and the diversity of functional isoforms. However, little is known about the effects of indels, in particular the ones involving core secondary structures, on the folding of protein structures. The long term goal of our study is to accurately predict the protein AS isoform structures. As a first step towards this goal, we performed a systematic analysis on the structural changes caused by short internal indels through mining highly homologous proteins in Protein Data Bank (PDB).

Results

We compiled a non-redundant dataset of short internal indels (2-40 amino acids) from highly homologous protein pairs and analyzed the sequence and structural features of the indels. We found that about one third of indel residues are in disordered state and majority of the residues are exposed to solvent, suggesting that these indels are generally located on the surface of proteins. Though naturally occurring indels are fewer than engineered ones in the dataset, there are no statistically significant differences in terms of amino acid frequencies and secondary structure types between the "Natural" indels and "All" indels in the dataset. Structural comparisons show that all the protein pairs with short internal indels in the dataset preserve the structural folds and about 85% of protein pairs have global RMSDs (root mean square deviations) of 2Å or less, suggesting that protein structures tend to be conserved and can tolerate short insertions and deletions. A few pairs with high RMSDs are results of relative domain positions of the proteins, probably due to the intrinsically dynamic nature of the proteins.

Conclusions

The analysis demonstrated that protein structures have the "plasticity" to tolerate short indels. This study can provide valuable guides in modeling protein AS isoform structures and homologous proteins with indels through placing the indels at the right locations since the accuracy of sequence alignments dictate model qualities in homology modeling.

Determinants of affinity and activity of the anti-sigma factor AsiA

The AsiA protein is a T4 bacteriophage early gene product that regulates transcription of host and viral genes. Monomeric AsiA binds tightly to the sigma(70) subunit of Escherichia coli RNA polymerase, thereby inhibiting transcription from bacterial promoters and phage early promoters and coactivating transcription from phage middle promoters. Results of structural studies have identified amino acids at the protomer-protomer interface in dimeric AsiA and at the monomeric AsiA-sigma(70) interface and demonstrated substantial overlap in the sets of residues that comprise each. Here we evaluate the contributions of individual interfacial amino acid side chains to protomer-protomer affinity in AsiA homodimers, to monomeric AsiA affinity for sigma(70), and to AsiA function in transcription. Sedimentation equilibrium, dynamic light scattering, electrophoretic mobility shift, and transcription activity measurements were used to assess affinity and function of site-specific AsiA mutants. Alanine substitutions for solvent-inaccessible residues positioned centrally in the protomer-protomer interface of the AsiA homodimer, V14, I17, and I40, resulted in the largest changes in free energy of dimer association, whereas alanine substitutions at other interfacial positions had little effect. These residues also contribute significantly to AsiA-dependent regulation of RNA polymerase activity, as do additional residues positioned at the periphery of the interface (K20 and F21). Notably, the relative contributions of a given amino acid side chain to RNA polymerase inhibition and activation (MotA-independent) by AsiA are very similar in most cases. The mainstay for intermolecular affinity and AsiA function appears to be I17. Our results define the core interfacial residues of AsiA, establish roles for many of the interfacial amino acids, are in agreement with the tenets underlying protein-protein interactions and interfaces, and will be beneficial for a general, comprehensive understanding of the mechanistic underpinnings of bacterial RNA polymerase regulation.

Concentration of specific amino acids at the catalytic/active centers of highly-conserved "housekeeping" enzymes of central metabolism in archaea, bacteria and Eukaryota: is there a widely conserved chemical signal of prebiotic assembly?

In alignments of 1969 protein sequences the amino acid glycine and others were found concentrated at most-conserved sites within approximately 15 A of catalytic/active centers (C/AC) of highly conserved kinases, dehydrogenases or lyases of Archaea, Bacteria and Eukaryota. Lysine and glutamic acid were concentrated at least-conserved sites furthest from their C/ACs. Logistic-regression analyses corroborated the "movement" of glycine towards and lysine away from their C/ACs: the odds of a glycine occupying a site were decreased by 19%, while the odds for a lysine were increased by 53%, for every 10 A moving away from the C/AC. Average conservation of MSA consensus sites was highest surrounding the C/AC and directly decreased in transition toward model's peripheries. Findings held with statistical confidence using sequences restricted to individual Domains or enzyme classes or to both. Our data describe variability in the rate of mutation and likelihoods for phylogenetic trees based on protein sequence data and endorse the extension of substitution models by incorporating data on conservation and distance to C/ACs rather than only using cumulative levels. The data support the view that in the most-conserved environment immediately surrounding the C/AC of taxonomically distant and highly conserved essential enzymes of central metabolism there are amino acids whose identity and degree of occupancy is similar to a proposed amino acid set and frequency associated with prebiotic evolution.

Critical interactions in the stability control region of tropomyosin

Our laboratory has recently described a stability control region in the two-stranded alpha-helical coiled-coil alpha-tropomyosin that accounts for overall protein stability but is not required for folding (Hodges et al., 2009). We have used a synthetic peptide approach to investigate three stability control sites within the stability control region (residues 97-118). Two of the sites, electrostatic cluster 1 (97-104, EELDRAQE) and electrostatic cluster 2 (112-118, KLEEAEK), feature sequences with unusually high charge density and the potential to form multiple intrachain and interchain salt bridges (ionic attractions). A third site (105-111, RLATALQ) features an e position Leu residue, an arrangement known previously to enhance coiled-coil stability modestly. A native peptide and seven peptide analogs of the tropomyosin sequence 85-119 were prepared by Fmoc solid-phase peptide synthesis. Thermal stability measurements by circular dichroism (CD) spectroscopy revealed the following T(m) values for the native peptide and three key analogs: 52.9 degrees C (Native), 46.0 degrees C (R101A), 45.3 degrees C (K112A/K118A), and 27.9 degrees C (L110A). The corresponding DeltaT(m) values for the analogs, relative to the native peptide, are -6.9 degrees C, -7.6 degrees C, and -25.0 degrees C, respectively. The dramatic contribution to stability made by L110e is three times greater than the contribution of either electrostatic cluster 1 or 2, likely resulting from a novel hydrophobic interaction not previously observed. These thermal stability results were corroborated by temperature profiling analyses using reversed-phase high-performance liquid chromatography (RP-HPLC). We believe that the combined contributions of the interactions within the three stability control sites are responsible for the effect of the stability control region in tropomyosin, with the Leu110e contribution being most critical.

Depletion of highly abundant proteins in blood plasma by hydrophobic interaction chromatography for proteomic analysis

The proteomic analysis of plasma is extremely complex due to the presence of few highly abundant proteins. These proteins have to be depleted in order to detect low abundance proteins, which are likely to be of biomedical interest. In this work it was investigated the applicability of hydrophobic interaction chromatography (HIC) as a plasma fractionation method prior to two-dimensional gel electrophoresis (2DGE). The average hydrophobicity of the 56 main plasma proteins was calculated. Plasma proteins were classified as low, medium and highly hydrophobic through a cluster analysis. The highly abundant proteins showed a medium hydrophobicity, and therefore a HIC step was designed to deplete them from plasma. HIC performance was assessed by 2DGE, and it was compared to that obtained by a commercial immuno-affinity (IA) column for albumin depletion. Both methods showed similar reproducibility. HIC allowed partially depleting alpha-1-antitrypsin and albumin, and permitted to detect twice the number of spots than IA. Since albumin depletion by HIC was incomplete, it should be further optimized for its use as a complementary or alternative method to IA.

The subunit interfaces of weakly associated homodimeric proteins

We analyzed subunit interfaces in 315 homodimers with an X-ray structure in the Protein Data Bank, validated by checking the literature for data that indicate that the proteins are dimeric in solution and that, in the case of the "weak" dimers, the homodimer is in equilibrium with the monomer. The interfaces of the 42 weak dimers, which are smaller by a factor of 2.4 on average than in the remainder of the set, are comparable in size with antibody-antigen or protease-inhibitor interfaces. Nevertheless, they are more hydrophobic than in the average transient protein-protein complex and similar in amino acid composition to the other homodimer interfaces. The mean numbers of interface hydrogen bonds and hydration water molecules per unit area are also similar in homodimers and transient complexes. Parameters related to the atomic packing suggest that many of the weak dimer interfaces are loosely packed, and we suggest that this contributes to their low stability. To evaluate the evolutionary selection pressure on interface residues, we calculated the Shannon entropy of homologous amino acid sequences at 60% sequence identity. In 93% of the homodimers, the interface residues are better conserved than the residues on the protein surface. The weak dimers display the same high degree of interface conservation as other homodimers, but their homologs may be heterodimers as well as homodimers. Their interfaces may be good models in terms of their size, composition, and evolutionary conservation for the labile subunit contacts that allow protein assemblies to share and exchange components, allosteric proteins to undergo quaternary structure transitions, and molecular machines to operate in the cell.

A computational pathway for bracketing native-like structures fo small alpha helical globular proteins

Impressive advances in the applications of bioinformatics for protein structure prediction coupled with growing structural databases on one hand and the insurmountable time-scale problem with ab initio computational methods on the other continue to raise doubts whether a computational solution to the protein folding problem--categorized as an NP-hard problem--is within reach in the near future. Combining some specially designed biophysical filters and vector algebra tools with ab initio methods, we present here a promising computational pathway for bracketing native-like structures of small alpha helical globular proteins departing from secondary structural information. The automated protocol is initiated by generating multiple structures around the loops between secondary structural elements. A set of knowledge-based biophysical filters namely persistence length and radius of gyration, developed and calibrated on approximately 1000 globular proteins, is introduced to screen the trial structures to filter out improbable candidates for the native and reduce the size of the library of probable structures. The ensemble so generated encompasses a few structures with native-like topology. Monte Carlo optimizations of the loop dihedrals are then carried out to remove steric clashes. The resultant structures are energy minimized and ranked according to a scoring function tested previously on a series of decoy sets vis-a-vis their corresponding natives. We find that the 100 lowest energy structures culled from the ensemble of energy optimized trial structures comprise at least a few to within 3-5 angstroms of the native. Thus the formidable "needle in a haystack" problem is narrowed down to finding an optimal solution amongst a computationally tractable number of alternatives. Encouraging results obtained on twelve small alpha helical globular proteins with the above outlined pathway are presented and discussed.

Surface-energy dependent contact activation of blood factor XII

Contact activation of blood factor XII (FXII, Hageman factor) in neat-buffer solution exhibits a parabolic profile when scaled as a function of silanized-glass-particle activator surface energy (measured as advancing water adhesion tension tau(a)(o)=gamma(lv)(o)cos theta in dyne/cm, where gamma(lv)(o) is water interfacial tension in dyne/cm and theta is the advancing contact angle). Nearly equal activation is observed at the extremes of activator water-wetting properties -36<tau(a)(o)<72 dyne/cm (0 degrees <or=theta<120 degrees), falling sharply through a broad minimum within the 20<tau(a)(o)<40 dyne/cm (55 degrees <theta<75 degrees) range over which activation yield (putatively FXIIa) rises just above detection limits. Activation is very rapid upon contact with all activators tested and did not significantly vary over 30 min of continuous FXII-procoagulant contact. Results suggest that materials falling within the 20<tau(a)(o)<40 dyne/cm surface-energy range should exhibit minimal activation of blood-plasma coagulation through the intrinsic pathway. Surface chemistries falling within this range are, however, a perplexingly difficult target for surface engineering because of the critical balance that must be struck between hydrophobicity and hydrophilicity. Results are interpreted within the context of blood plasma coagulation and the role of water and proteins at procoagulant surfaces.

Volumetric interpretation of protein adsorption: kinetics of protein-adsorption competition from binary solution

The standard solution-depletion method is implemented with SDS-gel electrophoresis as a multiplexing, separation-and-quantification tool to measure competition between two proteins (i and j) for adsorption to the same hydrophobic adsorbent particles (either octyl sepharose or silanized glass) immersed in binary-protein solutions. Adsorption kinetics reveals an unanticipated slow protein-size-dependent competition that controls steady-state adsorption selectivity. Two sequential pseudo-steady-state adsorption regimes (State 1 and State 2) are frequently observed depending on i, j solution concentrations. State 1 and State 2 are connected by a smooth transition, giving rise to sigmoidally-shaped adsorption-kinetic profiles with a downward inflection near 60 min of solution/adsorbent contact. Mass ratio of adsorbed i, j proteins (m(i)/m(j)) remains nearly constant between States 1 and 2, even though both m(i) and m(j) decrease in the transition between states. State 2 is shown to be stable for 24 h of continuous-adsorbent contact with stagnant solution whereas State 2 is eliminated by continuous mixing of adsorbent with solution. In sharp contrast to binary-competition results, adsorption to hydrophobic adsorbent particles from single-protein solutions (pure i or j) exhibits no detectable kinetics within the timeframe of experiment from either stagnant or continuously mixed solution, quickly achieving a single steady-state value in proportion to solution concentration. Comparison of binary competition between dissimilarly-sized protein pairs chosen to span a broad molecular-weight (MW) range demonstrates that selectivity between i and j scales with MW ratio that is proportional to protein-volume ratio (ubiquitin, Ub, MW=10.7 kDa; human serum albumin, HSA, MW=66.3 kDa; prothrombin, FII, 72 kDa; immunoglobulin G, IgG, MW=160 kDa; fibrinogen, Fib, MW=341 kDa). Results are interpreted in terms of a kinetic model of adsorption that has protein molecules rapidly diffusing into an inflating interphase that is spontaneously formed by bringing a protein solution into contact with a physical surface (State 1). State 2 follows by rearrangement of proteins within this interphase to achieve the maximum interphase concentration (dictated by energetics of interphase dehydration) within the thinnest (lowest volume) interphase possible by ejection of interphase water and initially-adsorbed proteins. Implications for understanding biocompatibility are discussed using a computational example relevant to the problem of blood-plasma coagulation.

A structure filter for the Eukaryotic Linear Motif Resource

Background

Many proteins are highly modular, being assembled from globular domains and segments of natively disordered polypeptides. Linear motifs, short sequence modules functioning independently of protein tertiary structure, are most abundant in natively disordered polypeptides but are also found in accessible parts of globular domains, such as exposed loops. The prediction of novel occurrences of known linear motifs attempts the difficult task of distinguishing functional matches from stochastically occurring non-functional matches. Although functionality can only be confirmed experimentally, confidence in a putative motif is increased if a motif exhibits attributes associated with functional instances such as occurrence in the correct taxonomic range, cellular compartment, conservation in homologues and accessibility to interacting partners. Several tools now use these attributes to classify putative motifs based on confidence of functionality.

Results

Current methods assessing motif accessibility do not consider much of the information available, either predicting accessibility from primary sequence or regarding any motif occurring in a globular region as low confidence. We present a method considering accessibility and secondary structural context derived from experimentally solved protein structures to rectify this situation. Putatively functional motif occurrences are mapped onto a representative domain, given that a high quality reference SCOP domain structure is available for the protein itself or a close relative. Candidate motifs can then be scored for solvent-accessibility and secondary structure context. The scores are calibrated on a benchmark set of experimentally verified motif instances compared with a set of random matches. A combined score yields 3-fold enrichment for functional motifs assigned to high confidence classifications and 2.5-fold enrichment for random motifs assigned to low confidence classifications. The structure filter is implemented as a pipeline with both a graphical interface via the ELM resource and through a Web Service protocol.

Conclusion

New occurrences of known linear motifs require experimental validation as the bioinformatics tools currently have limited reliability. The ELM structure filter will aid users assessing candidate motifs presenting in globular structural regions. Most importantly, it will help users to decide whether to expend their valuable time and resources on experimental testing of interesting motif candidates.

Volumetric interpretation of protein adsorption: capacity scaling with adsorbate molecular weight and adsorbent surface energy

Silanized-glass-particle adsorbent capacities are extracted from adsorption isotherms of human serum albumin (HSA, 66 kDa), immunoglobulin G (IgG, 160 kDa), fibrinogen (Fib, 341 kDa), and immunoglobulin M (IgM, 1000 kDa) for adsorbent surface energies sampling the observable range of water wettability. Adsorbent capacity expressed as either mass-or-moles per-unit-adsorbent-area increases with protein molecular weight (MW) in a manner that is quantitatively inconsistent with the idea that proteins adsorb as a monolayer at the solution-material interface in any physically-realizable configuration or state of denaturation. Capacity decreases monotonically with increasing adsorbent hydrophilicity to the limit-of-detection (LOD) near tau(o) = 30 dyne/cm (theta approximately 65 degrees) for all protein/surface combinations studied (where tau(o) identical with gamma(lv)(o) costheta is the water adhesion tension, gamma(lv)(o) is the interfacial tension of pure-buffer solution, and theta is the buffer advancing contact angle). Experimental evidence thus shows that adsorbent capacity depends on both adsorbent surface energy and adsorbate size. Comparison of theory to experiment implies that proteins do not adsorb onto a two-dimensional (2D) interfacial plane as frequently depicted in the literature but rather partition from solution into a three-dimensional (3D) interphase region that separates the physical surface from bulk solution. This interphase has a finite volume related to the dimensions of hydrated protein in the adsorbed state (defining "layer" thickness). The interphase can be comprised of a number of adsorbed-protein layers depending on the solution concentration in which adsorbent is immersed, molecular volume of the adsorbing protein (proportional to MW), and adsorbent hydrophilicity. Multilayer adsorption accounts for adsorbent capacity over-and-above monolayer and is inconsistent with the idea that protein adsorbs to surfaces primarily through protein/surface interactions because proteins within second (or higher-order) layers are too distant from the adsorbent surface to be held surface bound by interaction forces in close proximity. Overall, results are consistent with the idea that protein adsorption is primarily controlled by water/surface interactions.

End-to-end and end-to-middle interhelical interactions: new classes of interacting helix pairs in protein structures

Helix-helix interactions are important for the structure, stability and function of alpha-helical proteins. Helices that either cross in the middle or show extensive contacts between each other, such as coiled coils, have been investigated in previous studies. Interactions between two helices can also occur only at the terminal regions or between the terminal region of one helix and the middle region of another helix. Examples of such helix pairs are found in aquaporin, H(+)/Cl(-) transporter and Bcl-2 proteins. The frequency of the occurrence of such ;end-to-end' (EE) and ;end-to-middle' (EM) helix pairs in protein structures is not known. Questions regarding the residue preferences in the interface and the mode of interhelical interactions in such helix pairs also remain unanswered. In this study, high-resolution structures of all-alpha proteins from the PDB have been systematically analyzed and the helix pairs that interact only in EE or EM fashion have been extracted. EE and EM helix pairs have been categorized into five classes (N-N, N-C, C-C, N-MID and C-MID) depending on the region of interaction. Nearly 13% of 5725 helix pairs belonged to one of the five classes. Analysis of single-residue propensities indicated that hydrophobic and polar residues prefer to occur in the C-terminal and N-terminal regions, respectively. Hydrophobic C-terminal interacting residues and polar N-terminal interacting residues are also highly conserved. A strong correlation exists between some of the residue properties (surface area/volume and length of side chains) and their preferences for occurring in the interface of EE and EM helix pairs. In contrast to interacting non-EE/EM helix pairs, helices in EE and EM pairs are farther apart. In these helix pairs, residues with large surface area/volume and longer side chains are preferred in the interfacial region.

Solvent accessible surface area of amino acid residues in globular proteins: correlation of apparent transfer free energies with experimental hydrophobicity scales

It is known that the distribution of amino acid residues in globular proteins between surface and interior is in certain correlation with various experimental scales based on partitioning of amino acids or their analogs between water and organic solvents. These scales are often used in various quantitative structure-activity relationship (QSAR) studies as well as for evaluation of stability of proteins. In this work we have analyzed the distribution of residues based on their solvent accessible surface area in more than 8000 protein structures. Using extensive statistical sampling, we have computed residue apparent free energies of transfer between protein interior and surface applying various criteria for classifying residues as exposed or buried. The correlation of these statistical energies with several experimental hydrophobicity scales is discussed. We propose three types of statistical apparent transfer free energy scales and show that each of these scales is in better correlation with one of the experimental hydrophobicity scales (water/vapor, water/cyclohexane, and water/octanol transfer scales). The data are interpreted through the application of theoretical considerations by Finkelstein et al. (Protein Struct. Funct. Genet. 1995, 23, 142) based on random energy model of heteropolymer globules. The deviation of apparent transfer free energies from experimental scales is discussed and analyzed. The variations of amino acid distribution in proteins with the size of protein structure is discussed and the final protein set is chosen to minimize these variations.

Water molecule contributions to proton spin-lattice relaxation in rotationally immobilized proteins

Spin-lattice relaxation rates of protein and water protons in dry and hydrated immobilized bovine serum albumin were measured in the range of (1)H Larmor frequency from 10 kHz to 30 MHz at temperatures from 154 to 302 K. The water proton spin-lattice relaxation reports on that of protein protons, which causes the characteristic power law dependence on the magnetic field strength. Isotope substitution of deuterium for hydrogen in water and studies at different temperatures expose three classes of water molecule dynamics that contribute to the spin-lattice relaxation dispersion profile. At 185 K, a water (1)H relaxation contribution derives from reorientation of protein-bound molecules that are dynamically uncoupled from the protein backbone and is characterized by a Lorentzian function. Bound-water-molecule motions that can be dynamically uncoupled or coupled to the protein fluctuations make dominant contributions at higher temperatures as well. Surface water translational diffusion that is magnetically two-dimensional makes relaxation contributions at frequencies above 10 MHz. It is shown using isotope substitution that the exponent of the power law of the water signal in hydrated immobilized protein systems is the same as that for protons in lyophilized proteins over four orders of magnitude in the Larmor frequency, which implies that changes in the protein structure associated with hydration do not affect the (1)H spin relaxation.

A sequence and structure based method to predict putative substrates, functions and regulatory networks of endo proteases

BACKGROUND

Proteases play a central role in cellular homeostasis and are responsible for the spatio-temporal regulation of function. Many putative proteases have been recently identified through genomic approaches, leading to a surge in global profiling attempts to characterize their function. Through such efforts and others it has become evident that many proteases play non-traditional roles. Accordingly, the number and the variety of the substrate repertoire of proteases are expected to be much larger than previously assumed. In line with such global profiling attempts, we present here a method for the prediction of natural substrates of endo proteases (human proteases used as an example) by employing short peptide sequences as specificity determinants.

METHODOLOGY/PRINCIPAL FINDINGS

Our method incorporates specificity determinants unique to individual enzymes and physiologically relevant dual filters namely, solvent accessible surface area--a parameter dependent on protein three-dimensional structure and subcellular localization. By incorporating such hitherto unused principles in prediction methods, a novel ligand docking strategy to mimic substrate binding at the active site of the enzyme, and GO functions, we identify and perform subjective validation on putative substrates of matriptase and highlight new functions of the enzyme. Using relative solvent accessibility to rank order we show how new protease regulatory networks and enzyme cascades can be created.

CONCLUSION

We believe that our physiologically relevant computational approach would be a very useful complementary method in the current day attempts to profile proteases (endo proteases in particular) and their substrates. In addition, by using functional annotations, we have demonstrated how normal and unknown functions of a protease can be envisaged. We have developed a network which can be integrated to create a proteolytic world. This network can in turn be extended to integrate other regulatory networks to build a system wide knowledge of the proteome.

Structure of the Notch1-negative regulatory region: implications for normal activation and pathogenic signaling in T-ALL

Proteolytic resistance of Notch prior to ligand binding depends on the structural integrity of a negative regulatory region (NRR) of the receptor that immediately precedes the transmembrane segment. The NRR includes the 3 Lin12/Notch repeats and the juxtamembrane heterodimerization domain, the region of Notch1 most frequently mutated in T-cell acute lymphoblastic leukemia lymphoma (T-ALL). Here, we report the x-ray structure of the Notch1 NRR in its autoinhibited conformation. A key feature of the Notch1 structure that maintains its closed conformation is a conserved hydrophobic plug that sterically occludes the metalloprotease cleavage site. Crystal packing interactions involving a highly conserved, exposed face on the third Lin12/Notch repeat suggest that this site may normally be engaged in intermolecular or intramolecular protein-protein interactions. The majority of known T-ALL-associated point mutations map to residues in the hydrophobic interior of the Notch1 NRR. A novel mutation (H1545P), which alters a residue at the crystal-packing interface, leads to ligand-independent increases in signaling in reporter gene assays despite only mild destabilization of the NRR, suggesting that it releases the autoinhibitory clamp on the heterodimerization domain imposed by the Lin12/Notch repeats. The Notch1 NRR structure should facilitate a search for antibodies or compounds that stabilize the autoinhibited conformation.

A four-column theory for the origin of the genetic code: tracing the evolutionary pathways that gave rise to an optimized code

Background

The arrangement of the amino acids in the genetic code is such that neighbouring codons are assigned to amino acids with similar physical properties. Hence, the effects of translational error are minimized with respect to randomly reshuffled codes. Further inspection reveals that it is amino acids in the same column of the code (i.e. same second base) that are similar, whereas those in the same row show no particular similarity. We propose a 'four-column' theory for the origin of the code that explains how the action of selection during the build-up of the code leads to a final code that has the observed properties.

Results

The theory makes the following propositions. (i) The earliest amino acids in the code were those that are easiest to synthesize non-biologically, namely Gly, Ala, Asp, Glu and Val. (ii) These amino acids are assigned to codons with G at first position. Therefore the first code may have used only these codons. (iii) The code rapidly developed into a four-column code where all codons in the same column coded for the same amino acid: NUN = Val, NCN = Ala, NAN = Asp and/or Glu, and NGN = Gly. (iv) Later amino acids were added sequentially to the code by a process of subdivision of codon blocks in which a subset of the codons assigned to an early amino acid were reassigned to a later amino acid. (v) Later amino acids were added into positions formerly occupied by amino acids with similar properties because this can occur with minimal disruption to the proteins already encoded by the earlier code. As a result, the properties of the amino acids in the final code retain a four-column pattern that is a relic of the earliest stages of code evolution.

Conclusion

The driving force during this process is not the minimization of translational error, but positive selection for the increased diversity and functionality of the proteins that can be made with a larger amino acid alphabet. Nevertheless, the code that results is one in which translational error is minimized. We define a cost function with which we can compare the fitness of codes with varying numbers of amino acids, and a barrier function, which measures the change in cost immediately after addition of a new amino acid. We show that the barrier is positive if an amino acid is added into a column with dissimilar properties, but negative if an amino acid is added into a column with similar physical properties. Thus, natural selection favours the assignment of amino acids to the positions that they occupy in the final code.

Reviewers

This article was reviewed by David Ardell, Eugene Koonin and Stephen Freeland (nominated by Laurence Hurst)

Accessibility and partner number of protein residues, their relationship and a webserver, ContPlot for their display

Background

Depending on chemical features residues have preferred locations – interior or exterior – in protein structures, which also determine how many other residues are found around them. The close packing of residues is the hallmark of protein interior and protein-protein interaction sites.

Results

The average values of accessible surface area (ASA) and partner number (PN, the number of other residues within a distance of 4.5 Å from any atom of a given residue) of different residues have been determined and a webserver, ContPlot has been designed to display these values (relative to the average values) along the protein sequence. This would be useful to visually identify residues that are densely packed, or those involved in protein-protein interactions. The skewness observed in the distribution of PNs is indicative of the hydrophobic or hydrophilic nature of the residue. The variation of ASA with PN can be analytically expressed in terms of a cubic equation. These equations (one for each residue) can be used to estimate the ASA of a polypeptide chain using the PNs of the individual residues in the structure.

Conclusion

The atom-based PNs (obtained by counting surrounding atoms) are highly correlated to the residue-based PN, indicating that the latter can adequately capture the atomic details of packing. The average values of ASA and PN associated with each residue should be useful in protein structure prediction or fold-recognition algorithm. ContPlot would provide a handy tool to assess the importance of a residue in the protein structure or interaction site.

ProtSA: a web application for calculating sequence specific protein solvent accessibilities in the unfolded ensemble

Background

The stability of proteins is governed by the heat capacity, enthalpy and entropy changes of folding, which are strongly correlated to the change in solvent accessible surface area experienced by the polypeptide. While the surface exposed in the folded state can be easily determined, accessibilities for the unfolded state at the atomic level cannot be obtained experimentally and are typically estimated using simplistic models of the unfolded ensemble. A web application providing realistic accessibilities of the unfolded ensemble of a given protein at the atomic level will prove useful.

Results

ProtSA, a web application that calculates sequence-specific solvent accessibilities of the unfolded state ensembles of proteins has been developed and made freely available to the scientific community. The input is the amino acid sequence of the protein of interest. ProtSA follows a previously published calculation protocol which uses the Flexible-Meccano algorithm to generate unfolded conformations representative of the unfolded ensemble of the protein, and uses the exact analytical software ALPHASURF to calculate atom solvent accessibilities, which are averaged on the ensemble.

Conclusion

ProtSA is a novel tool for the researcher investigating protein folding energetics. The sequence specific atom accessibilities provided by ProtSA will allow obtaining better estimates of the contribution of the hydrophobic effect to the free energy of folding, will help to refine existing parameterizations of protein folding energetics, and will be useful to understand the influence of point mutations on protein stability.

Surface location of individual residues of SlpA provides insight into the Lactobacillus brevis S-layer

Bacterial surface layer (S-layer) proteins are excellent candidates for in vivo and in vitro nanobiotechnological applications because of their ability to self-assemble into two-dimensional lattices that form the outermost layer of many Eubacteria and most Archaea species. Despite this potential, knowledge about their molecular architecture is limited. In this study, we investigated SlpA, the S-layer protein of the potentially probiotic bacterium Lactobacillus brevis ATCC 8287 by cysteine-scanning mutagenesis and chemical modification. We generated a series of 46 mutant proteins by replacing single amino acids with cysteine, which is not present in the wild-type protein. Most of the replaced amino acids were located in the self-assembly domain (residues 179 to 435) that likely faces the outer surface of the lattice. As revealed by electron microscopy, all the mutant proteins were able to form self-assembly products identical to that of the wild type, proving that this replacement does not dramatically alter the protein conformation. The surface accessibility of the sulfhydryl groups introduced was studied with two maleimide-containing marker molecules, TMM(PEG)(12) (molecular weight [MW], 2,360) and AlexaFluor488-maleimide (MW = 720), using both monomeric proteins in solution and proteins allowed to self-assemble on cell wall fragments. Using the acquired data and available domain information, we assigned the mutated residues into four groups according to their location in the protein monomer and lattice structure: outer surface of the lattice (9 residues), inner surface of the lattice (9), protein interior (12), and protein-protein interface/pore regions (16). This information is essential, e.g., in the development of therapeutic and other health-related applications of Lactobacillus S-layers.

Putative re-entrant loop 1 of AE2 transmembrane domain has a major role in acute regulation of anion exchange by pH

Normal pH sensitivity of the SLC4A2/AE2 anion exchanger requires transmembrane domain (TMD) amino acid (aa) residues not conserved in the homologous but relatively pH-insensitive SLC4A1/AE1 polypeptide. We tested the hypothesis that the nonconserved aa cluster 1075DKPK1078 within the first putative re-entrant loop (RL1) of AE2 TMD contributes to pH sensor function by studying anion exchange function of AE2 mutants in which these and other RL1 aa were systematically substituted with corresponding RL1 aa from AE1. Regulation of Cl-/Cl- and Cl-/HCO(-)3 exchange by intracellular pH (pHi) or extracellular pH (pHo) was measured as 4,4'-di-isothiocyanatostilbene-2,2' disulfonic acid-sensitive 36Cl- efflux from Xenopus oocytes. AE2 RL1 mutants 1075AAAQ1078 and 1075AAAQN1079 showed reduced pHi sensitivity and pHo sensitivity was acid-shifted by approximately 1 pH unit. Individual mutants D1075A and P1077A exhibited moderately altered pH sensitivity, whereas a range of substitutions at conserved AE2 Ile-1079 substantially altered sensitivity to pHo and/or pHi. Substitution of the complete AE1 RL1 with AE2 RL1 failed to confer AE2-like pH sensitivity onto AE1. Replacement, however, of AE1 RL1 763SGPGAAAQ770 with AE2 1071VAPGDKPK1078 restored pHi sensitivity to the chimera AE2(1-920)/AE1(613-929) without affecting its low sensitivity to pHo. The results show that acute regulation of AE2 by pH requires RL1 of the TMD. We propose that critical segments of RL1 constitute part of an AE2 pH sensor that, together with residues within the N-terminal half of the TMD, constrain the AE2 polypeptide in a conformation required for regulation of anion exchange by pHi.

Joint evolutionary trees: a large-scale method to predict protein interfaces based on sequence sampling

The Joint Evolutionary Trees (JET) method detects protein interfaces, the core residues involved in the folding process, and residues susceptible to site-directed mutagenesis and relevant to molecular recognition. The approach, based on the Evolutionary Trace (ET) method, introduces a novel way to treat evolutionary information. Families of homologous sequences are analyzed through a Gibbs-like sampling of distance trees to reduce effects of erroneous multiple alignment and impacts of weakly homologous sequences on distance tree construction. The sampling method makes sequence analysis more sensitive to functional and structural importance of individual residues by avoiding effects of the overrepresentation of highly homologous sequences and improves computational efficiency. A carefully designed clustering method is parametrized on the target structure to detect and extend patches on protein surfaces into predicted interaction sites. Clustering takes into account residues' physical-chemical properties as well as conservation. Large-scale application of JET requires the system to be adjustable for different datasets and to guarantee predictions even if the signal is low. Flexibility was achieved by a careful treatment of the number of retrieved sequences, the amino acid distance between sequences, and the selective thresholds for cluster identification. An iterative version of JET (iJET) that guarantees finding the most likely interface residues is proposed as the appropriate tool for large-scale predictions. Tests are carried out on the Huang database of 62 heterodimer, homodimer, and transient complexes and on 265 interfaces belonging to signal transduction proteins, enzymes, inhibitors, antibodies, antigens, and others. A specific set of proteins chosen for their special functional and structural properties illustrate JET behavior on a large variety of interactions covering proteins, ligands, DNA, and RNA. JET is compared at a large scale to ET and to Consurf, Rate4Site, siteFiNDER|3D, and SCORECONS on specific structures. A significant improvement in performance and computational efficiency is shown.

Information obtained on the structure of macromolecular complexes is important for identifying functionally important partners but also for determining how such interactions will be perturbed by natural or engineered site mutations. Hence, to fully understand or control biological processes we need to predict in the most accurate manner protein interfaces for a protein structure, possibly without knowing its partners. Joint Evolutionary Trees (JET) is a method designed to detect very different types of interactions of a protein with another protein, ligands, DNA, and RNA. It uses a carefully designed sampling method, making sequence analysis more sensitive to the functional and structural importance of individual residues, and a clustering method parametrized on the target structure for the detection of patches on protein surfaces and their extension into predicted interaction sites. JET is a large-scale method, highly accurate and potentially applicable to search for protein partners.

A solution NMR investigation into the early events of amelogenin nanosphere self-assembly initiated with sodium chloride or calcium chloride

Using solution-state NMR spectroscopy, new insights into the early events governing amelogenin supramolecular self-assembly have been identified using sodium chloride and calcium chloride to trigger the association. Two-dimensional 1H-15N HSQC spectra were recorded for 15N- and 13C-labeled murine amelogenin as a function of increasing NaCl and CaCl2 concentration beginning with solution conditions of 2% acetic acid at pH 3.0, where amelogenin was monomeric. Residue specific changes in molecular dynamics, manifested by the reduction in intensity and disappearance of 1H-15N HSQC cross-peaks, were observed with the addition of either salt to the protein. With increasing NaCl concentrations, residues between T21 and R31 near the N-terminus were affected first, suggesting that these residues may initiate amelogenin dimerization, the first step in nanosphere assembly. At higher NaCl concentrations, more residues near the N-terminus (Y12-I51) were affected, and with further additions of NaCl, residues near the C-terminus (L141-T171) began to show a similar change in molecular dynamics. With increasing CaCl2 concentrations, a similar stepwise change in molecular dynamics involving essentially the same set of amelogenin residues was observed. As the concentration of either salt was increased, a concomitant increase in the estimated overall rotational correlation time (tau(c)) was observed, consistent with assembly. Self-assembly into a dimer or trimer was established with dynamic light scattering studies under similar conditions that showed an increase in diameter of the smallest species from 4.1 nm in the absence of salt to 10 nm in the presence of salt. These results suggest a possible stepwise interaction mechanism, starting with the N-terminus and followed by the C-terminus, leading to amelogenin nanosphere assembly.

Methods of calculating protein hydrophobicity and their application in developing correlations to predict hydrophobic interaction chromatography retention

Hydrophobic interaction chromatography (HIC) is a key technique for protein separation and purification. Different methodologies to estimate the hydrophobicity of a protein are reviewed, which have been related to the chromatographic behavior of proteins in HIC. These methodologies consider either knowledge of the three-dimensional structure or the amino acid composition of proteins. Despite some restrictions; they have proven to be useful in predicting protein retention time in HIC.

Evolutionary modeling of rate shifts reveals specificity determinants in HIV-1 subtypes

A hallmark of the human immunodeficiency virus 1 (HIV-1) is its rapid rate of evolution within and among its various subtypes. Two complementary hypotheses are suggested to explain the sequence variability among HIV-1 subtypes. The first suggests that the functional constraints at each site remain the same across all subtypes, and the differences among subtypes are a direct reflection of random substitutions, which have occurred during the time elapsed since their divergence. The alternative hypothesis suggests that the functional constraints themselves have evolved, and thus sequence differences among subtypes in some sites reflect shifts in function. To determine the contribution of each of these two alternatives to HIV-1 subtype evolution, we have developed a novel Bayesian method for testing and detecting site-specific rate shifts. The RAte Shift EstimatoR (RASER) method determines whether or not site-specific functional shifts characterize the evolution of a protein and, if so, points to the specific sites and lineages in which these shifts have most likely occurred. Applying RASER to a dataset composed of large samples of HIV-1 sequences from different group M subtypes, we reveal rampant evolutionary shifts throughout the HIV-1 proteome. Most of these rate shifts have occurred during the divergence of the major subtypes, establishing that subtype divergence occurred together with functional diversification. We report further evidence for the emergence of a new sub-subtype, characterized by abundant rate-shifting sites. When focusing on the rate-shifting sites detected, we find that many are associated with known function relating to viral life cycle and drug resistance. Finally, we discuss mechanisms of covariation of rate-shifting sites.

The AIDS epidemic, inflicted by the human immunodeficiency virus (HIV), has already claimed 25 million lives, thus posing a global threat. Since its discovery, several HIV subtypes have emerged, characterized by distinct genomic sequences and variable geographic locations. Here, we investigate the nature of the genetic differences among the subtypes. The neutral theory of evolution suggests that most genetic differences marginally affect the function of the encoded proteins (hence neutral) and thus occur randomly. Alternatively, changes in protein function are reflected by a pattern of nonrandom genetic differences. To address this issue, we developed a computational method, which studies the differences between sequences of different HIV subtypes, and estimates which of the explanations is more likely. Using a large sample of HIV protein sequences, we discovered that part of the variability among the subtypes is not random and possibly reflects different functional constraints imposed on the subtypes during the course of their evolution. An in-depth inspection of these nonrandom changes revealed a correlation with biological traits, such as drug resistance and mechanisms facilitating viral entry into the host cell. Interestingly, nonrandom changes are also characteristic of a viral strain that recently emerged in the former Soviet Union.