首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
The two most buried carboxyl groups in ribonuclease Sa (RNase Sa) are Asp33 (99% buried; pK 2.4) and Asp79 (85% buried; pK 7.4). Above these pK values, the stability of the D33A variant is 6kcal/mol less than wild-type RNase Sa, and the stability of the D79A variant is 3.3kcal/mol greater than wild-type RNase Sa. The key structural difference between the carboxyl groups is that Asp33 forms three intramolecular hydrogen bonds, and Asp79 forms no intramolecular hydrogen bond. Here, we focus on Asp79 and describe studies of 11 Asp79 variants. Most of the variants were at least 2kcal/mol more stable than wild-type RNase Sa, and the most interesting was D79F. At pH 3, below the pK of Asp79, RNase Sa is 0.3kcal/mol more stable than the D79F variant. At pH 8.5, above the pK of Asp79, RNase Sa is 3.7kcal/mol less stable than the D79F variant. The unfavorable contribution of Asp79 to the stability appears to result from the Born self-energy of burying the charge and, more importantly, from unfavorable charge-charge interactions. To counteract the effect of the negative charge on Asp79, we prepared the Q94K variant and the crystal structure showed that the amino group of the Lys formed a hydrogen-bonded ion pair (distance, 2.71A; angle, 100 degrees ) with the carboxyl group of Asp79. The stability of the Q94K variant was about the same as the wild-type at pH 3, where Asp79 is uncharged, but 1kcal/mol greater than that of wild-type RNase Sa at pH 8.5, where Asp79 is charged. Differences in hydrophobicity, steric strain, Born self-energy, and electrostatic interactions all appear to contribute to the range of stabilities observed in the variants. When it is possible, replacing buried, non-hydrogen bonded, ionizable side-chains with non-polar side-chains is an excellent means of increasing protein stability.  相似文献   

2.
The aim of this study was to gain a better understanding of the contribution of hydrogen bonds by tyrosine -OH groups to protein stability. The amino acid sequences of RNases Sa and Sa3 are 69 % identical and each contains eight Tyr residues with seven at equivalent structural positions. We have measured the stability of the 16 tyrosine to phenylalanine mutants. For two equivalent mutants, the stability increases by 0.3 kcal/mol (RNase Sa Y30F) and 0.5 kcal/mol (RNase Sa3 Y33F) (1 kcal=4.184 kJ). For all of the other mutants, the stability decreases with the greatest decrease being 3.6 kcal/mol for RNase Sa Y52F. Seven of the 16 tyrosine residues form intramolecular hydrogen bonds and the average decrease in stability for these is 2.0(+/-1.0) kcal/mol. For the nine tyrosine residues that do not form intramolecular hydrogen bonds, the average decrease in stability is 0.4(+/-0.6) kcal/mol. Thus, most tyrosine -OH groups contribute favorably to protein stability even if they do not form intramolecular hydrogen bonds. Generally, the stability changes for equivalent positions in the two proteins are remarkably similar. Crystal structures were determined for two of the tyrosine to phenylalanine mutants of RNase Sa: Y80F (1.2 A), and Y86F (1.7 A). The structures are very similar to that of wild-type RNase Sa, and the hydrogen bonding partners of the tyrosine residues always form intermolecular hydrogen bonds to water in the mutants. These results provide further evidence that the hydrogen bonding and van der Waals interactions of polar groups in the tightly packed interior of folded proteins are more favorable than similar interactions with water in the unfolded protein, and that polar group burial makes a substantial contribution to protein stability.  相似文献   

3.
Our goal was to gain a better understanding of how protein stability can be increased by improving β‐turns. We studied 22 β‐turns in nine proteins with 66–370 residues by replacing other residues with proline and glycine and measuring the stability. These two residues are statistically preferred in some β‐turn positions. We studied: Cold shock protein B (CspB), Histidine‐containing phosphocarrier protein, Ubiquitin, Ribonucleases Sa2, Sa3, T1, and HI, Tryptophan synthetase α‐subunit, and Maltose binding protein. Of the 15 single proline mutations, 11 increased stability (Average = 0.8 ± 0.3; Range = 0.3–1.5 kcal/mol), and the stabilizing effect of double proline mutants was additive. On the basis of this and our previous work, we conclude that proteins can generally be stabilized by replacing nonproline residues with proline residues at the i + 1 position of Type I and II β‐turns and at the i position in Type II β‐turns. Other turn positions can sometimes be used if the φ angle is near ?60° for the residue replaced. It is important that the side chain of the residue replaced is less than 50% buried. Identical substitutions in β‐turns in related proteins give similar results. Proline substitutions increase stability mainly by decreasing the entropy of the denatured state. In contrast, the large, diverse group of proteins considered here had almost no residues in β‐turns that could be replaced by Gly to increase protein stability. Improving β‐turns by substituting Pro residues is a generally useful way of increasing protein stability. Proteins 2009. © 2009 Wiley‐Liss, Inc.  相似文献   

