Brush effects on DNA chips: thermodynamics, kinetics, and design guidelines

In biology experiments, oligonucleotide microarrays are contacted with a solution of long nucleic acid targets. The hybridized probes thus carry long tails. When the surface density of the oligonucleotide probes is high enough, the progress of hybridization gives rise to a polyelectrolyte brush due to mutual crowding of the nucleic acid tails. The free-energy penalty associated with the brush modifies both the hybridization isotherms and the rate equations: the attainable hybridization is lowered significantly as is the hybridization rate. When the equilibrium hybridization fraction, x(eq), is low, the hybridization follows a Langmuir type isotherm, x(eq)/(1 - x(eq)) = c(t)K where c(t) is the target concentration and K is the equilibrium constant. K is smaller than its bulk value by a factor (n/N)(2/5) due to wall effects where n and N denote the number of bases in the probe and the target. At higher x(eq), when the brush is formed, the leading correction is x(eq)/(1 - x(eq)) = c(t)K exp - const'x(eq)(2/3) - x(B)(2/3) where x(B) corresponds to the onset of the brush regime. The denaturation rate constant in the two regimes is identical. However, the hybridization rate constant in the brush regime is lower, the leading correction being exp -const' x(2/3) - x(B)(2/3).     

Comparative Analysis of <Emphasis Type="Italic">Period</Emphasis> Genes in Teleost Fish Genomes

Period (Per) is a canonical circadian clock gene. The fruit fly, an invertebrate, has one per gene, while the human, a tetrapod vertebrate, has three Per genes. Per1, Per2, and Per3 of the tetrapods were generated from two rounds of ancient genome duplications from the ancestral chordate Per gene. Searching for five teleost fish genomes in a combination of phylogenetic, splicing site, and syntenic analyses revealed that zebrafish have two per1 genes, per1a and per1b, one per2, and one per3; medaka, fugu, and tetraodon each have two per2 genes, per2a and per2b, one per1, and one per3; sticklebacks also have per2a, per2b, and one per1 but lack per3; and per1a/per1b in zebrafish and per2a/per2b in madaka, fugu, tetraodon, and stickleback are ancient duplicates. While the dN/dS ratios of the five fish per duplicates are all <1, suggesting that they likely have been subject to purifying selection, the Tajima relative rate test showed that zebrafish per1a/per1b and fugu and medaka per2a/per2b have asymmetric evolutionary rates, implicating that one of these duplicates might have been under positive selection or relaxed functional constraint. Further, in situ hybridization showed that zebrafish per1a and per1b clearly have distinct patterns of temporal and spatial expression. These results support the notion that extra copies of teleost per genes were generated from the fish-specific genome duplication, and divergent resolution after the duplication resulted in retention of different per duplicates in different fish, most of which have diverged significantly.     

Effects of geometric isomerism in dinuclear platinum antitumor complexes on aquation reactions in the presence of perchlorate, acetate and phosphate

The aquation and subsequent reactions of the dinuclear Pt antitumor complexes [{trans-PtCl(NH(3))(2)}(2)(mu-NH(2)(CH(2))(6)NH(2))](2+) (1,1/t,t) and [{cis-PtCl(NH(3))(2)}(2)(mu-NH(2)(CH(2))(6)NH(2))](2+) (1,1/c,c) in 15 mM perchlorate, acetate or phosphate solutions were followed at 298 K by [(1)H,(15)N] HSQC 2D NMR spectroscopy. Rate and equilibrium constants for the initial reversible aquation and the subsequent reversible reaction with phosphate or acetate are reported. The rate constant for the first aquation step is two-fold lower for 1,1/c,c than 1,1/t,t but the anation rate constants are similar so that the equilibrium lies further towards the chloro form for the 1,1/c,c compound. A pK (a) value of 6.01+/-0.03 was determined for the diaquated species [{cis-Pt(NH(3))(2)(H(2)O)}(2)(mu-NH(2)(CH(2))(6)NH(2))](4+) (1,1/c,c-3) which is 0.4 units higher than that of the 1,1/t,t compound. The rate constants for the binding of acetate and phosphate to 1,1/t,t are similar, but the rate constant for the reverse reaction is close to ten-fold higher in the case of phosphate so that equilibrium conditions are attained more rapidly (12 h compared with 64 h). On the other hand, for 1,1/c,c the rate constants for the forward and reverse reactions with acetate and phosphate are quite similar so that equilibrium conditions are reached very slowly (80-100 h) and a greater proportion of phosphate-bound species are present. The reduced lability of the bound phosphate for 1,1/c,c is attributed to the formation of a macrochelate phosphate-bridged species which was characterized by (31)P NMR and ESI-MS. The speciation profiles of 1,1/t,t and 1,1/c,c under physiological conditions are explored.     

