Functional analysis of the large periplasmic loop of the Escherichia coli K‐12 WaaL O‐antigen ligase

WaaL is a membrane enzyme implicated in ligating undecaprenyl‐diphosphate (Und‐PP)‐linked O antigen to lipid A‐core oligosaccharide. We determined the periplasmic location of a large (EL5) and small (EL4) adjacent loops in the Escherichia coli K‐12 WaaL. Structural models of the EL5 from the K‐12, R1 and R4 E. coli ligases were generated by molecular dynamics. Despite the poor amino acid sequence conservation among these proteins, the models afforded similar folds consisting of two pairs of almost perpendicular α‐helices. One α‐helix in each pair contributes a histidine and an arginine facing each other, which are highly conserved in WaaL homologues. Mutations in either residue rendered WaaL non‐functional, since mutant proteins were unable to restore O antigen surface expression. Replacements of residues located away from the putative catalytic centre and non‐conserved residues within the centre itself did not affect ligation. Furthermore, replacing a highly conserved arginine in EL4 with various amino acids inactivates WaaL function, but functionality reappears when the positive charge is restored by a replacement with lysine. These results lead us to propose that the conserved amino acids in the two adjacent periplasmic loops could interact with Und‐PP, which is the common component in all WaaL substrates.     

Discarding functional residues from the substitution table improves predictions of active sites within three-dimensional structures

Substitutions of individual amino acids in proteins may be under very different evolutionary restraints depending on their structural and functional roles. The Environment Specific Substitution Table (ESST) describes the pattern of substitutions in terms of amino acid location within elements of secondary structure, solvent accessibility, and the existence of hydrogen bonds between side chains and neighbouring amino acid residues. Clearly amino acids that have very different local environments in their functional state compared to those in the protein analysed will give rise to inconsistencies in the calculation of amino acid substitution tables. Here, we describe how the calculation of ESSTs can be improved by discarding the functional residues from the calculation of substitution tables. Four categories of functions are examined in this study: protein–protein interactions, protein–nucleic acid interactions, protein–ligand interactions, and catalytic activity of enzymes. Their contributions to residue conservation are measured and investigated. We test our new ESSTs using the program CRESCENDO, designed to predict functional residues by exploiting knowledge of amino acid substitutions, and compare the benchmark results with proteins whose functions have been defined experimentally. The new methodology increases the Z-score by 98% at the active site residues and finds 16% more active sites compared with the old ESST. We also find that discarding amino acids responsible for protein–protein interactions helps in the prediction of those residues although they are not as conserved as the residues of active sites. Our methodology can make the substitution tables better reflect and describe the substitution patterns of amino acids that are under structural restraints only.     

The conservation profile of a protein bears the imprint of the molecule that is evolutionarily coupled to the protein

The conservation profile of a protein is a curve of the conservation levels of amino acids along the sequence. Biologists are usually more interested in individual points on the curve (namely, the conserved amino acids) than the overall shape of the curve. Here, we show that the conservation curves of proteins bear the imprints of molecules that are evolutionarily coupled to the proteins. Our method is based on recent studies that a sequence conservation profile is quantitatively linked to its structural packing profile. We find that the conservation profiles of nucleic acid (NA) binding proteins are better correlated with the packing profiles of the protein–NA complexes than those of the proteins alone. This indicates that a nucleic acid binding protein evolves to accommodate the nucleic acid in such a way that the residues involved in binding have their conservation levels closely coupled with the specific nucleotides. Proteins 2015; 83:1407–1413. © 2015 Wiley Periodicals, Inc.     

Distinguishing structural and functional restraints in evolution in order to identify interaction sites

Structural genomics projects are producing many three-dimensional structures of proteins that have been identified only from their gene sequences. It is therefore important to develop computational methods that will predict sites involved in productive intermolecular interactions that might give clues about functions. Techniques based on evolutionary conservation of amino acids have the advantage over physiochemical methods in that they are more general. However, the majority of techniques neither use all available structural and sequence information, nor are able to distinguish between evolutionary restraints that arise from the need to maintain structure and those that arise from function. Three methods to identify evolutionary restraints on protein sequence and structure are described here. The first identifies those residues that have a higher degree of conservation than expected: this is achieved by comparing for each amino acid position the sequence conservation observed in the homologous family of proteins with the degree of conservation predicted on the basis of amino acid type and local environment. The second uses information theory to identify those positions where environment-specific substitution tables make poor predictions of the overall amino acid substitution pattern. The third method identifies those residues that have highly conserved positions when three-dimensional structures of proteins in a homologous family are superposed. The scores derived from these methods are mapped onto the protein three-dimensional structures and contoured, allowing identification clusters of residues with strong evolutionary restraints that are sites of interaction in proteins involved in a variety of functions. Our method differs from other published techniques by making use of structural information to identify restraints that arise from the structure of the protein and differentiating these restraints from others that derive from intermolecular interactions that mediate functions in the whole organism.     

