Start both systems from zero and compare their times to reach half of their own steady concentrations. Without regulation,
In the strong logic approximation for gene regulation, production is initially fully on and the concentration rises towards the repression threshold. The regulated level is pinned near , so it reaches before switching:
The faster relative response comes from allowing high initial production while keeping the final abundance low. Negative autoregulation can reduce response delays and sensitivity to production fluctuations, and can reduce unnecessary accumulation of a protein. It does not suppress every possible noise source at every frequency.
A perfectly discontinuous on/off right-hand side has a threshold-switching ambiguity rather than an ordinary smooth equilibrium at . Here threshold pinning means the limiting rapid-switching or steep smooth-response idealization. A continuous repression law resolves that idealization. Comparing half of the same prescribed concentration, rather than half of each model's final concentration, would be a different response-time comparison.
For the strong Hill autorepression response, in the requested approximation, . Put . Then
so the formal zero-initial-concentration solution is
Setting gives
The formal ratio decreases with positive ; for it is about , and for about . Its large- asymptote is . This is a nonlinear half-response time, not simply the near-equilibrium relaxation time .
Approximation qualification. The omitted initial layer is precisely the one the source permits us to ignore: the approximate derivative diverges at , whereas the full Hill equation has initial slope . The approximation is not uniform as at fixed . Strong repression already at half of the final concentration needs . The actual Hill equation obeys , so . Its steep-switch limit approaches the logic-response scale in a; the formal vanishing large- result above cannot be extrapolated past its regime of validity.
For the parameter sensitivity of autorepressed expression, unregulated production has , so its logarithmic sensitivity is one. In the strong Hill approximation, , with held fixed. Hence
The ideal logic level instead gives , because changing production primarily changes the switching duty fraction, not the threshold. Finite Hill repression has a residual production sensitivity.
The exact Hill steady-state relation makes the comparison continuous. Write ; then and . Thus
It tends to one for weak repression and to for strong repression. The calculation concerns a local relative parameter change, not the response to an arbitrary large production perturbation.
Use steady or initial-rate conditions with effectively clamped substrate concentrations. For either substrate , put and . The intermediate balances are
The product flux is therefore
The free enzyme concentration is common to the two competing pathways, so it cancels even when enzyme sequestration changes its magnitude. Consequently
Each post-binding factor is the ratio of probabilities of proceeding rather than dissociating at that stage. The expression assumes positive common forward rates and no significant product back reaction. Without the steady/intermediate-balance or initial-rate assumption, the ratio can depend on intermediate populations and observation time; binding and dissociation constants alone would not determine a unique transient ratio.
With equal substrate abundance and association rates, let and . Then
Both factors increase with their forward rate: for example . Thus maximum discrimination requires enough waiting at each stage for the less stable incorrect complex to be rejected. In the slow-forward limit,
The two rejection stages multiply their discrimination. Fast modification or fast product formation bypasses the corresponding dissociation test and drives that factor towards one.
There is a kinetic proofreading fidelity-throughput tradeoff: taking either forward rate literally to zero gives zero product flux, not a useful finite-output optimum. If dissociation rates are well separated, a practical high-yield discrimination window is and . Correct complexes then usually pass both stages while incorrect complexes usually dissociate; its fidelity is not the mathematical infimum obtained with vanishing forward rates. The irreversible modification represents a driven, energy-consuming step, permitting successive nonequilibrium checks rather than merely repeating one equilibrium binding comparison. No finite positive optimum can be selected from minimization of alone without a speed or yield requirement.
A skeletal-muscle fibre contains parallel myofibrils assembled from repeating sarcomeres between Z discs. Thin actin filaments extend from the Z discs; bipolar thick myosin II filaments occupy the central region. Overlap between them allows many myosin heads to interact with actin in parallel.
A rise in cytosolic calcium exposes actin binding sites through the troponin–tropomyosin regulatory system. ATP binding releases a strongly attached myosin head; hydrolysis primes its lever arm, attachment and product release permit a force-generating stroke, and another ATP allows detachment. The cycle converts chemical free energy into mechanical work rather than consuming actin or shortening the filaments themselves.
Relative filament sliding draws the Z discs together. The I band and H zone shorten as overlap increases, whereas the thick-filament A-band length is approximately unchanged. Calcium removal restores inhibition and permits relaxation. The many-head ensemble sustains force even though individual cross-bridges repeatedly attach and detach. This sliding-filament description concerns myosin II in muscle; the processive two-head walking of Myosin V in a is a related but distinct motor arrangement.
In spatially gated myosin attachment, for a freshly available head encountering an isolated site, take before the attachment window. Converting time derivatives to position derivatives with gives in the binding window, in the protected bound region, and after detachment becomes allowed. Continuity gives
For , the plateau exceeds and approaches one for much slower motion. For , the binding probability is lower, becoming at very high speed. The spatial detachment tail has length , so faster sliding makes a longer low-probability tail. The following is an original calculation of those profiles; the force panel is derived in c.
Figure 1.
Spatial attachment probability and force in the prescribed-velocity myosin model
. Original profiles compare velocities below and above the binding-window crossing-rate scale. Binding rises only inside the narrow window, remains constant up to the detachment threshold, and then decays. The second panel shows the resulting signed mean force in the narrow-window approximation. Details.
Reset condition. Geometric separation does not alone guarantee detachment before the next site: independently available encounters also need . Without that extra timescale separation the head may still be attached to an earlier site.
For completeness, the periodic steady-state correction can be written explicitly. Put , , and
Here is the probability of entering the new window still attached to an earlier site. In the window the new-site probability is ; after it, use its value in the plateau and tail above. For , replace by . This follows from the old-site probability and new attachment rate . The isolated-site result is recovered as .
For the work of a spatially gated myosin head, the mechanical force exerted by an attached head in the increasing- direction is . It assists motion before and resists it afterwards. For one site encounter the expected work is
In the narrow binding-window limit, retain the capture probability but neglect the attachment displacement compared with . The positive work is . After , use to obtain the negative contribution . Hence
An independent energy calculation gives the same factor: binding near zero loads energy , while an exponentially distributed post-threshold attachment time has , leaving mean elastic energy to be lost upon detachment. The factor of two in that second moment is essential.
One site passes every distance , so average work per distance, or power divided by , gives
The force tends to as and changes sign at in this idealized model. Above that speed forced sliding is resisted on average. The formula uses the reset condition in b, as well as ; if reset fails, replace by the periodic capture probability to leading order in . For arbitrary window width, the isolated-site work is exactly . The loaded spring energy must be supplied by the binding/power-stroke chemistry; this model does not specify an ATP efficiency. Real force-dependent transition rates and energy budgets limit extrapolation of the toy force-velocity law.
In spectrin unfolding by force spectroscopy, the cantilever measures tension through its calibrated deflection, while displacement of the pulling point increases the separation of the protein attachments. The protein extension is obtained after allowing for cantilever compliance and any linker extension; it is not simply the commanded displacement.
A folded spectrin repeat stores much of its backbone in a compact domain. Between unfolding events, the already released chain is stretched: increasingly straight conformations have fewer accessible configurations, and the resisting entropic force rises nonlinearly. A peak is reached when a still-folded domain crosses an unfolding barrier on the timescale of the imposed loading. Its released backbone increases the available contour length abruptly. At nearly the same end-to-end distance the fractional extension is then smaller, so tension drops. Pulling resumes on the next, longer-chain elastic branch. Repeated release of domains produces the saw-tooth sequence.
The record is therefore evidence for sequential mechanical unfolding of modular structure, not repeated stretching of a fixed-length rod. For the nominal eight-domain construct, successive branches correspond to increasing numbers of unfolded repeats. The order and rupture forces can fluctuate because barrier crossing is thermally activated and depends on loading rate, structural heterogeneity and pulling geometry; equivalent repeats need not rupture at identical force.
The horizontal separation between worm-like-chain fits at the same force estimates the added contour length, whereas the vertical force drop at one extension reflects redistribution of tension. These are distinct measurements. A repeat increment near the specified is consistent with repeated similar structural units. The gain is a net release relative to the folded span and attachment geometry; is not automatically the universal contour length of a peptide residue.
During unloading, relaxation and possible refolding can restore compact domains. Hysteresis would indicate dissipative or kinetically delayed processes. Neither the peak force alone nor an arbitrary saw-tooth area directly equals an equilibrium unfolding free energy: one must separate chain elasticity, cantilever storage, irreversible loss and pulling-rate effects. The worm-like-chain branches describe released-chain elasticity, not the molecular transition state responsible for a rupture peak.
Write and . The worm-like chain force-extension interpolation is
A fit to an elastic rising branch estimates its contour length and, if not fixed independently, its persistence length . Comparing successive fitted contour lengths gives the length released by a domain and helps distinguish whole-repeat unfolding from partial events. estimates the flexibility of the unfolded chain and the bending rigidity ; at a stated temperature, .
For small extension, expansion gives , so the effective entropic spring constant is . Close to full extension, , giving strong strain stiffening. These regimes constrain the fit differently and correlate its parameters if the measured interval is short. The fit does not by itself identify a unique folding barrier or unfolding kinetic rate.
Persistence-length convention. In a three-dimensional isotropic worm-like chain the standard definition is . Thus at , . This is not literally the mean angle . For a simple counterexample, equally probable angles and have mean angle but . The source's angle description is shorthand for the tangent-correlation convention.
After a domain unfolds, its backbone joins the same continuous chain as previously unfolded material. The new branch differs mainly by an added contour length, while the local backbone bending statistics remain similar. It is therefore consistent to fit successive branches with different and one common . Tension is transmitted through serial portions of the same chain; one should not add their tensile forces as though they were independent parallel springs.
The persistence length is short compared with the contour increment per repeat, so each released domain contributes many flexible bending segments. Stretching suppresses those thermal modes and raises free energy through loss of entropy, which is precisely the worm-like-chain mechanism. The common is a property of the unfolded backbone under the measurement conditions, not the folded domain's axial length or the number of domains already released. With only as an illustrative value, .
A common persistence-length fit is evidence for approximately repeat-independent chain flexibility over the observed branches; it is not proof that a continuum inextensible model describes every molecular scale. Residual folded structure, heterogeneous linkers, force-induced backbone stretching, and very short-chain effects can cause deviations. The formal force divergence at signals the model's inextensibility assumption, not a physical infinite tensile force.
The sodium trace is a rapid inward current after depolarization, followed by loss of current even while voltage remains depolarized. A voltage-gated sodium channel therefore activates rapidly and then becomes unavailable through inactivation. Under ordinary cardiac ionic gradients, the test potential is below the sodium reversal potential, so opening carries inward positive charge. Rapid sodium entry supports the fast action-potential upstroke and propagation in working myocardium; inactivation and subsequent recovery contribute to refractoriness.
The potassium trace develops more slowly, is outward during the positive step, and remains approximately sustained until repolarization. It is consistent with a delayed outward potassium conductance: activation counteracts depolarization, and deactivation follows the return towards the holding voltage. Potassium currents contribute to repolarization and stable diastolic voltage. The trace does not identify a unique cardiac potassium-channel subtype, and not all potassium channels lack inactivation.
Together, fast inward activation initiates excitation and slower outward activation terminates it. Current direction is set by electrochemical driving force as well as gating: in an ohmic approximation. A voltage clamp fixes and reveals those currents; in a freely firing cell the currents instead jointly determine the voltage through the membrane-capacitance equation. Calcium current is also important for the cardiac plateau, so these two traces alone are not a complete ventricular action-potential model. The sodium-potassium pump maintains long-term gradients but is not the fast gating mechanism seen in these recordings.
Voltage gating and ion selectivity are different structural functions. A typical voltage-gated potassium channel has four pore-forming subunits; their voltage-sensing regions move charged residues in the membrane electric field and couple that motion to a pore gate. Carbonyl coordination in the selectivity filter compensates potassium dehydration in an ion-specific geometry. The smaller sodium ion is not automatically passed more readily: its hydration and coordination energetics differ. Structural work on a mammalian Kv channel resolves the sensor/pore architecture and an open activation gate.[ref]
A sodium-channel alpha subunit contains four homologous domains assembled into one pore, with separate voltage sensors. Rapid activation and fast inactivation are distinguishable conformational transitions; an intracellular inactivation apparatus stabilizes nonconducting states after activation. Repolarization favours recovery of availability. A channel can therefore be closed but activatable, or inactivated and temporarily unavailable, despite zero current in both cases.
Several types of evidence test these ideas. Ionic reversal potentials shift when ionic gradients are changed, establishing selectivity and driving force. Single-channel recordings show discrete conducting levels and stochastic dwell times; whole-cell current changes often reflect changed open probability rather than a continuously varying single-channel current. Small gating-charge transients can remain when ionic permeation is blocked, demonstrating sensor movement rather than ion flux. Conditioning and paired-pulse records separate activation, inactivation and recovery. Structural perturbations of voltage sensors and of inactivation machinery alter different kinetic properties, supporting their functional distinction.
These observations constrain models but do not make an individual macroscopic trace a unique structural reconstruction. Different distributions of closed and inactivated states can produce similar averages. In particular, the sustained potassium trace supports little inactivation on the displayed interval; it does not justify assigning every potassium channel a universal absence of an inactivation mechanism.
An increased late inward sodium current can delay repolarization if fast inactivation is incomplete or slowed, if channels reopen excessively, or if the voltage range of persistent/window current is enlarged. A larger brief peak current alone is not equivalent to a sustained depolarizing current. Conversely, reduced outward potassium current can prolong the action potential through fewer functional surface channels, reduced conductance, slower or less favourable activation, or altered deactivation/availability. Either change shifts the balance towards inward current during repolarization.
These are distinct biophysical routes to Long QT syndrome, not claims that every subtype has the same channel defect.[ref] The experimental association of sodium-channel abnormalities with inherited disease supports the inward-current mechanism.[ref]
At fixed voltage let and be opening and closing rates. The open probability obeys
This illustrates equal open probability with unequal ion-channel relaxation times: both channels have , but channel A has and channel B has . Starting closed, and . If channel number, unitary conductance and driving force agree, both have the same steady macroscopic current; B approaches it ten times faster.
An individual channel remains either open or closed. Its mean closed and open dwell times are and , respectively: both are in A and in B. Thus identical equilibrium open probability does not imply identical kinetics or noise correlation time. Ensemble probability relaxation depends on the sum of rates, while the mean dwell in a given state depends on its exit rate.
Figure 1.
Equal equilibrium opening probabilities can conceal different channel kinetics
. Original calculations show the specified two-state rates giving relaxation times of one and one tenth of a millisecond. Separate illustrative activation-inactivation and delayed-activation curves explain transient inward sodium current and sustained outward potassium current; those illustrative curves are not fits to the source recordings. Details.
A two-state closed–open scheme at fixed voltage relaxes monotonically and cannot by itself produce a sodium-current transient during a maintained step with constant driving force. A minimal sodium scheme therefore includes an unavailable state: , with voltage-dependent activation, inactivation and recovery; closed-state recovery pathways may be added when needed. Opening is initially rapid, occupancy then accumulates in , and recovery on repolarization restores an activatable population.
For the illustrated potassium current, a closed–open model gives sustained opening, but several activation steps can account for its initial delay. A conventional illustrative coarse-grained description is
with each gate following . Fast activation, subsequent decrease of , and slower increase of produce the two qualitative records. The powers are phenomenological gate assumptions, not values uniquely inferred from these traces or an exact structural count for every cardiac channel.
A 1-Hz cardiac rhythm has about a second between successive activations, whereas a 100-Hz neuronal rhythm has only ; the available resting interval is shorter still after subtracting the action potential. Recovery or deactivation with a timescale comparable to that interval can then leave residual inactivation or potassium activation, giving use-dependent changes and refractoriness. Millisecond or faster relaxation usually completes much more readily during a long diastolic interval. The particular two-state A and B in d would both nearly equilibrate over a full at their stipulated constant rates; their values alone do not demonstrate high-frequency failure. Real comparison requires voltage-dependent recovery rates and the different cardiac and neuronal waveforms. A frozen-rate model cannot quantitatively predict either firing pattern from frequency alone.
At initiation, RNA polymerase recognizes a promoter and forms a closed bound complex; bacterial promoter recognition involves its specificity factor. Local strand opening exposes the template in an open complex, after which initial RNA synthesis can precede promoter escape into elongation. Binding, opening and productive escape are distinct transitions.
For DNA unwinding measured by magnetic tweezers, magnetic tweezers report the mechanics of a torsionally constrained DNA tether through the position of an attached magnetic bead. Axial magnetic force extends the tether, while imposed rotation controls its linking constraint. Video localization measures extension changes. The important physical relation is
Opening a small helical region decreases twist, so the rest of the tether acquires compensating positive writhe. On a negatively supercoiled tether this removes some writhe magnitude and raises its extension; on a positively supercoiled tether it increases writhe magnitude and lowers extension. A change with opposite signs under the two torsional conditions therefore distinguishes unwinding from simple compaction, which would shorten both tethers. Calibrating extension against imposed turns converts a mechanical step into a change in helical winding. This is a conceptual readout of local opening, rather than a direct optical image of separated bases.[ref]
Figure 1.
Torsional constraint converts local DNA opening into a bead-extension signal
. An original conceptual magnetic-tweezers schematic is paired with illustrative extension-versus-turns curves. Local unwinding shifts the curve and produces opposite extension changes on its negatively and positively supercoiled branches. The diagram does not reproduce the source experiment or specify a laboratory protocol. Details.
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.
Write and . The given mean-search scaling is . Three-dimensional relocation sees an effective capture length of order , giving a bulk encounter time . A sliding excursion explores length of order in time of order ; sampling order contour segments gives a total sliding cost of order . Prefactors, correlations and repeated visits are coarse-grained in this scaling model.
Small sliding lengths make relocation inefficient; very long sliding wastes time revisiting a local DNA region rather than relocating to another. The optimum balances both contributions:
Since , it is the unique minimum for unrestricted positive . At the optimum, bulk and sliding times are equal. With ,
Figure 1.
Bulk relocation and local sliding impose opposing search-time costs
. The original dimensionless plot shows a decreasing inverse-sliding-length contribution and an increasing sliding contribution. Their sum has a minimum where the two costs agree, at a sliding length equal to the analytic optimum. Details.
The model is meaningful only over an accessible range between molecular steps and the available chromosome length. An optimum outside that range is replaced by a boundary choice; the unphysical divergences at arbitrarily small or large describe limits of the scaling expression. A common multiplicative searcher-number factor would not change the optimum, but an exact many-searcher first-arrival mean is not generally obtained by dividing an arbitrary single-searcher mean by their number without a waiting-time assumption.
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.
A single fluorophore creates a diffraction-broadened point-spread function, but its centre can be fitted much more precisely than that function's width. With Gaussian width and detected independent photons, the shot-noise limit is . Pixelation, background photons, dye orientation, motion blur, drift and label-to-target geometry impose additional uncertainty and systematic error. Localization precision is therefore not identical to either optical resolution or position accuracy.[ref]
In myosin-V tracking, localizing a label on one head or its lever arm distinguishes the paths of individual motor parts from the average motion of the whole motor. A trailing head that becomes the leading head travels roughly twice the stalk advance: about versus a roughly motor stride. A near-head label can therefore alternate a large movement with a pause while the partner head advances. Other lever-arm label positions give unequal alternating steps whose sum is approximately twice the stride. Those patterns support hand-over-hand walking rather than a simple inchworm model in which the same head remains in front. The label geometry is needed when interpreting individual step sizes.[ref]
Figure 1.
A fitted diffraction spot can locate alternating myosin heads beyond its optical width
. An original Gaussian spot illustrates the distinction between point-spread width and centre precision. A separate schematic shows a labelled head alternating a roughly 74-nanometre movement and a pause while the two-head motor advances in roughly 37-nanometre strides. It is a conceptual calculation, not copied tracking data. Details.
Super-resolution localization microscopy extends the same principle to a map containing many molecules by separating otherwise overlapping emitters in time. Sparse activation, blinking or transient binding yields individually fitted positions, which are accumulated. Dense sampling, registration and labelling specificity determine the final structural resolution; merely taking more photons from an unresolved simultaneous pair does not make it a single isolated emitter.
In PAINT, fluorescent probes bind transiently so that separated binding events can be localized and accumulated. DNA-PAINT uses reversible complementary-strand recognition at docking sites on a DNA nanostructure. The bound probe supplies a localized spot while freely diffusing probes contribute background; repeated visits build up a positional map, with replenishment rather than permanent reliance on one fluorophore. This is an imaging principle, not a requirement for a particular target sequence or preparation protocol.[ref]
Usable biological free energy has several interconvertible forms. At approximately fixed temperature and pressure, chemical reactions are described by chemical potentials, , and . ATP hydrolysis, for example, has under the specified biochemical standard-state convention. The activities and proton/ion convention matter. Breaking a bond alone costs energy; the net free-energy change includes the products, solvent and all bond changes, not a store released simply by snapping one bond.
Redox reactions provide free energy when electrons move between donors and acceptors with a favourable potential difference: . Electron-transfer chains can couple that drop to ion transport, building an electrochemical gradient. For transfer from an outer to an inner reservoir,
Both concentration difference and electrical potential contribute. ATP synthase couples downhill ion flow to nucleotide free energy; driven pumps can perform the reverse conversion and maintain gradients used by transporters and excitability.
Mechanical work appears as or torque times angle, and can be stored elastically, for example . Molecular motors convert chemical free energy into movement or force; mechanically driven reactions can conversely reshape chemical transitions. Electrical energy can also be stored in membrane capacitance, while bulk ionic concentration gradients store additional chemical free energy. Light supplies photon energy , which photosynthetic machinery can convert into redox and chemical potential differences.
Coupling is essential: an overall favourable reaction can drive an unfavourable one only when their processes are linked, rather than merely present in the same cell. At fixed conditions the extracted work is bounded by the consumed free-energy decrease. Dissipation produces heat and entropy, so energy conservation does not mean all energy remains available for work. Thermal fluctuations can drive random motion, but a persistent directional cycle requires a nonequilibrium free-energy source; an equilibrium asymmetric structure alone is not a free motor.
A point mutation changes a nucleotide at one genomic position. Replication errors and chemical damage can generate substitutions; repair and other proofreading processes influence the retained mutation rate. A coding change may be synonymous or alter an amino acid, while a noncoding change can affect regulation. Its evolutionary effect depends on context, so mutation, functional consequence and eventual fixation are distinct events.
Mutation supplies new variants, while genetic drift, selection and population history determine their persistence. Under a simple neutral model, increasing population size raises the number of new variants but lowers each one's fixation probability; the neutral substitution rate is approximately the mutation rate. Selective constraint, generation time and lineage-specific processes can therefore complicate a universal molecular clock. Molecular evolution is not directed towards anticipated needs.
To infer a phylogeny, compare homologous positions in appropriately chosen orthologous sequences and estimate the tree and branch lengths under a substitution model. Shared derived characters support common ancestry, but reversals and parallel substitutions mean raw mismatches underestimate evolutionary distance. In a simple equal-rates four-state illustration, a mismatch fraction gives corrected pairwise distance . The formula's strong divergence as reflects saturation and loss of reliable historical information, not arbitrarily precise evidence for an ancient split.
The molecular-clock calibration and rate-time identifiability issue is fundamental: branch lengths first measure expected substitutions per site. A clock relates them to time by ; for two equal-rate lineages separated for time , pairwise distance is , so . Absolute times require external information, such as independently dated fossils or geological constraints, or informative sampling dates when the data contain a temporal signal. Sequence differences alone generally identify only rate-times-time, not both quantities separately.
Time-calibrated inference should propagate uncertainty in topology, substitution rates, ancestral-state reconstruction and calibration bounds. Relaxed clocks permit different lineage rates. Fossil minima are not automatically exact split dates; orthology mistakes, incomplete lineage sorting, recombination, horizontal transfer and selection can cause discordance between one gene tree and species history. Multiple independent loci and well-justified calibrations make a stronger historical argument than a single uncalibrated sequence comparison. No unique numerical tree or date can be constructed here without actual sequence observations and calibration data.

Articles by others on the same topic (0)

There are currently no matching articles.