首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 281 毫秒
1.
With its simple glycotripeptide repeating structure the antifreeze glycoprotein of polar fish may be an especially simple conformational mode for mucin glycoproteins with similar but more complex structures. The fully assigned proton n.m.r. spectrum confirms the anomeric configurations of the hexapyranosidic sugars of the side chains and the coupling constants of the alpha GalNAc and the beta Gal residues show both to be in the expected 4C1 chair conformation. The assignment of a single resonance for each proton of the (Ala-Thr-Ala)n repeat unit coupled with the observation of long range nuclear Overhauser effects (n.O.e.) implies a three-fold repeating conformation. The resonances of the two alanines are distinct and can be assigned to their correct positions in the peptide sequence by n.O.e. observed at the amide proton resonances on saturation of the alpha proton signals. The amide proton coupling constants of all three peptide residues are similar and imply a limited range of peptide backbone torsion angles, phi CN. The large n.O.e. which has been observed between the amide proton and the alpha proton of the residue preceding it in the sequence implies large positive values for the peptide dihedral angle, psi CC. Limits are placed on possible values of side chain dihedral angles by the observation of the coupling constant between the alpha and beta protons of the threonyl residue. The observation of n.O.e. between the anomeric proton of GalNAc and the threonyl side chain protons gives information on the conformation of the alpha glycosidic linkage between the disaccharide and the peptide. n.O.e. observed between the protons of the beta glycosidic linkage indicates the conformation of the disaccharide and the large amide proton coupling constant of the GalNAc residue shows a trans proton relationship. The spectroscopically derived data have been combined with conformational energy calculations to give a conformational model for antifreeze glycoprotein in which the hydrophobic surfaces of the disaccharide side chains are wrapped closely against a three-fold left handed helical peptide backbone. The hydrophilic sides of the disaccharides are aligned so that they may bind to the ice crystal face, which is perpendicular to the fast growth axis inhibiting normal crystal growth.  相似文献   

2.
Abstract

The seven α-helical segments of Bacteriorhodopsin (BR) passing through the membrane are investigated for a continuous Hydrogen Bonded Chain (HBC). The study is carried out by computer modelling approach. It is assumed that the seven helices are placed as (AGFEDCB), which has been accepted as the best model by several groups. Helices A, D, E and G are considered to be present in right handed α-helical conformation. The inter-orientation of these helices are represented by Eulerian angles α, β and γ. For the helices B, C and F which contain Proline in the middle, several conformational possibilities were considered. In these cases apart from the Eulerian angles α, β and γ, the dihedral angles φp_1 and ψp_1 of the residues that are succeeded by Proline residue in the helical regions were also used in fixing the position of the helices with respect to each other. All these parameters were varied to fit with the top, middle and bottom distances reported by electron diffraction studies. Good fit was obtained for all right handed α-helical conformations and also for helices B, C and F with a left handed turn at the residue preceeding proline. Hence two structures were analysed for continuous HBC, Structure I which contained all the seven helices in right handed α-helical conformation and Structure II, which had the helices A, D, E and G in right handed conformation and the helices B, C and F in right handed α-helical conformation with a left handed turn at the residue preceeding proline. All possible staggered conformations were considered for the side chains and the inter atomic distances were analysed for Hydrogen bonds. It was possible to obtain a continuous chain in both the structures which includes most of the residues found to be important by the experiments. However Lys-216 has to be considered in two different conformations to connect the cytoplasmic side with the extra cellular side. The overall height spanned by HBC is about 25Å. The chains obtained by both the structures I and II are analysed in terms of the conformational parameters. It has also been possible to place the retinal in the region as predicted by the experiments. The Tryptophan residues which affect the spectral characterestics can be aligned on either side of the retinal.  相似文献   

3.
The accuracy and reliability of the recently proposed scaling-relaxation method for loop closure were examined by using extensive conformational sampling. For each of the eight heptapeptides chosen to represent a variety of protein conformations, 1,000–2,000 conformations were sampled. Each segment contained 14 rotatable backbone dihedral angles. The average root mean square deviations (RMSDs) between the predicted and the native conformations were 0.7 Å for the backbone and 1.2 Å for the side chain atoms. These predictions were substantially more accurate than the previous predictions (1.1 Å for the backbone and 2.2 Å for the side chain atoms) of the same eight protein segments based on limited conformational sampling (100 conformations for each segment). Large prediction errors mostly occurred at polar and surface side chains that are unlikely to have any meaningful conformation. Moreover, the reliability of seven of the eight predictions was demonstrated with their energy-RMSD and stability-RMSD correlations of the low-energy conformations, where the conformational stability was estimated by using the multiple copy simultaneous sampling method.  相似文献   

