A Pro to Gly mutation in the hinge of the arabinose-binding protein enhances binding and alters specificity. Sugar-binding and crystallographic studies

Laboratory Evolution and Reverse Engineering of Clostridium thermocellum for Growth on Glucose and Fructose

The native ability of Clostridium thermocellum to efficiently solubilize cellulose makes it an interesting platform for sustainable biofuel production through consolidated bioprocessing. Together with other improvements, industrial implementation of C. thermocellum, as well as fundamental studies into its metabolism, would benefit from improved and reproducible consumption of hexose sugars. To investigate growth of C. thermocellum on glucose or fructose, as well as the underlying molecular mechanisms, laboratory evolution was performed in carbon-limited chemostats with increasing concentrations of glucose or fructose and decreasing cellobiose concentrations. Growth on both glucose and fructose was achieved with biomass yields of 0.09 ± 0.00 and 0.18 ± 0.00 gbiomass gsubstrate-1, respectively, compared to 0.15 ± 0.01 gbiomass gsubstrate-1 for wild type on cellobiose. Single-colony isolates had no or short lag times on the monosaccharides, while wild type showed 42 ± 4 h on glucose and >80 h on fructose. With good growth on glucose, fructose, and cellobiose, the fructose isolates were chosen for genome sequence-based reverse metabolic engineering. Deletion of a putative transcriptional regulator (Clo1313_1831), which upregulated fructokinase activity, reduced lag time on fructose to 12 h with a growth rate of 0.11 ± 0.01 h-1 and resulted in immediate growth on glucose at 0.24 ± 0.01 h-1 Additional introduction of a G-to-V mutation at position 148 in cbpA resulted in immediate growth on fructose at 0.32 ± 0.03 h-1 These insights can guide engineering of strains for fundamental studies into transport and the upper glycolysis, as well as maximizing product yields in industrial settings.IMPORTANCEC. thermocellum is an important candidate for sustainable and cost-effective production of bioethanol through consolidated bioprocessing. In addition to unsurpassed cellulose deconstruction, industrial application and fundamental studies would benefit from improvement of glucose and fructose consumption. This study demonstrated that C. thermocellum can be evolved for reproducible constitutive growth on glucose or fructose. Subsequent genome sequencing, gene editing, and physiological characterization identified two underlying mutations with a role in (regulation of) transport or metabolism of the hexose sugars. In light of these findings, such mutations have likely (and unknowingly) also occurred in previous studies with C. thermocellum using hexose-based media with possible broad regulatory consequences. By targeted modification of these genes, industrial and research strains of C. thermocellum can be engineered to (i) reduce glucose accumulation, (ii) study cellodextrin transport systems in vivo, (iii) allow experiments at >120 g liter-1 soluble substrate concentration, or (iv) reduce costs for labeling studies.

Contribution of buried distal amino acid residues in horse liver alcohol dehydrogenase to structure and catalysis

The dynamics of enzyme catalysis range from the slow time scale (∼ms) for substrate binding and conformational changes to the fast time (∼ps) scale for reorganization of substrates in the chemical step. The contribution of global dynamics to catalysis by alcohol dehydrogenase was tested by substituting five different, conserved amino acid residues that are distal from the active site and located in the hinge region for the conformational change or in hydrophobic clusters. X-ray crystallography shows that the structures for the G173A, V197I, I220 (V, L, or F), V222I, and F322L enzymes complexed with NAD⁺ and an analogue of benzyl alcohol are almost identical, except for small perturbations at the sites of substitution. The enzymes have very similar kinetic constants for the oxidation of benzyl alcohol and reduction of benzaldehyde as compared to the wild-type enzyme, and the rates of conformational changes are not altered. Less conservative substitutions of these amino acid residues, such as G173(V, E, K, or R), V197(G, S, or T), I220(G, S, T, or N), and V222(G, S, or T) produced unstable or poorly expressed proteins, indicating that the residues are critical for global stability. The enzyme scaffold accommodates conservative substitutions of distal residues, and there is no evidence that fast, global dynamics significantly affect the rate constants for hydride transfers. In contrast, other studies show that proximal residues significantly participate in catalysis.

A new role for Escherichia coli DsbC protein in protection against oxidative stress

