Lipid Modulation of a Class B GPCR: Elucidating the Modulatory Role of PI(4,5)P2 Lipids.

Phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2) lipids have been shown to stabilize an active conformation of class A G-protein coupled receptors (GPCRs) through a conserved binding site, not present in class B GPCRs. For class B GPCRs, previous molecular dynamics (MD) simulation studies have shown PI(4,5)P2 interacting with the Glucagon receptor (GCGR), which constitutes an important target for diabetes and obesity therapeutics. In this work, we applied MD simulations supported by native mass spectrometry (nMS) to study lipid interactions with GCGR. We demonstrate how tail composition plays a role in modulating the binding of PI(4,5)P2 lipids to GCGR. Specifically, we find the PI(4,5)P2 lipids to have a higher affinity toward the inactive conformation of GCGR. Interestingly, we find that in contrast to class A GPCRs, PI(4,5)P2 appear to stabilize the inactive conformation of GCGR through a binding site conserved across class B GPCRs but absent in class A GPCRs. This suggests differences in the regulatory function of PI(4,5)P2 between class A and class B GPCRs.

General Protocol for Constructing Molecular Models of Nanodiscs.

Nanodisc technology is increasingly being applied for structural and biophysical studies of membrane proteins. In this work, we present a general protocol for constructing molecular models of nanodiscs for molecular dynamics simulations. The protocol is written in python and based on geometric equations, making it fast and easy to modify, enabling automation and customization of nanodiscs in silico. The novelty being the ability to construct any membrane scaffold protein (MSP) variant fast and easy given only an input sequence. We validated and tested the protocol by simulating seven different nanodiscs of various sizes and with different membrane scaffold proteins, both circularized and noncircularized. The structural and biophysical properties were analyzed and shown to be in good agreement with previously reported experimental data and simulation studies.

Deep neural networks identify signaling mechanisms of ErbB-family drug resistance from a continuous cell morphology space.

It is well known that the development of drug resistance in cancer cells can lead to changes in cell morphology. Here, we describe the use of deep neural networks to analyze this relationship, demonstrating that complex cell morphologies can encode states of signaling networks and unravel cellular mechanisms hidden to conventional approaches. We perform high-content screening of 17 cancer cell lines, generating more than 500 billion data points from ∼850 million cells. We analyze these data using a deep learning model, resulting in the identification of a continuous 27-dimension space describing all of the observed cell morphologies. From its morphology alone, we could thus predict whether a cell was resistant to ErbB-family drugs, with an accuracy of 74%, and predict the potential mechanism of resistance, subsequently validating the role of MET and insulin-like growth factor 1 receptor (IGF1R) as drivers of cetuximab resistance in in vitro models of lung and head/neck cancer.

Bowhead: Bayesian modelling of cell velocity during concerted cell migration.

Cell migration is a central biological process that requires fine coordination of molecular events in time and space. A deregulation of the migratory phenotype is also associated with pathological conditions including cancer where cell motility has a causal role in tumor spreading and metastasis formation. Thus cell migration is of critical and strategic importance across the complex disease spectrum as well as for the basic understanding of cell phenotype. Experimental studies of the migration of cells in monolayers are often conducted with 'wound healing' assays. Analysis of these assays has traditionally relied on how the wound area changes over time. However this method does not take into account the shape of the wound. Given the many options for creating a wound healing assay and the fact that wound shape invariably changes as cells migrate this is a significant flaw. Here we present a novel software package for analyzing concerted cell velocity in wound healing assays. Our method encompasses a wound detection algorithm based on cell confluency thresholding and employs a Bayesian approach in order to estimate concerted cell velocity with an associated likelihood. We have applied this method to study the effect of siRNA knockdown on the migration of a breast cancer cell line and demonstrate that cell velocity can track wound healing independently of wound shape and provides a more robust quantification with significantly higher signal to noise ratios than conventional analyses of wound area. The software presented here will enable other researchers in any field of cell biology to quantitatively analyze and track live cell migratory processes and is therefore expected to have a significant impact on the study of cell migration, including cancer relevant processes. Installation instructions, documentation and source code can be found at http://bowhead.lindinglab.science licensed under GPLv3.

Dynamic Rearrangement of Cell States Detected by Systematic Screening of Sequential Anticancer Treatments.

