The role of loop closure propensity in the refolding of Rop protein probed by molecular dynamics simulations

The Curious Case of A31P, a Topology-Switching Mutant of the Repressor of Primer Protein: A Molecular Dynamics Study of Its Folding and Misfolding

The effect of mutations on protein structures is usually rather localized and minor. Finding a mutation that can single-handedly change the fold and/or topology of a protein structure is a rare exception. The A31P mutant of the homodimeric Repressor of primer (Rop) protein is one such exception: This single mutation ─and as demonstrated by two independent crystal structure determinations─ can convert the canonical (left-handed/all-antiparallel) 4-α-helical bundle of Rop to a new form (right-handed/mixed parallel and antiparallel bundle) displaying a previously unobserved "bisecting U" topology. The main problem with understanding the dramatic effect of this mutation on the folding of Rop is to understand its very existence: Most computational methods appear to agree that the mutation should have had no appreciable effect, with the majority of energy minimization methods and protein structure prediction protocols indicating that this mutation is fully consistent with the native Rop structure, requiring only a local and minor change at the mutation site. Here we use two long (10 μs each) molecular dynamics simulations to compare the stability and dynamics of the native Rop versus a hypothetical structure that is identical with the native Rop but is carrying this single Alanine31 to Proline mutation. Comparative analysis of the two trajectories convincingly shows that, in contrast to the indications from energy minimization ─but in agreement with the experimental data─, this hypothetical native-like A31P structure is unstable, with its turn regions almost completely unfolding, even under the relatively mild 320 K NpT simulations that we have used for this study. We discuss the implication of these findings for the folding of the A31P mutant, especially with respect to the proposed model of a double-funneled energy landscape.

Structure and Thermal Stability of wtRop and RM6 Proteins through All-Atom Molecular Dynamics Simulations and Experiments

In the current work we study, via molecular simulations and experiments, the folding and stability of proteins from the tertiary motif of 4-α-helical bundles, a recurrent motif consisting of four amphipathic α-helices packed in a parallel or antiparallel fashion. The focus is on the role of the loop region in the structure and the properties of the wild-type Rop (wtRop) and RM6 proteins, exploring the key factors which can affect them, through all-atom molecular dynamics (MD) simulations and supporting by experimental findings. A detailed investigation of structural and conformational properties of wtRop and its RM6 loopless mutation is presented, which display different physical characteristics even in their native states. Then, the thermal stability of both proteins is explored showing RM6 as more thermostable than wtRop through all studied measures. Deviations from native structures are detected mostly in tails and loop regions and most flexible residues are indicated. Decrease of hydrogen bonds with the increase of temperature is observed, as well as reduction of hydrophobic contacts in both proteins. Experimental data from circular dichroism spectroscopy (CD), are also presented, highlighting the effect of temperature on the structural integrity of wtRop and RM6. The central goal of this study is to explore on the atomic level how a protein mutation can cause major changes in its physical properties, like its structural stability.

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.

Structural plasticity of 4-α-helical bundles exemplified by the puzzle-like molecular assembly of the Rop protein

The dimeric Repressor of Primer (Rop) protein, a widely used model system for the study of coiled-coil 4-α-helical bundles, is characterized by a remarkable structural plasticity. Loop region mutations lead to a wide range of topologies, folding states, and altered physicochemical properties. A protein-folding study of Rop and several loop variants has identified specific residues and sequences that are linked to the observed structural plasticity. Apart from the native state, native-like and molten-globule states have been identified; these states are sensitive to reducing agents due to the formation of nonnative disulfide bridges. Pro residues in the loop are critical for the establishment of new topologies and molten globule states; their effects, however, can be in part compensated by Gly residues. The extreme plasticity in the assembly of 4-α-helical bundles reflects the capacity of the Rop sequence to combine a specific set of hydrophobic residues into strikingly different hydrophobic cores. These cores include highly hydrated ones that are consistent with the formation of interchain, nonnative disulfide bridges and the establishment of molten globules. Potential applications of this structural plasticity are among others in the engineering of bio-inspired materials.