Search

Arrival of the novel SARS-CoV-2 has launched a worldwide effort to identify both pre-approved and novel therapeutics targeting the viral proteome, highlighting the urgent need for efficient drug discovery strategies. Even with effective vaccines, infection is possible, and at-risk populations would benefit from effective drug compounds that reduce the lethality and lasting damage of COVID-19 infection. The CoV-2 MacroD-like macrodomain (Mac1) is implicated in viral pathogenicity by disrupting host innate immunity through its mono (ADP-ribosyl) hydrolase activity, making it a prime target for antiviral therapy. We therefore solved the structure of CoV-2 Mac1 from non-structural protein 3 (Nsp3) and applied structural and sequence-based genetic tracing, including newly determined A. pompejana MacroD2 and GDAP2 amino acid sequences, to compare and contrast CoV-2 Mac1 with the functionally related human DNA-damage signaling factor poly (ADP-ribose) glycohydrolase (PARG). Previously, identified targetable features of the PARG active site allowed us to develop a pharmacologically useful PARG inhibitor (PARGi). Here, we developed a focused chemical library and determined 6 novel PARGi X-ray crystal structures for comparative analysis. We applied this knowledge to discovery of CoV-2 Mac1 inhibitors by combining computation and structural analysis to identify PARGi fragments with potential to bind the distal-ribose and adenosyl pockets of the CoV-2 Mac1 active site. Scaffold development of these PARGi fragments has yielded two novel compounds, PARG-345 and PARG-329, that crystallize within the Mac1 active site, providing critical structure-activity data and a pathway for inhibitor optimization. The reported structural findings demonstrate ways to harness our PARGi synthesis and characterization pipeline to develop CoV-2 Mac1 inhibitors targeting the ADP-ribose active site. Together, these structural and computational analyses reveal a path for accelerating development of antiviral therapeutics from pre-existing drug optimization pipelines.

Abstract Background Along with the fast development and urbanization in developing countries, the waterbodies aside the growing cities become heavily polluted and highly eutrophic, thus leading to the seasonal outbreak of cyanobacterial bloom. Systematic isolation and characterization of freshwater cyanophages might provide a biological solution to control the awful blooms. However, genomic sequences and related investigations on the freshwater cyanophages remain very limited to date. Results Following our recently reported five cyanophages Pam1~Pam5 from Lake Chaohu in China, here we isolated another five cyanophages, termed Pan1~Pan5, which infect the cyanobacterium Pseudanabaena sp. Chao 1811. Whole-genome sequencing showed that they all contain a double-stranded DNA genome of 37.2 to 72.0 kb in length, with less than half of the putative open reading frames annotated with known functions. Remarkably, the siphophage Pan1 encodes an auxiliary metabolic gene phoH and constitutes, together with the host, a complete queuosine modification pathway. Proteomic analyses revealed that although Pan1~Pan5 are distinct from each other in evolution, Pan1 and Pan3 are somewhat similar to our previously identified cyanophages Pam3 and Pam1 at the genomic level, respectively. Moreover, phylogenetic analyses suggested that Pan1 resembles the α-proteobacterial phage vB_DshS-R5C, revealing direct evidence for phage-mediated horizontal gene transfer between cyanobacteria and α-proteobacteria. Conclusion In addition to the previous reports of Pam1~Pam5, the present findings on Pan1~Pan5 largely enrich the library of reference freshwater cyanophages. The abundant genomic information provides a pool to identify novel genes and proteins of unknown function. Moreover, we found for the first time the evolutionary traces in the cyanophage that horizontal gene transfer might occur at the level of not only inter-species, but even inter-phylum. It indicates that the bacteriophage or cyanophage could be developed as a powerful tool for gene manipulation among various species or phyla.

Background: Along with the fast development and urbanization in developing countries, the waterbodies aside the growing cities become heavily polluted and highly eutrophic, thus leading to the seasonal outbreak of cyanobacterial bloom. Systematic isolation and characterization of freshwater cyanophages might provide a biological solution to control the awful blooms. However, genomic sequences and related investigations on the freshwater cyanophages remain very limited to date. Results: Following our recently reported five cyanophages Pam1~Pam5 from Lake Chaohu in China, here we isolated another five cyanophages, termed Pan1~Pan5, which infect the cyanobacterium Pseudanabaena sp. Chao 1811. Whole-genome sequencing showed that they all contain a double-stranded DNA genome of 37.2 to 72.0 kb in length, with less than half of the putative open reading frames annotated with known functions. Remarkably, the siphophage Pan1 encodes an auxiliary metabolic gene phoH and constitutes, together with the host, a complete queuosine modification pathway. Proteomic analyses revealed that although Pan1~Pan5 are distinct from each other in evolution, Pan1 and Pan3 are somewhat similar to our previously identified cyanophages Pam3 and Pam1 at the genomic level, respectively. Moreover, phylogenetic analyses suggested that Pan1 resembles the α-proteobacterial phage vB_DshS-R5C, revealing direct evidence for phage-mediated horizontal gene transfer between cyanobacteria and α-proteobacteria. Conclusion: In addition to the previous reports of Pam1~Pam5, the present findings on Pan1~Pan5 largely enrich the library of reference freshwater cyanophages. The abundant genomic information provides a pool to identify novel genes and proteins of unknown function. Moreover, we found for the first time the evolutionary traces in the cyanophage that horizontal gene transfer might occur at the level of not only inter-species, but even inter-phylum. It indicates that the bacteriophage or cyanophage could be developed as a powerful tool for gene manipulation among various species or phyla. [ABSTRACT FROM AUTHOR]