Signaling networks are nonlinear and complex, involving a large ensemble of dynamic interaction states that fluctuate in space and time. However, therapeutic strategies, such as combination chemotherapy, rarely consider the timing of drug perturbations. If we are to advance drug discovery for complex diseases, it will be essential to develop methods capable of identifying dynamic cellular responses to clinically relevant perturbations. Here, we present a Bayesian dose-response framework and the screening of an oncological drug matrix, comprising 10,000 drug combinations in melanoma and pancreatic cancer cell lines, from which we predict sequentially effective drug combinations. Approximately 23% of the tested combinations showed high-confidence sequential effects (either synergistic or antagonistic), demonstrating that cellular perturbations of many drug combinations have temporal aspects, which are currently both underutilized and poorly understood.

Computational Redesign of Thioredoxin Is Hypersensitive toward Minor Conformational Changes in the Backbone Template.

Despite the development of powerful computational tools, the full-sequence design of proteins still remains a challenging task. To investigate the limits and capabilities of computational tools, we conducted a study of the ability of the program Rosetta to predict sequences that recreate the authentic fold of thioredoxin. Focusing on the influence of conformational details in the template structures, we based our study on 8 experimentally determined template structures and generated 120 designs from each. For experimental evaluation, we chose six sequences from each of the eight templates by objective criteria. The 48 selected sequences were evaluated based on their progressive ability to (1) produce soluble protein in Escherichia coli and (2) yield stable monomeric protein, and (3) on the ability of the stable, soluble proteins to adopt the target fold. Of the 48 designs, we were able to synthesize 32, 20 of which resulted in soluble protein. Of these, only two were sufficiently stable to be purified. An X-ray crystal structure was solved for one of the designs, revealing a close resemblance to the target structure. We found a significant difference among the eight template structures to realize the above three criteria despite their high structural similarity. Thus, in order to improve the success rate of computational full-sequence design methods, we recommend that multiple template structures are used. Furthermore, this study shows that special care should be taken when optimizing the geometry of a structure prior to computational design when using a method that is based on rigid conformations.

Kinome-wide decoding of network-attacking mutations rewiring cancer signaling.

Cancer cells acquire pathological phenotypes through accumulation of mutations that perturb signaling networks. However, global analysis of these events is currently limited. Here, we identify six types of network-attacking mutations (NAMs), including changes in kinase and SH2 modulation, network rewiring, and the genesis and extinction of phosphorylation sites. We developed a computational platform (ReKINect) to identify NAMs and systematically interpreted the exomes and quantitative (phospho-)proteomes of five ovarian cancer cell lines and the global cancer genome repository. We identified and experimentally validated several NAMs, including PKCγ M501I and PKD1 D665N, which encode specificity switches analogous to the appearance of kinases de novo within the kinome. We discover mutant molecular logic gates, a drift toward phospho-threonine signaling, weakening of phosphorylation motifs, and kinase-inactivating hotspots in cancer. Our method pinpoints functional NAMs, scales with the complexity of cancer genomes and cell signaling, and may enhance our capability to therapeutically target tumor-specific networks.

Microbial biofilm as a smart material.

Microbial biofilm colonies will in many cases form a smart material capable of responding to external threats dependent on their size and internal state. The microbial community accordingly switches between passive, protective, or attack modes of action. In order to decide which strategy to employ, it is essential for the biofilm community to be able to sense its own size. The sensor designed to perform this task is termed a quorum sensor, since it only permits collective behaviour once a sufficiently large assembly of microbes have been established. The generic quorum sensor construct involves two genes, one coding for the production of a diffusible signal molecule and one coding for a regulator protein dedicated to sensing the signal molecules. A positive feedback in the signal molecule production sets a well-defined condition for switching into the collective mode. The activation of the regulator involves a slow dimerization, which allows low-pass filtering of the activation of the collective mode. Here, we review and combine the model components that form the basic quorum sensor in a number of Gram-negative bacteria, e.g., Pseudomonas aeruginosa.

Equilibrium simulations of proteins using molecular fragment replacement and NMR chemical shifts.