Selection and self-organization of self-reproducing macromolecules under the constraint of constant flux.

We investigate the dynamic behavior of a set of self-reproducing macromolecules (e.g., polynucleotides) under conditions such that the fluxes of all monomer units into the system are kept constant. Such conditions might prevail in an evolution reactor or in certain naturally occurring situations. A general set of equations is developed to describe the behavior of both the macromolecule and the monomer concentrations. The question of how the rate of macromolecule synthesis varies with the monomer levels is discussed briefly. With the help of several physically reasonable approximations, we obtain an exact solution for a simplified constant flux system. Comparison with the corresponding system under the constraint of constant overall organization reveals important similarities, most notably in the existance and composition of quasispecies. Given the same set of physical and chemical parameters, a system subject to constant flux will always evolve toward selective equilibrium more slowly than under the constraint of constant organization.     

Equilibrium studies on the binding of cadmium(II) to human serum transferrin

The binding of cadmium(II) to human serum transferrin in 0.01 M N-(2-hydroxyethyl)-piperazine-N'-2-ethanesulfonic acid with 5 mM bicarbonate at 25 degrees C has been evaluated by difference ultraviolet spectroscopy. Equilibrium constants were determined by competition versus three different low molecular weight chelating agents: nitrilotriacetic acid, ethylenediamine-N,N'-diacetic acid, and triethylenetetramine. Conditional equilibrium constants for the sequential binding of two cadmium ions to transferrin under the stated experimental conditions are log K1 = 5.95 +/- 0.10 and log K2 = 4.86 +/- 0.13. A linear free energy relationship for the complexation of cadmium and zinc has been prepared by using equilibrium data on 243 complexes of these metal ions with low molecular weight ligands. The transferrin binding constants for cadmium and zinc are in good agreement with this linear free energy relationship. This indicates that the larger size of the cadmium(II) ion does not significantly hinder its binding to the protein.     

Molecular cytogenetics of Leymus: Mapping the Ns genome-specific repetitive sequences

The objective of this paper is to summarize the work in my group on FISH (fluorescent in situ hybridization) mapping of Ns-specific repetitive DNA sequences fromLeymus and discuss the results in the context of classification based on the genome system currently used among Triticeae researchers. The key question here is whether the genome composition of a tetraploid Leymus species should be NsXm or NsNs (Ns₁Ns₂). Different types of Leymus-specific dispersed retroelement-like repeats have been isolated and characterized. Because the sequences occur in significantly high copy number in Leymus, based on strong hybridization signal in Southern blots, they are considered essentially specific to Leymus. They are also abundant in Psathyrostachys, the progenitor of Ns genome in Leymus. These dispersed repeats are found to distribute over the whole of all Leymus chromosomes, without any differentiation between chromosomes that have been suggested to be of different genomic origins, meaning that all genomes in Leymus are the same. GISH (genomic in situ hybridization) experiments on Leymus chromosomes using Psathyrostachys genomic DNA as probes further support the NsNs (Ns₁Ns₂) genome constitution for Leymus. The Xm genome of an unknown origin might have been there in the beginning of the allopolyploidization process, but the Ns genome-specific elements must have spread predominantly and rapidly across genomes, thus homogenizing the nuclear genomes of Leymus. I present here for the first time evidence that Ns-specific dispersed repeats can spread in a very short time, from Leymus over to wheat in Triticum × Leymus hybrids growing in artificial conditions.     

Acyl migration in 1,2-dipalmitoyl-sn-glycerol

