Mutational analysis of N381, a key trimer contact residue in Tsr, the Escherichia coli serine chemoreceptor

Chemoreceptors such as Tsr, the serine receptor, function in trimer-of-dimer associations to mediate chemotactic behavior in Escherichia coli. The two subunits of each receptor homodimer occupy different positions in the trimer, one at its central axis and the other at the trimer periphery. Residue N381 of Tsr contributes to trimer stability through interactions with its counterparts in a central cavity surrounded by hydrophobic residues at the trimer axis. To assess the functional role of N381, we created and characterized a full set of amino acid replacements at this Tsr residue. We found that every amino acid replacement at N381 destroyed Tsr function, and all but one (N381G) of the mutant receptors also blocked signaling by Tar, the aspartate chemoreceptor. Tar jamming reflects the formation of signaling-defective mixed trimers of dimers, and in vivo assays with a trifunctional cross-linking reagent demonstrated trimer-based interactions between Tar and Tsr-N381 mutants. Mutant Tsr molecules with a charged amino acid or proline replacement exhibited the most severe trimer formation defects. These trimer-defective receptors, as well as most of the trimer-competent mutant receptors, were unable to form ternary signaling complexes with the CheA kinase and with CheW, which couples CheA to receptor control. Some of the trimer-competent mutant receptors, particularly those with a hydrophobic amino acid replacement, may not bind CheW/CheA because they form conformationally frozen or distorted trimers. These findings indicate that trimer dynamics probably are important for ternary complex assembly and that N381 may not be a direct binding determinant for CheW/CheA at the trimer periphery.     

A model for switch-like phenomena in biological systems.

We present a model for the activity of protein clusters based on a simultaneous desorption of an activator (agonist, substrate molecule, etc.) and an inactivator (antagonist, inhibitor, etc.) caused by the collision or interaction between two effector molecules (e.g. receptors, enzymes). This model gives rise to switch-like dose-response curves, which are difficult to explain by ordinary co-operativity. It fits with recent experimental results obtained on single cells. Some other interesting aspects of the model are also pointed out. The model is similar to the model used to explain steep 'dose-response curves' in heterogeneous catalysis, caused by the reaction between two different molecules or atoms on the surface of the catalyst.     

Aggregation of cyclodextrins as an important factor to determine their complexation behavior

Here, we report a study on the complexation behavior of carotenoids with cyclodextrins (CDs) using solubility experiments and molecular-modelling methods. Carotenoids are an important group of naturally occurring dyes found in vegetables and fruits. Their antioxidant property has initiated investigations on their possible use as drugs. However, carotenoids are lipophilic molecules with very little inherent aqueous solubility. Cyclodextrin complexation has been widely used in order to increase the potential applications of hydrophobic compounds. Thus, the aim of our investigation was to design carotenoids with enhanced water solubility by cyclodextrin complexation. Molecular modelling of carotenoid-cyclodextrin complexes with a 1 : 1 stoichiometry successfully explained the experimentally observed capability of beta-cyclodextrins (beta-CDs) to form complexes with carotenoids as opposed to alpha-cyclodextrins (alpha-CDs) and gamma-cyclodextrins (gamma-CDs). Furthermore, molecular-dynamics calculations revealed that the aggregation properties of CD derivatives significantly influence their complexation behavior. Our docking calculations showed that RAMEB (random methylated beta-CD) is the beta-CD derivative that possesses the lowest tendency to aggregate. Solubility experiments yielded the same results, namely, RAMEB complexes possess the best water solubility. Our results showed that complexation of a ligand not buried inside of the CD cavity is dependent on two factors: i) the geometry of the inclusion part of the complex; ii) the self-aggregation property of the CD itself. The lower affinity the CDs possess for self-aggregation, the more likely are they involved in interactions with carotenoids. These results suggest that self-aggregation of CDs should be considered as an important parameter determining complexation in general.     

Crystal structure of beta-cyclodextrin-benzoic acid inclusion complex

