Imaging therapeutic peptide transport across intestinal barriers

Oral delivery is a highly preferred method for drug administration due to high patient compliance. However, oral administration is intrinsically challenging for pharmacologically interesting drug classes, in particular pharmaceutical peptides, due to the biological barriers associated with the gastrointestinal tract. In this review, we start by summarizing the pharmacological performance of several clinically relevant orally administrated therapeutic peptides, highlighting their low bioavailabilities. Thus, there is a strong need to increase the transport of peptide drugs across the intestinal barrier to realize future treatment needs and further development in the field. Currently, progress is hampered by a lack of understanding of transport mechanisms that govern intestinal absorption and transport of peptide drugs, including the effects of the permeability enhancers commonly used to mediate uptake. We describe how, for the past decades, mechanistic insights have predominantly been gained using functional assays with end-point read-out capabilities, which only allow indirect study of peptide transport mechanisms. We then focus on fluorescence imaging that, on the other hand, provides opportunities to directly visualize and thus follow peptide transport at high spatiotemporal resolution. Consequently, it may provide new and detailed mechanistic understanding of the interplay between the physicochemical properties of peptides and cellular processes; an interplay that determines the efficiency of transport. We review current methodology and state of the art in the field of fluorescence imaging to study intestinal barrier transport of peptides, and provide a comprehensive overview of the imaging-compatible in vitro, ex vivo, and in vivo platforms that currently are being developed to accelerate this emerging field of research.

Single-particle trajectories reveal two-state diffusion-kinetics of hOGG1 proteins on DNA

We reanalyze trajectories of hOGG1 repair proteins diffusing on DNA. A previous analysis of these trajectories with the popular mean-squared-displacement approach revealed only simple diffusion. Here, a new optimal estimator of diffusion coefficients reveals two-state kinetics of the protein. A simple, solvable model, in which the protein randomly switches between a loosely bound, highly mobile state and a tightly bound, less mobile state is the simplest possible dynamic model consistent with the data. It yields accurate estimates of hOGG1’s (i) diffusivity in each state, uncorrupted by experimental errors arising from shot noise, motion blur and thermal fluctuations of the DNA; (ii) rates of switching between states and (iii) rate of detachment from the DNA. The protein spends roughly equal time in each state. It detaches only from the loosely bound state, with a rate that depends on pH and the salt concentration in solution, while its rates for switching between states are insensitive to both. The diffusivity in the loosely bound state depends primarily on pH and is three to ten times higher than in the tightly bound state. We propose and discuss some new experiments that take full advantage of the new tools of analysis presented here.

How To Characterize Individual Nanosize Liposomes with Simple Self-Calibrating Fluorescence Microscopy

Nanosize lipid vesicles are used extensively at the interface between nanotechnology and biology, e.g., as containers for chemical reactions at minute concentrations and vehicles for targeted delivery of pharmaceuticals. Typically, vesicle samples are heterogeneous as regards vesicle size and structural properties. Consequently, vesicles must be characterized individually to ensure correct interpretation of experimental results. Here we do that using dual-color fluorescence labeling of vesicles-of their lipid bilayers and lumens, separately. A vesicle then images as two spots, one in each color channel. A simple image analysis determines the total intensity and width of each spot. These four data all depend on the vesicle radius in a simple manner for vesicles that are spherical, unilamellar, and optimal encapsulators of molecular cargo. This permits identification of such ideal vesicles. They in turn enable calibration of the dual-color fluorescence microscopy images they appear in. Since this calibration is not a separate experiment but an analysis of images of vesicles to be characterized, it eliminates the potential source of error that a separate calibration experiment would have been. Nonideal vesicles in the same images were characterized by how their four data violate the calibrated relationship established for ideal vesicles. In this way, our method yields size, shape, lamellarity, and encapsulation efficiency of each imaged vesicle. Applying this procedure to extruded samples of vesicles, we found that, contrary to common assumptions, only a fraction of vesicles are ideal.

Enrichment of megabase-sized DNA molecules for single-molecule optical mapping and next-generation sequencing

