Spontaneous noise reduction in a strongly cooperative reaction model

Noise characteristics of a simple multiparticle reaction model were investigated. The model can describe strongly cooperative reactions and is defined as follows. Particles A react in clusters of size k and each reaction forms a product molecule P. The back reaction is also allowed, and each cluster can dissociate into k reactants A. To describe a situation where the reaction is part of a pathway, the system is made open by assuming that particles A are injected and that particles P decay. This is a continuation study. The model is studied by using the same technique as employed previously, the pair approach reaction noise estimator (PARNES) method. Several new issues have been addressed. (i) In the previous work closeness to the Poisson distribution was used as a noise measure. In this work a more traditional noise measure, the ratio of the variance and the mean, was used to analyze stochastic features of the problem. (ii) The dependence of the new noise measure on k has been analyzed in detail, with an emphasis of investigating reactions with large k values. (iii) The previous study focused on understanding time-dependent issues, while this study focuses on describing the equilibrium state of the system. (iv) An exact solution (published elsewhere), available for closed system only, was used to re-investigate the validity of the PARNES method for describing reactions with large k. It was found that the PARNES method cannot describe fine details of the noise characteristics of such reactions.     

The synergetic enzyme theory of muscular contraction: II Relation between Hill's equations and functions of two-headed myosin.

Based on the assumption of nonidentical two heads of myosin it is pointed out that a strong motive force is generated in actomyosin pair only when ATP-decomposition occurs co-operatively at the both heads and that the tension-independent part of shortening heat is liberated when an ATP molecule is decomposed only at the burst head. These two actions of actomyosin pair are related to the two states of force-generator in Huxley-Simmons' model. Elementary cycles at different positions in a sarcomere are interacted each other through feedback loop via sliding motion of muscular filaments. Due to this synergetic interaction the rate constant for the rate-determining step of elementary cycle has a dependence on velocity v of shortening such as k = k ° + κv. From these functions and properties of actomyosin system in vivo, the following properties of muscle are explained consistently in a quantitative manner: (1) Hill's equation on the relationship between tension and velocity of shortening, (2) damped oscillations in tension and in muscular length around steady state, (3) Hill's energy equation improved in 1964, (4) the chemical equivalence of shortening heat, (5) the influence of tension on the incorporation of radio-active phosphate into ATP and (6) the asymmetric activation by actomyosin system only for the forward reaction, the decomposition of ATP.     

Photosystem I charge separation in the absence of centers A and B. II. ESR spectral characterization of center ‘X’ and correlation with optical signal ‘A2’

The Photosystem I electron acceptor complex was characterized by optical flash photolysis and electron spin resonance (ESR) spectroscopy after treatment of a subchloroplast particle with lithium dodecyl sulfate (LDS). The following properties were observed after 60 s of incubation with 1% LDS followed by rapid freezing. (i) ESR centers A and B were not observed during or after illumination of the sample at 19 K, although the P-700⁺ radical at g = 2.0026 showed a large, reversible light-minus-dark difference signal. (ii) Center ‘X’, characterized by g factors of 2.08, 1.88 and 1.78, exhibited reversible photoreduction at 8 K in the absence of reduced centers A and B. (iii) The backreaction kinetics at 8 K between P-700, observed at g = 2.0026, and center X, observed at g = 1.78, was 0.30 s. (iv) The amplitudes of the reversible g = 2.0026 radical observed at 19 K and the 1.2 ms optical 698 nm transient observed at 298 K were diminished to the same extent when treated with 1% LDS at room temperature for periods of 1 and 45 min. We interpret the strict correlation between the properties and lifetimes of the optical P-700⁺ A⁻₂ reaction pair and the ESR P-700⁺ center X⁻ reaction pair to indicate that signal A₂ and center X represent the same iron-sulfur center in Photosystem I.     

Electron transfer in reaction centers of Rhodopseudomonas sphaeroides. I. Determination of the charge recombination pathway of D+QAQ−B and free energy and kinetic relations between Q−AQB and QAQ−B

