The Role of Non-native Interactions in the Folding of Knotted Proteins

Stochastic simulations of coarse-grained protein models are used to investigate the propensity to form knots in early stages of protein folding. The study is carried out comparatively for two homologous carbamoyltransferases, a natively-knotted N-acetylornithine carbamoyltransferase (AOTCase) and an unknotted ornithine carbamoyltransferase (OTCase). In addition, two different sets of pairwise amino acid interactions are considered: one promoting exclusively native interactions, and the other additionally including non-native quasi-chemical and electrostatic interactions. With the former model neither protein show a propensity to form knots. With the additional non-native interactions, knotting propensity remains negligible for the natively-unknotted OTCase while for AOTCase it is much enhanced. Analysis of the trajectories suggests that the different entanglement of the two transcarbamylases follows from the tendency of the C-terminal to point away from (for OTCase) or approach and eventually thread (for AOTCase) other regions of partly-folded protein. The analysis of the OTCase/AOTCase pair clarifies that natively-knotted proteins can spontaneously knot during early folding stages and that non-native sequence-dependent interactions are important for promoting and disfavoring early knotting events.Comment: Accepted for publication on PLOS Computational Biolog

Folding Pathways of a Knotted Protein with a Realistic Atomistic Force Field

We report on atomistic simulation of the folding of a natively-knotted protein, MJ0366, based on a realistic force field. To the best of our knowledge this is the first reported effort where a realistic force field is used to investigate the folding pathways of a protein with complex native topology. By using the dominant-reaction pathway scheme we collected about 30 successful folding trajectories for the 82-amino acid long trefoil-knotted protein. Despite the dissimilarity of their initial unfolded configuration, these trajectories reach the natively-knotted state through a remarkably similar succession of steps. In particular it is found that knotting occurs essentially through a threading mechanism, involving the passage of the C-terminal through an open region created by the formation of the native $\beta$-sheet at an earlier stage. The dominance of the knotting by threading mechanism is not observed in MJ0366 folding simulations using simplified, native-centric models. This points to a previously underappreciated role of concerted amino acid interactions, including non-native ones, in aiding the appropriate order of contact formation to achieve knotting

Can Polarity-Inverted Surfactants Self-Assemble in Nonpolar Solvents