Next-generation sequencing (NGS) has caused a revolution, yet left a gap: long-range genetic information from native, non-amplified DNA fragments is unavailable. It might be obtained by optical mapping of megabase-sized DNA molecules. Frequently only a specific genomic region is of interest, so here we introduce a method for selection and enrichment of megabase-sized DNA molecules intended for single-molecule optical mapping: DNA from a human cell line is digested by the NotI rare-cutting enzyme and size-selected by pulsed-field gel electrophoresis. For demonstration, more than 600 sub-megabase- to megabase-sized DNA molecules were recovered from the gel and analysed by denaturation-renaturation optical mapping. Size-selected molecules from the same gel were sequenced by NGS. The optically mapped molecules and the NGS reads showed enrichment from regions defined by NotI restriction sites. We demonstrate that the unannotated genome can be characterized in a locus-specific manner via molecules partially overlapping with the annotated genome. The method is a promising tool for investigation of structural variants in enriched human genomic regions for both research and diagnostic purposes. Our enrichment method could potentially work with other genomes or target specified regions by applying other genomic editing tools, such as the CRISPR/Cas9 system.

Classification of DNA nucleotides with transverse tunneling currents

It has been theoretically suggested and experimentally demonstrated that fast and low-cost sequencing of DNA, RNA, and peptide molecules might be achieved by passing such molecules between electrodes embedded in a nanochannel. The experimental realization of this scheme faces major challenges, however. In realistic liquid environments, typical currents in tunneling devices are of the order of picoamps. This corresponds to only six electrons per microsecond, and this number affects the integration time required to do current measurements in real experiments. This limits the speed of sequencing, though current fluctuations due to Brownian motion of the molecule average out during the required integration time. Moreover, data acquisition equipment introduces noise, and electronic filters create correlations in time-series data. We discuss how these effects must be included in the analysis of, e.g., the assignment of specific nucleobases to current signals. As the signals from different molecules overlap, unambiguous classification is impossible with a single measurement. We argue that the assignment of molecules to a signal is a standard pattern classification problem and calculation of the error rates is straightforward. The ideas presented here can be extended to other sequencing approaches of current interest.

How to connect time-lapse recorded trajectories of motile microorganisms with dynamical models in continuous time

We provide a tool for data-driven modeling of motility, data being time-lapse recorded trajectories. Several mathematical properties of a model to be found can be gleaned from appropriate model-independent experimental statistics, if one understands how such statistics are distorted by the finite sampling frequency of time-lapse recording, by experimental errors on recorded positions, and by conditional averaging. We give exact analytical expressions for these effects in the simplest possible model for persistent random motion, the Ornstein-Uhlenbeck process. Then we describe those aspects of these effects that are valid for any reasonable model for persistent random motion. Our findings are illustrated with experimental data and Monte Carlo simulations.

"Calibration-on-the-spot": How to calibrate an EMCCD camera from its images

In order to count photons with a camera, the camera must be calibrated. Photon counting is necessary, e.g., to determine the precision of localization-based super-resolution microscopy. Here we present a protocol that calibrates an EMCCD camera from information contained in isolated, diffraction-limited spots in any image taken by the camera, thus making dedicated calibration procedures redundant by enabling calibration post festum, from images filed without calibration information.

How to determine local stretching and tension in a flow-stretched DNA molecule

We determine the nonuniform stretching of and tension in a mega base pairs-long fragment of deoxyribonucleic acid (DNA) that is flow stretched in a nanofluidic chip. We use no markers, do not know the contour length of the DNA, and do not have the full DNA molecule inside our field of view. Instead, we analyze the transverse thermal motion of the DNA. Tension at the center of the DNA adds up to 16 pN, giving almost fully stretched DNA. This method was devised for optical mapping of DNA, specifically, DNA denaturation patterns. It may be useful also for other studies, e.g., DNA-protein interactions, specifically, their tension dependence. Generally, wherever long strands of DNA-e.g., native DNA extracted from human cells or bacteria-must be stretched with ease for inspection, this method applies.

Transition state theory demonstrated at the micron scale with out-of-equilibrium transport in a confined environment

Transition state theory (TST) provides a simple interpretation of many thermally activated processes. It applies successfully on timescales and length scales that differ several orders of magnitude: to chemical reactions, breaking of chemical bonds, unfolding of proteins and RNA structures and polymers crossing entropic barriers. Here we apply TST to out-of-equilibrium transport through confined environments: the thermally activated translocation of single DNA molecules over an entropic barrier helped by an external force field. Reaction pathways are effectively one dimensional and so long that they are observable in a microscope. Reaction rates are so slow that transitions are recorded on video. We find sharp transition states that are independent of the applied force, similar to chemical bond rupture, as well as transition states that change location on the reaction pathway with the strength of the applied force. The states of equilibrium and transition are separated by micrometres as compared with angstroms/nanometres for chemical bonds.

