Bayesian inference of protein ensembles from SAXS data

The inherent flexibility of intrinsically disordered proteins (IDPs) and multi-domain proteins with intrinsically disordered regions (IDRs) presents challenges to structural analysis. These macromolecules need to be represented by an ensemble of conformations, rather than a single structure. Small-angle X-ray scattering (SAXS) experiments capture ensemble-averaged data for the set of conformations. We present a Bayesian approach to ensemble inference from SAXS data, called Bayesian ensemble SAXS (BE-SAXS). We address two issues with existing methods: the use of a finite ensemble of structures to represent the underlying distribution, and the selection of that ensemble as a subset of an initial pool of structures. This is achieved through the formulation of a Bayesian posterior of the conformational space. BE-SAXS modifies a structural prior distribution in accordance with the experimental data. It uses multi-step expectation maximization, with alternating rounds of Markov-chain Monte Carlo simulation and empirical Bayes optimization. We demonstrate the method by employing it to obtain a conformational ensemble of the antitoxin PaaA2 and comparing the results to a published ensemble.ISSN:1463-9084ISSN:1463-907

ENCORE:Software for Quantitative Ensemble Comparison

There is increasing evidence that protein dynamics and conformational changes can play an important role in modulating biological function. As a result, experimental and computational methods are being developed, often synergistically, to study the dynamical heterogeneity of a protein or other macromolecules in solution. Thus, methods such as molecular dynamics simulations or ensemble refinement approaches have provided conformational ensembles that can be used to understand protein function and biophysics. These developments have in turn created a need for algorithms and software that can be used to compare structural ensembles in the same way as the root-mean-square-deviation is often used to compare static structures. Although a few such approaches have been proposed, these can be difficult to implement efficiently, hindering a broader applications and further developments. Here, we present an easily accessible software toolkit, called ENCORE, which can be used to compare conformational ensembles generated either from simulations alone or synergistically with experiments. ENCORE implements three previously described methods for ensemble comparison, that each can be used to quantify the similarity between conformational ensembles by estimating the overlap between the probability distributions that underlie them. We demonstrate the kinds of insights that can be obtained by providing examples of three typical use-cases: comparing ensembles generated with different molecular force fields, assessing convergence in molecular simulations, and calculating differences and similarities in structural ensembles refined with various sources of experimental data. We also demonstrate efficient computational scaling for typical analyses, and robustness against both the size and sampling of the ensembles. ENCORE is freely available and extendable, integrates with the established MDAnalysis software package, reads ensemble data in many common formats, and can work with large trajectory files

Inference of structure ensembles of flexible biomolecules from sparse, averaged data

We present the theoretical foundations of a general principle to infer structure ensembles of flexible biomolecules from spatially and temporally averaged data obtained in biophysical experiments. The central idea is to compute the Kullback-Leibler optimal modification of a given prior distribution [Image: see text] with respect to the experimental data and its uncertainty. This principle generalizes the successful inferential structure determination method and recently proposed maximum entropy methods. Tractability of the protocol is demonstrated through the analysis of simulated nuclear magnetic resonance spectroscopy data of a small peptide

Bioinformatics analysis identifies several intrinsically disordered human E3 ubiquitin-protein ligases

The ubiquitin-proteasome system targets misfolded proteins for degradation. Since the accumulation of such proteins is potentially harmful for the cell, their prompt removal is important. E3 ubiquitin-protein ligases mediate substrate ubiquitination by bringing together the substrate with an E2 ubiquitin-conjugating enzyme, which transfers ubiquitin to the substrate. For misfolded proteins, substrate recognition is generally delegated to molecular chaperones that subsequently interact with specific E3 ligases. An important exception is San1, a yeast E3 ligase. San1 harbors extensive regions of intrinsic disorder, which provide both conformational flexibility and sites for direct recognition of misfolded targets of vastly different conformations. So far, no mammalian ortholog of San1 is known, nor is it clear whether other E3 ligases utilize disordered regions for substrate recognition. Here, we conduct a bioinformatics analysis to examine >600 human and S. cerevisiae E3 ligases to identify enzymes that are similar to San1 in terms of function and/or mechanism of substrate recognition. An initial sequence-based database search was found to detect candidates primarily based on the homology of their ordered regions, and did not capture the unique disorder patterns that encode the functional mechanism of San1. However, by searching specifically for key features of the San1 sequence, such as long regions of intrinsic disorder embedded with short stretches predicted to be suitable for substrate interaction, we identified several E3 ligases with these characteristics. Our initial analysis revealed that another remarkable trait of San1 is shared with several candidate E3 ligases: long stretches of complete lysine suppression, which in San1 limits auto-ubiquitination. We encode these characteristic features into a San1 similarity-score, and present a set of proteins that are plausible candidates as San1 counterparts in humans. In conclusion, our work indicates that San1 is not a unique case, and that several other yeast and human E3 ligases have sequence properties that may allow them to recognize substrates by a similar mechanism as San1

