首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
An automatic procedure for building a protein polyalanine backbone from C alpha positions and 'spare parts' retrieved from a data base of 66 high-resolution protein structures is described. Protein backbones are constructed from overlapping fragments of variable length, which allows the backbone of regular secondary structure elements to be built in one block. The procedure is shown to yield backbones which compare very favourably with those from highly refined X-ray structures (r.m.s. deviation between generated and crystal structures less than 1A). The method is furthermore quite insensitive to experimental errors in C alpha positions as well as to the size of the data base, and is seen to yield valuable insight into the relationships between sequence and 3-D structure: one example on triose phosphate isomerase, a beta-barrel protein, shows that beta alpha loops can be considered as structurally more uncommon than alpha beta loops. The 'spare parts' approach is also found to be useful for general-purpose modelling of local structural changes produced by insertion or deletion of residues. It should, however, be used with caution. Crude selection criteria based solely on fragment length and geometric fit to the loop base regions yield realistic backbones in about two-thirds of the test cases (r.m.s. deviations from refined crystal structure approximately 1A). In the remaining cases, sequence information, in particular the presence of glycine residues which tend to adopt more unusual backbone conformations, must be considered to obtain comparable results.  相似文献   

2.
Protein C alpha coordinates are used to accurately reconstruct complete protein backbones and side-chain directions. This work employs potentials of mean force to align semirigid peptide groups around the axes that connect successive C alpha atoms. The algorithm works well for all residue types and secondary structure classes and is stable for imprecise C alpha coordinates. Tests on known protein structures show that root mean square errors in predicted main-chain and C beta coordinates are usually less than 0.3 A. These results are significantly more accurate than can be obtained from competing approaches, such as modeling of backbone conformations from structurally homologous fragments.  相似文献   

3.
The crystal structure of beta-lactamase TEM1 from E. coli has been solved to 2.5 A resolution by X-ray diffraction methods and refined to a crystallographic R-factor of 22.7%. The structure was determined by multiple isomorphous replacement using four heavy atom derivatives. The solution from molecular replacement, using a polyalanine model constructed from the C alpha coordinates of S. Aureus PCl enzyme, provided a set of phases used for heavy atom derivatives analysis. The E. coli beta-lactamase TEM1 is made up of two domains whose topology is similar to that of the PCl enzyme. However, global superposition of the two proteins shows significant differences.  相似文献   

4.
Folding simulations of polyalanine peptides were carried out using an off-lattice Monte Carlo simulation technique. The peptide was represented as a chain of residues, each of which contains two interaction sites: one corresponding to the C(alpha) atom and the other to the side chain. A statistical potential was used to describe the interaction between these sites. The preferred conformations of the peptide chain on the energy surface, starting from several initial conditions, were searched by perturbations on its generalized coordinates with the Metropolis criterion. We observed that, at low temperatures, the effective energy was low and the helix content high. The calculated helix propagation (s) and nucleation (sigma) parameters of the Zimm-Bragg model were in reasonable agreement with the empirical data. Exploration of the energy surface of the alanine-based peptides (AAQAA)(3) and AAAAA(AAARA)(3)A demonstrated that their behavior is similar to that of polyalanine, in regard to their effective energy, helix content, and the temperature-dependence of their helicity. In contrast, stable secondary structures were not observed for (Gly)(20) at similar temperatures, which is consistent with the nonfolder nature of this peptide. The fluctuations in the slowest dynamics mode, which describe the elastic behavior of the chain, showed that as the temperature decreases, the polyalanine peptides become stiffer and retain conformations with higher helix content. Clustering of conformations during the folding phase implied that polyalanine folds into a helix through fewer numbers of intermediate conformations as the temperature decreases.  相似文献   

