Model-based fitting of single-channel dwell-time distributions

Single-channel recordings provide unprecedented resolutions on kinetics of conformational changes of ion channels. Several approaches exist for analysis of the data, including the dwell-time histogram fittings and the full maximal-likelihood approaches that fit either the idealized dwell-time sequence or more ambitiously the noisy data directly using hidden Markov modeling. Although the full maximum likelihood approaches are statistically advantageous, they can be time-consuming especially for large datasets and/or complex models. We present here an alternative approach for model-based fitting of one-dimensional and two-dimensional dwell-time histograms. To improve performance, we derived analytical expressions for the derivatives of one-dimensional and two-dimensional dwell-time distribution functions and employed the gradient-based variable metric method for fast search of optimal rate constants in a model. The algorithm also has the ability to allow for a first-order correction for the effects of missed events, global fitting across different experimental conditions, and imposition of typical constraints on rate constants including microscopic reversibility. Numerical examples are presented to illustrate the performance of the algorithm, and comparisons with the full maximum likelihood fitting are discussed.     

A general solution to the time interval omission problem applied to single channel analysis.

To obtain the open or closed time interval distributions of patch clamp signals, several workers have used a half-amplitude minimum time interval criterion. Within this framework, no transition between states of different conductance levels is considered to have taken place if it leads to a time interval smaller than a certain critical value. This procedure modifies substantially the open or closed time interval distribution of the random signal to be analyzed, since time intervals well above the time resolution of the recording system may be interrupted by short gaps that may or may not satisfy the minimum time interval criterion. We present here a general theoretical framework by means of which the effect of time interval omission on time interval distributions can be taken into account. Based on the mathematical formalism provided by the Kolmogorov forward equation, special matrix operators are first defined. The general solution to the time omission problem in its integral form is then derived. In view of the poor computational feasibility of the resulting solution, a first-order approximation is also presented. This approximation consists essentially in neglecting the contribution of the undetected gaps to the total length of the resulting time interval. The exact and approximate solutions are then applied to two special kinetic schemes commonly found in single-channel studies, namely the O-C and C-O-C models. The applicability of the proposed formalism to the time interval distribution problem of a damped random signal is finally discussed.     

Regulation by phosphorylation. Yet another twist in the WASP story

Cell migration and other complex cellular processes involve a variety of signaling molecules and require the integration of multiple signals into a coherent cytoskeletal response. Two papers in the May issue of Molecular Cell now demonstrate that phosphorylation plays a critical role in WASP function as a regulator of Arp2/3-mediated actin polymerization.     

Yet another numbering scheme for immunoglobulin variable domains: an automatic modeling and analysis tool

A common residue numbering scheme for all immunoglobulin variable domains (immunoglobulin light chain lambda (V(lambda)) and kappa (V(kappa)) variable domains, heavy chain variable domains (V(H)) and T-cell receptor alpha (V(alpha)), beta (V(beta)), gamma (V(gamma)) and delta (V(delta)) variable domains) has been devised. Based on the spatial alignment of known three-dimensional structures of immunoglobulin domains, it places the alignment gaps in a way that minimizes the average deviation from the averaged structure of the aligned domains. This residue numbering scheme was applied to the immunoglobulin variable domain structures in the PDB database to automate the extraction of information on structural variations in homologous positions of the different molecules. A number of methods are presented that allow the automated projection of information derived from individual structures or from the comparison of multi-structure alignments onto a graphical representation of the sequence alignment.     

Yet another improved silver staining method for the detection of proteins in polyacrylamide gels

Silver staining is very sensitive for detection of proteins in polyacrylamide gels and different procedures have been published. By combining and modifying some of the recipes, a very reproducible method, which is based upon staining with diamine complexes of silver, has been developed. The background staining is negligible and reduced silver does not precipitate on the gel surface. The technique works very well for sodium dodecyl sulfate-polyacrylamide gel electrophoresis in both homogeneous and in gradient gels as well as for two-dimensional (2-D) PAGE. It was possible to detect 1-10 ng of protein corresponding to approximately 50 pg/mm2, provided that a discontinuous buffer system was used, which gives sharp bands.     

