MapDisto: fast and efficient computation of genetic linkage maps

Several options are available to the scientific community for genetic map construction but few are simple to install and use. Available programs either lack intuitive interface or are commercial, expensive for many laboratories. We present MapDisto, a free, user-friendly and powerful program for constructing genetic maps from experimental segregating populations. MapDisto is freely available at http://mapdisto.free.fr/DL/. Current version: 1.7.5.     

The logistic transform for bounded outcome scores

The logistic transformation, originally suggested by Johnson (1949), is applied to analyze responses that are restricted to a finite interval (e.g. (0,1)), so-called bounded outcome scores. Bounded outcome scores often have a non-standard distribution, e.g. J- or U-shaped, precluding classical parametric statistical approaches for analysis. Applying the logistic transformation on a normally distributed random variable, gives rise to a logit-normal (LN) distribution. This distribution can take a variety of shapes on (0,1). Further, the model can be extended to correct for (baseline) covariates. Therefore, the method could be useful for comparative clinical trials. Bounded outcomes can be found in many research areas, e.g. drug compliance research, quality-of-life studies, and pain (and pain relief) studies using visual analog scores, but all these scores can attain the boundary values 0 or 1. A natural extension of the above approach is therefore to assume a latent score on 0,1) having a LN distribution. Two cases are considered: (a) the bounded outcome score is a proportion where the true probabilities have a LN distribution on (0,1) and (b) the bounded outcome score on [0,1] is a coarsened version of a latent score with a LN distribution on (0,1). We also allow the variance (on the transformed scale) to depend on treatment. The usefulness of our approach for comparative clinical trials will be assessed in this paper. It turns out to be important to distinguish the case of equal and unequal variances. For a bounded outcome score of the second type and with equal variances, our approach comes close to ordinal probit (OP) regression. However, ignoring the inequality of variances can lead to highly biased parameter estimates. A simulation study compares the performance of our approach with the two-sample Wilcoxon test and with OP regression. Finally, the different methods are illustrated on two data sets.     

The conformational analysis of peptides using fourier transform IR spectroscopy

Fourier transform infrared spectroscopy (FTIR) can be used for conformational analysis of peptides in a wide range of environments. Measurements can be performed in aqueous solution, organic solvents, detergent micelles as well as in phospholipid membranes. Information on the secondary structure of peptides can be derived from the analysis of the strong amide I band. Orientation of secondary structural elements within a lipid bilayer matrix can be determined by means of polarized attenuated total reflectance–FTIR spectroscopy. Hydrogen–deuterium exchange can be monitored by the analysis of the, amide II band. This review gives some example of peptide systems studied by FTIR spectroscopy. Studies on alamethicin and α-aminoisobutyric acid containing peptides have shown that FTIR spectroscopy is a sensitive tool for identifying 3₁₀-helical structures. Changes in the structure of the magainins upon interaction with charged lipids were detected using FTIR spectroscopy. Tachyplesin is an example of a β-sheet containing membrane active peptide. Polarized ir spectroscopy reveals that the antiparallel β-sheet structures of tachyplesin are oriented parallel to the membrane surface. Synthesis of peptides corresponding to functionally/structurally important regions of large proteins is becoming increasingly popular. FTIR spectroscopy has been used to analyze the structure of synthetic peptides corresponding to the ion-selective pore of the voltage-gated potassium channel. In biomembrane systems these peptides adopt a highly helical structure. Under conditions, where these peptides are aggregated the presence of some intermolecular β-sheet structure can also be detected. © 1994 John Wiley & Sons, Inc.     

Interaction of dipalmitoylphosphatidylcholine and dimyristoylphosphatidylcholine-d54 mixtures with glycophorin. A fourier transform infrared investigation

