Computational design of Candida boidinii xylose reductase for altered cofactor specificity

In this study we introduce a computationally‐driven enzyme redesign workflow for altering cofactor specificity from NADPH to NADH. By compiling and comparing data from previous studies involving cofactor switching mutations, we show that their effect cannot be explained as straightforward changes in volume, hydrophobicity, charge, or BLOSUM62 scores of the residues populating the cofactor binding site. Instead, we find that the use of a detailed cofactor binding energy approximation is needed to adequately capture the relative affinity towards different cofactors. The implicit solvation models Generalized Born with molecular volume integration and Generalized Born with simple switching were integrated in the iterative protein redesign and optimization (IPRO) framework to drive the redesign of Candida boidinii xylose reductase (CbXR) to function using the non‐native cofactor NADH. We identified 10 variants, out of the 8,000 possible combinations of mutations, that improve the computationally assessed binding affinity for NADH by introducing mutations in the CbXR binding pocket. Experimental testing revealed that seven out of ten possessed significant xylose reductase activity utilizing NADH, with the best experimental design (CbXR‐GGD) being 27‐fold more active on NADH. The NADPH‐dependent activity for eight out of ten predicted designs was either completely abolished or significantly diminished by at least 90%, yielding a greater than 10⁴‐fold change in specificity to NADH (CbXR‐REG). The remaining two variants (CbXR‐RTT and CBXR‐EQR) had dual cofactor specificity for both nicotinamide cofactors.     

Shape-specific nucleotide binding of single-stranded RNA by the GLD-1 STAR domain

Proteins containing the STAR RNA-binding domain fulfill vital roles in RNA biogenesis, yet a detailed understanding of STAR domain RNA binding specificity is lacking. In Caenorhabditis elegans, the STAR protein GLD-1 directly binds the 28 nucleotide recognition element TGE within the 3' untranslated region of tra-2 mRNA. The GLD-1:TGE interaction promotes translational silencing of tra-2 mRNA, marking a pivotal event in the spermatogenesis to oogenesis switch in C.elegans hermaphrodites. By measuring the binding affinities of both GLD-1 and TGE mutants, we have explored the molecular determinants of STAR domain specificity. Site-directed GLD-1 mutants were guided by sequence homology with human splicing factor 1 (SF1), for which an RNA:protein complex structure is available in the work done by Liu et al. The RNA binding affinity of 11 mutant GLD-1 proteins was measured, and their binding specificity was assessed with a series of TGE RNAs containing natural or modified nucleotides. This combinatorial analysis of both RNA and protein mutants revealed a diverse array of specificities of individual nucleotide-binding pockets along the interface. At nucleotide position 18, adenosine appears to be specified by the overall shape of a pocket lined with aliphatic side-chains. At position 19, the high preference for cytidine is dependent on both the length of an amino acid side-chain and the identity of terminal functional groups. The nucleotide 21 binding pocket exhibits low discrimination for cytidine, and accommodates most nucleobases. The highly hydrophobic binding interface and apparent small number of hydrogen bonding read-out interactions at these positions is consistent with our finding that few amino acids seem to function individually in establishing binding specificity. Rather, specificity is conferred by the shape of the nucleotide-binding pocket. Our data provide the first detailed, quantitative analysis of the STAR domain, and highlight features of STAR:RNA recognition that are distinct among single-stranded RNA-binding proteins.     

Change in protein-ligand specificity through binding pocket grafting

Recognition and discrimination of small molecules are crucial for biological processes in living systems. Understanding the mechanisms that underlie binding specificity is of particular interest to synthetic biology, e.g. the engineering of biosensors with de novo ligand affinities. Promising scaffolds for such biosensors are the periplasmic binding proteins (PBPs) due to their ligand-mediated structural change that can be translated into a physically measurable signal.In this study we focused on the two homologous polyamine binding proteins PotF and PotD. Despite their structural similarity, PotF and PotD have different binding specificities for the polyamines putrescine and spermidine. To elucidate how specificity is determined, we grafted the binding site of PotD onto PotF. The introduction of 7 mutations in the first shell of the binding pocket leads to a swap in the binding profile as confirmed by isothermal titration calorimetry. Furthermore, the 1.7 Å crystal structure of the new variant complexed with spermidine reveals the interactions of the specificity determining residues including a defined water network. Altogether our study shows that specificity is encoded in the first shell residues of the PotF binding pocket and that transplantation of these residues allows the swap of the binding specificity.     

Deciphering metal ion preference and primary coordination sphere robustness of a designed zinc finger with high‐resolution mass spectrometry