Our understanding of enzymes with high substrate ambiguity remains limited because their large active sites allow substrate docking freedom to an extent that seems incompatible with stereospecificity. One possibility is that some of these enzymes evolved a set of evolutionarily fitted sequence positions that stringently allow switching substrate ambiguity and chiral specificity. To explore this hypothesis, we targeted for mutation a serine ester hydrolase (EH3) that exhibits an impressive 71-substrate repertoire but is not stereospecific (e.e. 50%). We used structural actions and the computational evolutionary trace method to explore specificity-swapping sequence positions and hypothesized that position I244 was critical. Driven by evolutionary action analysis, this position was substituted to leucine, which together with isoleucine appears to be the amino acid most commonly present in the closest homologous sequences (max. identity, ca. 67.1%), and to phenylalanine, which appears in distant homologues. While the I244L mutation did not have any functional consequences, the I244F mutation allowed the esterase to maintain a remarkable 53-substrate range while gaining stereospecificity properties (e.e. 99.99%). These data support the possibility that some enzymes evolve sequence positions that control the substrate scope and stereospecificity. Such residues, which can be evolutionarily screened, may serve as starting points for further designing substrate-ambiguous, yet chiral-specific, enzymes that are greatly appreciated in biotechnology and synthetic chemistry.

Aim: To reconstruct the ancestral sequence of human adenoviral hexon protein by combining sequence variations and structural information. And to provide a candidate hexon protein for developing new adenoviral vector capable of escaping the pre-existing immunity in healthy populations. Methods: The sequences of 74 adenovirus-type strains were used to predict the ancestral sequence of human adenovirus hexon protein using FastML and MEGA software. The three-dimensional structure model was built using homology modeling methods. The immunological features of ancestral loop 1 and loop 2 regions of sequences were tested using protein segments expressed in a prokaryotic expression system and polypeptides synthesized with human serum samples. Results: The tower region of the hexon protein had the highest sequence variability, while the neck and base regions remained constant among different types. The modern strains successfully predicted the common ancestral sequence of the human adenovirus hexon. The positive sera against neutralizing epitopes on the common ancestor of adenoviral hexon were relatively rare among healthy adults. Conclusion: The existing strains inferred the common ancestor of human adenoviruses, with epitopes never observed in the current human strains. The predicted common ancestor hexon is a good prospect in the improvement of adenovirus vectors. [ABSTRACT FROM AUTHOR]

Aim: To reconstruct the ancestral sequence of human adenoviral hexon protein by combining sequence variations and structural information. And to provide a candidate hexon protein for developing new adenoviral vector capable of escaping the pre-existing immunity in healthy populations.Methods: The sequences of 74 adenovirus-type strains were used to predict the ancestral sequence of human adenovirus hexon protein using FastML and MEGA software. The three-dimensional structure model was built using homology modeling methods. The immunological features of ancestral loop 1 and loop 2 regions of sequences were tested using protein segments expressed in a prokaryotic expression system and polypeptides synthesized with human serum samples.Results: The tower region of the hexon protein had the highest sequence variability, while the neck and base regions remained constant among different types. The modern strains successfully predicted the common ancestral sequence of the human adenovirus hexon. The positive sera against neutralizing epitopes on the common ancestor of adenoviral hexon were relatively rare among healthy adults.Conclusion: The existing strains inferred the common ancestor of human adenoviruses, with epitopes never observed in the current human strains. The predicted common ancestor hexon is a good prospect in the improvement of adenovirus vectors.

