Random locomotion and chemotaxis of human blood polymorphonuclear leukocytes (PMN) in the presence of EDTA: PMN in close quarters require neither leukocyte integrins nor external divalent cations

Divalent cations are thought essential for motile function of leukocytes in general, and for the function of critical adhesion molecules in particular. In the current study, under direct microscopic observation with concomitant time-lapse video recording, we examined the effects of 10 mM EDTA on locomotion of human blood polymorphonuclear leukocytes (PMN). In very thin slide preparations, EDTA did not impair either random locomotion or chemotaxis; motile behavior appeared to benefit from the close approximation of slide and coverslip (“chimneying”). In preparations twice as thick, PMN in EDTA first exhibited active deformability with little or no displacement, then rounded up and became motionless. However, on creation of a chemotactic gradient, the same cells were able to orient and make their way to the target, often, however, losing momentarily their purchase on the substrate. In either of these preparations without EDTA, specific antibodies to β2 integrins did not prevent random locomotion or chemotaxis, even when we added antibodies to β1 and αvβ3 integrins and to integrin-associated protein, and none of these antibodies added anything to the effects of EDTA. In the more turbulent environment of even more media, effects of anti-β2 integrins became evident: PMN still could locomote but adhered to substrate largely by their uropods and by uropod-associated filaments. We relate these findings to the reported independence from integrins of PMN in certain experimental and disease states. Moreover, we suggest that PMN locomotion in close quarters is not only integrin-independent, but independent of external divalent cations as well.

A ribozyme-mediated, gene “knockdown” strategy for the identification of gene function in zebrafish

The zebrafish system offers many unique opportunities for the study of molecular biology. To date, only random mutagenesis, and not directed gene knockouts, have been demonstrated in this system. To more fully develop the potential of the zebrafish system, an approach to effectively inhibit the expression of any targeted gene in the developing zebrafish embryo has been developed. This approach uses a transient, cytoplasmic, T7 expression system, injected into the fertilized zebrafish egg to rapidly produce high levels of a ribozyme directed against the mRNA encoded by the targeted gene to inhibit its expression. In a demonstration of this strategy, expression of the recessive dominant zebrafish no tail gene was effectively inhibited by using this strategy to yield a phenotype identical to that resulting from a known defective mutation in this same gene. This, ribozyme-mediated, message deletion strategy may have use in determining the function of genetic coding sequences of unknown function.

Identification of two short internal ribosome entry sites selected from libraries of random oligonucleotides

Sequences that control translation of mRNA may play critical roles in regulating protein levels. One such element is the internal ribosome entry site (IRES). We previously showed that a 9-nt segment in the 5′ leader sequence of the mRNA encoding Gtx homeodomain protein could function as an IRES. To identify other short sequences with similar properties, we designed a selection procedure that uses a retroviral vector to express dicistronic mRNAs encoding enhanced green and cyan fluorescent proteins as the first and second cistrons, respectively. Expression of the second cistron was dependent upon the intercistronic sequences and was indicative of IRES activity. B104 cells were infected with two retroviral libraries that contained random sequences of 9 or 18 nt in the intercistronic region. Cells expressing both cistrons were sorted, and sequences recovered from selected cells were reassayed for IRES activity in a dual luciferase dicistronic mRNA. Two novel IRESes were identified by this procedure, and both contained segments with complementarity to 18S rRNA. When multiple copies of either segment were linked together, IRES activities were dramatically enhanced. Moreover, these synthetic IRESes were differentially active in various cell types. These properties are similar to those of the previously identified 9-nt IRES module from Gtx mRNA. These results provide further evidence that short nucleotide sequences can function as IRESes and support the idea that some cellular IRESes may be composed of shorter functional modules. The ability to identify IRES modules with specific expression properties may be useful in the design of vectors for biotechnology and gene therapy.

An essential amino acid residue in the protein translocation channel revealed by targeted random mutagenesis of SecY

The SecY/Sec61α family of membrane proteins are the central subunits of the putative protein translocation channel. We introduced random mutations into a segment of Escherichia coli SecY within its cytoplasmic domain 5, which was shown previously to be important for the SecA-dependent translocation activity. Mutations were classified into those retaining function and those gaining a dominant-interfering ability caused by a loss of function. These analyses showed that Arg-357, Pro-358, Gly-359, and Thr-362 are functionally important; Arg-357, conserved in almost all organisms, was identified as an indispensable residue.

Advances in random matrix theory, zeta functions, and sphere packing

