iGNM 2.0: the Gaussian network model database for biomolecular structural dynamics

Ensemble-based modeling and rigidity decomposition of allosteric interaction networks and communication pathways in cyclin-dependent kinases: Differentiating kinase clients of the Hsp90-Cdc37 chaperone

The overarching goal of delineating molecular principles underlying differentiation of protein kinase clients and chaperone-based modulation of kinase activity is fundamental to understanding activity of many oncogenic kinases that require chaperoning of Hsp70 and Hsp90 systems to attain a functionally competent active form. Despite structural similarities and common activation mechanisms shared by cyclin-dependent kinase (CDK) proteins, members of this family can exhibit vastly different chaperone preferences. The molecular determinants underlying chaperone dependencies of protein kinases are not fully understood as structurally similar kinases may often elicit distinct regulatory responses to the chaperone. The regulatory divergences observed for members of CDK family are of particular interest as functional diversification among these kinases may be related to variations in chaperone dependencies and can be exploited in drug discovery of personalized therapeutic agents. In this work, we report the results of a computational investigation of several members of CDK family (CDK5, CDK6, CDK9) that represented a broad repertoire of chaperone dependencies—from nonclient CDK5, to weak client CDK6, and strong client CDK9. By using molecular simulations of multiple crystal structures we characterized conformational ensembles and collective dynamics of CDK proteins. We found that the elevated dynamics of CDK9 can trigger imbalances in cooperative collective motions and reduce stability of the active fold, thus creating a cascade of favorable conditions for chaperone intervention. The ensemble-based modeling of residue interaction networks and community analysis determined how differences in modularity of allosteric networks and topography of communication pathways can be linked with the client status of CDK proteins. This analysis unveiled depleted modularity of the allosteric network in CDK9 that alters distribution of communication pathways and leads to impaired signaling in the client kinase. According to our results, these network features may uniquely define chaperone dependencies of CDK clients. The perturbation response scanning and rigidity decomposition approaches identified regulatory hotspots that mediate differences in stability and cooperativity of allosteric interaction networks in the CDK structures. By combining these synergistic approaches, our study revealed dynamic and network signatures that can differentiate kinase clients and rationalize subtle divergences in the activation mechanisms of CDK family members. The therapeutic implications of these results are illustrated by identifying structural hotspots of pathogenic mutations that preferentially target regions of the increased flexibility to enable modulation of activation changes. Our study offers a network-based perspective on dynamic kinase mechanisms and drug design by unravelling relationships between protein kinase dynamics, allosteric communications and chaperone dependencies.

PEBP1 Wardens Ferroptosis by Enabling Lipoxygenase Generation of Lipid Death Signals

Ferroptosis is a form of programmed cell death that is pathogenic to several acute and chronic diseases and executed via oxygenation of polyunsaturated phosphatidylethanolamines (PE) by 15-lipoxygenases (15-LO) that normally use free polyunsaturated fatty acids as substrates. Mechanisms of the altered 15-LO substrate specificity are enigmatic. We sought a common ferroptosis regulator for 15LO. We discovered that PEBP1, a scaffold protein inhibitor of protein kinase cascades, complexes with two 15LO isoforms, 15LO1 and 15LO2, and changes their substrate competence to generate hydroperoxy-PE. Inadequate reduction of hydroperoxy-PE due to insufficiency or dysfunction of a selenoperoxidase, GPX4, leads to ferroptosis. We demonstrated the importance of PEBP1-dependent regulatory mechanisms of ferroptotic death in airway epithelial cells in asthma, kidney epithelial cells in renal failure, and cortical and hippocampal neurons in brain trauma. As master regulators of ferroptotic cell death with profound implications for human disease, PEBP1/15LO complexes represent a new target for drug discovery.

Gaussian network model can be enhanced by combining solvent accessibility in proteins

Gaussian network model (GNM), regarded as the simplest and most representative coarse-grained model, has been widely adopted to analyze and reveal protein dynamics and functions. Designing a variation of the classical GNM, by defining a new Kirchhoff matrix, is the way to improve the residue flexibility modeling. We combined information arising from local relative solvent accessibility (RSA) between two residues into the Kirchhoff matrix of the parameter-free GNM. The undetermined parameters in the new Kirchhoff matrix were estimated by using particle swarm optimization. The usage of RSA was motivated by the fact that our previous work using RSA based linear regression model resulted out higher prediction quality of the residue flexibility when compared with the classical GNM and the parameter free GNM. Computational experiments, conducted based on one training dataset, two independent datasets and one additional small set derived by molecular dynamics simulations, demonstrated that the average correlation coefficients of the proposed RSA based parameter-free GNM, called RpfGNM, were significantly increased when compared with the parameter-free GNM. Our empirical results indicated that a variation of the classical GNMs by combining other protein structural properties is an attractive way to improve the quality of flexibility modeling.