Methods of protein structure determination based on NMR chemical shifts are becoming increasingly common. The most widely used approaches adopt the molecular fragment replacement strategy, in which structural fragments are repeatedly reassembled into different complete conformations in molecular simulations. Although these approaches are effective in generating individual structures consistent with the chemical shift data, they do not enable the sampling of the conformational space of proteins with correct statistical weights. Here, we present a method of molecular fragment replacement that makes it possible to perform equilibrium simulations of proteins, and hence to determine their free energy landscapes. This strategy is based on the encoding of the chemical shift information in a probabilistic model in Markov chain Monte Carlo simulations. First, we demonstrate that with this approach it is possible to fold proteins to their native states starting from extended structures. Second, we show that the method satisfies the detailed balance condition and hence it can be used to carry out an equilibrium sampling from the Boltzmann distribution corresponding to the force field used in the simulations. Third, by comparing the results of simulations carried out with and without chemical shift restraints we describe quantitatively the effects that these restraints have on the free energy landscapes of proteins. Taken together, these results demonstrate that the molecular fragment replacement strategy can be used in combination with chemical shift information to characterize not only the native structures of proteins but also their conformational fluctuations.

Combining experiments and simulations using the maximum entropy principle.

A key component of computational biology is to compare the results of computer modelling with experimental measurements. Despite substantial progress in the models and algorithms used in many areas of computational biology, such comparisons sometimes reveal that the computations are not in quantitative agreement with experimental data. The principle of maximum entropy is a general procedure for constructing probability distributions in the light of new data, making it a natural tool in cases when an initial model provides results that are at odds with experiments. The number of maximum entropy applications in our field has grown steadily in recent years, in areas as diverse as sequence analysis, structural modelling, and neurobiology. In this Perspectives article, we give a broad introduction to the method, in an attempt to encourage its further adoption. The general procedure is explained in the context of a simple example, after which we proceed with a real-world application in the field of molecular simulations, where the maximum entropy procedure has recently provided new insight. Given the limited accuracy of force fields, macromolecular simulations sometimes produce results that are at not in complete and quantitative accordance with experiments. A common solution to this problem is to explicitly ensure agreement between the two by perturbing the potential energy function towards the experimental data. So far, a general consensus for how such perturbations should be implemented has been lacking. Three very recent papers have explored this problem using the maximum entropy approach, providing both new theoretical and practical insights to the problem. We highlight each of these contributions in turn and conclude with a discussion on remaining challenges.

Robust Estimation of Diffusion-Optimized Ensembles for Enhanced Sampling.

The multicanonical, or flat-histogram, method is a common technique to improve the sampling efficiency of molecular simulations. The idea is that free-energy barriers in a simulation can be removed by simulating from a distribution where all values of a reaction coordinate are equally likely, and subsequently reweight the obtained statistics to recover the Boltzmann distribution at the temperature of interest. While this method has been successful in practice, the choice of a flat distribution is not necessarily optimal. Recently, it was proposed that additional performance gains could be obtained by taking the position-dependent diffusion coefficient into account, thus placing greater emphasis on regions diffusing slowly. Although some promising examples of applications of this approach exist, the practical usefulness of the method has been hindered by the difficulty in obtaining sufficiently accurate estimates of the diffusion coefficient. Here, we present a simple, yet robust solution to this problem. Compared to current state-of-the-art procedures, the new estimation method requires an order of magnitude fewer data to obtain reliable estimates, thus broadening the potential scope in which this technique can be applied in practice.

Probabilistic Determination of Native State Ensembles of Proteins.

The motions of biological macromolecules are tightly coupled to their functions. However, while the study of fast motions has become increasingly feasible in recent years, the study of slower, biologically important motions remains difficult. Here, we present a method to construct native state ensembles of proteins by the combination of physical force fields and experimental data through modern statistical methodology. As an example, we use NMR residual dipolar couplings to determine a native state ensemble of the extensively studied third immunoglobulin binding domain of protein G (GB3). The ensemble accurately describes both local and nonlocal backbone fluctuations as judged by its reproduction of complementary experimental data. While it is difficult to assess precise time-scales of the observed motions, our results suggest that it is possible to construct realistic conformational ensembles of biomolecules very efficiently. The approach may allow for a dramatic reduction in the computational as well as experimental resources needed to obtain accurate conformational ensembles of biological macromolecules in a statistically sound manner.

Size of quorum sensing communities.

