首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 93 毫秒
1.
Liang S  Liu Z  Li W  Ni L  Lai L 《Biopolymers》2000,54(7):515-523
We have developed a strategy for grafting a protein-protein interface based on the known crystal structure of a native ligand and receptor proteins in a complex. The key interaction residues at the ligand protein binding interface are transferred onto a scaffold protein so that the mutated scaffold protein will bind the receptor protein in the same manner as the ligand protein. First, our method identifies key residues and atoms in the ligand protein, which strongly interact with the receptor protein. Second, this method searches the scaffold protein for combinations of candidate residues, among which the distance between any two candidate residues is similar to that between relevant key interaction residues in the ligand protein. These candidate residues are mutated to key interaction residues in the ligand protein respectively. The scaffold protein is superposed onto the ligand protein based upon the coordinates of corresponding atoms, which are assumed to strongly interact with the receptor protein. Complementarity between scaffold and receptor proteins is evaluated. Scaffold proteins with a low superposing rms difference and high complementary score are accepted for further analysis. Then, the relative position of the scaffold protein is adjusted so that the interfaces between the scaffold and receptor proteins have a reasonable packing density. Other mutations are also considered to reduce the desolvation energy or bad steric contacts. Finally, the scaffold protein is cominimized with the receptor protein and evaluated. To test the method, the binding interface of barstar, the inhibitor of barnase, was grafted onto small proteins. Four scaffold proteins with high complementary scores are accepted.  相似文献   

2.
Relaxin, an emerging pharmaceutical treatment for acute heart failure, activates the relaxin family peptide receptor (RXFP1), which is a class A G-protein-coupled receptor. In addition to the classic transmembrane (TM) domain, RXFP1 possesses a large extracellular domain consisting of 10 leucine-rich repeats and an N-terminal low density lipoprotein class A (LDLa) module. Relaxin-mediated activation of RXFP1 requires multiple coordinated interactions between the ligand and various receptor domains including a high affinity interaction involving the leucine-rich repeats and a predicted lower affinity interaction involving the extracellular loops (ELs). The LDLa is essential for signal activation; therefore the ELs/TM may additionally present an interaction site to facilitate this LDLa-mediated signaling. To overcome the many challenges of investigating relaxin and the LDLa module interactions with the ELs, we engineered the EL1 and EL2 loops onto a soluble protein scaffold, mapping specific ligand and loop interactions using nuclear magnetic resonance spectroscopy. Key EL residues were subsequently mutated in RXFP1, and changes in function and relaxin binding were assessed alongside the RXFP1 agonist ML290 to monitor the functional integrity of the TM domain of these mutant receptors. The outcomes of this work make an important contribution to understanding the mechanism of RXFP1 activation and will aid future development of small molecule RXFP1 agonists/antagonists.  相似文献   

3.
Pheromone‐binding proteins (PBPs) are distributed widely on the antennae of insects, and they are believed to be involved in the process of chemical signal transduction, but their interaction with chemicals is largely unknown. Here, we present our findings on the key amino acid residues of PBPs in the gypsy moth, Lymantria dispar. Potential key residues were screened with the Calculate Mutation Energy program and molecular docking methods. Mutated proteins were obtained by mutating residues to alanine via site‐directed mutagenesis. Circular dichroism (CD) spectroscopy showed that the mutated proteins formed α‐helix, and the stability of protein structure was influenced due to mutations. Fluorescence binding assays were further conducted with the mutated proteins, sex pheromones and analogues. Results showed that to PBP 1, tyrosine at position 41 and phenylalanine at position 76 could be the key amino acid residues influencing the stability of structure; in addition, phenylalanine at 36 and lysine at position 94 could be key amino acid residues interacting with chemicals. To PBP 2, glycine at position 49, phenylalanine at position 76 and lysine at position 121 could be the key amino acid residues in the structural stability. These results shed light on the relationship between the specific amino acids and functions of PBPs in transmitting the chemical signals.  相似文献   