Glycophorin from the human erythrocyte membrane has been isolated in pure form and reconstituted into large unilamellar vesicles comprised of binary mixtures of 1,2-dipalmitoyl-3-sn-phosphatidylcholine (DPPC) and chain perdeuterated 1,2-dimyristoyl-3-sn-phosphatidylcholine (DMPC-d54). The effect of temperature and protein on lipid structure and mixing was monitored by using Fourier transform infrared spectroscopy; deuteration of one of the components of the mixture permits observation of the protein interaction with each lipid species. The melting curves were analyzed by assuming that each lipid chain can exist in one of two physical states (i.e., gel or liquid crystalline), characterized by a temperature-dependent Lorentzian distribution for the line shape of the C-H or C-D stretching vibrations. The fraction of each lipid component melted at temperatures within the two-phase region of the phase diagram was calculated and approximate phase diagrams were constructed. Addition of protein lowers the liquidus line of the phase diagram while leaving the solidus line essentially unchanged. No lipid phase separation is observed. The effect of protein is more pronounced on the DPPC component than on the DMPC-d54. The former is significantly more disordered and/or fluidized at all lipid mole fractions in the ternary system than in the binary phospholipid mixture.     

Cluster analysis of protein fourier transform infrared spectra

Cluster analysis techniques were used to examine a set of Fourier transform infrared (FT-IR) spectra of bovine serum albumin (BSA) in the adsorbed and nonadsorbed states. The region from 1480 to 1600 cm^?1, comprising the amide II band, was used. Spectra were preprocessed to compensate for linear baseline variation, and the single linkage method of cluster analysis was applied. As expected, the spectra of adsorbed and nonadsorbed BSA fell into two distinct clusters. However, no further clustering was observed among the adsorbed BSA spectra on the basis of surface type, suggesting that surface specificity of the spectral changes induced in BSA by adsorption is not detectable above experimental variation. This work illustrates the value of using cluster analysis in the FT-IR study of proteins as a complement to other data analysis methods.     

Complexes of aluminium with peptide ligands: A fourier transform IR spectroscopic study

Aluminium has been recognized to be a neurotoxic agent and a risk factor in Alzheimer's disease and other neuronal dysfunctions. CD spectroscopic studies on two synthetic fragments of the human neurofilament protein midsized subunit (NF-M), and their alanine-for-serine-substituled and /or serine-phosphorylated derivatives showed the formation of stable, citric acid resistant complexes of Al³⁺with peptide ligands [M. Hollósi, Z.M. Shen, A. Perczel, and G.D. Fasman (1994) Proc. Natl. Acad. Sci. USA, vol. 9 , pp.4902-4906]. In the case of Ser-phosphorylated fragments, aβ-sheet inducing effect of Ca²⁺ and Al³⁺ ions was observed. However, the serine-containing parent peptides, NF-M 13 (KSPVPKSPVEEKG) and NF-M 17 (EEKGKSPVPKSPVEEKG), failed to show CD spectral changes reflecting β-sheet formation upon addition of Al³⁺ ions. On the basis of the amide I region of the Fourier transform ir spectra, in triftuoroethanol, the peptide backbone of NF-M17 and NF-M17 (A⁶A¹¹) shows marked changes in the presence ofAl³⁺. The most significant spectral differences are seen in the car-boxyl region (> 1700 cm^?l). The high-frequency component bands above 1760 cm^?1 in both spectra belong to the C? O of undissociated CF₃COOH. Another strong band at 1710 cm^?1 which appears only in the spectrum of NF-Ml 7 (A⁶A¹¹)(NF-M17 with Ser⁶ andSer¹¹ replaced by Ala) can be assigned to the side chain or C-terminal COOH groups. The differential proton-ation state of the carboxyl groups in the two peptides suggests the format ion ofAl³⁺ complexes of different structure and stability. The Al³⁺ complex ofNF-Ml 7 (A⁶A¹¹) is likely less stable, or one or more of the carboxylates are not coordinated to the Al³⁺ and thus can serve as a base to bind the liberated protons. In NF-M17 the OH groups of serines facilitate the formation of type [Al-pep(H_-1)] complexes with the involvement of all carboxylategroups in the molecule. The relevance of intramolecular and intermolecular Al³⁺ binding to the controversial biological role of aluminium is also discussed. © 1995 John Wiley & Sons, Inc.     