Over four hundred years ago, Sir Walter Raleigh asked his mathematical assistant to find formulas for the number of cannonballs in regularly stacked piles. These investigations aroused the curiosity of the astronomer Johannes Kepler and led to a problem that has gone centuries without a solution: why is the familiar cannonball stack the most efficient arrangement possible? Here we discuss the solution that Hales found in 1998. Almost every part of the 282-page proof relies on long computer verifications. Random matrix theory was developed by physicists to describe the spectra of complex nuclei. In particular, the statistical fluctuations of the eigenvalues (“the energy levels”) follow certain universal laws based on symmetry types. We describe these and then discuss the remarkable appearance of these laws for zeros of the Riemann zeta function (which is the generating function for prime numbers and is the last special function from the last century that is not understood today.) Explaining this phenomenon is a central problem. These topics are distinct, so we present them separately with their own introductory remarks.

Random mutagenesis of Thermus aquaticus DNA polymerase I: concordance of immutable sites in vivo with the crystal structure.

Expression of Thermus aquaticus (Taq) DNA polymerase I (pol I) in Escherichia, coli complements the growth defect caused by a temperature-sensitive mutation in the host pol I. We replaced the nucleotide sequence encoding amino acids 659-671 of the O-helix of Taq DNA pol I, corresponding to the substrate binding site, with an oligonucleotide containing random nucleotides. Functional Taq pol I mutants were selected based on colony formation at the nonpermissive temperature. By using a library with 9% random substitutions at each of 39 positions, we identified 61 active Taq pol I mutants, each of which contained from one to four amino acid substitutions. Some amino acids, such as alanine-661 and threonine-664, were tolerant of several or even many diverse replacements. In contrast, no replacements or only conservative replacements were identified at arginine-659, lysine-663, and tyrosine-671. By using a library with totally random nucleotides at five different codons (arginine-659, arginine-660, lysine-663, phenylalanine-667, and glycine-668), we confirmed that arginine-659 and lysine-663 were immutable, and observed that only tyrosine substituted for phenylalanine-667. The two immutable residues and the two residues that tolerate only highly conservative replacements lie on the side of O-helix facing the incoming deoxynucleoside triphosphate, as determined by x-ray analysis. Thus, we offer a new approach to assess concordance of the active conformation of an enzyme, as interpreted from the crystal structure, with the active conformation inferred from in vivo function.

Novel human DNA alkyltransferases obtained by random substitution and genetic selection in bacteria.

DNA repair alkyltransferases protect organisms against the cytotoxic, mutagenic, and carcinogenic effects of alkylating agents by transferring alkyl adducts from DNA to an active cysteine on the protein, thereby restoring the native DNA structure. We used random sequence substitutions to gain structure-function information about the human O6-methylguanine-DNA methyltransferase (EC 2.1.1.63), as well as to create active mutants. Twelve codons surrounding but not including the active cysteine were replaced by a random nucleotide sequence, and the resulting random library was selected for the ability to provide alkyltransferase-deficient Escherichia coli with resistance to the methylating agent N-methyl-N'-nitro-N-nitrosoguanidine. Few amino acid changes were tolerated in this evolutionarily conserved region of the protein. One mutation, a valine to phenylalanine change at codon 139 (V139F), was found in 70% of the selected mutants; in fact, this mutant was selected much more frequently than the wild type. V139F provided alkyltransferase-deficient bacteria with greater protection than the wild-type protein against both the cytotoxic and mutagenic effects of N-methyl-N'-nitro-N-nitrosoguanidine, increasing the D37 over 4-fold and reducing the mutagenesis rate 2.7-5.5-fold. This mutant human alkyltransferase, or others similarly created and selected, could be used to protect bone marrow cells from the cytotoxic side effects of alkylation-based chemotherapeutic regimens.

Active barnase variants with completely random hydrophobic cores.