Solution structure of a conserved domain of antizyme: a protein regulator of polyamines

Antizyme and its isoforms are members of an unusual yet broadly conserved family of proteins, with roles in regulating polyamine levels within cells. Antizyme has the ability to bind and inhibit the enzyme ornithine decarboxylase (ODC), targeting it for degradation at the proteasome; antizyme is also known to affect the transport of polyamines and interact with the antizyme inhibitor protein (AZI), as well as the cell-cycle protein cyclin D1. In the present work, NMR methods were used to determine the solution structure of a stable, folded domain of mammalian antizyme isoform-1 (AZ-1), consisting of amino acid residues 87-227. The protein was found to contain eight beta strands and two alpha helices, with the strands forming a mixed parallel and antiparallel beta sheet. At the level of primary sequence, antizyme is not similar to any protein of known structure, and results show that antizyme exhibits a novel arrangement of its strands and helices. Interestingly, however, the fold of antizyme is similar to that found in a family of acetyl transferases, as well as translation initiation factor IF3, despite a lack of functional relatedness between these proteins. Structural results, combined with amino acid sequence comparisons, were used to identify conserved features among the various homologues of antizyme and their isoforms. Conserved surface residues, including a cluster of acidic amino acids, were found to be located on a single face of antizyme, suggesting this surface is a possible site of interaction with target proteins such as ODC. This structural model provides an essential framework for an improved future understanding of how the different parts of antizyme play their roles in polyamine regulation.     

Identification of residues in a hydrophilic loop of the Papaver rhoeas S protein that play a crucial role in recognition of incompatible pollen.

The self-incompatibility response involves S allele-specific recognition between stigmatic S proteins and incompatible pollen. This response results in pollen inhibition. Defining the amino acid residues within the stigmatic S proteins that participate in S allele-specific inhibition of incompatible pollen is essential for the elucidation of the molecular basis of the self-incompatibility response. We have constructed mutant derivatives of the S1 protein from Papaver rhoeas by using site-directed mutagenesis and have tested their biological activity. This has enabled us to identify amino acid residues in the stigmatic S proteins of P. rhoeas that are required for S-specific inhibition of incompatible pollen. We report here the identification of several amino acid residues in the predicted hydrophilic loop 6 of the P. rhoeas stigmatic S1 protein that are involved in the inhibition of S1 pollen. Mutation of the only hypervariable amino acid, which is situated in this loop, resulted in the complete loss of ability of the S protein to inhibit S1 pollen. This clearly demonstrates that this residue plays a crucial role in pollen recognition and may also participate in defining allelic specificity. We have also established the importance of highly conserved amino acids adjacent to this hypervariable site. Our studies demonstrate that both variable and conserved amino acids in the region of the S protein corresponding to surface loop 6 are key elements that play a role in the recognition and inhibition of incompatible pollen in the pollen-pistil self-incompatibility reaction.     

In Silico Analysis of Prion Protein Mutants: A Comparative Study by Molecular Dynamics Approach

Polymorphisms in the human prion proteins lead to amino acid substitutions by the conversion of PrPC to PrPSc and amyloid formation, resulting in prion diseases such as familial Creutzfeldt–Jakob disease, Gerstmann–Straussler–Scheinker disease and fatal familial insomnia. Cation–π interaction is a non-covalent binding force that plays a significant role in protein stability. Here, we employ a novel approach by combining various in silico tools along with molecular dynamics simulation to provide structural and functional insight into the effect of mutation on the stability and activity of mutant prion proteins. We have investigated impressions of prevalent mutations including 1E1S, 1E1P, 1E1U, 1E1P, 1FKC and 2K1D on the human prion proteins and compared them with wild type. Structural analyses of the models were performed with the aid of molecular dynamics simulation methods. According to our results, frequently occurred mutations were observed in conserved sequences of human prion proteins and the most fluctuation values appear in the 2K1D mutant model at around helix 4 with residues ranging from 190 to 194. Our observations in this study could help to further understand the structural stability of prion proteins.     

DNA binding analysis of glucocorticoid receptor specificity mutants.