Tryptophan perturbation in the L intermediate of bacteriorhodopsin: fourier transform infrared analysis with indole-15N shift.

In the photoreaction of bacteriorhodopsin, the L intermediate shows an intense band at 3486 cm-1 which is unaffected by 2H2O (Maeda, A., Sasaki, J., Shichida, Y., & Yoshizawa, T. (1992) Biochemistry 31, 462-467]. This band is shifted to 3477 cm-1 by [indole-15N]tryptophan substitution and therefore is assigned to the N-H stretching vibration of the indole of tryptophan. Free indole in carbon tetrachloride shows its N-H stretching vibration at 3491 cm-1 [Fuson, N., Josien, M.-L., Powell, R. L., & Utterback, E. (1952) J. Chem. Phys. 20, 145-152]. Thus, it is suggested that at least one tryptophan residue in the L intermediate is not hydrogen bonded.     

fjoin: simple and efficient computation of feature overlaps.

Sets of biological features with genome coordinates (e.g., genes and promoters) are a particularly common form of data in bioinformatics today. Accordingly, an increasingly important processing step involves comparing coordinates from large sets of features to find overlapping feature pairs. This paper presents fjoin, an efficient, robust, and simple algorithm for finding these pairs, and a downloadable implementation. For typical bioinformatics feature sets, fjoin requires O(n log(n)) time (O(n) if the inputs are sorted) and uses O(1) space. The reference implementation is a stand-alone Python program; it implements the basic algorithm and a number of useful extensions, which are also discussed in this paper.     

The application of fourier transform infrared transmission spectroscopy to the study of model and natural membranes

Fourier transform infrared transmission spectroscopy is presented as a technique with great potential for the study of aqueous membrane preparations. The methodology of sample preparation, spectra recording and data reduction is outlined. Spectral parameters are derived from FT-IR difference spectra which provide an extremely sensitive means to monitor the temperature-dependent behavior of individual vibrational modes in model and natural membranes.     

Effect of anion binding on iodopsin studied by low-temperature fourier transform infrared spectroscopy.

The effect of anion binding on iodopsin, the chicken red-sensitive cone visual pigment, was studied by measurements of the Fourier transform infrared spectra of chloride- and nitrate-bound forms of iodopsin at 77 K. In addition to the blue shift of the absorption maximum upon substituting nitrate for chloride, the C=C stretching vibrations of iodopsin and its photoproducts were upshifted 5-6 cm(-)(1). The C=NH and C=ND stretching vibrations were the same in wavenumber between the chloride- and nitrate-bound forms, indicating that the binding of either chloride or nitrate has no effect on the interaction between the protonated Schiff base and the counterion. The vibrational bands of iodopsin in the fingerprint and the hydrogen out-of-plane wagging regions were insensitive to anion substitution, suggesting that local chromophore interactions with the anions are not crucial for the absorption spectral shift. In contrast, bathoiodopsin in the chloride-bound form exhibited an intense C(14)H wagging mode, whose intensity was considerably weakened upon substitution of nitrate for chloride. These results suggest that binding of chloride changes the environment near the C(14) position of the chromophore, which could be one of the factors in the thermal reverse reaction of bathoiodopsin to iodopsin in the chloride-bound form.     

On the computation of molecular surface correlations for protein docking using fourier techniques

The computation of surface correlations using a variety of molecular models has been applied to the unbound protein docking problem. Because of the computational complexity involved in examining all possible molecular orientations, the fast Fourier transform (FFT) (a fast numerical implementation of the discrete Fourier transform (DFT)) is generally applied to minimize the number of calculations. This approach is rooted in the convolution theorem which allows one to inverse transform the product of two DFTs in order to perform the correlation calculation. However, such a DFT calculation results in a cyclic or "circular" correlation which, in general, does not lead to the same result as the linear correlation desired for the docking problem. In this work, we provide computational bounds for constructing molecular models used in the molecular surface correlation problem. The derived bounds are then shown to be consistent with various intuitive guidelines previously reported in the protein docking literature. Finally, these bounds are applied to different molecular models in order to investigate their effect on the correlation calculation.     