The inclusion complex of beta-cyclodextrin (beta-CD) with benzoic acid (BA) has been characterized crystallographically. Two beta-CDs cocrystallize with two BAs, 0.7 ethanol and 20.65 water molecules [2(C(6)H(10)O(5))(7).2(C(7)H(6)O(2)).0.7(C(2)H(6)O).20.65H(2)O] in the triclinic space group P1 with unit cell constants: a=15.210(1), b=15.678(1), c=15.687(1) A, alpha=89.13(1), beta=74.64(1), gamma=76.40(1) degrees. The anisotropic refinement of 1840 atomic parameters against 16,201 X-ray diffraction data converged at R=0.078. In the crystal lattice, beta-CD forms dimers stabilized by direct O-2(m)_1/O-3(m)_1...O-2(n)_2/O-3(n)_2 hydrogen bonds (intradimer) and by indirect O-6(m)_1...,O-6(n)_2 hydrogen bonds with one or two bridging water molecules joined in between (interdimer). These dimers are stacked like coins in a roll constructing endless channels where the guest molecules are included. The BA molecules protrude with their COOH groups at the beta-CD O-6-sides and are maintained in positions by hydrogen bonding to the surrounding O-6-H groups and water molecules. Water molecules (20.65) are distributed over 30 positions in the interstices between beta-CD molecules, except the water sites W-1, W-2 that are located in the channel of the beta-CD dimer. Water site W-2 is hydrogen bonded to the disordered ethanol molecule (occupancy 0.7).     

Monte Carlo simulation of 4‐α‐glucanotransferase reaction

4‐α‐Glucanotransferase (GTase, D ‐enzyme) catalyzes disproportionation between two short polymers of maltooligosaccharides linked by α‐1,4‐glucoside bonds. Using action modes of the potato GTase for the donor and acceptor substrates, the Monte Carlo method was applied to simulate the GTase reaction. The simulation starts from a single enzyme molecule and a finite number (10⁵) of substrate molecules. All selection processes were performed using random numbers produced by computer. The initial substrates were from trimer to 10‐mer. In every case, the final stage was the steady‐state distribution of polymers. The steady‐state distribution by the potato GTase reaction was different from those by the hypothetical random disproportionation reaction. The simulated data from the reaction of potato GTase and trimer almost quantitatively agreed with experimental data. The mechanism of the GTase reaction was accumulation of probabilistic processes and was well simulated by the Monte Carlo method. GTase randomizes the overall distribution of chain length of the substrate. Therefore the GTase reaction is an entropy‐driven process. © 1999 John Wiley & Sons, Inc. Biopoly 50: 145–151, 1999     

Structure of the complex of beta-cyclodextrin with beta-naphthyloxyacetic acid in the solid state and in aqueous solution

The structure of the complex of beta-cyclodextrin (cyclomaltoheptaose) with beta-naphthyloxyacetic acid was studied in solid state by X-ray diffraction and in aqueous solution by 1H NMR spectroscopy. The complex crystallizes in the channel mode, space group C2, with a stoichiometry of 2:1; two beta-cyclodextrin molecules related by a twofold crystal axis form dimers, in the cavity of which one guest molecule is found on average. The above stoichiometry indicates one guest per beta-CD dimer statistically oriented over two positions or two guest molecules in pi-pi interactions in half of the beta-CD dimers and the rest of the beta-CD dimers empty. In addition, occupancy of 0.5 for the guest per every beta-CD dimer is in accord with the occupancy of the two disordered primary hydroxyls. These two hydroxyl groups, to which the carboxylic oxygen atoms of the guest are hydrogen bonded, point towards the interior of the beta-CD cavity. In aqueous solution, the 1H NMR spectroscopic study indicated that there is a mixture of complexes with host-guest stoichiometries both 1:1 and 2:1.     

Zinc nanoparticles interact with olfactory receptor neurons

Zinc metal nanoparticles strongly enhance odorant responses of olfactory receptor neurons. Olfactory receptors belong to the large superfamily of G-protein coupled receptors. A theoretical model based on experimental results explains a stoichiometry of metal nanoparticles receptor interaction. The model is similar to that used by A.V. Hill for the binding reaction between hemoglobin and oxygen. The model predicted that one metal nanoparticle binds two receptor molecules to create a dimer. This result is consistent with the evidence that many G-protein-coupled receptors form dimers or larger oligomers.     

肿瘤坏死因子与其受体相互作用的计算机模拟研究

