The evolution of protein conformational dynamics contains important information about protein function and regulation. Here, we describe an approach to dynamical-evolution analysis based on multiple microsecond molecular dynamics simulations and residue− residue contact analysis. We illustrate our approach by comparing three human cyclophilin isoforms, cyclophilin A, D, and E, which belong to a family of enzymes catalyzing peptidyl-prolyl cis−trans isomerization. Our results reveal that despite distinct overall equilibrium conformations between cyclophilins under substrate-f ree conditions, functional dynamical changes resembling substrate-binding and catalytic processes tend to be conserved. Key residues displaying either concerted or specific dynamical changes among isoforms during the reactions are identified, which delineate two distinct allosteric pathways for cyclophilin function consistent with recent nuclear magnetic resonance experiments. A sequence-based coevolution analysis is also employed for further understanding dynamical consequences. Our results collectively provide a framework where both common and specific functional mechanisms of a protein family can be elucidated.

The local conformation of DNA plays a critical role in the recognition-andbinding process of transcription factors. In particular, the most common form of the DNA double helix, B-DNA, exists as a conformational equilibrium between BI and BII substates that are associated with changes in major and minor groove dimensions. It has been hypothesized that cytosine methylation, the most common epigenetic modification and one that is often found in so-called CpG islands, modulates transcription factor binding affinity by altering the BI-BII equilibrium in the neighborhood of methylation sites. Indeed, previous molecular dynamics studies have suggested that cytosine methylation leads to a uniform stabilization of the BI substate in CpG islands. In this work we first benchmark the ability of the latest Amber DNA force fields and a variety of three-point water models to reproduce BI-BII equilibrium in the Dickerson dodecamer and reduce terminal base pair fraying. Using the best-performing combinations, we then re-examine the conformational effects of cytosine methylation in a prototypical CpG island, (GC) 5 . We find that the BI substate is stabilized for GpC steps in such sequences, but that the BII substate is stabilized for the CpG steps. More interestingly, the Amber ff99bsc0 DNA force field used in previous studies yields opposite conclusions. These findings highlight the need for continued fixed-charge force field development for nucleic acids, as well as solution state experimental datasets involving epigenetic modification that can serve as benchmarks for simulation accuracy.