Small zinc finger (ZnF) motifs are promising molecular scaffolds for protein design owing to their structural robustness and versatility. Moreover, their characterization provides important insights into protein folding in general. ZnF motifs usually possess an exceptional specificity and high affinity towards Zn(II) ion to drive folding. While the Zn(II) ion is canonically coordinated by two cysteine and two histidine residues, many other coordination spheres also exist in small ZnFs, all having four amino acid ligands. Here we used high‐resolution mass spectrometry to study metal ion binding specificity and primary coordination sphere robustness of a designed zinc finger, named MM1. Based on the results, MM1 possesses high specificity for zinc with sub‐micromolar binding affinity. Surprisingly, MM1 retains metal ion binding affinity even in the presence of selective alanine mutations of the primary zinc coordinating amino acid residues.     

Computational active site analysis of molecular pathways to improve functional classification of enzymes

This study describes a method to computationally assess the function of homologous enzymes through small molecule binding interaction energy. Three experimentally determined X-ray structures and four enzyme models from ornithine cyclo-deaminase, alanine dehydrogenase, and mu-crystallin were used in combination with nine small molecules to derive a function score (FS) for each enzyme-model combination. While energy values varied for a single molecule-enzyme combination due to differences in the active sites, we observe that the binding energies for the entire pathway were proportional for each set of small molecules investigated. This proportionality of energies for a reaction pathway appears to be dependent on the amino acids in the active site and their direct interactions with the small molecules, which allows a function score (FS) to be calculated to assess the specificity of each enzyme. Potential of mean force (PMF) calculations were used to obtain the energies, and the resulting FS values demonstrate that a measurement of function may be obtained using differences between these PMF values. Additionally, limitations of this method are discussed based on: (a) larger substrates with significant conformational flexibility; (b) low homology enzymes; and (c) open active sites. This method should be useful in accurately predicting specificity for single enzymes that have multiple steps in their reactions and in high throughput computational methods to accurately annotate uncharacterized proteins based on active site interaction analysis.     

Xenopus embryonic poly(A) binding protein 2 (ePABP2) defines a new family of cytoplasmic poly(A) binding proteins expressed during the early stages of vertebrate development

We describe a new RNA binding protein from Xenopus we have named ePABP2 (embryonic poly(A) binding protein type II). Based on amino acid similarity, ePABP2 is closely related to the ubiquitously expressed nuclear PABP2 protein that directs the elongation of mRNA poly(A) tails during pre-mRNA processing. However, in contrast to known PABP2 proteins, Xenopus ePABP2 is a cytoplasmic protein that is predominantly expressed during the early stages of Xenopus development and in adult ovarian tissue. Biochemical experiments indicate ePABP2 binds poly(A) with specificity and that this binding requires the RRM domain. Mouse and human ePABP2 proteins were also identified and mouse ePABP2 expression is also confined to the earliest stages of mouse development and adult ovarian tissue. We propose that Xenopus ePABP2 is the founding member of a new class of poly(A) binding proteins expressed in vertebrate embryos. Possible roles for this protein in regulating mRNA function in early vertebrate development are discussed.     

Crucial Roles of Single Residues in Binding Affinity,Specificity, and Promiscuity in the Cellulosomal Cohesin-Dockerin Interface

Interactions between cohesin and dockerin modules play a crucial role in the assembly of multienzyme cellulosome complexes. Although intraspecies cohesin and dockerin modules bind in general with high affinity but indiscriminately, cross-species binding is rare. Here, we combined ELISA-based experiments with Rosetta-based computational design to evaluate the contribution of distinct residues at the Clostridium thermocellum cohesin-dockerin interface to binding affinity, specificity, and promiscuity. We found that single mutations can show distinct and significant effects on binding affinity and specificity. In particular, mutations at cohesin position Asn³⁷ show dramatic variability in their effect on dockerin binding affinity and specificity: the N37A mutant binds promiscuously both to cognate (C. thermocellum) as well as to non-cognate Clostridium cellulolyticum dockerin. N37L in turn switches binding specificity: compared with the wild-type C. thermocellum cohesin, this mutant shows significantly increased preference for C. cellulolyticum dockerin combined with strongly reduced binding to its cognate C. thermocellum dockerin. The observation that a single mutation can overcome the naturally observed specificity barrier provides insights into the evolutionary dynamics of this system that allows rapid modulation of binding specificity within a high affinity background.     

Intracellular cholesterol transporter StarD4 binds free cholesterol and increases cholesteryl ester formation