利用同源模建方法,以肿瘤坏死因子(TNF)R55受体胞外区的晶体结构为参考模板,预测了TNFR75受体胞外区Cys18-Phe147片段的三维结构。根据R55受体胞外区与LT相结合的复合物的晶体结构,预测了TNF与R55及R75胞外区的复合物的三维结构,模拟了TNF与受体之间的相互作用。由于TNF与受体的作用形式是三聚体对三聚体,因此在模拟TNF与受体相互作用时选择了包括一个非对称的TNF三聚体和一个受体(R55或R75)单体的模拟系统。结合已有的突变体实验结果,利用计算机分析手段,发现了一些TNF突变体之所以具有受体选择性的三维结构基础和发挥了关键作用的氨基酸残基以及这些残基之间的主要作用形式。研究深化了对已有的突变体实验结果的认识,建立了不同的实验结果之间的内在关联,为以后有目的的新型突变体设计和实验研究打下了基础。     

Protein components of the nonactivated glucocorticoid receptor.

The nonactivated glucocorticoid receptor (Mr approximately 350,000) of WEHI-7 mouse lymphoma cells was investigated with respect to the stoichiometry of protein subunits. Cross-linking patterns obtained by affinity labeling and denaturing gel electrophoresis revealed a heterotetramer consisting of one receptor polypeptide in association with two 90- and one approximately 50-kDa subunits. The receptor stabilized by molybdate, disulfide bond formation, or chemical cross-linking was purified roughly 6000-fold by immunoaffinity chromatography and analyzed by gel electrophoresis and immunoblotting. The 90-kDa component was consistently detected in a 2:1 ratio with respect to the receptor polypeptide and was identified as the 90-kDa heat shock protein, hsp90. A 70-kDa heat shock protein was found in both stabilized and nonstabilized receptors and bound to the immunomatrix independent of receptor. The additional receptor subunit was unequivocally identified as the 59-kDa protein previously described (Tai, P.-K. K., Maeda, Y., Nakao, K., Wakim, N. G., Duhring, J. L., and Faber, L. E. (1986) Biochemistry 25, 5269-5275). This component was found only in complexes cross-linked via amino groups. It was removed from the molybdate-stabilized receptor under our purification conditions, thus leaving behind a trimer composed of the receptor polypeptide and two molecules of hsp90. In the absence of hormone, the receptor had the same subunit composition as in its presence.     

Crystal structures of the signal transducing protein GlnK from Thermus thermophilus HB8

The Thermus thermophilus HB8 genome encodes a signal transducing PII protein, GlnK. The crystal structures of GlnK have been determined in two different space groups, P2(1)2(1)2(1) and P3(1)21. The PII protein has the T-loop, which is essential for interactions with receptor proteins. In both crystal forms, three GlnK molecules form a trimer in the asymmetric unit. In one P2(1)2(1)2(1) crystal form, the three T-loops in the trimer are disordered, while in another P2(1)2(1)2(1) crystal form, the T-loop from one molecule in the trimer is ordered. In the P3(1)21 crystal, one T-loop is ordered while the other two T-loops are disordered. The conformations of the ordered T-loops significantly differ between the two crystal forms; one makes the alpha-helix in the middle of the T-loop, while the other has an extension of the beta-hairpin. Two different conformations are captured by the crystal contacts. The observation of multiple T-loop conformations suggests that the T-loop could potentially exhibit "polysterism," which would be important for interactions with receptor proteins. The crystal structures of the nucleotide-bound forms, GlnK.ATP and GlnK.ADP, have also been determined. ATP/ADP binding within a cleft at the interface of two adjacent T. thermophilus GlnK monomers might affect the conformation of the T-loop.     

Crystalline beta-cyclodextrin.12H2O reversibly dehydrates to beta-cyclodextrin.10.5 H2O under ambient conditions.

In contact with mother liquor, crystalline beta-cyclodextrin (beta-CD) hydrate has composition approximately beta-CD.12H2O. If crystals are dried at ambient conditions (18 degrees C, approximately 50% humidity), the unit cell volume diminishes approximately 30 to 50 A3. X-ray structure analysis of a dry crystal (0.89 A resolution, 4617 data, R = 0.059) showed the composition beta-CD.10.5 H2O, with approximately 5.5 water molecules in the beta-CD cavity (7 partially and 2 fully occupied sites) and approximately 5.0 between the beta-CD molecules. The positions of the beta-CD host and of most of the hydration waters are conserved during dehydration, but the occupancies of the waters in the beta-CD cavity diminish. Dry crystals put into solvent re-hydrate to the original form. The mechanism of de- and re-hydration is not evident.     