We report a new function for Escherichia coli DsbC, a protein best known for disulfide bond isomerization in the periplasm. We found that DsbC regulates the redox state of the single cysteine of the L-arabinose-binding protein AraF. This cysteine, which can be oxidized to a sulfenic acid, mediates the formation of a disulfide-linked homodimer under oxidative stress conditions, preventing L-arabinose binding. DsbC, unlike the homologous protein DsbG, reduces the intermolecular disulfide, restoring AraF binding properties. Thus, our results reveal a new link between oxidative protein folding and the defense mechanisms against oxidative stress.

Allosteric signaling and dynamics of the clamshell-like NMDA receptor GluN1 N-terminal domain

N-methyl-D-aspartate receptors (NMDARs), neuronal glutamate-gated ion channels, are obligatory heterotetramers composed of GluN1 and GluN2 subunits. Each subunit contains two extracellular clamshell-like domains with an agonist-binding domain and a distal N-terminal domain (NTD). The GluN2 NTDs form mobile regulatory domains. In contrast, the dynamics of GluN1 NTD and its contribution to NMDAR function remain poorly understood. Here we show that GluN1 NTD is neither static nor functionally silent. Perturbing the conformation of GluN1 NTD affects both receptor gating and pharmacological properties. GluN1 NTD undergoes structural rearrangements that involve hinge bending and large twisting and untwisting motions, allowing for new intra- and intersubunit contacts. GluN1 NTD acts in trans with GluN2 NTD to influence binding of glutamate but, notably, not of GluN1 coagonist glycine. Our work uncovers a dynamic role of GluN1 NTD in controlling NMDAR function through new interdomain allosteric interactions.

Insights derived from molecular dynamics simulation into the molecular motions of serine protease proteinase K

Serine protease proteinase K, a member of the subtilisin family of enzymes, is of significant industrial, agricultural and biotechnological importance. Despite the wealth of structural information about proteinase K provided by static X-ray structures, a full understanding of the enzymatic mechanism requires further insight into the dynamic properties of this enzyme. Molecular dynamics simulations and essential dynamics (ED) analysis were performed to investigate the molecular motions in proteinase K. The results indicate that the internal core of proteinase K is relatively rigid, whereas the surface-exposed loops, most notably the substrate-binding regions, exhibit considerable conformational fluctuations. Further ED analysis reveals that the large concerted motions in the substrate-binding regions cause opening/closing of the substrate-binding pockets, thus supporting the proposed induced-fit mechanism of substrate binding. The distinct electrostatic/hydrogen-bonding interactions between Asp39 and His69 and between His69 and Ser224 within the catalytic triad lead to different thermal motions and orientations of these three catalytic residues, which can be related to their different functional roles in the catalytic process. Statistical analyses of the geometrical/functional properties as well as evolutionary conservation of the glycines in proteinase K-like proteins reveal that glycines may play an important role in determining the folding architecture and structural flexibility of this class of enzymes. Our simulation study complements the biochemical and structural studies and provides new insights into the dynamic structural basis of the functional properties of this class of enzymes.

Biological roles of specific peptides in enzymes

It has recently been shown (Kunik et al., PLOS Comput Biol 2007;3(8):e167) that the occurrence of specific peptides (SPs) on sequences of enzymes allows for accurate EC classification of enzymes. We inquire whether these SPs play important roles in bringing about the enzymatic function. This is assessed by cross-checking the occurrence of SPs on enzymes with Swiss-Prot annotations and PDB spatial structures of enzymes. Analyzing the coverage of functional annotations of enzymes, we demonstrate that SPs contain major fractions of all annotated features. This result is statistically highly significant and associates over 10% of all SPs with important biological markers. Concentrating on DNA binding regions, relevant to LexA repressor enzymes, we find interesting coverage patterns. Moreover, for the same data, we demonstrate that SPs allow for subclassification of the relevant bacteria into phylogenetic classes. An analysis of mutagen annotations on SPs appearing on all enzymes leads to the conclusion that mutations on SPs tend to damage the enzymatic function much more than expected from a background model, hence SPs are of high importance to enzymatic functions. SPs that lie in 3D pockets that are shared by active and binding sites, are shown to be significantly enriched by glycine, leading to the hypothesis that they are responsible for conformational plasticity. Finally we show that SPs can partially resolve outstanding difficult problems of convergent evolution by representing correctly enzyme functions in spite of remote homologies in sequence and in structure.

