TADOSS: computational estimation of tandem domain swap stability.
Ontology highlight
ABSTRACT:
SUBMITTER: Lafita A
PROVIDER: S-EPMC6612889 | biostudies-other | 2018 Nov
REPOSITORIES: biostudies-other
Publications
TADOSS: computational estimation of tandem domain swap stability.
Lafita Aleix A Tian Pengfei P Best Robert B RB Bateman Alex A
Bioinformatics (Oxford, England) 20190701 14
<h4>Summary</h4>Proteins with highly similar tandem domains have shown an increased propensity for misfolding and aggregation. Several molecular explanations have been put forward, such as swapping of adjacent domains, but there is a lack of computational tools to systematically analyze them. We present the TAndem DOmain Swap Stability predictor (TADOSS), a method to computationally estimate the stability of tandem domain-swapped conformations from the structures of single domains, based on prev ...[more]
Similar Datasets
Project description:Three-dimensional (3D) domain swapping creates a bond between two or more protein molecules as they exchange their identical domains. Since the term '3D domain swapping' was first used to describe the dimeric structure of diphtheria toxin, the database of domain-swapped proteins has greatly expanded. Analyses of the now about 40 structurally characterized cases of domain-swapped proteins reveal that most swapped domains are at either the N or C terminus and that the swapped domains are diverse in their primary and secondary structures. In addition to tabulating domain-swapped proteins, we describe in detail several examples of 3D domain swapping which show the swapping of more than one domain in a protein, the structural evidence for 3D domain swapping in amyloid proteins, and the flexibility of hinge loops. We also discuss the physiological relevance of 3D domain swapping and a possible mechanism for 3D domain swapping. The present state of knowledge leads us to suggest that 3D domain swapping can occur under appropriate conditions in any protein with an unconstrained terminus. As domains continue to swap, this review attempts not only a summary of the known domain-swapped proteins, but also a framework for understanding future findings of 3D domain swapping.
| S-EPMC2373619 | biostudies-literature
Project description:AGel amyloidosis is a genetic degenerative disease characterized by the deposition of insoluble gelsolin protein aggregates in different tissues. Until recently, this disease was associated with two mutations of a single residue (Asp187 to Asn/Tyr) in the second domain of the protein. The general opinion is that pathogenic variants are not per se amyloidogenic but rather that the mutations trigger an aberrant proteolytic cascade, which results in the production of aggregation prone fragments. Here, we report the crystal structure of the second domain of gelsolin carrying the recently identified Gly167Arg mutation. This mutant dimerizes through a three-dimensional domain swapping mechanism, forming a tight but flexible assembly, which retains the structural topology of the monomer. To date, such dramatic conformational changes of this type have not been observed. Structural and biophysical characterizations reveal that the Gly167Arg mutation alone is responsible for the monomer to dimer transition and that, even in the context of the full-length protein, the pathogenic variant is prone to form dimers. These data suggest that, in addition to the well-known proteolytic-dependent mechanism, an alternative oligomerization pathway may participate in gelsolin misfolding and aggregation. We propose to integrate this alternative pathway into the current model of the disease that may also be relevant for other types of AGel amyloidosis, and other related diseases with similar underlying pathological mechanisms.
| S-EPMC5886171 | biostudies-literature
Project description:Confinement of the X chromosome into a territory for dosage compensation (DC) is a prime example of subnuclear compartmentalization for transcription regulation at the megabase scale. In D. melanogaster, two sex-specific non-coding (nc) RNAs roX1 and roX2 are transcribed from the X chromosome. They associate with the Male-specific lethal (MSL) complex, which by acetylating histone H4 lysine 16 (H4K16ac) confers approximately 2-fold upregulated expression of male X-linked genes. Current models explain X-over-autosome specificity based on the MSL2 subunit recognizing cis-regulatory DNA high-affinity sites (HAS). However, HAS motifs are also found on autosomes, indicating that additional factors have evolved to confer stable association of the MSL complex with the X. Here, we show that the low-complexity C-terminal domain (CTD) of MSL2 renders its recruitment to the X chromosome sensitive to roX ncRNAs. roX-MSL2-CTD form a stably condensated state, and functional analyses in Drosophila and mammalian cells reveal the critical importance of their interplay for DC in vivo. Replacing the CTD of mammalian MSL2 with that from Drosophila and expressing roX in cis is sufficient to nucleate ectopic DC in mammalian cells. Thus, the condensating nature of roX-MSL2CTD is the primary determinant for specific compartmentalization of the X in Drosophila.
2020-08-07 | GSE129831 | GEO
Project description:Intracellular signal transduction proteins typically utilize multiple interaction domains for proper targeting, and thus a broad diversity of distinct signaling complexes may be assembled. Considering the coordination of only two such domains, as in tandem Src homology 2 (SH2) domain constructs, gives rise to a kinetic scheme that is not adequately described by simple models used routinely to interpret in vitro binding measurements. To analyze the interactions between tandem SH2 domains and bisphosphorylated peptides, we formulated detailed kinetic models and applied them to the phosphoinositide 3-kinase p85 regulatory subunit/platelet-derived growth factor beta-receptor system. Data for this system from different in vitro assay platforms, including surface plasmon resonance, competition binding, and isothermal titration calorimetry, were reconciled to estimate the magnitude of the cooperativity characterizing the sequential binding of the high and low affinity SH2 domains (C-SH2 and N-SH2, respectively). Compared with values based on an effective volume approximation, the estimated cooperativity is 3 orders of magnitude lower, indicative of significant structural constraints. Homodimerization of full-length p85 was found to be an alternative mechanism for high avidity binding to phosphorylated platelet-derived growth factor receptors, which would render the N-SH2 domain dispensable for receptor binding.
| S-EPMC2276335 | biostudies-literature
Project description:TIS11d is a member of the CCCH-type family of tandem zinc finger (TZF) proteins; the TZF domain of TIS11d (residues 151-220) is sufficient to bind and destabilize its target mRNAs with high specificity. In this study, the TZF domain of TIS11d is simulated in an aqueous environment in both the free and RNA-bound states. Multiple nanosecond timescale molecular dynamics trajectories of TIS11d wild-type and E157R/E195K mutant with different RNA sequences were performed to investigate the molecular basis for RNA binding specificities of this TZF domain. A variety of measures of the protein structure, fluctuations, and dynamics were used to analyze the trajectories. The results of this study support the following conclusions: (1) the structure of the two fingers is maintained in the free state but a global reorientation occurs to yield a more compact structure; (2) mutation of the glutamate residues at positions 157 and 195 to arginine and lysine, respectively, affects the RNA recognition by this TIS11d mutant in agreement with the findings of Pagano et al. (J Biol Chem 2007; 282:8883-8894); and (3) we predict that the E157R/E195K mutant will present a more relaxed RNA binding specificity relative to wild-type TIS11d based on the more favorable nonsequence-specific Coulomb interaction of the two positively charged residues at positions 157 and 195 with the RNA backbone, which compensates for a partial loss of the stacking interaction of aromatic side chains with the RNA bases.
| S-EPMC2895246 | biostudies-literature
Project description:Methods for efficient and accurate prediction of RNA structure are increasingly valuable, given the current rapid advances in understanding the diverse functions of RNA molecules in the cell. To enhance the accuracy of secondary structure predictions, we developed and refined optimization techniques for the estimation of energy parameters. We build on two previous approaches to RNA free-energy parameter estimation: (1) the Constraint Generation (CG) method, which iteratively generates constraints that enforce known structures to have energies lower than other structures for the same molecule; and (2) the Boltzmann Likelihood (BL) method, which infers a set of RNA free-energy parameters that maximize the conditional likelihood of a set of reference RNA structures. Here, we extend these approaches in two main ways: We propose (1) a max-margin extension of CG, and (2) a novel linear Gaussian Bayesian network that models feature relationships, which effectively makes use of sparse data by sharing statistical strength between parameters. We obtain significant improvements in the accuracy of RNA minimum free-energy pseudoknot-free secondary structure prediction when measured on a comprehensive set of 2518 RNA molecules with reference structures. Our parameters can be used in conjunction with software that predicts RNA secondary structures, RNA hybridization, or ensembles of structures. Our data, software, results, and parameter sets in various formats are freely available at http://www.cs.ubc.ca/labs/beta/Projects/RNA-Params.
| S-EPMC2995392 | biostudies-literature
Project description:Computational design of surface charge-charge interactions has been demonstrated to be an effective way to increase both the thermostability and the stability of proteins. To test the robustness of this approach for proteins with predominantly beta-sheet secondary structure, the chicken isoform of the Fyn SH3 domain was used as a model system. Computational analysis of the optimal distribution of surface charges showed that the increase in favorable energy per substitution begins to level off at five substitutions; hence, the designed Fyn sequence contained four charge reversals at existing charged positions and one introduction of a new charge. Three additional variants were also constructed to explore stepwise contributions of these substitutions to Fyn stability. The thermodynamic stabilities of the variants were experimentally characterized using differential scanning calorimetry and far-UV circular dichroism spectroscopy and are in very good agreement with theoretical predictions from the model. The designed sequence was found to have increased the melting temperature, DeltaT (m) = 12.3 +/- 0.2 degrees C, and stability, DeltaDeltaG(25 degrees C) = 7.1 +/- 2.2 kJ/mol, relative to the wild-type protein. The experimental data suggest that a significant increase in stability can be achieved through a very small number of amino acid substitutions. Consistent with a number of recent studies, the presented results clearly argue for a seminal role of surface charge-charge interactions in determining protein stability and suggest that the optimization of surface interactions can be an attractive strategy to complement algorithms optimizing interactions in the protein core to further enhance protein stability.
| S-EPMC2222822 | biostudies-literature
Project description:A central challenge in computational modeling of biological systems is the determination of the model parameters. Typically, only a fraction of the parameters (such as kinetic rate constants) are experimentally measured, while the rest are often fitted. The fitting process is usually based on experimental time course measurements of observables, which are used to assign parameter values that minimize some measure of the error between these measurements and the corresponding model prediction. The measurements, which can come from immunoblotting assays, fluorescent markers, etc., tend to be very noisy and taken at a limited number of time points. In this work we present a new approach to the problem of parameter selection of biological models. We show how one can use a dynamic recursive estimator, known as extended Kalman filter, to arrive at estimates of the model parameters. The proposed method follows. First, we use a variation of the Kalman filter that is particularly well suited to biological applications to obtain a first guess for the unknown parameters. Secondly, we employ an a posteriori identifiability test to check the reliability of the estimates. Finally, we solve an optimization problem to refine the first guess in case it should not be accurate enough. The final estimates are guaranteed to be statistically consistent with the measurements. Furthermore, we show how the same tools can be used to discriminate among alternate models of the same biological process. We demonstrate these ideas by applying our methods to two examples, namely a model of the heat shock response in E. coli, and a model of a synthetic gene regulation system. The methods presented are quite general and may be applied to a wide class of biological systems where noisy measurements are used for parameter estimation or model selection.
| S-EPMC2832681 | biostudies-other
Project description:MotivationProtein domain duplications are a major contributor to the functional diversification of protein families. These duplications can occur one at a time through single domain duplications, or as tandem duplications where several consecutive domains are duplicated together as part of a single evolutionary event. Existing methods for inferring domain-level evolutionary events are based on reconciling domain trees with gene trees. While some formulations consider multiple domain duplications, they do not explicitly model tandem duplications; this leads to inaccurate inference of which domains duplicated together over the course of evolution.ResultsHere, we introduce a reconciliation-based framework that considers the relative positions of domains within extant sequences. We use this information to uncover tandem domain duplications within the evolutionary history of these genes. We devise an integer linear programming approach that solves our problem exactly, and a heuristic approach that works well in practice. We perform extensive simulation studies to demonstrate that our approaches can accurately uncover single and tandem domain duplications, and additionally test our approach on a well-studied orthogroup where lineage-specific domain expansions exhibit varying and complex domain duplication patterns.Availability and implementationCode is available on github at https://github.com/Singh-Lab/TandemDuplications.Supplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC8275333 | biostudies-literature
Project description:Cytochromes P450 (P450, CYP) metabolize a wide variety of endogenous and exogenous lipophilic molecules, including most drugs. Sterol 14α-demethylase (CYP51) is a target for antifungal drugs known as conazoles. Using X-ray crystallography, we have discovered a domain-swap homodimerization mode in CYP51 from a human pathogen, Acanthamoeba castellanii CYP51 (AcCYP51). Recombinant AcCYP51 with a truncated transmembrane helix was purified as a heterogeneous mixture corresponding to the dimer and monomer units. Spectral analyses of these two populations have shown that the CO-bound ferrous form of the dimeric protein absorbed at 448 nm (catalytically competent form), whereas the monomeric form absorbed at 420 nm (catalytically incompetent form). AcCYP51 dimerized head-to-head via N-termini swapping, resulting in formation of a nonplanar protein-protein interface exceeding 2000 Å2 with a total solvation energy gain of -35.4 kcal/mol. In the dimer, the protomers faced each other through the F and G α-helices, thus blocking the substrate access channel. In the presence of the drugs clotrimazole and isavuconazole, the AcCYP51 drug complexes crystallized as monomers. Although clotrimazole-bound AcCYP51 adopted a typical CYP monomer structure, isavuconazole-bound AcCYP51 failed to refold 74 N-terminal residues. The failure of AcCYP51 to fully refold upon inhibitor binding in vivo would cause an irreversible loss of a structurally aberrant enzyme through proteolytic degradation. This assumption explains the superior potency of isavuconazole against A. castellanii The dimerization mode observed in this work is compatible with membrane association and may be relevant to other members of the CYP family of biologic, medical, and pharmacological importance. SIGNIFICANCE STATEMENT: We investigated the mechanism of action of antifungal drugs in the human pathogen Acanthamoeba castellanii. We discovered that the enzyme target [Acanthamoeba castellanii sterol 14α-demethylase (AcCYP51)] formed a dimer via an N-termini swap, whereas drug-bound AcCYP51 was monomeric. In the AcCYP51-isavuconazole complex, the protein target failed to refold 74 N-terminal residues, suggesting a fundamentally different mechanism of AcCYP51 inactivation than only blocking the active site. Proteolytic degradation of a structurally aberrant enzyme would explain the superior potency of isavuconazole against A. castellanii.
| S-EPMC7674934 | biostudies-literature