Connecting Protein Conformational Dynamics with Catalytic Function As Illustrated in Dihydrofolate Reductase

Combined quantum mechanics/molecular mechanics molecular dynamics simulations reveal that the M20 loop conformational dynamics of dihydrofolate reductase (DHFR) is severely restricted at the transition state of the hydride transfer as a result of the M42W/G121V double mutation. Consequently, the double-mutant enzyme has a reduced entropy of activation, i.e., increased entropic barrier, and altered temperature dependence of kinetic isotope effects in comparison with those of wild-type DHFR. Interestingly, in both wild-type DHFR and the double mutant, the average donor–acceptor distances are essentially the same in the Michaelis complex state (∼3.5 Å) and the transition state (2.7 Å). It was found that an additional hydrogen bond is formed to stabilize the M20 loop in the closed conformation in the M42W/G121V double mutant. The computational results reflect a similar aim designed to knock out precisely the dynamic flexibility of the M20 loop in a different double mutant, N23PP/S148A

Conformational Equilibrium of N‑Myristoylated cAMP-Dependent Protein Kinase A by Molecular Dynamics Simulations

The catalytic subunit of protein kinase A (PKA-C) is subject to several post- or cotranslational modifications that regulate its activity both spatially and temporally. Among those, N-myristoylation increases the kinase affinity for membranes and might also be implicated in substrate recognition and allosteric regulation. Here, we investigated the effects of N-myristoylation on the structure, dynamics, and conformational equilibrium of PKA-C using atomistic molecular dynamics simulations. We found that the myristoyl group inserts into the hydrophobic pocket and leads to a tighter packing of the A-helix against the core of the enzyme. As a result, the conformational dynamics of the A-helix are reduced and its motions are more coupled with the active site. Our simulations suggest that cation−π interactions among W30, R190, and R93 are responsible for coupling these motions. Two major conformations of the myristoylated N-terminus are the most populated: a long loop (LL conformation), similar to Protein Data Bank (PDB) entry 1CMK, and a helix–turn–helix structure (HTH conformation), similar to PDB entry 4DFX, which shows stronger coupling between the conformational dynamics observed at the A-helix and active site. The HTH conformation is stabilized by S10 phosphorylation of the kinase via ionic interactions between the protonated amine of K7 and the phosphate group on S10, further enhancing the dynamic coupling to the active site. These results support a role of N-myristoylation in the allosteric regulation of PKA-C