The acyl migration of 1,2-dipalmitoyl-sn-glycerol (1,2-DPG) to 1,3-dipalmitoylglycerol (1,3-DPG) in different states, neat, in the presence of egg yolk lecithin (sonicated and unsonicated) and on silica gel was studied. The isomerization was quantitated by scanning densitometry of charred TLC plates, at different temperatures and for varying periods of time. At equilibrium the amount of 1,3-DPG was found to be 56%. The rates of initial isomerization, and the time required to isomerize to half the equilibrium quantity (i.e., t1/2 eq. = 1,3-DPG 28%) under the above conditions was estimated. In the case of neat melt at 74 degrees C and in an organic solvent the time required to t1/2 eq. is 18 h and a few days, respectively. However, at 62 degrees C in the presence of a polar solvent (sodium phosphate buffer at pH 7.0) the t1/2 eq. is 1-2 h. On dry silica gel (TLC plate) at 24 degrees C the t1/2 eq. is reached in less than 1 h.     

The nonclassical MHC class I molecule Qa-1 forms unstable peptide complexes

The MHC class Ib molecule Qa-1 is the primary ligand for mouse CD94/NKG2A inhibitory receptors expressed on NK cells, in addition to presenting Ags to a subpopulation of T cells. CD94/NKG2A receptors specifically recognize Qa-1 bound to the MHC class Ia leader sequence-derived peptide Qdm. Qdm is the dominant peptide loaded onto Qa-1 under physiological conditions and this peptide has an optimal sequence for binding to Qa-1. Peptide dissociation experiments demonstrated that Qdm dissociates from soluble or cell surface Qa-1(b) molecules with a t(1/2) of approximately 1.5 h at 37 degrees C. In comparison, complexes of an optimal peptide (SIINFEKL) bound to the MHC class Ia molecule H-2K(b) dissociated with a t(1/2) in the range from 11 to 31 h. In contrast to K(b), the stability of cell surface Qa-1(b) molecules was independent of bound peptides, and several observations suggested that empty cell surface Qa-1(b) molecules might be unusually stable. Consistent with the rapid dissociation rate of Qdm from Qa-1(b), cells become susceptible to lysis by CD94/NKG2A(+) NK cells under conditions in which new Qa-1(b)/Qdm complexes cannot be continuously generated at the cell surface. These results support the hypothesis that Qa-1 has been selected as a specialized MHC molecule that is unable to form highly stable peptide complexes. We propose that the CD94/NKG2A-Qa-1/Qdm recognition system has evolved as a rapid sensor of the integrity of the MHC class I biosynthesis and Ag presentation pathway.     

The importance of thermodynamic equilibrium for high throughput gene expression arrays

We present an analysis of physical chemical constraints on the accuracy of DNA micro-arrays under equilibrium and nonequilibrium conditions. At the beginning of the article we describe an algorithm for choosing a probe set with high specificity for targeted genes under equilibrium conditions. The algorithm as well as existing methods is used to select probes from the full Saccharomyces cerevisiae genome, and these probe sets, along with a randomly selected set, are used to simulate array experiments and identify sources of error. Inasmuch as specificity and sensitivity are maximum at thermodynamic equilibrium, we are particularly interested in the factors that affect the approach to equilibrium. These are analyzed later in the article, where we develop and apply a rapidly executable method to simulate the kinetics of hybridization on a solid phase support. Although the difference between solution phase and solid phase hybridization is of little consequence for specificity and sensitivity when equilibrium is achieved, the kinetics of hybridization has a pronounced effect on both. We first use the model to estimate the effects of diffusion, crosshybridization, relaxation time, and target concentration on the hybridization kinetics, and then investigate the effects of the most important kinetic parameters on specificity. We find even when using probe sets that have high specificity at equilibrium that substantial crosshybridization is present under nonequilibrium conditions. Although those complexes that differ from perfect complementarity by more than a single base do not contribute to sources of error at equilibrium, they slow the approach to equilibrium dramatically and confound interpretation of the data when they dissociate on a time scale comparable to the time of the experiment. For the best probe set, our simulation shows that steady-state behavior is obtained in a relaxation time of approximately 12-15 h for experimental target concentrations approximately (10(-13) - 10(-14))M, but the time is greater for lower target concentrations in the range (10(-15)-10(-16))M. The result points to an asymmetry in the accuracy with which up- and downregulated genes are identified.     

