Grid computing solutions for distributed repositories of protein folding and unfolding simulations

The P-found protein folding and unfolding simulation repository is designed to allow scientists to perform analyses across large, distributed simulation data sets. There are two storage components in P-found: a primary repository of simulation data and a data warehouse. Here we demonstrate how grid technologies can support multiple, distributed P-found installations. In particular we look at two aspects, first how grid data management technologies can be used to access the distributed data warehouses; and secondly, how the grid can be used to transfer analysis programs to the primary repositories --- this is an important and challenging aspect of P-found because the data volumes involved are too large to be centralised. The grid technologies we are developing with the P-found system will allow new large data sets of protein folding simulations to be accessed and analysed in novel ways, with significant potential for enabling new scientific discoveries.

P-found: grid-enabling distributed repositories of protein folding and unfolding simulations for data mining

The P-found protein folding and unfolding simulation repository is designed to allow scientists to perform data mining and other analyses across large, distributed simulation data sets. There are two storage components in P-found: a primary repository of simulation data that is used to populate the second component, and a data warehouse that contains important molecular properties. These properties may be used for data mining studies. Here we demonstrate how grid technologies can support multiple, distributed P-found installations. In particular, we look at two aspects: firstly, how grid data management technologies can be used to access the distributed data warehouses; and secondly, how the grid can be used to transfer analysis programs to the primary repositories — this is an important and challenging aspect of P-found, due to the large data volumes involved and the desire of scientists to maintain control of their own data. The grid technologies we are developing with the P-found system will allow new large data sets of protein folding simulations to be accessed and analysed in novel ways, with significant potential for enabling scientific discovery.

Mechanical and chemical unfolding of a single protein: A comparison

Is the mechanical unraveling of protein domains by atomic force microscopy (AFM) just a technological feat or a true measurement of their unfolding? By engineering a protein made of tandem repeats of identical Ig modules, we were able to get explicit AFM data on the unfolding rate of a single protein domain that can be accurately extrapolated to zero force. We compare this with chemical unfolding rates for untethered modules extrapolated to 0 M denaturant. The unfolding rates obtained by the two methods are the same. Furthermore, the transition state for unfolding appears at the same position on the folding pathway when assessed by either method. These results indicate that mechanical unfolding of a single protein by AFM does indeed reflect the same event that is observed in traditional unfolding experiments. The way is now open for the extensive use of AFM to measure folding reactions at the single-molecule level. Single-molecule AFM recordings have the added advantage that they define the reaction coordinate and expose rare unfolding events that cannot be observed in the absence of chemical denaturants.

Two-dimensional IR correlation spectroscopy: Sequential events in the unfolding process of the λ Cro-V55C repressor protein

A question often posed in protein folding/unfolding studies is whether the process is fully cooperative or whether it contains sequential elements. To address this question, one needs tools capable of resolving different events. It seems that, at least in certain cases, two-dimensional (2D) IR correlation spectroscopy can provide answers to this question. To illustrate this point, we have turned to the Cro-V55C dimer of the λ Cro repressor, a protein known to undergo thermal unfolding in two discrete steps through a stable equilibrium intermediate. The secondary structure of this intermediate is compatible with that of a partially unfolded protein and involves a reorganization of the N terminus, whereas the antiparallel β-ribbon formed by the C-terminal part of each subunit remains largely intact. To establish whether the unfolding process involves sequential events, we have performed a 2D correlation analysis of IR spectra recorded over the temperature range of 20–95°C. The 2D IR correlation analysis indeed provides evidence for a sequential formation of the stable intermediate, which is created in three (closely related) steps. A first step entails the unfolding of the short N-terminal β-strand, followed by the unfolding of the α-helices in a second step, and the third step comprises the reorganization of the remaining β-sheet and of some unordered segments in the protein. The complete unfolding of the stable intermediate at higher temperatures also undergoes sequential events that ultimately end with the breaking of the H bonds between the two β-strands at the dimer interface.

Pressure-induced protein-folding/unfolding kinetics

We use an off-lattice minimalist model to describe the effects of pressure in slowing down the folding/unfolding kinetics of proteins when subjected to increasingly larger pressures. The potential energy function used to describe the interactions between beads in the model includes the effects of pressure on the pairwise interaction of hydrophobic groups in water. We show that pressure affects the participation of contacts in the transition state. More significantly, pressure exponentially decreases the chain reconfigurational diffusion coefficient. These results are consistent with experimental results on the kinetics of pressure-denaturation of staphylococcal nuclease.