5.
The problem of constructing all-atom model co-ordinates of a protein from an outline of the polypeptide chain is encountered in protein structure determination by crystallography or nuclear magnetic resonance spectroscopy, in model building by homology and in protein design. Here, we present an automatic procedure for generating full protein co-ordinates (backbone and, optionally, side-chains) given the C alpha trace and amino acid sequence. To construct backbones, a protein structure database is first scanned for fragments that locally fit the chain trace according to distance criteria. A best path algorithm then sifts through these segments and selects an optimal path with minimal mismatch at fragment joints. In blind tests, using fully known protein structures, backbones (C alpha, C, N, O) can be reconstructed with a reliability of 0.4 to 0.6 A root-mean-square position deviation and not more than 0 to 5% peptide flips. This accuracy is sufficient to identify possible errors in protein co-ordinate sets. To construct full co-ordinates, side-chains are added from a library of frequently occurring rotamers using a simple and fast Monte Carlo procedure with simulated annealing. In tests on X-ray structures determined at better than 2.5 A resolution, the positions of side-chain atoms in the protein core (less than 20% relative accessibility) have an accuracy of 1.6 A (r.m.s. deviation) and 70% of chi 1 angles are within 30 degrees of the X-ray structure. The computer program MaxSprout is available on request.  相似文献   

6.
The building of protein structures from alpha-carbon coordinates   总被引:3,自引:0,他引:3  
P E Correa 《Proteins》1990,7(4):366-377
A procedure for the construction of complete protein structures from only alpha-carbon coordinates is described. This involves building the backbone by sequential addition of Pro, Gly, or Ala residues. This main chain structure is then refined using molecular dynamics. Side chains are constructed by sequential addition of atoms with intermediate molecular dynamics refinement. For alpha lytic protease (a structure that is mostly beta sheet) a backbone root mean square deviation (RMSD) of 0.19 A and an overall RMSD of 1.24 A from the crystallographic coordinates are attained. For troponin C (67% alpha-helix), where the coordinates are available only for the alpha-carbons, a backbone RMSD of 0.41 A and an overall RMSD of 1.68 A are attained (fits kindly provided by Dr. Michael James and Natalie Strynadka). For flavodoxin a backbone RMSD of 0.49 A and an overall RMSD of 1.64 A were attained.  相似文献   

7.
We present an approach to the analysis of low-frequency (0-200 cm?1) α-helix vibrations in molecular dynamics simulation. The approach employs the P-Curves algorithm [H. Sklenar, C. Etchebest, and R. Lavery, (1989) Proteins: Structure, Function and Genetics, Vol. 6. pp. 46–60] to determine the helical axis and a set of helicoidal parameters describing the axis curvature and the position of the repealing units with respect to the axis and each other. The vibrations are analyzed in terms of time correlation functions of the fluctuations of P-Curves parameters and their Fourier transforms. Simulations of polyalanine and myoglobin are analyzed. For polyalanine, global twisting, bending, and stretching vibrations are found at 11, 20, and 40 cm?1, respectively. In myoglobin, the spectra of the global helix vibrations are qualitatively different from those of polyalanine and considerably more complicated. Local vibrations of individual amino acid units in the helix backbones are also analyzed with P-Curves and compared. © 1995 John Wiley & Sons, Inc.  相似文献   

8.
There are many kinds of silks from silkworms and spiders with different structures and properties, and thus, silks are suitable to study the structure-property relationship of fibrous proteins. Silk fibroin from a wild silkworm, Samia cynthia ricini, mainly consists of the repeated similar sequences by about 100 times where there are alternative appearances of the polyalanine (Ala)(12-13) region and the Gly-rich region. In this paper, a sequential model peptide, GGAGGGYGGDGG(A)(12)GGAGDGYGAG, which is a typical sequence of the silk fibroin, was synthesized, and the atomic-level conformations of Gly residues at the N- and C-terminal ends of the polyalanine region were determined as well as that of the central Ala residue using (13)C 2D spin diffusion solid-state nuclear magnetic resonance (NMR) under off-magic angle spinning. In the model peptide with alpha-helical conformation, the torsion angle of the central Ala residue, the 19th Ala, was determined to be (phi, psi) = (-60 degrees, -50 degrees ), which was a typical alpha-helical structure, but the torsion angles of two Gly residues, the 12th and 25th Gly residues, which are located at the N- and C-terminal ends of the polyalanine region, were determined to be (phi, psi) = (-70 degrees, -30 degrees ) and (phi, psi) = (-70 degrees, -20 degrees ), respectively. Thus, it was observed that the turns at both ends of polyalanine with alpha-helix conformation in the model peptide are tightly wound.  相似文献   