4.
Conformational changes upon protein-protein association are the key element of the binding mechanism. The study presents a systematic large-scale analysis of such conformational changes in the side chains. The results indicate that short and long side chains have different propensities for the conformational changes. Long side chains with three or more dihedral angles are often subject to large conformational transition. Shorter residues with one or two dihedral angles typically undergo local conformational changes not leading to a conformational transition. A relationship between the local readjustments and the equilibrium fluctuations of a side chain around its unbound conformation is suggested. Most of the side chains undergo larger changes in the dihedral angle most distant from the backbone. The frequencies of the core-to-surface interface transitions of six nonpolar residues and Tyr are larger than the frequencies of the opposite surface-to-core transitions. The binding increases both polar and nonpolar interface areas. However, the increase of the nonpolar area is larger for all considered classes of protein complexes, suggesting that the protein association perturbs the unbound interfaces to increase the hydrophobic contribution to the binding free energy. To test modeling approaches to side-chain flexibility in protein docking, conformational changes in the X-ray set were compared with those in the docking decoy sets. The results lead to a better understanding of the conformational changes in proteins and suggest directions for efficient conformational sampling in docking protocols.  相似文献   

5.
Proline-induced constraints in alpha-helices   总被引:9,自引:0,他引:9  
L Piela  G Némethy  H A Scheraga 《Biopolymers》1987,26(9):1587-1600
The disrupting effect of a prolyl residue on an α-helix has been analyzed by means of conformational energy computations. In the preferred, nearly α-helical conformations of Ac-Ala4-Pro-NHMe and of Ac-Ala7-Pro-Ala7-NHMe, only the residue preceding Pro is not α-helical, while all other residues can occur in the α-helical A conformation; i.e., it is sufficient to introduce a conformational change of only one residue in order to accommodate proline in a distorted α-helix. Other low-energy conformations exist in which the conformational state of three residues preceding proline is altered considerably; on the other hand, another conformation in which these three residues retain the near-α-helical A-conformational state (with up to 26° changes of their dihedral angles ? and ψ, and a 48° change in one ω from those of the ideal α-helix) has a considerably higher energy. These conclusions are not altered by the substitution of other residues in the place of the Ala preceding Pro. The conformations of the peptide chain next to prolyl residues in or near an α-helix have been analyzed in 58 proteins of known structure, based on published atomic coordinates. Of 331 α-helices, 61 have a Pro at or next to their N-terminus, 21 have a Pro next to their C-terminus, and 30 contain a Pro inside the helix. Of the latter, 16 correspond to a break in the helix, 9 are located inside distorted first turns of the helix, and 5 are parts of irregular helices. Thus, the reported occurrence of prolyl residues next to or inside observed α-helices in proteins is consistent with the computed steric and energetic requirements of prolyl peptides.  相似文献   

6.
β‐Peptides are analogs of natural α‐peptides and form a variety of remarkably stable structures. Having an additional carbon atom in the backbone of each residue, their folded conformation is not only influenced by the side‐chain sequence but also and foremost by their substitution pattern. The precise mechanism by which the side chains interact with the backbone is, however, hitherto not completely known. To unravel the various effects by which the side chains influence the backbone conformation, we quantify to which extent the dihedral angles of a β3‐substited peptide with an additional methyl group on the central Cα‐atom can be regarded as independent degrees of freedom and analyze the distributions of these dihedral angles. We also selectively capture the steric effect of substituents on the Cα‐ and Cβ‐atoms of the central residue by alchemically changing them into dummy atoms, which have no nonbonded interactions. We find that the folded state of the β3‐peptide is primarily stabilized by a steric exclusion of large parts of the unfolded state (entropic effect) and only subsequently by mutual dependence of the ψ‐dihedral angles (enthalpic effect). The folded state of β‐peptides is stabilized by a different mechanism than that of α‐peptides. Proteins 2010. © 2010 Wiley‐Liss, Inc.  相似文献   

