Loop propensity of the sequence YKGQP from staphylococcal nuclease: implications for the folding of nuclease

Conformational dynamics of Ribonuclease Sa and its S48P mutant: Implications for the stability of the mutant protein

The optimization of β-turns has been used as a strategy to increase protein thermal stability. One example is the S48P mutation in Ribonuclease Sa, introduced to optimize a β-turn, which increases the stability of the protein as determined experimentally. Here, we have studied ⁴⁸SYGY⁵¹ β-turn and its S48P mutant from RNase Sa, as a peptide and as part of the protein, using molecular dynamics simulations. The turn propensity of the region ⁴⁸SYGY⁵¹ shows an increase in both the peptide and protein models on S48P mutation. The mutant protein shows an overall decrease in conformational dynamics and a decrease in conformational heterogeneity as compared to the wildtype protein. A comparatively restricted sampling of the φ-ψ region of GLN47, a pre PRO48 residue, in the mutant protein and some local changes in hydrogen bonding patterns involving residues 20-24 might be contributing to the mutant protein stability. In addition, some long-range hydrogen bonding interactions involving the 60s loop and the salt-bridge interaction involving ASP17-ARG63 could also be contributing to the increase in rigidity and stability of the mutant protein.

Impact of Turn Propensity on the Folding Rates of Z34C Protein: Implications for the Folding of Helix-Turn-Helix Motif

The rate-limiting step for the folding of the helix-turn-helix (HTH) protein, Z34C, involves β-turn region ²⁰DPNL²³. This reverse turn has been observed to be part of the transition state in the folding process for Z34C, influencing its folding rates. Molecular dynamics simulations were performed on this turn peptide and its two mutants, D20A and P21A, to study turn formation using GROMOS54A7 force field. We find that this region has a turn propensity of its own, and the highest turn propensity is observed for the wild-type, which correlates well with available experimental results. We also find that a slight unfavorable change in ΔG _{turn folding} causes a drastic change in the folding rates of HTH motif and a mechanistic interpretation is given. Implications of these observations for the folding of the HTH protein Z34C are discussed.

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.

Molecular dynamics simulations of certain mutant peptide models from staphylococcal nuclease reveal that initial hydrophobic collapse associated with turn propensity drives β-hairpin folding

An important nucleation event during the folding of staphylococcal nuclease involves the formation of a β-hairpin by the sequence (21) DTVKLMYKGQPMTFR(35) . Earlier studies show that the turn sequence 'YKGQP' has an important role in the folding of this β-hairpin. To understand the active or passive nature of the turn sequence 'YKGQP' in the folding of the aforementioned β-hairpin sequence, we studied glycine mutant peptides Ac-(2) DTVKLMYGGQPMTFR(16) -NMe (K9G:15), Ac-(2) DTVKLMYKGGPMTFR(16) -NMe (Q11G:15), Ac-(2) DTVKLMYGGGPMTFR(16) -NMe (K9G/Q11G:15), and Ac-(2) DTVKLMGGGGGMTFR(16) -NMe (penta-G:15) by using molecular dynamics simulations, starting with two different unfolded states, polyproline II and extended conformational forms. Further, 5mer mutant turn peptides Ac-(2) YGGQP(6) -NMe (K3G:5), Ac-(2) YKGGP(6) -NMe (Q5G:5), Ac-(2) YGGGP(6) -NMe (K3G/Q5G:5), and Ac-(2) GGGGG(6) -NMe (penta-G:5) were also studied individually. Our results show that an initial hydrophobic collapse and loop closure occurs in all 15mer mutants, but only K9G:15 mutant forms a stable native-like β-hairpin. In the other 15mer mutants, the hydrophobic collapsed state would not proceed to β-hairpin formation. Of the different simulations performed for the penta-G:15 mutant, in only one simulation a nonnative β-hairpin conformation is sampled with highly flexible loop region ((8) GGGGG(12) ), which has no specific conformational preference as a 5mer. While the sequence 'YGGQP' in the K3G:5 simulation shows relatively higher β-turn propensity, the presence of this sequence in K9G:15 peptide seems to be driving the β-hairpin formation. Thus, these results seem to suggest that for the formation of a stable β-hairpin, the initial hydrophobic collapse is to be assisted by a turn propensity. Initial hydrophobic collapse alone is not sufficient to guide β-hairpin formation.

A shorter peptide model from staphylococcal nuclease for the folding-unfolding equilibrium of a beta-hairpin shows that unfolded state has significant contribution from compact conformational states

It is important to understand the conformational features of the unfolded state in equilibrium with folded state under physiological conditions. In this paper, we consider a short peptide model LMYKGQPM from staphylococcal nuclease to model the conformational equilibrium between a hairpin conformation and its unfolded state using molecular dynamics simulation under NVT conditions at 300K using GROMOS96 force field. The free energy landscape has overall funnel-like shape with hairpin conformations sampling the minima. The "unfolded" state has a higher free energy of approximately 12kJ/mol with respect to native hairpin minimum and occupies a plateau region. We find that the unfolded state has significant contributions from compact conformations. Many of these conformations have hairpin-like topology. Further, these compact conformational forms are stabilized by hydrophobic interactions. Conversion between native and non-native hairpins occurs via unfolded states. Frequent conversions between folded and unfolded hairpins are observed with single exponential kinetics. We compare our results with the emerging picture of unfolded state from both experimental and theoretical studies.