Evolutionary Trace Analysis Research Articles

The structural and functional determinants of the Axin and Dishevelled DIX domains

BackgroundThe dishevelled and axin genes encode multi-domain proteins that play key roles in WNT signalling. Dishevelled prevents β-catenin degradation by interfering with the interaction of β-catenin with the degradation-mediating Axin-APC-GSK3β complex. This interference leads to an accumulation of cytoplasmic β-catenin, which enters the nucleus and interacts with transcription factors that induce expression of Wnt-target genes. Axin, as a component of the degradation-mediating complex, is a potent negative regulator of Wnt signalling, whereas Dishevelled is a potent activator. Both Dishevelled and Axin possess a DIX (Dishevelled/Axin) domain, which mediates protein-protein interactions, specifically homodimerization.ResultsAn evolutionary trace analysis of DIX domains identified conserved residues which, when mapped onto the crystal structure of the Axin DIX domain and a comparative model of the Dishevelled DIX domain, allow their categorization as residues of either structural or functional importance. We identify residues that are structural and functional determinants of the DIX domain fold, as well as those that are specific to homodimerization of Axin and Dishevelled.ConclusionThis report provides the first explanation of the mutant phenotypes caused by non-synonymous substitutions in the Dishevelled and Axin DIX domain by correlating their presumed functional significance with molecular structure.

Gene expression, antiparasitic activity, and functional evolution of the drosomycin family

Drosophila employs various antimicrobial peptides as effective weapons to defend against diverse pathogens. Drosomycin is an inducible antifungal peptide initially isolated from the Drosophila melanogaster haemolymph. Here we report the expression pattern of seven drosomycin genes in four different developmental stages (egg, larva, pupa and adult). Results show that drosomycin and drosomycin-2 are expressed in larva, pupa and adult, whereas drosomycin-1 and drosomycin-6 were not detected in all the stages. Moreover, all the seven drosomycin genes are shut off in egg. Functional comparison of recombinant drosomycin and drosomycin-2, both with identical expression pattern, produced from Escherichia coli, revealed their significant differences in potency against a specific fungal species. In addition, we found for the first time that drosomycin and drosomycin-2 both are antiparasitic peptides which show inhibitory effect on the ookinete development of the parasite Plasmodium berghei with differential potenc y. Functional differentiation between them was further evaluated by evolutionary trace analysis which identified two evolutionary epitopes (named α- and γ-patch, respectively) and an important site in the m-loop. Substitutions in these regions are possibly associated with the antifungal and antiparasitic potency difference among members of the drosomycin family.

Evolutionary trace analysis at the ligand binding site of laccase

Laccase belongs to the family of blue multi-copper oxidases and are capable of oxidizing a wide range of aromatic compounds. Laccases have industrial applications in paper pulping or bleaching and hydrocarbon bioremediation as a biocatalyst. We describe the design of a laccase with broader substrate spectrum in bioremediation. The application of evolutionary trace (ET) analysis of laccase at the ligand binding site for optimal design of the enzyme is described. In this attempt, class specific sites from ET analysis were mapped onto known crystal structure of laccase. The analysis revealed 162PHE as a critical residue in structure function relationship studies.

Docking of protein molecular surfaces with evolutionary trace analysis

We have developed a new method to predict protein- protein complexes based on the shape complementarity of the molecular surfaces, along with sequence conservation obtained by evolutionary trace (ET) analysis. The docking is achieved by optimization of an object function that evaluates the degree of shape complementarity weighted by the conservation of the interacting residues. The optimization is carried out using a genetic algorithm in combination with Monte Carlo sampling. We applied this method to CAPRI targets and evaluated the performance systematically. Consequently, our method could achieve native-like predictions in several cases. In addition, we have analyzed the feasibility of the ET method for docking simulations, and found that the conservation information was useful only in a limited category of proteins (signal related proteins and enzymes).

Solution Structure of Inhibitor-Free Human Metalloelastase (MMP-12) Indicates an Internal Conformational Adjustment