Staphylococcal protein A: Unfolding pathways, unfolded states, and differences between the B and E domains

We present multiple native and denaturation simulations of the B and E domains of the three-helix bundle protein A, totaling 60 ns. The C-terminal helix (H3) consistently denatures later than either of the other two helices and contains residual helical structure in the denatured state. These results are consistent with experiments suggesting that the isolated H3 fragment is more stable than H1 and H2 and that H3 forms early in folding. Interestingly, the denatured state of the B domain is much more compact than that of the E domain. This sequence-dependent effect on the dimensions of the denatured state and the lack of correlation with structure suggest that the radius of gyration can be a misleading reaction coordinate for unfolding/folding. Various unfolding and refolding events are observed in the denaturation simulations. In some cases, the transitions are facilitated through interactions with other portions of the protein—contact-assisted helix formation. In the native simulations, the E domain is very stable: after 6 ns, the Cα root-mean-square deviation from the starting structure is less than 1.4 Å. In contrast, the native state of the B domain deviates more and its inter-helical angles fluctuate. In apparent contrast, we note that the B domain is thermodynamically more stable than the E domain. The simulations suggest that the increased stability of the B domain may be due to heightened mobility, and therefore entropy, in the native state and decreased mobility/entropy in the more compact denatured state.

A scrapie-like unfolding intermediate of the prion protein domain PrP(121–231) induced by acidic pH

The infectious agent of transmissible spongiform encephalopathies is believed to consist of an oligomeric isoform, PrPSc, of the monomeric cellular prion protein, PrPC. The conversion of PrPC to PrPSc is characterized by a decrease in α-helical structure, an increase in β-sheet content, and the formation of PrPSc amyloid. Whereas the N-terminal part of PrPC comprising residues 23–120 is flexibly disordered, its C-terminal part, PrP(121–231), forms a globular domain with three α-helices and a small β-sheet. Because the segment of residues 90–231 is protease-resistant in PrPSc, it is most likely structured in the PrPSc form. The conformational change of the segment containing residues 90–120 thus constitutes the minimal structural difference between PrPC and a PrPSc monomer. To test whether PrP(121–231) is also capable to undergo conformational transitions, we analyzed its urea-dependent unfolding transitions at neutral and acidic pH. We identified an equilibrium unfolding intermediate of PrP(121–231) that is exclusively populated at acidic pH and shows spectral characteristics of a β-sheet protein. The intermediate is in rapid equilibrium with native PrP(121–231), significantly populated in the absence of urea at pH 4.0, and may have important implications for the presumed formation of PrPSc during endocytosis.

Evaluation of novel Streptococcus pyogenes vaccine candidates incorporating multiple conserved sequences from the C-repeat region of the M-protein

A major challenge for Streptococcus pyogenes vaccine development is the identification of epitopes that confer protection from infection by multiple S. pyogenes M-types. Here we have identified and characterised the distribution of common variant sequences from individual repeat units of the C-repeat region (CRR) of M-proteins representing 77 different M-types. Three polyvalent fusion vaccine candidates (SV1, SV2 and SV3) incorporating the most common variants were subsequently expressed and purified, and demonstrated to be alpha-helical by Circular Dichroism (CD), a secondary conformational characteristic of the CRR in the M-protein. Antibodies raised against each of these constructs recognise M-proteins that vary in their CRR, and bind to the surface of multiple S. pyogenes isolates. Antibodies raised against SV1, containing five variant sequences, also kill heterologous S. pyogenes isolates in in vitro bactericidal assays. Further structural characterisation of this construct demonstrated the conformation of SV1 was stable at different pHs, and thermal unfolding of SV1 a reversible process. Our findings demonstrate that linkage of multiple variant sequences into a single recombinant construct overcomes the need to embed the variant sequences in foreign helix promoting flanking sequences for conformational stability, and demonstrates the viability of the polyvalent candidates as global S. pyogenes vaccine candidates.

Effect of Crowding Agents, Signal Peptide, and Chaperone SecB on the Folding and Aggregation of E. coli Maltose Binding Protein