Wind Turbine Adaptive Blade Integrated Design and Analysis

This project aims to develop efficient and robust tools for optimal design of wind turbine adaptive blades. In general, wind turbine adaptive blade design is an aero-structure coupled design process, in which, the evaluation of aerodynamic performance cannot be carried out precisely without structural deformation analysis of the adaptive blade. However, employing finite element analysis (FEA) based structural analysis commercial packages as part of the aerodynamic objective evaluation process has been proven time consuming and it results in inefficient and redundant design optimisation of adaptive blades caused by elastic-coupled (bend-twist or stretch-twist) iteration. In order to achieve the goal of wind turbine adaptive blade integrated design and analysis, this project is carried out from three aspects. Firstly, a general geometrically linear model for thin-walled composite beams with multi-cell, non-uniform cross-section and arbitrary lay-ups under various types of loadings is developed for implementing structural deformation analysis. After that, this model is validated by a simple box-beam, single- and multi-cell wind turbine blades. Through validation, it denotes that this thin-walled composite beam model is efficient and accurate for predicting the structural deformations compared to FEA based commercial packages (ANSYS). This developed beam model thus provides more probabilities for further investigations of dynamic performance of adaptive blades. Secondly in order to investigate the effects of aero elastic tailoring and implanting elastic coupling on aerodynamic performance of adaptive blades, auxiliary software tools with graphical interfaces are developed via MATLAB codes. Structural/material characteristics and configurations of adaptive blades (i.e. elastic coupling topology, layup configuration and material properties of blade) are defined by these auxiliary software tools. By interfacing these software tools to the structural analysers based on the developed thin-walled composite beam model to an aerodynamic performance evaluator, an integrated design environment is developed. Lastly, by using the developed thin-walled composite beam model as a search platform, the application of the decoupled design method, a method of design of smart aero-structures based on the concept of variable state design parameter, is also extended

Genome-wide interaction study of a proxy for stress-sensitivity and its prediction of major depressive disorder

Individual response to stress is correlated with neuroticism and is an important predictor of both neuroticism and the onset of major depressive disorder (MDD). Identification of the genetics underpinning individual differences in response to negative events (stress-sensitivity) may improve our understanding of the molecular pathways involved, and its association with stress-related illnesses. We sought to generate a proxy for stress-sensitivity through modelling the interaction between SNP allele and MDD status on neuroticism score in order to identify genetic variants that contribute to the higher neuroticism seen in individuals with a lifetime diagnosis of depression compared to unaffected individuals. Meta-analysis of genome-wide interaction studies (GWIS) in UK Biobank (N = 23,092) and Generation Scotland: Scottish Family Health Study (N = 7,155) identified no genome-wide significance SNP interactions. However, gene-based tests identified a genome-wide significant gene, ZNF366, a negative regulator of glucocorticoid receptor function implicated in alcohol dependence (p = 1.48x10-7; Bonferroni-corrected significance threshold p < 2.79x10-6). Using summary statistics from the stress-sensitivity term of the GWIS, SNP heritability for stress-sensitivity was estimated at 5.0%. In models fitting polygenic risk scores of both MDD and neuroticism derived from independent GWAS, we show that polygenic risk scores derived from the UK Biobank stress-sensitivity GWIS significantly improved the prediction of MDD in Generation Scotland. This study may improve interpretation of larger genome-wide association studies of MDD and other stress-related illnesses, and the understanding of the etiological mechanisms underpinning stress-sensitivity

Hair Cortisol in Twins : Heritability and Genetic Overlap with Psychological Variables and Stress-System Genes