Macrophage metalloelastase or matrix metalloproteinase-12 (MMP-12) appears to exacerbate atherosclerosis, emphysema, aortic aneurysm, rheumatoid arthritis, and inflammatory bowel disease. An inactivating E219A mutation, validated by crystallography and NMR spectra, prevents autolysis of MMP-12 and allows us to determine its NMR structure without an inhibitor. The structural ensemble of the catalytic domain without an inhibitor is based on 2813 nuclear Overhauser effects (NOEs) and has an average RMSD to the mean structure of 0.25 Å for the backbone and 0.61 Å for all heavy atoms for residues Trp109-Gly263. Compared to crystal structures of MMP-12, helix B (hB) at the active site is unexpectedly more deeply recessed under the β-sheet. This opens a pocket between hB and β-strand IV in the active-site cleft. Both hB and an internal cavity are shifted toward β-strand I, β-strand III, and helix A on the back side of the protease. About 25 internal NOE contacts distinguish the inhibitor-free solution structure and indicate hB's greater depth and proximity to the sheet and helix A. Line broadening and multiplicity of amide proton NMR peaks from hB are consistent with hB undergoing a slow conformational exchange among subtly different environments. Inhibitor-binding-induced perturbations of the NMR spectra of MMP-1 and MMP-3 map to similar locations across MMP-12 and encompass the internal conformational adjustments. Evolutionary trace analysis suggests a functionally important network of residues that encompasses most of the locations adjusting in conformation, including 18 residues with NOE contacts unique to inhibitor-free MMP-12. The conformational change, sequence analysis, and inhibitor perturbations of NMR spectra agree on the network they identify between structural scaffold and the active site of MMPs.

Homology modeling and evolutionary trace analysis of superoxide dismutase from extremophile Acidithiobacillus ferrooxidans

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.

Cell surface proteins in archaeal and bacterial genomes comprising “LVIVD”, “RIVW” and “LGxL” tandem sequence repeats are predicted to fold as β-propeller

Proteins that share even low sequence homologies are known to adopt similar folds. The β-propeller structural motif is one such example. Identifying sequences that adopt a β-propeller fold is useful to annotate protein structure and function. Often, tandem sequence repeats provide the necessary signal for identifying β-propellers in proteins. In our recent analysis to identify cell surface proteins in archaeal and bacterial genomes, we identified some proteins that contain novel tandem repeats “LVIVD”, “RIVW” and “LGxL”. In this work, based on protein fold predictions and three-dimensional comparative modeling methods, we predicted that these repeat types fold as β-propeller. Further, the evolutionary trace analysis of all proteins constituting amino acid sequence repeats in β-propellers suggest that the novel repeats have diverged from a common ancestor.

Phosphoprotein and Phosphopeptide Interactions with the FHA Domain from Arabidopsis Kinase-Associated Protein Phosphatase

FHA domains are phosphoThr recognition modules found in diverse signaling proteins, including kinase-associated protein phosphatase (KAPP) from Arabidopsis thaliana. The kinase-interacting FHA domain (KI-FHA) of KAPP targets it to function as a negative regulator of some receptor-like kinase (RLK) signaling pathways important in plant development and environmental responses. To aid in the identification of potential binding sites for the KI-FHA domain, we predicted (i) the structure of a representative KAPP-binding RLK, CLAVATA1, and (ii) the functional surfaces of RLK kinase domains using evolutionary trace analysis. We selected phosphopeptides from KAPP-binding Arabidopsis RLKs for in vitro studies of association with KI-FHA from KAPP. Three phosphoThr peptide fragments from the kinase domain of CLV1 or BAK1 were found to bind KI-FHA with KD values of 8-20 microM, by NMR or titration calorimetry. Their affinity is driven by favorable enthalpy and solvation entropy gain. Mutagenesis of these three threonine sites suggests Thr546 in the C-lobe of the BAK1 kinase domain to be a principal but not sole site of KI-FHA binding in vitro. The brassinosteroid receptor BRI1 and KAPP are shown to associate in vivo and in vitro. Further genetic studies indicate that KAPP may be a negative regulator of the BRI1 signaling transduction pathway. 15N-Labeled KI-FHA was titrated with the GST-BRI1 kinase domain and monitored by NMR. BRI1 interacts with the same 3/4, 4/5, 6/7, 8/9, and 10/11 recognition loops of KI-FHA, with similar affinity as the phosphoThr peptides.

Evolutionary trace analysis of plant haemoglobins: implications for site-directed mutagenesis

