Characterization of the impact of alternative splicing on protein dynamics: the cases of glutathione S-transferase and ectodysplasin-A isoforms

Cognate RNA-Binding Modes by the Alternative-Splicing Regulator MBNL1 Inferred from Molecular Dynamics

The muscleblind-like protein family (MBNL) plays a prominent role in the regulation of alternative splicing. Consequently, the loss of MBNL function resulting from sequestration by RNA hairpins triggers the development of a neuromuscular disease called myotonic dystrophy (DM). Despite the sequence and structural similarities between the four zinc-finger domains that form MBNL1, recent studies have revealed that the four binding domains have differentiated splicing activity. The dynamic behaviors of MBNL1 ZnFs were simulated using conventional molecular dynamics (cMD) and steered molecular dynamics (sMD) simulations of a structural model of MBNL1 protein to provide insights into the binding selectivity of the four zinc-finger (ZnF) domains toward the GpC steps in YGCY RNA sequence. In accordance with previous studies, our results suggest that both global and local residue fluctuations on each domain have great impacts on triggering alternative splicing, indicating that local motions in RNA-binding domains could modulate their affinity and specificity. In addition, all four ZnF domains provide a distinct RNA-binding environment in terms of structural sampling and mobility that may be involved in the differentiated MBNL1 splicing events reported in the literature.

Quantitative Characterization of Binding Pockets and Binding Complementarity by Means of Zernike Descriptors

In this work, we describe the application of the Zernike formalism to quantitatively characterize the binding pockets of two sets of biologically relevant systems. Such an approach, when applied to molecular dynamics trajectories, is able to pinpoint the subtle differences between very similar molecular regions and their impact on the local propensity to ligand binding, allowing us to quantify such differences. The statistical robustness of our procedure suggests that it is very suitable to describe protein binding sites and protein-ligand interactions within a rigorous and well-defined framework.

Molecular dynamics study of naturally existing cavity couplings in proteins

Couplings between protein sub-structures are a common property of protein dynamics. Some of these couplings are especially interesting since they relate to function and its regulation. In this article we have studied the case of cavity couplings because cavities can host functional sites, allosteric sites, and are the locus of interactions with the cell milieu. We have divided this problem into two parts. In the first part, we have explored the presence of cavity couplings in the natural dynamics of 75 proteins, using 20 ns molecular dynamics simulations. For each of these proteins, we have obtained two trajectories around their native state. After applying a stringent filtering procedure, we found significant cavity correlations in 60% of the proteins. We analyze and discuss the structure origins of these correlations, including neighbourhood, cavity distance, etc. In the second part of our study, we have used longer simulations (≥100 ns) from the MoDEL project, to obtain a broader view of cavity couplings, particularly about their dependence on time. Using moving window computations we explored the fluctuations of cavity couplings along time, finding that these couplings could fluctuate substantially during the trajectory, reaching in several cases correlations above 0.25/0.5. In summary, we describe the structural origin and the variations with time of cavity couplings. We complete our work with a brief discussion of the biological implications of these results.

Computational extraction of a neural molecular network through alternative splicing

BACKGROUND

Generally, the results of high throughput analyses contain information about gene expressions, and about exon expressions. Approximately 90% of primary protein-coding transcripts undergo alternative splicing in mammals. However, changes induced by alternative exons have not been properly analyzed for their impact on important molecular networks or their biological events. Even when alternative exons are identified, they are usually subjected to bioinformatics analysis in the same way as the gene ignoring the possibility of functionality change because of the alteration of domain caused by alternative exon. Here, we reveal an effective computational approach to explore an important molecular network based on potential changes of functionality induced by alternative exons obtained from our comprehensive analysis of neuronal cell differentiation.

RESULTS

From our previously identified 262 differentially alternatively spliced exons during neuronal cell differentiations, we extracted 241 sets that changed the amino acid sequences between the alternatively spliced sequences. Conserved domain searches indicated that annotated domain(s) were changed in 128 sets. We obtained 49 genes whose terms overlapped between domain description and gene annotation. Thus, these 49 genes have alternatively differentially spliced in exons that affect their main functions. We performed pathway analysis using these 49 genes and identified the EGFR (epidermal growth factor receptor) and mTOR (mammalian target of rapamycin) signaling pathway as being involved frequently. Recent studies reported that the mTOR pathway is associated with neuronal cell differentiation, vindicating that our approach extracted an important molecular network successfully.

CONCLUSIONS

Effective informatics approaches for exons should be more complex than those for genes, because changes in alternative exons affect protein functions via alterations of amino acid sequences and functional domains. Our method extracted alterations of functional domains and identified key alternative splicing events. We identified the EGFR and mTOR signaling pathway as the most affected pathway. The mTOR pathway is important for neuronal differentiation, suggesting that this in silico extraction of alternative splicing networks is useful. This preliminary analysis indicated that automated analysis of the effects of alternative splicing would provide a rich source of biologically relevant information.

The relationship between gene isoform multiplicity, number of exons and protein divergence

At present we know that phenotypic differences between organisms arise from a variety of sources, like protein sequence divergence, regulatory sequence divergence, alternative splicing, etc. However, we do not have yet a complete view of how these sources are related. Here we address this problem, studying the relationship between protein divergence and the ability of genes to express multiple isoforms. We used three genome-wide datasets of human-mouse orthologs to study the relationship between isoform multiplicity co-occurrence between orthologs (the fact that two orthologs have more than one isoform) and protein divergence. In all cases our results showed that there was a monotonic dependence between these two properties. We could explain this relationship in terms of a more fundamental one, between exon number of the largest isoform and protein divergence. We found that this last relationship was present, although with variations, in other species (chimpanzee, cow, rat, chicken, zebrafish and fruit fly). In summary, we have identified a relationship between protein divergence and isoform multiplicity co-occurrence and explained its origin in terms of a simple gene-level property. Finally, we discuss the biological implications of these findings for our understanding of inter-species phenotypic differences.