The electron-transfer reactions and thermodynamic equilibria involving the quinone acceptor complex in bacterial reaction centers from R. sphaeroides were investigated. The reactions are described by the scheme: We found that the charge recombination pathway of D⁺Q_AQ^?_B proceeds via the intermediate state D⁺Q^?_AQ_B, the direct pathway contributing less than approx. 5% to the observed recombination rate. The method used to obtain this result was based on a comparison of the kinetics predicted for the indirect pathway (given by the product k_AD-times the fraction of reaction centers in the Q^?_AQ_B state) with the observed recombination rate, k^obs_{D⁺ →D}. The kinetic measurements were used to obtain the pH dependence (6.1 ? pH ? 11.7) of the free energy difference between the states Q^?_AQ_B and Q_AQ^?_B. At low pH (less than 9) Q_AQ^?_B is stabilized relative to Q^?_AQ_B by 67 meV, whereas at high pH Q^?_AQ_B is energetically favored. Both Q^?_A and Q^?_B associate with a proton, with pK values of 9.8 and 11.3, respectively. The stronger interaction of the proton with Q^?_B provides the driving force for the forward electron transfer.     

Expression of the I-E target antigen for T-cell killing requires two genes

TheH?2I ^k region encodes at least two different target antigens for unrestricted T-cell mediated killing. The first is controlled by theI?A region alone and the second depends on a pair of alleles, one located to the left ofI?B (presumably inI?A) and the other to the right ofI?J (presumably inI?E). Hence, effector cells nominally specific for a product of theI?E region do not kill target cells with the sameI?E region as the stimulator unless theI?A region is also shared. Some effectors specific forH?2I ^k , such as A.TH anti-A.TL and B10.A(4R) anti-B10.A(2R), cross-react with B10.A(3R) and B10.A(5R) target cells. A product of theH?2 ^b haplotype was shown to complement products of theH?2 ^d orH?2 ^k haplotypes in forming this cross-reactive determinant. The results are consistent with recent biochemical data on the component chains of Ia antigens.     

Stochastic Oscillations Induced by Intrinsic Fluctuations in a Self-Repressing Gene

Biochemical reaction networks are subjected to large fluctuations attributable to small molecule numbers, yet underlie reliable biological functions. Thus, it is important to understand how regularity can emerge from noise. Here, we study the stochastic dynamics of a self-repressing gene with arbitrarily long or short response time. We find that when the mRNA and protein half-lives are approximately equal to the gene response time, fluctuations can induce relatively regular oscillations in the protein concentration. To gain insight into this phenomenon at the crossroads of determinism and stochasticity, we use an intermediate theoretical approach, based on a moment-closure approximation of the master equation, which allows us to take into account the binary character of gene activity. We thereby obtain differential equations that describe how nonlinearity can feed-back fluctuations into the mean-field equations to trigger oscillations. Finally, our results suggest that the self-repressing Hes1 gene circuit exploits this phenomenon to generate robust oscillations, inasmuch as its time constants satisfy precisely the conditions we have identified.     

SENSORY ADAPTATION AND THE STATIONARY STATE

1. Experiments are described which measure the sensitivity of animals exposed to continued illumination to which they have become adapted. It is shown that the amount of outside light energy necessary to stimulate an adapted animal increases with the intensity of the adapting illumination. 2. The data are analyzed quantitatively in terms of the reversible reaction S ⇌ P + A shown previously to account for the photic sensitivity of these animals. This analysis demonstrates that, though the amount of incident energy necessary for a minimal response varies with the adapting intensity, the actual amount of photochemical decomposition required to set off the sensory mechanism is a constant quantity. 3. The ability of these animals to come into sensory equilibrium with any sustained illumination is accounted for quantitatively by the presence of a stationary state in the reversible photochemical reaction S ⇌ P + A during which the concentrations of the three components are constant. 4. It is shown that the concentrations of these substances at the stationary state are automatically controlled by the outside intensity. Therefore, given the sensory mechanism as a basis, the adaptation of the animals to light and the consequent changes in sensitivity, are determined entirely by the light to which the animals are exposed. 5. Because of the properties of the stationary state, and of the constancy of photochemical decomposition for a minimal effect, it is suggested that the sensory system is not only the traditional receptor system, but is also a protecting layer which stabilizes and buffers the relation between the nervous system and the environment.     

Reaction of the Co(II)-substrate radical pair catalytic intermediate in coenzyme B12-dependent ethanolamine ammonia-lyase in frozen aqueous solution from 190 to 217 K