StarD4 protein is a member of the StarD4 subfamily of steroidogenic acute regulatory-related lipid transfer (START) domain proteins that includes StarD5 and StarD6, proteins whose functions remain poorly defined. The objective of this study was to isolate and characterize StarD4's sterol binding and to determine in a hepatocyte culture model its sterol transport capabilities. Utilizing purified full-length StarD4, in vitro binding assays demonstrated a concentration-dependent binding of [(14)C]cholesterol by StarD4 similar to that of the cholesterol binding START domain proteins StarD1 and StarD5. Other tested sterols showed no detectable binding to StarD4, except for 7alpha-hydroxycholesterol, for which StarD4 demonstrated weak binding on lipid protein overlay assays. Subsequently, an isolated mouse hepatocyte model was used to study the ability of StarD4 to bind/mobilize/distribute cellular cholesterol. Increased expression of StarD4 in primary mouse hepatocytes led to a marked increase in the intracellular cholesteryl ester concentration and in the rates of bile acid synthesis. The ability and specificity of StarD4 to bind cholesterol and, as a function of its level of expression, to direct endogenous cellular cholesterol suggest that StarD4 plays an important role as a directional cholesterol transporter in the maintenance of cellular cholesterol homeostasis.     

Fungal cell wall phosphomannans facilitate the toxic activity of a plant PR-5 protein

Osmotin is a plant PR-5 protein. It has a broad spectrum of antifungal activity, yet also exhibits specificity for certain fungal targets. The structural bases for this specificity remain unknown. We show here that full sensitivity of Saccharomyces cerevisiae cells to the PR-5 protein osmotin is dependent on the function of MNN2, MNN4 and MNN6. MNN2 is an alpha-1, 2-mannosyltransferase catalyzing the addition of the first mannose to the branches on the poly l,6-mannose backbone of the outer chain of cell wall N-linked mannans. MNN4 and MNN6 are required for the transfer of mannosylphosphate to cell wall mannans. Null mnn2, mnn4 or mnn6 mutants lack phosphomannans and are defective in binding osmotin to the fungal cell wall. Both antimannoprotein antibody and the cationic dye alcian blue protect cells against osmotin cytotoxicity. MNN1 is an alpha-1,3-mannosyltransferase that adds the terminal mannose to the outer chain branches of N-linked mannan, masking mannosylphosphate. Null mnn1 cells exhibit enhanced osmotin binding and sensitivity. Several cell wall mannoproteins can bind to immobilized osmotin, suggesting that their polysaccharide constituent determines osmotin binding. Our results demonstrating a causal relationship between cell surface phosphomannan and the susceptibility of a yeast strain to osmotin suggest that cell surface polysaccharides of invading pathogens control target specificity of plant PR-5 proteins.     

The Tiam1 PDZ Domain Couples to Syndecan1 and Promotes Cell-Matrix Adhesion

The T-cell lymphoma invasion and metastasis gene 1 (Tiam1) is a guanine exchange factor (GEF) for the Rho-family GTPase Rac1 that is crucial for the integrity of adherens junctions, tight junctions, and cell-matrix interactions. This GEF contains several protein-protein interaction domains, including a PDZ domain. Earlier studies identified a consensus PDZ-binding motif and a synthetic peptide capable of binding to the Tiam1 PDZ domain, but little is known about its ligand specificity and physiological role in cells. Here, we investigated the structure, specificity, and function of the Tiam1 PDZ domain. We determined the crystal structures of the Tiam1 PDZ domain free and in complex with a “model” peptide, which revealed the structural basis for ligand specificity. Protein database searches using the consensus PDZ-binding motif identified two eukaryotic cell adhesion proteins, Syndecan1 and Caspr4, as potential Tiam1 PDZ domain binding proteins. Equilibrium binding experiments confirmed that C-terminal peptides derived from Syndecan1 and Caspr4 bound the Tiam1 PDZ domain. NMR chemical shift perturbation experiments indicated that the Tiam1 PDZ/Syndecan1 and PDZ/Caspr4 complexes were structurally distinct and identified key residues likely to be responsible for ligand selectivity. Moreover, cell biological analysis established that Syndecan1 is a physiological binding partner of Tiam1 and that the PDZ domain has a function in cell-matrix adhesion and cell migration. Collectively, our data provide insight into the structure, specificity, and function of the Tiam1 PDZ domain. Importantly, our data report on a physiological role for the Tiam1 PDZ domain and establish a novel link between two previously unrelated signal transduction pathways, both of which are implicated in cancer.     

The structural basis of receptor-binding by Escherichia coli associated with diarrhea and septicemia

GafD in Escherichia coli G (F17) fimbriae is associated with diarrheal disease, and the structure of the ligand-binding domain, GafD1-178, has been determined at 1.7A resolution in the presence of the receptor sugar N-acetyl-D-glucosamine. The overall fold is a beta-barrel jelly-roll fold. The ligand-binding site was identified and localized to the side of the molecule. Receptor binding is mediated by side-chain as well main-chain interactions. Ala43-Asn44, Ser116-Thr117 form the sugar acetamide specificity pocket, while Asp88 confers tight binding and Trp109 appears to position the ligand. There is a disulfide bond that rigidifies the acetamide specificity pocket. The three fimbrial lectins, GafD, FimH and PapG share similar beta-barrel folds but display different ligand-binding regions and disulfide-bond patterns. We suggest an evolutionary path for the evolution of the very diverse fimbrial lectins from a common ancestral fold.     