7.
In this study, the transition process of the helix‐sense inversion of poly(β‐phenethyl‐L‐aspartate) was investigated by Raman scattering and 2‐dimensional correlation spectroscopy. Temperature‐dependent Raman spectra were obtained during the helix‐sense inversion. The results of 2‐dimensional correlation analysis in the spectral regions of 1600‐1800 and 3200‐3400 cm?1 showed that the intensity changes of the side‐chain ester C═O stretching bands occurred prior to those of amide A and amide I bands in the unwinding process of αR‐helix on heating. The sequential order of the intensity changes for amide A, amide I, and the side‐chain ester C═O stretching bands during the inversion process was determined. It was found that the conformation change of the side chain occurred prior to that of the main chain for the αR‐helix on heating. Thus, we concluded that the transformation of the backbone chain from right‐handed to left‐handed is triggered by the conformational change of the side chains.  相似文献   

8.
The conformational preferences of azaphenylalanine-containing peptide were investigated using a model compound, Ac-azaPhe-NHMe with ab initio method at the HF/3-21G and HF/6-31G(*) levels, and the seven minimum energy conformations with trans orientation of acetyl group and the 4 minimum energy conformations with cis orientation of acetyl group were found at the HF/6-31G(*) level if their mirror images were not considered. An average backbone dihedral angle of the 11 minimum energy conformations is phi=+/-91 degrees +/-24 degrees , psi =+/-18 degrees +/-10 degrees (or +/-169 degrees +/-8 degrees ), corresponding to the i+2 position of beta-turn (delta(R)) or polyproline II (beta(P)) structure, respectively. The chi(1) angle in the aromatic side chain of azaPhe residue adopts preferentially between +/-60 degrees and +/-130 degrees, which reflect a steric hindrance between the N-terminal carbonyl group or the C-terminal amide group and the aromatic side chain with respect to the configuration of the acetyl group. These conformational preferences of Ac-azaPhe-NHMe predicted theoretically were compared with those of For-Phe-NHMe to characterize the structural role of azaPhe residue. Four tripeptides containing azaPhe residue, Boc-Xaa-azaPhe-Ala-OMe [Xaa=Gly(1), Ala(2), Phe(3), Asn(4)] were designed and synthesized to verify whether the backbone torsion angles of azaPhe reside are still the same as compared with theoretical conformations and how the preceding amino acids of azaPhe residue perturb the beta-turn skeleton in solution. The solution conformations of these tripeptide models containing azaPhe residue were determined in CDCl(3) and DMSO solvents using NMR and molecular modeling techniques. The characteristic NOE patterns, the temperature coefficients of amide protons and small solvent accessibility for the azapeptides 1-4 reveal to adopt the beta-turn structure. The structures of azapeptides containing azaPhe residue from a restrained molecular dynamics simulation indicated that average dihedral angles [(phi(1), psi(1)), (phi(2), psi(2))] of Xaa-azaPhe fragment in azapeptide, Boc-Xaa-azaPhe-Ala-OMe were [(-68 degrees, 135 degrees ), (116 degrees, -1 degrees )], and this implies that the intercalation of an azaPhe residue in tripeptide induces the betaII-turn conformation, and the volume change of a preceding amino acid of azaPhe residue in tripeptides would not perturb seriously the backbone dihedral angle of beta-turn conformation. We believe such information could be critical in designing useful molecules containing azaPhe residue for drug discovery and peptide engineering.  相似文献   

9.
V N Damle 《Biopolymers》1970,9(8):937-954
Absorption, circular dichroism (CD), and optical rotatory dispersion (ORD) measurements were carried out on poly-L -tyrosine in trimethyl phosphate solution over the spectral range 185–600 mμ. There is evidence in the CD spectrum for side chain-side chain interactions in poly-L -tyrosine. ORD and CD data in dimethylformamide and pyridine closely parallel those in trimethyl phosphate, indicating a similarity in conformation of the polymer in all three solvents. In the polarized infrared spectrum both position and polarization of amide A, I, and II bands are characteristic of α-helical polypeptides. Bands corresponding to side chain also exhibit dichroism, suggesting that the side chains are not randomly oriented. Viscosity and light-scattering studies are consistent with α-helical structure for the polymer that, remains rigid over a temperature range of 15–50°C and becomes somewhat flexible at higher temperatures. Optical rotatory properties were found to vary gradually and continuously with temperature over the range of ?30 to +100°C. This suggested that all three electronic transitions of tyrosyl side chain are optically active, and that the side chains have some freedom of motion that decreases with decreasing temperature, disappearing only at about ?30°C.  相似文献   

