Algorithms and software for support of gene identification experiments

MOTIVATION: Gene annotation is the final goal of gene prediction algorithms. However, these algorithms frequently make mistakes and therefore the use of gene predictions for sequence annotation is hardly possible. As a result, biologists are forced to conduct time-consuming gene identification experiments by designing appropriate PCR primers to test cDNA libraries or applying RT-PCR, exon trapping/amplification, or other techniques. This process frequently amounts to 'guessing' PCR primers on top of unreliable gene predictions and frequently leads to wasting of experimental efforts. RESULTS: The present paper proposes a simple and reliable algorithm for experimental gene identification which bypasses the unreliable gene prediction step. Studies of the performance of the algorithm on a sample of human genes indicate that an experimental protocol based on the algorithm's predictions achieves an accurate gene identification with relatively few PCR primers. Predictions of PCR primers may be used for exon amplification in preliminary mutation analysis during an attempt to identify a gene responsible for a disease. We propose a simple approach to find a short region from a genomic sequence that with high probability overlaps with some exon of the gene. The algorithm is enhanced to find one or more segments that are probably contained in the translated region of the gene and can be used as PCR primers to select appropriate clones in cDNA libraries by selective amplification. The algorithm is further extended to locate a set of PCR primers that uniformly cover all translated regions and can be used for RT-PCR and further sequencing of (unknown) mRNA.      

Molecular evolution of a repetitive region within the per gene of Drosophila

The clock gene period (per) controls a number of biological rhythms in Drosophila. In D. melanogaster, per has a repetitive region that encodes a number of alternating threonine-glycine residues. We sequenced and compared this region from several different Drosophila species belonging to various groups within the Drosophila and Sophophora subgenera. This part of per shows a great variability in both DNA sequence and length. Furthermore, analysis of the data suggests that changes in the length of this variable region might be associated with amino acid replacements in the more conserved flanking sequences.      

Anatomical and chemical characteristics of the seed coat of <i>Medicago sativa</i> L. (alfalfa) cv. Baralfa 85 seeds and their association with seed dormancy

Seeds of alfalfa (Medicago sativa L.) can exhibit seedcoat imposed dormancy, which produces hard seeds within a seed lot. These seeds do not germinate because they do not imbibe water due to a barrier to water entry in the seed coat. The aim of this work was to analyze the anatomical and chemical characteristics of the testa of alfalfa seeds with respect to water permeability levels. The anatomy of seeds of the cv. Baralfa 85 was studied and structural substances, polyphenols, tannins and cutin present in the testa of seeds of different water permeability levels were determined. The anatomical characteristics of the seed coat and the proportions of components were found to determine the permeability level of the seed coat, an aspect that is associated with the physical seed dormancy level. Anatomically, increased thickness of the testa was associated with a lower permeability level. The difference may be attributed to the variation in cuticle thickness, length of macrosclereids and thickness of the cell wall, and presence and development of osteosclereids. From the physiological and chemical points of view, the mechanism of physical dormancy of the testa is explained by a greater amount of components that repel water and cement the cell wall, such as polyphenols, lignins, condensed tannins, pectic substances, and a lower proportion of cellulose and hemicellulose.     

The genetic architecture of fitness in a seed beetle: assessing the potential for indirect genetic benefits of female choice

Background  

Quantifying the amount of standing genetic variation in fitness represents an empirical challenge. Unfortunately, the shortage of detailed studies of the genetic architecture of fitness has hampered progress in several domains of evolutionary biology. One such area is the study of sexual selection. In particular, the evolution of adaptive female choice by indirect genetic benefits relies on the presence of genetic variation for fitness. Female choice by genetic benefits fall broadly into good genes (additive) models and compatibility (non-additive) models where the strength of selection is dictated by the genetic architecture of fitness. To characterize the genetic architecture of fitness, we employed a quantitative genetic design (the diallel cross) in a population of the seed beetle Callosobruchus maculatus, which is known to exhibit post-copulatory female choice. From reciprocal crosses of inbred lines, we assayed egg production, egg-to-adult survival, and lifetime offspring production of the outbred F1 daughters (F1 productivity).     

Mapping the nucleotide and isoform-dependent structural and dynamical features of Ras proteins

Ras GTPases are conformational switches controlling cell proliferation, differentiation, and development. Despite their prominent role in many forms of cancer, the mechanism of conformational transition between inactive GDP-bound and active GTP-bound states remains unclear. Here we describe a detailed analysis of available experimental structures and molecular dynamics simulations to quantitatively assess the structural and dynamical features of active and inactive states and their interconversion. We demonstrate that GTP-bound and nucleotide-free G12V H-ras sample a wide region of conformational space, and show that the inactive-to-active transition is a multiphase process defined by the relative rearrangement of the two switches and the orientation of Tyr32. We also modeled and simulated N- and K-ras proteins and found that K-ras is more flexible than N- and H-ras. We identified a number of isoform-specific, long-range side chain interactions that define unique pathways of communication between the nucleotide binding site and the C terminus.     