Haemoglobins are found ubiquitously in eukaryotes and many bacteria. In plants, haemoglobins were first identified in species, which can fix nitrogen via symbiosis with bacteria. Recent findings suggest that another class of haemoglobins termed as nonsymbiotic haemoglobins are present through out the plant kingdom and are expressed differentially during plant development. Limited data available suggests that non-symbiotic haemoglobins are involved in hypoxic stress and oversupply of nutrients. Due to lack of information on structurally conserved, functionally important residues in non-symbiotic haemoglobins, further studies to elucidate the molecular mechanisms underlying the biological role are hampered. To determine functionally important residues in non-symbiotic haemoglobins, I have analyzed a number of sequences from plant haemoglobin family, in the context of the known crystal structures of plant by evolutionary trace method. Results indicate that the, evolutionary trace method like conventional phylogentic analysis, could resolve phylogentic relationships between plant haemoglobin family. Evolutionary trace analysis has identified candidate functional (trace) residues that uniquely characterize the heme-binding pocket, dimer interface and possible novel functional surfaces. Such residues from specific three-dimensional clusters might be of functional importance in nonsymbiotic haemoglobins. These data, together with our improved knowledge of possible functional residues, can be used in future structure-function analysis experiments.

A chimeric cysteine protease of Plasmodium berghei engineered to resemble the Plasmodium falciparum protease falcipain-2

The cysteine proteases falcipain-2 and falcipain-3 are hemoglobinases and potential targets for chemotherapy directed against Plasmodium falciparum, the most important human malaria parasite. Most in vivo evaluations of candidate antimalarials are conducted in murine malaria models, and falcipain homologs from rodent malaria parasites differ importantly from falcipain-2 and falcipain-3. We expressed berghepain-2, the single homolog of falcipain-2 and falcipain-3 of the rodent parasite P. berghei, in Escherichia coli, and characterized the refolded active enzyme. Berghepain-2 was biochemically very similar to the previously characterized rodent plasmodial protease vinckepain-2, but differed from falcipain-2 and falcipain-3 in its fine substrate and inhibitor specificity. We then used homology modeling and evolutionary trace analysis to predict key amino acids that mediate functional differences between falcipain-2 and berghepain-2. Thirteen amino acids were sequentially altered to replace berghepain-2 residues with those in falcipain-2. Mutant enzymes varied in activity and sensitivity to inhibitors. A berghepain-2 mutant with eight substitutions retained good activity and demonstrated fine substrate and inhibitor sensitivity more similar to that of falcipain-2 than berghepain-2. These results suggest that, to facilitate drug discovery, we can produce mutant animal model malaria parasites with biochemical properties more like those of the key drug target, P. falciparum.

Ancestral sequence evolutionary trace and crystal structure analyses of alkaline α‐amylase from Bacillus sp. KSM‐1378 to clarify the alkaline adaptation process of proteins

The crystal structure of alkaline liquefying alpha-amylase (AmyK) from the alkaliphilic Bacillus sp. KSM-1378 was determined at 2.1 A resolution. The AmyK structure belongs to the GH13 glycoside hydrolase family, which consists of three domains, and bound three calcium and one sodium ions. The alkaline adaptation mechanism of AmyK was investigated by the ancestral sequence evolutionary trace method and by extensive comparisons between alkaline and nonalkaline enzyme structures, including three other protein families: protease, cellulase, and phosphoserine aminotransferase. The consensus change for the alkaline adaptation process was a decrease in the Lys content. The loss of a Lys residue is associated with ion pair remodeling, which mainly consists of the loss of Lys-Asp/Glu ion pairs and the acquisition of Arg ion pairs, preferably Arg-Glu. The predicted replacements of the positively charged amino acids were often, although not always, used for ion pair remodeling.

Essential Helix Interactions in the Anion Transporter Domain of Prestin Revealed by Evolutionary Trace Analysis

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.

Loss of neurturin in frog—Comparative genomics study of GDNF family ligand-receptor pairs

Four different GDNF family ligand (GFL)-receptor (GFRα) binding pairs exist in mammals, and they all signal via the RET receptor tyrosine kinase. However, the evolution of these molecules is poorly understood. We identified orthologs of all four GFRα receptors and GRAL (GDNF Receptor Alpha-Like) in all vertebrate classes, and a predicted GFR-like protein in several invertebrates. In addition, Gas1 (growth arrest-specific 1), a distant member of the GFR-superfamily, is present in both vertebrates and invertebrates. Analysis of exon structures suggests a common origin of GFR-superfamily proteins and early divergence of Gas1 from the common ancestor. Bony fishes have orthologs of all four mammalian GFLs, consistent with genome duplications in early vertebrates. Surprisingly, the clawed frog and chicken have only three GFLs: synteny analysis indicates loss of neurturin in frog and of persephin in chicken. Evolutionary trace analysis and protein structure homology modeling points at GDNF as the endogenous ligand of frog GFRα2.