10.
Like all other complex biological systems, proteins exhibit properties not found in free amino acids (i.e., emergent properties). Here, we explore top-down constraints experienced by the residue side chains in proteins compared to amino acids in increasingly complex molecular environments: free amino acids, end-capped amino acids, and the central residue in an alpha-helical nonapeptide. The crystalline structure of the contractile protein profilin Ib and the enzyme trypsin were chosen as objects of study, and submitted to 10 ns molecular dynamics (MD) simulations. The results revealed increased conformational constraints on the side chains when going from the simpler to the more complex compounds. A Shannon entropy (SE) analysis of the conformational behavior of the side chains showed in most cases a progressive and marked decrease in the SE of the chi1 and chi2 dihedral angles. This is equivalent to stating that conformational constraints on the side chain of residues increase their information content and, hence, recognition specificity compared to free amino acids. In other words, the vastly increased information content of a protein relative to its free monomers is embedded not only in the tertiary structure of the backbone, but also in the conformational behavior of the side chains. The postulated implication is that both backbone and side chains, by virtue of being conformationally constrained, contribute to the protein's recognition specificity toward other macromolecules and ligands.  相似文献   

11.
Conformational energy computations were carried out on collagenlike triple-stranded conformations of several poly(tripeptide)s with the general structure CH3CO? (Gly? X? Y)3? NHCH3. The sequences considered had various amino acid residues in position X or Y of the central tripeptide, with either Pro or Ala as a neighbor, i.e., Gly-X-Pro, Gly-X-Ala, Gly-Pro-Y, and Gly-Ala-Y. Minimum-energy conformations were computed for the side chains, and their distributions were compared for the four sequences. The residues used were Abu (= α-aminobutyric acid), Leu, Phe, Ser, Asp, Asn, Val, Ile, and Thr. The conformational energy of a ? Ch2? CH3 side chain in Abu was mapped as a function of the dihedral angle χ1. Intrastrand interactions with neighboring residues do not affect the conformations of a side chain in position Y, and they have a minor effect on it in the X-Ala sequence, but they strongly restrict the conformational freedom of the side chain in the X-Pro sequence. Conversely, interstrand interactions do not affect side chains in position X, but they strongly restrict the conformational freedom of a side chain in position Y if there is a nearby Pro residue in a neighboring strand. Hydrogen bonds with the backbone can be formed in some conformations of long polar side chains, such as Asp, Asn, or Gln. All amino acid residues can be accommodated in collagen. Because of the interactions mentioned above, steric and energetic constraints can be correlated with observed preferences of certain amino acids for positions X or Y in collagen. Hence, these preferences may be explained, in part, in terms of differences in the conformational freedom of the side chains in the triple-stranded structure.  相似文献   

12.
13.
A method for generating a complete polypeptide backbone structure from a set of Cα coordinates is presented. Initial trial values of ? and ψ for a selected residue are chosen (essentially from an identification of the conformational region of the virtual-bond backbone, e.g., and α-helical region), and values of ? and ψ for the remaining residues (both towards the N- and C-terminus) are then computed, subject to the constraint that the chain have the same virtual-bond angles and virtual-bond dihedral angles as the given set of Cα coordinates. The conversion from Cα coordinates to full backbone dihedral angles (?,ψ) involves the solution of a set of algebraic equations relating the virtual-bond angles and virtual-bond dihedral angles to standard peptide geometry and backbone dihedral angles. The procedure has been tested successfully on Cα coordinates taken from standard-geometry full-atom structures of bovine pancreatic trypsin inhibitor (BPTI). Some difficulty was encountered with error-sensitive residues, but on the whole the backbone generation was successful. Application of the method to Cα coordinates for BPTI derived from simplified model calculations (involving nonstandard geometry) showed that such coordinates may be inconsistent with the requirement that ?Pro be near ?75°. In such a case, i.e., for residues for which the algebraic method failed, a leastsquares minimizer was then used in conjunction with the algebraic method; the mean-square deviation of the calculated Cα coordinates from the given ones was minimized by varying the backbone dihedral angles. Thus, these inconsistencies were circumvented and a full backbone structure whose Cα coordinates had an rms deviation of 0.26 Å from the given set of Cα coordinates was obtained.  相似文献   

