A practical guide to the simultaneous determination of protein structure and dynamics using metainference

Accurate protein structural ensembles can be determined with metainference, a Bayesian inference method that integrates experimental information with prior knowledge of the system and deals with all sources of uncertainty and errors as well as with system heterogeneity. Furthermore, metainference can be implemented using the metadynamics approach, which enables the computational study of complex biological systems requiring extensive conformational sampling. In this chapter, we provide a step-by-step guide to perform and analyse metadynamic metainference simulations using the ISDB module of the open-source PLUMED library, as well as a series of practical tips to avoid common mistakes. Specifically, we will guide the reader in the process of learning how to model the structural ensemble of a small disordered peptide by combining state-of-the-art molecular mechanics force fields with nuclear magnetic resonance data, including chemical shifts, scalar couplings and residual dipolar couplings.Comment: 49 pages, 9 figure

Identifying Alternative Hyper-Splicing Signatures in MG-Thymoma by Exon Arrays

BACKGROUND: The vast majority of human genes (>70%) are alternatively spliced. Although alternative pre-mRNA processing is modified in multiple tumors, alternative hyper-splicing signatures specific to particular tumor types are still lacking. Here, we report the use of Affymetrix Human Exon Arrays to spot hyper-splicing events characteristic of myasthenia gravis (MG)-thymoma, thymic tumors which develop in patients with MG and discriminate them from colon cancer changes. METHODOLOGY/PRINCIPAL FINDINGS: We combined GO term to parent threshold-based and threshold-independent ad-hoc functional statistics with in-depth analysis of key modified transcripts to highlight various exon-specific changes. These denote alternative splicing in MG-thymoma tumors compared to healthy human thymus and to in-house and Affymetrix datasets from colon cancer and healthy tissues. By using both global and specific, term-to-parent Gene Ontology (GO) statistical comparisons, our functional integrative ad-hoc method allowed the detection of disease-relevant splicing events. CONCLUSIONS/SIGNIFICANCE: Hyper-spliced transcripts spanned several categories, including the tumorogenic ERBB4 tyrosine kinase receptor and the connective tissue growth factor CTGF, as well as the immune function-related histocompatibility gene HLA-DRB1 and interleukin (IL)19, two muscle-specific collagens and one myosin heavy chain gene; intriguingly, a putative new exon was discovered in the MG-involved acetylcholinesterase ACHE gene. Corresponding changes in spliceosome composition were indicated by co-decreases in the splicing factors ASF/SF(2) and SC35. Parallel tumor-associated changes occurred in colon cancer as well, but the majority of the apparent hyper-splicing events were particular to MG-thymoma and could be validated by Fluorescent In-Situ Hybridization (FISH), Reverse Transcription-Polymerase Chain Reaction (RT-PCR) and mass spectrometry (MS) followed by peptide sequencing. Our findings demonstrate a particular alternative hyper-splicing signature for transcripts over-expressed in MG-thymoma, supporting the hypothesis that alternative hyper-splicing contributes to shaping the biological functions of these and other specialized tumors and opening new venues for the development of diagnosis and treatment approaches

Dimer Formation Enhances Structural Differences between Amyloid β-Protein (1–40) and (1–42): An Explicit-Solvent Molecular Dynamics Study

Amyloid -protein (A) is central to the pathology of Alzheimer's disease. A 5% difference in the primary structure of the two predominant alloforms, A and A, results in distinct assembly pathways and toxicity properties. Discrete molecular dynamics (DMD) studies of A and A assembly resulted in alloform-specific oligomer size distributions consistent with experimental findings. Here, a large ensemble of DMD–derived A and A monomers and dimers was subjected to fully atomistic molecular dynamics (MD) simulations using the OPLS-AA force field combined with two water models, SPCE and TIP3P. The resulting all-atom conformations were slightly larger, less compact, had similar turn and lower -strand propensities than those predicted by DMD. Fully atomistic A and A monomers populated qualitatively similar free energy landscapes. In contrast, the free energy landscape of A dimers indicated a larger conformational variability in comparison to that of A dimers. A dimers were characterized by an increased flexibility in the N-terminal region D1-R5 and a larger solvent exposure of charged amino acids relative to A dimers. Of the three positively charged amino acids, R5 was the most and K16 the least involved in salt bridge formation. This result was independent of the water model, alloform, and assembly state. Overall, salt bridge propensities increased upon dimer formation. An exception was the salt bridge propensity of K28, which decreased upon formation of A dimers and was significantly lower than in A dimers. The potential relevance of the three positively charged amino acids in mediating the A oligomer toxicity is discussed in the light of available experimental data