A new approach to the species classification problem in floristic analysis

A new method of species (inverse) classification of vegetation data, i.e. classification of species into groups with similar ecological tolerances, is presented which overcomes the problems of species abundance distorting the results. The algorithm TWO-STEP is based on the use of an asymmetric measure of dissimilarity: where i, j are species, h is the stand, n is the total number of stands, and x_ih is the amount of species i in stand h. The algorithm uses the rows of the asymmetric dissimilarity matrix generated as above to form a second symmetric dissimilarity matrix using the measure: where m is the number of species and k the species. Flexible sorting is applied to generate a species classification. Comparison of results after applying the TWO-STEP algorithm and a standard alternative to an artificial data set demonstrates its efficacy. TWO-STEP also shows considerable advantages over previous analyses for a Queensland rainforest data set (quantitative) and an English heath (qualitative) data set. Normalization of species data appears advantageous for quantitative data only.     

Optical single-channel analysis of the aerolysin pore in erythrocyte membranes.

Scanning microphotolysis (Scamp), a recently developed photobleaching technique, was used to analyze the transport of two small organic anions and one inorganic cation through single pores formed in human erythrocyte membranes by the channel-forming toxin aerolysin secreted by Aeromonas species. The transport rate constants of erythrocyte ghosts carrying a single aerolysin pore were determined to be (1.83 +/- 0.43) x 10(-3) s-1 for Lucifer yellow, (0.33 +/- 0.10) x 10(-3) s-1 for carboxyfluorescein, and (8.20 +/- 2.30) x 10(-3) s-1 for Ca2+. The radius of the aerolysin pore was derived from the rate constants to be 19-23 A, taking steric hindrance and viscous drag into account. The size of the Ca2+ rate constant implies that at physiological extracellular Ca2+ concentrations (> 1 mM) the intracellular Ca2+ concentration would be elevated to the critical level of > 1 microM in much less than a second after formation of a single aerolysin pore in the plasma membrane. Thus changes in the levels of Ca2+ or other critical intracellular components may be more likely to cause cell death than osmotic imbalance.     

Regulation of USP7/HAUSP in response to DNA damage: Yet another role for ATM

Comment on: Khoronenkova SV, et al. Mol Cell 2012; 45:801-13.     

Regulation of USP7/HAUSP in response to DNA damage: Yet another role for ATM

Comment on: Khoronenkova SV, et al. Mol Cell 2012; 45:801-13.     

Resistance to levamisole resolved at the single-channel level.

Levamisole is commonly used to treat nematode parasite infections but therapy is limited by resistance. The purpose of this study was to determine the mechanism of resistance to this selective nicotinic drug. Levamisole receptor channel currents in muscle patches from levamisole-sensitive and levamisole-resistant isolates of the parasitic nematode Oesophagostomum dentatum were compared. The number of channels present in patches of sensitive and resistant isolates was similar at 10 microM levamisole, but at 30 microM and 100 microM the resistant isolate contained fewer active patches, suggesting desensitization. Mean Po and open times were reduced in resistant isolates. The distribution of conductances of channels in the sensitive isolate revealed a heterogeneous receptor population and the presence of G25, G35, G40, and G45 subtypes. A G35 subtype was missing in the resistant isolate. Resistance to levamisole was produced by changes in the averaged properties of the levamisole receptor population, with some receptors from sensitive and resistant isolates having indistinguishable characteristics.     

Topology fingerprint approach to the inverse protein folding problem.

We describe the most general solution to date of the problem of matching globular protein sequences to the appropriate three-dimensional structures. The screening template, against which sequences are tested, is provided by a protein "structural fingerprint" library based on the contact map and the buried/exposed pattern of residues. Then, a lattice Monte Carlo algorithm validates or dismisses the stability of the proposed fold. Examples of known structural similarities between proteins having weakly or unrelated sequences such as the globins and phycocyanins, the eight-member alpha/beta fold of triose phosphate isomerase and even a close structural equivalence between azurin and immunoglobulins are found.