New technologies for DNA analysis--a review of the READNA Project

The REvolutionary Approaches and Devices for Nucleic Acid analysis (READNA) project received funding from the European Commission for 41/2 years. The objectives of the project revolved around technological developments in nucleic acid analysis. The project partners have discovered, created and developed a huge body of insights into nucleic acid analysis, ranging from improvements and implementation of current technologies to the most promising sequencing technologies that constitute a 3(rd) and 4(th) generation of sequencing methods with nanopores and in situ sequencing, respectively.

Configurational Statistics of Magnetic Bead Detection with Magnetoresistive Sensors

Magnetic biosensors detect magnetic beads that, mediated by a target, have bound to a functionalized area. This area is often larger than the area of the sensor. Both the sign and magnitude of the average magnetic field experienced by the sensor from a magnetic bead depends on the location of the bead relative to the sensor. Consequently, the signal from multiple beads also depends on their locations. Thus, a given coverage of the functionalized area with magnetic beads does not result in a given detector response, except on the average, over many realizations of the same coverage. We present a systematic theoretical analysis of how this location-dependence affects the sensor response. The analysis is done for beads magnetized by a homogeneous in-plane magnetic field. We determine the expected value and standard deviation of the sensor response for a given coverage, as well as the accuracy and precision with which the coverage can be determined from a single sensor measurement. We show that statistical fluctuations between samples may reduce the sensitivity and dynamic range of a sensor significantly when the functionalized area is larger than the sensor area. Hence, the statistics of sampling is essential to sensor design. For illustration, we analyze three important published cases for which statistical fluctuations are dominant, significant, and insignificant, respectively.

Optimized measurements of separations and angles between intra-molecular fluorescent markers

We demonstrate a novel, yet simple tool for the study of structure and function of biomolecules by extending two-colour co-localization microscopy to fluorescent molecules with fixed orientations and in intra-molecular proximity. From each colour-separated microscope image in a time-lapse movie and using only simple means, we simultaneously determine both the relative (x,y)-separation of the fluorophores and their individual orientations in space with accuracy and precision. The positions and orientations of two domains of the same molecule are thus time-resolved. Using short double-stranded DNA molecules internally labelled with two fixed fluorophores, we demonstrate the accuracy and precision of our method using the known structure of double-stranded DNA as a benchmark, resolve 10-base-pair differences in fluorophore separations, and determine the unique 3D orientation of each DNA molecule, thereby establishing short, double-labelled DNA molecules as probes of 3D orientation of anything to which one can attach them firmly.

Concentrating genomic length DNA in a microfabricated array

We demonstrate that a microfabricated bump array can concentrate genomic-length DNA molecules efficiently at continuous, high flow velocities, up to 40 μm/s, if the single-molecule DNA globule has a sufficiently large shear modulus. Increase in the shear modulus is accomplished by compacting the DNA molecules to minimal coil size using polyethylene glycol (PEG) derived depletion forces. We map out the sweet spot, where concentration occurs, as a function of PEG concentration and flow speed using a combination of theoretical analysis and experiment. Purification of DNA from enzymatic reactions for next-generation DNA-sequencing libraries will be an important application of this development.

Thermophoretic forces on DNA measured with a single-molecule spring balance

We stretch a single DNA molecule with thermophoretic forces and measure these forces with a spring balance: the DNA molecule itself. It is an entropic spring which we calibrate, using as a benchmark its Brownian motion in the nanochannel that contains and prestretches it. This direct measurement of the thermophoretic force in a static configuration finds forces up to 130 fN. This is eleven times stronger than the force experienced by the same molecule in the same thermal gradient in bulk, where the molecule shields itself. Our stronger forces stretch the middle of the molecule up to 80% of its contour length. We find the Soret coefficient per unit length of DNA at various ionic strengths. It agrees, with novel precision, with results obtained in bulk for DNA too short to shield itself and with the thermodynamic model of thermophoresis.

Optimal estimation of diffusion coefficients from single-particle trajectories