A Gibbs free energy correlation for automated docking of carbohydrates

Thermodynamic information can be inferred from static atomic configurations. To model the thermodynamics of carbohydrate binding to proteins accurately, a large binding data set has been assembled from the literature. The data set contains information from 262 unique protein-carbohydrate crystal structures for which experimental binding information is known. Hydrogen atoms were added to the structures and training conformations were generated with the automated docking program AutoDock 3.06, resulting in a training set of 225,920 all-atom conformations. In all, 288 formulations of the AutoDock 3.0 free energy model were trained against the data set, testing each of four alternate methods of computing the van der Waals, solvation, and hydrogen-bonding energetic components. The van der Waals parameters from AutoDock 1 produced the lowest errors, and an entropic model derived from statistical mechanics produced the only models with five physically and statistically significant coefficients. Eight models predict the Gibbs free energy of binding with an error of less than 40% of the error of any similar models previously published.

Hexose/Pentose and Hexitol/Pentitol Metabolism

Escherichia coli and Salmonella enterica serovar Typhimurium exhibit a remarkable versatility in the usage of different sugars as the sole source of carbon and energy, reflecting their ability to make use of the digested meals of mammalia and of the ample offerings in the wild. Degradation of sugars starts with their energy-dependent uptake through the cytoplasmic membrane and is carried on further by specific enzymes in the cytoplasm, destined finally for degradation in central metabolic pathways. As variant as the different sugars are, the biochemical strategies to act on them are few. They include phosphorylation, keto-enol isomerization, oxido/reductions, and aldol cleavage. The catabolic repertoire for using carbohydrate sources is largely the same in E. coli and in serovar Typhimurium. Nonetheless, significant differences are found, even among the strains and substrains of each species. We have grouped the sugars to be discussed according to their first step in metabolism, which is their active transport, and follow their path to glycolysis, catalyzed by the sugar-specific enzymes. We will first discuss the phosphotransferase system (PTS) sugars, then the sugars transported by ATP-binding cassette (ABC) transporters, followed by those that are taken up via proton motive force (PMF)-dependent transporters. We have focused on the catabolism and pathway regulation of hexose and pentose monosaccharides as well as the corresponding sugar alcohols but have also included disaccharides and simple glycosides while excluding polysaccharide catabolism, except for maltodextrins.

Characterization of potential stress responses in ancient Siberian permafrost psychroactive bacteria

Past studies of cold-acclimated bacteria have focused primarily on organisms not capable of sub-zero growth. Siberian permafrost isolates Exiguobacterium sp. 255-15 and Psychrobacter sp. 273-4, which grow at subzero temperatures, were used to study cold-acclimated physiology. Changes in membrane composition and exopolysaccharides were defined as a function of growth at 24, 4 and -2.5 degrees C in the presence and absence of 5% NaCl. As expected, there was a decrease in fatty acid saturation and chain length at the colder temperatures and a further decrease in the degree of saturation at higher osmolarity. A shift in carbon source utilization and antibiotic resistance occurred at 4 versus 24 degrees C growth, perhaps due to changes in the membrane transport. Some carbon substrates were used uniquely at 4 degrees C and, in general, increased antibiotic sensitivity was observed at 4 degrees C. All the permafrost strains tested were resistant to long-term freezing (1 year) and were not particularly unique in their UVC tolerance. Most of the tested isolates had moderate ice nucleation activity, and particularly interesting was the fact that the Gram-positive Exiguobacterium showed some soluble ice nucleation activity. In general the features measured suggest that the Siberian organisms have adapted to the conditions of long-term freezing at least for the temperatures of the Kolyma region which are -10 to -12 degrees C where intracellular water is likely not frozen.

A single nucleotide polymorphism alters the activity of the renal Na+:Cl- cotransporter and reveals a role for transmembrane segment 4 in chloride and thiazide affinity