Different signaling roles of two conserved residues in the cytoplasmic hairpin tip of Tsr, the Escherichia coli serine chemoreceptor

Bacterial chemoreceptors form ternary signaling complexes with the histidine kinase CheA through the coupling protein CheW. Receptor complexes in turn cluster into cellular arrays that produce highly sensitive responses to chemical stimuli. In Escherichia coli, receptors of different types form mixed trimer-of-dimers signaling teams through the tips of their highly conserved cytoplasmic domains. To explore the possibility that the hairpin loop at the tip of the trimer contact region might promote interactions with CheA or CheW, we constructed and characterized mutant receptors with amino acid replacements at the two nearly invariant hairpin charged residues of Tsr: R388, the most tip-proximal trimer contact residue, and E391, the apex residue of the hairpin turn. Mutant receptors were subjected to in vivo tests for the assembly and function of trimers, ternary complexes, and clusters. All R388 replacements impaired or destroyed Tsr function, apparently through changes in trimer stability or geometry. Large-residue replacements locked R388 mutant ternary complexes in the kinase-off (F, H) or kinase-on (W, Y) signaling state, suggesting that R388 contributes to signaling-related conformational changes in the trimer. In contrast, most E391 mutants retained function and all formed ternary signaling complexes efficiently. Hydrophobic replacements of any size (G, A, P, V, I, L, F, W) caused a novel phenotype in which the mutant receptors produced rapid switching between kinase-on and -off states, indicating that hairpin tip flexibility plays an important role in signal state transitions. These findings demonstrate that the receptor determinants for CheA and CheW binding probably lie outside the hairpin tip of the receptor signaling domain.     

Thiol-directed immobilization of recombinant IgG-binding receptors

A genetic engineering approach is described for directed immobilization of IgG-binding receptors to a thiol-containing matrix using a single cysteine residue. The cysteine residue is introduced into the C-terminal part of receptors based on staphylococcal protein A. Receptors containing one, two or five IgG-binding domains were produced in Escherichia coli and subsequently immobilized on thiopropyl-Sepharose. A high amount, 5 mumol/ml gel (45 mg/ml), of recombinant receptor could be rapidly immobilized to the solid support and both the gel and the immobilized receptor could be regenerated by reduction and oxidation reactions. The gel was used for affinity purification of human IgG and analysis of IgG-binding capacity at different amounts of immobilized recombinant protein show the same maximal IgG-binding capacity (20-25 mg/ml) for all three immobilized receptors. However, at low substitution grade of receptors, the immobilized receptor molecules were shown to bind one (Z-Cys) and two (ZZ-Cys) IgG molecules, respectively. These results demonstrate that the immobilized protein molecules are in a functionally active form and that a two-domain receptor can bind two molecules of IgG without steric hindrance. Interestingly, the five-domain receptor (ZV-Cys), with a structure similar to native protein A, can only bind approximately two IgG molecules, despite the five-domain structure of the molecule.     

Addition of oligonucleotides to the 5''-terminus of DNA by T4 RNA ligase.

Bacteriophage T4-induced RNA ligase catalyzes the controlled template-independent addition of RNA to the 5'-phosphoryl end of large DNA molecules. Restriction enzyme-generated fragments of Co1E1 DNA with completely basepaired or cohesive ends and linear single-stranded ?X174 viral DNA were all good substrates. DNA molecules from 10 to 6000 nucleotides long were quantitatively joined in an hour to a number of different RNA homopolymers. The most efficient of these was A(pA)5; I(pI)5 and C(pC)5 were also utilized while U(pU)5 was not. The optimum ribohomopolymer length was six nucleotides. Joining of ribohomopolymers between 10 and 20 nucleotides long occurred at approximately 1/2 the maximal rate and a trimer was the shortest substrate. Thus RNA ligase provides a method for generating extensions of predetermined length and base composition at the 5'-end of large DNA molecules that complements the available procedures for extending the 3'-hydroxyl terminus of DNA.     