4.
Kinase Suppressor of Ras (KSR) is a molecular scaffold that interacts with the core kinase components of the ERK cascade, Raf, MEK, ERK to provide spatial and temporal regulation of Ras-dependent ERK cascade signaling. Interruption of this mechanism can have a high influence in inhibiting the downstream signaling of the mutated tyrosine kinase receptor kinase upon ligand binding. Still none of the studies targeted to prevent the binding of Raf, MEK binding on kinase suppressor of RAS. In that perspective the cysteine rich C1 domain of scaffold proteins kinase suppressor of Ras-1 was targeted rather than its ATP binding site with small ligand molecules like flavones and anthocyanidins and analyzed through insilico docking studies. The binding energy evaluation shows the importance of hydroxyl groups at various positions on the flavone and anthocyanidin nucleus. Over all binding interaction shows these ligands occupied the potential sites of cysteine rich C1 domain of scaffold protein KSR.  相似文献   

5.
Hydrogen bond, hydrophobic and vdW interactions are the three major non-covalent interactions at protein–protein interfaces. We have developed a method that uses only these properties to describe interactions between proteins, which can qualitatively estimate the individual contribution of each interfacial residue to the binding and gives the results in a graphic display way. This method has been applied to analyze alanine mutation data at protein–protein interfaces. A dataset containing 13 protein–protein complexes with 250 alanine mutations of interfacial residues has been tested. For the 75 hot-spot residues (G1.5 kcal mol-1), 66 can be predicted correctly with a success rate of 88%. In order to test the tolerance of this method to conformational changes upon binding, we utilize a set of 26 complexes with one or both of their components available in the unbound form. The difference of key residues exported by the program is 11% between the results using complexed proteins and those from unbound ones. As this method gives the characteristics of the binding partner for a particular protein, in-depth studies on protein–protein recognition can be carried out. Furthermore, this method can be used to compare the difference between protein–protein interactions and look for correlated mutation. Figure Key interaction grids at the interface between barnase and barstar. Key interaction grid for barnase and barstar are presented in one figure according to their coordinates. In order to distinguish the two proteins, different icons were assigned. Crosses represent key grids for barstar and dots represent key grids for barnase. The four residues in ball and stick are Asp40 in barstar and Arg83, Arg87, His102 in barnase.  相似文献   

6.
The ciliate Euplotes raikovi produces a family of diffusible signal proteins (pheromones) that function as prototypic growth factors. They may either promote cell growth, by binding to pheromone receptors synthesized by the same cells from which they are secreted (autocrine activity), or induce a temporary cell shift from the growth stage to a mating (sexual) one by binding to pheromone receptors of other, conspecific cells (paracrine activity). In cells constitutively secreting the pheromone Er-1, it was first observed that the expression of the Er-1 receptor "p15," a type II membrane protein of 130 amino acids, is quantitatively correlated with the extracellular concentration of secreted pheromone. p15 expression on the cell surface rapidly and markedly increased after the removal of secreted Er-1 and gradually decreased in parallel with new Er-1 secretion. It was then shown that p15 is internalized through endocytic vesicles following Er-1 binding and that the internalization of p15/Er-1 complexes is specifically blocked by the paracrine p15 binding of Er-2, a pheromone structurally homologous to, and thus capable of fully antagonizing, Er-1. Based on previous findings that the p15 pheromone-binding site is structurally equivalent to Er-1 and that Er-1 molecules polymerize in crystals following a pattern of cooperative interaction, it was proposed that p15/Er-1 complexes are internalized as a consequence of their unique property (not shared by p15/Er-2 complexes) of undergoing clustering.  相似文献   

7.
We used a novel charge optimization technique to study the small ribonuclease barnase and to analyze its interaction with a natural tight binding inhibitor, the protein barstar. The approach uses a continuum model to explicitly determine the charge distributions that lead to the most favorable electrostatic contribution to binding when competing desolvation and interaction effects are included. Given its backbone fold, barstar is electrostatically optimized for tight binding to barnase when compared with mutants where residues have been substituted with one of the 20 common amino acids. Natural proteins thus appear to use optimization of electrostatic interactions as one strategy for achieving tight binding.  相似文献   