How does one optimally determine the diffusion coefficient of a diffusing particle from a single-time-lapse recorded trajectory of the particle? We answer this question with an explicit, unbiased, and practically optimal covariance-based estimator (CVE). This estimator is regression-free and is far superior to commonly used methods based on measured mean squared displacements. In experimentally relevant parameter ranges, it also outperforms the analytically intractable and computationally more demanding maximum likelihood estimator (MLE). For the case of diffusion on a flexible and fluctuating substrate, the CVE is biased by substrate motion. However, given some long time series and a substrate under some tension, an extended MLE can separate particle diffusion on the substrate from substrate motion in the laboratory frame. This provides benchmarks that allow removal of bias caused by substrate fluctuations in CVE. The resulting unbiased CVE is optimal also for short time series on a fluctuating substrate. We have applied our estimators to human 8-oxoguanine DNA glycolase proteins diffusing on flow-stretched DNA, a fluctuating substrate, and found that diffusion coefficients are severely overestimated if substrate fluctuations are not accounted for.

Assessment of automated analyses of cell migration on flat and nanostructured surfaces

Motility studies of cells often rely on computer software that analyzes time-lapse recorded movies and establishes cell trajectories fully automatically. This raises the question of reproducibility of results, since different programs could yield significantly different results of such automated analysis. The fact that the segmentation routines of such programs are often challenged by nanostructured surfaces makes the question more pertinent. Here we illustrate how it is possible to track cells on bright field microscopy images with image analysis routines implemented in an open-source cell tracking program, PACT (Program for Automated Cell Tracking). We compare the automated motility analysis of three cell tracking programs, PACT, Autozell, and TLA, using the same movies as input for all three programs. We find that different programs track overlapping, but different subsets of cells due to different segmentation methods. Unfortunately, population averages based on such different cell populations, differ significantly in some cases. Thus, results obtained with one software package are not necessarily reproducible by other software.

Cell motility, morphology, viability and proliferation in response to nanotopography on silicon black

Knowledge of cells' interactions with nanostructured materials is fundamental for bio-nanotechnology. We present results for how individual mouse fibroblasts from cell line NIH3T3 respond to highly spiked surfaces of silicon black that were fabricated by maskless reactive ion etching (RIE). We did standard measurements of cell viability, proliferation, and morphology on various surfaces. We also analyzed the motility of cells on the same surfaces, as recorded in time lapse movies of sparsely populated cell cultures. We find that motility and morphology vary strongly with nano-patterns, while viability and proliferation show little dependence on substrate type. We conclude that motility analysis can show a wide range of cell responses e.g. over a factor of two in cell speed to different nano-topographies, where standard assays, such as viability or proliferation, in the tested cases show much less variation of the order 10-20%.

Why is the microtubule lattice helical?

Microtubules polymerize from identical tubulin heterodimers, which form a helical lattice pattern that is the microtubule. This pattern always has left-handed chirality, but it is not known why. But as tubulin, similar to other proteins, evolved for a purpose, the question of the title of this artcile appears to be meaningful. In a computer simulation that explores the 'counterfactual biology' of microtubules without helicity, we demonstrate that these have the same mechanical properties as Nature's microtubules with helicity. Thus only a dynamical reason for helicity is left as potential explanation. We find that helicity solves 'the problem of the blind mason', i.e. how to correctly build a structure, guided only by the shape of the bricks. This answer in turn raises some new questions for researchers to address.

Intracellular signaling by diffusion: can waves of hydrogen peroxide transmit intracellular information in plant cells?

Amplitude- and frequency-modulated waves of Ca(2+) ions transmit information inside cells. Reactive Oxygen Species (ROS), specifically hydrogen peroxide, have been proposed to have a similar role in plant cells. We consider the feasibility of such an intracellular communication system in view of the physical and biochemical conditions in plant cells. As model system, we use a H(2)O(2) signal originating at the plasma membrane (PM) and spreading through the cytosol. We consider two maximally simple types of signals, isolated pulses and harmonic oscillations. First we consider the basic limits on such signals as regards signal origin, frequency, amplitude, and distance. Then we establish the impact of ROS-removing enzymes on the ability of H(2)O(2) to transmit signals. Finally, we consider to what extent cytoplasmic streaming distorts signals. This modeling allows us to predict the conditions under which diffusion-mediated signaling is possible. We show that purely diffusive transmission of intracellular information by H(2)O(2) over a distance of 1 μm (typical distance between organelles, which may function as relay stations) is possible at frequencies well above 1 Hz, which is the highest frequency observed experimentally. This allows both frequency and amplitude modulation of the signal. Signaling over a distance of 10 μm (typical distance between the PM and the nucleus) may be possible, but requires high signal amplitudes or, equivalently, a very low detection threshold. Furthermore, at this longer distance a high rate of enzymatic degradation is required to make signaling at frequencies above 0.1 Hz possible. In either case, cytoplasmic streaming does not seriously disturb signals. We conclude that although purely diffusion-mediated signaling without relaying stations is theoretically possible, it is unlikely to work in practice, since it requires a much faster enzymatic degradation and a much lower cellular background concentration of H(2)O(2) than observed experimentally.