4.
Gaining a better understanding of the denatured state ensemble of proteins is important for understanding protein stability and the mechanism of protein folding. We studied the folding kinetics of ribonuclease Sa (RNase Sa) and a charge-reversal variant (D17R). The refolding kinetics are similar, but the unfolding rate constant is 10-fold greater for the variant. This suggests that charge-charge interactions in the denatured state and the transition state ensembles are more favorable in the variant than in RNase Sa, and shows that charge-charge interactions can influence the kinetics and mechanism of protein folding.  相似文献   

5.
Escherichia coli RNase HI is a well-characterized model system for protein folding and stability. Controlling protein stability is critical for both natural proteins and for the development of engineered proteins that function under extreme conditions. We have used native-state hydrogen exchange on a variant containing the stabilizing mutation Asp10 to alanine in order to determine its residue-specific stabilities. On average, the DeltaG(unf) value for each residue was increased by 2-3 kcal/mol, resulting in a lower relative population of partially unfolded forms. Though increased in stability by a uniform factor, D10A shows a distribution of stabilities in its secondary structural units that is similar to that of E. coli RNase H, but not the closely related protein from Thermus thermophilus. Hence, the simple mutation used to stabilize the enzyme does not recreate the balance of conformational flexibility evolved in the thermophilic protein.  相似文献   

6.
It is known that several naturally occurring substances known as osmolytes increase the conformational stability of proteins. Bolen and co-worker proposed the osmophobic theory, which asserts the osmolyte effect occurs because of an unfavorable interaction of osmolytes mainly with the protein backbone, based on the results on the transfer Gibbs energy of amino acids (Deltag) [Bolen and Baskakov (2001) J. Mol. Biol. 310, 955-963]. In this paper, we report the effect of sarcosine on the conformational stability (DeltaG) of RNase Sa (96 residues and one disulfide bond) and four mutant proteins. The thermal denaturation curves for RNase Sa in sarcosine fitted a two-state model on nonlinear least-squares analysis. All the RNase Sa proteins were stabilized by sarcosine. For example, the increase in stability of the wild-type protein in 4 M sarcosine due to the osmolyte effect (Delta(o)DeltaG) is 3.2 kcal/mol. Mutational analysis of the osmolyte effect indicated that the changed Delta(o)DeltaG values upon mutation (Delta(m)Delta(o)DeltaG), as estimated from the Deltag values, are similar to the experimental values. Structural-based analysis of the osmolyte effect was also performed using model denatured structures: (a) a fully extended model (single chain) with no disulfide bond, (b) two-part, unfolded models (two chains) with a disulfide bond constructed through molecular dynamic (MD) simulation, and (c) a two-part, folded model (two chains). The two-part, unfolded models were expected to be more suitable as denatured structures. The Delta(o)DeltaG values calculated using the two-part, unfolded models were more consistent with experimental values than those calculated using the fully extended and two-part, folded models. This suggests that MD simulation is useful for testing denatured structures. These results indicate that the osmophobic theory can explain the osmolyte effect on protein stability.  相似文献   