High reproducibility using sodium hydroxide-stripped long oligonucleotide DNA microarrays

Recently, long oligonucleotide (60- to 70-mer) microarrays for two-color experiments have been developed and are gaining widespread use. In addition, when there is limited availability of mRNA from tissue sources, RNA amplification can and is being used to produce sufficient quantities of cRNA for microarray hybridization. Taking advantage of the selective degradation of RNA under alkaline conditions, we have developed a method to "strip" glass-based oligonucleotide microarrays that use fluorescent RNA in the hybridization, while leaving the DNA oligonucleotide probes intact and usable for a second experiment. Replicate microarray experiments conducted using stripped arrays showed high reproducibility, however, we found that arrays could only be stripped and reused once without compromising data quality. The intraclass correlation (ICC) between a virgin array and a stripped array hybridized with the same sample showed a range of 0.90-0.98, which is comparable to the ICC of two virgin arrays hybridized with the same sample. Using this method, once-stripped oligonucleotide microarrays are usable, reliable, and help to reduce costs.     

Graviresponse of cress roots under varying gravitational forces

During the spacelab mission IML-2 threshold values concerning gravity controlled growth processes have been estimated in order to test the reciprocity law (dose = stimulus x time = constant) for the first time under exact physiological conditions. Cress seedlings have been cultivated from dry seeds under conditions of microgravity and on a 1 x g-centrifuge in the ESA-BIORACK. With the help of NIZEMI--the slow rotating centrifuge microscope--these seedlings have been stimulated by different doses ranging from 12 to 60 x g x s. Two different values of acceleration--0.1 x g and 1 x g--have been used. Graviresponses of the roots have been documented by video recording for 60 min under conditions of microgravity. The response of roots to accelerations of 0.1 x g was remarkably less than to 1 x g in spite of the same doses being applied to the seedlings. Roots cultivated under conditions of microgravity showed a higher sensitivity than those grown on the 1 x g-centrifuge. Displacement of statoliths in gravity perceiving cells was mainly less than 1 micron under the different stimulation procedures. These results together with results from former space flights do not confirm the validity of the reciprocity law. They indicate that transformation of the gravistimulus has to occur in close vicinity to the statoliths, probably mediated by the ground cytoplasm and the cytoskeleton suspended therein.     

Binding of colchicine to purified microtubule protein.

The binding of colchicine to tubulin, purified by two cycles of assembly-disassembly, has been studied. Equilibrium studies indicated a dissociation constant which declined during incubation approaching a minimum value of approximately 0.30 times 10- minus 6 M after 13 hours of incubation. Because tubulin is unstable during prolonged incubation (t1/2 of 5.2 hours for free tubulin, t1/2 of 12.5 hours for tubulin bound to colchicine), the equilibrium Kd was felt to be an overestimation of the true Kd. The rate constant of dissociation (k-1 equal to 0.009 hour- minus 1 hour- minus 1) and the rate constant of association (k1 equal to 0.37 times 10-6 M-minus 1) were measured under conditions designed to circumvent or correct for tubulin instability. The dissociation constant determined by the ratio k-1/k1 was 0.024 times -minus 6 M. To determine whether the discrepancy between the "equilibrium" and "kinetic" determined dissociation constants could be accounted for on the basis of tubulin instability, the binding reaction was computer-simulated using the measured association and dissociation rate constants and the rate constants for decay of bound and free tubulin. Computer simulation was in close agreement with the experimentally determined behavior of the reaction during a 13-hour incubation. It is concluded that the Kd determined by equilibrium methodology results in a considerable overestimation due to the instability of tubulin, and that the best estimate for the Kd of the colchicine-tubulin equilibrium is the value determined by the ratio of the rate constants.     

Development of a fast DNA-DNA hybridization method based on melting profiles in microplates