The central structural feature of natural proteins is a tightly packed and highly ordered hydrophobic core. If some measure of exquisite, native-like core packing is necessary for enzymatic function, this would constitute a significant obstacle to the development of novel enzymes, either by design or by natural or experimental evolution. To test the minimum requirements for a core to provide sufficient structural integrity for enzymatic activity, we have produced mutants of the ribonuclease barnase in which 12 of the 13 core residues have together been randomly replaced by hydrophobic alternatives. Using a sensitive biological screen, we find that a strikingly high proportion of these mutants (23%) retain enzymatic activity in vivo. Further substitution at the 13th core position shows that a similar proportion of completely random hydrophobic cores supports enzyme function. Of the active mutants produced, several have no wild-type core residues. These results imply that hydrophobicity is nearly a sufficient criterion for the construction of a functional core and, in conjunction with previous studies, that refinement of a crudely functional core entails more stringent sequence constraints than does the initial attainment of crude core function. Since attainment of crude function is the critical initial step in evolutionary innovation, the relatively scant requirements contributed by the hydrophobic core would greatly reduce the initial hurdle on the evolutionary pathway to novel enzymes. Similarly, experimental development of novel functional proteins might be simplified by limiting core design to mere specification of hydrophobicity and using iterative mutation-selection to optimize core structure.

On ANOVA Decompositions of Kernels and Gaussian Random Field Paths

The FANOVA (or “Sobol’-Hoeffding”) decomposition of multivariate functions has been used for high-dimensional model representation and global sensitivity analysis. When the objective function f has no simple analytic form and is costly to evaluate, computing FANOVA terms may be unaffordable due to numerical integration costs. Several approximate approaches relying on Gaussian random field (GRF) models have been proposed to alleviate these costs, where f is substituted by a (kriging) predictor or by conditional simulations. Here we focus on FANOVA decompositions of GRF sample paths, and we notably introduce an associated kernel decomposition into 4 d 4d terms called KANOVA. An interpretation in terms of tensor product projections is obtained, and it is shown that projected kernels control both the sparsity of GRF sample paths and the dependence structure between FANOVA effects. Applications on simulated data show the relevance of the approach for designing new classes of covariance kernels dedicated to high-dimensional kriging.

Evaluation of random plasma cortisol and the low dose corticotropin test as indicators of adrenal secretory capacity in critically ill patients: A prospective study

It is unclear whether a random plasma cortisol measurement and the corticotropin (ACTH) test adequately reflect glucocorticoid secretory capacity in critical illness. This study aimed to determine whether these tests provide information representative of the 24 hour period. Plasma cortisol was measured hourly for 24 hours in 21 critically ill septic patients followed by a corticotropin test with 1 μ g dose administered intravenously. Serum and urine were analysed for ACTH and free cortisol respectively. Marked hourly variability in plasma cortisol was evident (coefficient of variation 8-30%) with no demonstrable circadian rhythm. The individual mean plasma cortisol concentrations ranged from 286 59 nmol/l to 796 &PLUSMN; 83 nmol/l. The 24 hour mean plasma cortisol was strongly correlated with both random plasma cortisol (r(2) 0.9, P< 0.0001) and the cortisol response to corticotropin (r(2) 0.72, P< 0.001). Only nine percent of patients increased their plasma cortisol by 250 nmol/l after corticotropin (euadrenal response). However, 35% of non-responders had spontaneous hourly rises > 250 nmol/l thus highlighting the limitations of a single point corticotropin test. Urinary free cortisol was elevated (865&PLUSMN; 937 nmol) in both corticotropin responders and non-responders suggesting elevated plasma free cortisol. No significant relationship was demonstrable between plasma cortisol and ACTH. We conclude that although random cortisol measurements and the low dose corticotropin tests reliably reflect the 24 hour mean cortisol in critical illness, they do not take into account the pulsatile nature of cortisol secretion. Consequently, there is the potential for erroneous conclusions about adrenal function based on a single measurement. We suggest that caution be exercised when drawing conclusions on the adequacy of adrenal function based on a single random plasma cortisol or the corticotropin test.

A Mixture model with random-effects components for clustering correlated gene-expression profiles

Motivation: The clustering of gene profiles across some experimental conditions of interest contributes significantly to the elucidation of unknown gene function, the validation of gene discoveries and the interpretation of biological processes. However, this clustering problem is not straightforward as the profiles of the genes are not all independently distributed and the expression levels may have been obtained from an experimental design involving replicated arrays. Ignoring the dependence between the gene profiles and the structure of the replicated data can result in important sources of variability in the experiments being overlooked in the analysis, with the consequent possibility of misleading inferences being made. We propose a random-effects model that provides a unified approach to the clustering of genes with correlated expression levels measured in a wide variety of experimental situations. Our model is an extension of the normal mixture model to account for the correlations between the gene profiles and to enable covariate information to be incorporated into the clustering process. Hence the model is applicable to longitudinal studies with or without replication, for example, time-course experiments by using time as a covariate, and to cross-sectional experiments by using categorical covariates to represent the different experimental classes. Results: We show that our random-effects model can be fitted by maximum likelihood via the EM algorithm for which the E(expectation) and M(maximization) steps can be implemented in closed form. Hence our model can be fitted deterministically without the need for time-consuming Monte Carlo approximations. The effectiveness of our model-based procedure for the clustering of correlated gene profiles is demonstrated on three real datasets, representing typical microarray experimental designs, covering time-course, repeated-measurement and cross-sectional data. In these examples, relevant clusters of the genes are obtained, which are supported by existing gene-function annotation. A synthetic dataset is considered too.