Molecular dynamics simulation of tau peptides for the investigation of conformational changes induced by specific phosphorylation patterns

The Tau protein plays an important role due to its biomolecular interactions in neurodegenerative diseases. The lack of stable structure and various posttranslational modifications such as phosphorylation at various sites in the Tau protein pose a challenge for many experimental methods that are traditionally used to study protein folding and aggregation. Atomistic molecular dynamics (MD) simulations can help around deciphering relationship between phosphorylation and various intermediate and stable conformations of the Tau protein which occur on longer timescales. This chapter outlines protocols for the preparation, execution, and analysis of all-atom MD simulations of a 21-amino acid-long phosphorylated Tau peptide with the aim of generating biologically relevant structural and dynamic information. The simulations are done in explicit solvent and starting from nearly extended configurations of the peptide. The scaled MD method implemented in AMBER14 was chosen to achieve enhanced conformational sampling in addition to a conventional MD approach, thereby allowing the characterization of folding for such an intrinsically disordered peptide at 293 K. Emphasis is placed on the analysis of the simulation trajectories to establish correlations with NMR data (i.e., chemical shifts and NOEs). Finally, in-depth discussions are provided for commonly encountered problems

A Practical Guide to the Simultaneous Determination of Protein Structure and Dynamics Using Metainference

Accurate protein structural ensembles can be determined with metainference, a Bayesian inference method that integrates experimental information with prior knowledge of the system and deals with all sources of uncertainty and errors as well as with system heterogeneity. Furthermore, metainference can be implemented using the metadynamics approach, which enables the computational study of complex biological systems requiring extensive conformational sampling. In this chapter, we provide a step-by-step guide to perform and analyse metadynamic metainference simulations using the ISDB module of the open-source PLUMED library, as well as a series of practical tips to avoid common mistakes. Specifically, we will guide the reader in the process of learning how to model the structural ensemble of a small disordered peptide by combining state-of-the-art molecular mechanics force fields with nuclear magnetic resonance data, including chemical shifts, scalar couplings and residual dipolar couplings

Global ubiquitination analysis reveals extensive modification and proteasomal degradation of cowpox virus proteins, but preservation of viral cores

The emergence of Variola virus-like viruses by natural evolution of zoonotic Orthopoxviruses, like Cowpox virus (CPXV), is a global health threat. The proteasome is essential for poxvirus replication, making the viral components interacting with the ubiquitin-proteasome system attractive antiviral targets. We show that proteasome inhibition impairs CPXV replication by prevention of uncoating, suggesting that uncoating is mediated by proteasomal degradation of viral core proteins. Although Orthopoxvirus particles contain considerable amounts of ubiquitin, distinct modification sites are largely unknown. Therefore, for the first time, we analyzed globally ubiquitination sites in CPXV mature virion proteins using LC-MS/MS. Identification of 137 conserved sites in 54 viral proteins among five CPXV strains revealed extensive ubiquitination of structural core proteins. Moreover, since virions contained primarily K48-linked polyubiquitin, we hypothesized that core proteins are modified accordingly. However, quantitative analysis of ubiquitinated CPXV proteins early in infection showed no proteasomal degradation of core proteins. Instead, our data indicate that the recently suggested proteasomal regulation of the uncoating factor E5 is a prerequisite for uncoating. Expanding our understanding of poxvirus uncoating and elucidating a multitude of novel ubiquitination sites in poxvirus proteins, the present study verifies the major biological significance of ubiquitin in poxvirus infection