It has become commonplace to employ principal component analysis to reveal the most important motions in proteins. This method is more commonly known by its acronym, PCA. While most popular molecular dynamics packages inevitably provide PCA tools to analyze protein trajectories, researchers often make inferences of their results without having insight into how to make interpretations, and they are often unaware of limitations and generalizations of such analysis. Here we review best practices for applying standard PCA, describe useful variants, discuss why one may wish to make comparison studies, and describe a set of metrics that make comparisons possible. In practice, one will be forced to make inferences about the essential dynamics of a protein without having the desired amount of samples. Therefore, considerable time is spent on describing how to judge the significance of results, highlighting pitfalls. The topic of PCA is reviewed from the perspective of many practical considerations, and useful recipes are provided.

We used FRET to examine the kinetics and thermodynamics of structural changes associated with ADP release in myosin V, which is thought to be a strain sensitive step in many muscle and non-muscle myosins. We also explore essential dynamics using FIRST/FRODA starting with three different myosin V X-ray crystal structures to examine intrinsic flexibility and correlated motions. Our steady-state and time resolved FRET analysis demonstrates a temperature dependent reversible conformational change in the nucleotide binding pocket. Our kinetic results demonstrate that the nucleotide binding pocket goes from a closed to an open conformation prior to the release of ADP while the actin binding cleft remains closed. Interestingly, we find that the temperature dependence of the maximum actin-activated myosin V ATPase rate is similar to the pocket opening step, demonstrating this is the rate limiting structural transition in the ATPase cycle. Thermodynamic analysis demonstrates the transition from the open to closed nucleotide binding pocket conformation is unfavorable because of a decrease in entropy. The intrinsic flexibility analysis is consistent with conformational entropy playing a role in this transition as the MV.ADP structure is highly flexible compared to the MV.APO structure. Our experimental and modeling studies support the conclusion of a novel post-power-stroke actomyosin.ADP state in which the nucleotide binding pocket and actin binding cleft are closed. The novel state may be important for strain sensitivity as the transition from the closed to open nucleotide binding pocket conformation may be altered by lever arm position.

Conserved active-site elements in myosins and other P-loop NTPases play critical roles in nucleotide binding and hydrolysis; however, the mechanisms of allosteric communication among these mechanoenzymes remain unresolved. In this work we introduced the E442A mutation, which abrogates a salt-bridge between switch I and switch II, and the G440A mutation, which abolishes a main-chain hydrogen bond associated with the interaction of switch II with the γ phosphate of ATP, into myosin V. We used fluorescence resonance energy transfer between mant-labeled nucleotides or IAEDANS-labeled actin and FlAsH-labeled myosin V to examine the conformation of the nucleotide- and actin-binding regions, respectively. We demonstrate that in the absence of actin, both the G440A and E442A mutants bind ATP with similar affinity and result in only minor alterations in the conformation of the nucleotide-binding pocket (NBP). In the presence of ADP and actin, both switch II mutants disrupt the formation of a closed NBP actomyosin.ADP state. The G440A mutant also prevents ATP-induced opening of the actin-binding cleft. Our results indicate that the switch II region is critical for stabilizing the closed NBP conformation in the presence of actin, and is essential for communication between the active site and actin-binding region.

To clarify the extent structure plays in determining protein dynamics, a comparative study is made using three models that characterize native state dynamics of single domain proteins starting from known structures taken from four distinct SCOP classifications. A geometrical simulation using the Framework Rigidity Optimized Dynamics Algorithm (FRODA) based on rigid cluster decomposition is compared to the commonly employed elastic network model (specifically the Anisotropic Network Model ANM) and molecular dynamics (MD) simulation. The essential dynamics are quantified by a mode subspace constructed from ANM and a principal component analysis (PCA) on FRODA and MD trajectories. Aggregate conformational ensembles are constructed to provide a basis for quantitative comparisons between FRODA runs using different parameter settings to critically assess how the predictions of essential dynamics depend on a priori arbitrary user-defined distance constraint rules. We established a range of physicality for these parameters. Surprisingly, FRODA maintains greater intra-consistent results than obtained from MD trajectories, comparable to ANM. Additionally, a mode subspace is constructed from PCA on an exemplar set of myoglobin structures from the Protein Data Bank. Significant overlap across the three model subspaces and the experimentally derived subspace is found. While FRODA provides the most robust sampling and characterization of the native basin, all three models give similar dynamical information of a native state, further demonstrating that structure is the key determinant of dynamics.