首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
    
Wang P  Klimov DK 《Proteins》2008,70(3):925-937
We use lattice protein models and Monte Carlo simulations to study cotranslational folding of small single domain proteins. We show that the assembly of native structure begins during late extrusion stages, but final formation of native state occurs during de novo folding, when all residues are extruded. There are three main results in our study. First, for the sequences displaying two-state refolding mechanism de novo cotranslational folding pathway differs from that sampled in in vitro refolding. The change in folding pathways is due to partial assembly of native interactions during extrusion that results in different starting conditions for in vitro refolding and for de novo cotranslational folding. For small single domain proteins cotranslational folding is slower than in vitro refolding, but is generally fast enough to be completed before the release from a ribosome. Second, we found that until final stages of biosynthesis cotranslational folding is essentially equilibrium. This observation is explained by low stability of structured states for partially extruded chains. Finally, our data suggest that the proteins, which refold in vitro slowly via intermediates, complete their de novo folding after the release from a ribosome. Comparison of our lattice cotranslational simulations with recent experimental and computational studies is discussed.  相似文献   

2.
    
The hinge-bending proteins provide the most pronounced example of the large-amplitude slow motions in a number of proteins, which are critical for their functioning. They are often used as the test ground for developing novel approaches to the simulation of slow protein dynamics. In the present study, we present the algorithm, which allows physically-consistent simulations of slow protein dynamics in globular proteins. Our algorithm is based on the hierarchical clustering of the correlation patterns (HCCP) technique of domain identification, which allows subdividing the protein into the hierarchy of the rigid-body-like clusters. The clusters are allowed to rotate relative to one another on the automatically identified hinges. The clusters are found in the course of automated, objective and well-tested procedure. In the present communication, our technique is applied to 10 hinge-bending proteins. For each of the proteins, we performed the blind search for the closed conformation, staring from the open one. Resulting closed conformations are compared with the closed states observed in crystallographic structures. It is shown that our technique produces realistic closed conformations for 8 out of 10 studied proteins. This demonstrates that HCCP technique can be used for finding alternative protein conformations and for sampling the slow motions in proteins.  相似文献   

3.
Domains are basic evolutionary units of proteins and most proteins have more than one domain. Advances in domain modeling and collection are making it possible to annotate a large fraction of known protein sequences by a linear ordering of their domains, yielding their architecture. Protein domain architectures link evolutionarily related proteins and underscore their shared functions. Here, we attempt to better understand this association by identifying the evolutionary pathways by which extant architectures may have evolved. We propose a model of evolution in which architectures arise through rearrangements of inferred precursor architectures and acquisition of new domains. These pathways are ranked using a parsimony principle, whereby scenarios requiring the fewest number of independent recombination events, namely fission and fusion operations, are assumed to be more likely. Using a data set of domain architectures present in 159 proteomes that represent all three major branches of the tree of life allows us to estimate the history of over 85% of all architectures in the sequence database. We find that the distribution of rearrangement classes is robust with respect to alternative parsimony rules for inferring the presence of precursor architectures in ancestral species. Analyzing the most parsimonious pathways, we find 87% of architectures to gain complexity over time through simple changes, among which fusion events account for 5.6 times as many architectures as fission. Our results may be used to compute domain architecture similarities, for example, based on the number of historical recombination events separating them. Domain architecture "neighbors" identified in this way may lead to new insights about the evolution of protein function.  相似文献   

4.
Domain insertions and deletions lead to variations in the domain architectures of the proteins from their common ancestor. In this work, we investigated four groups of the RhoGEF-containing proteins from different organisms with domain architectures RhoGEF-PH-SH3, SH3-RhoGEF-PH, RhoGEF-PH, and SH3-RhoGEF defined in the Pfam database. The phylogenetic trees were constructed using each individual domain and/or the combinations of all the domains. The phylogenetic analysis suggests that RhoGEF-PH-SH3 and SH3-RhoGEF-PH might have evolved from RhoGEF-PH through the insertion of SH3 independently, while SH3- RhoGEF of proteins in fruit fly might have evolved from SH3-RhoGEF-PH by the degeneration of PH domain.  相似文献   