The thiazide-sensitive Na+:Cl- cotransporter is the major salt transport pathway in the distal convoluted tubule of the kidney, and a role of this cotransporter in blood pressure homeostasis has been defined by physiological studies on pressure natriuresis and by its involvement in monogenic diseases that feature arterial hypotension or hypertension. Data base analysis revealed that 135 single nucleotide polymorphisms along the human SLC12A3 gene that encodes the Na+:Cl- cotransporter have been reported. Eight are located within the coding region, and one results in a single amino acid change; the residue glycine at the position 264 is changed to alanine (G264A). This residue is located within the fourth transmembrane domain of the predicted structure. Because Gly-264 is a highly conserved residue, we studied the functional properties of this polymorphism by using in vitro mutagenesis and the heterologous expression system in Xenopus laevis oocytes. G264A resulted in a significant and reproducible reduction ( approximately 50%) in (22)Na+ uptake when compared with the wild type cotransporter. The affinity for extracellular Cl- and for thiazide diuretics was increased in G264A. Western blot analysis showed similar immunoreactive bands between the wild type and the G264A cotransporters, and confocal images of oocytes injected with enhanced green fluorescent protein-tagged wild type and G264A cotransporter showed no differences in the protein surface expression level. These observations suggest that the G264A polymorphism is associated with reduction in the substrate translocation rate of the cotransporter, due to a decrease in the intrinsic activity. Our study also reveals a role of the transmembrane segment 4 in defining the affinity for extracellular Cl- and thiazide diuretics.

BHB: a simple knowledge-based scoring function to improve the efficiency of database screening

A new knowledge-based scoring function was developed in this work to facilitate the rapid ranking of ligands in databases. The acronym of the method is BHB based on the descriptors it utilizes: buriedness, hydrogen bonding, and binding energy. Receptor buriedness is a measure of how well molecules occupy the binding pocket in comparison to known high-affinity ligands or, alternatively, whether they have contact with identified residues in the pocket. The possibility of hydrogen bond formation is checked for selected residues that are recognized as being important in the binding of known ligands. The approximate binding energy is calculated from the thermodynamic cycle using the optimized bound and free solvent conformations of the ligand-receptor system. The information necessary for the scoring function can ideally be gleaned from the 3D structure of the receptor-ligand complex. Alternatively, the descriptors can be derived from the 3D structure of the unbound receptor, provided this receptor has a known ligand that binds to the given site with nanomolar activity. We show that the new scoring functions provide up to 12 times improvement in enrichment compared to the popular commercial docking program GOLD.

Agonist-induced isomerization in a glutamate receptor ligand-binding domain. A kinetic and mutagenetic analysis

Agonist binding to glutamate receptor ion channels occurs within an extracellular domain (S1S2) that retains ligand affinity when expressed separately. S1S2 is homologous to periplasmic binding proteins, and it has been proposed that a Venus flytrap-style cleft closure triggers opening of glutamate receptor ion channels. Here we compare the kinetics of S1S2-agonist binding to those of the periplasmic binding proteins and show that the reaction involves an initial rapid association, followed by slower conformational changes that stabilize the complex: "docking" followed by "locking." The motion detected here reflects the mechanism by which the energy of glutamate binding is converted into protein conformational changes within S1S2 alone. In the intact channel, these load-free conformational changes are harnessed and possibly modified as the agonist binding reaction is used to drive channel opening and subsequent desensitization. Using mutagenesis, key residues in each step were identified, and their roles were interpreted in light of a published S1S2 crystal structure. In contrast to the Venus flytrap proposal, which focuses on motion between the two lobes as the readout for agonist binding, we argue that smaller, localized conformational rearrangements allow agonists to bridge the cleft, consistent with published hydrodynamic measurements.

Specificity mutants of the binding protein of the oligopeptide transport system of Lactococcus lactis

The kinetic properties of wild-type and mutant oligopeptide binding proteins of Lactococcus lactis were determined. To observe the properties of the mutant proteins in vivo, the oppA gene was deleted from the chromosome of L. lactis to produce a strain that was totally defective in oligopeptide transport. Amplified expression of the oppA gene resulted in an 8- to 12-fold increase in OppA protein relative to the wild-type level. The amplified expression was paralleled by increased bradykinin binding activity, but had relatively little effect on the overall transport of bradykinin via Opp. Several site-directed mutants were constructed on the basis of a comparison of the primary sequences of OppA from Salmonella enterica serovar Typhimurium and L. lactis, taking into account the known structure of the serovar Typhimurium protein. Putative peptide binding-site residues were mutated. All the mutant OppA proteins exhibited a decreased binding affinity for the high-affinity peptide bradykinin. Except for OppA(D471R), the mutant OppA proteins displayed highly defective bradykinin uptake, whereas the transport of the low-affinity substrate KYGK was barely affected. Cells expressing OppA(D471R) had a similar K(m) for transport, whereas the V(max) was increased more than twofold as compared to the wild-type protein. The data are discussed in the light of a kinetic model and imply that the rate of transport is determined to a large extent by the donation of the peptide from the OppA protein to the translocator complex.