The decay kinetics of the aminoethanol-generated Co^II-substrate radical pair catalytic intermediate in ethanolamine ammonia-lyase from Salmonella typhimurium have been measured on timescales of <10⁵ s in frozen aqueous solution from 190 to 217 K. X-band continuous-wave electron paramagnetic resonance (EPR) spectroscopy of the disordered samples has been used to continuously monitor the full radical pair EPR spectrum during progress of the decay after temperature step reaction initiation. The decay to a diamagnetic state is complete and no paramagnetic intermediate states are detected. The decay exhibits three kinetic regimes in the measured temperature range, as follows. i), Low temperature range, 190 ≤ T ≤ 207 K: the decay is biexponential with constant fast (0.57 ± 0.04) and slow (0.43 ± 0.04) phase amplitudes. ii), Transition temperature range, 207 < T < 214 K: the amplitude of the slow phase decreases to zero with a compensatory rise in the fast phase amplitude, with increasing temperature. iii), High temperature range, T ≥ 214 K: the decay is monoexponential. The observed first-order rate constants for the monoexponential (k_obs,m) and the fast phase of the biexponential decay (k_obs,f) adhere to the same linear relation on an lnk versus T⁻¹ (Arrhenius) plot. Thus, k_obs,m and k_obs,f correspond to the same apparent Arrhenius prefactor and activation energy (logA_app,f (s⁻¹) = 13.0, E_a,app,f = 15.0 kcal/mol), and therefore, a common decay mechanism. We propose that k_obs,m and k_obs,f represent the native, forward reaction of the substrate through the radical rearrangement step. The slow phase rate constant (k_obs,s) for 190 ≤ T ≤ 207 K obeys a different linear Arrhenius relation (logA_app,s (s⁻¹) = 13.9, E_a,app,s = 16.6 kcal/mol). In the transition temperature range, k_obs,s displays a super-Arrhenius increase with increasing temperature. The change in E_a,app,s with temperature and the narrow range over which it occurs suggest an origin in a liquid/glass or dynamical transition. A discontinuity in the activation barrier for the chemical reaction is not expected in the transition temperature range. Therefore, the transition arises from a change in the properties of the protein. We propose that a protein dynamical contribution to the reaction, which is present above the transition temperature, is lost below the transition temperature, owing to an increase in the activation energy barrier for protein motions that are coupled to the reaction. For both the fast and slow phases of the low temperature decay, the dynamical transition in protein motions that are obligatorily coupled to the reaction of the Co^II-substrate radical pair lies below 190 K.     

Enhanced Regulatory Sequence Prediction Using Gapped k-mer Features

Oligomers of length k, or k-mers, are convenient and widely used features for modeling the properties and functions of DNA and protein sequences. However, k-mers suffer from the inherent limitation that if the parameter k is increased to resolve longer features, the probability of observing any specific k-mer becomes very small, and k-mer counts approach a binary variable, with most k-mers absent and a few present once. Thus, any statistical learning approach using k-mers as features becomes susceptible to noisy training set k-mer frequencies once k becomes large. To address this problem, we introduce alternative feature sets using gapped k-mers, a new classifier, gkm-SVM, and a general method for robust estimation of k-mer frequencies. To make the method applicable to large-scale genome wide applications, we develop an efficient tree data structure for computing the kernel matrix. We show that compared to our original kmer-SVM and alternative approaches, our gkm-SVM predicts functional genomic regulatory elements and tissue specific enhancers with significantly improved accuracy, increasing the precision by up to a factor of two. We then show that gkm-SVM consistently outperforms kmer-SVM on human ENCODE ChIP-seq datasets, and further demonstrate the general utility of our method using a Naïve-Bayes classifier. Although developed for regulatory sequence analysis, these methods can be applied to any sequence classification problem.     

Kinetic analysis of biphasic protein modification reactions cooperative effects

A mathematical treatment of protein modification reactions is presented, and it is shown thai in these cases protein modification is described by a summation of exponential functions of reaction time, the number of exponentials being equal to the number of modified protein species. It is shown that in cases of protein modification cooperativity, there is a strict dependence of the coefficients of the multiexponential modification equation on the constants of the same equation. The conditions necessary for a reduction of a multiexponential protein modification equation to one of a summation of two exponentials only are examined. The possible formulae for the coefficients of a two-exponential-summation equation, used to describe the modification of protein models with two, three or four modifiable residues (as well as some aspects of models with five and six modifiable residues) per protein molecule are derived. It is seen that the number of such coefficients is severely limited. The most frequently obtained formula for the lower stoichiomelric coefficient of a 'wo-exponential-summation equation is Ak_a/(k_a-k_b). where k_b and k_b are the constants of the two exponentials of the equation, and A is a constant. The value most frequently arrived at for A is (n?1)/n, where n is the number of modifiable residues per protein molecule, while values such as 1/n, or a/n (where a is an integer, and also where a < n) are also possible. In most of the cooperative protein modification models worked out, k_a is identical with k_n, viz., k_a is identical with the rate constant for the first stoichiometric protein modification.     

