Integrating protein structural dynamics and evolutionary analysis with Bio3D

Abstract Background Popular 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. Results Here, 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. Conclusions The 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/ .http://deepblue.lib.umich.edu/bitstream/2027.42/109747/1/12859_2014_Article_399.pd

Neck linker docking is critical for Kinesin-1 force generation in cells but at a cost to motor speed and processivity.

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

Abdominal Computed Tomography Imaging Findings in Hospitalized COVID-19 Patients: A Year-Long Experience and Associations Revealed by Explainable Artificial Intelligence.

The aim of this retrospective study is to assess any association between abdominal CT findings and the radiological stage of COVID-19 pneumonia, pulmonary embolism and patient outcomes. We included 158 adult hospitalized COVID-19 patients between 1 March 2020 and 1 March 2021 who underwent 206 abdominal CTs. Two radiologists reviewed all CT images. Pathological findings were classified as acute or not. A subset of patients with inflammatory pathology in ACE2 organs (bowel, biliary tract, pancreas, urinary system) was identified. The radiological stage of COVID pneumonia, pulmonary embolism, overall days of hospitalization, ICU admission and outcome were registered. Univariate statistical analysis coupled with explainable artificial intelligence (AI) techniques were used to discover associations between variables. The most frequent acute findings were bowel abnormalities

Mapping the processivity determinants of the kinesin-3 motor domain

Kinesin superfamily members play important roles in many diverse cellular processes, including cell motility, cell division, intracellular transport, and regulation of the microtubule cytoskeleton. How the properties of the family-defining motor domain of distinct kinesins are tailored to their different cellular roles remains largely unknown. Here, we employed molecular-dynamics simulations coupled with energetic calculations to infer the family-specific interactions of kinesin-1 and kinesin-3 motor domains with microtubules in different nucleotide states. We then used experimental mutagenesis and single-molecule motility assays to further assess the predicted residue-wise determinants of distinct kinesin-microtubule binding properties. Collectively, our results identify residues in the L8, L11, and α6 regions that contribute to family-specific microtubule interactions and whose mutation affects motor-microtubule complex stability and processive motility (the ability of an individual motor to take multiple steps along its microtubule filament). In particular, substitutions of prominent kinesin-3 residues with those found in kinesin-1, namely, R167S/H171D, K266D, and R346M, were found to decrease kinesin-3 processivity 10-fold and thus approach kinesin-1 levels

Mapping the Structural and Dynamical Features of Kinesin Motor Domains

<div><p>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.</p></div

Results of loop13 G325A/G326A simulations.

<p>(<b>A</b>) Correlation of atomic displacement for all residue pairs for Eg5 mutant G325A/G326A. The color scale runs from pink (for values ranging between −1 to −0.5), through white (−0.5 to 0.5) to cyan (0.5 to 1). Negative values are indicative of displacements along opposite directions, namely anticorrelated motions, whereas positive values depict correlated motions occurring along the same direction. Major secondary structure elements are labeled and indicated schematically with helices in black and strands in gray. (<b>B</b>) Molecular representation of dynamic community partitioning for G325A/G326A mutant.</p

Simulations performed.

<p>Simulations were run in pairs, differing only in the assignment of initial velocities, to improve sampling <a href="http://www.ploscompbiol.org/article/info:doi/10.1371/journal.pcbi.1003329#pcbi.1003329-Gorfe1" target="_blank">[33]</a>, <a href="http://www.ploscompbiol.org/article/info:doi/10.1371/journal.pcbi.1003329#pcbi.1003329-Caves1" target="_blank">[49]</a>. RMSD values averaged over production phase dynamics for each simulation were calculated based on the Cα atom subset used for experimental structure PCA (see <a href="http://www.ploscompbiol.org/article/info:doi/10.1371/journal.pcbi.1003329#s3" target="_blank">Methods</a>).</p