9.
Here we perform a systematic exploration of the use of distance constraints derived from small angle X-ray scattering (SAXS) measurements to filter candidate protein structures for the purpose of protein structure prediction. This is an intrinsically more complex task than that of applying distance constraints derived from NMR data where the identity of the pair of amino acid residues subject to a given distance constraint is known. SAXS, on the other hand, yields a histogram of pair distances (pair distribution function), but the identities of the pairs contributing to a given bin of the histogram are not known. Our study is based on an extension of the Levitt-Hinds coarse grained approach to ab initio protein structure prediction to generate a candidate set of C(alpha) backbones. In spite of the lack of specific residue information inherent in the SAXS data, our study shows that the implementation of a SAXS filter is capable of effectively purifying the set of native structure candidates and thus provides a substantial improvement in the reliability of protein structure prediction. We test the quality of our predicted C(alpha) backbones by doing structural homology searches against the Dali domain library, and find that the results are very encouraging. In spite of the lack of local structural details and limited modeling accuracy at the C(alpha) backbone level, we find that useful information about fold classification can be extracted from this procedure. This approach thus provides a way to use a SAXS data based structure prediction algorithm to generate potential structural homologies in cases where lack of sequence homology prevents identification of candidate folds for a given protein. Thus our approach has the potential to help in determination of the biological function of a protein based on structural homology instead of sequence homology.  相似文献   

10.
Dihedral probability grid Monte Carlo (DPG-MC) is a general-purpose method of conformational sampling that can be applied to many problems in peptide and protein modeling. Here we present the DPG-MC method and apply it to predicting complete protein structures from C alpha coordinates. This is useful in such endeavors as homology modeling, protein structure prediction from lattice simulations, or fitting protein structures to X-ray crystallographic data. It also serves as an example of how DPG-MC can be applied to systems with geometric constraints. The conformational propensities for individual residues are used to guide conformational searches as the protein is built from the amino-terminus to the carboxyl-terminus. Results for a number of proteins show that both the backbone and side chain can be accurately modeled using DPG-MC. Backbone atoms are generally predicted with RMS errors of about 0.5 A (compared to X-ray crystal structure coordinates) and all atoms are predicted to an RMS error of 1.7 A or better.  相似文献   

11.
We propose a coarse-grained (CG) model to study the native structure and physical properties of helical membrane proteins (HMPs) using off-lattice computer simulations. Instead of considering sequence heterogeneity explicitly, we model its effect on the packing of helices by employing a mean packing parameter r(0), which is calculated from an all-atom (AA) model. Specifically, this CG model is applied to investigate the packing of helices in bacteriorhodopsin (BR), and predicts the seven helix bundle structure of BR with a root mean square deviation (RMSD) in coordinates of helix backbone atoms (N, C, C(alpha)) of 3.99 A from its crystal structure. This predicted structure is further refined in an AA model by Amber and the refined structure has a RMSD (in coordinates of helix backbone atoms) of 2.64 A. The predicted packing position, tilting angle, and orientation angle of each helix in the refined structure are consistent with experimental data and their physical origins can be well understood in our model. Our results show that a reasonably good structure of BR can be predicted by using such a dual-scale approach, provided that its secondary structure is known. Starting from a random initial configuration, the folded structure can be obtained in days using a regular desktop computer. Various thermodynamic properties of helix packing of BR are also investigated in this CG model.  相似文献   