Ensembles of bacteria are able to coordinate their phenotypic behavior in accordance with the size, density, and growth state of the ensemble. This is achieved through production and exchange of diffusible signal molecules in a cell-cell regulatory system termed quorum sensing. In the generic quorum sensor a positive feedback in the production of signal molecules defines the conditions at which the collective behavior switches on. In spite of its conceptual simplicity, a proper measure of biofilm colony "size" appears to be lacking. We establish that the cell density multiplied by a geometric factor which incorporates the boundary conditions constitutes an appropriate size measure. The geometric factor is the square of the radius for a spherical colony or a hemisphere attached to a reflecting surface. If surrounded by a rapidly exchanged medium, the geometric factor is divided by three. For a disk-shaped biofilm the geometric factor is the horizontal dimension multiplied by the height, and the square of the height of the biofilm if there is significant flow above the biofilm. A remarkably simple factorized expression for the size is obtained, which separates the all-or-none ignition caused by the positive feedback from the smoother activation outside the switching region.

Formulation of probabilistic models of protein structure in atomic detail using the reference ratio method.

We propose a method to formulate probabilistic models of protein structure in atomic detail, for a given amino acid sequence, based on Bayesian principles, while retaining a close link to physics. We start from two previously developed probabilistic models of protein structure on a local length scale, which concern the dihedral angles in main chain and side chains, respectively. Conceptually, this constitutes a probabilistic and continuous alternative to the use of discrete fragment and rotamer libraries. The local model is combined with a nonlocal model that involves a small number of energy terms according to a physical force field, and some information on the overall secondary structure content. In this initial study we focus on the formulation of the joint model and the evaluation of the use of an energy vector as a descriptor of a protein's nonlocal structure; hence, we derive the parameters of the nonlocal model from the native structure without loss of generality. The local and nonlocal models are combined using the reference ratio method, which is a well-justified probabilistic construction. For evaluation, we use the resulting joint models to predict the structure of four proteins. The results indicate that the proposed method and the probabilistic models show considerable promise for probabilistic protein structure prediction and related applications.

Kinetic model for signal binding to the Quorum sensing regulator LasR.

We propose a kinetic model for the activation of the las regulon in the opportunistic pathogen Pseudomonas aeruginosa. The model is based on in vitro data and accounts for the LasR dimerization and consecutive activation by binding of two OdDHL signal molecules. Experimentally, the production of the active LasR quorum-sensing regulator was studied in an Escherichia coli background as a function of signal molecule concentration. The functional activity of the regulator was monitored via a GFP reporter fusion to lasB expressed from the native lasB promoter. The new data shows that the active form of the LasR dimer binds two signal molecules cooperatively and that the timescale for reaching saturation is independent of the signal molecule concentration. This favors a picture where the dimerized regulator is protected against proteases and remains protected as it is activated through binding of two successive signal molecules. In absence of signal molecules, the dimerized regulator can dissociate and degrade through proteolytic turnover of the monomer. This resolves the apparent contradiction between our data and recent reports that the fully protected dimer is able to "degrade" when the induction of LasR ceases.

PHAISTOS: a framework for Markov chain Monte Carlo simulation and inference of protein structure.

We present a new software framework for Markov chain Monte Carlo sampling for simulation, prediction, and inference of protein structure. The software package contains implementations of recent advances in Monte Carlo methodology, such as efficient local updates and sampling from probabilistic models of local protein structure. These models form a probabilistic alternative to the widely used fragment and rotamer libraries. Combined with an easily extendible software architecture, this makes PHAISTOS well suited for Bayesian inference of protein structure from sequence and/or experimental data. Currently, two force-fields are available within the framework: PROFASI and OPLS-AA/L, the latter including the generalized Born surface area solvent model. A flexible command-line and configuration-file interface allows users quickly to set up simulations with the desired configuration. PHAISTOS is released under the GNU General Public License v3.0. Source code and documentation are freely available from http://phaistos.sourceforge.net. The software is implemented in C++ and has been tested on Linux and OSX platforms.

An efficient null model for conformational fluctuations in proteins.