Lipid-Sorting Specificity Encoded in K-Ras Membrane Anchor Regulates Signal Output

1. Download high-res image (273KB)
2. Download full-size image

     

Alternative splicing and protein function

Background  

Alternative splicing is a major mechanism of generating protein diversity in higher eukaryotes. Although at least half, and probably more, of mammalian genes are alternatively spliced, it was not clear, whether the frequency of alternative splicing is the same in different functional categories. The problem is obscured by uneven coverage of genes by ESTs and a large number of artifacts in the EST data.     

Feedback equilibrium control during human standing

Equilibrium maintenance during standing in humans was investigated with a 3-joint (ankle, knee and hip) sagittal model of body movement. The experimental paradigm consisted of sudden perturbations of humans in quiet stance by backward displacements of the support platform. Data analysis was performed using eigenvectors of motion equation. The results supported three conclusions. First, independent feedback control of movements along eigenvectors (eigenmovements) can adequately describe human postural responses to stance perturbations. This conclusion is consistent with previous observations (Alexandrov et al., 2001b) that these same eigenmovements are also independently controlled in a feed-forward manner during voluntary upper-trunk bending. Second, independent feedback control of each eigenmovement is sufficient to provide its stability. Third, the feedback loop in each eigenmovement can be modeled as a linear visco-elastic spring with delay. Visco-elastic parameters and time-delay values result from the combined contribution of passive visco-elastic mechanisms and sensory systems of different modalities     

Salt bridges destabilize a leucine zipper designed for maximized ion pairing between helices

Interhelical salt bridges are common in leucine zippers and are thought to stabilize the coiled coil conformation. Here we present a detailed thermodynamic investigation of the designed, disulfide-linked leucine zipper AB(SS) whose high-resolution NMR structure shows six interhelical ion pairs between heptad positions g of one helix and e' of the other helix but no ion pairing within single helices. The average pK(a) value of the Glu side chain carboxyl groups of AB(SS) is slightly higher than the pK(a) of a freely accessible Glu in an unfolded peptide [Marti, D. N., Jelesarov, I., and Bosshard, H. R. (2000) Biochemistry 39, 12804-12818]. This indicates that the salt bridges are destabilizing, a prediction we now have confirmed by determining the pH +/- stability profile of AB(SS). Circular dichroism-monitored unfolding by urea and by heating and differential scanning calorimetry show that the coiled coil conformation is approximately 5 kJ/mol more stable when salt bridges are broken by protonation of the carboxyl side chains. Using guanidinium chloride as the denaturant, the increase in the free energy of unfolding on protonation of the carboxyl side chains is larger, approximately 17 kJ/mol. The discrepancy between urea and guanidinium chloride unfolding can be ascribed to the ionic nature of guanidinium chloride, which screens charge-charge interactions. This work demonstrates the difficulty of predicting the energetic contribution of salt bridges from structural data alone even in a case where the ion pairs are seen in high-resolution NMR structures. The reason is that the contribution to stability results from a fine balance between energetically favorable Coulombic attractions and unfavorable desolvation of charges and conformational constraints of the residues involved in ion pairing. The apparent discrepancy between the results presented here and mutational studies indicating stabilization by salt bridges is discussed and resolved. An explanation is proposed for why interhelical salt bridges are frequently found in natural coiled coils despite evidence that they do not directly contribute to stability.     

The PINE study: rationale and design of a randomised comparison of epidural injection of local anaesthetics and steroids versus care-as-usual to prevent postherpetic neuralgia in the elderly [ISRCTN32866390]

BACKGROUND: Postherpetic neuralgia (PHN) is by far the most common complication of herpes zoster (HZ) and one of the most intractable pain disorders. Since PHN is seen most often in the elderly, the number of patients with this disorder is expected to increase in our ageing society. PHN may last for months to years and has a high impact on the quality of life. The results of PHN treatment are rather disappointing. Epidural injection of local anaesthetics and steroids in the acute phase of HZ is a promising therapy for the prevention of PHN. Since randomised trials on the effectiveness of this intervention are lacking, the PINE (Prevention by epidural Injection of postherpetic Neuralgia in the Elderly) study was set up. The PINE study compares the effectiveness and cost-effectiveness of a single epidural injection of local anaesthetics and steroids during the acute phase of HZ with that of care-as-usual (i.e. antivirals and analgesics) in preventing PHN in elderly patients. METHODS / DESIGN: The PINE study is an open, multicenter clinical trial in which 550 elderly (age >/= 50 yr.) patients who consult their general practitioner in the acute phase of HZ (rash < 7 days) are randomised to one of the treatment groups. The primary clinical endpoint is the presence of HZ-related pain one month after the onset of the rash. Secondary endpoints include duration and severity of pain, re-interventions aiming to treat the existing pain, side effects, quality of life, and cost-effectiveness. CONCLUSION: The PINE study is aimed to quantify the (cost-) effectiveness of a single epidural injection during the acute phase of HZ on the prevention of PHN.