8.
Birch-allergic patients often experience oral allergy syndrome upon ingestion of vegetables and fruits, most prominently apple, that is caused by antibody cross-reactivity of the IgE antibodies in patients to proteins sharing molecular surface structures with the major birch pollen group 1 allergen from Betula verrucosa (Bet v 1). Still, to what extent two molecular surfaces need to be similar for clinically relevant antibody cross-reactivity to occur is unknown. Here, we describe the grafting of a defined conformational antibody epitope from Bet v 1 onto the surface of the homologous apple allergen Malus domestica group 1 (Mal d 1). Engineering of the epitope was accomplished by genetic engineering substituting amino acid residues in Mal d 1 differing between Bet v 1 and Mal d 1 within the epitope defined by the mAb BV16. The kinetic parameters characterizing the antibody binding interaction to Bet v 1 and to the mutated Mal d 1 variant, respectively, were assessed by Biacore experiments demonstrating indistinguishable binding kinetics. This demonstrates that a conformational epitope defined by a high affinity antibody-allergen interaction can successfully be grafted onto a homologous scaffold molecule without loss of epitope functionality. Furthermore, we show that increasing surface similarity to Bet v 1 of Mal d 1 variants by substitution of 6-8 residues increased the ability to trigger basophil histamine release with blood from birch-allergic patients not responding to natural Mal d 1. Conversely, reducing surface similarity to Bet v 1 of a Mal d 1 variant by substitution of three residues abolished histamine release in one patient reacting to Mal d 1.  相似文献   

9.
The fibroblast growth factor receptor 1 (flg) contains eight acidic amino acids between the first and second immunoglobulin domain. This report examines the role of the acidic domain in the interaction of the flg receptor with its ligands. We observed a marked inhibition of binding of bFGF to the receptor when the acidic domain was completely deleted, but mutants with two and four amino acids deleted (flgΔ2 and flgΔ4, respectively) still bound the ligand. After addition of a bifunctional cross-linking reagent, cross-linked complexes (between bFGF and receptor) with the expected size were observed in cells expressing mutants lacking two or four acidic residues, but not in cells expressing mutants lacking six or eight acidic residues. Immunoprecipitation with anti-flg antibody followed by electrophoresis produced a band of 90 Kd in tunicamycin-treated cells expressing the mutant as well as the wild-type receptors, indicating that the inhibition of binding was not due to defective expression of the protein. The ability of flgΔ8 to mediate a mitogenic response to FGFs was also greatly reduced when this mutated receptor was expressed in receptor-negative cells. The effect of replacing the acidic amino acids with lysine residues was also studied. Binding of bFGF to cells transfected with a plasmid encoding a mutated protein with four amino acid substitutions was totally inhibited, but an eight amino acid substitution did not alter ligand binding to the receptor. In this case the mutation with four amino acids substitution caused a drastic impairment of protein expression. Thus the acidic domain of the FGFR-1 plays an essential role in receptor function, either because it is important for a stable protein configuration or for ligand-receptor interaction. © 1993 Wiley-Liss, Inc.  相似文献   

10.
Monocyte chemoattractant protein-1 (MCP-1) is a chemotactic cytokine mainly acting on monocytes and T cells that elicits its biological effects by interacting with the seven-transmembrane helix receptor CCR2B. The vaccinia virus strain Lister and many other poxviruses express soluble proteins (vCCI) that bind MCP-1 and other CC chemokines and inhibit their function. In order to define the interaction site of MCP-1 with vCCI from vaccinia, surface exposed residues of MCP-1 were identified and mutated to alanine. The MCP-1 variants were expressed, purified, and their interaction with vCCI was characterized. The site on MCP-1 for vCCI binding is dominated by arginine 18 with important additional contributions from tyrosine 13 and arginine 24. These residues define a binding site that largely overlaps with the CCR2B receptor interaction site. The viral chemokine-binding protein vCCI thus inhibits the biological function of MCP-1 by directly masking its CCR2B receptor-binding site.  相似文献   