Dynamic Modulation of Binding Affinity as a Mechanism for Regulating Interferon Signaling

How structural dynamics affects cytokine signaling is under debate. Here, we investigated the dynamics of the type I interferon (IFN) receptor, IFNAR1, and its effect on signaling upon binding IFN and IFNAR2 using a combination of structure-based mechanistic studies, in situ binding, and gene induction assays. Our study reveals that IFNAR1 flexibility modulates ligand-binding affinity, which, in turn, regulates biological signaling. We identified the hinge sites and key interactions implicated in IFNAR1 inter-subdomain (SD1-SD4) movements. We showed that the predicted cooperative movements are essential to accommodate intermolecular interactions. Engineered disulfide bridges, computationally predicted to interfere with IFNAR1 dynamics, were experimentally confirmed. Notably, introducing disulfide bonds between subdomains SD2 and SD3 modulated IFN binding and activity in accordance with the relative attenuation of cooperative movements with varying distance from the hinge center, whereas locking the SD3-SD4 interface flexibility in favor of an extended conformer increased activity.

Chromosomal dynamics predicted by an elastic network model explains genome-wide accessibility and long-range couplings

Understanding the three-dimensional (3D) architecture of chromatin and its relation to gene expression and regulation is fundamental to understanding how the genome functions. Advances in Hi-C technology now permit us to study 3D genome organization, but we still lack an understanding of the structural dynamics of chromosomes. The dynamic couplings between regions separated by large genomic distances (>50 Mb) have yet to be characterized. We adapted a well-established protein-modeling framework, the Gaussian Network Model (GNM), to model chromatin dynamics using Hi-C data. We show that the GNM can identify spatial couplings at multiple scales: it can quantify the correlated fluctuations in the positions of gene loci, find large genomic compartments and smaller topologically-associating domains (TADs) that undergo en bloc movements, and identify dynamically coupled distal regions along the chromosomes. We show that the predictions of the GNM correlate well with genome-wide experimental measurements. We use the GNM to identify novel cross-correlated distal domains (CCDDs) representing pairs of regions distinguished by their long-range dynamic coupling and show that CCDDs are associated with increased gene co-expression. Together, these results show that GNM provides a mathematically well-founded unified framework for modeling chromatin dynamics and assessing the structural basis of genome-wide observations.

Intrinsic Dynamics Analysis of a DNA Octahedron by Elastic Network Model

DNA is a fundamental component of living systems where it plays a crucial role at both functional and structural level. The programmable properties of DNA make it an interesting building block for the construction of nanostructures. However, molecular mechanisms for the arrangement of these well-defined DNA assemblies are not fully understood. In this paper, the intrinsic dynamics of a DNA octahedron has been investigated by using two types of Elastic Network Models (ENMs). The application of ENMs to DNA nanocages include the analysis of the intrinsic flexibilities of DNA double-helices and hinge sites through the calculation of the square fluctuations, as well as the intrinsic collective dynamics in terms of cross-collective map calculation coupled with global motions analysis. The dynamics profiles derived from ENMs have then been evaluated and compared with previous classical molecular dynamics simulation trajectories. The results presented here revealed that ENMs can provide useful insights into the intrinsic dynamics of large DNA nanocages and represent a useful tool in the field of structural DNA nanotechnology.

The 2016 database issue of Nucleic Acids Research and an updated molecular biology database collection

The 2016 Database Issue of Nucleic Acids Research starts with overviews of the resources provided by three major bioinformatics centers, the U.S. National Center for Biotechnology Information (NCBI), the European Bioinformatics Institute (EMBL-EBI) and Swiss Institute for Bioinformatics (SIB). Also included are descriptions of 62 new databases and updates on 95 databases that have been previously featured in NAR plus 17 previously described elsewhere. A number of papers in this issue deal with resources on nucleic acids, including various kinds of non-coding RNAs and their interactions, molecular dynamics simulations of nucleic acid structure, and two databases of super-enhancers. The protein database section features important updates on the EBI's Pfam, PDBe and PRIDE databases, as well as a variety of resources on pathways, metabolomics and metabolic modeling. This issue also includes updates on popular metagenomics resources, such as MG-RAST, EBI Metagenomics, and probeBASE, as well as a newly compiled Human Pan-Microbe Communities database. A significant fraction of the new and updated databases are dedicated to the genetic basis of disease, primarily cancer, and various aspects of drug research, including resources for patented drugs, their side effects, withdrawn drugs, and potential drug targets. A further six papers present updated databases of various antimicrobial and anticancer peptides. The entire Database Issue is freely available online on the Nucleic Acids Research website (http://nar.oxfordjournals.org/). The NAR online Molecular Biology Database Collection, http://www.oxfordjournals.org/nar/database/c/, has been updated with the addition of 88 new resources and removal of 23 obsolete websites, which brought the current listing to 1685 databases.