Group I chaperonins are a highly conserved family of essential proteins that self-assemble into molecular nanoboxes that mediate the folding of cytoplasmic proteins in bacteria and organelles. GroEL, the chaperonin of Escherichia coli, is the archetype of the family. Protein folding-independent functions have been described for numerous chaperonins, including HtpB, the chaperonin of the bacterial pathogen Legionella pneumophila. Several protein folding-independent functions attributed to HtpB are not shared by GroEL, suggesting that differences in the amino acid (aa) sequence between these two proteins could correlate with functional differences. GroEL and HtpB differ in 137 scattered aa positions. Using the Evolutionary Trace (ET) bioinformatics method, site-directed mutagenesis, and a functional reporter test based upon a yeast-two-hybrid interaction with the eukaryotic protein ECM29, it was determined that out of those 137 aa, ten (M68, M212, S236, K298, N507 and the cluster AEHKD in positions 471-475) were involved in the interaction of HtpB with ECM29. GroEL was completely unable to interact with ECM29, but when GroEL was modified at those 10 aa positions, to display the HtpB aa, it acquired a weak ability to interact with ECM29. This constitutes proof of concept that the unique functional abilities of HtpB can be mapped to specific aa positions.

Many computational approaches estimate the effect of coding variants, but their predictions often disagree with each other. These contradictions confound users and raise questions regarding reliability. Performance assessments can indicate the expected accuracy for each method and highlight advantages and limitations. The Critical Assessment of Genome Interpretation (CAGI) community aims to organize objective and systematic assessments: They challenge predictors on unpublished experimental and clinical data and assign independent assessors to evaluate the submissions. We participated in CAGI experiments as predictors, using the Evolutionary Action (EA) method to estimate the fitness effect of coding mutations. EA is untrained, uses homology information, and relies on a formal equation: The fitness effect equals the functional sensitivity to residue changes multiplied by the magnitude of the substitution. In previous CAGI experiments (between 2011 and 2016), our submissions aimed to predict the protein activity of single mutants. In 2018 (CAGI5), we also submitted predictions regarding clinical associations, folding stability, and matching genomic data with phenotype. For all these diverse challenges, we used EA to predict the fitness effect of variants, adjusted to specifically address each question. Our submissions had consistently good performance, suggesting that EA predicts reliably the effects of genetic variants. [ABSTRACT FROM AUTHOR]

β-propeller domains composed of WD repeats are highly ubiquitous and typically used as multi-site docking platforms to coordinate and integrate the activities of groups of proteins. Here, we have used extensive homology modelling of the WD40-repeat family of seven-bladed β-propellers coupled with subsequent structural classification and clustering of these models to define subfamilies of β-propellers with common structural, and probable, functional characteristics. We show that it is possible to assign seven-bladed WD β-propeller proteins into functionally different groups based on the information gained from homology modelling. We examine general structural diversity within the WD40-repeat family of seven-bladed β-propellers and demonstrate that seven-bladed β-propellers composed of WD-repeats are structurally distinct from other seven-bladed β-propellers. We further provide some insights into the multifunctional diversity of the seven-bladed WD β-propeller surfaces. This report once again reinforces the importance of structural data and the usefulness of homology models in functional classification.

The structural basis of allosteric signaling in G protein-coupled receptors (GPCRs) is important in guiding design of therapeutics and understanding phenotypic consequences of genetic variation. The Evolutionary Trace (ET) algorithm previously proved effective in redesigning receptors to mimic the ligand specificities of functionally distinct homologs. We now expand ET to consider mutual information, with validation in GPCR structure and dopamine D2 receptor (D2R) function. The new algorithm, called ET-MIp, identifies evolutionarily relevant patterns of amino acid covariations. The improved predictions of structural proximity and D2R mutagenesis demonstrate that ET-MIp predicts functional interactions between residue pairs, particularly potency and efficacy of activation by dopamine. Remarkably, although most of the residue pairs chosen for mutagenesis are neither in the binding pocket nor in contact with each other, many exhibited functional interactions, implying at-a-distance coupling. The functional interaction between the coupled pairs correlated best with the evolutionary coupling potential derived from dopamine receptor sequences rather than with broader sets of GPCR sequences. These data suggest that the allosteric communication responsible for dopamine responses is resolved by ET-MIp and best discerned within a short evolutionary distance.Most double mutants restored dopamine response to wild-type levels, also suggesting that tight regulation of the response to dopamine drove the coevolution and intramolecular communications between coupled residues. Our approach provides a general tool to identify evolutionary covariation patterns in small sets of close sequence homologs and to translate them into functional linkages between residues. [ABSTRACT FROM AUTHOR]

