The Bio3D packages for structural bioinformatics.

Bio3D is a family of R packages for the analysis of biomolecular sequence, structure, and dynamics. Major functionality includes biomolecular database searching and retrieval, sequence and structure conservation analysis, ensemble normal mode analysis, protein structure and correlation network analysis, principal component, and related multivariate analysis methods. Here, we review recent package developments, including a new underlying segregation into separate packages for distinct analysis, and introduce a new method for structure analysis named ensemble difference distance matrix analysis (eDDM).

Author Correction: A structural model for microtubule minus-end recognition and protection by CAMSAP proteins.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Comparative Protein Structure Analysis with Bio3D-Web.

Bio3D-web is an online application for the interactive analysis of sequence-structure-dynamics relationships in user-defined protein structure sets. Major functionality includes structure database searching, sequence and structure conservation assessment, inter-conformer relationship mapping and clustering with principal component analysis (PCA), and flexibility prediction and comparison with ensemble normal mode analysis (eNMA). Collectively these methods allow users to start with a single sequence or structure and characterize the structural, conformational, and internal dynamic properties of homologous proteins for which there are high-resolution structures available.

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.

Comparative structural dynamic analysis of GTPases.

GTPases regulate a multitude of essential cellular processes ranging from movement and division to differentiation and neuronal activity. These ubiquitous enzymes operate by hydrolyzing GTP to GDP with associated conformational changes that modulate affinity for family-specific binding partners. There are three major GTPase superfamilies: Ras-like GTPases, heterotrimeric G proteins and protein-synthesizing GTPases.

A posttranslational modification of the mitotic kinesin Eg5 that enhances its mechanochemical coupling and alters its mitotic function.

Numerous posttranslational modifications have been described in kinesins, but their consequences on motor mechanics are largely unknown. We investigated one of these-acetylation of lysine 146 in Eg5-by creating an acetylation mimetic lysine to glutamine substitution (K146Q). Lysine 146 is located in the α2 helix of the motor domain, where it makes an ionic bond with aspartate 91 on the neighboring α1 helix.

Altered chemomechanical coupling causes impaired motility of the kinesin-4 motors KIF27 and KIF7.

Kinesin-4 motors play important roles in cell division, microtubule organization, and signaling. Understanding how motors perform their functions requires an understanding of their mechanochemical and motility properties. We demonstrate that KIF27 can influence microtubule dynamics, suggesting a conserved function in microtubule organization across the kinesin-4 family.

A structural model for microtubule minus-end recognition and protection by CAMSAP proteins.

CAMSAP and Patronin family members regulate microtubule minus-end stability and localization and thus organize noncentrosomal microtubule networks, which are essential for cell division, polarization and differentiation. Here, we found that the CAMSAP C-terminal CKK domain is widely present among eukaryotes and autonomously recognizes microtubule minus ends. Through a combination of structural approaches, we uncovered how mammalian CKK binds between two tubulin dimers at the interprotofilament interface on the outer microtubule surface.

Navigating the conformational landscape of G protein-coupled receptor kinases during allosteric activation.

G protein-coupled receptors (GPCRs) are essential for transferring extracellular signals into carefully choreographed intracellular responses controlling diverse aspects of cell physiology. The duration of GPCR-mediated signaling is primarily regulated via GPCR kinase (GRK)-mediated phosphorylation of activated receptors. Although many GRK structures have been reported, the mechanisms underlying GRK activation are not well-understood, in part because it is unknown how these structures map to the conformational landscape available to this enzyme family.

Investigating Protein Sequence-structure-dynamics Relationships with Bio3D-web.

We demonstrate the usage of Bio3D-web for the interactive analysis of biomolecular structure data. The Bio3D-web application provides online functionality for: (1) The identification of related protein structure sets to user specified thresholds of similarity; (2) Their multiple alignment and structure superposition; (3) Sequence and structure conservation analysis; (4) Inter-conformer relationship mapping with principal component analysis, and (5) comparison of predicted internal dynamics via ensemble normal mode analysis. This integrated functionality provides a complete online workflow for investigating sequence-structure-dynamic relationships within protein families and superfamilies.

Structural and Molecular Mechanisms of Cytokine-Mediated Endocrine Resistance in Human Breast Cancer Cells.

Human breast cancers that exhibit high proportions of immune cells and elevated levels of pro-inflammatory cytokines predict poor prognosis. Here, we demonstrate that treatment of human MCF-7 breast cancer cells with pro-inflammatory cytokines results in ERα-dependent activation of gene expression and proliferation, in the absence of ligand or presence of 4OH-tamoxifen (TOT). Cytokine activation of ERα and endocrine resistance is dependent on phosphorylation of ERα at S305 in the hinge domain.

Computational and biochemical characterization of two partially overlapping interfaces and multiple weak-affinity K-Ras dimers.

Recent studies found that membrane-bound K-Ras dimers are important for biological function. However, the structure and thermodynamic stability of these complexes remained unknown because they are hard to probe by conventional approaches. Combining data from a wide range of computational and experimental approaches, here we describe the structure, dynamics, energetics and mechanism of assembly of multiple K-Ras dimers.