The insulin receptor tyrosine kinase domain in a chimaeric epidermal growth factor-insulin receptor generates Ca2+ signals through the PLC-gamma1 pathway.

The receptors for insulin (IR) and epidermal growth factor (EGFR) are members of the tyrosine kinase receptor (TKR) family. Despite homology of their cytosolic TK domains, both receptors induce different cellular responses. Tyrosine phosphorylation of insulin receptor substrate (IRS) molecules is a specific IR post-receptor response. The EGFR specifically activates phospholipase C-gamma1 (PLC-gamma1). Recruitment of substrate molecules with Src homology 2 (SH2) domains or phosphotyrosine binding (PTB) domains to phosphotyrosines in the receptor is one of the factors creating substrate specificity. In addition, it has been shown that the TK domains of the IR and EGFR show preferences to phosphorylate distinct peptides in vitro, suggesting additional mechanisms of substrate recognition. We have examined to what extent the substrate preference of the TK domain contributes to the specificity of the receptor in vivo. For this purpose we determined whether the IR TK domain, in situ, is able to tyrosine-phosphorylate substrates normally used by the EGFR. A chimaeric receptor, consisting of an EGFR in which the juxtamembrane and tyrosine kinase domains were exchanged by their IR counterparts, was expressed in CHO-09 cells lacking endogenous EGFR. This receptor was found to activate PLC-gamma1, indicating that the IR TK domain, in situ, is able to tyrosine phosphorylate substrates normally used by the EGFR. These findings suggest that the IR TK domain, in situ, has a low specificity for selection and phosphorylation of non-cognate substrates.     

Coarse-Grained Molecular Simulation of Epidermal Growth Factor Receptor Protein Tyrosine Kinase Multi-Site Self-Phosphorylation

Upon the ligand-dependent dimerization of the epidermal growth factor receptor (EGFR), the intrinsic protein tyrosine kinase (PTK) activity of one receptor monomer is activated, and the dimeric receptor undergoes self-phosphorylation at any of eight candidate phosphorylation sites (P-sites) in either of the two C-terminal (CT) domains. While the structures of the extracellular ligand binding and intracellular PTK domains are known, that of the ∼225-amino acid CT domain is not, presumably because it is disordered. Receptor phosphorylation on CT domain P-sites is critical in signaling because of the binding of specific signaling effector molecules to individual phosphorylated P-sites. To investigate how the combination of conventional substrate recognition and the unique topological factors involved in the CT domain self-phosphorylation reaction lead to selectivity in P-site phosphorylation, we performed coarse-grained molecular simulations of the P-site/catalytic site binding reactions that precede EGFR self-phosphorylation events. Our results indicate that self-phosphorylation of the dimeric EGFR, although generally believed to occur in trans, may well occur with a similar efficiency in cis, with the P-sites of both receptor monomers being phosphorylated to a similar extent. An exception was the case of the most kinase-proximal P-site-992, the catalytic site binding of which occurred exclusively in cis via an intramolecular reaction. We discovered that the in cis interaction of P-site-992 with the catalytic site was facilitated by a cleft between the N-terminal and C-terminal lobes of the PTK domain that allows the short CT domain sequence tethering P-site-992 to the PTK core to reach the catalytic site. Our work provides several new mechanistic insights into the EGFR self-phosphorylation reaction, and demonstrates the potential of coarse-grained molecular simulation approaches for investigating the complexities of self-phosphorylation in molecules such as EGFR (HER/ErbB) family receptors and growth factor receptor PTKs in general.     

Towards a general model for protein-substrate stereoselectivity

Protein-substrate interactions in enzymatic, neurological, and immunological systems are typically characterized by a high degree of stereoselectivity towards complex substrates. We propose a novel stereocenter-recognition (SR) model for stereoselectivity of proteins (or receptors in general) towards substrates that have multiple stereocenters, based on the topology of substrate stereocenters. The model provides the minimum number of substrate locations that need to enter into binding, nonbinding, or repulsive interactions with receptor sites, for stereoselectivity to occur. According to this model, a substrate location may interact with multiple receptor sites, or multiple substrate locations may interact with a single receptor site, but a stereoselective receptor has to offer, in the correct geometry, at least as many interactions as the required minimum number of substrate locations. The SR model predicts that stereoselectivity towards an acyclic substrate with N stereocenters distributed along a single chain requires interactions involving a minimum of N + 2 substrate locations, distributed over all stereocenters in the substrate, such that effectively three locations exist per stereocenter. Thus, enantioselective recognition of molecules with one chiral center requires a protein to interact with a minimum of three substrate locations, while stereoselectivity towards substrates with two or three stereocenters requires interactions with a minimum of four or five substrate locations, respectively, and so on. We demonstrate the general applicability of this model to protein-substrate interactions by interpreting several previous experimental observations.     

