The Effect of Newly Developed OPLS-AA Alanyl Radical Parameters on Peptide Secondary Structure

Activation of Glycyl Radical Enzymes─Multiscale Modeling Insights into Catalysis and Radical Control in a Pyruvate Formate-Lyase-Activating Enzyme

Pyruvate formate-lyase (PFL) is a glycyl radical enzyme (GRE) playing a pivotal role in the metabolism of strict and facultative anaerobes. Its activation is carried out by a PFL-activating enzyme, a member of the radical S-adenosylmethionine (rSAM) superfamily of metalloenzymes, which introduces a glycyl radical into the Gly radical domain of PFL. The activation mechanism is still not fully understood and is structurally based on a complex with a short model peptide of PFL. Here, we present extensive molecular dynamics simulations in combination with quantum mechanics/molecular mechanics (QM/MM)-based kinetic and thermodynamic reaction evaluations of a more complete activation model comprising the 49 amino acid long C-terminus region of PFL. We reveal the benefits and pitfalls of the current activation model, providing evidence that the bound peptide conformation does not resemble the bound protein-protein complex conformation with PFL, with implications for the activation process. Substitution of the central glycine with (S)- and (R)-alanine showed excellent binding of (R)-alanine over unstable binding of (S)-alanine. Radical stabilization calculations indicate that a higher radical stability of the glycyl radical might not be the sole origin of the evolutionary development of GREs. QM/MM-derived radical formation kinetics further demonstrate feasible activation barriers for both peptide and C-terminus activation, demonstrating why the crystalized model peptide system is an excellent inhibitory system for natural activation. This new evidence supports the theory that GREs converged on glycyl radical formation due to the better conformational accessibility of the glycine radical loop, rather than the highest radical stability of the formed peptide radicals.

Protein Stability and Unfolding Following Glycine Radical Formation

Glycine (Gly) residues are particularly susceptible to hydrogen abstraction; which results in the formation of the capto-dative stabilized C_α-centered Gly radical (GLR) on the protein backbone. We examined the effect of GLR formation on the structure of the Trp cage; tryptophan zipper; and the villin headpiece; three fast-folding and stable miniproteins; using all-atom (OPLS-AA) molecular dynamics simulations. Radicalization changes the conformation of the GLR residue and affects both neighboring residues but did not affect the stability of the Trp zipper. The stability of helices away from the radical center in villin were also affected by radicalization; and GLR in place of Gly15 caused the Trp cage to unfold within 1 µs. These results provide new evidence on the destabilizing effects of protein oxidation by reactive oxygen species.

Radical Formation Initiates Solvent-Dependent Unfolding and β-sheet Formation in a Model Helical Peptide

We examined the effects of Cα-centered radical formation on the stability of a model helical peptide, N-Ac-KK(AL)10KK-NH2. Three, 100 ns molecular dynamics simulations using the OPLS-AA force field were carried out on each α-helical peptide in six distinct binary TIP4P water/2,2,2-trifluoroethanol (TFE) mixtures. The α-helicity was at a maximum in 20% TFE, which was inversely proportional to the number of H-bonds between water molecules and the peptide backbone. The radial distribution of TFE around the peptide backbone was highest in 20% TFE, which enhanced helix stability. The Cα-centered radical initiated the formation of a turn within 5 ns, which was a smaller kink at high TFE concentrations, and a loop at lower TFE concentrations. The highest helicity of the peptide radical was measured in 100% TFE. The formation of hydrogen bonds between the peptide backbone and water destabilized the helix, whereas the clustering of TFE molecules around the radical center stabilized the helix. Following radical termination, the once helical structure converted to a β-sheet rich state in 100% water only, and this transition did not occur in the nonradical control peptide. This study gives evidence on how the formation of peptide radicals can initiate α-helical to β-sheet transitions under oxidative stress conditions.