7.
Characterizing the denatured state ensemble is crucial to understanding protein stability and the mechanism of protein folding. The aim of this research was to see if fluorescence could be used to gain new information on the denatured state ensemble. Ribonuclease Sa (RNase Sa) contains no Trp residues. We made five variants of RNase Sa by adding Trp residues at locations where they are found in other members of the microbial ribonuclease family. To better understand the protein denatured state, we also studied the fluorescence properties of the following peptides: N-acetyl-Trp-amide (NATA), N-acetyl-Ala-Trp-Ala-amide (AWA), N-acetyl-Ala-Ala-Trp-Ala-Ala-amide (AAWAA), and the five pentapeptides with the same sequence as the Trp substitution sites in RNase Sa. The major conclusions are: 1), the wavelength of maximum fluorescence intensity, λmax, does not differ significantly for the peptides and the denatured proteins; 2), the fluorescence intensity at λmax, IF, differs significantly for the five Trp containing variants of RNase Sa; 3), the IF differences for the denatured proteins are mirrored in the peptides, showing that the short-range effects giving rise to the IF differences in the peptides are also present in the proteins; 4) the IF values for the denatured proteins are more than 30% greater than for the peptides, showing the presence of long-range effects in the proteins; 5), fluorescence quenching of Trp by acrylamide and iodide is more than 50% greater in the peptides than in the denatured proteins, showing that long-range effects limit the accessibility of the quenchers to the Trp side chains in the proteins; and 6), these results show that nonlocal effects in the denatured states of proteins influence Trp fluorescence and accessibility significantly.  相似文献   

8.
We previously suggested that proteins gain more stability from the burial and hydrogen bonding of polar groups than from the burial of nonpolar groups (Pace, C. N. (2001) Biochemistry 40, 310-313). To study this further, we prepared eight Thr-to-Val mutants of RNase Sa, four in which the Thr side chain is hydrogen-bonded and four in which it is not. We measured the stability of these mutants by analyzing their thermal denaturation curves. The four hydrogen-bonded Thr side chains contribute 1.3 +/- 0.9 kcal/mol to the stability; those that are not still contribute 0.4 +/- 0.9 kcal/mol to the stability. For 40 Thr-to-Val mutants of 11 proteins, the average decrease in stability is 1.0 +/- 1.0 kcal/mol when the Thr side chain is hydrogen-bonded and 0.0 +/- 0.5 kcal/mol when it is not. This is clear evidence that hydrogen bonds contribute favorably to protein stability. In addition, we prepared four Val-to-Thr mutants of RNase Sa, measured their stability, and determined their crystal structures. In all cases, the mutants are less stable than the wild-type protein, with the decreases in stability ranging from 0.5 to 4.4 kcal/mol. For 41 Val-to-Thr mutants of 11 proteins, the average decrease in stability is 1.8 +/- 1.3 kcal/mol and is unfavorable for 40 of 41 mutants. This shows that placing an [bond]OH group at a site designed for a [bond]CH3 group is very unfavorable. So, [bond]OH groups can contribute favorably to protein stability, even if they are not hydrogen-bonded, if the site was selected for an [bond]OH group, but they will make an unfavorable contribution to stability, even if they are hydrogen-bonded, when they are placed at a site selected for a [bond]CH3 group. The contribution that polar groups make to protein stability depends strongly on their environment.  相似文献   

9.
10.
Abstract

Recent site-directed mutagenesis and thermodynamic studies have shown that the V74I mutant of Escherichia coli ribonuclease HI (RNase HI) is more stable than the wild type protein [Ishikawa et al., Biochemistry 32, 6171 (1993)]. In order to clarify the stabilization mechanism of this mutant, we calculated the free energy change due to the mutation Val 74→Ile in both the native and denatured states by free energy perturbations based on molecular dynamics (MD) simulations. We carried out inclusive MD simulations for the protein in water; i.e., fully solvated, no artificial constraints applied, and all long-range Coulomb interactions included. We found that the free energy of the mutant increased slightly relative to the wild type, in the native state by 1.60 kcal/mol, and in the denatured state by 2.25 kcal/mol. The unfolding free energy increment of the mutant (0.66 ± 0.19 kcal/mol) was in good agreement with the experimental value (0.6 kcal/mol). The hysteresis error in the free energy calculations, i.e., forward and reverse perturbations, was only ±0.19 kcal/mol. These results show that the V74I mutant is stabilized relative to the wild type by the increased free energy of the denatured state and not by a decrease in the free energy of the native state as had been proposed earlier based on the mutant X-ray structure. It was found that the stabilization was caused by a loss of solvation energy in the mutant denatured state and not by improved packing interactions inside the native protein.  相似文献   