11.
Computational grafting of functional motifs onto scaffold proteins is a promising way to engineer novel proteins with pre-specified functionalities. Typically, protein grafting involves the transplantation of protein side chains from a functional motif onto structurally homologous regions of scaffold proteins. Using this approach, we previously transplanted the human immunodeficiency virus 2F5 and 4E10 epitopes onto heterologous proteins to design novel "epitope-scaffold" antigens. However, side-chain grafting is limited by the availability of scaffolds with compatible backbone for a given epitope structure and offers no route to modify backbone structure to improve mimicry or binding affinity. To address this, we report here a new and more aggressive computational method-backbone grafting of linear motifs-that transplants the backbone and side chains of linear functional motifs onto scaffold proteins. To test this method, we first used side-chain grafting to design new 2F5 epitope scaffolds with improved biophysical characteristics. We then independently transplanted the 2F5 epitope onto three of the same parent scaffolds using the newly developed backbone grafting procedure. Crystal structures of side-chain and backbone grafting designs showed close agreement with both the computational models and the desired epitope structure. In two cases, backbone grafting scaffolds bound antibody 2F5 with 30- and 9-fold higher affinity than corresponding side-chain grafting designs. These results demonstrate that flexible backbone methods for epitope grafting can significantly improve binding affinities over those achieved by fixed backbone methods alone. Backbone grafting of linear motifs is a general method to transplant functional motifs when backbone remodeling of the target scaffold is necessary.  相似文献   

12.
Binding of inositol 1,4,5-trisphosphate (IP(3)) to the amino-terminal region of IP(3) receptor promotes Ca(2+) release from the endoplasmic reticulum. Within the amino terminus, the first 220 residues directly preceding the IP(3) binding core domain play a key role in IP(3) binding suppression and regulatory protein interaction. Here we present a crystal structure of the suppressor domain of the mouse type 1 IP(3) receptor at 1.8 A. Displaying a shape akin to a hammer, the suppressor region contains a Head subdomain forming the beta-trefoil fold and an Arm subdomain possessing a helix-turn-helix structure. The conserved region on the Head subdomain appeared to interact with the IP(3) binding core domain and is in close proximity to the previously proposed binding sites of Homer, RACK1, calmodulin, and CaBP1. The present study sheds light onto the mechanism underlying the receptor's sensitivity to the ligand and its communication with cellular signaling proteins.  相似文献   

13.
Recombinant ligands derived from small protein scaffolds show promise as robust research and diagnostic reagents and next generation protein therapeutics. Here, we derived high-affinity binders of human interferon gamma (hIFNγ) from the three helix bundle scaffold of the albumin-binding domain (ABD) of protein G from Streptococcus G148. Computational interaction energy mapping, solvent accessibility assessment, and in silico alanine scanning identified 11 residues from the albumin-binding surface of ABD as suitable for randomization. A corresponding combinatorial ABD scaffold library was synthesized and screened for hIFNγ binders using in vitro ribosome display selection, to yield recombinant ligands that exhibited K(d) values for hIFNγ from 0.2 to 10 nM. Molecular modeling, computational docking onto hIFNγ, and in vitro competition for hIFNγ binding revealed that four of the best ABD-derived ligands shared a common binding surface on hIFNγ, which differed from the site of human IFNγ receptor 1 binding. Thus, these hIFNγ ligands provide a proof of concept for design of novel recombinant binding proteins derived from the ABD scaffold.  相似文献   

14.
The alpha/beta‐hydrolases are a family of acid‐base‐nucleophile catalytic triad enzymes with a common fold, but using a wide variety of substrates, having different pH optima, catalyzing unique catalytic reactions and often showing improved chemical and thermo stability. The ABH enzymes are prime targets for protein engineering. Here, we have classified active sites from 51 representative members of 40 structural ABH fold families into eight distinct conserved geometries. We demonstrate the occurrence of a common structural motif, the catalytic acid zone, at the catalytic triad acid turn. We show that binding of an external ligand does not change the structure of the catalytic acid zone and both the ligand‐free and ligand‐bound forms of the protein belong to the same catalytic acid zone subgroup. We also show that the catalytic acid zone coordinates the position of the catalytic histidine loop directly above its plane, and consequently, fixes the catalytic histidine in a proper position near the catalytic acid. Finally, we demonstrate that the catalytic acid zone plays a key role in multi‐subunit complex formation in ABH enzymes, and is involved in interactions with other proteins. As a result, we speculate that each of the catalytic triad residues has its own supporting structural scaffold, similar to the catalytic acid zone, described above, which together form the extended catalytic triad motif. Each scaffold coordinates the function of its respective catalytic residue, and can even compensate for the loss of protein function, if the catalytic amino acid is mutated.  相似文献   