Electron transfer after flash photolysis of mixed-valence carboxycytochrome c oxidase

The light-induced difference spectra of the fully reduced (a³⁺a²⁺₃-CO) complex and the mixed-valence carboxycytochrome c oxidase (a³⁺a²⁺₃-CO) during steady-state illumination and after flash photolysis showed marked differences. The differences appear to be due to electron transfer between the redox centres in the enzyme. The product of the absorbance coefficient and the quantum yield was found to be equal in both enzyme species, both when determined from the rates of photolysis and from the values of the dissociation constants of the cytochrome a²⁺₃-CO complex. This would confirm that the spectral properties of cytochrome a₃ are not affected by the redox state of cytochrome a and Cu_A. When the absorbance changes after photolysis of cytochrome a²⁺₃-CO with a laser flash were followed on a time scale from 1 μs to 1 s in the fully reduced carboxycytochrome c oxidase, only the CO recombination reaction was observed. However, in the mixed-valence enzyme an additional fast absorbance change (k = 7·10³s^?1) was detected. The kinetic difference spectrum of this fast change showed a peak at 415 nm and a trough at 445 nm, corresponding to oxidation of cytochrome a₃. Concomitantly, a decrease of the 830 nm band was observed due to reduction of Cu_A. This demonstrates that in the partially reduced enzyme a pathway is present between Cu_A and the cytochrome a₃-Cu_B pair, via which electrons are transferred rapidly.     

Free energy change accompanying electron transfer from P870 to quinone in Rhodopseudomonas sphaeroides chromatophores

Delayed fluorescence of chromatophores of Rhodopseudomonas sphaeroides was measured to estimate the standard free energy change accompanying the electron transfer from the bacteriochlorophyll dimer (P) to the primary acceptor quinone (Q_A). The chromatophores emitted delayed fluorescence with a lifetime of about 60 ms in the presence of o-phenanthroline. By comparing the intensity of the delayed fluorescence with that of the prompt fluorescence, the standard free energy of the P⁺Q_A^? radical pair was evaluated. It was about 0.87 eV below the level of excited singlet state, P1Q_A, or 0.51 eV above the ground state, PQ_A, independent of pH.     

Reaction of turbulence at the edge and in the center of the plasma column to pulsed impurity injection caused by the sputtering of the wall coating in L-2M stellarator

Impurity injection into plasma caused by the sputtering of the wall coating in the L-2M stellarator during auxiliary electron cyclotron resonance heating leads to a change in the level of plasma density fluctuations with frequencies above 0.25 MHz: suppression of long-wavelength (k _⊥ = 2 cm^–1) density fluctuations in the edge plasma, intensification of short-wavelength (k _⊥ = 30 cm^–1) and long-wavelength (k _⊥ = 1 cm^–1) fluctuations at the midradius of the plasma column, and intensification of short-wavelength fluctuations (k _⊥ = 20 cm^–1) in the plasma center (including the gyroresonance region). At the same time, the level of fluctuations with frequencies below 0.25 MHz remains unchanged. In the edge plasma, a decrease in the plasma potential and suppression of its fluctuations is observed during impurity injection, which also causes an increase in MHD activity.     

Exploration of the conserved A+C wobble pair within the ribosomal peptidyl transferase center using affinity purified mutant ribosomes

Protein synthesis in the ribosome's large subunit occurs within an active site comprised exclusively of RNA. Mutational studies of rRNA active site residues could provide valuable insight into the mechanism of peptide bond formation, but many of these mutations cause a dominant lethal phenotype, which prevents production of the homogeneous mutant ribosomes needed for analysis. We report a general method to affinity purify in vivo assembled 50S ribosomal subunits containing lethal active site mutations via a U1A protein-binding tag inserted onto the 23S rRNA. The expected pH-dependent formation of the A2450+C2063 wobble pair has made it a potential candidate for the pH-dependent conformational change that occurs within the ribosomal active site. Using this approach, the active site A2450+C2063 pair was mutated to the isosteric, but pH-independent, G2450•U2063 wobble pair, and 50S subunits containing the mutations were affinity purified. The G•U mutation caused the adjacent A2451 to become hyper-reactive to dimethylsulfate (DMS) modification in a pH-independent manner. Furthermore, the G•U mutation decreased both the rate of peptide bond formation and the affinity of the post-translocation complex for puromycin. The reaction rate (k_pep) was reduced ~200-fold for both puromycin and the natural aminoacyl-tRNA A-site substrate. The mutations also substantially altered the pH dependence of the reaction. Mutation of this base pair has significant deleterious effects upon peptidyl transferase activity, but because G•U mutation disrupts several tertiary contacts with the wobble pair, the assignment of A2450 as the active site residue with the neutral pK_a important for the peptidyl transferase reaction cannot be fully supported or excluded based upon these data.     