11.
Recent work on the thermodynamics of protein denatured states is providing insight into the stability of residual structure and the conformational constraints that affect the disordered states of proteins. Current data from native state hydrogen exchange and the pH dependence of protein stability indicate that residual structure can modulate the stability of the denatured state by up to 4 kcal mol(-1). NMR structural data have emphasized the role of hydrophobic clusters in stabilizing denatured state residual structures, however recent results indicate that electrostatic interactions, both favorable and unfavorable, are also important modulators of the stability of the denatured state. Thermodynamics methods that take advantage of histidine-heme ligation chemistry have also been developed to probe the conformational constraints that act on denatured states. These methods have provided insights into the role of excluded volume, chain stiffness, and loop persistence in modulating the conformational preferences of highly disordered proteins. New insights into protein folding and novel methods to manipulate protein stability are emerging from this work.  相似文献   

12.
It is difficult to increase protein stability by adding hydrogen bonds or burying nonpolar surface. The results described here show that reversing the charge on a side chain on the surface of a protein is a useful way of increasing stability. Ribonuclease T1 is an acidic protein with a pI approximately 3.5 and a net charge of approximately -6 at pH 7. The side chain of Asp49 is hyperexposed, not hydrogen bonded, and 8 A from the nearest charged group. The stability of Asp49Ala is 0.5 kcal/mol greater than wild-type at pH 7 and 0.4 kcal/mol less at pH 2.5. The stability of Asp49His is 1.1 kcal/mol greater than wild-type at pH 6, where the histidine 49 side chain (pKa = 7.2) is positively charged. Similar results were obtained with ribonuclease Sa where Asp25Lys is 0.9 kcal/mol and Glu74Lys is 1.1 kcal/mol more stable than the wild-type enzyme. These results suggest that protein stability can be increased by improving the coulombic interactions among charged groups on the protein surface. In addition, the stability of RNase T1 decreases as more hydrophobic aromatic residues are substituted for Ala49, indicating a reverse hydrophobic effect.  相似文献   

13.
A Giletto  C N Pace 《Biochemistry》1999,38(40):13379-13384
The side-chain carboxyl of Asp 76 in ribonuclease T1 (RNase T1) is buried, charged, non-ion-paired, and forms three good intramolecular hydrogen bonds (2.63, 2.69, and 2.89 A) and a 2.66 A hydrogen bond to a buried, conserved water molecule. When Asp 76 was replaced by Asn, Ser, and Ala, the conformational stability of the protein decreased by 3.1, 3.2, and 3.7 kcal/mol, respectively. The stability was measured as a function of pH for wild-type RNase T1 and the D76N mutant and showed that the pH dependence below pH 3 was almost entirely due to Asp 76. The pK of Asp 76 is 0.5 in the native state and 3.7 in the denatured state. Thus, the hydrogen bonding of the carboxyl group of Asp 76 contributes more than half of the net stability of RNase T1 at pH 7. In addition, the charged carboxyl of Asp 76 stabilizes structure in the denatured states of RNase T1 that is not present in D76N, D76S, and D76A.  相似文献   

14.
Our goal was to gain a better understanding of the contribution of the burial of polar groups and their hydrogen bonds to the conformational stability of proteins. We measured the change in stability, Δ(ΔG), for a series of hydrogen bonding mutants in four proteins: villin headpiece subdomain (VHP) containing 36 residues, a surface protein from Borrelia burgdorferi (VlsE) containing 341 residues, and two proteins previously studied in our laboratory, ribonucleases Sa (RNase Sa) and T1 (RNase T1). Crystal structures were determined for three of the hydrogen bonding mutants of RNase Sa: S24A, Y51F, and T95A. The structures are very similar to wild type RNase Sa and the hydrogen bonding partners form intermolecular hydrogen bonds to water in all three mutants. We compare our results with previous studies of similar mutants in other proteins and reach the following conclusions. (1) Hydrogen bonds contribute favorably to protein stability. (2) The contribution of hydrogen bonds to protein stability is strongly context dependent. (3) Hydrogen bonds by side chains and peptide groups make similar contributions to protein stability. (4) Polar group burial can make a favorable contribution to protein stability even if the polar groups are not hydrogen bonded. (5) The contribution of hydrogen bonds to protein stability is similar for VHP, a small protein, and VlsE, a large protein.  相似文献   

