Biomedical impact of splicing mutations revealed through exome sequencing

Splicing is a cellular mechanism, which dictates eukaryotic gene expression by removing the noncoding introns and ligating the coding exons in the form of a messenger RNA molecule. Alternative splicing (AS) adds a major level of complexity to this mechanism and thus to the regulation of gene expression. This widespread cellular phenomenon generates multiple messenger RNA isoforms from a single gene, by utilizing alternative splice sites and promoting different exon-intron inclusions and exclusions. AS greatly increases the coding potential of eukaryotic genomes and hence contributes to the diversity of eukaryotic proteomes. Mutations that lead to disruptions of either constitutive splicing or AS cause several diseases, among which are myotonic dystrophy and cystic fibrosis. Aberrant splicing is also well established in cancer states. Identification of rare novel mutations associated with splice-site recognition, and splicing regulation in general, could provide further insight into genetic mechanisms of rare diseases. Here, disease relevance of aberrant splicing is reviewed, and the new methodological approach of starting from disease phenotype, employing exome sequencing and identifying rare mutations affecting splicing regulation is described. Exome sequencing has emerged as a reliable method for finding sequence variations associated with various disease states. To date, genetic studies using exome sequencing to find disease-causing mutations have focused on the discovery of nonsynonymous single nucleotide polymorphisms that alter amino acids or introduce early stop codons, or on the use of exome sequencing as a means to genotype known single nucleotide polymorphisms. The involvement of splicing mutations in inherited diseases has received little attention and thus likely occurs more frequently than currently estimated. Studies of exome sequencing followed by molecular and bioinformatic analyses have great potential to reveal the high impact of splicing mutations underlying human disease.     

Sequence, structure, and evolution of a complete human olfactory receptor gene cluster

The olfactory receptor (OR) gene cluster on human chromosome 17p13.3 was subjected to mixed shotgun automated DNA sequencing. The resulting 412 kb of genomic sequence include 17 OR coding regions, 6 of which are pseudogenes. Six of the coding regions were discovered only upon genomic sequencing, while the others were previously reported as partial sequences. A comparison of DNA sequences in the vicinity of the OR coding regions revealed a common gene structure with an intronless coding region and at least one upstream noncoding exon. Potential gene control regions including specific pyrimidine:purine tracts and Olf-1 sites have been identified. One of the pseudogenes apparently has evolved into a CpG island. Four extensive CpG islands can be discerned within the cluster, not coupled to specific OR genes. The cluster is flanked at its telomeric end by an unidentified open reading frame (C17orf2) with no significant similarity to any known protein. A high proportion of the cluster sequence (about 60%) belongs to various families of interspersed repetitive elements, with a clear predominance of LINE repeats. The OR genes in the cluster belong to two families and seven subfamilies, which show a relatively high degree of intermixing along the cluster, in seemingly random orientations. This genomic organization may be best accounted for by a complex series of evolutionary events.     

Reciprocal effects in anther cultures of wheat hybrids

This study was conducted to determine the reciprocal effects for anther culture response in wheat (Triticum aestivum L.) using a set of 4 × 4 full diallel crosses. Both reciprocal and nuclear genetic effects were highly significant for anther culture response and useful for selection and breeding purposes. General combining ability (GCA) effects were predominant for all investigated anther culture traits. Also, significant differences for specific combining ability (SCA) effects were detected between reciprocal crosses. Although significant reciprocal differences for responding anther, callus number and green plant regeneration were recorded in some reciprocal crosses, there were no significant reciprocal differences for albino plant regeneration. The use of one parent as male or female could lead to change at the production of green plants from the F₁ hybrids and screening of inbred lines for response to anther culture, without reciprocal effects, could decrease the utilization of breeding material.     

Mechanism of signal propagation upon retinal isomerization: insights from molecular dynamics simulations of rhodopsin restrained by normal modes

As one of the best studied members of the pharmaceutically relevant family of G-protein-coupled receptors, rhodopsin serves as a prototype for understanding the mechanism of G-protein-coupled receptor activation. Here, we aim at exploring functionally relevant conformational changes and signal transmission mechanisms involved in its photoactivation brought about through a cis-trans photoisomerization of retinal. For this exploration, we propose a molecular dynamics simulation protocol that utilizes normal modes derived from the anisotropic network model for proteins. Deformations along multiple low-frequency modes of motion are used to efficiently sample collective conformational changes in the presence of explicit membrane and water environment, consistent with interresidue interactions. We identify two highly stable regions in rhodopsin, one clustered near the chromophore, the other near the cytoplasmic ends of transmembrane helices H1, H2, and H7. Due to redistribution of interactions in the neighborhood of retinal upon stabilization of the trans form, local structural rearrangements in the adjoining H3-H6 residues are efficiently propagated to the cytoplasmic end of these particular helices. In the structures obtained by our simulations, all-trans retinal interacts with Cys¹⁶⁷ on H4 and Phe²⁰³ on H5, which were not accessible in the dark state, and exhibits stronger interactions with H5, while some of the contacts made (in the cis form) with H6 are lost.     

Vaccinia Virus Proteins A52 and B14 Share a Bcl-2–Like Fold but Have Evolved to Inhibit NF-κB rather than Apoptosis