Restriction of bacteriophage lambda by Escherichia coli K

Derivatives of phage lambda, for which the numbers and positions of the recognition sites for endonuclease R. Eco_k are known, were used as substrates for the Escherichia coli K restriction system in vivo and in vitro. A single unmodified recognition site was sufficient for a DNA molecule to be bound and broken by the K restriction enzyme. Although discrete fragments of DNA were not produced, the breaks were made preferentially in the proximity of the recognition site. Breakage of a DNA molecule with only one recognition site required a 10 to 40-fold higher concentration of restriction enzyme than breakage of a DNA molecule with two or more recognition sites, but these substrates were all equally effective in a binding assay for the enzyme.The polynucleotide kinase reaction provided no evidence for new 5′-terminal sequences generated by restriction in vitro; the 5′ termini were either refractory to the polynucleotide kinase reaction or had no sequence specificity.     

A slow fluctuation in the molecular conformation of a soya bean trypsin inhibitor

The kinetics of the hydrogen-deuterium exchange reaction in a trypsin inhibitor (Kunitz) from soya bean have been followed by infrared absorption measurements in aqueous solutions at various temperatures and pH values. It was found that, in every case, 49% of the total peptide hydrogen atoms exchange relatively slowly. This amount corresponds to 83 peptide groups per molecule, and this is considered to be equal to the number of peptide NH groups involved in hydrogen bonds with the carbonyls of other peptide groups in the protein molecule in its native form. Each rate constant (k) determined at pH 2.75 for this category of the NH groups is in good agreement with the value expected from an idea that the breaking of the peptide-peptide hydrogen bonds takes place very slowly, and that this is the rate-determining process in the hydrogen-deuterium exchange reaction. Thus, by ultraviolet absorption measurements at 297 nm, the equilibrium constant of the native and denatured forms has been determined in the temperature range from 42 to 53.5 °C, as well as the reaction rate of reaching equilibrium from an off-equilibrium state. From these data the rate constant (k₁) of the denaturation reaction is determined, and the k₁ value is found to be practically equal to the hydrogen exchange rate constant (k). The Arrhenius plot of this rate constant (k) gives a straight line in the 25 to 55 °C region, and this gives a value of 48.6 kcal/mol for the activation energy of the denaturation reaction. The rate of this reaction is found to be very low at 25 °C; its half-life is about eleven days. Infrared absorption spectra observed in the amide I region suggest that the very slow denaturation of this protein is accompanied by a conformation change from an α-helix to a β-form. The number of the peptide groups involved in this α → β change is estimated to be 9 ± 3.     

Distinct effects of the C-terminal octapeptide of cholecystokinin and of a cholera toxin pretreatment on the kinetics of rat pancreatic adenylate cyclase activity