ABSTRACT In Pseudomonas aeruginosa, the algH gene regulates the cellular concentrations of a number of enzymes and the production of several virulence factors, and is suggested to serve a global regulatory function. The precise mechanism by which the algH gene product, the AlgH protein, functions is unknown. The same is true for AlgH family members from other bacteria. In order to lay the groundwork for understanding the physical underpinnings of AlgH function, we examined the structure and physical properties of AlgH in solution. Under reducing conditions, results of NMR, electrophoretic mobility, and sedimentation equilibrium experiments indicate AlgH is predominantly monomeric and monodisperse in solution. Under nonreducing conditions intra and intermolecular disulfide bonds form, the latter promoting AlgH oligomerization. The high-resolution solution structure of AlgH reveals alpha/beta-sandwich architecture fashioned from ten beta strands and seven alpha helices. Comparison with available structures of orthologues indicates conservation of overall structural topology. The region of the protein most strongly conserved structurally also shows the highest amino acid sequence conservation and, as revealed by hydrogen-deuterium exchange studies, is also the most stable. In this region, evolutionary trace analysis identifies two clusters of amino acid residues with the highest evolutionary importance relative to all other AlgH residues. These frame a partially solvent exposed shallow hydrophobic cleft, perhaps identifying a site for intermolecular interactions. The results establish a physical foundation for understanding the structure and function of AlgH and AlgH family proteins and should be of general importance for further investigations of these and related proteins. Proteins 2015; 83:1137-1150. © 2015 Wiley Periodicals, Inc. [ABSTRACT FROM AUTHOR]

Understanding the molecular basis of protein function remains a central goal of biology, with the hope to elucidate the role of human genes in health and in disease, and to rationally design therapies through targeted molecular perturbations. We review here some of the computational techniques and resources available for characterizing a critical aspect of protein function – those mediated by protein–protein interactions (PPI). We describe several applications and recent successes of the Evolutionary Trace (ET) in identifying molecular events and shapes that underlie protein function and specificity in both eukaryotes and prokaryotes. ET is a part of analytical approaches based on the successes and failures of evolution that enable the rational control of PPI. [ABSTRACT FROM AUTHOR]

Abstract: Leprosy remains a major public health problem, since single and multi-drug resistance has been reported worldwide over the last two decades. In the present study, we report the novel multi-targeted therapy for leprosy to overcome multi drug resistance and to improve therapeutic efficacy. If multiple enzymes of an essential metabolic pathway of a bacterium were targeted, then the therapy would become more effective and can prevent the occurrence of drug resistance. The MurC, MurD, MurE and MurF enzymes of peptidoglycan biosynthetic pathway were selected for multi targeted therapy. The conserved or class specific active site residues important for function or stability were predicted using evolutionary trace analysis and site directed mutagenesis studies. Ten such residues which were present in at least any three of the four Mur enzymes (MurC, MurD, MurE and MurF) were identified. Among the ten residues G125, K126, T127 and G293 (numbered based on their position in MurC) were found to be conserved in all the four Mur enzymes of the entire bacterial kingdom. In addition K143, T144, T166, G168, H234 and Y329 (numbered based on their position in MurE) were significant in binding substrates and/co-factors needed for the functional events in any three of the Mur enzymes. These are the probable residues for designing newer anti-leprosy drugs in an attempt to reduce drug resistance. [Copyright &y& Elsevier]

Abstract: Seven transmembrane (7TM) or G protein-coupled receptors constitute a large superfamily of cell surface receptors sharing a structural motif of seven transmembrane spanning alpha helices. Their activation mechanism most likely involves concerted movements of the transmembrane helices, but remains to be completely resolved. Evolutionary Trace (ET) analysis is a computational method, which identifies clusters of functionally important residues by integrating information on evolutionary important residue variations with receptor structure. Combined with known mutational data, ET predicted a patch of residues in the cytoplasmic parts of TM2, TM3, and TM6 to form an activation switch that is common to all family A 7TM receptors. We tested this hypothesis in the rat Angiotensin II (Ang II) type 1a (AT1a) receptor. The receptor has important roles in the cardiovascular system, but has also frequently been applied as a model for 7TM receptor activation and signaling. Six mutations: F66A, L67R, L70R, L119R, D125A, and I245F were targeted to the putative switch and assayed for changes in activation state by their ligand binding, signaling, and trafficking properties. All but one receptor mutant (that was not expressed well) displayed phenotypes associated with changed activation state, such as increased agonist affinity or basal activity, promiscuous activation, or constitutive internalization highlighting the importance of testing different signaling pathways. We conclude that this evolutionary important patch mediates interactions important for maintaining the inactive state. More broadly, these observations in the AT1 receptor are consistent with computational predictions of a generic role for this patch in 7TM receptor activation. [Copyright &y& Elsevier]