5.
    
An interaction between a pair of proteins unique for a particular tissue is denoted as a tissue-specific interaction (TSI). Tissue-specific (TS) proteins always perform TSIs with a limited number of interacting partners. However, it has been claimed that housekeeping (HK) proteins frequently take part in TSIs. This is actually an unusual phenomenon. How a single HK protein mediates TSIs – remains an interesting yet an unsolved question. We have hypothesized that HK proteins have attained a high degree of structural flexibility to modulate TSIs efficiently. We have observed that HK proteins are selected to be intrinsically disordered compared to TS proteins. Therefore, the purposeful adaptation of structural disorder brings out special advantages for HK proteins compared to TS proteins. We have demonstrated that TSIs may play vital roles in shaping the molecular adaptation of disordered regions within HK proteins. We also have noticed that HK proteins, mediating a huge number of TSIs, have a greater portion of their interacting interfaces overlapped with the adjacent disordered segment. Moreover, these HK proteins, mediating TSIs, preferably adapt single domain (SD). We have concluded that HK proteins adapt a high degree of structural flexibility to mediate TSIs. Besides, having a SD along with structural flexibility is more economic than maintaining multiple domains with a rigid structure. This assists them in attaining various structural conformations upon binding to their partners, thereby designing an economically optimum molecular system.  相似文献   

6.
    
Bond‐orientational correlations for finite‐length homopolypeptides and a selected group of denatured proteins are obtained by numerical simulations using a polypeptide model with a potential of mean force. These correlations characterize the stiffness of the polypeptide backbone and are generally described by either an exponential or a power‐law decay in the asymptotic limit. However, for finite length polypeptides and unfolded proteins the correlations significantly deviate from either single exponential or power‐law behavior. A heuristic model is developed to analyze the correlations of homopolypeptides, which depends on the chain length and the side‐chain properties. The model contains power‐law and multi‐exponential terms, the latter which could be interpreted as local persistence lengths. In the asymptotic limit, the model reduces to the expected power‐law behavior. Simulations of denatured proteins show that the power‐law behavior of the correlations is significantly suppressed and only the multi‐exponential term is needed to model the correlations. In addition, average persistence lengths (ranging from 2.0 to 2.5 nm) are obtained from the correlations by fitting single exponentials and shown to be in general agreement with experiments, which also assume single exponential decay. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 312–323, 2016.  相似文献   

7.
    
La Penna G  Furlan S  Perico A 《Biopolymers》2006,83(2):135-147
Molecular dynamics computer simulations were performed for the 25-residue N-terminal tail of the H3 histone protein in the proximity of a DNA segment of 10 base pairs (bp), representing a model for the linker DNA in chromatin. Several least biased configurations were used as initial configurations. The secondary structure content of the protein was increased by the presence of DNA close to it, but the locations of the secondary motifs were different for different initial orientations of the DNA grooves with respect to the protein. As a common feature to all simulations, the electrostatic attraction between negatively charged DNA and positively charged protein was screened by the water solvent and counterbalanced by the intrinsic compaction of the protein due to hydrophobic effects. The protein secondary structure limited the covering of DNA by the protein to 4-5 bp. The degree of compaction and charge density of the bound protein suggests a possible role of H3 tail in a nonspecific bending and plasticity of the linker DNA when the protein is located in the crowded dense chromatin.  相似文献   

8.
    