A molecular envelope of the ligand-binding domain of a glutamate receptor in the presence and absence of agonist

Solution scattering studies were performed on a ligand-binding domain (S1S2) of a glutamate receptor ion channel (GluR) in order to study GluR-binding and signal-transduction mechanisms. The core of the ligand-binding domain is homologous to prokaryotic periplasmic binding proteins (PBP), whose binding mechanism involves a dramatic cleft closure: the "Venus flytrap". Several models of GluR function have proposed that a similar cleft closure is induced by agonist binding. We have directly tested this putative functional homology by measuring the radius of gyration of S1S2 in the presence and absence of saturating concentrations of agonists. In contrast to the PBP, S1S2 shows no reduction in radius of gyration upon agonist binding, excluding a comparably large conformational change. Furthermore, we determined an ab initio molecular envelope for our S1S2 construct, which also contains the peptides that connect the PBP homology core to the three transmembrane domains and to an N-terminal domain. By fitting an atomic model of the ligand-binding domain core to the envelope of our extended construct, we were able to establish the likely position of these connecting peptides. Their positions relative to one another and to the expected sites of an agonist-induced conformational change suggest that ion channel gating and desensitization may involve more subtle and complex mechanisms than have been assumed based on the structural homology to the PBP.

Molecular dynamics simulations of protein-tyrosine phosphatase 1B. I. ligand-induced changes in the protein motions

Activity of enzymes, such as protein tyrosine phosphatases (PTPs), is often associated with structural changes in the enzyme, resulting in selective and stereospecific reactions with the substrate. To investigate the effect of a substrate on the motions occurring in PTPs, we have performed molecular dynamics simulations of PTP1B and PTP1B complexed with a high-affinity peptide DADEpYL, where pY stands for phosphorylated tyrosine. The peptide sequence is derived from the epidermal growth factor receptor (EGFR988-993). Simulations were performed in water for 1 ns, and the concerted motions in the protein were analyzed using the essential dynamics technique. Our results indicate that the predominately internal motions in PTP1B occur in a subspace of only a few degrees of freedom. Upon substrate binding, the flexibility of the protein is reduced by approximately 10%. The largest effect is found in the protein region, where the N-terminal of the substrate is located, and in the loop region Val198-Gly209. Displacements in the latter loop are associated with the motions in the WPD loop, which contains a catalytically important aspartic acid. Estimation of the pKa of the active-site cysteine along the trajectory indicates that structural inhomogeneity causes the pKa to vary by approximately +/-1 pKa unit. In agreement with experimental observations, the active-site cysteine is negatively charged at physiological pH.

mRNA cap recognition: dominant role of enhanced stacking interactions between methylated bases and protein aromatic side chains

We have determined, by high resolution x-ray analysis, 10 structures comprising the mRNA cap-specific methyltransferase VP39 or specific mutants thereof in the presence of methylated nucleobase analogs (N1-methyladenine, N3-methyladenine, N1-methylcytosine, N3-methylcytosine) and their unmethylated counterparts, or nucleoside N7-methylguanosine. Together with solution affinity studies and previous crystallographic data for N7-methylguanosine and its phosphorylated derivatives, these data demonstrate that only methylated, positively charged bases are bound, indicating that their enhanced stacking with two aromatic side chains of VP39 (Tyr 22 and Phe 180) plays a dominant role in cap recognition. Four key features characterize this stacking interaction: (i) near perfect parallel alignment between the sandwiched methylated bases and aromatic side chains, (ii) substantial areas of overlap in the two-stacked rings, (iii) a 3.4-A interplanar spacing within the overlapping region, and (iv) positive charge in the heterocyclic nucleobase.

Analysis of the dynamics of rhizomucor miehei lipase at different temperatures