Group I chaperonins are a highly conserved family of essential proteins that self-assemble into molecular nanoboxes that mediate the folding of cytoplasmic proteins in bacteria and organelles. GroEL, the chaperonin of Escherichia coli, is the archetype of the family. Protein folding-independent functions have been described for numerous chaperonins, including HtpB, the chaperonin of the bacterial pathogen Legionella pneumophila. Several protein folding-independent functions attributed to HtpB are not shared by GroEL, suggesting that differences in the amino acid (aa) sequence between these two proteins could correlate with functional differences. GroEL and HtpB differ in 137 scattered aa positions. Using the Evolutionary Trace (ET) bioinformatics method, site-directed mutagenesis, and a functional reporter test based upon a yeast-two-hybrid interaction with the eukaryotic protein ECM29, it was determined that out of those 137 aa, ten (M68, M212, S236, K298, N507 and the cluster AEHKD in positions 471-475) were involved in the interaction of HtpB with ECM29. GroEL was completely unable to interact with ECM29, but when GroEL was modified at those 10 aa positions, to display the HtpB aa, it acquired a weak ability to interact with ECM29. This constitutes proof of concept that the unique functional abilities of HtpB can be mapped to specific aa positions. [ABSTRACT FROM AUTHOR]

β-propeller domains composed of WD repeats are highly ubiquitous and typically used as multi-site docking platforms to coordinate and integrate the activities of groups of proteins. Here, we have used extensive homology modelling of the WD40-repeat family of seven-bladed β-propellers coupled with subsequent structural classification and clustering of these models to define subfamilies of β-propellers with common structural, and probable, functional characteristics. We show that it is possible to assign seven-bladed WD β-propeller proteins into functionally different groups based on the information gained from homology modelling. We examine general structural diversity within the WD40-repeat family of seven-bladed β-propellers and demonstrate that seven-bladed β-propellers composed of WD-repeats are structurally distinct from other seven-bladed β-propellers. We further provide some insights into the multifunctional diversity of the seven-bladed WD β-propeller surfaces. This report once again reinforces the importance of structural data and the usefulness of homology models in functional classification. [ABSTRACT FROM AUTHOR]

The gene sod in Acidithiobacillus ferrooxidans may play a crucial role in its tolerance to the extremely acidic, toxic and oxidative environment of bioleaching. For insight into the anti-toxic mechanism of the bacteria, a three-dimensional (3D) molecular structure of the protein encoded by this gene was built by homology modeling techniques, refined by molecular dynamics simulations, assessed by PROFILE-3D and PROSTAT programs and its key residues were further detected by evolutionary trace analysis. Through these procedures, some trace residues were identified and spatially clustered. Among them, the residues of Asn38, Gly103 and Glu161 are randomly scattered throughout the mapped structure; interestingly, the other residues are all distinctly clustered in a subgroup near Fe atom. From these results, this gene can be confirmed at 3D level to encode the Fe-depending superoxide dismutase and subsequently play an anti-toxic role. Furthermore, the detected key residues around Fe binding site can be conjectured to be directly responsible for Fe binding and catalytic function. [ABSTRACT FROM AUTHOR]

Abstract: We discuss the question of which residues are sufficiently important for protein–protein interaction to be under notable evolutionary pressure. Its interest stems from the applicability of this knowledge in the reverse direction, to detect a protein–protein interface on a single protomer, starting from the rate of mutation of participating residues. Using the analysis of trajectories produced by the molecular dynamics simulations, we suggest that, in the case of water soluble proteins, a large fraction of evolutionarily privileged residues can be found by considering the dynamic behavior of the protein interface and by looking for residues which exchange water molecules with the bulk of the solvent outstandingly slowly (tentatively termed “dry residues”). We show that the dry interface residues are better conserved across homologues than the generic “geometric footprint” and can be quite reliably detected through comparative analysis of protein homologues, without strong dependence on the choice of method. Furthermore, we show that dry residues distinguish themselves through a set of biophysical properties consistent with the known mechanisms of protein oligomerization: their compositional shift toward nonpolar, overlap, and co-location with residues exhibiting low mobility, their two- to threefold increased propensity over the rest of the geometric footprint to form hydrogen bonds, and four- to almost tenfold increased likelihood to participate in formation of salt bridges. These properties, consistently, help understand the observed increase in the evolutionary pressure that dry residues experience. [Copyright &y& Elsevier]

Prestin, a member of the SLC26A family of anion transporters, is a polytopic membrane protein found in outer hair cells (OHCs) of the mammalian cochlea. Prestin is an essential component of the membrane-based motor that enhances electromotility of OHCs and contributes to frequency sensitivity and selectivity in mammalian hearing. Mammalian cells expressing prestin display a nonlinear capacitance (NLC), widely accepted as the electrical signature of electromotility. The associated charge movement requires intracellular anions reflecting the membership of prestin in the SLC26A family. We used the computational approach of evolutionary trace analysis to identify candidate functional (trace) residues in prestin for mutational studies. We created a panel of mutations at each trace residue and determined membrane expression and nonlinear capacitance associated with each mutant. We observe that several residue substitutions near the conserved sulfate transporter domain of prestin either greatly reduce or eliminate NLC, and the effect is dependent on the size of the substituted residue. These data suggest that packing of helices and interactions between residues surrounding the "sulfate transporter motif" is essential for normal prestin activity. [ABSTRACT FROM AUTHOR]