SecB, a soluble cytosolic chaperone component of the Secexport pathway, binds to newly synthesized precursor proteins and prevents their premature aggregation and folding and subsequently targets them to the translocation machinery on the membrane. PreMBP, the precursor form of maltose binding protein, has a 26-residue signal sequence attached to the N-terminus of MBP and is a physiological substrate of SecB. We examine the effect of macromolecular crowding and SecB on the stability and refolding of denatured preMBP and MBP. PreMBP was less stable than MBP (ΔTm =7( 0.5 K) in both crowded and uncrowded solutions. Crowding did not cause any substantial changes in the thermal stability ofMBP(ΔTm=1(0.4 K) or preMBP (ΔTm=0(0.6 K), as observed in spectroscopically monitored thermal unfolding experiments. However, both MBP and preMBP were prone to aggregation while refolding under crowded conditions. In contrast to MBP aggregates, which were amorphous, preMBP aggregates form amyloid fibrils.Under uncrowded conditions, a molar excess of SecB was able to completely prevent aggregation and promote disaggregation of preformed aggregates of MBP. When a complex of the denatured protein and SecB was preformed, SecB could completely prevent aggregation and promote folding of MBP and preMBP even in crowded solution. Thus, in addition to maintaining substrates in an unfolded, export-competent conformation, SecB also suppresses the aggregation of its substrates in the crowded intracellular environment. SecB is also able to promote passive disaggregation of macroscopic aggregates of MBP in the absence of an energy source such as ATP or additional cofactors. These experiments also demonstrate that signal peptide can reatly influence protein stability and aggregation propensity.

Adsorption and Unfolding of Lysozyme at a Polarized Aqueous-Organic Liquid Interface

The adsorption of proteins at the interface between two immiscible electrolyte solutions has been found to be key to their bioelectroactivity at such interfaces. Combined with interfacial complexation of organic phase anions by cationic proteins, this adsorption process may be exploited to achieve nanomolar protein detection. In this study, replica exchange molecular dynamics simulations have been performed to elucidate for the first time the molecular mechanism of adsorption and subsequent unfolding of hen egg white lysozyme at low pH at a polarized 1,2-dichloroethane/water interface. The unfolding of lysozyme was observed to occur as soon as it reaches the organic−aqueous interface,which resulted in a number of distinct orientations at the interface. In all cases, lysozyme interacted with the organic phase through regions rich in nonpolar amino acids, such that the side chains are directed toward the organic phase, whereas charged and polar residues were oriented toward the aqueous phase. By contrast, as expected, lysozyme in neat water at low pH does not exhibit significant structural changes. These findings demonstrate the key influence of the organic phase upon adsorption of lysozyme under the influence of an electric field, which results in the unfolding of its structure.

Thermodynamics of metal ion binding and denaturation of a calcium binding protein from Entamoeba histolytica

The thermodynamics of tie binding of calcium and magnesium ions to a calcium binding protein from Entamoeba histolytica was investigated by isothermal titration calorimetry (ITC) in 20 mM MOPS buffer (pH 7.0) at 20 degrees C. Enthalpy titration curves of calcium show the presence of four Ca2+ binding sites, There exist two low-affinity sites for Ca2+, both of which are exothermic in nature and with positive cooperative interaction between them. Two other high affinity sites for Ca2+ exist of which one is endothermic and the other exothermic, again with positive cooperative interaction. The binding constants for Ca2+ at the four sites have been verified by a competitive binding assay, where CaBP competes with a chromophoric chelator 5, 5'-Br-2 BAPTA to bind Ca2+ and a Ca2+ titration employing intrinsic tyrosine fluorescence of the protein, The enthalpy of titration of magnesium in the absence of calcium is single site and endothermic in nature. In the case of the titrations performed using protein presaturated with magnesium, the amount of heat produced is altered. Further, the interaction between the high-affinity sites changes to negative cooperativity. No exchange of heat was observed throughout the addition of magnesium in the presence of 1 mM calcium, Titrations performed on a cleaved peptide comprising the N-terminus and the central linker show the existence of two Ca2+ specific sites, These results indicate that this CaBP has one high-affinity Ca-Mg site, one high-affinity Ca-specific site, and two low-affinity Ca-specific sites. The thermodynamic parameters of the binding of these metal ions were used to elucidate the energetics at the individual site(s) and the interactions involved therein at various concentrations of the denaturant, guanidine hydrochloride, ranging from 0.05 to 6.5 M. Unfolding of the protein was also monitored by titration calorimetry as a function of the concentration of the denaturant. These data show that at a GdnHCl concentration of 0.25 M the binding affinity for the Mg2+ ion is lost and there are only two sites which can bind to Ca2+, with substantial loss cooperativity. At concentrations beyond 2.5 M GdnHCl, at which the unfolding of the tertiary structure of this protein is observed by near UV CD spectroscopy, the binding of Ca2+ ions is lost. We thus show that the domain containing the two low-affinity sites is the first to unfold in the presence of GdnHCl. Control experiments with change in ionic strength by addition of KCI in the range 0.25-1 M show the existence of four sites with altered ion binding parameters.

On the relationship of thermodynamic parameters with the buried surface area in protein-ligand complex formation

Prediction of thermodynamic parameters of protein-protein and antigen-antibody complex formation from high resolution structural parameters has recently received much attention, since an understanding of the contributions of different fundamental processes like hydrophobic interactions, hydrogen bonding, salt bridge formation, solvent reorganization etc. to the overall thermodynamic parameters and their relations with the structural parameters would lead to rational drug design. Using the results of the dissolution of hydrocarbons and other model compounds the changes in heat capacity (DeltaCp), enthalpy (DeltaH) and entropy (DeltaS) have been empirically correlated with the polar and apolar surface areas buried during the process of protein folding/unfolding and protein-ligand complex formation. In this regard, the polar and apolar surfaces removed from the solvent in a protein-ligand complex have been calculated from the experimentally observed values of changes in heat capacity (DeltaCp) and enthalpy (DeltaH) for protein-ligand complexes for which accurate thermodynamic and high resolution structural data are available, and the results have been compared with the x-ray crystallographic observations. Analyses of the available results show poor correlation between the thermodynamic and structural parameters. Probable reasons for this discrepancy are mostly related with the reorganization of water accompanying the reaction which is indeed proven by the analyses of the energetics of the binding of the wheat germ agglutinin to oligosaccharides.

Insight into the early stages of thermal unfolding of peanut agglutinin by molecular dynamics simulations

Peanut agglutinin is a homotetrameric nonglycosylated protein. The protein has a unique open quaternary structure. Molecular dynamics simulations have been employed follow the atomistic details of its unfolding at different temperatures. The early events of the deoligomerization of the protein have been elucidated in the present study. Simulation trajectories of the monomer as well as those of the tetramer have been compared and the tetramer is found to be substantially more stable than its monomeric counterpart. The tetramer shows retention of most of its.. secondary structure but considerable loss of the tertiary structure at high temperature. e generation of a This observation impies the molten globule-like intermediate in the later stages of deoligomerization. The quaternary structure of the protein has weakened to a large extent, but none of the subunits are separated. In addition, the importance of the metal-binding to the stability of the protein structure has also been investigated. Binding of the metal ions not only enhances the local stability of the metal-ion binding loop, but also imparts a global stability to the overall structure. The dynamics of different interfaces vary significantly as probed through interface clusters. The differences are substantially enhanced at higher temperatures. The dynamics and the stability of the interfaces have been captured mainly by cluster analysis, which has provided detailed information on the thermal deoligomerization of the protein.

Unfolding of Plasmodium falciparum Triosephosphate Isomerase in Urea and Guanidinium Chloride: Evidence for a Novel Disulfide Exchange Reaction in a Covalently Cross-Linked Mutant†

The conformational stability of Plasmodium falciparum triosephosphate isomerase (TIMWT) enzyme has been investigated in urea and guanidinium chloride (GdmCl) solutions using circular dichroism, fluorescence, and size-exclusion chromatography. The dimeric enzyme is remarkably stable in urea solutions. It retains considerable secondary, tertiary, and quaternary structure even in 8 M urea. In contrast, the unfolding transition is complete by 2.4 M GdmCl. Although the secondary as well as the tertiary interactions melt before the perturbation of the quaternary structure, these studies imply that the dissociation of the dimer into monomers ultimately leads to the collapse of the structure, suggesting that the interfacial interactions play a major role in determining multimeric protein stability. The C-m(urea)/C-m(GdmCl) ratio (where C-m is the concentration of the denaturant required at the transition midpoint) is unusually high for triosephosphate isomerase as compared to other monomeric and dimeric proteins. A disulfide crosslinked mutant protein (Y74C) engineered to form two disulfide cross-links across the interface (13-74') and (13'-74) is dramatically destablized in urea. The unfolding transition is complete by 6 M urea and involves a novel mechanism of dimer dissociation through intramolecular thiol-disulfide exchange.