15.
Ubiquitin, a post-translational protein modifier inside the cell, functions as a CXC chemokine receptor (CXCR) 4 agonist outside the cell. However, the structural determinants of the interaction between extracellular ubiquitin and CXCR4 remain unknown. Utilizing C-terminal truncated ubiquitin and ubiquitin mutants, in which surface residues that are known to interact with ubiquitin binding domains in interacting proteins are mutated (Phe-4, Leu-8, Ile-44, Asp-58, Val-70), we provide evidence that the ubiquitin-CXCR4 interaction follows a two-site binding mechanism in which the hydrophobic surfaces surrounding Phe-4 and Val-70 are important for receptor binding, whereas the flexible C terminus facilitates receptor activation. Based on these findings and the available crystal structures, we then modeled the ubiquitin-CXCR4 interface with the RosettaDock software followed by small manual adjustments, which were guided by charge complementarity and anticipation of a conformational switch of CXCR4 upon activation. This model suggests three residues of CXCR4 (Phe-29, Phe-189, Lys-271) as potential interaction sites. Binding studies with HEK293 cells overexpressing wild type and CXCR4 after site-directed mutagenesis confirm that these residues are important for ubiquitin binding but that they do not contribute to the binding of stromal cell-derived factor 1α. Our findings suggest that the structural determinants of the CXCR4 agonist activity of ubiquitin mimic the typical structure-function relationship of chemokines. Furthermore, we provide evidence for separate and specific ligand binding sites on CXCR4. As exogenous ubiquitin has been shown to possess therapeutic potential, our findings are expected to facilitate the structure-based design of new compounds with ubiquitin-mimetic actions on CXCR4.  相似文献   

16.
The cyclic AMP-dependent kinase-anchoring proteins (AKAPs) function as scaffolds for a wide-range of protein-protein interactions. The 250-kDa AKAP known as gravin plays a central role in organizing G-protein-coupled receptors to the protein kinases and phosphatases that regulate receptor function in desensitization, resensitization, and sequestration. Although gravin is critical for G-protein-linked receptor biology, the molecular features of the receptor necessary for interaction with this scaffold are not known. Herein, we map the regions of the beta(2)-adrenergic receptor that are required for binding to gravin. Intracellular loops 1, 2, and 3 appear not to participate in the binding of the receptor to the scaffold. In contrast, the C-terminal cytoplasmic region of the receptor (Arg-329 to Leu-413) competes readily for the binding of the beta(2)-adrenergic receptor by gravin, both using in vitro and in vivo assays. C-terminally truncated peptides with sequences ranging from Arg-329 to Leu-342 (13 aminoacyl residues), to Asn-352 (23 residues), to Tyr-366 (37 residues), to Asp-380 (51 residues), or to His-390 (61 residues), as well as N-terminally truncated peptides from Gln-391 to Leu-413 (23 residues) or Leu-381 to Leu-413 (33 residues) displayed no ability to block binding of receptor to gravin. The combination of Arg-329 to His-390 peptide and Gln-391 to Leu-413 peptide, however, reconstitutes a fragmented but full-length C-terminal region and also potently blocks the ability of gravin to bind the beta(2)-adrenergic receptor. The gravin-receptor interaction was examined in response to agonist by confocal microscopy. Remarkably, the association of the receptor with gravin was not disrupted during agonist-induced sequestration. The receptor-scaffold complex was maintained during agonist-induced sequestration. These data, in agreement with the biochemical data, reveal that gravin binds the receptor through the beta(2)-adrenergic receptor C-terminal cytoplasmic domain and that this interaction is maintained as the receptor is internalized. This is the first report of an AKAP scaffold protein translocating with its receptor, in this case a G-protein-coupled receptor.  相似文献   

17.
The 5-HT(3) receptor is a member of the Cys-loop family of ligand-gated ion channels. The extracellular domains of these proteins contain six amino acid loops (A-F) that converge to form the ligand binding site. In this study we have mutated 21 residues in or close to the 5-HT(3) receptor F-loop (Ile(192) to Gly(212)) to Ala or to a residue with similar chemical properties. Mutant receptors were expressed in HEK293 cells, and binding affinity was measured using [(3)H]granisetron. Two regions displayed decreases in binding affinity when mutated to Ala (Ile(192)-Arg(196) and Asp(204)-Ser(206)), but only one region was sensitive when mutated to chemically similar residues (Ile(192)-Val(201)). Homology modeling using acetylcholine-binding protein crystal structures with a variety of different bound ligands suggests there may be distinct movements of Trp(195) and Asp(204) upon ligand binding, indicating that these residues and their immediate neighbors have the ability to interact differently with different ligands. The models suggest predominantly lateral movement around Asp(204) and rotational movement around Trp(195), indicating the former is in a more flexible region. Overall our results are consistent with a flexible 5-HT(3) receptor F-loop with two regions that have specific but distinct roles in ligand binding.  相似文献   