Understanding the molecular mechanisms of the specific interaction between transcription factor proteins and DNA is key to comprehend the regulation of gene expression and to develop technologies to engineer transcription factors. Thus far, although there have been several attempts to elucidate protein–DNA interaction through amino acid–base recognition codes, sequence based profiles, or physical models of interaction, the greatest successes in engineering DNA binding specificity remain experimental. Here we present the first systematic evidence of correlated evolutionary pressure at interacting amino acid residues and DNA base-pairs in transcription factors, and show that it can be used to rationally engineer DNA binding specificity. The correlation is between the relative evolutionary importance of protein residues and DNA bases, measured, respectively, in terms of the Evolutionary Trace (ET) rank and information entropy. The evolutionarily most important residues interact with the most conserved base-pairs within the response element while residues of least importance interact with the most variable base-pairs. The correlation averages 0.74 over 12 unrelated families of transcriptional regulators, including nuclear hormone receptors, basic helix–loop–helix, ETS- and homeo-domain family. To test the predictive power of this correlation, we targeted a mutational swap of top-ranked ET residues in a transcription factor, LRH-1. This redirects LRH-1 binding as predicted and showed that, in this case, evolutionary importance and binding specificity are coupled sufficiently strongly for the Evolutionary Trace to guide the computational design of DNA binding specificity. This establishes the existence of evolutionary importance correlation at protein–DNA interfaces, and demonstrates that it is a useful principle for the rational engineering of binding specificity. [Copyright &y& Elsevier]

The visual receptor of rods and cones is a covalent complex of the apoprotein, opsin, and the light-sensitive chromophore, 11- cis-retinal. This pigment must fulfill many functions including photoactivation, spectral tuning, signal transmission, inactivation, and chromophore regeneration. Rod and cone photoreceptors employ distinct families of opsins. Although it is well known that these opsin families provide unique ranges in spectral sensitivity, it is unclear whether the families have additional functional differences. In this study, we use evolutionary trace (ET) analysis of 188 vertebrate opsin sequences to identify functionally important sites in each opsin family. We demonstrate the following results. (1) The available vertebrate opsin sequences produce a definitive description of all five vertebrate opsin families. This is the first demonstration of sequence saturation prior to ET analysis, which we term saturated ET (SET). (2) The cone opsin classes have class-specific sites compared to the rod opsin class. These sites reside in the transmembrane region and tune the spectral sensitivity of each opsin class to its characteristic wavelength range. (3) The cytoplasmic loops, primarily responsible for signal transmission and inactivation, are essentially invariant in rod versus cone opsins. This indicates that the electrophysiological differences between rod and cone photoreceptors cannot be ascribed to differences in the protein interaction regions of the opsins. SET shows that chromophore binding and regeneration are the only aspects of opsin structure likely to have functionally significant differences between rods and cones, whereas excitatory and adaptational properties of the opsin families appear to be functionally invariant. [ABSTRACT FROM AUTHOR]

The studies of novel inhibitors of DNA topoisomerase I (Topo I) have already become very promising in cancer chemotherapy. Identifying the new drug-binding residues is playing an important role in the design and optimization of Topo I inhibitors. The designed compounds may have novel scaffolds, thus will be helpful to overcome the toxicities of current camptothecin (CPT) drugs and may provide a solution to cross resistance with these drugs. Multiple sequence alignments were performed on eukaryotic DNA topoisomerase I superfamily and thus the evolutionary tree was constructed. The Evolutionary Trace method was applied to identify functionally important residues of human Topo I. It has been demonstrated that class-specific hydrophobic residues Ala351, Met428, Pro431 are located around the 7,9-position of CPT, indicating suitable substitution of hydrophobic group on CPT will increase antitumor activity. The conservative residue Lys436 in the superfamily is of particular interest and new CPT derivatives designed based on this residue may greatly increase water solubility of such drugs. It has also been demonstrated that the residues Asn352 and Arg364 were conservative in the superfamily, whose mutation will render CPT resistance. As our molecular docking studies demonstrated they did not make any direct interaction with CPT, they are important drug-binding site residues for future design of novel non-camptothecin lead compounds. This work provided a strong basis for the design and synthesis of novel highly potent CPT derivatives and virtual screening for novel lead compounds. [ABSTRACT FROM AUTHOR]