Systems Structural Biology Analysis of Ligand Effects on ERα Predicts Cellular Response to Environmental Estrogens and Anti-hormone Therapies.

Environmental estrogens and anti-hormone therapies for breast cancer have diverse tissue- and signaling-pathway-selective outcomes, but how estrogen receptor alpha (ERα) mediates this phenotypic diversity is poorly understood. We implemented a statistical approach to allow unbiased, parallel analyses of multiple crystal structures, and identified subtle perturbations of ERα structure by different synthetic and environmental estrogens. Many of these perturbations were in the sub-Å range, within the noise of the individual structures, but contributed significantly to the activities of synthetic and environmental estrogens.

Online interactive analysis of protein structure ensembles with Bio3D-web.

Unlabelled: Bio3D-web is an online application for analyzing the sequence, structure and conformational heterogeneity of protein families. Major functionality is provided for identifying protein structure sets for analysis, their alignment and refined structure superposition, sequence and structure conservation analysis, mapping and clustering of conformations and the quantitative comparison of their predicted structural dynamics.

Availability: Bio3D-web is based on the Bio3D and Shiny R packages.

Rapid Characterization of Allosteric Networks with Ensemble Normal Mode Analysis.

Allosteric regulation is a primary means of controlling protein function. By definition, allostery involves the propagation of structural dynamic changes between distal protein sites that yields a functional change. Gaining improved knowledge of these fundamental mechanisms is important for understanding many biomolecular processes and for guiding protein engineering and drug design efforts.

Dynamic Coupling and Allosteric Networks in the α Subunit of Heterotrimeric G Proteins.

G protein α subunits cycle between active and inactive conformations to regulate a multitude of intracellular signaling cascades. Important structural transitions occurring during this cycle have been characterized from extensive crystallographic studies. However, the link between observed conformations and the allosteric regulation of binding events at distal sites critical for signaling through G proteins remain unclear.

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.

Enrichment of druggable conformations from apo protein structures using cosolvent-accelerated molecular dynamics.

Here we describe the development of an improved workflow for utilizing experimental and simulated protein conformations in the structure-based design of inhibitors for anti-apoptotic Bcl-2 family proteins. Traditional structure-based approaches on similar targets are often constrained by the sparsity of available structures and difficulties in finding lead compounds that dock against flat, flexible protein-protein interaction surfaces. By employing computational docking of known small molecule inhibitors, we have demonstrated that structural ensembles derived from either accelerated MD (aMD) or MD in the presence of an organic cosolvent generally give better scores than those assessed from analogous conventional MD.

Integrating protein structural dynamics and evolutionary analysis with Bio3D.

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.

Kinesin-5 allosteric inhibitors uncouple the dynamics of nucleotide, microtubule, and neck-linker binding sites.

Kinesin motor domains couple cycles of ATP hydrolysis to cycles of microtubule binding and conformational changes that result in directional force and movement on microtubules. The general principles of this mechanochemical coupling have been established; however, fundamental atomistic details of the underlying allosteric mechanisms remain unknown. This lack of knowledge hampers the development of new inhibitors and limits our understanding of how disease-associated mutations in distal sites can interfere with the fidelity of motor domain function.

Mapping the structural and dynamical features of kinesin motor domains.

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.

Domain-opening and dynamic coupling in the α-subunit of heterotrimeric G proteins.

Heterotrimeric G proteins are conformational switches that turn on intracellular signaling cascades in response to the activation of G-protein-coupled receptors. Receptor activation by extracellular stimuli promotes a cycle of GTP binding and hydrolysis on the G protein α-subunit (Gα). Important conformational transitions occurring during this cycle have been characterized from extensive crystallographic studies of Gα.

LeuT conformational sampling utilizing accelerated molecular dynamics and principal component analysis.

Monoamine transporters (MATs) function by coupling ion gradients to the transport of dopamine, norepinephrine, or serotonin. Despite their importance in regulating neurotransmission, the exact conformational mechanism by which MATs function remains elusive. To this end, we have performed seven 250 ns accelerated molecular dynamics simulations of the leucine transporter, a model for neurotransmitter MATs.

Calcium binding and allosteric signaling mechanisms for the sarcoplasmic reticulum Ca²+ ATPase.

The sarcoplasmic reticulum Ca²⁺ ATPase (SERCA) is a membrane-bound pump that utilizes ATP to drive calcium ions from the myocyte cytosol against the higher calcium concentration in the sarcoplasmic reticulum. Conformational transitions associated with Ca²⁺-binding are important to its catalytic function. We have identified collective motions that partition SERCA crystallographic structures into multiple catalytically-distinct states using principal component analysis.

Identification of potential small molecule binding pockets on Rho family GTPases.

Rho GTPases are conformational switches that control a wide variety of signaling pathways critical for eukaryotic cell development and proliferation. They represent attractive targets for drug design as their aberrant function and deregulated activity is associated with many human diseases including cancer. Extensive high-resolution structures (>100) and recent mutagenesis studies have laid the foundation for the design of new structure-based chemotherapeutic strategies.