'Dicty dynamics': Dictyostelium motility as persistent random motion

We model the motility of Dictyostelium cells in a systematic data-driven manner. We deduce a minimal dynamical model that reproduces the statistical features of experimental trajectories. These are trajectories of the centroid of the cell perimeter, which is more sensitive to pseudopod activity than the usual tracking by centroid or nucleus. Our data account for cell individuality and dictate a model that extends the cell-type specific models recently derived for mammalian cells. Two generalized Langevin equations model stochastic periodic pseudopod motion parallel and orthogonal to the amoeba's direction of motion. This motion propels the amoeba with a random periodic left-right waddle in a direction that has a long persistence time. The model fully accounts for the statistics of the experimental trajectories, including velocity power spectra and auto-correlations, non-Gaussian velocity distributions, and multiplicative noise. Thus, we find neither need nor place in our data for an interpretation in terms of anomalous diffusion. The model faithfully captures cell individuality as different parameter values in the model, and serves as a basis for integrating the local mechanics of cell motion with our observed long-term behavior.

Pressure-driven DNA in nanogroove arrays: complex dynamics leads to length- and topology-dependent separation

The motion of linear and circular DNA molecules is studied under pressure driven buffer flow in a 50 nm slit channel with arrays of transverse 150 nm deep nanogrooves. Transport occurs through two states of propagation unique to this nanogroove geometry, a slow, stepwise groove-to-groove translation called the "sidewinder" and a fast, continuous tumbling across the grooves called the "tumbleweed". Dynamical transitions between the two states are observed at fixed buffer velocity. Molecules exhibit size- and topology-dependent velocities.

Harmonic oscillator in heat bath: exact simulation of time-lapse-recorded data and exact analytical benchmark statistics

The stochastic dynamics of the damped harmonic oscillator in a heat bath is simulated with an algorithm that is exact for time steps of arbitrary size. Exact analytical results are given for correlation functions and power spectra in the form they acquire when computed from experimental time-lapse recordings. Three applications are discussed: (i) The effects of finite sampling rate and time, described exactly here, are similar for other stochastic dynamical systems--e.g., motile microorganisms and their time-lapse-recorded trajectories. (ii) The same statistics is satisfied by any experimental system to the extent that it is interpreted as a damped harmonic oscillator at finite temperature-such as an AFM cantilever. (iii) Three other models of fundamental interest are limiting cases of the damped harmonic oscillator at finite temperature; it consequently bridges their differences and describes the effects of finite sampling rate and sampling time for these models as well.

Integrative analysis correlates donor transcripts to recipient autoantibodies in primary graft dysfunction after lung transplantation

Up to one in four lung-transplanted patients develop pulmonary infiltrates and impaired oxygenation within the first days after lung transplantation. Known as primary graft dysfunction (PGD), this condition increases mortality significantly. Complex interactions between donor lung and recipient immune system are the suspected cause. We took an integrative, systems-level approach by first exploring whether the recipient's immune response to PGD includes the development of long-lasting autoreactivity. We next explored whether proteins displaying such differential autoreactivity also display differential gene expression in donor lungs that later develop PGD compared with those that did not. We evaluated 39 patients from whom autoantibody profiles were already available for PGD based on chest radiographs and oxygenation data. An additional nine patients were evaluated for PGD based on their medical records and set aside for validation. From two recent donor lung gene expression studies, we reanalysed and paired gene profiles with autoantibody profiles. Primary graft dysfunction can be distinguished by a profile of differentially reactive autoantibodies binding to 17 proteins. Functional analysis showed that 12 of these proteins are part of a protein-protein interaction network (P=3 x 10⁻⁶) involved in proliferative processes. A nearest centroid classifier assigned correct PGD grades to eight out of the nine patients in the validation cohort (P=0·048). We observed significant positive correlation (r=0·63, P=0·011) between differences in IgM reactivity and differences in gene expression levels. This connection between donor lung gene expression and long-lasting recipient IgM autoantibodies towards a specific set of proteins suggests a mechanism for the development of autoimmunity in PGD.