REPuter: fast computation of maximal repeats in complete genomes.

SUMMARY: A software tool was implemented that computes exact repeats and palindromes in entire genomes very efficiently. AVAILABILITY: Via the Bielefeld Bioinformatics Server (http://bibiserv.techfak.uni-bielefeld.de/rep uter/).     

Lipid reorganization in biological membranes. A study by fourier transform infrared difference spectroscopy

The first application of infrared difference spectroscopy to the study of a natural biological membrane is described. Perdeuterated palmitic acid was incorporated biosynthetically into the lipids of the plasma membrane of Acholeplasma laidlawii and the temperature-induced structural rearrangement of the endogenous lipids monitored via their C²H vibrational modes. Changes in infrared parameters were studied between 0 and 50°C and contrasted with those occurring in the model membrane system of $\text{1,2-}\text{diperdeuteropalmitoyl}\text{-sn-}\text{glycero-3-phosphocholine}$. The phase transition of the biomembrane occurs over a 20°C range with the temperature of the maximum rate of change of absorbance coinciding with that of the sharp phase transition of the model membrane.     

A note on efficient computation of haplotypes via perfect phylogeny.

The problem of inferring haplotype phase from a population of genotypes has received a lot of attention recently. This is partly due to the observation that there are many regions on human genomic DNA where genetic recombination is rare (Helmuth, 2001; Daly et al., 2001; Stephens et al., 2001; Friss et al., 2001). A Haplotype Map project has been announced by NIH to identify and characterize populations in terms of these haplotypes. Recently, Gusfield introduced the perfect phylogeny haplotyping problem, as an algorithmic implication of the no-recombination in long blocks observation, together with the standard population-genetic assumption of infinite sites. Gusfield's solution based on matroid theory was followed by direct theta(nm2) solutions that use simpler techniques (Bafna et al., 2003; Eskin et al., 2003), and also bound the number of solutions to the PPH problem. In this short note, we address two questions that were left open. First, can the algorithms of Bafna et al. (2003) and Eskin et al. (2003) be sped-up to O(nm + m2) time, which would imply an O(nm) time-bound for the PPH problem? Second, if there are multiple solutions, can we find one that is most parsimonious in terms of the number of distinct haplotypes. We give reductions that suggests that the answer to both questions is "no." For the first problem, we show that computing the output of the first step (in either method) is equivalent to Boolean matrix multiplication. Therefore, the best bound we can presently achieve is O(nm(omega-1)), where omega < or = 2.52 is the exponent of matrix multiplication. Thus, any linear time solution to the PPH problem likely requires a different approach. For the second problem of computing a PPH solution that minimizes the number of distinct haplotypes, we show that the problem is NP-hard using a reduction from Vertex Cover (Garey and Johnson, 1979).     

Low-frequency fourier transform infrared spectroscopy of the oxygen-evolving complex in Photosystem II*

In this communication, we report our progress on the development of low-frequency Fourier transform infrared (FTIR) spectroscopic techniques to study metal-substrate and metal-ligand vibrational modes in the Photosystem II/oxygen-evolving complex (PS II/OEC). This information will provide important structural and mechanistic insight into the OEC. Strong water absorption in the low-frequency region (below 1000 cm⁻¹), a lack of suitable materials, and temperature control problems have limited previous FTIR spectroscopic studies of the OEC to higher frequencies (>1000 cm⁻¹). We have overcome these technical difficulties that have blocked access to the low-frequency region and have developed successive instruments that allow us to move deeper into the low-frequency region (down to 350 cm⁻¹), while increasing both data accumulation efficiency and S/N ratio. We have detected several low-frequency modes in the S₂/S₁spectrum that are specifically associated with these two states. Our results demonstrate the utility of FTIR techniques in accessing low-frequency modes in Photosystem II and in proteins generally. This revised version was published online in June 2006 with corrections to the Cover Date.