Evolutionary trace (ET) and entropy are two related methods for analyzing a multiple sequence alignment to determine functionally important residues in proteins. In this article, these methods have been enhanced with a view to reinvestigate the issue of GPCR dimerization and oligomerization. In particular, cluster analysis has replaced the subjective visual analysis element of the original ET method. Previous applications of the ET method predicted two dimerization interfaces on the external transmembrane lipid-facing region of GPCRs; these were discussed in terms of dimerization and linear oligomers. Removing the subjective element of the ET method gives rise to the prediction of functionally important residues on the external face of each transmembrane helix for a large number of class A GPCRs. These results are consistent with a growing body of experimental information that, taken over many receptor subtypes, has implicated each transmembrane helix in dimeric interactions. In this application, entropy gave superior results to those obtained from the ET method in that its use gives rise to higher z-scores and fewer instances of z-scores below 3. [ABSTRACT FROM AUTHOR]

The α-D-phosphohexomutase superfamily is composed of four related enzymes that catalyze a reversible, intramolecular phosphoryl transfer on their sugar substrates. The enzymes in this superfamily play important and diverse roles in carbohydrate metabolism in organisms from bacteria to humans. Recent structural and mechanistic studies of one member of this superfamily, phosphomannomutase/phosphoglucomutase (PMM/PGM) from Pseudomonas aeruginosa, have provided new insights into enzyme mechanism and substrate recognition. Here we use sequence-sequence and sequence-structure comparisons via evolutionary trace analysis to examine 71 members of the α-D-phosphohexomutase superfamily. These analyses show that key residues in the active site, including many of those involved in substrate contacts in the P. aeruginosa PMM/PGM complexes, are conserved throughout the enzyme family. Several important regions show class-specific differences in sequence that appear to be correlated with differences in substrate specificity exhibited by subgroups of the family. In addition, we describe the translocation of a 20-residue segment containing the catalytic phosphoserine of phosphoacetylglucosamine mutase, which uniquely identifies members of this subgroup. [ABSTRACT FROM AUTHOR]

In order to identify the amino acids that determine protein structure and function it is useful to rank them by their relative importance. Previous approaches belong to two groups; those that rely on statistical inference, and those that focus on phylogenetic analysis. Here, we introduce a class of hybrid methods that combine evolutionary and entropic information from multiple sequence alignments. A detailed analysis in insulin receptor kinase domain and tests on proteins that are well-characterized experimentally show the hybrids'' greater robustness with respect to the input choice of sequences, as well as improved sensitivity and specificity of prediction. This is a further step toward proteome scale analysis of protein structure and function. [Copyright &y& Elsevier]

Using the tumor suppressor VHL protein as an example, we show that detailed analysis of conservation versus variation pattern in the multiple alignment can be coupled with the genomic pathway/complex conservation analysis to provide a more complete picture of the entire interaction/regulatory network. Results from the present study have allowed us to hypothesize that two additional proteins are involved in the VHL-mediated regulation of angiogenesis. Detailed modeling also has led to a prediction of the possible interaction mode between the known and the proposed parts of the VHL complex. To aid in an analysis of the VHL protein regulation of HIF-1α degradation, an important and only partially understood process that directly influences angiogenesis, we performed a comprehensive search for the orthologs of the VHL as well as for VHL-interacting proteins in all the available eukaryotic genomes. Analysis of a multiple alignment of thus identified VHL orthologs reveals an unusually high degree of conservation of the surface amino acid residues that almost exactly correspond to positions mutated in the VHL disease-associated tumors. In addition, these positions form well-defined clusters in three-dimensional space, and presence or absence of individual clusters correlates with the presence or absence of pathway elements in different genomes. We have also shown that relation trees derived from the multiple sequence alignment, functional surface-mapping, and HIF-1α degradation pathway structure are in complete agreement, linking the functional and structural evolution of the VHL protein and VHL-dependent HIF-1α degradation complex. [ABSTRACT FROM AUTHOR]

Starting from the hypothesis that evolutionarily important residues form a spatially limited cluster in a protein''s native fold, we discuss the possibility of detecting a non-native structure based on the absence of such clustering. The relevant residues are determined using the Evolutionary Trace method. We propose a quantity to measure clustering of the selected residues on the structure and show that the exact values for its average and variance over several ensembles of interest can be found. This enables us to study the behavior of the associated z-scores. Since our approach rests on an analytic result, it proves to be general, customizable, and computationally fast. We find that clustering is indeed detectable in a large representative protein set. Furthermore, we show that non-native structures tend to achieve lower residue-clustering z-scores than those attained by the native folds. The most important conclusion that we draw from this work is that consistency between structural and evolutionary information, manifested in clustering of key residues, imposes powerful constraints on the conformational space of a protein. [Copyright &y& Elsevier]