The glucocorticoid receptor (GR) DNA binding domain consists of several conserved amino acids and folds into two zinc finger-like structures. Previous transactivation experiments indicated that three amino acids residing in this region, Gly, Ser and Val, appear to be critical for target-site discrimination. Based on the solved crystal structure, these residues are at the beginning of an amphipathic alpha-helix that interacts with the DNA's major groove; of these, only valine, however, contacts DNA. In order to examine their functional role directly, we have substituted these residues for the corresponding amino acids from the estrogen receptor (ER), overexpressed and purified the mutant proteins, and assayed their binding specificity and affinity by gel mobility shifts using glucocorticoid or estrogen response elements (GRE or ERE, respectively) as DNA probes. We find that all three residues are indeed required to fully switch GR's specificity to an ERE. The contacting valine in GR is of primary importance. The corresponding residue in ER, alanine, is less important for specificity, while glutamic acid, four amino acids towards the N-terminus, is most critical for ER discrimination. Finally, we show that the GR DNA binding domain carrying all three ER-specific mutations has a significantly higher affinity for an ERE than the ER DNA binding domain itself. We interpret these results in the context of both the data presented here and the crystal structure of the GR DNA binding domain complexed to a GRE.     

Efficient prediction of nucleic acid binding function from low-resolution protein structures

Structural genomics projects as well as ab initio protein structure prediction methods provide structures of proteins with no sequence or fold similarity to proteins with known functions. These are often low-resolution structures that may only include the positions of C alpha atoms. We present a fast and efficient method to predict DNA-binding proteins from just the amino acid sequences and low-resolution, C alpha-only protein models. The method uses the relative proportions of certain amino acids in the protein sequence, the asymmetry of the spatial distribution of certain other amino acids as well as the dipole moment of the molecule. These quantities are used in a linear formula, with coefficients derived from logistic regression performed on a training set, and DNA-binding is predicted based on whether the result is above a certain threshold. We show that the method is insensitive to errors in the atomic coordinates and provides correct predictions even on inaccurate protein models. We demonstrate that the method is capable of predicting proteins with novel binding site motifs and structures solved in an unbound state. The accuracy of our method is close to another, published method that uses all-atom structures, time-consuming calculations and information on conserved residues.     

Triggering loops and enzyme function: identification of loops that trigger and modulate movements

Enzyme function often involves a conformational change. There is a general agreement that loops play a vital role in correctly positioning the catalytically important residues. Nevertheless, predicting the functional loops and most importantly their role in enzyme function remains a difficult task. A major reason for this difficulty is that loops that undergo conformational change are frequently not well conserved in their primary sequence. beta1,4-Galactosyltransferase is one such enzyme. There, the amino acid sequence of a long loop that undergoes a large conformational change upon substrate binding is not well conserved. Our molecular dynamics simulations show that the large conformational change in the long loop is brought about by a second, interacting loop. Interestingly, while the structural change of the second loop is much smaller than that of the long loop, its sequence (particularly glycine residues) is highly conserved. We further examine the generality of the proposition that there are loops that trigger movements but nevertheless show little or no structural changes in crystals. We focus on two other enzymes, enolase and lipase. We chose these enzymes, since they too undergo conformational change upon ligand binding, however, they have different folds and different functions. Through multiple sets of simulations we show that the conformational change of the functional loop(s) is brought about through communication of flexibility by triggering loops that have several glycine residues. We further propose that similar to the conservation of common favorable fold types and structural motifs, evolution has also conserved common "skillful" mechanisms. Mechanisms may be conserved across different folds, sequences and functions, with adaptation to specific enzymatic roles.     

Helicase motifs: the engine that powers DNA unwinding

Helicases play essential roles in nearly all DNA metabolic transactions and have been implicated in a variety of human genetic disorders. A hallmark of these enzymes is the existence of a set of highly conserved amino acid sequences termed the 'helicase motifs' that were hypothesized to be critical for helicase function. These motifs are shared by another group of enzymes involved in chromatin remodelling. Numerous structure-function studies, targeting highly conserved residues within the helicase motifs, have been instrumental in uncovering the functional significance of these regions. Recently, the results of these mutational studies were augmented by the solution of the three-dimensional crystal structure of three different helicases. The structural model for each helicase revealed that the conserved motifs are clustered together, forming a nucleotide-binding pocket and a portion of the nucleic acid binding site. This result is gratifying, as it is consistent with structure-function studies suggesting that all the conserved motifs are involved in the nucleotide hydrolysis reaction. Here, we review helicase structure-function studies in the light of the recent crystal structure reports. The current data support a model for helicase action in which the conserved motifs define an engine that powers the unwinding of duplex nucleic acids, using energy derived from nucleotide hydrolysis and conformational changes that allow the transduction of energy between the nucleotide and nucleic acid binding sites. In addition, this ATP-hydrolysing engine is apparently also associated with proteins involved in chromatin remodelling and provides the energy required to alter protein-DNA structure, rather than duplex DNA or RNA structure.     