15.
Ribonuclease Sa (RNase Sa) contains no tryptophan (Trp) residues. We have added single Trp residues to RNase Sa at sites where Trp is found in four other microbial ribonucleases, yielding the following variants of RNase Sa: Y52W, Y55W, T76W, and Y81W. We have determined crystal structures of T76W and Y81W at 1.1 and 1.0 A resolution, respectively. We have studied the fluorescence properties and stabilities of the four variants and compared them to wild-type RNase Sa and the other ribonucleases on which they were based. Our results should help others in selecting sites for adding Trp residues to proteins. The most interesting findings are: 1), Y52W is 2.9 kcal/mol less stable than RNase Sa and the fluorescence intensity emission maximum is blue-shifted to 309 nm. Only a Trp in azurin is blue-shifted to a greater extent (308 nm). This blue shift is considerably greater than observed for Trp71 in barnase, the Trp on which Y52W is based. 2), Y55W is 2.1 kcal/mol less stable than RNase Sa and the tryptophan fluorescence is almost completely quenched. In contrast, Trp59 in RNase T1, on which Y55W is based, has a 10-fold greater fluorescence emission intensity. 3), T76W is 0.7 kcal/mol more stable than RNase Sa, indicating that the Trp side chain has more favorable interactions with the protein than the threonine side chain. The fluorescence properties of folded Y76W are similar to those of the unfolded protein, showing that the tryptophan side chain in the folded protein is largely exposed to solvent. This is confirmed by the crystal structure of the T76W which shows that the side chain of the Trp is only approximately 7% buried. 4), Y81W is 0.4 kcal/mol less stable than RNase Sa. Based on the crystal structure of Y81W, the side chain of the Trp is 87% buried. Although all of the Trp side chains in the variants contribute to the unusual positive circular dichroism band observed near 235 nm for RNase Sa, the contribution is greatest for Y81W.  相似文献   

16.
A sequence alignment of mammalian cytochromes c with yeast iso-1-cytochrome c (y-cyt-c) shows that the yeast protein contains five extra N-terminal residues. We have been interested in understanding the question: What is the role of these five extra N-terminal residues in folding and stability of the protein? To answer this question we have prepared five deletants of y-cyt-c by sequentially removing these extra residues. During our studies on the wild type (WT) protein and its deletants, we observed that the amount of secondary structure in the guanidinium chloride (GdmCl)-induced denatured (D) state of each protein is different from that of the heat-induced denatured (H) state. This finding is confirmed by the observation of an additional cooperative transition curve of optical properties between H and D states on the addition of different concentrations of GdmCl to the already heat denatured WT y-cyt-c and its deletants at pH 6.0 and 68°C. For each protein, analysis of transition curves representing processes, native (N) state ? D state, N state ? H state, and H state ? D state, was done to obtain Gibbs free energy changes associated with all the three processes. This analysis showed that, for each protein, thermodynamic cycle accommodates Gibbs free energies associated with transitions between N and D states, N and H states, and H and D states, the characteristics required for a thermodynamic function. All these experimental observations have been supported by our molecular dynamics simulation studies.  相似文献   

17.
Are turns required for the folding of ribonuclease T1?   总被引:6,自引:5,他引:1       下载免费PDF全文
Ribonuclease T1 (RNase T1) is a small, globular protein of 104 amino acids for which extensive thermodynamic and structural information is known. To assess the specific influence of variations in amino acid sequence on the mechanism for protein folding, circularly permuted variants of RNase T1 were constructed and characterized in terms of catalytic activity and thermodynamic stability. The disulfide bond connecting Cys-2 and Cys-10 was removed by mutation of these residues to alanine (C2, 10A) to avoid potential steric problems imposed by the circular permutations. The original amino-terminus and carboxyl-terminus of the mutant (C2, 10A) were subsequently joined with a tripeptide linker to accommodate a reverse turn and new termini were introduced throughout the primary sequence in regions of solvent-exposed loops at Ser-35 (cp35S1), Asp-49 (cp49D1), Gly-70 (cp70G1), and Ser-96 (cp96S1). These circularly permuted RNase T1 mutants retained 35-100% of the original catalytic activity for the hydrolysis of guanylyl(3'-->5')cytidine, suggesting that the overall tertiary fold of these mutants is very similar to that of wild-type protein. Chemical denaturation curves indicated thermodynamic stabilities at pH 5.0 of 5.7, 2.9, 2.6, and 4.6 kcal/mol for cp35S1, cp49D1, cp70G1, and cp96S1, respectively, compared to a value of 10.1 kcal/mol for wild-type RNase T1 and 6.4 kcal/mol for (C2, 10A) T1. A fifth set of circularly permuted variants was attempted with new termini positioned in a tight beta-turn between Glu-82 and Gln-85. New termini were inserted at Asn-83 (cp83N1), Asn-84 (cp84N1), and Gln-85 (cp85Q1). No detectable amount of protein was ever produced for any of the mutations in this region, suggesting that this turn may be critical for the proper folding and/or thermodynamic stability of RNase T1.  相似文献   