We investigate the self-assembly process of a surfactant with inverted polarity in water and cyclohexane using both all-atom and coarse grained hybrid particle-field molecular dynamics simulations. Unlike conventional surfactants, the molecule under study, proposed in a recent experiment, is formed by a rigid and compact hydrophobic adamantane moiety, and a long and floppy triethylene glycol tail. In water, we report the formation of stable inverted micelles with the adamantane heads grouping together into a hydrophobic core, and the tails forming hydrogen bonds with water. By contrast, microsecond simulations do not provide evidence of stable micelle formation in cyclohexane. Validating the computational results by comparison with experimental diffusion constant and small-angle X-ray scattering intensity, we show that at laboratory thermodynamic conditions the mixture resides in the supercritical region of the phase diagram, where aggregated and free surfactant states co-exist in solution. Our simulations also provide indications about how to escape this region, to produce thermodynamically stable micellar aggregates.Comment: 14 pages, 10 Figures, accepted for publication (2020

Dominant Folding Pathways of a WW Domain

We investigate the folding mechanism of the WW domain Fip35 using a realistic atomistic force field by applying the Dominant Reaction Pathways (DRP) approach. We find evidence for the existence of two folding pathways, which differ by the order of formation of the two hairpins. This result is consistent with the analysis of the experimental data on the folding kinetics of WW domains and with the results obtained from large-scale molecular dynamics (MD) simulations of this system. Free-energy calculations performed in two coarse-grained models support the robustness of our results and suggest that the qualitative structure of the dominant paths are mostly shaped by the native interactions. Computing a folding trajectory in atomistic detail only required about one hour on 48 CPU's. The gain in computational efficiency opens the door to a systematic investigation of the folding pathways of a large number of globular proteins

Multidimensional severity assessment in bronchiectasis:An analysis of 7 European cohorts.

INTRODUCTION: Bronchiectasis is a multidimensional disease associated with substantial morbidity and mortality. Two disease-specific clinical prediction tools have been developed, the Bronchiectasis Severity Index (BSI) and the FACED score, both of which stratify patients into severity risk categories to predict the probability of mortality. METHODS: We aimed to compare the predictive utility of BSI and FACED in assessing clinically relevant disease outcomes across seven European cohorts independent of their original validation studies. RESULTS: The combined cohorts totalled 1612. Pooled analysis showed that both scores had a good discriminatory predictive value for mortality (pooled area under the curve (AUC) 0.76, 95% CI 0.74 to 0.78 for both scores) with the BSI demonstrating a higher sensitivity (65% vs 28%) but lower specificity (70% vs 93%) compared with the FACED score. Calibration analysis suggested that the BSI performed consistently well across all cohorts, while FACED consistently overestimated mortality in 'severe' patients (pooled OR 0.33 (0.23 to 0.48), p<0.0001). The BSI accurately predicted hospitalisations (pooled AUC 0.82, 95% CI 0.78 to 0.84), exacerbations, quality of life (QoL) and respiratory symptoms across all risk categories. FACED had poor discrimination for hospital admissions (pooled AUC 0.65, 95% CI 0.63 to 0.67) with low sensitivity at 16% and did not consistently predict future risk of exacerbations, QoL or respiratory symptoms. No association was observed with FACED and 6 min walk distance (6MWD) or lung function decline. CONCLUSION: The BSI accurately predicts mortality, hospital admissions, exacerbations, QoL, respiratory symptoms, 6MWD and lung function decline in bronchiectasis, providing a clinically relevant evaluation of disease severity

Application of Consensus Scoring and Principal Component Analysis for Virtual Screening against β-Secretase (BACE-1)

BACKGROUND: In order to identify novel chemical classes of β-secretase (BACE-1) inhibitors, an alternative scoring protocol, Principal Component Analysis (PCA), was proposed to summarize most of the information from the original scoring functions and re-rank the results from the virtual screening against BACE-1. METHOD: Given a training set (50 BACE-1 inhibitors and 9950 inactive diverse compounds), three rank-based virtual screening methods, individual scoring, conventional consensus scoring and PCA, were judged by the hit number in the top 1% of the ranked list. The docking poses were generated by Surflex, five scoring functions (Surflex_Score, D_Score, G_Score, ChemScore, and PMF_Score) were used for pose extraction. For each pose group, twelve scoring functions (Surflex_Score, D_Score, G_Score, ChemScore, PMF_Score, LigScore1, LigScore2, PLP1, PLP2, jain, Ludi_1, and Ludi_2) were used for the pose rank. For a test set, 113,228 chemical compounds (Sigma-Aldrich® corporate chemical directory) were docked by Surflex, then ranked by the same three ranking methods motioned above to select the potential active compounds for experimental test. RESULTS: For the training set, the PCA approach yielded consistently superior rankings compared to conventional consensus scoring and single scoring. For the test set, the top 20 compounds according to conventional consensus scoring were experimentally tested, no inhibitor was found. Then, we relied on PCA scoring protocol to test another different top 20 compounds and two low micromolar inhibitors (S450588 and 276065) were emerged through the BACE-1 fluorescence resonance energy transfer (FRET) assay. CONCLUSION: The PCA method extends the conventional consensus scoring in a quantitative statistical manner and would appear to have considerable potential for chemical screening applications

Syndecan-4 Is Essential for Development of Concentric Myocardial Hypertrophy via Stretch-Induced Activation of the Calcineurin-NFAT Pathway

Sustained pressure overload leads to compensatory myocardial hypertrophy and subsequent heart failure, a leading cause of morbidity and mortality. Further unraveling of the cellular processes involved is essential for development of new treatment strategies. We have investigated the hypothesis that the transmembrane Z-disc proteoglycan syndecan-4, a co-receptor for integrins, connecting extracellular matrix proteins to the cytoskeleton, is an important signal transducer in cardiomyocytes during development of concentric myocardial hypertrophy following pressure overload. Echocardiographic, histochemical and cardiomyocyte size measurements showed that syndecan-4−/− mice did not develop concentric myocardial hypertrophy as found in wild-type mice, but rather left ventricular dilatation and dysfunction following pressure overload. Protein and gene expression analyses revealed diminished activation of the central, pro-hypertrophic calcineurin-nuclear factor of activated T-cell (NFAT) signaling pathway. Cardiomyocytes from syndecan-4−/−-NFAT-luciferase reporter mice subjected to cyclic mechanical stretch, a hypertrophic stimulus, showed minimal activation of NFAT (1.6-fold) compared to 5.8-fold increase in NFAT-luciferase control cardiomyocytes. Accordingly, overexpression of syndecan-4 or introducing a cell-permeable membrane-targeted syndecan-4 polypeptide (gain of function) activated NFATc4 in vitro. Pull-down experiments demonstrated a direct intracellular syndecan-4-calcineurin interaction. This interaction and activation of NFAT were increased by dephosphorylation of serine 179 (pS179) in syndecan-4. During pressure overload, phosphorylation of syndecan-4 was decreased, and association between syndecan-4, calcineurin and its co-activator calmodulin increased. Moreover, calcineurin dephosphorylated pS179, indicating that calcineurin regulates its own binding and activation. Finally, patients with hypertrophic myocardium due to aortic stenosis had increased syndecan-4 levels with decreased pS179 which was associated with increased NFAT activation. In conclusion, our data show that syndecan-4 is essential for compensatory hypertrophy in the pressure overloaded heart. Specifically, syndecan-4 regulates stretch-induced activation of the calcineurin-NFAT pathway in cardiomyocytes. Thus, our data suggest that manipulation of syndecan-4 may provide an option for therapeutic modulation of calcineurin-NFAT signaling