14.
The Langevin dynamics method and statistical correlation analysis were used to study the α-helical structure folding dynamics of the (Ala)50, (AlaGly)25, and (AlaGly)75 polypeptides depending on the viscosity of the medium. Friction forces that arise when the effective viscosity of the medium is similar to the viscosity of water were found to result in strong correlations between the backbone torsion angles. The polypeptides under study folded mainly to produce α-helical structures. A structure of two contacting α-helices that were approximately equal in length and had a loop between them was observed for a longer chain of 150 residues. A method to visualize the correlation matrix of the dihedral angles of a polypeptide chain was developed for analyzing the effects of the dynamic correlation of conformational degrees of freedom. The analysis of the dynamics of the correlation matrix showed that rotations involving angles of the same type (φ–φ and ψ–ψ) occur predominantly in one direction. Rotations invoving different angles (φ–ψ) occur predominantly in opposite directions, so that the total macromolecule does not rotate. A significant reduction in the effective viscosity of the medium disrupts the correlation and makes the rotations stochastic, thus distorting the formation of the regular (helical) structure. The effects of correlated conformational motions are consequences of viscous friction forces. This conclusion agrees with our previous results that outlined the principle of the minimum rate of energy dissipation and the equipartition of energy dissipation rate between conformational degrees of freedom.  相似文献   

15.
Potential energies of conformation of a dipeptide unit with butyl, seryl, threonyl, eysteinyl, and valyl side groups have been computed by using classical energy expressions. The presence of a γ-atom introduces characteristic restrictions on the backbone rotational angles ? and ψ the γ-atom itself is restricted to three staggered positions about the Cα—Cβ bond. The important results are that a γ-carbon in position I (χ1 ? 60°) cannot be accommodated in the standard right-and left-handed α-helices, whereas a γ-oxygen or sulfur could easily be accommodated in the right-handed α-helix. Further, a γ-carbon or a heteroatom in position II (χ1 ? 180°) does not favor a conformation ψ ? 180°, compared to two other positions. The valyl side group significantly reduces the allowed ? and ψ values and energetically prefers a β-conformation compared to right-or left-handed α-helical conformations. The less favorable α-helical conformation is possible only for γ (III, II) combination of the valyl residue. The observed ?, ψ, and χ1 values of all the amino acid residues in the three protein molecules, lysozyme, myoglobin, and chymotrypsin are compared with the theoretical predictions and the agreement is excellent. The results bring out the important fact that even in large molecules, the conformation of local segments are predominantly governed by the short-range intramolecular interactions.  相似文献   

16.
B Mao 《Biophysical journal》1991,60(3):611-622
Atomic motions in protein molecules have been studied by molecular dynamics (MD) simulations; dynamics simulation methods have also been employed in conformational studies of polypeptide molecules. It was found that when atomic masses are weighted, the molecular dynamics method can significantly increase the sampling of dihedral conformation space in such studies, compared to a conventional MD simulation of the same total simulation time length. Herein the theoretical study of molecular conformation sampling by the molecular dynamics-based simulation method in which atomic masses are weighted is reported in detail; moreover, a numerical scheme for analyzing the extensive conformational sampling in the simulation of a tetrapeptide amide molecule is presented. From numerical analyses of the mass-weighted molecular dynamics trajectories of backbone dihedral angles, low-resolution structures covering the entire backbone dihedral conformation space of the molecule were determined, and the distribution of rotationally stable conformations in this space were analyzed quantitatively. The theoretical analyses based on the computer simulation and numerical analytical methods suggest that distinctive regimes in the conformational space of the peptide molecule can be identified.  相似文献   