Residues critical for retroviral integrative recombination in a region that is highly conserved among retroviral/retrotransposon integrases and bacterial insertion sequence transposases.

Our comparison of deduced amino acid sequences for retroviral/retrotransposon integrase (IN) proteins of several organisms, including Drosophila melanogaster and Schizosaccharomyces pombe, reveals strong conservation of a constellation of amino acids characterized by two invariant aspartate (D) residues and a glutamate (E) residue, which we refer to as the D,D(35)E region. The same constellation is found in the transposases of a number of bacterial insertion sequences. The conservation of this region suggests that the component residues are involved in DNA recognition, cutting, and joining, since these properties are shared among these proteins of divergent origin. We introduced amino acid substitutions in invariant residues and selected conserved and nonconserved residues throughout the D,D(35)E region of Rous sarcoma virus IN and in human immunodeficiency virus IN and assessed their effect upon the activities of the purified, mutant proteins in vitro. Changes of the invariant and conserved residues typically produce similar impairment of both viral long terminal repeat (LTR) oligonucleotide cleavage referred to as the processing reaction and the subsequent joining of the processed LTR-based oligonucleotides to DNA targets. The severity of the defects depended upon the site and the nature of the amino acid substitution(s). All substitutions of the invariant acidic D and E residues in both Rous sarcoma virus and human immunodeficiency virus IN dramatically reduced LTR oligonucleotide processing and joining to a few percent or less of wild type, suggesting that they are essential components of the active site for both reactions.(ABSTRACT TRUNCATED AT 250 WORDS)     

Structural basis for chitotetraose coordination by CGL3, a novel galectin-related protein from Coprinopsis cinerea

Recent advances in genome sequencing efforts have revealed an abundance of novel putative lectins. Among these, many galectin-related proteins, characterized by many conserved residues but intriguingly lacking critical amino acids, have been found in all corners of the eukaryotic superkingdom. Here we present a structural and biochemical analysis of one representative, the galectin-related lectin CGL3 found in the inky cap mushroom Coprinopsis cinerea. This protein contains all but one conserved residues known to be involved in β-galactoside binding in galectins. A Trp residue strictly conserved among galectins is changed to an Arg in CGL3 (R81). Accordingly, the galectin-related protein is not able to bind lactose. Screening of a glycan array revealed that CGL3 displays preference for oligomers of β1-4-linked N-acetyl-glucosamines (chitooligosaccharides) and GalNAcβ1-4GlcNAc (LacdiNAc). Carbohydrate-binding affinity of this novel lectin was quantified using isothermal titration calorimetry, and its mode of chitooligosaccharide coordination not involving any aromatic amino acid residues was studied by X-ray crystallography. Structural information was used to alter the carbohydrate-binding specificity and substrate affinity of CGL3. The importance of residue R81 in determining the carbohydrate-binding specificity was demonstrated by replacing this Arg with a Trp residue (R81W). This single-amino-acid change led to a lectin that failed to bind chitooligosaccharides but gained lactose binding. Our results demonstrate that, similar to the legume lectin fold, the galectin fold represents a conserved structural framework upon which dramatically altered specificities can be grafted by few alterations in the binding site and that, in consequence, many metazoan galectin-related proteins may represent lectins with novel carbohydrate-binding specificities.     

Conserved repeat motifs and glucan binding by glucansucrases of oral streptococci and Leuconostoc mesenteroides

Glucansucrases of oral streptococci and Leuconostoc mesenteroides have a common pattern of structural organization and characteristically contain a domain with a series of tandem amino acid repeats in which certain residues are highly conserved, particularly aromatic amino acids and glycine. In some glucosyltransferases (GTFs) the repeat region has been identified as a glucan binding domain (GBD). Such GBDs are also found in several glucan binding proteins (GBP) of oral streptococci that do not have glucansucrase activity. Alignment of the amino acid sequences of 20 glucansucrases and GBP showed the widespread conservation of the 33-residue A repeat first identified in GtfI of Streptococcus downei. Site-directed mutagenesis of individual highly conserved residues in recombinant GBD of GtfI demonstrated the importance of the first tryptophan and the tyrosine-phenylalanine pair in the binding of dextran, as well as the essential contribution of a basic residue (arginine or lysine). A microplate binding assay was developed to measure the binding affinity of recombinant GBDs. GBD of GtfI was shown to be capable of binding glucans with predominantly alpha-1,3 or alpha-1,6 links, as well as alternating alpha-1,3 and alpha-1,6 links (alternan). Western blot experiments using biotinylated dextran or alternan as probes demonstrated a difference between the binding of streptococcal GTF and GBP and that of Leuconostoc glucansucrases. Experimental data and bioinformatics analysis showed that the A repeat motif is distinct from the 20-residue CW motif, which also has conserved aromatic amino acids and glycine and which occurs in the choline-binding proteins of Streptococcus pneumoniae and other organisms.     