Random vibration of the functionally graded laminates in thermal environments

This work deals with the random free vibration of functionally graded laminates with general boundary conditions and subjected to a temperature change, taking into account the randomness in a number of independent input variables such as Young's modulus, Poisson's ratio and thermal expansion coefficient of each constituent material. Based on third-order shear deformation theory, the mixed-type formulation and a semi-analytical approach are employed to derive the standard eigenvalue problem in terms of deflection, mid-plane rotations and stress function. A mean-centered first-order perturbation technique is adopted to obtain the second-order statistics of vibration frequencies. A detailed parametric study is conducted, and extensive numerical results are presented in both tabular and graphical forms for laminated plates that contain functionally graded material which is made of aluminum and zirconia, showing the effects of scattering in thermo-clastic material constants, temperature change, edge support condition, side-to-thickness ratio, and plate aspect ratio on the stochastic characteristics of natural frequencies. (c) 2005 Elsevier B.V. All rights reserved.

Emulation of random output simulators

Computer models, or simulators, are widely used in a range of scientific fields to aid understanding of the processes involved and make predictions. Such simulators are often computationally demanding and are thus not amenable to statistical analysis. Emulators provide a statistical approximation, or surrogate, for the simulators accounting for the additional approximation uncertainty. This thesis develops a novel sequential screening method to reduce the set of simulator variables considered during emulation. This screening method is shown to require fewer simulator evaluations than existing approaches. Utilising the lower dimensional active variable set simplifies subsequent emulation analysis. For random output, or stochastic, simulators the output dispersion, and thus variance, is typically a function of the inputs. This work extends the emulator framework to account for such heteroscedasticity by constructing two new heteroscedastic Gaussian process representations and proposes an experimental design technique to optimally learn the model parameters. The design criterion is an extension of Fisher information to heteroscedastic variance models. Replicated observations are efficiently handled in both the design and model inference stages. Through a series of simulation experiments on both synthetic and real world simulators, the emulators inferred on optimal designs with replicated observations are shown to outperform equivalent models inferred on space-filling replicate-free designs in terms of both model parameter uncertainty and predictive variance.

A randomised controlled trial investigating the effect of lutein and antioxidant dietary supplementation on visual function in healthy eyes

Background & aims It has been suggested that retinal lutein may improve visual acuity for images that are illuminated by white light. Our aim was to determine the effect of a lutein and antioxidant dietary supplement on visual function. Methods A prospective, 9- and 18-month, double-masked randomised controlled trial. For the 9-month trial, 46 healthy participants were randomised (using a random number generator) to placebo (n=25) or active (n=21) groups. Twenty-nine of these subjects went on to complete 18 months of supplementation, 15 from the placebo group, and 14 from the active group. The active group supplemented daily with 6mg lutein combined with vitamins and minerals. Outcome measures were distance and near visual acuity, contrast sensitivity, and photostress recovery time. The study had 80% power at the 5% significance level for each outcome measure. Data were collected at baseline, 9, and 18 months. Results There were no statistically significant differences between groups for any of the outcome measures over 9 or 18 months. Conclusion There was no evidence of effect of 9 or 18 months of daily supplementation with a lutein-based nutritional supplement on visual function in this group of people with healthy eyes. ISRCTN78467674.

Noisy random Boolean formulae:a statistical physics perspective

Properties of computing Boolean circuits composed of noisy logical gates are studied using the statistical physics methodology. A formula-growth model that gives rise to random Boolean functions is mapped onto a spin system, which facilitates the study of their typical behavior in the presence of noise. Bounds on their performance, derived in the information theory literature for specific gates, are straightforwardly retrieved, generalized and identified as the corresponding macroscopic phase transitions. The framework is employed for deriving results on error-rates at various function-depths and function sensitivity, and their dependence on the gate-type and noise model used. These are difficult to obtain via the traditional methods used in this field.