18.
The introduction of non‐natural modules could provide unprecedented control over folding/unfolding behavior, conformational stability, and biological function of proteins. Success requires the interrogation of candidate modules in natural contexts. Here, expressed protein ligation is used to replace a reverse turn in bovine pancreatic ribonuclease (RNase A) with a synthetic β‐dipeptide: β2‐homoalanine–β3‐homoalanine. This segment is known to adopt an unnatural reverse‐turn conformation that contains a 10‐membered ring hydrogen bond, but one with a donor–acceptor pattern opposite to that in the 10‐membered rings of natural reverse turns. The RNase A variant has intact enzymatic activity, but unfolds more quickly and has diminished conformational stability relative to native RNase A. These data indicate that hydrogen‐bonding pattern merits careful consideration in the selection of beneficial reverse‐turn surrogates.  相似文献   

19.
We have examined the chemical denaturations of the Klentaq and Klenow large-fragment domains of the Type 1 DNA polymerases from Thermus aquaticus (Klentaq) and Escherichia coli (Klenow) under identical solution conditions in order to directly compare the stabilization energetics of the two proteins. The high temperature stability of Taq DNA polymerase is common knowledge, and is the basis of its use in the polymerase chain reaction. This study, however, is aimed at understanding the thermodynamic basis for this high-temperature stability. Chemical denaturations with guanidine hydrochloride report a folding free energy (DeltaG) for Klentaq that is over 20 kcal/mol more favorable than that for Klenow under the conditions examined. This difference between the stabilization free energies of a homologous mesophilic-thermophilic protein pair is significantly larger than generally observed. This is due in part to the fact that the stabilization free energy for Klentaq polymerase, at 27.5 kcal/mol, is one of the largest ever determined for a monomeric protein. Large differences in the chemical midpoints of the unfolding (Cm) and the dependences of the unfolding free energy on denaturant concentration in the transition region (m-value) between the two proteins are also observed. Measurements of the sedimentation coefficients of the two proteins in the native and denatured states report that both proteins approximately double in hydrodynamic size upon denaturation, but that Klentaq expands somewhat more than Klenow.  相似文献   

20.
We used site-specific mutagenesis by targeting E179 and F190 on the structure of photoprotein Mnemiopsin 2 (Mn2) from Mnemiopsis leidyi. The tertiary structure of E179S and F190L mutants was made by the MODELLER program. Far-ultraviolet circular dichroism data showed that the overall secondary structural content of photoprotein is not changed upon mutation, however the helicity and stabilizing interactions in helical structure decreases in mutants as compared with the wild-type (WT) photoprotein. Fluorescence spectra data revealed that the tertiary structure of the mutants is more compact than that of WT Mn2. According to the heat-induced denaturation experiments data, the melting temperature (Tm) for the unfolding of tertiary structure of the F190L variant increases by 3°C compared with that of the WT and E179S mutant. Interestingly, the conformational enthalpy of the F190L mutant (86 kcal mol−1) is considerably lower than those in the WT photoprotein (102 kcal mol−1) and E179S mutant (106 kcal mol−1). The significant difference in the enthalpy of the thermal unfolding process could be explained by considering that the thermally denatured state of the F190L mutant is structurally less expanded than the WT and E179S variants. Bioluminescence activity data showed that the maximum characteristic wavelengths of the mutants undergo blue shift as compared with the WT protein. Initial intensity of the F190L and E179S variants was recorded to be 137.5% and 55.9% of the WT protein, respectively.  相似文献   

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

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