Coupling of Agonist Binding to Effector Domain Activation in Metabotropic Glutamate-like Receptors

Many membrane receptors are made of a ligand binding domain and an effector domain mediating intracellular signaling. This is the case for the metabotropic glutamate-like G-protein-coupled receptors. How ligand binding leads to the active conformation of the effector domain in such receptors is largely unknown. Here, we used an evolutionary trace analysis and mutagenesis to identify critical residues involved in the allosteric coupling between the Venus flytrap ligand binding domain (VFT) and the heptahelical G-protein activating domain of the metabotropic glutamate-like receptors. We have shown that a conserved interdomain disulfide bridge is required for this allosteric interaction. Taking into account that these receptors are homodimers, this finding provides important new information explaining how the different conformations of the dimer of VFT lead to different signaling of such dimeric receptors.

Evolutionary identification of a subtype specific functional site in the ligand binding domain of steroid receptors

Nuclear receptors are ubiquitous eukaryotic ligand-activated transcription factors that modulate gene expression through varied interactions. However, the highly conserved functional sites known today seem insufficient to explain receptor specific recruitment of different coactivator and corepressor proteins and regulation of transcription. To search for new receptor-subtype specific functional sites, we applied difference evolutionary trace (difference ET) analysis to the ligand binding domain of steroid receptors, a subgroup of the nuclear receptor (NR) family. This computational approach identified a new functional site located on a surface opposite to currently known protein-protein interaction sites and distinct from the ligand binding pocket. Strikingly, the literature shows that in vivo variations at residues in the new site are linked to androgen resistance and leukemia, and our own targeted mutations to this site lower but do not eradicate transcriptional activation by estrogen receptor alpha (ERalpha), with reduced ligand binding affinity and SRC-1 interaction. Thus, these data demonstrate that this evolutionary important surface can function as an allosteric site that modulates some but not all receptor binding interactions. Evolutionary analysis further shows that this allosteric regulatory site is shared among all NRs from groups 2 (HNF4-like) and 4 (NGFIB-like), suggesting a role among many nuclear receptors. Its concave structure, hydrophobic composition, and residue variability among nuclear receptors further suggest that it would be amenable for specific drug design. This highlights the power of evolutionary information for the identification of new functional sites even in a protein family as well studied as NRs.

Hsc70 Contacts Helix III of the J Domain from Polyomavirus T Antigens: Addressing a Dilemma in the Chaperone Hypothesis of How They Release E2F from pRb

Hsc70's expected binding site on helix II of the J domain of T antigens appears to be blocked in its structure bound to tumor suppressor pRb. We used NMR to map where mammalian Hsc70 binds the J domain of murine polyomavirus T antigens (PyJ). The ATPase domain of Hsc70 unexpectedly has its biggest effects on the NMR peak positions of the C-terminal end of helix III of PyJ. The Hsc70 ATPase domain protects the C-terminal end of helix III of PyJ from an uncharged paramagnetic probe of chelated Gd(III), clearly suggesting the interface. Effects on the conserved HPD loop and helix II of PyJ are smaller. The NMR results are supported by a novel assay of Hsc70's ATP hydrolysis showing that mutations of surface residues in PyJ helix III impair PyJ-dependent stimulation of Hsc70 activity. Evolutionary trace analysis of J domains suggests that helix III usually may join helix II in contributing specificities for cognate hsp70s. Our novel evidence implicating helix III differs from evidence that Escherichia coli DnaK primarily affects helix II and the HPD loop of DnaJ. We find the pRb-binding fragment of E2F1 to be intrinsically unfolded and a good substrate for Hsc70 in vitro. This suggests that E2F1 could be a substrate for Hsc70 recruited by T antigen to an Rb family member. Importantly, our results strengthen the chaperone hypothesis for E2F release from an Rb family member by Hsc70 recruited by large T antigen. That is, it now appears that Hsc70 can freely access helix III and the HPD motif of large T antigen bound to an Rb family member.