Single-molecule denaturation mapping of DNA in nanofluidic channels

Here we explore the potential power of denaturation mapping as a single-molecule technique. By partially denaturing YOYO-1-labeled DNA in nanofluidic channels with a combination of formamide and local heating, we obtain a sequence-dependent "barcode" corresponding to a series of local dips and peaks in the intensity trace along the extended molecule. We demonstrate that this structure arises from the physics of local denaturation: statistical mechanical calculations of sequence-dependent melting probability can predict the barcode to be observed experimentally for a given sequence. Consequently, the technique is sensitive to sequence variation without requiring enzymatic labeling or a restriction step. This technique may serve as the basis for a new mapping technology ideally suited for investigating the long-range structure of entire genomes extracted from single cells.

Power spectrum analysis with least-squares fitting: amplitude bias and its elimination, with application to optical tweezers and atomic force microscope cantilevers

Optical tweezers and atomic force microscope (AFM) cantilevers are often calibrated by fitting their experimental power spectra of Brownian motion. We demonstrate here that if this is done with typical weighted least-squares methods, the result is a bias of relative size between -2/n and +1/n on the value of the fitted diffusion coefficient. Here, n is the number of power spectra averaged over, so typical calibrations contain 10%-20% bias. Both the sign and the size of the bias depend on the weighting scheme applied. Hence, so do length-scale calibrations based on the diffusion coefficient. The fitted value for the characteristic frequency is not affected by this bias. For the AFM then, force measurements are not affected provided an independent length-scale calibration is available. For optical tweezers there is no such luck, since the spring constant is found as the ratio of the characteristic frequency and the diffusion coefficient. We give analytical results for the weight-dependent bias for the wide class of systems whose dynamics is described by a linear (integro)differential equation with additive noise, white or colored. Examples are optical tweezers with hydrodynamic self-interaction and aliasing, calibration of Ornstein-Uhlenbeck models in finance, models for cell migration in biology, etc. Because the bias takes the form of a simple multiplicative factor on the fitted amplitude (e.g. the diffusion coefficient), it is straightforward to remove and the user will need minimal modifications to his or her favorite least-squares fitting programs. Results are demonstrated and illustrated using synthetic data, so we can compare fits with known true values. We also fit some commonly occurring power spectra once-and-for-all in the sense that we give their parameter values and associated error bars as explicit functions of experimental power-spectral values.

Optimized localization analysis for single-molecule tracking and super-resolution microscopy

We optimally localize isolated fluorescent beads and molecules imaged as diffraction-limited spots, determine the orientation of molecules, and present reliable formulae for the precisions of various localization methods. For beads, theory and experimental data both show that unweighted least-squares fitting of a Gaussian squanders one third of the available information, a popular formula for its precision exaggerates beyond Fisher's information limit, and weighted least-squares may do worse, while maximum likelihood fitting is practically optimal.

Chronic rejection of a lung transplant is characterized by a profile of specific autoantibodies

SUMMARY

Obliterative bronchiolitis (OB) continues to be the major limitation to long-term survival after lung transplantation. The specific aetiology and pathogenesis of OB are not well understood. To explore the role of autoreactivity in OB, we spotted 751 different self molecules onto glass slides, and used these antigen microarrays to profile 48 human serum samples for immunoglobulin G (IgG) and IgM autoantibodies; 27 patients showed no or mild bronchiolitis obliterans syndrome (BOS; a clinical correlate of OB) and 15 patients showed medium to severe BOS. We now report that these BOS grades could be differentiated by a profile of autoantibodies binding to 28 proteins or their peptides. The informative autoantibody profile included down-regulation as well as up-regulation of both IgM and IgG specific reactivities. This profile was evaluated for robustness using a panel of six independent test patients. Analysis of the functions of the 28 informative self antigens showed that eight of them are connected in an interaction network involved in apoptosis and protein metabolism. Thus, a profile of autoantibodies may reflect pathological processes in the lung allograft, suggesting a role for autoimmunity in chronic rejection leading to OB.

