BackgroundPopular bioinformatics approaches for studying protein functional dynamics include comparisons of crystallographic structures, molecular dynamics simulations and normal mode analysis. However, determining how observed displacements and predicted motions from these traditionally separate analyses relate to each other, as well as to the evolution of sequence, structure and function within large protein families, remains a considerable challenge. This is in part due to the general lack of tools that integrate information of molecular structure, dynamics and evolution.ResultsHere, we describe the integration of new methodologies for evolutionary sequence, structure and simulation analysis into the Bio3D package. This major update includes unique high-throughput normal mode analysis for examining and contrasting the dynamics of related proteins with non-identical sequences and structures, as well as new methods for quantifying dynamical couplings and their residue-wise dissection from correlation network analysis. These new methodologies are integrated with major biomolecular databases as well as established methods for evolutionary sequence and comparative structural analysis. New functionality for directly comparing results derived from normal modes, molecular dynamics and principal component analysis of heterogeneous experimental structure distributions is also included. We demonstrate these integrated capabilities with example applications to dihydrofolate reductase and heterotrimeric G-protein families along with a discussion of the mechanistic insight provided in each case.ConclusionsThe integration of structural dynamics and evolutionary analysis in Bio3D enables researchers to go beyond a prediction of single protein dynamics to investigate dynamical features across large protein families. The Bio3D package is distributed with full source code and extensive documentation as a platform independent R package under a GPL2 license from http://thegrantlab.org/bio3d/.Electronic supplementary materialThe online version of this article (doi:10.1186/s12859-014-0399-6) contains supplementary material, which is available to authorized users.

CAMSAP and Patronin family members regulate microtubule minus-end stability and localization and thus organize non-centrosomal microtubule networks, which are essential for cell division, polarization and differentiation. Here, we show that the C-terminal CKK domain of CAMSAPs is widely present among eukaryotes and autonomously recognizes microtubule minus ends. Through a combination of structural approaches, we uncover how mammalian CKK binds between two tubulin dimers at the inter-protofilament interface on the outer microtubule surface. In vitro reconstitution assays combined with high resolution fluorescence microscopy and cryo-electron tomography suggest that CKK preferentially associates with the transition zone between curved protofilaments and the regular microtubule lattice. We propose that minus-end-specific features of the inter-protofilament interface at this site form the basis for CKK’s minus-end preference. The steric clash between microtubule-bound CKK and kinesin motors explains how CKK protects microtubule minus ends against kinesin-13-induced depolymerization and thus controls the stability of free microtubule minus ends.

An increasing number of functional studies of proteins have shown that sequence and structural similarities alone may not be sufficient for reliable prediction of their interaction properties. This is particularly true for proteins recognizing specific antibodies, where the prediction of antibody-binding sites, called epitopes, has proven challenging. The antibody-binding properties of an antigen depend on its structure and related dynamics. Aiming to predict the antibody-binding regions of a protein, we investigate a new approach based on the integrated analysis of the dynamical and energetic properties of antigens, to identify nonoptimized, low-intensity energetic interaction networks in the protein structure isolated in solution. The method is based on the idea that recognition sites may correspond to localized regions with low-intensity energetic couplings with the rest of the protein, which allows them to undergo conformational changes, to be recognized by a binding partner, and to tolerate mutations with minimal energetic expense. Upon analyzing the results on isolated proteins and benchmarking against antibody complexes, it is found that the method successfully identifies binding sites located on the protein surface that are accessible to putative binding partners. The combination of dynamics and energetics can thus discriminate between epitopes and other substructures based only on physical properties. We discuss implications for vaccine design.

Kinesin force generation involves ATP-induced docking of the neck linker (NL) along the motor core. However, the roles of the proposed steps of NL docking, cover-neck bundle (CNB) and asparagine latch (N-latch) formation, during force generation are unclear. Furthermore, the necessity of NL docking for transport of membrane-bound cargo in cells has not been tested. We generated kinesin-1 motors impaired in CNB and/or N-latch formation based on molecular dynamics simulations. The mutant motors displayed reduced force output and inability to stall in optical trap assays but exhibited increased speeds, run lengths, and landing rates under unloaded conditions. NL docking thus enhances force production but at a cost to speed and processivity. In cells, teams of mutant motors were hindered in their ability to drive transport of Golgi elements (high-load cargo) but not peroxisomes (low-load cargo). These results demonstrate that the NL serves as a mechanical element for kinesin-1 transport under physiological conditions.

Kinesin motor proteins drive intracellular transport by coupling ATP hydrolysis to conformational changes that mediate directed movement along microtubules. Characterizing these distinct conformations and their interconversion mechanism is essential to determining an atomic-level model of kinesin action. Here we report a comprehensive principal component analysis of 114 experimental structures along with the results of conventional and accelerated molecular dynamics simulations that together map the structural dynamics of the kinesin motor domain. All experimental structures were found to reside in one of three distinct conformational clusters (ATP-like, ADP-like and Eg5 inhibitor-bound). These groups differ in the orientation of key functional elements, most notably the microtubule binding α4–α5, loop8 subdomain and α2b-β4-β6-β7 motor domain tip. Group membership was found not to correlate with the nature of the bound nucleotide in a given structure. However, groupings were coincident with distinct neck-linker orientations. Accelerated molecular dynamics simulations of ATP, ADP and nucleotide free Eg5 indicate that all three nucleotide states could sample the major crystallographically observed conformations. Differences in the dynamic coupling of distal sites were also evident. In multiple ATP bound simulations, the neck-linker, loop8 and the α4–α5 subdomain display correlated motions that are absent in ADP bound simulations. Further dissection of these couplings provides evidence for a network of dynamic communication between the active site, microtubule-binding interface and neck-linker via loop7 and loop13. Additional simulations indicate that the mutations G325A and G326A in loop13 reduce the flexibility of these regions and disrupt their couplings. Our combined results indicate that the reported ATP and ADP-like conformations of kinesin are intrinsically accessible regardless of nucleotide state and support a model where neck-linker docking leads to a tighter coupling of the microtubule and nucleotide binding regions. Furthermore, simulations highlight sites critical for large-scale conformational changes and the allosteric coupling between distal functional sites.