Lessons on protein structure from interleukin-4: All disulfides are not created equal.

Interleukin-4 (IL-4) is a hematopoietic cytokine composed by a four-helix bundle stabilized by an antiparallel beta-sheet and three disulfide bonds: Cys3-Cys127, Cys24-Cys65, and Cys46-Cys99. IL-4 is involved in several immune responses associated to infection, allergy, autoimmunity, and cancer. Besides its physiological relevance, IL-4 is often used as a "model" for protein design and engineering. Hence, to understand the role of each disulfide in the structure and dynamics of IL-4, we carried out several spectroscopic analyses (circular dichroism [CD], fluorescence, nuclear magnetic resonance [NMR]), and molecular dynamics (MD) simulations on wild-type IL-4 and four IL-4 disulfide mutants. All disulfide mutants showed loss of structure, altered interhelical angles, and looser core packings, showing that all disulfides are relevant for maintaining the overall fold and stability of the four-helix bundle motif, even at very low pH. In the absence of the disulfide connecting both protein termini Cys3-Cys127, C3T-IL4 showed a less packed protein core, loss of secondary structure (~9%) and fast motions on the sub-nanosecond time scale (lower S2 order parameters and larger τc correlation time), especially at the two protein termini, loops, beginning of helix A and end of helix D. In the absence of Cys24-Cys65, C24T-IL4 presented shorter alpha-helices (14% loss in helical content), altered interhelical angles, less propensity to form the small anti-parallel beta-sheet and increased dynamics. Simultaneously deprived of two disulfides (Cys3-Cys127 and Cys24-Cys65), IL-4 formed a partially folded "molten globule" with high 8-anilino-1-naphtalenesulphonic acid-binding affinity and considerable loss of secondary structure (~50%decrease), as shown by the far UV-CD, NMR, and MD data.

Potentially amyloidogenic conformational intermediates populate the unfolding landscape of transthyretin: insights from molecular dynamics simulations.

Protein aggregation into insoluble fibrillar structures known as amyloid characterizes several neurodegenerative diseases, including Alzheimer's, Huntington's and Creutzfeldt-Jakob. Transthyretin (TTR), a homotetrameric plasma protein, is known to be the causative agent of amyloid pathologies such as FAP (familial amyloid polyneuropathy), FAC (familial amyloid cardiomiopathy) and SSA (senile systemic amyloidosis). It is generally accepted that TTR tetramer dissociation and monomer partial unfolding precedes amyloid fibril formation. To explore the TTR unfolding landscape and to identify potential intermediate conformations with high tendency for amyloid formation, we have performed molecular dynamics unfolding simulations of WT-TTR and L55P-TTR, a highly amyloidogenic TTR variant. Our simulations in explicit water allow the identification of events that clearly discriminate the unfolding behavior of WT and L55P-TTR. Analysis of the simulation trajectories show that (i) the L55P monomers unfold earlier and to a larger extent than the WT; (ii) the single alpha-helix in the TTR monomer completely unfolds in most of the L55P simulations while remain folded in WT simulations; (iii) L55P forms, early in the simulations, aggregation-prone conformations characterized by full displacement of strands C and D from the main beta-sandwich core of the monomer; (iv) L55P shows, late in the simulations, severe loss of the H-bond network and consequent destabilization of the CBEF beta-sheet of the beta-sandwich; (v) WT forms aggregation-compatible conformations only late in the simulations and upon extensive unfolding of the monomer. These results clearly show that, in comparison with WT, L55P-TTR does present a much higher probability of forming transient conformations compatible with aggregation and amyloid formation.

Enthalpic and entropic contributions mediate the role of disulfide bonds on the conformational stability of interleukin-4.