Analysis of the subunit stoichiometries in viral entry

Virions of the Human Immunodeficiency Virus (HIV) infect cells by first attaching with their surface spikes to the CD4 receptor on target cells. This leads to conformational changes in the viral spikes, enabling the virus to engage a coreceptor, commonly CCR5 or CXCR4, and consecutively to insert the fusion peptide into the cellular membrane. Finally, the viral and the cellular membranes fuse. The HIV spike is a trimer consisting of three identical heterodimers composed of the gp120 and gp41 envelope proteins. Each of the gp120 proteins in the trimer is capable of attaching to the CD4 receptor and the coreceptor, and each of the three gp41 units harbors a fusion domain. It is still under debate how many of the envelope subunits within a given trimer have to bind to the CD4 receptors and to the coreceptors, and how many gp41 protein fusion domains are required for fusion. These numbers are referred to as subunit stoichiometries. We present a mathematical framework for estimating these parameters individually by analyzing infectivity assays with pseudotyped viruses. We find that the number of spikes that are engaged in mediating cell entry and the distribution of the spike number play important roles for the estimation of the subunit stoichiometries. Our model framework also shows why it is important to subdivide the question of the number of functional subunits within one trimer into the three different subunit stoichiometries. In a second step, we extend our models to study whether the subunits within one trimer cooperate during receptor binding and fusion. As an example for how our models can be applied, we reanalyze a data set on subunit stoichiometries. We find that two envelope proteins have to engage with CD4-receptors and coreceptors and that two fusion proteins must be revealed within one trimer for viral entry. Our study is motivated by the mechanism of HIV entry but the experimental technique and the model framework can be extended to other viral systems as well.     

Decarboxylation of oxalacetate by pyruvate carboxylase

The decarboxylation of oxalacetate by pyruvate carboxylase in the absence of ADP and Pi is stimulated 400-fold by the presence of oxamate, which is an inhibitory analogue of pyruvate. The observation of substrate inhibition when either oxamate or oxalacetate is varied at a fixed concentration of the other indicates that both molecules bind at the same site on the enzyme. The pH profiles for this reaction show no evidence of the involvement of an enzymic acid-base catalyst, suggesting that the proton and CO2 units may be exchanged directly between the reactants (although CO2 sequestered in the active site may be an intermediate in the process). The pH profiles of the full reverse reaction of pyruvate carboxylase in which oxalacetate decarboxylation is coupled to ATP formation and where Pi is the variable substrate do, however, indicate that such an acid-base catalyst is involved in the other partial reaction of the enzyme in proton transfer to and from biotin. The enzyme also displays two oxamate-independent oxalacetate decarboxylating activities, one of which is biotin-dependent and the other is independent of biotin.     

Atomic force microscopy imaging demonstrates that P2X2 receptors are trimers but that P2X6 receptor subunits do not oligomerize

P2X receptors are cation-selective channels activated by extracellular ATP. The architecture of these receptors is still not completely clear. Here we have addressed this issue by both chemical cross-linking and direct imaging of individual receptors by atomic force microscopy (AFM). Cross-linking of the P2X(2) receptor produced higher order adducts, consistent with the presence of trimers. The mean molecular volume of the receptor determined by AFM (409 nm(3)) also points to a trimeric structure. P2X(2) receptors bearing His(6) epitope tags were incubated with anti-His(6) antibodies, and the resultant complexes were imaged by AFM. For receptors with two bound antibodies, the mean angle between the antibodies was 123 degrees , again indicating that the receptor is a trimer. In contrast, cross-linking of the P2X(6) receptor did not produce higher order adducts, and the mean molecular volume of the receptor was 145 nm(3). We conclude that P2X(2) receptors are trimers, whereas the P2X(6) receptor subunits do not form stable oligomers.