17.
Backbone mimicry by the formation of closed-loop C7, C10 and C13 (mimics of gamma-, beta- and alpha-turns) conformations through side chain-main chain hydrogen bonds by polar groups is a frequent observation in protein structures. A data set of 250 non-homologous and high-resolution protein crystal structures was used to analyze these conformations for their characteristic features. Seven out of the nine polar residues (Ser, Thr, Asn, Asp, Gln, Glu and His) have hydrogen bonding groups in their side chains which can participate in such mimicry and as many as 15% of all these polar residues engage in such conformations. The distributions of dihedral angles of these mimics indicate that only certain combinations of the dihedral angles involved aid the formation of these mimics. The observed examples were categorized into various classes based on these combinations, resulting in well defined motifs. Asn and Asp residues show a very high capability to perform such backbone secondary structural mimicry. The most highly mimicked backbone structure is of the C10 conformation by the Asx residues. The mimics formed by His, Ser, Thr and Glx residues are also discussed. The role of such conformations in initiating the formation of regular secondary structures during the course of protein folding seems significant.  相似文献   

18.
In order to search for probable conformations of the peptide, the amino acid side chain, and the carbohydrate linkage in glycoproteins, conformational energy surfaces of glycopeptide model compounds were studied by Monte Carlo methods using the Metropolis algorithm. The potential energies were composed of empirical energy functions which include nonbonded interactions, electrostatics, hydrogen bonding, and torsional energies specified by parameters which have been used for peptides. Calculations were performed on 1-N-acetyl-2-acetamido-beta-D-glucopyranosyl amine and the glycosylated dipeptide N-acetyl-delta-N-(2-acetamido-beta-D-glucopyranosyl)-L-asparaginyl-N'-methyl amide as models for N-glycosylated peptides and on methyl-2-acetamido-alpha-D-galactopyranoside as well as the glycosylated dipeptides N-acetyl-gamma-O-(2-acetamido-alpha-D-galactopyranosyl)-L-threonyl-N'-methyl amide and its seryl analog as models for O-glycosylated glycoproteins. The probable conformations of these compounds were analyzed by single-angle probability tables and by two-dimensional conformation density maps projected from the Markov chains which contained up to six independently varied conformational dihedral angles. The presence of high barriers to rotation required the use of search strategies which resulted in a rather low acceptance rate for new conformations in the Metropolis algorithm in order to avoid trapping of the Markov chain in local energy minima. This problem contributed to the failure of these calculations to attain complete convergence to the thermodynamic limit for the glycosylated dipeptide models in which six dihedral angles were independently varied. Analysis of the results shows that the conformational space available to the highly crowded axial glycosides of the alpha-O-GalNAc type is much more restricted than that for the N-asparaginyl glycopeptides. The most probable conformation for the O-glycosylated peptides is is a beta-turn while N-glycosylated peptides may be either in a beta-turn or an extended conformation.  相似文献   

19.
B N Rao  V K Dua  C A Bush 《Biopolymers》1985,24(12):2207-2229
Spectroscopic data and conformational energy calculations are reported for eight oligosaccharides from ovarian cyst mucins and from human milk, the nonreducing terminals of which have fucose (α1 → 2)galactose linked either (β1 → 3) (type I) or (β1 → 4)(type II) to N-acetylglucosamine or in (β1 → 3) linkage to galactosaminitol. The fully assigned proton nmr spectra are reported along with nuclear Overhauser enhancement (NOE) data. Amide proton coupling constants and vacuum-uv CD spectra provide information on the amide plane orientation and amide environment. Our results imply that the fucosidic dihedral angles are similar for all three cases and that the substantial differences in the chemical shifts of the fucosyl protons of type I, type II, and 3-ol chains result from different perturbations by the amide group of the residue to which the β-galactose is linked. Stereopair diagrams of conformational models for both type I and II H chains are presented that are consistent with NOE, coupling constants, conformational energy calculations, and the CD data. While the temperature dependence of the observed NOE of penta- and hexasaccharides indicates that their rotational correlation times are strongly temperature dependent, we conclude that the conformations are essentially independent of temperature.  相似文献   

20.
Broadening of the infrared amide A, amide I and amide II bands of α-helical polypeptides has been observed for thermodynamically unstable α-helices. This spectroscopic fact can be explained now by the geometrical distortions of the backbone of the helical structure. Two models for distorted helices which include regular or irregular distortions of the angles of internal rotation of the main polypeptide chain have been considered. It is pointed out that the instability of α-helix is associated with irregular distortions of the polypeptide backbone.  相似文献   

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

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