The role of disulfide bridges in the structure, stability, and folding pathways of proteins has been the subject of wide interest in the fields of protein design and engineering. However, the relative importance of entropic and enthalpic contributions for the stabilization of proteins provided by disulfides is not always clear. Here, we perform a detailed analysis of the role of disulfides in the conformational stability of human Interleukin-4 (IL4), a four-helix bundle protein. In order to evaluate the contribution of two out of the three disulfides to the structure and stability of IL4, two IL4 mutants, C3T-IL4 and C24T-IL4, were used. NMR and ANS binding experiments were compatible with altered dynamics and an increase of the nonpolar solvent-accessible surface area of the folded state of the mutant proteins. Chemical and thermal unfolding experiments followed by fluorescence and circular dichroism revealed that both mutant proteins have lower conformational stability than the wild-type protein. Transition temperatures of unfolding decreased 14 degrees C for C3T-IL4 and 10 degrees C for C24T-IL4, when compared to WT-IL4, and the conformational stability, at 25 degrees C, decreased 4.9 kcal/mol for C3T-IL4 and 3.2 kcal/mol for C24T-IL4. Interestingly, both the enthalpy and the entropy of unfolding, at the transition temperature, decreased in the mutant proteins. Moreover, a smaller change in heat capacity of unfolding was also observed for the mutants. Thus, disulfide bridges in IL4 play a critical role in maintaining the thermodynamic stability and core packing of the helix bundle.

A structural model of an amyloid protofilament of transthyretin.

A docking-and-alignment protocol was devised in order to build amyloid protofilaments of Transthyretin (TTR), starting from partially disrupted TTR monomeric subunits and based on experimentally available information. The docking approach is driven by a combination of shape complementarity and energetic criteria, and uses constraints derived from experimental data obtained for the fibrillar state. The dimeric structures obtained were then subjected to an alignment scheme followed by clustering analysis, producing a collection of protofilaments with distinct geometric properties. The selected protofilament model presented here does agree with known experimental data and general amyloid properties; it is formed by two extended continuous beta-sheets with the beta-strands perpendicular to the main axis of the protofilament and a helical twist with a period of approximately 48 beta-strands. This TTR proto-filament model may be an important step in the understanding of the molecular mechanisms of TTR aggregation, as well as, a valuable instrument in drug design strategies against amyloid diseases.

Towards data warehousing and mining of protein unfolding simulation data.

OBJECTIVES: The prediction of protein structure and the precise understanding of protein folding and unfolding processes remains one of the greatest challenges in structural biology and bioinformatics. Computer simulations based on molecular dynamics (MD) are at the forefront of the effort to gain a deeper understanding of these complex processes. Currently, these MD simulations are usually on the order of tens of nanoseconds, generate a large amount of conformational data and are computationally expensive. More and more groups run such simulations and generate a myriad of data, which raises new challenges in managing and analyzing these data. Because the vast range of proteins researchers want to study and simulate, the computational effort needed to generate data, the large data volumes involved, and the different types of analyses scientists need to perform, it is desirable to provide a public repository allowing researchers to pool and share protein unfolding data. METHODS: To adequately organize, manage, and analyze the data generated by unfolding simulation studies, we designed a data warehouse system that is embedded in a grid environment to facilitate the seamless sharing of available computer resources and thus enable many groups to share complex molecular dynamics simulations on a more regular basis. RESULTS: To gain insight into the conformational fluctuations and stability of the monomeric forms of the amyloidogenic protein transthyretin (TTR), molecular dynamics unfolding simulations of the monomer of human TTR have been conducted. Trajectory data and meta-data of the wild-type (WT) protein and the highly amyloidogenic variant L55P-TTR represent the test case for the data warehouse. CONCLUSIONS: Web and grid services, especially pre-defined data mining services that can run on or 'near' the data repository of the data warehouse, are likely to play a pivotal role in the analysis of molecular dynamics unfolding data.

Protein folding and unfolding simulations: a new challenge for data mining.

One of the unsolved paradigms in molecular biology is the protein folding problem. In recent years, with the identification of several diseases as protein folding disorders and with the explosion of genome information and the need for efficient ways to predict protein structure, protein folding became a central issue in molecular sciences research. Using molecular dynamics unfolding simulations of an amyloidogenic protein--transthyretin--as an example, we put forward a series of ideas on how simulations of this type may be used to infer rules and unfolding behavior in amyloidogenic proteins, and to extrapolate rules for protein folding in different structural classes of proteins. These, in turn, could help in the development of protein structure prediction methods. The need to analyse different proteins and to run multiple simulations creates a huge amount of data which has to be stored, managed, analyzed and shared (database and Grid technology; data mining). Once the data is captured, the next challenge is to find meaningful patterns (associations, correlations, clusters, rules, relationships) among molecular properties, or their relative importance at different stages of the folding or unfolding processes. This clearly puts new and interesting challenges to the bioinformatics community.