Bioinformatic analysis of subfamily-specific regions in 3D-structures of homologs to study functional diversity and conformational plasticity in protein superfamilies

The Role of a Loop in the Non-catalytic Domain B on the Hydrolysis/Transglycosylation Specificity of the 4-α-Glucanotransferase from Thermotoga maritima

The mechanism by which glycoside hydrolases control the reaction specificity through hydrolysis or transglycosylation is a key element embedded in their chemical structures. The determinants of reaction specificity seem to be complex. We looked for structural differences in domain B between the 4-α-glucanotransferase from Thermotoga maritima (TmGTase) and the α-amylase from Thermotoga petrophila (TpAmylase) and found a longer loop in the former that extends towards the active site carrying a W residue at its tip. Based on these differences we constructed the variants W131G and the partial deletion of the loop at residues 120-124/128-131, which showed a 11.6 and 11.4-fold increased hydrolysis/transglycosylation (H/T) ratio relative to WT protein, respectively. These variants had a reduction in the maximum velocity of the transglycosylation reaction, while their affinity for maltose as the acceptor was not substantially affected. Molecular dynamics simulations allow us to rationalize the increase in H/T ratio in terms of the flexibility near the active site and the conformations of the catalytic acid residues and their associated pKas.

Beyond sequence: Structure-based machine learning

Recent breakthroughs in protein structure prediction demarcate the start of a new era in structural bioinformatics. Combined with various advances in experimental structure determination and the uninterrupted pace at which new structures are published, this promises an age in which protein structure information is as prevalent and ubiquitous as sequence. Machine learning in protein bioinformatics has been dominated by sequence-based methods, but this is now changing to make use of the deluge of rich structural information as input. Machine learning methods making use of structures are scattered across literature and cover a number of different applications and scopes; while some try to address questions and tasks within a single protein family, others aim to capture characteristics across all available proteins. In this review, we look at the variety of structure-based machine learning approaches, how structures can be used as input, and typical applications of these approaches in protein biology. We also discuss current challenges and opportunities in this all-important and increasingly popular field.

Guide tree optimization with genetic algorithm to improve multiple protein 3D-structure alignment

MOTIVATION

With the increasing availability of 3D-data, the focus of comparative bioinformatic analysis is shifting from protein sequence alignments toward more content-rich 3D-alignments. This raises the need for new ways to improve the accuracy of 3D-superimposition.

RESULTS

We proposed guide tree optimization with genetic algorithm (GA) as a universal tool to improve the alignment quality of multiple protein 3D-structures systematically. As a proof of concept, we implemented the suggested GA-based approach in popular Matt and Caretta multiple protein 3D-structure alignment (M3DSA) algorithms, leading to a statistically significant improvement of the TM-score quality indicator by up to 220-1523% on 'SABmark Superfamilies' (in 49-77% of cases) and 'SABmark Twilight' (in 59-80% of cases) datasets. The observed improvement in collections of distant homologies highlights the potentials of GA to optimize 3D-alignments of diverse protein superfamilies as one plausible tool to study the structure-function relationship.

AVAILABILITY AND IMPLEMENTATION

The source codes of patched gaCaretta and gaMatt programs are available open-access at https://github.com/n-canter/gamaps.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Computer simulation of the Receptor-Ligand Interactions of Mannose Receptor CD206 in Comparison with the Lectin Concanavalin A Model

Computer modeling of complexation of mono- and oligosaccharide ligands with the main (fourth) carbohydrate-binding domain of the mannose receptor CD206 (CRD4), as well as with the model receptor concanavalin A (ConA), was carried out for the first time, using methods of molecular dynamics and neural network analysis. ConA was shown to be a relevant model of CD206 (CRD4) due to similarity of the structural organization of the binding sites and high correlation of the values of free energies of complexation between the literature data and computer modeling (r > 0.9). Role of the main factors affecting affinity of the ligand–receptor interactions is discussed: the number and nature of carbohydrate residues, presence of Me-group in the O1 position, type of the glycoside bond in dimannose. Complexation of ConA and CD206 with ligands is shown to be energetically caused by electrostatic interactions (E) of the charged residues (Asn, Asp, Arg) with oxygen and hydrogen atoms in carbohydrates; contributions of hydrophobic and van der Waals components is lower. Possibility of the additional stabilization of complexes due to the CH–π stacking interactions of Tyr with the Man plane is discussed. The role of calcium and manganese ions in binding ligands has been studied. The values of free energies of complexation calculated in the course of molecular dynamics simulation correlate with experimental data (published for the model ConA): correlation coefficient r = 0.68. The Pafnucy neural network was trained based on the set of PDBbind2020 ligand–receptor complexes, which significantly increased accuracy of the energy predictions to r = 0.8 and 0.82 for CD206 and ConA receptors, respectively. A model of normalization of the complexation energy values for calculating the relevant values of ΔG_bind, K_d is proposed. Based on the developed technique, values of the dissociation constants of a series of CD206 complexes with nine carbohydrate ligands of different structures were determined, which were not previously known. The obtained data open up possibilities for using computer modeling for the development of optimal drug carriers capable of active macrophage targeting, and also determine the limits of applicability of using ConA as a relevant model for studying parameters of the CD206 binding to various carbohydrate ligands.