12.
Bacterial aromatic polyketides such as tetracycline and doxorubicin are a medicinally important class of natural products produced as secondary metabolites by actinomyces bacteria. Their backbones are derived from malonyl-CoA units by polyketide synthases (PKSs). The nascent polyketide chain is synthesized by the minimal PKS, a module consisting of four dissociated enzymes. Although the biosynthesis of most aromatic polyketide backbones is initiated through decarboxylation of a malonyl building block (which results in an acetate group), some polyketides, such as the estrogen receptor antagonist R1128, are derived from nonacetate primers. Understanding the mechanism of nonacetate priming can lead to biosynthesis of novel polyketides that have improved pharmacological properties. Recent biochemical analysis has shown that nonacetate priming is the result of stepwise activity of two dissociated PKS modules with orthogonal molecular recognition features. In these PKSs, an initiation module that synthesizes a starter unit is present in addition to the minimal PKS module. Here we describe a general method for the engineered biosynthesis of regioselectively modified aromatic polyketides. When coexpressed with the R1128 initiation module, the actinorhodin minimal PKS produced novel hexaketides with propionyl and isobutyryl primer units. Analogous octaketides could be synthesized by combining the tetracenomycin minimal PKS with the R1128 initiation module. Tailoring enzymes such as ketoreductases and cyclases were able to process the unnatural polyketides efficiently. Based upon these findings, hybrid PKSs were engineered to synthesize new anthraquinone antibiotics with predictable functional group modifications. Our results demonstrate that (i) bimodular aromatic PKSs present a general mechanism for priming aromatic polyketide backbones with nonacetate precursors; (ii) the minimal PKS controls polyketide chain length by counting the number of atoms incorporated into the backbone rather than the number of elongation cycles; and (iii) in contrast, auxiliary PKS enzymes such as ketoreductases, aromatases, and cyclases recognize specific functional groups in the backbone rather than overall chain length. Among the anthracyclines engineered in this study were compounds with (i) more superior activity than R1128 against the breast cancer cell line MCF-7 and (ii) inhibitory activity against glucose-6-phosphate translocase, an attractive target for the treatment of Type II diabetes.  相似文献   

13.
The interaction of purified alpha alpha cross-linked hemoglobin (alpha alpha Hb) with a pentaacylated mutant lipopolysaccharide (pLPS) and the corresponding lipid A (pLA) was studied biophysically and the effects correlated with data from biological assays, i.e. cytokine induction (tumor necrosis factor-alpha) in human mononuclear cells and the Limulus amebocyte lysate assay. Fourier transform infrared spectroscopic and Zeta-Sizer experiments indicated an electrostatic as well as a non-electrostatic binding of alpha alpha Hb to the hydrophilic and to the hydrophobic moieties of the endotoxins with an increase of the inclination angle of the pLA backbone, with respect to the membrane surface, from 25 degrees to more than 50 degrees. Small angle synchrotron radiation x-ray diffraction measurements indicated a reorientation of the lipid A aggregates from a multilamellar into a cubic structure as a result of alpha alpha Hb interaction. Thus, in the absence of alpha alpha Hb, the molecular shape of the pentaacyl samples was cylindrical with a moderate inclination of the diglucosamine backbone, whereas, in the presence of the protein, the shape was conical, and the inclination angle was high. The cytokine-inducing capability in human mononuclear cells, negligible for the pure pentaacylated compounds, increased markedly in the presence of alpha alpha Hb in a concentration-dependent manner. In the Limulus assay, the pentaacylated samples were active a priori, and their activity was enhanced following binding to alpha alphaHb, at least at the highest protein concentrations. The data can be understood in the light of a reaggregation of the endotoxins because of alpha alpha Hb binding, with the endotoxin backbones then readily accessible for serum and membrane proteins. By using fluorescence resonance energy transfer spectroscopy, an uptake of the endotoxin-Hb complex into phospholipid liposomes was observed, which provides a basis for cell activation.  相似文献   

14.
Musingarimi P  Plumb ME  Sodetz JM 《Biochemistry》2002,41(37):11255-11260
Human C8 is one of five complement components (C5b, C6, C7, C8, and C9) that interact to form the cytolytic membrane attack complex (MAC). It is an oligomeric protein composed of a disulfide-linked C8alpha-gamma heterodimer and a noncovalently associated C8beta chain. C8alpha and C8beta are homologous; both contain an N-terminal thrombospondin type 1 (TSP1) module, a low-density lipoprotein receptor class A (LDLRA) module, an extended central segment referred to as the membrane attack/perforin (MACPF) domain, an epidermal growth factor (EGF) module, and a second TSP1 module at the C-terminus. In this study, the segment of C8beta that confers binding specificity toward C8alpha-gamma was identified using recombinant C8beta constructs in which the N- and/or C-terminal modules were deleted or exchanged with those from C8alpha. Constructs were tested for their ability to bind C8alpha-gamma in solution and express C8 hemolytic activity. Binding to C8alpha-gamma was found to be dependent on the TSP1 + LDLRA + MACPF segment of C8beta. Within this segment, the TSP1 module and MACPF domain are principally involved and act cooperatively to mediate binding. Results from activity assays suggest that residues within this segment also mediate binding and incorporation of C8 into the MAC.  相似文献   

