Oseen approximation 2026-10-14
The Oseen approximation retains advection of a small disturbance by a uniform incident velocity , but neglects advection of the disturbance by itself. Its steady equations are
A Stokes disturbance of a sphere decays as , so the ratio of background advection to viscous diffusion grows as . Thus the approximation resolves the outer region , where the Stokes flow expansion is nonuniform, while self-convection remains smaller by .
Let , , and , so the single-bead mobility is . Assume an incompressible solvent, constant , negligible bead inertia, and a symmetric positive-definite . The last requirement makes the resistance invertible and Brownian diffusion physically admissible. At a given position and time, depends on the moment but not on the individual configuration .
Use the source's gradient convention , and write for the component-first gradient. Then . The spring forces on the two beads are and . Their overdamped probability currents are
Here the Brownian term written as is a force weighted by probability density. Its corresponding conditional osmotic force is ; treating the given term as an unweighted deterministic force per bead would not produce the stated Fokker–Planck equation.
With midpoint , the derivatives are and . Subtracting the bead currents gives the configuration current
The midpoint current is after neglecting the specified spatial Brownian diffusion. Configuration probability conservation therefore gives
Since has components , this agrees with the source's . The factor two in both spring relaxation and configuration diffusivity comes from relative motion of both beads.
The second moment is the unnormalized conformation tensor. Now multiply by and integrate over configuration space. Assume sufficiently decaying tails so every boundary term at infinite separation vanishes. The affine current gives ; spring drift gives ; two integrations by parts of diffusion give . Hence
Because , is a polynomial in and commutes with it. Thus
This is the Hookean dumbbell conformation equation with anisotropic mobility, an exact second-moment closure within the stated model; no Gaussian approximation to was needed.
The upper-convected derivative of the identity is . Applying it to the polymer stress therefore gives
Consequently the coefficients and the closed stress equation are
where in the printed convention. Here is a polymer relaxation time, unrelated to the inverse-length parameter of question c; is the polymer viscosity scale. The coefficient has inverse-stress units. In the dimensional mobility parameter in the Giesekus model, the dimensionless factor is , not the dimensional printed in this problem.
The upper-convected rate removes the affine stretching and rotation of material connectors. If , then a tensor transported purely as has . Under a time-dependent rigid rotation , the rate transforms as : this is material frame indifference, not an ordinary componentwise time derivative. The resulting equation is the Giesekus model in these dimensional conventions.
The isotropic subtraction is fixed by zero polymer extra stress at thermal equilibrium. With no flow and positive mobility, the stationary distribution is
so . The constitutive equation correspondingly has a zero-stress quiescent state and a Newtonian low-rate limit . Without subtracting exactly , the spring contribution would retain an equilibrium isotropic stress and the equation would acquire corresponding reference terms. Because mobility itself depends on , that reference cannot simply be shifted while retaining the same model. The separate fluid pressure enforcing incompressibility is not included in this polymer extra stress.
For steady uniaxial extension of the Giesekus model, the component-first velocity gradient is diagonal, . For a steady spatially uniform diagonal stress, the material derivative vanishes and the convected terms are and . The off-diagonal equations vanish identically, and symmetry permits equal transverse components. The remaining equations are
Set , and , with . For , the branches continuously connected to zero stress at zero rate are
A sufficient nondegenerate range for physical solutions at every nonnegative rate is . Indeed with , so positive conformation implies positive mobility. The longitudinal quadratic equation is nonpositive at and has a nonnegative root of a polynomial. The transverse quadratic equation has values at and at , so its larger root of a polynomial lies in . Thus the displayed branches satisfy , giving . The other algebraic roots do not provide a physical branch connected to the equilibrium distribution. For constant positive diagonal , the associated stationary configuration distribution is a Gaussian with covariance .
For positive on the physical branch, the stretching stress remains finite through . At that rate,
At large extension rate,
Thus the stress continues to grow, but the uniaxial extensional viscosity approaches the finite plateau rather than developing a finite-rate pole. The transverse stress approaches as its conformation tensor variance becomes small. At small rate the two stresses are and , so .
When , the physical steady branch instead is
The stretching stress diverges as from below. Above that threshold its algebraic continuation has and cannot be a conformation covariance. There is no physical finite steady stretching state. The transient equation
gives linear growth at and exponential growth at larger rates, the Hookean coil-stretch failure of the upper-convected Maxwell limit.
Figure 1.
Anisotropic mobility rounds the Hookean extensional-stress singularity
. Dimensionless longitudinal and transverse stresses are plotted against the extension-rate parameter W. Positive dimensionless mobility factors 0.1 and 0.5 give regular equilibrium-connected branches. The isotropic-mobility longitudinal branch is shown only below its physical steady-state threshold W equal to one half. Details.
Parameter qualification. Positivity of alone, without a restriction on , does not ensure real physical steady solutions for the entire tensor. For example , makes the transverse equation , which has no real root. Its diffusion/resistance interpretation is already inadmissible for some conformations. The physically positive-mobility regime stated above avoids that counterexample; the source's positive- case must not be taken as a guarantee for arbitrary parameters.
The membrane supplies G proteins with an address, an orienting surface, a reaction environment and a route for spatial redistribution. Its role therefore goes beyond retaining a protein near a receptor. Lipid composition, covalent anchors, membrane organization and trafficking can alter which partners meet, which interfaces are exposed and how long a signal persists.
For a heterotrimeric G protein, an activated G protein-coupled receptor promotes release of GDP from the subunit and binding of GTP. The nucleotide switch changes interaction with the complex and downstream effectors; GTP hydrolysis favors termination and reassociation. This is not phosphorylation of GDP by the receptor. Both the subunit and the complex can participate in downstream signaling, and their activation need not entail complete physical separation in every receptor–subunit combination. The relevant question is whether the effector-binding surfaces and local interactions change.
Covalent protein lipidation helps place this cycle at a membrane. Many subunits undergo S-palmitoylation, some also undergo N-myristoylation, and the subunit commonly undergoes prenylation near its C-terminus. These modifications attach chemically different hydrophobic groups; their distributions depend on the subunit, so no single anchor recipe describes every G protein. Additional basic residues, contacts among subunits and membrane-facing surfaces stabilize association. The membrane anchor on helps position the whole complex. Because one attachment can be individually weak or nonspecific, combined signals can determine both residence time and preference for a particular membrane.
small GTPases illustrate the same logic in a distinct signaling architecture. Ras isoforms have a prenylation anchor and additional targeting information. Palmitoylation contributes to targeting of some isoforms, while the polybasic region of KRAS4B promotes interaction with negatively charged inner-leaflet lipids. Reversible acylation and trafficking can therefore change localization without altering the nucleotide-binding core. Mutating an anchoring site should be interpreted with expression, processing and nucleotide-loading controls: an inactive cytosolic mutant does not by itself distinguish loss of a necessary location from a direct catalytic defect.
Membrane tethering can enhance productive encounters by restricting motion to a surface and increasing local availability near a membrane-associated partner. It also sets orientation and limits the orientations accessible to the protein. More surface density does not guarantee faster signaling: trapping in an inappropriate region, unfavorable orientation or crowding can reduce access. A clean experimental comparison is to retarget an anchor-deficient protein with an independent membrane tether, confirming comparable surface density. Rescue would support the importance of localization; failure could indicate the original anchor also specifies orientation, lipid preference or trafficking, or that the replacement tether introduces an artifact.
The membrane can regulate even a correctly localized G protein through contacts with its folded domain. For KRAS4B, NMR spectroscopy on lipid bilayer nanodiscs demonstrated membrane-associated orientations in which an effector-binding surface can be occluded. Some disease-associated mutations changed this interaction and favored effector accessibility.[ref] Thus a protein can remain attached while its signaling competence changes. The inference is specific to the tested constructs and lipid environments; it does not prove that every Ras isoform or every active state has one universal membrane orientation.
Lipids can also act at the receptor–G protein interface. Phosphatidylinositol 4,5-bisphosphate is a negatively charged inner-leaflet signaling lipid, not simply an inert bulk solvent. Native mass spectrometry and complementary experiments showed that it can stabilize active G protein-coupled receptor complexes and influence G protein coupling selectivity.[ref] Specific lipid binding, general surface electrostatics and changes in membrane mechanics are competing or cooperating explanations in any new system. Vary the relevant lipid while holding total charge and other physical properties as constant as possible; compare mutations in the candidate binding site with control surface mutations. A bound lipid detected after purification is useful evidence but needs functional corroboration in a lipid bilayer.
Bulk composition further affects hydrophobic mismatch, membrane thickness, local order and the conformational cost of receptor activation. Cholesterol and sphingomyelin can alter the organization of a lipid bilayer; acyl-chain diversity can affect diffusion and anchor partitioning. Distinct membrane nanodomains may enrich some receptors, G proteins and effectors while separating others. These are dynamic, context-dependent organizations. Detergent-resistant extraction does not demonstrate a pre-existing domain of a particular size in a living biological cell, and diffraction-limited colocalization cannot establish molecular interaction. Use live-cell tracking, proximity readouts and controlled membrane reconstitution, with caution that cholesterol depletion also changes many unrelated membrane processes.
The membrane is also a chemical participant downstream. Activation of phospholipase C hydrolyzes phosphatidylinositol 4,5-bisphosphate, producing membrane-retained diacylglycerol and soluble inositol trisphosphate. This both generates distinct signals and changes the abundance of a receptor-associated lipid. Some ion channels are regulated by and by phosphatidylinositol 4,5-bisphosphate, linking local protein interactions to local lipid availability. Consequently, removing this lipid could impair receptor coupling, an effector, or both. Separate receptor-mediated nucleotide exchange from effector activity in a reconstituted assay before assigning the step affected.
Localization changes in time. Receptor activation can redistribute particular combinations from the plasma membrane to internal membranes, creating opportunities for signaling at the Golgi apparatus and other sites. Experiments on receptor-induced movement linked this redistribution to Golgi apparatus structure and secretion.[ref] Subunit-dependent anchor interactions help determine the kinetics. Receptor trafficking can additionally move signaling platforms to endosomes: conformation-selective antibody probes have detected active beta-2 adrenergic receptors and their G proteins there, with internalized receptors contributing to the cellular response.[ref] A signal measured from the whole biological cell may therefore persist after its surface component declines. These possibilities require compartment-resolved measurements rather than an assumption that all signaling occurs at the surface.
To separate these roles, reconstitute a receptor and a defined heterotrimeric G protein in lipid bilayer nanodiscs or lipid vesicles of controlled composition, verifying receptor orientation and active fraction. Measure nucleotide exchange, GTP hydrolysis and effector response separately, then compare live-cell localization and signaling biosensors. Quantify expression, anchor processing, membrane integrity and sensor specificity; include rescue of a depleted lipid and an independent method of perturbing the same variable. The same G protein nucleotide state can generate different outputs because membrane attachment, orientation, lipid contacts and compartmental residence determine whether a productive signaling complex can form.
Lipid diversity allows a membrane to remain a selective, mechanically coherent barrier while retaining the molecular mobility required for function. Stability and fluidity are different properties: a lipid bilayer can be well ordered yet permit rapid lateral diffusion, or be mobile yet leaky and vulnerable to rupture. The effects depend on composition, temperature, leaflet organization and interactions among components rather than on one universal ranking of “rigid” and “fluid” lipids.
The basic stability comes from amphiphilic organization. Hydrophobic chains avoid exposure to water, while polar headgroups remain hydrated. A pore exposing chains has an energetic cost, encouraging local defects to close, but stress, chain disorder or incompatible molecular shapes can make defects easier to form. Useful measurements therefore include membrane permeability, leakage and mechanical failure as well as rotational mobility, lateral diffusion and acyl-chain order. Calling all of these “fluidity” conceals distinctions relevant to cell survival.
Acyl-chain length changes the contact area between neighboring chains. Other things equal, longer chains generally strengthen attractive contacts, increase bilayer thickness and raise the tendency to form an ordered phase. Shorter chains tend to lower the chain-melting transition and can increase permeability. A double bond between carbon atoms in the cis configuration restricts the chain geometry and disrupts close packing, generally lowering the chain-melting temperature relative to a comparable saturated chain. A double bond in the trans configuration has a different geometry and commonly permits closer packing than the corresponding cis configuration. Branching also changes packing. These comparisons hold the headgroup and other variables fixed; they are not predictions that all unsaturated lipids behave identically in a complex mixture.
An immediate consequence is that composition has to be interpreted relative to the growth temperature. A membrane made mostly of long saturated chains can become too ordered on cooling, slowing membrane protein conformational changes and molecular redistribution. Greater unsaturation can preserve mobility. On heating, the converse adjustments can limit excessive disorder and leakage. The original homeoviscous adaptation experiments found temperature-dependent changes in the saturated and long-chain fatty acid content of Escherichia coli, using spin-label electron paramagnetic resonance to compare membrane-lipid viscosity at the cells' growth temperatures.[ref] This establishes an adaptive relationship rather than a claim that every organism maintains a single invariant value by the same biosynthetic route.
Headgroup diversity changes both local interactions and preferred geometry. Phosphatidylcholine is approximately cylindrical in many common bilayer compositions, whereas phosphatidylethanolamine has a smaller headgroup relative to its chains and can favor negative monolayer spontaneous curvature. A single-chain lysophospholipid often favors curvature of the opposite sign. Mixing these shapes can support bending and membrane-fusion intermediates, while too much of a nonbilayer-forming component can destabilize a planar lipid bilayer. The amphiphile packing parameter is a useful qualitative guide, but hydrated headgroup size and chain volume depend on conditions; a molecular cartoon does not supply a quantitative phase diagram.
Headgroup charge adds another control. phosphatidylserine and phosphatidylinositol 4,5-bisphosphate contribute negative inner-leaflet charge, recruiting basic membrane-binding regions of proteins and changing local electrostatics. Ionic strength, pH and divalent ions can screen, alter or bridge those charges. A lipid-induced change in domain organization could therefore arise from ion-mediated interactions rather than acyl-chain order alone. Sphingomyelin can pack tightly through its chain structure and headgroup/backbone interactions, influencing thickness and local order. A comparison of headgroups must specify the attached chains, because changing a headgroup while simultaneously changing chain saturation confounds the interpretation.
Cholesterol is particularly effective because its small polar head and rigid ring system fit between phospholipid chains. In a relatively disordered bilayer, it often increases chain order, condenses packing and reduces permeability to small solutes. In a tendency toward a tightly packed gel, it can disrupt regular packing and inhibit the cooperative freezing transition. Thus “cholesterol decreases fluidity” is incomplete: it can reduce some local motions while preventing an immobile gel and retaining lateral mobility. Its effects depend on concentration, temperature and the surrounding lipids; excessive amounts can exceed the mixture's capacity to accommodate it.
Mixtures of saturated lipids, unsaturated lipids and cholesterol can form coexisting liquid-ordered phase and liquid-disordered phase. The former combines greater acyl-chain order with lateral mobility; it is distinct from a solid gel. fluorescence microscopy of defined giant lipid vesicles directly demonstrated domains and mapped composition- and temperature-dependent coexistence.[ref] The result proves that lipid interactions can organize a membrane without a protein scaffold. It does not establish that living plasma membranes contain large equilibrium domains of the same sizes. Cellular crowding, cytoskeleton coupling, asymmetry, trafficking and active turnover can change or suppress those patterns.
Diversity also works across the membrane. The outer and inner membrane leaflets differ in headgroups and chain distributions; equal total composition does not imply equal surface properties. Lipidomic and packing measurements found greater unsaturation and looser packing in the cytoplasmic leaflet of the examined mammalian plasma membranes.[ref] Such asymmetry can couple to membrane protein shape and binding preferences. flippases and scramblases change transverse lipid distribution; because spontaneous translocation of many polar phospholipids is slow, this organization can be maintained and regulated. Leaflet imbalance creates curvature stress, but the result depends on lipid exchange and mechanical constraints rather than on headgroup charge alone.
Membrane composition also sets the energetic cost of deforming the surface around a membrane protein. Longer ordered chains may increase thickness and stiffness; mismatch with a transmembrane protein's hydrophobic span can favor local deformation, changed association or altered activity. Curvature-preferring lipids can lower the cost of particular conformational or membrane-shape transitions. Greater mechanical stiffness can resist bending but need not mean a larger rupture threshold. Oxidation of polyunsaturated chains can change molecular geometry and produce defects or reactive breakdown products, making chemical stability another tradeoff alongside mobility. Appropriate experiments should measure oxidation state and leakage rather than infer stability from saturation alone.
A good causal experiment uses compositionally defined lipid vesicles, independently varying chain saturation, chain length, headgroup and cholesterol. Keep vesicle size, temperature, ionic strength and osmotic conditions controlled. differential scanning calorimetry identifies cooperative transitions; spin-label electron paramagnetic resonance or calibrated fluorescence anisotropy probes local order and rotational freedom; fluorescence recovery after photobleaching measures lateral redistribution. Recovery also depends on geometry, immobile fractions and probe partitioning, so distinguish these from the diffusion coefficient. Test barrier stability by entrapped-dye leakage with dye-quenching and complete-release controls, and mechanical behavior with defined membrane-tension measurements. Repeat key comparisons with different probes because a probe can partition preferentially into one phase. In living biological cells, combine lipidomics with targeted biosynthetic perturbation, rescue and viability measurements to separate compositional effects from secondary stress responses.
Lipid diversity stabilizes useful combinations of properties: a sealed barrier, adaptable curvature, controlled local order and sufficient mobility. It also allows different compartments and leaflets to satisfy different physical demands; no single lipid composition maximizes all of those properties simultaneously.
Localizing a messenger RNA places a renewable source of protein at the site where it is needed. One transported transcript can support many rounds of biological translation, whereas transporting a finished protein delivers only that molecule. Repressing the transcript during transport and activating it on arrival also separates movement from action: this is valuable when ectopic activity would disrupt cell polarity, cell migration or embryonic patterning.
The advantage is especially clear in neurons. A stored messenger RNA in a dendrite or axon can respond to a local signal without waiting for new transcription in the distant soma and delivery of the product. Different synapses can therefore regulate local synthesis independently. Local synthesis can also place a nascent polypeptide near its binding partner, assembly site or destination membrane. Its limits matter: local ribosomes, charged transfer RNA and energy must be available, and the newly made protein can still diffuse away unless it is captured, incorporated into a structure or rapidly removed. Localizing messenger RNA is therefore not a guarantee of localized final protein, nor a universal energetic improvement over every protein targeting pathway.
Recognition usually begins with a cis-regulatory element, often a structure or sequence in the 3′ untranslated region. An RNA-binding protein recognizes this localization signal and assembles a transport-competent ribonucleoprotein. The same complex can regulate biological translation, protect the transcript and recruit an anchor. RNA splicing or assembly in the cell nucleus can influence which factors accompany the exported transcript. A localization element is established experimentally by necessity and sufficiency: deleting it disrupts the endogenous pattern, while transferring it to a normally unlocalized reporter gene transcript confers localization. Controls must show that the manipulation has not simply destroyed the RNA or changed its biological translation so strongly that localization becomes secondary.
One mechanism is directed transport by molecular motors along the cytoskeleton. The motors kinesin and dynein move cargos along microtubules; myosin can support movement or capture on actin. Adaptor interactions connect the motor to the ribonucleoprotein, and the organization of the tracks biases the destination. Individual particles can move in both directions, pause and reverse. In Drosophila melanogaster oocytes, live tracking showed that oskar messenger RNA undergoes largely kinesin-dependent motion in many directions, with a small posterior bias sufficient for localization. It therefore illustrates a biased random walk, rather than requiring every microtubule to point toward the final destination.[ref]
A distinct mechanism is diffusion followed by selective capture or anchoring. Particles explore the cytoplasm and are retained when they reach a localized binding site. Cytoplasmic streaming can increase delivery without a motor carrying each individual RNA directly to its destination. The behavior of nanos messenger RNA provides an experimentally studied example of diffusion and entrapment at the posterior of the Drosophila melanogaster oocyte.[ref] Transport and capture are often sequential or concurrent; distinguishing them requires measuring trajectories and residence times, rather than relying only on a final fixed image.
A third mechanism is spatially selective stability. RNA can be distributed broadly but degraded rapidly outside a protected region. The localized population then reflects longer residence through slower decay, even if delivery itself is unbiased. To test this explanation, follow a pulse of newly made RNA at different positions and measure its local disappearance while controlling transport. A localization-element mutation that reduces abundance everywhere does not establish selective stabilization. Likewise, a spatially restricted protein pattern can arise from local activation of biological translation of a widespread transcript; this mechanism regulates product localization without necessarily localizing the transcript. Measuring both species distinguishes these possibilities.
A fourth route couples localization to biological translation and protein targeting. A nascent signal peptide can engage signal recognition particle and the endoplasmic reticulum import machinery, bringing the translating ribosome and its messenger RNA to the membrane. Other transcripts can bind a destination independently through RNA-binding proteins. Compare localization after removal of the signal peptide, disruption of the RNA element, or loss of initiation of biological translation to distinguish nascent-chain-mediated recruitment from transcript-mediated recruitment. Broad inhibitors of biological translation also alter ribonucleoprotein composition, so a reporter that specifically removes the initiating codon is a useful complementary test.
Localization must be coordinated with activation. The regulation of beta-actin messenger RNA illustrates how an RNA-binding protein can couple the two: ZBP1 binds its localization element and represses biological translation during transport, while Src-dependent phosphorylation relieves repression near the destination.[ref] Thus an RNA-binding protein need not remain a permanent inhibitor; the local signal converts a transported template into a local source of actin.
Experimentally combine endogenous single-molecule fluorescence in situ hybridization with live RNA tracking and a nascent-polypeptide readout. Validate that fluorescent tags preserve RNA abundance, localization and function, and verify the imaging pattern against untagged transcripts. Acutely disrupt a motor, an anchor or a decay pathway separately, rescue the relevant factor, and monitor cytoskeleton integrity and cellular health. Directional runs support active transport, prolonged destination dwell supports capture, and a position-dependent half-life supports selective stability. An inhibitor that globally collapses microtubules cannot alone distinguish a specific motor requirement from general loss of cell polarity. The persuasive explanation connects delivery, retention and local synthesis to the measured biological outcome.
For a fixed perfectly absorbing spherical target, steady diffusion outside radius obeys , , and . Spherical symmetry gives , hence . The inward diffusion flux integrated over the target is . Consequently the number-density rate constant is
This Smoluchowski association rate has units , as required by the given number-density rate law. If both partners moved independently, their relative diffusion coefficient would be the sum; the chromosomal target is treated as fixed here.
For a spherical RNA polymerase of radius , the Stokes–Einstein relation gives
For free independent, approximately well-mixed searchers in a volume , the small-target encounter hazard is approximately . Thus
The numerical coefficient explicitly uses a typical illustrative cell volume and an illustrative effective capture radius, not dimensions supplied by the source. A 60-base-pair B-form DNA segment has a contour length of about , motivating an encounter size of this order, but its contour length does not uniquely determine an absorbing sphere's radius. Taking instead doubles the estimate. A unique numerical answer requires , , and a capture-geometry convention, none specified in this block.
The estimate also treats all 6000 molecules as freely available and every encounter as productive. Crowding, nonspecific binding, orientation constraints, imperfect reaction, and finite-domain initial transients can alter it. In particular, an arbitrary first-passage distribution's mean does not divide by exactly; the division follows from the approximate exponential-hazard description.[ref]
Suppose initially . Positive ions initially move by diffusion from side 1 to side 2. Since negative counterions cannot cross, side 1 acquires a negative interfacial charge and side 2 a positive one. The resulting electrical force opposes further positive-ion transfer. The two contributions balance when the electrochemical potentials are equal; individual ions still cross in both directions, but there is no net electric current of the permeant chemical species.
A simple capacitor model makes the feedback quantitative. If moles of monovalent positive ions move from initially neutral side 1 to side 2, the charge on side 1 is , hence . The free-energy stationarity condition is
As increases, both the concentration imbalance and the voltage drive toward reversal diminish the original diffusion tendency. The equilibrium lies before the bulk concentrations become equal. This demonstrates why a nonzero Nernst potential can coexist with very little change in bulk composition and approximate bulk electroneutrality. A concentration difference alone does not maintain an electrical voltage if both counterions can freely follow each other; selective permeability is essential to this model.
The bacterium samples its surroundings in time. A sustained increase in chemoattractant makes a favorable run persist by reducing the chance of tumbling; declining chemoattractant makes reorientation more likely. Individual tumbles need not point up the gradient. Biasing run durations can generate net migration even if reorientation itself is random. This is temporal comparison, not a direct measurement of a concentration gradient across the short cell body.
Compare the swimming time over length , , with diffusion over the same scale, . Their ratio is the Péclet number . The crossover estimate is
Order-one factors depend on the chosen diffusion-distance convention. A hypothetical stirrer of extent must instead have of order one, giving
These are transport-scale estimates, not a proof of a particular increase in uptake. Around a one-micrometre body, : local diffusive supply is already rapid. Actual uptake enhancement depends on the nutrient field and the cell's uptake boundary condition, neither supplied. The long run can help the organism reach a better nutrient region without implying that it greatly increases the flux from a uniform bath. E. coli normally uses bacterial flagella; the cilium in this comparison is a hypothetical generic stirring appendage, not its actual locomotor organelle.
Model the target as a perfectly absorbing sphere of reaction radius , with a fixed target and a dilute far-field polymerase density . Steady radial diffusion satisfies , , , hence
Therefore, for a rate law multiplying number densities,
For the earliest diffusive target encounter among independent searchers, with independent well-mixed searchers in volume , the approximate first-arrival rate is , giving . The target's effective reaction radius includes encounter geometry and is not uniquely determined by its contour length.
A numerical estimate needs temperature, cell volume and a capture-geometry convention, none of which is fixed numerically in the source. State representative estimates and . The Stokes–Einstein estimate for a sphere of radius in the stipulated viscosity is
Using the target contour scale as an order-of-magnitude gives and
This is an idealized -second-scale estimate, not a unique geometric prediction. Finite cell shape, crowded motion, accessibility and orientational restrictions can reduce the ideal capture rate.
If nonspecific binding merely sequesters searchers and their free fraction is , the three-dimensional capture hazard is reduced fivefold: for those estimates. If bound molecules instead search productively by sliding, a free-time fraction alone cannot determine the total arrival time or its sign of change. That is the distinct facilitated-search mechanism treated in c.
A conceptual single-molecule measurement follows a labelled bound molecule along an independently mapped DNA contour. Project fitted positions onto contour distance and estimate the mean-squared displacement at several lag times. For unbiased one-dimensional diffusion,
Independent localization errors of variance add an intercept . Finite exposure adds motion blur, and drift, binding transitions and the finite DNA contour must be distinguished from free sliding. A slope fit with those limitations can estimate ; simply fitting all observed displacements to a noiseless line can bias it.
Retain the printed coefficient. The original PDF really prints , not a negative exponent. In sliding diffusivity with per-step dissociation, for a basepair step , an unbiased nearest-neighbour walk has step interval . If the dissociation probability per step is , the mean number of steps is approximately , giving
For the printed values, and . The latter is a random-walk displacement, not the contour distance travelled in a million steps. The step interval is about , even shorter than the light-travel time across a basepair by many orders of magnitude. This is an unphysical extrapolation for a protein, not a meaningful experimentally measurable diffusivity.
Consequently no fluorescence camera could temporally resolve trajectories under the stated formal values, regardless of the spatial excursion's length. An actual corrected exponent cannot be inferred uniquely from the paper. If a physically plausible were supplied instead, resolving sliding would require exposures substantially shorter than and diffusion displacements larger than localization error, with enough detected photons. The dependence on is explicit above, so no guessed source correction is needed.
PALM separates emitters in time. A sparse subset of photoactivatable or photoconvertible labels is switched on so that its diffraction-broadened images can be fitted individually; repeated activation, imaging and localization builds a map from many isolated detections. The positional uncertainty can be much smaller than the point spread function width because many photons estimate one image centre: ideally for image width and detected photon number .
That is fluorescence localization precision, not automatically the structural resolution of the reconstructed cell. Pixelation, background, imperfect optics, drift, label linkage, blinking and labelling density matter. Repeated localizations of one molecule cannot fill unlabelled gaps; a crowded frame with overlapping images cannot simply be treated as many isolated emitters. For a static structural map, the structure must remain sufficiently stable while the many frames are accumulated, or its time dependence must be reconstructed explicitly.
For mobility in living cells, keep activated proteins observable for successive frames and link isolated positions into trajectories before bleaching. Repeating this for successive sparse subsets yields many trajectories without making the entire protein population simultaneously fluorescent. This is the sptPALM strategy.[ref]
For free isotropic two-dimensional diffusion with per-coordinate localization error and instantaneous exposures,
Fit several lags or a full trajectory likelihood to infer D and its uncertainty. Finite exposure causes motion blur; for a uniform shutter of duration , a simple correction is . Confinement, directed motion and switching binding states change the lag dependence. Short tracks, missed detections and false linking in crowded regions can mimic apparent heterogeneity, so mobility classes require explicit measurement-error and tracking controls. A trajectory map recovers dynamics; it is not identical to a static PALM density map.
The given sphere has radius . The Stokes–Einstein relation at an illustrative gives
For well-mixed independent searchers in cell volume , the target-arrival rate is . The corresponding mean first arrival time is
The source does not fix , the reaction radius or temperature. For an order-of-magnitude estimate, a bacterial volume around and reaction radius – give about –. A 20-base-pair B-form site is about long, supporting a nanometre-scale capture radius but not uniquely defining its spherical value. The estimate is the ideal diffusion-limited time, not a measured intracellular search time: nonspecific binding, orientational matching, crowding and limited accessibility are excluded.
The term is the bulk-diffusion search cost: increasing the effective length of DNA searched after each encounter enlarges the useful capture region and reduces this contribution. The term is the accumulated sliding cost. Approximately DNA landing-and-sliding episodes, each taking order , give . These are scaling estimates; geometric constants and correlated revisits are omitted by the supplied model.
Thus the facilitated target search is slow as , because too little sequence is examined per encounter, and slow at very large , because sliding and repeated local exploration take too long. A finite intermediate minimum balances 1D scanning against 3D relocation. The right panel of the figure in solution shows the two contributions and their sum. The model requires a useful range of sliding lengths below the full chromosome length, finite positive diffusion coefficients and effective relocation between episodes; it is not a claim that arbitrarily large is physically available.
Take forward as increasing , and specify the sinusoidal amplitude convention by
Thus the unforced barrier is , not . Let and . Stationary points exist because . A minimum is at and its forward saddle at . The distances and activation barriers are
Use an Arrhenius estimate with the supplied distance-squared waiting-time scaling, . Calibrate at with and rate in each direction. At the explicit room-temperature estimate , , giving
Applying the same distance rule separately backwards gives and hence the literal scaling-model estimate
Thermodynamic qualification: for isothermal diffusion in this translated tilted potential, adjacent wells differ by , so local detailed balance instead requires . The two separate distance prefactors introduce an extra factor , inconsistent with that requirement. Keeping the forward estimate while enforcing the proper ratio gives
Both interpretations are shown because the supplied distance heuristic is not, by itself, a thermodynamically consistent prescription for both directional prefactors. Additional hidden chemical driving would require its own affinity, not an unexplained prefactor asymmetry. A common transition-state prefactor, as in a consistent Kramers model, preserves the Boltzmann rate ratio. Moreover the forward barrier is only about , so rare-escape exponential dwell kinetics and the numerical rate should be regarded as rough estimates. If “amplitude” meant the entire unforced barrier height, the maximum restoring force would be only and no dwell wells would survive at ; the amplitude convention is therefore essential.
A charged biological molecule attracts counterions and repels ions of the same sign. Thermal diffusion opposes complete localization, producing a mobile ionic atmosphere whose charge partly cancels the molecule's electric field at larger distances. Solvent polarization also reduces the field through its permittivity, but that dielectric response is distinct from the mobile-ion screening cloud.
For an ideal electrolyte in equilibrium at temperature , an ion of valence has local number density
For , linearization about an electroneutral bulk gives
This is the Debye length. Here counts each species separately and is the absolute solvent permittivity. For a point charge , outside its immediate molecular environment the linearized solution is
It is a screened potential, not a claim that the molecular charge has disappeared.
In terms of ionic strength with in , the density sum is , so
In water near , : about at and at . Increasing ion density or valence shortens the range; temperature and solvent permittivity also matter. The displayed temperature dependence holds at fixed permittivity, which itself changes with temperature in water.
The derivation assumes dilute mobile ions, a continuum uniform dielectric and weak potentials. Near highly charged DNA or protein surfaces, nonlinear Poisson-Boltzmann equation behaviour, counterion condensation, finite ion sizes, specific binding and correlations can matter. At physiological salt the predicted length approaches molecular scales, so the formula is a useful far-field scale with limited microscopic accuracy, rather than a complete account of the first hydration layer.
The small binding RNA can serve conceptually as an RNA aptamer recognition probe. Coupling recognition to an optical label would make a specifically bound virus appear as a moving fluorescent spot. Single-particle tracking estimates successive spot positions and fits their displacement statistics. Specificity is assessed by comparison with a nonbinding recognition probe or a sample without the target, while the much faster diffusion of unbound probes distinguishes freely moving label from target-associated signal. Subnanomolar dissociation constant supports appreciable occupancy at low probe abundance, but affinity alone does not determine the dissociation lifetime because . This describes detection principles, not an operational virus preparation or labelling protocol.
For a dilute spherical particle of radius in bulk water, the Stokes–Einstein relation gives, at an explicitly assumed ,
A free probe would have , roughly thirty times larger. The recognition probe changes the target's hydrodynamic size only modestly in this estimate.
For isotropic Brownian motion, the mean vector displacement is zero. In , distinct measures of how far it moves are
A two-dimensional image instead gives projected root-mean-square displacement and mean projected distance . The projected mean squared displacement is ; independent localization errors of variance in each coordinate add . Finite exposure blurs motion, and nearby walls reduce diffusion; those effects must be separated from particle size when interpreting a trajectory.
A small swimming member of the bacteria samples too little spatial separation across its own body to compare weak chemical gradients accurately. Temporal sensing in bacterial chemotaxis instead compares recent receptor activity with an adapted memory. Straight runs carry it between locations; tumbles reorient it. When conditions improve, tumbling becomes less frequent and favorable runs last longer. An unfavorable direction is abandoned sooner. Repetition creates a biased random walk without requiring the cell to identify the best direction at each tumble. Rotational diffusion and environmental obstacles eventually destroy the benefit of maintaining a direction indefinitely.
An order-of-magnitude transport estimate compares the time to swim a distance with the time for an attractant to diffuse that distance. Their crossover is a Péclet number of order one:
This estimates a useful run scale: over much shorter distances, nutrient diffusion rapidly replenishes the local environment compared with the organism's displacement. Runs of order a second also permit comparison along a substantial path before reorientation. The given and determine this scale, not an exact biological mean run length or a complete sensory optimization. If a one-coordinate root-mean-square diffusion length is equated with , the resulting prefactor is two; a three-dimensional root-mean-square convention gives six. Such constants cannot be fixed by a scaling argument alone.
For a cilium of length , compare its advective time with the local diffusion time . Significant stirring relative to diffusion requires , so
This is about times the stated swimming speed. If the tip travelled on a circular path of radius , that speed would correspond to an illustrative frequency . That geometric conversion is not a specified beat pattern. At microscopic scales the flow has small Reynolds number; useful transport depends on the actual nonreciprocal motion and resulting flow, and rapid motion alone does not guarantee increased uptake. The scaling explains why simply stirring with a very short appendage is demanding when molecular diffusion is already fast.
Treat the target as an absorbing sphere of effective capture radius , surrounded by a uniform dilute reservoir of RNA polymerase with number density . The steady diffusion equation, spherical symmetry, and boundary conditions give
By Fick's laws, the inward capture rate per target is . Multiplying by the target number density and comparing with the bimolecular rate law therefore yields the Smoluchowski association rate
Its units here are , since the concentrations count particles rather than moles. The absorbing boundary assumes every encounter is productive; orientation constraints would reduce the effective association rate constant.
The Stokes–Einstein relation, for RNA polymerase radius at an illustrative temperature , gives
Neither cell volume nor the effective is specified, so a unique numerical time is unavailable. A representative Escherichia coli volume is . The -base pair target has a DNA contour extent of about ; using this as an order-of-magnitude capture radius gives . For independent searchers, the target's approximately constant encounter hazard is , hence
Using half the target's contour extent for instead gives about ; the spherical capture model is not exact target geometry. One searcher's mean would be about under the first convention, whereas the first arrival among is much faster.
If only a fraction of each searcher's time is available for three-dimensional diffusion, and nonspecific DNA binding neither scans nor delivers the target, the effective searching population is . Thus , about under the stated estimate. This assumes stationarity, well-mixed searchers and no additional target-access route; nonspecific binding alone is sequestration, not facilitated diffusion.
Write and . The search time is : it decreases from infinity at small positive , reaches one minimum, and increases without bound at large within this idealized model.
During a facilitated diffusion search, roughly separate DNA landing-and-sliding episodes sample a chromosome. A three-dimensional encounter with any chromosomal landing region has characteristic time , so the total bulk-search contribution is . Each one-dimensional diffusion excursion lasts order ; multiplying by gives . Longer sliding therefore reduces the number of bulk relocations but increases time spent repeatedly exploring local DNA. Numerical constants and overlap corrections are already omitted in the given model.
Differentiating and using the positive second derivative gives
At the two time contributions are equal. The optimum is meaningful only if lies between the microscopic capture scale and the available chromosomal contour length.
A decrease in DNA concentration means a decrease in , but its effect on the time is not determined by that ratio alone. At fixed chromosome length , dilution by increasing increases at fixed and increases as ; bulk relocation becomes slower. If instead decreases at fixed , while the target is retained, the stated model predicts a smaller sliding contribution and a smaller optimized time, proportional to . In both cases . This distinction prevents attributing a unique time trend to concentration without specifying how it changes.
Figure 1.
Search-time optimum and position-dependent molecular motor binding and force
. Details.
The left panel normalizes the search law to , . The other panels illustrate b and c.