The dynamics of Rhizomucor miehei lipase has been studied by molecular dynamics simulations at temperatures ranging from 200-500K. Simulations carried out in periodic boundary conditions and using explicit water molecules were performed for 400 ps at each temperature. Our results indicate that conformational changes and internal motions in the protein are significantly influenced by the temperature increase. With increasing temperature, the number of internal hydrogen bonds decreases, while surface accessibility, radius of gyration and the number of residues in random coil conformation increase. In the temperature range studied, the motions can be described in a low dimensional subspace, whose dimensionality decreases with increasing temperature. Approximately 80% of the total motion is described by the first (i) 80 eigenvectors at T=200K, (ii) 30 eigenvectors at T=300K and (iii) 10 eigenvectors at T=400K. At high temperature, the alpha-helix covering the active site in the native Rhizomucor miehei lipase, the helix at which end the active site is located, and in particular, the loop (Gly35-Lys50) show extensive flexibility.

Dominant role of local dipolar interactions in phosphate binding to a receptor cleft with an electronegative charge surface: equilibrium, kinetic, and crystallographic studies

Stringent specificity and complementarity between the receptor, a periplasmic phosphate-binding protein (PBP) with a two-domain structure, and the completely buried and dehydrated phosphate are achieved by hydrogen bonding or dipolar interactions. We recently found that the surface charge potential of the cleft between the two domains that contains the anion binding site is intensely electronegative. This novel finding prompted the study reported here of the effect of ionic strength on the equilibrium and rapid kinetics of phosphate binding. To facilitate this study, Ala197, located on the edge of the cleft, was replaced by a Trp residue (A197W PBP) to generate a fluorescence reporter group. The A197W PBP-phosphate complex retains wild-type Kd and X-ray structure beyond the replacement residue. The Kd (0.18 microM) at no salt is increased by 20-fold at greater than 0.30 M NaCl. Stopped-flow fluorescence kinetic studies indicate a two-step binding process: (1) The phosphate (L) binds, at near diffusion-controlled rate, to the open cleft form (Po) of PBP to produce an intermediate, PoL. This rate decreases with increasing ionic strength. (2) The intermediate isomerizes to the closed-conformation form, PcL. The results indicate that the high specificity, affinity, and rate of phosphate binding are not influenced by the noncomplementary electronegative surface potential of the cleft. That binding depends almost entirely on local dipolar interactions with the receptor has important ramification in electrostatic interactions in protein structures and in ligand recognition.

Structure-based redesign of corepressor specificity of the Escherichia coli purine repressor by substitution of residue 190

Guanine or hypoxanthine, physiological corepressors of the Escherichia coli purine repressor (PurR), promote formation of the ternary PurR-corepressor-operator DNA complex that functions to repress pur operon gene expression. Structure-based predictions on the importance of Arg190 in determining 6-oxopurine specificity and corepressor binding affinity were tested by mutagenesis, analysis of in vivo function, and in vitro corepressor binding measurements. Replacements of Arg190 with Ala or Gln resulted in functional repressors in which binding of guanine and hypoxanthine was retained but specificity was relaxed to permit binding of adenine. X-ray structures were determined for ternary complexes of mutant repressors with purines (adenine, guanine, hypoxanthine, and 6-methylpurine) and operator DNA. These structures indicate that R190A binds guanine, hypoxanthine, and adenine with nearly equal, albeit reduced, affinity in large part because of a newly made compensatory hydrogen bond between the rotated hydroxyl side chain of Ser124 and the exocyclic 6 positions of the purines. Through direct and water-mediated contacts, the R190Q protein binds adenine with a nearly 75-fold higher affinity than the wild type repressor while maintaining wild type affinity for guanine and hypoxanthine. The results establish at the atomic level the basis for the critical role of Arg190 in the recognition of the exocyclic 6 position of its purine corepressors and the successful redesign of corepressor specificity.

The role of the variable region in Tet repressor for inducibility by tetracycline

A set of deletions and substitutions to alanine was introduced into the loop separating helices alpha8 and alpha9 of Tn10 Tet repressor (TetR). This region appears as an unstructured loop in the crystal structure of the TetR(D).([Mg-tc]+)2 complex and is the only internal segment of variable length in an alignment of Tet repressors from seven different resistance determinants. In vivo analysis of 10 mutants shows that this loop is important for inducibility by tetracycline (tc), whereas DNA binding is not or only marginally affected. All deletions have an induction-deficient TetRS phenotype, but the corresponding substitutions do not or only slightly affect inducibility. The purified mutant TetR proteins have a reduced affinity for tc in vitro that correlates with their lack of inducibility. The association rate of [Mg-tc]+ to the TetR mutants is enhanced. Since none of the mutated residues contacts tc directly in the crystal structure, we propose that the length of the loop is important for the structural transition between a closed, tc binding and an open, operator binding conformation of TetR. We propose that the deletions in the loop shift the equilibrium between both forms toward the open, operator binding conformation.