(1) The kinetic parameters of rat pancreatic adenylate cyclase were evaluated, using GTP, p[NH]ppG or GTPγS as nucleotide activator, cholecystokinin as peptide hormone, and GDPβS and dibutyryl cyclic GMP as inhibitors of guanosine triphosphate and CCK-8, respectively. The time courses of activation and the degree of activation at steady state (E_A/E_TOT) were compatible with a simple two-state model of activation-deactivation based on a pseudo-monomolecular activation process (rate constant k₊₂, and a deactivation process (rate constant k_off) that included, depending on the activating nucleotide, the hydrolysis of GTP (rate constant k₂) and/or the dissociation of the intact nucleotide (rate constant k_?1), so that E_A/E_TOT = k₊₁/(k₊₁ + k₂ + k_?a). (2) The hormone CCK-8 increased the value of k₊₁ with GTP dose-dependently, from 0.2 to 10.9 min^?1. The value of k_?1 increased 0.01 to 0.3 min^?1 but the value of k₂ was unaltered at 7 min^?1, so that E_A/E_TOT increased 15-fold, from 4% to 61%. (3) A cholera toxin pretreatment at 30 μg/ml allowed also a large increase in E_A/E_TOT with GTP (up to 51%) but the underlying mechanism was different. It consisted of a 14-fold decrease in the k_off value of the GTP-activated enzyme (from 7 min^? to 0.5 min^?1) that corresponded to a reduction in GTPase activity. When testing the system with p[NH]ppG, two added effects of the cholera toxin pretreatment were observed: a 4-fold increase in the value of k₊₁ (from 0.2 to 0.8 min^?1) and the occurrence of a significant 0.3 min^?1 value for k_?1.     

Evidence for a new early acceptor in Photosystem I of plants. An ESR investigation of reaction center triplet yield and of the reduced intermediary acceptors

The yield of the triplet state of the primary electron donor of Photosystem I of photosynthesis (P^T-700) and the characteristic parameters (g value, line shape, saturation behavior) of the ESR signal of the photoaccumulated intermediary acceptor A have been measured for two types of Photosystem I subchloroplast particles: Triton particles (TSF 1, about 100 chlorophyll molecules per P-700) that contain the iron-sulfur acceptors F_X, F_B and F_A, and lithium dodecyl sulfate (LDS) particles (about 40 chlorophyll molecules per P-700) that lack these iron-sulfur acceptors. The results are: (i) In Triton particles the yield of P^T-700 upon illumination is independent of the redox state of A and of F_X,B,A and is maximally about 5% of the active reaction centers at 5 K. The molecular sublevel decay rates are k_x = 1100 s^?1 ± 10%, k_y = 1300 s^?1 ± 10% and k_z = 83 s^?1 ± 20%. In LDS particles the triplet yield decreases linearly with concentration of reduced intermediary acceptors, the maximal yield being about 4% at 5 K assuming full P-700 activity. (ii) In Triton particles the acceptor complex A consists of two acceptors A₀ and A₁, with A₀ preceding A₁. In LDS particles at temperatures below ?30°C only A₀ is photoactive. (iii) The spin-polarized ESR signal found in the time-resolved ESR experiments with Triton particles is attributed to a polarized P-700-A^?₁ spectrum. The decay kinetics are complex and are influenced by transient nutation effects, even at low microwave power. It is concluded that the lifetime at 5 K of P-700A₀A^?₁ must exceed 5 ms. We conclude that P^T-700 originates from charge recombination of P-700A^?₀, and that in Triton particles A₀ and A₁ are both photoaccumulated upon cooling at low redox potential in the light. Since the state P-700AF^?_X does not give rise to triplet formation the 5% triplet yield in Triton particles is probably due to centers with damaged electron transport.     

A Population Based Study of the Genetic Association between Catecholamine Gene Variants and Spontaneous Low-Frequency Fluctuations in Reaction Time

The catecholamines dopamine and noradrenaline have been implicated in spontaneous low-frequency fluctuations in reaction time, which are associated with attention deficit hyperactivity disorder (ADHD) and subclinical attentional problems. The molecular genetic substrates of these behavioral phenotypes, which reflect frequency ranges of intrinsic neuronal oscillations (Slow-4: 0.027-0.073 Hz; Slow-5: 0.010-0.027 Hz), have not yet been investigated. In this study, we performed regression analyses with an additive model to examine associations between low-frequency fluctuations in reaction time during a sustained attention task and genetic markers across 23 autosomal catecholamine genes in a large young adult population cohort (n = 964), which yielded greater than 80% power to detect a small effect size (f² = 0.02) and 100% power to detect a small/medium effect size (f² = 0.15). At significance levels corrected for multiple comparisons, none of the gene variants were associated with the magnitude of low-frequency fluctuations. Given the study’s strong statistical power and dense coverage of the catecholamine genes, this either indicates that associations between low-frequency fluctuation measures and catecholamine gene variants are absent or that they are of very small effect size. Nominally significant associations were observed between variations in the alpha-2A adrenergic receptor gene (ADRA2A) and the Slow-5 band. This is in line with previous reports of an association between ADRA2A gene variants and general reaction time variability during response selection tasks, but the specific association of these gene variants and low-frequency fluctuations requires further confirmation. Pharmacological challenge studies could in the future provide convergent evidence for the noradrenergic modulation of both general and time sensitive measures of intra-individual variability in reaction time.