The estimation of statistical parameters for local alignment score distributions

The distribution of optimal local alignment scores of random sequences plays a vital role in evaluating the statistical significance of sequence alignments. These scores can be well described by an extreme-value distribution. The distribution’s parameters depend upon the scoring system employed and the random letter frequencies; in general they cannot be derived analytically, but must be estimated by curve fitting. For obtaining accurate parameter estimates, a form of the recently described ‘island’ method has several advantages. We describe this method in detail, and use it to investigate the functional dependence of these parameters on finite-length edge effects.

Creating local bibliographic databases: new tools for evidence-based health care*

The Internet has created new opportunities for librarians to present literature search results to clinicians. In order to take full advantage of these opportunities, libraries need to create locally maintained bibliographic databases. A simple method of creating a local bibliographic database and publishing it on the Web is described. The method uses off-the-shelf software and requires minimal programming. A hedge search strategy for outcome studies of clinical process interventions is created, and Ovid is used to search MEDLINE. The search results are saved and imported into EndNote libraries. The citations are modified, exported to a Microsoft Access database, and published on the Web. Clinicians can use a Web browser to search the database. The bibliographic database contains 13,803 MEDLINE citations of outcome studies. Most searches take between four and ten seconds and retrieve between ten and 100 citations. The entire cost of the software is under $900. Locally maintained bibliographic databases can be created easily and inexpensively. They significantly extend the evidence-based health care services that libraries can offer to clinicians.

A fast marching level set method for monotonically advancing fronts.

A fast marching level set method is presented for monotonically advancing fronts, which leads to an extremely fast scheme for solving the Eikonal equation. Level set methods are numerical techniques for computing the position of propagating fronts. They rely on an initial value partial differential equation for a propagating level set function and use techniques borrowed from hyperbolic conservation laws. Topological changes, corner and cusp development, and accurate determination of geometric properties such as curvature and normal direction are naturally obtained in this setting. This paper describes a particular case of such methods for interfaces whose speed depends only on local position. The technique works by coupling work on entropy conditions for interface motion, the theory of viscosity solutions for Hamilton-Jacobi equations, and fast adaptive narrow band level set methods. The technique is applicable to a variety of problems, including shape-from-shading problems, lithographic development calculations in microchip manufacturing, and arrival time problems in control theory.

cAMP compartmentation is responsible for a local activation of cardiac Ca2+ channels by beta-adrenergic agonists.

The role of cAMP subcellular compartmentation in the progress of beta-adrenergic stimulation of cardiac L-type calcium current (ICa) was investigated by using a method based on the use of whole-cell patch-clamp recording and a double capillary for extracellular microperfusion. Frog ventricular cells were sealed at both ends to two patch-clamp pipettes and positioned approximately halfway between the mouths of two capillaries that were separated by a 5-micron thin wall. ICa could be inhibited in one half or the other by omitting Ca2+ from one solution or the other. Exposing half of the cell to a saturating concentration of isoprenaline (ISO, 1 microM) produced a nonmaximal increase in ICa (347 +/- 70%; n = 4) since a subsequent application of ISO to the other part induced an additional effect of nearly similar amplitude to reach a 673 +/- 130% increase. However, half-cell exposure to forskolin (FSK, 30 microM) induced a maximal stimulation of ICa (561 +/- 55%; n = 4). This effect was not the result of adenylyl cyclase activation due to FSK diffusion in the nonexposed part of the cell. To determine the distant effects of ISO and FSK on ICa, the drugs were applied in a zero-Ca solution. Adding Ca2+ to the drug-containing solutions allowed us to record the local effect of the drugs. Dose-response curves for the local and distant effects of ISO and FSK on ICa were used as an index of cAMP concentration changes near the sarcolemma. We found that ISO induced a 40-fold, but FSK induced only a 4-fold, higher cAMP concentration close to the Ca2+ channels, in the part of the cell exposed to the drugs, than it did in the rest of the cell. cAMP compartmentation was greatly reduced after inhibition of phosphodiesterase activity with 3-isobutyl-methylxanthine, suggesting the colocalization of enzymes involved in the cAMP cascade. We conclude that beta-adrenergic receptors are functionally coupled to nearby Ca2+ channels via local elevations of cAMP.

Thermodynamic parameters for the helix-coil transition of oligopeptides: molecular dynamics simulation with the peptide growth method.

The helix-coil transition equilibrium of polypeptides in aqueous solution was studied by molecular dynamics simulation. The peptide growth simulation method was introduced to generate dynamic models of polypeptide chains in a statistical (random) coil or an alpha-helical conformation. The key element of this method is to build up a polypeptide chain during the course of a molecular transformation simulation, successively adding whole amino acid residues to the chain in a predefined conformation state (e.g., alpha-helical or statistical coil). Thus, oligopeptides of the same length and composition, but having different conformations, can be incrementally grown from a common precursor, and their relative conformational free energies can be calculated as the difference between the free energies for growing the individual peptides. This affords a straightforward calculation of the Zimm-Bragg sigma and s parameters for helix initiation and helix growth. The calculated sigma and s parameters for the polyalanine alpha-helix are in reasonable agreement with the experimental measurements. The peptide growth simulation method is an effective way to study quantitatively the thermodynamics of local protein folding.