Single amino acid residue mediates reciprocal specificity in two mosquito odorant receptors

The southern house mosquito, Culex quinquefasciatus, utilizes two odorant receptors, CquiOR10 and CquiOR2, narrowly tuned to oviposition attractants and well conserved among mosquito species. They detect skatole and indole, respectively, with reciprocal specificity. We swapped the transmembrane (TM) domains of CquiOR10 and CquiOR2 and identified TM2 as a specificity determinant. With additional mutations, we showed that CquiOR10A73L behaved like CquiOR2. Conversely, CquiOR2L74A recapitulated CquiOR10 specificity. Next, we generated structural models of CquiOR10 and CquiOR10A73L using RoseTTAFold and AlphaFold and docked skatole and indole using RosettaLigand. These modeling studies suggested space-filling constraints around A73. Consistent with this hypothesis, CquiOR10 mutants with a bulkier residue (Ile, Val) were insensitive to skatole and indole, whereas CquiOR10A73G retained the specificity to skatole and showed a more robust response than the wildtype receptor CquiOR10. On the other hand, Leu to Gly mutation of the indole receptor CquiOR2 reverted the specificity to skatole. Lastly, CquiOR10A73L, CquiOR2, and CquiOR2L74I were insensitive to 3-ethylindole, whereas CquiOR2L74A and CquiOR2L74G gained activity. Additionally, CquiOR10A73G gave more robust responses to 3-ethylindole than CquiOR10. Thus, we suggest the specificity of these receptors is mediated by a single amino acid substitution, leading to finely tuned volumetric space to accommodate specific oviposition attractants.

Modulating Glycoside Hydrolase Activity between Hydrolysis and Transfer Reactions Using an Evolutionary Approach

The proteins within the CAZy glycoside hydrolase family GH13 catalyze the hydrolysis of polysaccharides such as glycogen and starch. Many of these enzymes also perform transglycosylation in various degrees, ranging from secondary to predominant reactions. Identifying structural determinants associated with GH13 family reaction specificity is key to modifying and designing enzymes with increased specificity towards individual reactions for further applications in industrial, chemical, or biomedical fields. This work proposes a computational approach for decoding the determinant structural composition defining the reaction specificity. This method is based on the conservation of coevolving residues in spatial contacts associated with reaction specificity. To evaluate the algorithm, mutants of α-amylase (TmAmyA) and glucanotransferase (TmGTase) from Thermotoga maritima were constructed to modify the reaction specificity. The K98P/D99A/H222Q variant from TmAmyA doubled the transglycosydation/hydrolysis (T/H) ratio while the M279N variant from TmGTase increased the hydrolysis/transglycosidation ratio five-fold. Molecular dynamic simulations of the variants indicated changes in flexibility that can account for the modified T/H ratio. An essential contribution of the presented computational approach is its capacity to identify residues outside of the active center that affect the reaction specificity.

Computational approaches to predict protein functional families and functional sites

Understanding the mechanisms of protein function is indispensable for many biological applications, such as protein engineering and drug design. However, experimental annotations are sparse, and therefore, theoretical strategies are needed to fill the gap. Here, we present the latest developments in building functional subclassifications of protein superfamilies and using evolutionary conservation to detect functional determinants, for example, catalytic-, binding- and specificity-determining residues important for delineating the functional families. We also briefly review other features exploited for functional site detection and new machine learning strategies for combining multiple features.