Demerdash ON  Mitchell JC 《Proteins》2012,80(7):1766-1779
Normal mode analysis has emerged as a useful technique for investigating protein motions on long time scales. This is largely due to the advent of coarse-graining techniques, particularly Hooke's Law-based potentials and the rotational-translational blocking (RTB) method for reducing the size of the force-constant matrix, the Hessian. Here we present a new method for domain decomposition for use in RTB that is based on hierarchical clustering of atomic density gradients, which we call Density-Cluster RTB (DCRTB). The method reduces the number of degrees of freedom by 85-90% compared with the standard blocking approaches. We compared the normal modes from DCRTB against standard RTB using 1-4 residues in sequence in a single block, with good agreement between the two methods. We also show that Density-Cluster RTB and standard RTB perform well in capturing the experimentally determined direction of conformational change. Significantly, we report superior correlation of DCRTB with B-factors compared with 1-4 residue per block RTB. Finally, we show significant reduction in computational cost for Density-Cluster RTB that is nearly 100-fold for many examples.  相似文献   

9.
    
Cytochrome c (cyt c) is an electron‐transfer protein in the respiratory chain of mitochondria. It is known to form polymers, but its polymerization mechanism is still unknown. Dimeric and trimeric cyt c from horse were successfully crystallized by the sitting‐drop vapour‐diffusion method using polyethylene glycol as a precipitating reagent. The crystal of dimeric cyt c belonged to space group P1, with unit‐cell parameters a = 41.8, b = 56.3, c = 60.8 Å, α = 66.3, β = 89.9, γ = 73.7°, whereas that of trimeric cyt c belonged to space group P212121, with unit‐cell parameters a = 57.2, b = 95.7, c = 130.9 Å. Initial structure models showed that the crystals of dimeric and trimeric cyt c contained two dimers and two trimers, respectively, in the asymmetric unit.  相似文献   

10.
    
The 30 kDa membrane‐binding domain of protein 4.1 from human erythrocytes has been expressed in Escherichia coli and crystallized in a form suitable for X‐ray crystallographic study. Crystals were grown using a salting‐in technique. Crystals have a tetragonal plate shape and belong to the C2 space group, with unit‐cell parameters a = 163.9, b = 106.5, c = 93.5 Å, β = 95.5°. The crystals diffract to 2.8 Å resolution.  相似文献   

11.
    
Guardiani C  Cecconi F  Livi R 《Proteins》2008,70(4):1313-1322
Thermal folding molecular dynamics simulations of the domain C5 of Myosin binding protein C were performed using a native-centric model to study the role of three mutations related to Familial Hypertrophic Cardiomyopathy. Mutation of Asn755 causes the largest shift of the folding temperature, and the residue is located in the CFGA' beta-sheet featuring the highest phi-values. The mutation thus appears to reduce the thermodynamic stability in agreement with experimental data. The mutations on Arg654 and Arg668, conversely, cause little change in the folding temperature and they reside in the low phi-value BDE beta-sheet, so that their pathological role cannot be related to impairment of the folding process but possibly to the binding with target molecules. As the typical signature of Domain C5 is the presence of a longer and destibilizing CD-loop with respect to the other Ig-like domains, we completed the work with a bioinformatic analysis of this loop showing a high density of negative charge and low hydrophobicity. This indicates the CD-loop as a natively unfolded sequence with a likely coupling between folding and ligand binding.  相似文献   

12.
Generic features associated with the adsorption of proteins on solid surfaces are reviewed within the framework of the hydrophobic-polar (HP) lattice protein model. The thermodynamic behaviour and structural properties of various HP protein sequences interacting with attractive surfaces have been studied using extensive Wang–Landau sampling with different types of surfaces, each of which attracts either: all monomers, only hydrophobic (H) monomers or only polar (P) monomers, respectively. Consequently, different types of folding behaviour occur for varied surface strengths. Analysis of the combined patterns of various structural observables, e.g. the derivatives of the number of interaction contacts, together with the specific heat, leads to the identification of fundamental categories of folding and transition hierarchies. We also inferred a connection between the transition categories and the relative surface strengths, i.e. the ratios of the surface attractive strengths to the intra-chain attraction among H monomers. Thus, we believe that the folding hierarchies and identification scheme are generic for different HP sequences interacting with attractive surfaces, regardless of the chain length, sequence or surface attraction.  相似文献   

13.
    