Protein dynamics play a crucial role in function, catalytic activity, and pathogenesis. Consequently, there is great interest in computational methods that probe the conformational fluctuations of a protein. However, molecular dynamics simulations are computationally costly and therefore are often limited to comparatively short timescales. TYPHON is a probabilistic method to explore the conformational space of proteins under the guidance of a sophisticated probabilistic model of local structure and a given set of restraints that represent nonlocal interactions, such as hydrogen bonds or disulfide bridges. The choice of the restraints themselves is heuristic, but the resulting probabilistic model is well-defined and rigorous. Conceptually, TYPHON constitutes a null model of conformational fluctuations under a given set of restraints. We demonstrate that TYPHON can provide information on conformational fluctuations that is in correspondence with experimental measurements. TYPHON provides a flexible, yet computationally efficient, method to explore possible conformational fluctuations in proteins.

Subtle Monte Carlo Updates in Dense Molecular Systems.

Although Markov chain Monte Carlo (MC) simulation is a potentially powerful approach for exploring conformational space, it has been unable to compete with molecular dynamics (MD) in the analysis of high density structural states, such as the native state of globular proteins. Here, we introduce a kinetic algorithm, CRISP, that greatly enhances the sampling efficiency in all-atom MC simulations of dense systems. The algorithm is based on an exact analytical solution to the classic chain-closure problem, making it possible to express the interdependencies among degrees of freedom in the molecule as correlations in a multivariate Gaussian distribution. We demonstrate that our method reproduces structural variation in proteins with greater efficiency than current state-of-the-art Monte Carlo methods and has real-time simulation performance on par with molecular dynamics simulations. The presented results suggest our method as a valuable tool in the study of molecules in atomic detail, offering a potential alternative to molecular dynamics for probing long time-scale conformational transitions.

Generative probabilistic models extend the scope of inferential structure determination.

Conventional methods for protein structure determination from NMR data rely on the ad hoc combination of physical forcefields and experimental data, along with heuristic determination of free parameters such as weight of experimental data relative to a physical forcefield. Recently, a theoretically rigorous approach was developed which treats structure determination as a problem of Bayesian inference. In this case, the forcefields are brought in as a prior distribution in the form of a Boltzmann factor. Due to high computational cost, the approach has been only sparsely applied in practice. Here, we demonstrate that the use of generative probabilistic models instead of physical forcefields in the Bayesian formalism is not only conceptually attractive, but also improves precision and efficiency. Our results open new vistas for the use of sophisticated probabilistic models of biomolecular structure in structure determination from experimental data.

Neural mechanism of activity spread in the cat motor cortex and its relation to the intrinsic connectivity.

Motor cortical points are linked by intrinsic horizontal connections having a recurrent network topology. However, it is not known whether neural activity can propagate over the area covered by these intrinsic connections and whether there are spatial anisotropies of synaptic strength, as opposed to synaptic density. Moreover, the mechanisms by which activity spreads have yet to be determined. To address these issues, an 8 × 8 microelectrode array was inserted in the forelimb area of the cat motor cortex (MCx). The centre of the array had a laser etched hole ∼500 µm in diameter. A microiontophoretic pipette, with a tip diameter of 2-3 µm, containing bicuculline methiodide (BIC) was inserted in the hole and driven to a depth of 1200-1400 µm from the cortical surface. BIC was ejected for ∼2min from the tip of the micropipette with positive direct current ranging between 20 and 40 nA in different experiments. This produced spontaneous nearly periodic bursts (0.2-1.0 Hz) of multi-unit activity in a radius of about 400 µm from the tip of the micropipette. The bursts of neural activity spread at a velocity of 0.11-0.24 ms⁻¹ (mean=0.14 mm ms⁻¹, SD=0.05)with decreasing amplitude.The area activated was on average 7.22 mm² (SD=0.91 mm²), or ∼92% of the area covered by the recording array. The mode of propagation was determined to occur by progressive recruitment of cortical territory, driven by a central locus of activity of some 400 µm in radius. Thus, activity did not propagate as a wave. Transection of the connections between the thalamus and MCx did not significantly alter the propagation velocity or the size of the recruited area, demonstrating that the bursts spread along the routes of intrinsic cortical connectivity. These experiments demonstrate that neural activity initiated within a small motor cortical locus (≤ 400 µm in radius) can recruit a relatively large neighbourhood in which a variety of muscles acting at several forelimb joints are represented. These results support the hypothesis that the MCx controls the forelimb musculature in an integrated and anticipatory manner based on a recurrent network topology