Crystal structure of PotD, the primary receptor of the polyamine transport system in Escherichia coli

PotD protein is a periplasmic binding protein and the primary receptor of the polyamine transport system, which regulates the polyamine content in Escherichia coli. The crystal structure of PotD in complex with spermidine has been solved at 2.5-A resolution. The PotD protein consists of two domains with an alternating beta-alpha-beta topology. The polyamine binding site is in a central cleft lying in the interface between the domains. In the cleft, four acidic residues recognize the three positively charged nitrogen atoms of spermidine, while five aromatic side chains anchor the methylene backbone by van der Waals interactions. The overall fold of PotD is similar to that of other periplasmic binding proteins, and in particular to the maltodextrin-binding protein from E. coli, despite the fact that sequence identity is as low as 20%. The comparison of the PotD structure with the two maltodextrin-binding protein structures, determined in the presence and absence of the substrate, suggests that spermidine binding rearranges the relative orientation of the PotD domains to create a more compact structure.

Reactivity and structure of a mouse anti-F(ab')2 IgM. Comparison of its variable region sequences with those of a structurally close polyreactive natural IgM

IE12 is a monoclonal IgM with strong anti-F(ab')2 activity that inhibits the binding of normal mouse IgG to self antigens. In this study, we found that this IgM was also reactive with several monoclonal antibodies (mAbs), myeloma proteins and B lymphocytes from normal BALB/c mouse. The nucleotide sequences of the variable region of the heavy and light chains of IE12 were determined, and compared to those of another mAb already described in the literature. This mAb uses the same light chain and also the same VH, D and JH segments, but unlike IE12, is polyreactive. The comparison of the amino acid composition of these two mAbs and of the computer predictions for their structure and hydrophilicity indicated that the most striking difference between them was located in the third complementarity determining region (CDR3) of the heavy chain. Indeed, they used the same D segment but translated in two different reading frames, leading to different amino acid compositions. The CDR3 of IE12 contains aliphatic amino acids, while that of the polyreactive IgM does not. In addition, IE12 has two prolines, one at each at each extremity of its D segment, that could confer a certain rigidity to this region. Finally, the CDR3 of IE12 is predicted to be hydrophobic, while the one of the polyreactive IgM is predicted to be hydrophilic and more flexible, suggesting that the hydrophilicity and the flexibility of this region might be critical for polyreactivity.

Critical amino acids responsible for converting specificities of proteins and for enhancing enzyme evolution are located around beta-turn potentials: data-based prediction

Various reports have described that amino acid substitutions can alter substrate, positional, inhibitory, and target gene specificities of proteins. By using the method of Chou and Fasman, the present work predicts that critical amino acids for converting these substrate specificities of trypsin, L-lactate dehydrogenase, aspartate aminotransferase, beta-lactamase, and cytochrome P-450 are found to exist within regions predicted as beta-turns. The ratios of hydroxylation and oxygenation positions of substrates by cytochrome P-450 and lipoxygenase, respectively, are varied by changes of the protein structures, probably around turn conformations. Inhibitory specificities of bovine pancreatic trypsin inhibitor and alpha 1-antitrypsin and target gene specificity of glucocorticoid receptor are converted by changing turn structures. Occurrence of beta-turn probabilities can be predicted around the amino acid alteration positions of an evolutionally antecedent protein of a nylon degradation enzyme. These findings will have relevance to work on protein engineering and enzyme evolution.

Inversion of receptor binding preferences by mutagenesis: free energy thermodynamic integration studies on sugar binding to L-arabinose binding proteins

The Escherichia coli L-arabinose-binding protein (ABP) participates as a specific receptor in the transport of L-arabinose, D-fucose, and D-galactose through the periplasmic space. The wild-type protein binds L-arabinose about 40 times more strongly than D-fucose. A mutation of the protein at position 108 (Met-->Leu) causes a specificity change. The Met108Leu ABP slightly prefers binding of D-fucose over L-arabinose. Molecular dynamics (MD) and thermodynamic integration (TI) computer simulations were performed to study the mechanism of sugar discrimination and specificity change based on the known high-resolution X-ray structures. The specificity change was evaluated by calculating the difference in free energy of L-arabinose versus D-fucose bound to wild-type and Met108Leu ABP. The calculated free energy differences are consistent with the experimentally observed specificity of wild-type and Met108Leu ABP. The simulations indicate that the specificity change of Met108Leu ABP is accomplished mainly by reduced Lennard-Jones interactions of residue 108 with L-arabinose and improved interactions with D-fucose. In addition to MD/TI calculations on sugar binding, finite difference Poisson-Boltzmann calculations were performed to identify the most stable ionization state of buried ionizable residues in ABP.

Crystallographic evidence of a large ligand-induced hinge-twist motion between the two domains of the maltodextrin binding protein involved in active transport and chemotaxis

The periplasmic maltodextrin binding protein of Escherichia coli serves as an initial receptor for the active transport of and chemotaxis toward maltooligosaccharides. The three-dimensional structure of the binding protein complexed with maltose has been previously reported [Spurlino, J. C., Lu, G.-Y., & Quiocho, F. A. (1991) J. Biol. Chem. 266, 5202-5219]. Here we report the structure of the unliganded form of the binding protein refined to 1.8-A resolution. This structure, combined with that for the liganded form, provides the first crystallographic evidence that a major ligand-induced conformational change occurs in a periplasmic binding protein. The unliganded structure shows a rigid-body "hinge-bending" between the two globular domains by approximately 35 degrees, relative to the maltose-bound structure, opening the sugar binding site groove located between the two domains. In addition, there is an 8 degrees twist of one domain relative to the other domain. The conformational changes observed between this structure and the maltose-bound structure are consistent with current models of maltose/maltodextrin transport and maltose chemotaxis and solidify a mechanism for receptor differentiation between the ligand-free and ligand-bound forms in signal transduction.

Protein-ligand energetics assessed using deoxy and fluorodeoxy sugars in equilibrium binding and high resolution crystallographic studies

Hydrogen bonding interactions are one of the most important single factors in protein-ligand interactions and molecular recognition. To probe the energetics of the interactions, we have analyzed the binding of 1-deoxy-, 2-deoxy- and 6-fluoro-6-deoxy- analogues of D-galactose (Gal) to a primary high-affinity periplasmic receptor for monosaccharide active transport. Kd values and atomic structures refined at 1.81 to 1.45 A resolution of the complexes have been determined and compared with those of Gal binding. With binding site residues and the bound modified sugars in nearly identical positions as found in the complex with Gal, the binding of 1-deoxy-Gal or 2-deoxy-Gal reflects the overall contribution of 1.8 kcal mol-1 per hydrogen bond (neutral-charge type) to the affinity of Gal. Neglected in these estimates is the contribution of van der Waals' forces that accompany the formation of hydrogen bonds with each sugar hydroxyl. Contrary to expectations, the 6-fluoro-6-deoxy analogue proved to be an inadequate probe of Gal OH6 as a hydrogen bond donor due to the binding of a new water molecule and structural changes arising from the electronegative fluoro group. This study sheds new light on the energetics of protein-ligand interactions and the use of engineered ligands in assessing these interactions.

Interdomain salt bridges modulate ligand-induced domain motion of the sulfate receptor protein for active transport

The refined crystal structure of the liganded form of the Salmonella typhimurium sulfate-binding protein, a periplasmic receptor of active transport, is made up of two globular domains bisected by a deep cleft wherein the dehydrated sulfate is completely engulfed and bound by hydrogen bonds and van der Waals' forces. Two salt bridges (between Glu15 and Arg174 and between Asp68 and Arg134) span the cleft opening. To elucidate the role of the inter-domain salt bridges in the ligand-induced domain motion, the acidic residues were changed (singly and together) to their corresponding amide side-chains by site-directed mutagenesis of the recombinant Escherichia coli sulfate-binding protein. Rapid kinetics and equilibrium measurements of sulfate binding to the purified mutant proteins demonstrate that these salt bridges stabilize the closed liganded form of the receptor and modulate the rate of cleft opening. Our results have new implications in understanding the dynamics of many other multidomain proteins that undergo similar large-scale domain motions.