Recently Bekker et al. [Bekker G‐J et al. Protein Sci. 2019;28:429–438.] described a computational strategy of applying molecular‐dynamics simulations to estimate the relative stabilities of single‐domain antibodies, and utilized their method to design changes with the aim of increasing the stability of a single‐domain antibody with a known crystal structure. The structure from which they generated potentially stabilizing mutations is an anti‐cholera toxin single domain antibody selected from a naïve library which has relatively low thermal stability, reflected by a melting point of 48°C. Their work was purely theoretical, so to examine their predictions, we prepared the parental and predicted stabilizing mutant single domain antibodies and examined their thermal stability, ability to refold and affinity. We found that the mutation that improved stability the most (~7°C) was one which changed an amino acid in CDR1 from an asparagine to an aspartic acid. This change unfortunately was also accompanied by a reduction in affinity. Thus, while their modeling did appear to successfully predict stabilizing mutations, introducing mutations in the binding regions is problematic. Of further interest, the mutations selected via their high temperature simulations, did improve refolding, suggesting that they were successful in stabilizing the structure at high temperatures and thereby decrease aggregation. Our result should permit them to reassess and refine their model and may one day lead to a usefulin silico approach to protein stabilization.  相似文献   

14.
    
The protein folding problem was first articulated as question of how order arose from disorder in proteins: How did the various native structures of proteins arise from interatomic driving forces encoded within their amino acid sequences, and how did they fold so fast? These matters have now been largely resolved by theory and statistical mechanics combined with experiments. There are general principles. Chain randomness is overcome by solvation-based codes. And in the needle-in-a-haystack metaphor, native states are found efficiently because protein haystacks (conformational ensembles) are funnel-shaped. Order-disorder theory has now grown to encompass a large swath of protein physical science across biology.  相似文献   

15.
Molecular graphics and molecular mechanics techniques have been used to study the mode of ligand binding and mechanism of action of the enzyme phospholipase A2. A substrate-enzyme complex was constructed based on the crystal structure of the apoenzyme. The complex was minimized to relieve initial strain, and the structural and energetic features of the resultant complex analyzed in detail, at the molecular and residue level. The minimized complex was then used as a basis for examining the action of the enzyme on modified substrates, binding of inhibitors to the enzyme, and possible reaction intermediate complexes. The model is compatible with the suggested mechanism of hydrolysis and with experimental data about stereoselectivity, efficiency of hydrolysis of modified substrates, and inhibitor potency. In conclusion, the model can be used as a tool in evaluating new ligands as possible substrates and in the rational design of inhibitors, for the therapeutic treatment of diseases such as rheumatoid arthritis, atherosclerosis, and asthma.  相似文献   

16.
We use an integrated computational approach to reconstruct accurately the transition state ensemble (TSE) for folding of the src-SH3 protein domain. We first identify putative TSE conformations from free energy surfaces generated by importance sampling molecular dynamics for a fully atomic, solvated model of the src-SH3 protein domain. These putative TSE conformations are then subjected to a folding analysis using a coarse-grained representation of the protein and rapid discrete molecular dynamics simulations. Those conformations that fold to the native conformation with a probability (P(fold)) of approximately 0.5, constitute the true transition state. Approximately 20% of the putative TSE structures were found to have a P(fold) near 0.5, indicating that, although correct TSE conformations are populated at the free energy barrier, there is a critical need to refine this ensemble. Our simulations indicate that the true TSE conformations are compact, with a well-defined central beta sheet, in good agreement with previous experimental and theoretical studies. A structured central beta sheet was found to be present in a number of pre-TSE conformations, however, indicating that this element, although required in the transition state, does not define it uniquely. An additional tight cluster of contacts between highly conserved residues belonging to the diverging turn and second beta-sheet of the protein emerged as being critical elements of the folding nucleus. A number of commonly used order parameters to identify the transition state for folding were investigated, with the number of native Cbeta contacts displaying the most satisfactory correlation with P(fold) values.  相似文献   