Vaccinia virus (VACV), the prototype poxvirus, encodes numerous proteins that modulate the host response to infection. Two such proteins, B14 and A52, act inside infected cells to inhibit activation of NF-κB, thereby blocking the production of pro-inflammatory cytokines. We have solved the crystal structures of A52 and B14 at 1.9 Å and 2.7 Å resolution, respectively. Strikingly, both these proteins adopt a Bcl-2–like fold despite sharing no significant sequence similarity with other viral or cellular Bcl-2–like proteins. Unlike cellular and viral Bcl-2–like proteins described previously, A52 and B14 lack a surface groove for binding BH3 peptides from pro-apoptotic Bcl-2–like proteins and they do not modulate apoptosis. Structure-based phylogenetic analysis of 32 cellular and viral Bcl-2–like protein structures reveals that A52 and B14 are more closely related to each other and to VACV N1 and myxoma virus M11 than they are to other viral or cellular Bcl-2–like proteins. This suggests that a progenitor poxvirus acquired a gene encoding a Bcl-2–like protein and, over the course of evolution, gene duplication events have allowed the virus to exploit this Bcl-2 scaffold for interfering with distinct host signalling pathways.     

Structure and function of A41, a vaccinia virus chemokine binding protein

The vaccinia virus (VACV) A41L gene encodes a secreted 30 kDa glycoprotein that is nonessential for virus replication but affects the host response to infection. The A41 protein shares sequence similarity with another VACV protein that binds CC chemokines (called vCKBP, or viral CC chemokine inhibitor, vCCI), and strains of VACV lacking the A41L gene induced stronger CD8+ T-cell responses than control viruses expressing A41. Using surface plasmon resonance, we screened 39 human and murine chemokines and identified CCL21, CCL25, CCL26 and CCL28 as A41 ligands, with Kds of between 8 nM and 118 nM. Nonetheless, A41 was ineffective at inhibiting chemotaxis induced by these chemokines, indicating it did not block the interaction of these chemokines with their receptors. However the interaction of A41 and chemokines was inhibited in a dose-dependent manner by heparin, suggesting that A41 and heparin bind to overlapping sites on these chemokines. To better understand the mechanism of action of A41 its crystal structure was solved to 1.9 A resolution. The protein has a globular beta sandwich structure similar to that of the poxvirus vCCI family of proteins, but there are notable structural differences, particularly in surface loops and electrostatic charge distribution. Structural modelling suggests that the binding paradigm as defined for the vCCI-chemokine interaction is likely to be conserved between A41 and its chemokine partners. Additionally, sequence analysis of chemokines binding to A41 identified a signature for A41 binding. The biological and structural data suggest that A41 functions by forming moderately strong (nM) interactions with certain chemokines, sufficient to interfere with chemokine-glycosaminoglycan interactions at the cell surface (microM-nM) and thereby to destroy the chemokine concentration gradient, but not strong enough to disrupt the (pM) chemokine-chemokine receptor interactions.     

Analysis of correlated mutations in HIV-1 protease using spectral clustering

MOTIVATION: The ability of human immunodeficiency virus-1 (HIV-1) protease to develop mutations that confer multi-drug resistance (MDR) has been a major obstacle in designing rational therapies against HIV. Resistance is usually imparted by a cooperative mechanism that can be elucidated by a covariance analysis of sequence data. Identification of such correlated substitutions of amino acids may be obscured by evolutionary noise. RESULTS: HIV-1 protease sequences from patients subjected to different specific treatments (set 1), and from untreated patients (set 2) were subjected to sequence covariance analysis by evaluating the mutual information (MI) between all residue pairs. Spectral clustering of the resulting covariance matrices disclosed two distinctive clusters of correlated residues: the first, observed in set 1 but absent in set 2, contained residues involved in MDR acquisition; and the second, included those residues differentiated in the various HIV-1 protease subtypes, shortly referred to as the phylogenetic cluster. The MDR cluster occupies sites close to the central symmetry axis of the enzyme, which overlap with the global hinge region identified from coarse-grained normal-mode analysis of the enzyme structure. The phylogenetic cluster, on the other hand, occupies solvent-exposed and highly mobile regions. This study demonstrates (i) the possibility of distinguishing between the correlated substitutions resulting from neutral mutations and those induced by MDR upon appropriate clustering analysis of sequence covariance data and (ii) a connection between global dynamics and functional substitution of amino acids.     

Coupling between global dynamics and signal transduction pathways: a mechanism of allostery for chaperonin GroEL

Despite significant efforts toward understanding the molecular basis of allosteric communication, the mechanisms by which local energetic and conformational changes cooperatively diffuse from ligand-binding sites to distal regions across the 3-dimensional structure of allosteric proteins remain to be established. Recent experimental and theoretical evidence supports the view that allosteric communication is facilitated by the intrinsic ability of the biomolecules to undergo collective changes in structure, triggered by ligand binding. Two groups of studies recently proved to provide insights into such intrinsic, structure-induced effects: elastic network models that permit us to visualize the cooperative changes in conformation that are most readily accessible near native state conditions, and information-theoretic approaches that elucidate the most efficient pathways of signal transmission favored by the overall architecture. Using a combination of these two approaches, we highlight, by way of application to the bacterial chaperonin complex GroEL-GroES, how the most cooperative modes of motion play a role in mediating the propagation of allosteric signals. A functional coupling between the global dynamics sampled under equilibrium conditions and the signal transduction pathways inherently favored by network topology appears to control allosteric effects.