DNA-DNA hybridization is still the “gold standard” for the genotypic delineation of bacterial species. However, it is not widely used because traditional DNA-DNA hybridization techniques are rather time-consuming and not easy to perform in routine laboratories. In the present study, DNA of reference strains was digested with Sau3A, ligated with linker oligonucleotides S1/2 and in vitro amplified. The amplified DNA fragments were immobilized on MaxiSorb 96-well plates. DNA isolated from target strains was also digested with Sau3A, ligated with linker oligonuleotides P1/2 and in vitro amplified in the presence of digoxygenin modified dUTP. The labeled amplificate was hybridized to the immobilized reference DNA under isothermal conditions. Thermal denaturation curves of the DNA-DNA hybrids were obtained by using washing solutions of increasing stringency. Remaining hybrids were colorimetrically detected with anti-digoxygenin-horseradish peroxidase anti-bodies. The new method was validated with strains of the genus Pedioccocus for which DNA-DNA similarities have also been determined by the filter hybridization method. In addition, DNA-DNA hybridizations were performed with genotypically defined Enterobacter species.     

Kinetic study of the reduction mechanism for Desulfovibrio gigas cytochrome c3.

The kinetic aspects of the reduction process in cytochrome c3 from Desulfovibrio gigas have been investigated over a wide range of pH values ranging between pH 5.8 and pH 9.8. The data have been analyzed in the framework of an I2H4 interaction network coupled to a proton-linked equilibrium between two tertiary structures (Cornish-Bowden, A. & Koshland, D.E. Jr (1970) J. Biol. Chem. 245, 6241-6250). The kinetic rate constants for the reduction of the four hemes for the two tertiary conformations have been characterized in the framework of the thermodynamic network obtained from the equilibrium analysis (Coletta, M., Catarino, T., LeGall, J.J. & Xavier, A.V. (1991) Eur. J. Biochem. 202, 1101-1106). The intrinsic reduction rate constants determined by reaction with sodium dithionite for two hemes (namely heme 4 and heme 1) are significantly faster than those for the other two heme residues. In view of the equilibrium redox properties, heme 4 (with the fastest reduction rate) may then work as the kinetic electron-capturing site for the electrons from sodium dithionite. The transfer to hemes 2 and 3 then occurs by virtue of their free-energy levels at equilibrium. At our experimental conditions, there is also transfer of electrons to hemes 2 and 3 from heme 1, which is reduced at a slower rate than heme 4, thus contributing to the biphasic kinetics observed for the overall process. The kinetic parameters obtained are discussed in terms of the mechanism proposed for the coupling between the electron and proton transfer, as induced by the heme/heme cooperativity network.     

Introgression in natural populations of bioindicators: a case study of Carabus splendens and Carabus punctatoauratus

The evolutionary importance of hybridization in wild plants and animals has become increasingly widely recognized in the last decade. In practical terms, hybridization provides an exceptionally tough set of problems for conservation biologists. We illustrate this in a case study of two Carabidae species widely used to evaluate the impact of human activities on biodiversity. These two species live in a complex mosaic of sympatry/allopatry and are known to hybridize in controlled conditions. Hybridization has not been quantified in natural populations to date due to the lack of a simple set of phenotypic traits for identifying hybrids. We thus screened for hybrids in natural populations, by multilocus genotyping at nine microsatellite loci. A high level of genetic differentiation between these two taxa was observed, as shown by allelic frequency distributions. Two Bayesian assignment procedures without obligatory pure taxon references were used to infer different classes of hybrids (F(1), F(2) and backcrosses) and mixture proportions between the two species. A low level of hybridization (F(1) genotypes) was observed in natural populations, contrasting with results obtained in controlled conditions. A high level of introgression was, however, detected at three of 12 sites, as revealed by the detection of backcrossed genotypes. This interspecific gene flow was detected in a limited zone of the common geographical range of the two species and was not related to the pattern of sympatry/allopatry. We then considered the origin and repercussions of this introgression, based on intraspecific genetic diversity and geographical structure.     

Optimization of DNA blot hybridization conditions for various types of membranes]

