Molecular dynamics study of an insertion/duplication mutant of bacteriophage T4 lysozyme reveals the nature of α→β transition in full protein context

Molecular dynamics simulations of an engineered T4 lysozyme exclude helix to sheet transition, and provide insights into long distance, intra-protein switchable motion

An engineered variant of T4 lysozyme serves as a model for studying induced remote conformational changes in a full protein context. The design involves a duplicated surface helix, flanked by two loops, that switches between two different conformations spanning about 20 Å. Molecular dynamics simulations of the engineered protein, up to 1 μs, rule out α-helix to β-sheet transitions within the duplicated helix as suggested by others. These simulations highlight how the use of different force fields can lead to radical differences in the structure of the protein. In addition, Markov state modeling and transition path theory were employed to map a 6.6 μs simulation for possible early intermediate states and to provide insights into the onset of the switching motion. The putative intermediates involve the folding of one helical turn in the C-terminal loop through energy driven, sequential rearrangement of nearby salt bridges around the key residue Arg63. These results provide a first step towards understanding the energetics and dynamics of a rather complicated intra-protein motion.

Complex kinetics and residual structure in the thermal unfolding of yeast triosephosphate isomerase

BACKGROUND

Saccharomyces cerevisiae triosephosphate isomerase (yTIM) is a dimeric protein that shows noncoincident unfolding and refolding transitions (hysteresis) in temperature scans, a phenomenon indicative of the slow forward and backward reactions of the native-unfolded process. Thermal unfolding scans suggest that no stable intermediates appear in the unfolding of yTIM. However, reported evidence points to the presence of residual structure in the denatured monomer at high temperature.

RESULTS

Thermally denatured yTIM showed a clear trend towards the formation of aggregation-prone, β-strand-like residual structure when pH decreased from 8.0 to 6.0, even though thermal unfolding profiles retained a simple monophasic appearance regardless of pH. However, kinetic studies performed over a relatively wide temperature range revealed a complex unfolding mechanism comprising up to three observable phases, with largely different time constants, each accompanied by changes in secondary structure. Besides, a simple sequential mechanism is unlikely to explain the observed variation of amplitudes and rate constants with temperature. This kinetic complexity is, however, not linked to the appearance of residual structure. Furthermore, the rate constant for the main unfolding phase shows small, rather unvarying values in the pH region where denatured yTIM gradually acquires a β-strand-like conformation. It appears, therefore, that the residual structure has no influence on the kinetic stability of the native protein. However, the presence of residual structure is clearly associated with increased irreversibility.

CONCLUSIONS

The slow temperature-induced unfolding of yeast TIM shows three kinetic phases. Rather than a simple sequential pathway, a complex mechanism involving off-pathway intermediates or even parallel pathways may be operating. β-strand-type residual structure, which appears below pH 8.0, is likely to be associated with increased irreversible aggregation of the unfolded protein. However, this denatured form apparently accelerates the refolding process.

Environmental polarity induces conformational transitions in a helical peptide sequence from bacteriophage T4 lysozyme and its tandem duplicate: a molecular dynamics simulation study

Our recent molecular dynamics (MD) simulation of an insertion/duplication mutant 'L20' of bacteriophage T4 lysozyme demonstrated a solvent induced α→β transition in a loosely held duplicate helical region, while α-helical conformation in the parent region was relatively stabilized by its tertiary interactions with the neighboring residues. The solution NMR of the parent helical sequence, sans its protein context, showed no inherent tendency to adopt a particular secondary structure in pure water but showed α-helical propensity in TFE/water and SDS micelles. In this study we investigate the conformational preference of the 'parent' and 'duplicate' sequences, sans the protein context, in pure water and an apolar TFE/water solution. Apolar TFE/water solution is a model for non-polar protein context. We performed MD simulations of the two peptides, in explicit water and 80% (v/v) TFE/water, using GROMOS 53a6 force field, at 300 K and 1 bar (under NPT conditions). We show that in TFE/water mixture, salt bridges are stabilized by apolar TFE molecules and main chain-main chain hydrogen bonds promote the α-helical conformation, particularly in the duplicate peptide. Solvent exposure, in pure water, resulted in an α→β transition to form a triple stranded β-sheet structure in the 'duplicate' sequence, with a rare psi-loop topology, while a mixture of turn/bend conformations were adopted by the 'parent' sequence. Thus the differences in conformational preference of the parent and duplicate sequence sans protein context, in pure water and TFE/water, implicate the importance of the environment polarity in dictating the peptide conformation. Mechanism of folding of the observed psi-loop in the duplicate sequence gives insights into folding of this rare β-sheet topology.