Value of eight-amino-acid matches in predicting the allergenicity status of proteins: an empirical bioinformatic investigation

The use of biotechnological techniques to introduce novel proteins into food crops (transgenic or GM crops) has motivated investigation into the properties of proteins that favor their potential to elicit allergic reactions. As part of the allergenicity assessment, bioinformatic approaches are used to compare the amino-acid sequence of candidate proteins with sequences in a database of known allergens to predict potential cross reactivity between novel food proteins and proteins to which people have become sensitized. Two criteria commonly used for these queries are searches over 80-amino-acid stretches for >35% identity, and searches for 8-amino-acid contiguous matches. We investigated the added value provided by the 8-amino-acid criterion over that provided by the >35%-identity-over-80-amino-acid criterion, by identifying allergens pairs that only met the former criterion, but not the latter criterion. We found that the allergen-sequence pairs only sharing 8-amino-acid identity, but not >35% identity over 80 amino acids, were unlikely to be cross reactive allergens. Thus, the common search for 8-amino-acid identity between novel proteins and known allergens appears to be of little additional value in assessing the potential allergenicity of novel proteins

Interconversion of Functional Motions between Mesophilic and Thermophilic Adenylate Kinases

Dynamic properties are functionally important in many proteins, including the enzyme adenylate kinase (AK), for which the open/closed transition limits the rate of catalytic turnover. Here, we compare our previously published coarse-grained (double-well Gō) simulation of mesophilic AK from E. coli (AKmeso) to simulations of thermophilic AK from Aquifex aeolicus (AKthermo). In AKthermo, as with AKmeso, the LID domain prefers to close before the NMP domain in the presence of ligand, but LID rigid-body flexibility in the open (O) ensemble decreases significantly. Backbone foldedness in O and/or transition state (TS) ensembles increases significantly relative to AKmeso in some interdomain backbone hinges and within LID. In contact space, the TS of AKthermo has fewer contacts at the CORE-LID interface but a stronger contact network surrounding the CORE-NMP interface than the TS of AKmeso. A “heated” simulation of AKthermo at 375K slightly increases LID rigid-body flexibility in accordance with the “corresponding states” hypothesis. Furthermore, while computational mutation of 7 prolines in AKthermo to their AKmeso counterparts produces similar small perturbations, mutation of these sites, especially positions 8 and 155, to glycine is required to achieve LID rigid-body flexibility and hinge flexibilities comparable to AKmeso. Mutating the 7 sites to proline in AKmeso reduces some hinges' flexibilities, especially hinge 2, but does not reduce LID rigid-body flexibility, suggesting that these two types of motion are decoupled in AKmeso. In conclusion, our results suggest that hinge flexibility and global functional motions alike are correlated with but not exclusively determined by the hinge residues. This mutational framework can inform the rational design of functionally important flexibility and allostery in other proteins toward engineering novel biochemical pathways

Nonlinear elasticity, proteinquakes, and the energy landscapes of functional transitions in proteins

Large-scale motions of biomolecules involve linear elastic deformations along low-frequency normal modes, but for function nonlinearity is essential. In addition, unlike macroscopic machines, biological machines can locally break and then reassemble during function. We present a model for global structural transformations, such as allostery, that involve large-scale motion and possible partial unfolding, illustrating the method with the conformational transition of adenylate kinase. Structural deformation between open and closed states occurs via low-frequency modes on separate reactant and product surfaces, switching from one state to the other when energetically favorable. The switching model is the most straightforward anharmonic interpolation, which allows the barrier for a process to be estimated from a linear normal mode calculation, which by itself cannot be used for activated events. Local unfolding, or cracking, occurs in regions where the elastic stress becomes too high during the transition. Cracking leads to a counterintuitive catalytic effect of added denaturant on allosteric enzyme function. It also leads to unusual relationships between equilibrium constant and rate like those seen recently in single-molecule experiments of motor proteins