A common difficulty in post genomics biology is that large-scale techniques of data collection often strip away information on the biological context of these data. The result is a massive number of disconnected observations on sequence, structure, and function from which underlying patterns and biological meaning are obscured. One solution is to build computational filters that pick out sufficiently few facts, relevant to a query, that their relationship is immediately apparent and experimentally testable. Typically, these filters rely on mathematics and statistics, and on first principles from physics and chemistry. We show here that evolution itself can be used to filter sequence and structure data in order to identify evolutionarily important amino acids. A general property of these residues is that they form clusters in native protein structures and point to regions where mutations have the greatest biological impact. The result is an accurate method of functional site annotation that is scalable for structural proteomics. [ABSTRACT FROM AUTHOR]

Group I chaperonins are a highly conserved family of essential proteins that self-assemble into molecular nanoboxes that mediate the folding of cytoplasmic proteins in bacteria and organelles. GroEL, the chaperonin of Escherichia coli , is the archetype of the family. Protein folding-independent functions have been described for numerous chaperonins, including HtpB, the chaperonin of the bacterial pathogen Legionella pneumophila . Several protein folding-independent functions attributed to HtpB are not shared by GroEL, suggesting that differences in the amino acid (aa) sequence between these two proteins could correlate with functional differences. GroEL and HtpB differ in 137 scattered aa positions. Using the Evolutionary Trace (ET) bioinformatics method, site-directed mutagenesis, and a functional reporter test based upon a yeast-two-hybrid interaction with the eukaryotic protein ECM29, it was determined that out of those 137 aa, ten (M68, M212, S236, K298, N507 and the cluster AEHKD in positions 471-475) were involved in the interaction of HtpB with ECM29. GroEL was completely unable to interact with ECM29, but when GroEL was modified at those 10 aa positions, to display the HtpB aa, it acquired a weak ability to interact with ECM29. This constitutes proof of concept that the unique functional abilities of HtpB can be mapped to specific aa positions.

Our understanding of enzymes with high substrate ambiguity remains limited because their large active sites allow substrate docking freedom to an extent that seems incompatible with stereospecificity. One possibility is that some of these enzymes evolved a set of evolutionarily fitted sequence positions that stringently allow switching substrate ambiguity and chiral specificity. To explore this hypothesis, we targeted for mutation a serine ester hydrolase (EH 3 ) that exhibits an impressive 71-substrate repertoire but is not stereospecific ( e.e. 50%). We used structural actions and the computational evolutionary trace method to explore specificity-swapping sequence positions and hypothesized that position I244 was critical. Driven by evolutionary action analysis, this position was substituted to leucine, which together with isoleucine appears to be the amino acid most commonly present in the closest homologous sequences (max. identity, ca . 67.1%), and to phenylalanine, which appears in distant homologues. While the I244L mutation did not have any functional consequences, the I244F mutation allowed the esterase to maintain a remarkable 53-substrate range while gaining stereospecificity properties ( e.e. 99.99%). These data support the possibility that some enzymes evolve sequence positions that control the substrate scope and stereospecificity. Such residues, which can be evolutionarily screened, may serve as starting points for further designing substrate-ambiguous, yet chiral-specific, enzymes that are greatly appreciated in biotechnology and synthetic chemistry., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2021 The Author(s).)

The contribution of newly designed salt bridges to protein stabilization remains controversial even today. In order to solve this problem, we investigated salt bridges from two aspects: spatial distribution and evolutionary characteristics of salt bridges. Firstly, we analyzed spatial distribution of salt bridges in proteins, elucidating the basic requirements of forming salt bridges. Then, from an evolutionary point of view, the evolutionary characteristics of salt bridges as well as their neighboring residues were investigated in our study. The results demonstrate that charged residues appear more frequently than other neutral residues at certain positions of sequence even under evolutionary pressure, which are able to form electrostatic interactions that could increase the evolutionary stability of corresponding amino acid regions, enhancing their importance to stability of proteins. As a corollary, we conjectured that the newly designed salt bridges with more contribution to proteins, not only, are qualified spatial distribution of salt bridges, but also, are needed to further increase the evolutionary stability of corresponding amino acid regions. Based on analysis, the 8 mutations were accordingly constructed in the 1,4-α-glucan branching enzyme (EC 2.4.1.18, GBE) from Geobacillus thermoglucosidans STB02, of which 7 mutations improved thermostability of GBE. The enhanced thermostability of 7 mutations might be a result of additional salt bridges on residue positions that at least one of amino acids positions is conservative, improving their contribution of stabilization to proteins.