Combining single-molecule optical trapping and small-angle x-ray scattering measurements to compute the persistence length of a protein ER/K alpha-helix

A relatively unknown protein structure motif forms stable isolated single alpha-helices, termed ER/K alpha-helices, in a wide variety of proteins and has been shown to be essential for the function of some molecular motors. The flexibility of the ER/K alpha-helix determines whether it behaves as a force transducer, rigid spacer, or flexible linker in proteins. In this study, we quantify this flexibility in terms of persistence length, namely the length scale over which it is rigid. We use single-molecule optical trapping and small-angle x-ray scattering, combined with Monte Carlo simulations to demonstrate that the Kelch ER/K alpha-helix behaves as a wormlike chain with a persistence length of approximately 15 nm or approximately 28 turns of alpha-helix. The ER/K alpha-helix length in proteins varies from 3 to 60 nm, with a median length of approximately 5 nm. Knowledge of its persistence length enables us to define its function as a rigid spacer in a translation initiation factor, as a force transducer in the mechanoenzyme myosin VI, and as a flexible spacer in the Kelch-motif-containing protein.

Nanoconfinement-enhanced conformational response of single DNA molecules to changes in ionic environment

We show that the ionic environment plays a critical role in determining the configurational properties of DNA confined in silica nanochannels. The extension of DNA in the nanochannels increases as the ionic strength is reduced, almost tripling over two decades in ionic strength for channels around 100 x 100 nm in dimension. Surprisingly, we find that the variation of the persistence length alone with ionic strength is not enough to explain our results. The effect is due mainly to increasing self-avoidance created by the reduced screening of electrostatic interactions at low ionic strength. To quantify the increase in self-avoidance, we introduce a new parameter into the de Gennes theory: an effective DNA width that gives the increase in the excluded volume due to electrostatic repulsion.

Strong physical constraints on sequence-specific target location by proteins on DNA molecules

Sequence-specific binding to DNA in the presence of competing non-sequence-specific ligands is a problem faced by proteins in all organisms. It is akin to the problem of parking a truck at a loading bay by the side of a road in the presence of cars parked at random along the road. Cars even partially covering the loading bay prevent correct parking of the truck. Similarly on DNA, non-specific ligands interfere with the binding and function of sequence-specific proteins. We derive a formula for the probability that the loading bay is free from parked cars. The probability depends on the size of the loading bay and allows an estimation of the size of the footprint on the DNA of the sequence-specific protein by assaying protein binding or function in the presence of increasing concentrations of non-specific ligand. Assaying for function gives an 'activity footprint'; the minimum length of DNA required for function rather than the more commonly measured physical footprint. Assaying the complex type I restriction enzyme, EcoKI, gives an activity footprint of approximately 66 bp for ATP hydrolysis and 300 bp for the DNA cleavage function which is intimately linked with translocation of DNA by EcoKI. Furthermore, considering the coverage of chromosomal DNA by proteins in vivo, our theory shows that the search for a specific DNA sequence is very difficult; most sites are obscured by parked cars. This effectively rules out any significant role in target location for mechanisms invoking one-dimensional, linear diffusion along DNA.

The free NADH concentration is kept constant in plant mitochondria under different metabolic conditions

The reduced coenzyme NADH plays a central role in mitochondrial respiratory metabolism. However, reports on the amount of free NADH in mitochondria are sparse and contradictory. We first determined the emission spectrum of NADH bound to proteins using isothermal titration calorimetry combined with fluorescence spectroscopy. The NADH content of actively respiring mitochondria (from potato tubers [Solanum tuberosum cv Bintje]) in different metabolic states was then measured by spectral decomposition analysis of fluorescence emission spectra. Most of the mitochondrial NADH is bound to proteins, and the amount is low in state 3 (substrate + ADP present) and high in state 2 (only substrate present) and state 4 (substrate + ATP). By contrast, the amount of free NADH is low but relatively constant, even increasing a little in state 3. Using modeling, we show that these results can be explained by a 2.5- to 3-fold weaker average binding of NADH to mitochondrial protein in state 3 compared with state 4. This indicates that there is a specific mechanism for free NADH homeostasis and that the concentration of free NADH in the mitochondrial matrix per se does not play a regulatory role in mitochondrial metabolism. These findings have far-reaching consequences for the interpretation of cellular metabolism.