Structural complementarity of distance constraints obtained from chemical cross-linking and amino acid coevolution

The analysis of amino acid coevolution has emerged as a practical method for protein structural modeling by providing structural contact information from alignments of amino acid sequences. In parallel, chemical cross-linking/mass spectrometry (XLMS) has gained attention as a universally applicable method for obtaining low-resolution distance constraints to model the quaternary arrangements of proteins, and more recently even protein tertiary structures. Here, we show that the structural information obtained by XLMS and coevolutionary analysis are effectively complementary: the distance constraints obtained by each method are almost exclusively associated with non-coincident pairs of residues, and modeling results obtained by the combination of both sets are improved relative to considering the same total number of constraints of a single type. The structural rationale behind the complementarity of the distance constraints is discussed and illustrated for a representative set of proteins with different sizes and folds.     

Evolution of amino acid frequencies in proteins over deep time: inferred order of introduction of amino acids into the genetic code

To understand more fully how amino acid composition of proteins has changed over the course of evolution, a method has been developed for estimating the composition of proteins in an ancestral genome. Estimates are based upon the composition of conserved residues in descendant sequences and empirical knowledge of the relative probability of conservation of various amino acids. Simulations are used to model and correct for errors in the estimates. The method was used to infer the amino acid composition of a large protein set in the Last Universal Ancestor (LUA) of all extant species. Relative to the modern protein set, LUA proteins were found to be generally richer in those amino acids that are believed to have been most abundant in the prebiotic environment and poorer in those amino acids that are believed to have been unavailable or scarce. It is proposed that the inferred amino acid composition of proteins in the LUA probably reflects historical events in the establishment of the genetic code.     

Nearest-neighbor classifier as a tool for classification of protein families

Knowledge about protein function is essential in understanding the biological processes. A specific class or family of protein shares common structural and chemical properties amongst its member sequences. The set of properties that display its unique characteristics for clearly classifying a protein sequence into its corresponding protein family needs to be studied. Our study of these important properties conducted on four major classes of proteins namely Globins, Homeoboxes, Heat Shock proteins (HSP) and Kinase have shown that frequency of twenty naturally occurring amino acids, hydrophobic content of protein, molecular weight of protein, isoelectric point of protein, secondary structure composition of amino acid residues as helices, coils and sheets and the composition of helices, coils and sheets in the secondary structure topology plays a significant role in correctly classifying the protein into its corresponding class or family as indicated by the overall efficiency of Nearest Neighbor Classifier as 84.92%.     

Amino acids conserved in interleukin-1 receptors (IL-1Rs) and the Drosophila toll protein are essential for IL-1R signal transduction.

The cytoplasmic domain of the human T cell-type interleukin-1 receptor (hIL-1R) is not involved in the binding, internalization, or nuclear localization of interleukin-1 (IL-1), but is essential for signal transduction. We have previously localized a 50-amino acid region (residues 477-527) critical for IL-1-mediated activation of the interleukin-2 promoter in T cells. This region displays a striking degree of amino acid conservation in human, murine, and chicken IL-1Rs. Here we report the results of a site-directed mutational analysis of the cytoplasmic domain of the hIL-1R. We have introduced single-amino acid substitutions at positions conserved in all three receptors and at nonconserved positions and identified key amino acids for IL-1R function in signal transduction. Three basic (Arg431, Lys515, and Arg518) and 3 aromatic (Phe513, Trp514, and Tyr519) amino acids that are conserved in human, murine, and chicken IL-1Rs could not be replaced without abolishing IL-1R-mediated signal transduction. A substitution at another conserved position (Pro521) reduces significantly the ability of the IL-1R to transmit the IL-1 signal. Nonconserved residues could be replaced without affecting signal transduction. The cytoplasmic domain of the IL-1R is related to that of the Drosophila Toll protein, with a 26% identity and a 43% similarity in amino acid sequence. The amino acids shown to be essential for IL-1R function are conserved in the Toll protein. Our experimental data indicate that the amino acid sequence similarity between the IL-1R and the Drosophila toll protein reflects a functional homology between the two proteins.