Monobodies and other synthetic binding proteins for expanding protein science

Synthetic binding proteins are constructed using nonantibody molecular scaffolds. Over the last two decades, in‐depth structural and functional analyses of synthetic binding proteins have improved combinatorial library designs and selection strategies, which have resulted in potent platforms that consistently generate binding proteins to diverse targets with affinity and specificity that rival those of antibodies. Favorable attributes of synthetic binding proteins, such as small size, freedom from disulfide bond formation and ease of making fusion proteins, have enabled their unique applications in protein science, cell biology and beyond. Here, we review recent studies that illustrate how synthetic binding proteins are powerful probes that can directly link structure and function, often leading to new mechanistic insights. We propose that synthetic proteins will become powerful standard tools in diverse areas of protein science, biotechnology and medicine.     

Efficient T-cell receptor signaling requires a high-affinity interaction between the Gads C-SH3 domain and the SLP-76 RxxK motif

The relationship between the binding affinity and specificity of modular interaction domains is potentially important in determining biological signaling responses. In signaling from the T-cell receptor (TCR), the Gads C-terminal SH3 domain binds a core RxxK sequence motif in the SLP-76 scaffold. We show that residues surrounding this motif are largely optimized for binding the Gads C-SH3 domain resulting in a high-affinity interaction (K(D)=8-20 nM) that is essential for efficient TCR signaling in Jurkat T cells, since Gads-mediated signaling declines with decreasing affinity. Furthermore, the SLP-76 RxxK motif has evolved a very high specificity for the Gads C-SH3 domain. However, TCR signaling in Jurkat cells is tolerant of potential SLP-76 crossreactivity, provided that very high-affinity binding to the Gads C-SH3 domain is maintained. These data provide a quantitative argument that the affinity of the Gads C-SH3 domain for SLP-76 is physiologically important and suggest that the integrity of TCR signaling in vivo is sustained both by strong selection of SLP-76 for the Gads C-SH3 domain and by a capacity to buffer intrinsic crossreactivity.     

Mutation of active site residues in the chitin-binding domain ChBDChiA1 from chitinase A1 of Bacillus circulans alters substrate specificity: use of a green fluorescent protein binding assay

A fluorescent binding assay was developed to investigate the effects of mutagenesis on the binding affinity and substrate specificity of the chitin-binding domain of chitinase A1 from Bacillus circulans WL-12. The chitin-binding domain was genetically fused to the N-terminus of a green fluorescent protein, and the polyhistidine-tagged hybrid protein was expressed in Escherichia coli. Residues likely to be involved in the binding site were mutated and their contributions to binding and substrate specificity were evaluated by affinity electrophoresis and depletion assays. The experimental binding isotherms were analyzed by non-linear regression using a modified Langmuir equation. Non-conservative substitution of tryptophan residue (W687) nearly abolished chitin-binding affinity and dramatically lowered chitosan binding while retaining the original level of curdlan binding. Double mutation E668K/P689A had altered specificity for several substrates and also impaired chitin binding significantly. Other substitutions in the binding site altered substrate specificity but had little effect on overall affinity for chitin. Interestingly, mutation T682A led to a higher specificity towards chitinous substrates than the wildtype. Furthermore, the ChBD-GFP hybrid protein was tested for use in diagnostic staining of cell walls of fungi and yeast and for the detection of fungal infections in tissue samples.     

Efficient recognition of protein fold at low sequence identity by conservative application of Psi-BLAST: application

Based on a study involving structural comparisons of proteins sharing 25% or less sequence identity, three rounds of Psi-BLAST appear capable of identifying remote evolutionary homologs with greater than 95% confidence provided that more than 50% of the query sequence can be aligned with the target sequence. Since it seems that more than 80% of all homologous protein pairs may be characterized by a lack of significant sequence similarity, the experimental biologist is often confronted with a lack of guidance from conventional homology searches involving pair-wise sequence comparisons. The ability to disregard levels of sequence identity and expect value in Psi-BLAST if at least 50% of the query sequence has been aligned allows for generation of new hypotheses by consideration of matches that are conventionally disregarded. In one example, we suggest a possible evolutionary linkage between the cupredoxin and immunoglobulin fold families. A thermostable hypothetical protein of unknown function may be a circularly permuted homolog to phosphotriesterase, an enzyme capable of detoxifying organophosphate nerve agents. In a third example, the amino acid sequence of another hypothetical protein of unknown function reveals the ATP binding-site, metal binding site, and catalytic sidechain consistent with kinase activity of unknown specificity. This approach significantly expands the utility of existing sequence data to define the primary structure degeneracy of binding sites for substrates, cofactors and other proteins.