A non-Gaussian distribution quantifies distances measured with fluorescence localization techniques

When single-molecule fluorescence localization techniques are pushed to their lower limits in attempts to measure ever-shorter distances, measurement errors become important to understand. Here we describe the non-Gaussian distribution of measured distances that is the key to proper interpretation of distance measurements. We test it on single-molecule high-resolution colocalization data for a known distance, 10 nm, and find that it gives the correct result, whereas interpretation of the same data with a Gaussian distribution gives a result that is systematically too large.

Cell motility as persistent random motion: theories from experiments

Experimental time series for trajectories of motile cells may contain so much information that a systematic analysis will yield cell-type-specific motility models. Here we demonstrate how, using human keratinocytes and fibroblasts as examples. The two resulting models reflect the cells' different roles in the organism, it seems, and show that a cell has a memory of past velocities. They also suggest how to distinguish quantitatively between various surfaces' compatibility with the two cell types.

Microtubule's conformational cap

The molecular mechanisms that allow elongation of the unstable microtubule lattice remain unclear. It is usually thought that the GDP-liganded tubulin lattice is capped by a small layer of GTP- or GDP-Pi-liganded molecules, the so called "GTP-cap". Here, we point-out that the elastic properties of the microtubule lattice cause a difference in stability between the elongating tubulin sheet and the completed microtubule wall. The implications of our observations for microtubule structure and dynamics are discussed.

Modelling NADH turnover in plant mitochondria

NADH is central to the functioning of mitochondrial respiration. It is produced by enzymes in, or associated with, the tricarboxylic acid cycle in the matrix, and it is oxidized by two respiratory chain enzymes in the inner membrane, the rotenone-sensitive complex I and the rotenone-insensitive internal NADH dehydrogenase (ND(in)). A simplified kinetic model for NADH turnover in the matrix of plant mitochondria is presented. Only the two main NADH-producing enzymes, NAD-malate dehydrogenase [EC 1.1.1.37] (MDH) and NAD-malic enzyme [EC 1.1.1.39] (ME), are considered. This model reproduces the complex behaviour of malate oxidation by isolated mitochondria in response to additions of ADP (state 3/state 4), NAD(+) and/or rotenone, as well as to changes in pH. It is found that MDH always operates at or close to equilibrium. Changes in the activity of complex I, ND(in), or ME are predicted to cause clear changes in the pattern of malate oxidation. In general, the model predicts high sensitivity to changes in the ME activity. In contrast, MDH activity can be reduced 100-fold without detectable changes in malate oxidation. It is demonstrated that it is not the high activity, but the equilibrium properties of MDH that are important for the redox-buffering function of MDH in the mitochondrial matrix. Binding of NAD(+) and NADH in the matrix reduces the concentrations of free NAD(+) and NADH, depending on the concentration of binding sites and the binding strength. On the basis of the modelling results it is estimated that a significant proportion of the mitochondrial NAD is bound.

Structural microtubule cap: stability, catastrophe, rescue, and third state

Microtubules polymerize from GTP-liganded tubulin dimers, but are essentially made of GDP-liganded tubulin. We investigate the tug-of-war resulting from the fact that GDP-liganded tubulin favors a curved configuration, but is forced to remain in a straight one when part of a microtubule. We point out that near the end of a microtubule, the proximity of the end shifts the balance in this tug-of-war, with some protofilament bending as result. This somewhat relaxes the microtubule lattice near its end, resulting in a structural cap. This structural cap thus is a simple mechanical consequence of two well-established facts: protofilaments made of GDP-liganded tubulin have intrinsic curvature, and microtubules are elastic, made from material that can yield to forces, in casu its own intrinsic forces. We explore possible properties of this structural cap, and demonstrate 1) how it allows both polymerization from GTP-liganded tubulin and rapid depolymerization in its absence; 2) how rescue can occur; 3) how a third, meta-stable intermediate state is possible and can explain some experimental results; and 4) how the tapered tips observed at polymerizing microtubule ends are stabilized during growth, though unable to accommodate a lateral cap. This scenario thus supports the widely accepted GTP-cap model by suggesting a stabilizing mechanism that explains the many aspects of dynamic instability.