A set of experiments has been conducted to choose the optimal conditions for DNA transfer and fixation on two types of the nitrocellulose, three types of nylon membranes and on capron filters. The buffer capillary transfer systems, electroblotting and gel hybridization are analyzed. Two techniques for DNA binding have been tested under different transfer conditions for all the membrane types: a vacuum fixation at 80 degrees C and a UV-exposure. The results indicate the critical dependence of the efficiency of blot-hybridization on the conditions of UV-treatment. The UV-exposure longer or shorter than the optimal one resulted in a loss of the hybridization efficiency. The optimal DNA-transfer and fixation conditions are recommended for all the membranes tested. The dependence of the optimal transfer and binding conditions on the specific characteristics of different membrane types was demonstrated for maximal sensitivity of the blot-hybridization.     

Assignment of the human gene for β-microseminoprotein (MSMB) to chromosome 10 and demonstration of related genes in other vertebrates

The gene for β-microseminoprotein MSMB has been studied by DNA hybridization and molecular cloning techniques. Comparative analysis of restriction endonuclease digests of the cloned gene and of leukocyte DNA strongly suggested that the gene is present in a single copy in the haploid human genome. By Southern blot analysis of DNA from somatic cell hybrids, the gene was assigned to chromosome 10. The coding nucleotides of the human gene are separated into four exons by relatively large introns. A related gene might be present in other mammals, birds, and amphibians as revealed by DNA hybridization under conditions of low stringency.     

DNA hybridization on microparticles: determining capture-probe density and equilibrium dissociation constants

Many DNA-probe assays utilize oligonucleotide-coated microparticles for capture of complementary nucleic acids from solution. During development of these assays, as well as in other particle-based nucleic acid applications, it is useful to know both the amount of duplex formation expected under various experimental conditions and the coating density of the capture oligonucleotide on the particle surface. We examined the simplest form of a DNA-probe microparticle assay: hybridization of a particle-bound capture oligonucleotide to its solution-phase complement. Fluorescein-labeled solution-phase oligonucleotide was hybridized to varying amounts of particles, and the amount of labeled oligonucleotide remaining in solution at equilibrium was measured. We present a simple two-state, all-or-none model for bimolecular hybridization of non-self-complementary sequences that can be used to calculate the equilibrium dissociation constant ( Kd ) from hybridization data. With experimental conditions where both the Kd value and the concentration of capture probe in the reaction are small relative to the concentration of labeled complementary oligonucleotide in the reaction, density of the capture probe on the particle's surface can also be determined. Kd values for particle-based hybridization were different from those obtained from solution-phase thermodynamic parameters. At higher temperatures, hybridization on particles was more efficient than hybridization in solution.     

Environment-Related Adaptive Changes of Gut Commensal Microbiota Do not Alter Colonic Toll-Like Receptors but Modulate the Local Expression of Sensory-Related Systems in Rats

Pathogenic and protective roles have been attributed to gut commensal microbiota (GCM) in gastrointestinal inflammatory and functional disorders. We have shown that the adaptation to a new environment implies specific changes in the composition of GCM. Here we assessed if environment-related adaptive changes of GCM modulate the expression of colonic Toll-like receptors (TLRs) and sensory-related systems in rats. Adult male SD rats were maintained under different environmental conditions: barrier-breed-and-maintained, barrier-breed adapted to conventional conditions or conventional-breed-and-maintained. Fluorescent in situ hybridization and real-time quantitative PCR (qPCR) were used to characterize luminal ceco-colonic microbiota. Colonic expression of TLR2, TLR4, TLR5, and TLR7, cannabinoid receptors (CB1/CB2), μ-opioid receptor (MOR), transient receptor potential vanilloid (TRPV1, TRPV3, and TRPV4), protease-activated receptor 2 (PAR-2), and calcitonin gene-related peptide were quantified by RT-qPCR. CB1, CB2 and MOR expression, was evaluated also by immunohistochemistry. In rats, housing-related environmental conditions induce specific changes of GCM, without impact on the expression of TLR-dependent bacterial recognition systems. Expression of sensory-related markers (MOR, TRPV3, PAR-2, and CB2) decreased with the adaptation to a conventional environment, correlating with changes in Bacteroides spp., Lactobacillus spp., and Bifidobacterium spp. counts. This suggests an interaction between GCM and visceral sensory mechanisms, which might be part of the mechanisms underlying the beneficial effects of some bacterial groups on functional and inflammatory gastrointestinal disorders.