17.
The graph theory was combined with fluctuation dynamics to investigate the structural communication in four small G proteins, Arf1, H-Ras, RhoA, and Sec4. The topology of small GTPases is such that it requires the presence of the nucleotide to acquire a persistent structural network. The majority of communication paths involves the nucleotide and does not exist in the unbound forms. The latter are almost devoid of high-frequency paths. Thus, small Ras GTPases acquire the ability to transfer signals in the presence of nucleotide, suggesting that it modifies the intrinsic dynamics of the protein through the establishment of regions of hyperlinked nodes with high occurrence of correlated motions. The analysis of communication paths in the inactive (SGDP) and active (SGTP) states of the four G proteins strengthened the separation of the Ras-like domain into two dynamically distinct lobes, i.e. lobes 1 and 2, representing, respectively, the N-terminal and C-terminal halves of the domain. In the framework of this separation, interfunctional states and interfamily differences could be inferred. The structure network undergoes a reshaping depending on the bound nucleotide. Nucleotide-dependent divergences in structural communication reach the maximum in Arf1 and the minimum in RhoA. In Arf1, the nucleotide-dependent paths essentially express a communication between the G box 4 (G4) and distal portions of lobe 1. In the SGDP state, the G4 communicates with the N-term, while, in the SGTP state, the G4 communicates with the switch II. Clear differences could be also found between Arf1 and the other three G proteins. In Arf1, the nucleotide tends to communicate with distal portions of lobe 1, whereas in H-Ras, RhoA, and Sec4 it tends to communicate with a cluster of aromatic/hydrophobic amino acids in lobe 2. These differences may be linked, at least in part, to the divergent membrane anchoring modes that would involve the N-term for the Arf family and the C-term for the Rab/Ras/Rho families.  相似文献   

18.
19.
The insecticidal bacterium Bacillus thuringiensis synthesizes endotoxin Cry proteins of two size classes, 135 and 70 kDa, and both form crystalline inclusions in cells after synthesis. Crystallization of 135-kDa proteins is due to intermolecular attraction of regions in the C-terminal half of the molecule, and the N-terminal half fails to crystallize when synthesized in vivo. Alternatively, endotoxins of the 70-kDa class such as Cry2A and Cry3A, which correspond to the N-terminal half of 135-kDa molecules, crystallize readily after synthesis. Cry molecules of this size class consist of three principal domains, but the domains responsible for crystallization are not known. To identify these domains, chimeric proteins were constructed in which Cry3A Domains I or III, or I and III were substituted for the corresponding domains in truncated Cry1C molecules. Cry1C molecules with only Cry3A Domain III did not crystallize, whereas when Cry3A Domains I and III, or Domain I alone, were substituted, large inclusions were obtained. Except for the chimera consisting of Cry3A Domains I and III and Cry1C Domain II, most chimeras were not as stable as wild-type Cry3A or truncated Cry1C. These results show that Cry3A Domain I plays an important role in its crystallization in vivo.  相似文献   

20.
    
Sikora M  Cieplak M 《Proteins》2011,79(6):1786-1799
We estimate the size of mechanostability for 318 multidomain proteins which are single-chain and contain up to 1021 amino acids. We predict existence of novel types of mechanical clamps in which interdomain contacts play an essential role. Mechanical clamps are structural regions which are the primary source of a protein's resistance to pulling. Among these clamps there is one that opposes tensile stress due to two domains swinging apart. This movement strains and then ruptures the contacts that hold the two domains together. Another clamp also involves tensile stress but it originates from an immobilization of a structural region by a surrounding knot-loop (without involving any disulfide bonds). Still another mechanism involves shear between helical regions belonging to two domains. We also consider the amyloid-prone cystatin C which provides an example of a two-chain 3D domain-swapped protein. We predict that this protein should withstand remarkably large stress, perhaps of order 800 pN, when inducing a shearing strain. The survey is generated through molecular dynamics simulations performed within a structure-based coarse grained model.  相似文献   

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

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