A. Palotie on työryhmän jäsen.Hair cortisol concentration (HCC) is a promising measure of long-term hypothalamus-pituitary-adrenal (HPA) axis activity. Previous research has suggested an association between HCC and psychological variables, and initial studies of inter-individual variance in HCC have implicated genetic factors. However, whether HCC and psychological variables share genetic risk factors remains unclear. The aims of the present twin study were to: (i) assess the heritability of HCC; (ii) estimate the phenotypic and genetic correlation between HPA axis activity and the psychological variables perceived stress, depressive symptoms, and neuroticism; using formal genetic twin models and molecular genetic methods, i.e. polygenic risk scores (PRS). HCC was measured in 671 adolescents and young adults. These included 115 monozygotic and 183 dizygotic twin-pairs. For 432 subjects PRS scores for plasma cortisol, major depression, and neuroticism were calculated using data from large genome wide association studies. The twin model revealed a heritability for HCC of 72%. No significant phenotypic or genetic correlation was found between HCC and the three psychological variables of interest. PRS did not explain variance in HCC. The present data suggest that HCC is highly heritable. However, the data do not support a strong biological link between HCC and any of the investigated psychological variables.Peer reviewe

Evidence for Increased Genetic Risk Load for Major Depression in Patients Assigned to Electroconvulsive Therapy

Electroconvulsive therapy (ECT) is the treatment of choice for severe and treatment-resistant depression; disorder severity and unfavorable treatment outcomes are shown to be influenced by an increased genetic burden for major depression (MD). Here, we tested whether ECT assignment and response/nonresponse are associated with an increased genetic burden for major depression (MD) using polygenic risk score (PRS), which summarize the contribution of diseaserelated common risk variants. Fifty-one psychiatric inpatients suffering from a major depressive episode underwent ECT. MD-PRS were calculated for these inpatients and a separate population-based sample (n = 3,547 healthy; n = 426 self-reported depression) based on summary statistics from the Psychiatric Genomics Consortium MDD-working group (Cases: n = 59,851; Controls: n = 113,154). MD-PRS explained a significant proportion of disease status between ECT patients and healthy controls (p = .022, R2 = 1.173%); patients showed higher MD-PRS. MD-PRS in population-based depression self-reporters were intermediate between ECT patients and controls (n.s.). Significant associations between MD-PRS and ECT response (50% reduction in Hamilton depression rating scale scores) were not observed. Our findings indicate that ECT cohorts show an increased genetic burden for MD and are consistent with the hypothesis that treatment-resistant MD patients represent a subgroup with an increased genetic risk for MD. Larger samples are needed to better substantiate these findings

Association of whole-genome and NETRIN1 signaling pathway-derived polygenic risk scores for Major Depressive Disorder and white matter microstructure in UK Biobank

Background: Major depressive disorder is a clinically heterogeneous psychiatric disorder with a polygenic architecture. Genome-wide association studies have identified a number of risk-associated variants across the genome and have reported growing evidence of NETRIN1 pathway involvement. Stratifying disease risk by genetic variation within the NETRIN1 pathway may provide important routes for identification of disease mechanisms by focusing on a specific process, excluding heterogeneous risk-associated variation in other pathways. Here, we sought to investigate whether major depressive disorder polygenic risk scores derived from the NETRIN1 signaling pathway (NETRIN1-PRSs) and the whole genome, excluding NETRIN1 pathway genes (genomic-PRSs), were associated with white matter microstructure. Methods: We used two diffusion tensor imaging measures, fractional anisotropy (FA) and mean diffusivity (MD), in the most up-to-date UK Biobank neuroimaging data release (FA: n = 6401; MD: n = 6390). Results: We found significantly lower FA in the superior longitudinal fasciculus (β = −.035, p =.029) and significantly higher MD in a global measure of thalamic radiations (β =.029, p =.021), as well as higher MD in the superior (β =.034, p =.039) and inferior (β =.029, p =.043) longitudinal fasciculus and in the anterior (β =.025, p =.046) and superior (β =.027, p =.043) thalamic radiation associated with NETRIN1-PRS. Genomic-PRS was also associated with lower FA and higher MD in several tracts. Conclusions: Our findings indicate that variation in the NETRIN1 signaling pathway may confer risk for major depressive disorder through effects on a number of white matter tracts