Binding site prediction of galanin peptide using evolutionary trace method

Galanin is a neuropeptide with aminoacid length ranging from 29 to 31 is widely distributed in central and peripheral nervous system. Galanin controls various psychological processes such as sensation of pain, learning, feeding, and sexual behaviour. The N-terminal region of this neuropeptide has highly conserved 15 amino acids, which is triggered by galanin receptors. We performed evolutionary trace analysis for galanin sequences to gather information about functional residues. The consensus pattern given by the evolutionary trace (ET) analysis is supported by CLUSTALW and WEBLOGO results. Our observations strongly suggest the presence of functional residues in the N-terminal region of galanin for agonist-receptor binding.

Solution structure of the antifreeze-like domain of human sialic acid synthase.

The structure of the C-terminal antifreeze-like (AFL) domain of human sialic acid synthase was determined by NMR spectroscopy. The structure comprises one alpha- and two single-turn 3(10)-helices and two beta-strands, and is similar to those of the type III antifreeze proteins. Evolutionary trace analyses of the type III antifreeze protein family suggested that the class-specific residues in the human and bacterial AFL domains are important for their substrate binding, while the class-specific residues of the fish antifreeze proteins are gathered on the ice-binding surface.

Hsc70 Contacts Helix III of DnaJ‐like Domain from Polyomavirus T Antigens

The J domain of large T antigen, Hsc70, and its hydrolysis of ATP were found necessary for polyomavirus infection to release E2Fs from pRb family members. Hsc70's expected binding site on helix II of the DnaJ-like domain of T antigens is blocked in the crystal structure of a T antigen fragment bound to tumor suppressor pRb. We used NMR to find recognition sites for mammalian Hsc70 upon the J domain of murine polyomavirus T antigens (PyJ). The Hsc70 ATPase domain unexpectedly shifts most of the amide NMR peaks of the C-terminal end of helix III of PyJ. The Hsc70 ATPase domain buries the C-terminal end of helix III of PyJ, protecting it from an uncharged paramagnetic probe. This region corresponds to the QKRAA motif of orthodox J domains. Evolutionary trace analysis suggests helix III of J domains to be of wide importance to specificity for hsp70s. The NMR-detected Hsc70 contacts map to a PyJ surface homologous to a solvent-accessible surface of the SV40 T antigen fragment bound to pRb. This surface of helix III and the HPD loop may allow Hsc70 to form a ternary complex with T antigen and pRb, avoiding obstructions to helix II. The helix III docking site is consistent with large T antigen being able to free an E2F from pRb by positioning the substrate-binding domain of Hsc70 for direct contact with the E2F. NMR spectra show the pRb-binding region of E2F1 (aa243–437) to be unfolded in the absence of DP1. Hsc70 is known to bind unfolded proteins. In the absence of DP1, E2F1(aa243–437) acts as a substrate for Hsc70 in vitro, stimulating its ATPase activity in a manner dependent on PyJ from T antigen. If segments of an E2F are unstructured in vivo, it is conceivable that Hsc70, recruited by large T antigen to a complex with pRb, might act on the unstructured portion of the E2F when liberating it from the pRb family member.

Structural Characterization of the UL25 DNA-Packaging Protein from Herpes Simplex Virus Type 1

Herpesviruses replicate their double stranded DNA genomes as high-molecular-weight concatemers which are subsequently cleaved into unit-length genomes by a complex mechanism that is tightly coupled to DNA insertion into a preformed capsid structure, the procapsid. The herpes simplex virus type 1 UL25 protein is incorporated into the capsid during DNA packaging, and previous studies of a null mutant have demonstrated that its function is essential at the late stages of the head-filling process, either to allow packaging to proceed to completion or for retention of the viral genome within the capsid. We have expressed and purified an N-terminally truncated form of the 580-residue UL25 protein and have determined the crystallographic structure of the region corresponding to amino acids 134 to 580 at 2.1-Angstroms resolution. This structure, the first for any herpesvirus protein involved in processing and packaging of viral DNA, reveals a novel fold, a distinctive electrostatic distribution, and a unique "flexible" architecture in which numerous flexible loops emanate from a stable core. Evolutionary trace analysis of UL25 and its homologues in other herpesviruses was used to locate potentially important amino acids on the surface of the protein, leading to the identification of four putative docking regions for protein partners.