15.
Human interferon-alpha 8 (HuIFN alpha 8), a type I interferon (IFN), is a cytokine belonging to the hematopoietic super-family that includes human growth hormone (HGH). Recent data identified two human type I IFN receptor components. One component (p40) was purified from human urine by its ability to bind to immobilized type I IFN. A second receptor component (IFNAR), consisting of two cytokine receptor-like domains (D200 and D200'), was identified by expression cloning. Murine cells transfected with a gene encoding this protein were able to produce an antiviral response to human IFN alpha 8. Both of these receptor proteins have been identified as members of the immunoglobulin superfamily of which HGH receptor is a member. The cytokine receptor-like structural motifs present in p40 and IFNAR were modeled based on the HGH receptor X-ray structure. Models of the complexes of HuIFN alpha 8 with the receptor subunits were built by superpositioning the conserved C alpha backbone of the HuIFN alpha 8 and receptor subunit models with HGH and its receptor complex. The HuIFN alpha 8 model was constructed from the C alpha coordinates of murine interferon-beta crystal structure. Electrostatic potentials and hydrophobic interactions appear to favor the model of HuIFN alpha 8 interacting with p40 at site 1 and the D200' domain of IFNAR at site 2 because there are regions of complementary electrostatic potential and hydrophobic interactions at both of the proposed binding interfaces. Some of the predicted receptor binding residues within HuIFN alpha 8 correspond to functionally important residues determined previously for human IFN alpha 1, IFN alpha 2, and IFN alpha 4 subtypes by site-directed mutagenesis studies. The models predict regions of interaction between HuIFN alpha 8 and each of the receptor proteins, and provide insights into interactions between other type I IFNs (IFN-alpha subtypes and IFN-beta) and their respective receptor components.  相似文献   

16.
The laminin alpha3 chain LG4 module (alpha3LG4 module) has cell adhesion, heparin binding, migration, and neurite outgrowth activities. The LG4 module consists of a 14-stranded beta-sheet (A-N) sandwich structure. Previously, we identified the A3G756 sequence (KNSFMALYLSKGRLVFALG in the human laminin alpha3 chain 1411-1429) as a biologically active site in the alpha3LG4 module. The A3G756 sequence is located on the E and F strands based on a crystal structure-based sequence alignment. The Lys1421 and Arg1423 residues, critical amino acids for the biological activity of A3G756, are located on the E-F connecting loop region as a KGR sequence. In this study, we focused on the KGR sequence and investigated the structural requirements of the E-F connecting loop region in the alpha3LG4 module. We synthesized three linear peptides containing the KGR sequence at the middle and the N and C termini and also prepared three cyclic analogues corresponding to the linear peptides. cyclo-hEF3A (CLYLSKGRLVFAC), which is a cyclic peptide containing the KGR sequence at the middle, showed the strongest inhibitory effect on both the heparin binding and the cell attachment to the recombinant alpha3LG4 module protein. The cyclo-hEF3A peptide was more active for syndecan-4 binding and neurite outgrowth than the linear form. Furthermore, we found that the structure of cyclo-hEF3A is similar to that of the connecting E-F loop region in human laminin alpha3LG4 module by structural analysis using molecular dynamics simulations. These results suggest that the loop structure of the E-F connecting region of the alpha3LG4 module is important for its biological activities. The cyclo-hEF3A peptide may be useful for the development of therapeutic reagents especially for wound healing and nerve regeneration.  相似文献   

17.
The C alpha backbones of the glucose isomerase molecules of Streptomyces rubiginosus and Arthrobacter have been determined by X-ray crystallography and compared. Each molecule is a tetramer of eight-stranded alpha/beta barrels, and the mode of association of the tetramers is identical in each case. The Arthrobacter electron density shows four additional amino acids at the carboxyl terminus. There is also an insertion of six amino acids at position 277, and two individual insertions at about positions 348 and 357 (numbering according to the Streptomyces structure). There is a close structural homology throughout the whole molecule, which is most accurate up to position 325. The r.m.s. displacement for 315 homologous C alpha positions up to this position is 0.92 A.  相似文献   

18.
It is known that folding a protein chain into a cubic lattice is an NP-complete problem. We consider a seemingly easier problem: given a three-dimensional (3D) fold of a protein chain (coordinates of its C(alpha) atoms), we want to find the closest lattice approximation of this fold. This problem has been studied under names such as "lattice approximation of a protein chain", "the protein chain fitting problem", and "building of protein lattice models". We show that this problem is NP-complete for the cubic lattice with side close to 3.8 A and coordinate root mean square deviation.  相似文献   

19.
The global fold of maltose-binding protein in complex with the substrate beta-cyclodextrin was determined by solution NMR methods. The two-domain protein is comprised of a single polypeptide chain of 370 residues, with a molecular mass of 42 kDa. Distance information in the form of H(N)-H(N), H(N)-CH(3) and CH(3)-CH(3) NOEs was recorded on (15)N, (2)H and (15)N, (13)C, (2)H-labeled proteins with methyl protonation in Val, Leu, and Ile (C(delta1) only) residues. Distances to methyl protons, critical for the structure determination, comprised 77 % of the long-range restraints. Initial structures were calculated on the basis of 1943 NOEs, 48 hydrogen bond and 555 dihedral angle restraints. A global pair-wise backbone rmsd of 5.5 A was obtained for these initial structures with rmsd values for the N and C domains of 2.4 and 3.8 A, respectively. Direct refinement against one-bond (1)H(N)-(15)N, (13)C(alpha)-(13)CO, (15)N-(13)CO, two-bond (1)H(N)-(13)CO and three-bond (1)H(N)-(13)C(alpha) dipolar couplings resulted in structures with large numbers of dipolar restraint violations. As an alternative to direct refinement against measured dipolar couplings we have developed an approach where discrete orientations are calculated for each peptide plane on the basis of the dipolar couplings described above. The orientation which best matches that in initial NMR structures calculated from NOE and dihedral angle restraints exclusively is used to refine further the structures using a new module written for CNS. Modeling studies from four different proteins with diverse structural motifs establishes the utility of the methodology. When applied to experimental data recorded on MBP the precision of the family of structures generated improves from 5.5 to 2.2 A, while the rmsd with respect to the X-ray structure (1dmb) is reduced from 5.1 to 3.3 A.  相似文献   

20.
Wang Z  Plaxco KW  Makarov DE 《Biopolymers》2007,86(4):321-328
Although recent spectroscopic studies of chemically denatured proteins hint at significant nonrandom residual structure, the results of extensive small angle X-ray scattering studies suggest random coil behavior, calling for a coherent understanding of these seemingly contradicting observations. Here, we report the results of a Monte Carlo study of the effects of two types of local structures, alpha helix and Polyproline II (PPII) helix, on the dimensions of random coil polyalanine chains viewed as a model of highly denatured proteins. We find that although Flory's power law scaling, long regarded as a signature of random coil behavior, holds for chains containing up to 90% alpha or PPII helix, the absolute magnitude of the chain dimensions is sensitive to helix content. As residual alpha helix content increases, the chain contracts until it reaches a minimum radius at approximately 70% helix, after which the chain dimensions expand rapidly. With an alpha helix content of approximately 20%, corresponding to the Ramachandran probability of being in the helical basin, experimentally observed radii of gyration are recovered. Experimental radii are similarly recovered at an alpha helix content of approximately 87%, providing an explanation for the previously puzzling experimental finding that the dimensions of the highly helical methanol-induced unfolded state are experimentally indistinguishable from those of the helix-poor urea-unfolded state. In contrast, the radius of gyration increases monotonically with increasing PPII content, and is always more expanded than the dimensions observed experimentally. These results suggest that PPII is unlikely the sole, dominant preferred conformation for unfolded proteins.  相似文献   

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

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