Comparison of intron-containing and intron-lacking human genes elucidates putative exonic splicing enhancers.

Of the rules used by the splicing machinery to precisely determine intron-exon boundaries only a fraction is known. Recent evidence suggests that specific short sequences within exons help in defining these boundaries. Such sequences are known as exonic splicing enhancers (ESE). A possible bioinformatical approach to studying ESE sequences is to compare genes that harbor introns with genes that do not. For this purpose two non-redundant samples of 719 intron-containing and 63 intron-lacking human genes were created. We performed a statistical analysis on these datasets of intron-containing and intron-lacking human coding sequences and found a statistically significant difference (P = 0.01) between these samples in terms of 5-6mer oligonucleotide distributions. The difference is not created by a few strong signals present in the majority of exons, but rather by the accumulation of multiple weak signals through small variations in codon frequencies, codon biases and context-dependent codon biases between the samples. A list of putative novel human splicing regulation sequences has been elucidated by our analysis.

The genomic structure of c14orf1 is conserved across eukarya.

We have recently cloned the gene C14orf1, which is strongly expressed in normal testis and in several cancer cell lines and tumors. This gene maps to 14q24.3 and is interrupted by four introns. Two of them are also represented in the open reading frame of Schizosaccharomyces pombe in the same phase. In Arabidopsis taliana only the first of the two introns was found, in the same phase as the corresponding ones in S. pombe and human. Disruption of the ortholog in Saccharomyces cerevisiae (Yer044c) led to a severe growth defect, and C14orf1 failed to complement mutant yeast when put under the control of the natural Yer044c promoter. Further studies are needed to understand the causes underlying the high degree of conservation of the C14orf1 genomic structure.

Structure of the Sm binding site from human U4 snRNA derived from a 3 ns PME molecular dynamics simulation.

A molecular dynamics simulation of the Sm binding site from human U4 snRNA was undertaken to determine the conformational flexibility of this region and to identify RNA conformations that were important for binding of the Sm proteins. The RNA was fully-solvated (>9,000 water molecules) and charge neutralized by inclusion of potassium ions. A three nanosecond MD simulation was conducted using AMBER with long-range electrostatic forces considered using the particle mesh Ewald summation method. The initial model of the Sm binding site region had the central and 3' stem-loops that flanked the Sm site co-axial with one another, and with the single-stranded Sm binding site region ([I] conformation). During the course of the trajectory, the axes of the 3' stem-loop, and later the central stem-loop, became roughly orthogonal from their original anti-parallel orientation. As these conformational changes occurred, the snRNA adopted first an [L] conformation, and finally a [U] conformation. The [U] conformation was more stable than either the [I] or [L] conformations, and persisted for the final 1 ns of the trajectory. Analysis of the structure resulting from the MD simulations revealed the bulged nucleotide, U114, and the mismatched Ag91-G110 base pair provided distinctive structural features that may enhance Sm protein binding. Based on the results of the MD simulation and the available experimental data, we proposed a mechanism for the binding of the Sm protein sub-complexes to the snRNA. In this model, the D1/D2 and E/F/G Sm protein sub-complexes first bind the snRNA in the [U] conformation, followed by conformational re-arrangement to the [I] conformation and binding of the D3/B Sm protein sub-complex.

EID: the Exon-Intron Database-an exhaustive database of protein-coding intron-containing genes.

To aid studies of molecular evolution and to assist in gene prediction research, we have constructed an Exon-Intron Database (EID) in FASTA format. Currently, the database is derived from GenBank release 112, and it contains 51 289 protein-coding genes (287 209 exons) that harbor introns, along with extensive descriptions of each gene and its DNA and protein sequences, as well as splice motif information. There is 17% redundancy inherited from GenBank-a purge at the 99% identity level reduced the database to 42 460 genes (243 589 exons). We have created subdatabases of genes whose intron positions have been experimentally determined. One such database, constructed by comparing genomic and mRNA sequences, contains 11 242 genes (62 474 exons). A larger database of 22 196 genes (105 595 exons) was constructed by selecting on keywords to eliminate computer-predicted genes. By examining the two nucleotides adjacent to the intron boundary, we infer that there is a 2% rate of errors or other deviations from the standard GTellipsisAG motif in nuclear genes. This criterion can be used to eliminate 4921 genes from the overall database. Various tools are provided to enable generation of user-specific subsets of the EID. The EID distribution can be obtained from http://mcb.harvard.edu/gilbert/EID

Pathways of electron transfer in Escherichia coli DNA photolyase: Trp306 to FADH.

We describe the results of a series of theoretical calculations of electron transfer pathways between Trp306 and *FADH. in the Escherichia coli DNA photolyase molecule, using the method of interatomic tunneling currents. It is found that there are two conformationally orthogonal tryptophans, Trp359 and Trp382, between donor and acceptor that play a crucial role in the pathways of the electron transfer process. The pathways depend vitally on the aromaticity of tryptophans and the flavin molecule. The results of this calculation suggest that the major pathway of the electron transfer is due to a set of overlapping orthogonal pi-rings, which starts from the donor Trp306, runs through Trp359 and Trp382, and finally reaches the flavin group of the acceptor complex, FADH.

Effect of protein dynamics on biological electron transfer.

Computer simulations of the effect of protein dynamics on the long distance tunneling mediated by the protein matrix have been carried out for a Ru-modified (His 126) azurin molecule. We find that the tunneling matrix element is a sensitive function of the atomic configuration of the part of the protein matrix in which tunneling currents (pathways) are localized. Molecular dynamics simulations show that fluctuations of the matrix element can occur on a time scale as short as 10 fs. These short time fluctuations are an indication of a strong dynamic coupling of a tunneling electron to vibrational motions of the protein nuclear coordinates. The latter results in a modification of the conventional Marcus picture of electron transfer in proteins. The new element in the modified theory is that the tunneling electron is capable of emitting or absorbing vibrational energy (phonons) from the medium. As a result, some biological reactions may occur in an activationless fashion. An analytical theoretical model is proposed to account for thermal fluctuations of the medium in long distance electron transfer reactions. The model shows that, at long distances, the phonon-modified inelastic tunneling always dominates over the conventional elastic tunneling.