18.
Protein-protein interactions are very important in the function of a cell. Computational studies of these interactions have been of interest, but often they have utilized classical modelling techniques. In recent years, quantum mechanical (QM) treatment of entire proteins has emerged as a powerful approach to study biomolecular systems. Herein, we apply a semi-empirical divide and conquer (DC) methodology coupled with a dielectric continuum model for the solvent, to explore the contribution of electrostatics, polarization and charge transfer to the interaction energy between barnase and barstar in their complex form. Molecular dynamic (MD) simulation was performed to account for the dynamic behavior of the complex. The results show that electrostatics, charge transfer and polarization favor the formation of the complex. Our study shows that electrostatics dominates the interaction between barnase and barstar ( approximately 73%), while charge transfer and polarization are approximately 21% and approximately 6%, respectively. Close inspection of the polarization and charge-transfer effects on the charge distribution of the complex reveals the existence of two, well localized, regions in barstar. The first region includes the residues between P27 and Y47 and the second region is between N65 and D83. Since no such regions could be detected in barnase clearly suggests that barstar is well optimized for efficiently binding barnase. Furthermore, using our interaction energy decomposition scheme, we were able to identify all residues that have been experimentally determined to be important for the complex formation and to suggest other residues never have been investigated. This suggests that our approach will be useful as an aid in further understanding protein-protein contacts for the ultimate goal to produce successful inhibitors for protein complexes.  相似文献   

19.
Discrete peptide domains within the primary sequence of cell-surface receptor glycoproteins are believed to regulate not only their function but also their targeting to the cell membrane. To identify sequence elements required for intracellular transport and ligand binding by the human Tac interleukin-2 (IL-2) receptor, we prepared expression plasmids encoding a series of artificially mutated or naturally occurring variants of the Tac cDNA. In particular, we sought to further delineate the functional role of the sequences contributed by each of the eight exons that together encode the Tac protein. Deletion of exons 5 through 8 of the receptor had no detectable effect on IL-2 binding or intracellular transport of the Tac protein, and resulted in secreted forms of this IL-2-binding protein. Removal of sequences corresponding to all of exon 4 ablated IL-2 binding activity yet still permitted transport to the cell surface. In contrast, partial deletion of exon 4 sequences resulted in proteins that not only lacked IL-2 binding activity but also were sequestered within the endoplasmic reticulum. Removal of one or both of the N-linked glycosylation sites present in the Tac protein did not impair receptor transport or ligand binding. These results demonstrate that exon 4 of the Tac gene encodes amino acid residues that play an important role in regulating both the intracellular transport and function of this IL-2 receptor.  相似文献   

20.
We have defined regions in the Sonic hedgehog (Shh) molecule that are important for Patched (Ptc) receptor binding by targeting selected surface amino acid residues with probes of diverse sizes and shapes and assessing the effects of these modifications on function. Eleven amino acid residues that surround the surface of the protein were chosen for these studies and mutated to cysteine residues. These cysteines were then selectively modified with thiol-specific probes, and the modified proteins were tested for hedgehog receptor binding activity and their ability to induce differentiation of C3H10T1/2 cells into osteoblasts. Based on these analyses, approximately one-third of the Shh surface can be modified without effect on function regardless of the size of the attachment. These sites are located near to where the C terminus protrudes from the surface of the protein. All other sites were sensitive to modification, indicating that the interaction of Shh with its primary receptor Ptc is mediated over a large surface of the Shh protein. For sites Asn-50 and Ser-156, function was lost with the smallest of the probes tested, indicating that these residues are in close proximity to the Ptc-binding site. The epitope for the neutralizing mAb 5E1 mapped to a close but distinct region of the structure. The structure-activity data provide a unique view of the interactions between Shh and Ptc that is not readily attainable by conventional mapping strategies.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号