Experimental assessment of the importance of amino acid positions identified by an entropy-based correlation analysis of multiple-sequence alignments

A Single Mutation Increases the Thermostability and Activity of Aspergillus terreus Amine Transaminase

Enhancing the thermostability of (R)-selective amine transaminases (AT-ATA) will expand its application in the asymmetric synthesis of chiral amines. In this study, mutual information and coevolution networks of ATAs were analyzed by the Mutual Information Server to Infer Coevolution (MISTIC). Subsequently, the amino acids most likely to influence the stability and function of the protein were investigated by alanine scanning and saturation mutagenesis. Four stabilized mutants (L118T, L118A, L118I, and L118V) were successfully obtained. The best mutant, L118T, exhibited an improved thermal stability with a 3.7-fold enhancement in its half-life (t_1/2) at 40 °C and a 5.3 °C increase in T₅₀¹⁰ compared to the values for the wild-type protein. By the differential scanning fluorimetry (DSF) analysis, the best mutant, L118T, showed a melting temperature (T_m) of 46.4 °C, which corresponded to a 5.0 °C increase relative to the wild-type AT-ATA (41.4 °C). Furthermore, the most stable mutant L118T displayed the highest catalytic efficiency among the four stabilized mutants.

Correlated positions in protein evolution and engineering

Statistical analysis of a protein multiple sequence alignment can reveal groups of positions that undergo interdependent mutations throughout evolution. At these so-called correlated positions, only certain combinations of amino acids appear to be viable for maintaining proper folding, stability, catalytic activity or specificity. Therefore, it is often speculated that they could be interesting guides for semi-rational protein engineering purposes. Because they are a fingerprint from protein evolution, their analysis may provide valuable insight into a protein's structure or function and furthermore, they may also be suitable target positions for mutagenesis. Unfortunately, little is currently known about the properties of these correlation networks and how they should be used in practice. This review summarises the recent findings, opportunities and pitfalls of the concept.

Protein thermostability engineering

Using structure and sequence based analysis we can engineer proteins to increase their thermal stability.

Antibody Binding Selectivity: Alternative Sets of Antigen Residues Entail High-Affinity Recognition

Understanding the relationship between protein sequence and molecular recognition selectivity remains a major challenge. The antibody fragment scFv1F4 recognizes with sub nM affinity a decapeptide (sequence 6TAMFQDPQER15) derived from the N-terminal end of human papilloma virus E6 oncoprotein. Using this decapeptide as antigen, we had previously shown that only the wild type amino-acid or conservative replacements were allowed at positions 9 to 12 and 15 of the peptide, indicating a strong binding selectivity. Nevertheless phenylalanine (F) was equally well tolerated as the wild type glutamine (Q) at position 13, while all other amino acids led to weaker scFv binding. The interfaces of complexes involving either Q or F are expected to diverge, due to the different physico-chemistry of these residues. This would imply that high-affinity binding can be achieved through distinct interfacial geometries. In order to investigate this point, we disrupted the scFv-peptide interface by modifying one or several peptide positions. We then analyzed the effect on binding of amino acid changes at the remaining positions, an altered susceptibility being indicative of an altered role in complex formation. The 23 starting variants analyzed contained replacements whose effects on scFv1F4 binding ranged from minor to drastic. A permutation analysis (effect of replacing each peptide position by all other amino acids except cysteine) was carried out on the 23 variants using the PEPperCHIP® Platform technology. A comparison of their permutation patterns with that of the wild type peptide indicated that starting replacements at position 11, 12 or 13 modified the tolerance to amino-acid changes at the other two positions. The interdependence between the three positions was confirmed by SPR (Biacore® technology). Our data demonstrate that binding selectivity does not preclude the existence of alternative high-affinity recognition modes.

H2rs: deducing evolutionary and functionally important residue positions by means of an entropy and similarity based analysis of multiple sequence alignments

BACKGROUND

The identification of functionally important residue positions is an important task of computational biology. Methods of correlation analysis allow for the identification of pairs of residue positions, whose occupancy is mutually dependent due to constraints imposed by protein structure or function. A common measure assessing these dependencies is the mutual information, which is based on Shannon's information theory that utilizes probabilities only. Consequently, such approaches do not consider the similarity of residue pairs, which may degrade the algorithm's performance. One typical algorithm is H2r, which characterizes each individual residue position k by the conn(k)-value, which is the number of significantly correlated pairs it belongs to.

RESULTS

To improve specificity of H2r, we developed a revised algorithm, named H2rs, which is based on the von Neumann entropy (vNE). To compute the corresponding mutual information, a matrix A is required, which assesses the similarity of residue pairs. We determined A by deducing substitution frequencies from contacting residue pairs observed in the homologs of 35 809 proteins, whose structure is known. In analogy to H2r, the enhanced algorithm computes a normalized conn(k)-value. Within the framework of H2rs, only statistically significant vNE values were considered. To decide on significance, the algorithm calculates a p-value by performing a randomization test for each individual pair of residue positions. The analysis of a large in silico testbed demonstrated that specificity and precision were higher for H2rs than for H2r and two other methods of correlation analysis. The gain in prediction quality is further confirmed by a detailed assessment of five well-studied enzymes. The outcome of H2rs and of a method that predicts contacting residue positions (PSICOV) overlapped only marginally. H2rs can be downloaded from http://www-bioinf.uni-regensburg.de.

CONCLUSIONS

Considering substitution frequencies for residue pairs by means of the von Neumann entropy and a p-value improved the success rate in identifying important residue positions. The integration of proven statistical concepts and normalization allows for an easier comparison of results obtained with different proteins. Comparing the outcome of the local method H2rs and of the global method PSICOV indicates that such methods supplement each other and have different scopes of application.

Loop-loop interactions govern multiple steps in indole-3-glycerol phosphate synthase catalysis

Substrate binding, product release, and likely chemical catalysis in the tryptophan biosynthetic enzyme indole-3-glycerol phosphate synthase (IGPS) are dependent on the structural dynamics of the β1α1 active-site loop. Statistical coupling analysis and molecular dynamic simulations had previously indicated that covarying residues in the β1α1 and β2α2 loops, corresponding to Arg54 and Asn90, respectively, in the Sulfolobus sulfataricus enzyme (ssIGPS), are likely important for coordinating functional motions of these loops. To test this hypothesis, we characterized site mutants at these positions for changes in catalytic function, protein stability and structural dynamics for the thermophilic ssIGPS enzyme. Although there were only modest changes in the overall steady-state kinetic parameters, solvent viscosity and solvent deuterium kinetic isotope effects indicated that these amino acid substitutions change the identity of the rate-determining step across multiple temperatures. Surprisingly, the N90A substitution had a dramatic effect on the general acid/base catalysis of the dehydration step, as indicated by the loss of the descending limb in the pH rate profile, which we had previously assigned to Lys53 on the β1α1 loop. These changes in enzyme function are accompanied with a quenching of ps-ns and µs-ms timescale motions in the β1α1 loop as measured by nuclear magnetic resonance studies. Altogether, our studies provide structural, dynamic and functional rationales for the coevolution of residues on the β1α1 and β2α2 loops, and highlight the multiple roles that the β1α1 loop plays in IGPS catalysis. Thus, substitution of covarying residues in the active-site β1α1 and β2α2 loops of indole-3-glycerol phosphate synthase results in functional, structural, and dynamic changes, highlighting the multiple roles that the β1α1 loop plays in enzyme catalysis and the importance of regulating the structural dynamics of this loop through noncovalent interactions with nearby structural elements.