Atomically Detailed Description of the Unfolding of α-Lactalbumin by the Combined Use of Experiments and Simulations

Molecular dynamics simulations at high temperatures of the Aeropyrum pernix L7Ae thermostable protein: Insight into the unfolding pathway

Most life forms on earth live at temperatures below 50 °C. Within these organisms are proteins that form the three-dimensional structures essential to their biological activity and function. However, some thermophilic life forms can resist higher temperatures and have corresponding adaptations to preserve protein function at these high temperatures. Among the structural factors responsible for this resistance of thermophilic proteins to high temperatures is the presence of additional hydrogen bonds in the thermophilic proteins, which means that the structure of the protein is more resistant to unfolding. Similarly, thermostable proteins are rich in structure-stabilizing salt bridges and/or disulfide bridges. In this context, we perform multiple replica molecular dynamics simulations at different temperatures on the Aeropyrum pernix (L7Ae) protein (from the crenarchaeal species A. pernix), known for its high melting temperature, and this in the aim to elucidate the structural factors responsible for its high thermostability. The results reveal that between the most sensitive regions of the protein to the increase of temperature are the loops L1, and L5, which surround the hydrophobic core region of the protein, besides the loop L9, and the C-terminal α5 region. This latter is the longer alpha helix of the protein secondary structure motifs and it is the first to be denaturated at 450 K, while the rest of the protein secondary structure motifs at this temperature were intact. The mechanism of unfolding that follows this protein at 550 K is similar to other thermophile proteins found in literature, with the opening of the loops that surround the hydrophobic core of the protein. So, the latter is completely exposed to the solvent, and partially denatured. The total denaturation process of the protein takes an average time of 40 ns to be achieved. Our investigation also shows that all the calculated salt bridges, with distances less than or equal to 6 A°, are on the periphery part of the protein, exposed to the solvent. However, the hydrophobic core of the protein is not involved in the formation of salt bridges, but rather with formation of some important hydrogen bondings that still persist even at 450 K. So, optimizing hydrogen bonding, near or within the core region, at high temperatures is a strategy that follows this thermostable protein to protect its hydrophobic core from denaturation, and ensure the thermal stability of the protein.

Impacts of Ca2+ cation and temperature on bovine α-lactalbumin secondary structures and foamability - Insights from computational molecular dynamics

We used computational molecular dynamics (MD) to assess molecular conformations of apo- and holo-forms (respectively without and with Ca²⁺) of bovine α-lactalbumin (α-La) at different temperatures, and to correlate them with the protein's foaming properties. At 4 °C and 25 °C no major protein conformation changes occurred. At 75 °C, lots of changes were evidenced: the Ca²⁺ depletion triggered the complete loss of h2b, h3c helices and S1, S2 and S3 β-sheets, and partial losses of H1, H2 and H3 α-helices. The absence of Ca²⁺ in apo-α-La and its leaving from holo-α-La triggered electrostatic repulsion among Asp82, Asp84 and Asp87, leading to the formation of a hydrophobic cluster involving Phe9, Phe31, Ile1, Va42, Ile55, Phe80 and Leu81. These conformational changes were related to an interfacial tension decrease and to a foaming capacity increase, for both apo-α-La and holo-α-La. This study exemplifies how powerful MD is as a tool to provide a better understanding of the molecular origins of food proteins' techno-functionalities.

Sequential four-state folding/unfolding of goat α-lactalbumin and its N-terminal variants

Equilibria and kinetics of folding/unfolding of α-lactalbumin and its two N-terminal variants were studied by circular dichroism spectroscopy. The two variants were wild-type recombinant and Glu1-deletion (E1M) variants expressed in Escherichia coli. The presence of an extra methionine at the N terminus in recombinant α-lactalbumin destabilized the protein by 2 kcal/mol, while the stability was recovered in the E1M variant in which Glu1 was replaced by Met1. Kinetic folding/unfolding reactions of the proteins, induced by stopped-flow concentration jumps of guanidine hydrochloride, indicated the presence of a burst-phase in refolding, and gave chevron plots with significant curvatures in both the folding and unfolding limbs. The folding-limb curvature was interpreted in terms of accumulation of the burst-phase intermediate. However, there was no burst phase observed in the unfolding kinetics to interpret the unfolding-limb curvature. We thus assumed a sequential four-state mechanism, in which the folding from the burst-phase intermediate takes place via two transition states separated by a high-energy intermediate. We estimated changes in the free energies of the burst-phase intermediate and two transition states, caused by the N-terminal variations and also by the presence of stabilizing calcium ions. The Φ values at the N terminus and at the Ca(2+)-binding site thus obtained increased successively during folding, demonstrating the validity of the sequential mechanism. The stability and the folding behavior of the E1M variant were essentially identical to those of the authentic protein, allowing us to use this variant as a pseudo-wild-type α-lactalbumin in future studies.

Thermal unfolding simulations of bacterial flagellin: insight into its refolding before assembly

Flagellin is the subunit of the bacterial filament, the micrometer-long propeller of a bacterial flagellum. The protein is believed to undergo unfolding for transport through the channel of the filament and to refold in a chamber at the end of the channel before being assembled into the growing filament. We report a thermal unfolding simulation study of S. typhimurium flagellin in aqueous solution as an attempt to gain atomic-level insight into the refolding process. Each molecule comprises two filament-core domains {D0, D1} and two hypervariable-region domains {D2, D3}. D2 can be separated into subdomains D2a and D2b. We observed a similar unfolding order of the domains as reported in experimental thermal denaturation. D2a and D3 exhibited high thermal stability and contained persistent three-stranded beta-sheets in the denatured state which could serve as folding cores to guide refolding. A recent mutagenesis study on flagellin stability seems to suggest the importance of the folding cores. Using crude size estimates, our data suggests that the chamber might be large enough for either denatured hypervariable-region domains or filament-core domains, but not whole flagellin; this implicates a two-staged refolding process.

Unfolding pathways of goat alpha-lactalbumin as revealed in multiple alignment of molecular dynamics trajectories

Molecular dynamics simulations of protein unfolding were performed at an elevated temperature for the authentic and recombinant forms of goat alpha-lactalbumin. Despite very similar three-dimensional structures, the two forms have significantly different unfolding rates due to an extra N-terminal methionine in the recombinant protein. To identify subtle differences between the two forms in the highly stochastic kinetics of unfolding, we classified the unfolding trajectories using the multiple alignment method based on the analogy between the biological sequences and the molecular dynamics trajectories. A dendrogram derived from the multiple trajectory alignment revealed a clear difference in the unfolding pathways of the authentic and recombinant proteins, i.e. the former reached the transition state in an all-or-none manner while the latter unfolded less cooperatively. It was also found in the classification that the two forms of the protein shared a common transition state structure, which was in excellent agreement with the transition state structure observed experimentally in the Phi-value analysis.