Recent advances in microfluidic devices for single-cell cultivation: methods and applications

Single-cell analysis is essential to improve our understanding of cell functionality from cellular and subcellular aspects for diagnosis and therapy. Single-cell cultivation is one of the most important processes in single-cell analysis, which allows the monitoring of actual information of individual cells and provides sufficient single-cell clones and cell-derived products for further analysis. The microfluidic device is a fast-rising system that offers efficient, effective, and sensitive single-cell cultivation and real-time single-cell analysis conducted either on-chip or off-chip. Here, we introduce the importance of single-cell cultivation from the aspects of cellular and subcellular studies. We highlight the materials and structures utilized in microfluidic devices for single-cell cultivation. We further discuss biological applications utilizing single-cell cultivation-based microfluidics, such as cellular phenotyping, cell-cell interactions, and omics profiling. Finally, present limitations and future prospects of microfluidics for single-cell cultivation are also discussed.

A simple and reversible glass-glass bonding method to construct a microfluidic device and its application for cell recovery

Compared with polymer microfluidic devices, glass microfluidic devices have advantages for diverse lab-on-a-chip applications due to their rigidity, optical transparency, thermal stability, and chemical/biological inertness. However, the bonding process to construct glass microfluidic devices usually involves treatment(s) like high temperature over 400 °C, oxygen plasma or piranha solution. Such processes require special skill, apparatus or harsh chemicals, and destroy molecules or cells in microchannels. Here, we present a simple method for glass-glass bonding to easily form microchannels. This method consists of two steps: placing water droplets on a glass substrate cleaned by neutral detergent, followed by fixing a cover glass plate on the glass substrate by binding clips for a few hours at room temperature. Surface analyses showed that the glass surface cleaned by neutral detergent had a higher ratio of SiOH over SiO than glass surfaces prepared by other cleaning steps. Thus, the suggested method could achieve stronger glass-glass bonding via dehydration condensation due to the higher density of SiOH. The pressure endurance reached over 600 kPa within 6 h of bonding, which is sufficient for practical microfluidic applications. Moreover, by exploiting the reversibility of this bonding method, cell recoveries after cultivating cells in a microchannel were demonstrated. This new bonding method can significantly improve both the productivity and the usability of glass microfluidic devices and extend the possibility of glass microfluidic applications in future.

A sub-population of Dictyostelium discoideum cells shows extremely high sensitivity to cAMP for directional migration

The chemotaxis of Dictysotelium discoideum cells in response to a chemical gradient of cyclic adenosine 3',5'-monophosphate (cAMP) was studied using a newly designed microfluidic device. The device consists of 800 cell-sized channels in parallel, each 4 μm wide, 5 μm high, and 100 μm long, allowing us to prepare the same chemical gradient in all channels and observe the motility of 500-1000 individual cells simultaneously. The percentage of cells that exhibited directed migration was determined for various cAMP concentrations ranging from 0.1 pM to 10 μM. The results show that chemotaxis was highest at 100 nM cAMP, consistent with previous observations. At concentrations as low as 10 pM, about 16% of cells still exhibited chemotaxis, suggesting that the receptor occupancy of only 6 cAMP molecules/cell can induce chemotaxis in very sensitive cells. At 100 pM cAMP, chemotaxis was suppressed due to the self-production and secretion of intracellular cAMP induced by extracellular cAMP. Overall, systematic observations of a large number of individual cells under the same chemical gradients revealed the heterogeneity of chemotaxis responses in a genetically homogeneous cell population, especially the existence of a sub-population with extremely high sensitivity for chemotaxis.

The potential geographical distribution and phenology of Bemisia tabaci Middle East/Asia Minor 1, considering irrigation and glasshouse production

The Bemisia tabaci species complex is one of the most important pests of open field and protected cropping globally. Within this complex, one species (Middle East Asia Minor 1, B. tabaci MEAM1, formerly biotype B) has been especially problematic, invading widely and spreading a large variety of plant pathogens, and developing broad spectrum pesticide resistance. Here, we fit a CLIMEX model to the distribution records of B. tabaci MEAM1, using experimental observations to calibrate its temperature responses. In fitting the model, we consider the effects of irrigation and glasshouses in extending its potential range. The validated niche model estimates its potential distribution as being considerably broader than its present known distribution, especially in the Americas, Africa and Asia. The potential distribution of the fitted model encompasses the known distribution of B. tabaci sensu lato, highlighting the magnitude of the threat posed globally by this invasive pest species complex and the viruses it vectors to open field and protected agriculture.

Ultrasensitive Single Cell Metabolomics by Capillary Electrophoresis-Mass Spectrometry with a Thin-Walled Tapered Emitter and Large-Volume Dual Sample Preconcentration

Single cell metabolome analysis is essential for studying microscale life phenomena such as neuronal networks and tumor microenvironments. Capillary electrophoresis-mass spectrometry (CE-MS) is one of the most sensitive technologies; however, its sensitivity is still not enough for single cell analysis on general human cells such as HeLa. To address these issues, we first developed an efficient ionization emitter, named as a "nanoCESI" emitter, that had a thin-walled (∼10 μm) and tapered (5-10 μm) end. The thin conductive wall enabled sheathless ionization and minimized the flow rate of ionizing sample, and the tapered end efficiently ionized analytes via an electrospray ionization mechanism, providing up to 3.5-fold increase in sensitivity compared with a conventional sheathless emitter. Fifty repetitive analyses on 20 amino acids were successfully achieved with a nanoCESI emitter. Relative standard deviations of 50 analyses were 1.5%, 4.4%, and 6.8% for migration time, peak height, and peak area, respectively, where a limit of detection (LOD) of 170 pM (850 zmol) was achieved. Second, a sample enrichment method, large-volume dual preconcentration by isotachophoresis and stacking (LDIS), was applied to a newly designed protocol of nanoCESI-MS. This approach achieved up to 380-fold enhanced sensitivity and LOD of 450 fM. Compared with normal sheathless CE-MS, coupling of nanoCESI and LDIS provided up to 800-fold increase of sensitivity in total. Finally, metabolome analyses of single HeLa cells were performed, where 20 amino acids were successfully quantified with triple-quadrupole MS and 40 metabolites were identified with quadrupole-time-of-flight MS, as a promising analytical platform for microscale bioanalysis for the next generation.

Precision Analysis of the ^{136}Xe Two-Neutrino ββ Spectrum in KamLAND-Zen and Its Impact on the Quenching of Nuclear Matrix Elements

We present a precision analysis of the ^{136}Xe two-neutrino ββ electron spectrum above 0.8 MeV, based on high-statistics data obtained with the KamLAND-Zen experiment. An improved formalism for the two-neutrino ββ rate allows us to measure the ratio of the leading and subleading 2νββ nuclear matrix elements (NMEs), ξ_{31}^{2ν}=-0.26_{-0.25}^{+0.31}. Theoretical predictions from the nuclear shell model and the majority of the quasiparticle random-phase approximation (QRPA) calculations are consistent with the experimental limit. However, part of the ξ_{31}^{2ν} range allowed by the QRPA is excluded by the present measurement at the 90% confidence level. Our analysis reveals that predicted ξ_{31}^{2ν} values are sensitive to the quenching of NMEs and the competing contributions from low- and high-energy states in the intermediate nucleus. Because these aspects are also at play in neutrinoless ββ decay, ξ_{31}^{2ν} provides new insights toward reliable neutrinoless ββ NMEs.

Enhancement in acoustic focusing of micro and nanoparticles by thinning a microfluidic device

The manipulation of micro/nanoparticles has become increasingly important in biological and industrial fields. As a non-contact method for particle manipulation, acoustic focusing has been applied in sorting, enrichment and analysis of particles with microfluidic devices. Although the frequency and amplitude of acoustic waves and the dimensions of microchannels have been recognized as important parameters for acoustic focusing, the thickness of microfluidic devices has not been considered so far. Here, we report that thin glass microfluidic devices enhance acoustic focusing of micro/nanoparticles. It was found that the thickness of a microfluidic device strongly influences its ability to focus particles via acoustic radiation, because the energy propagation of acoustic waves is affected by the total mass of the device. Acoustic focusing of submicrometre polystyrene beads and Escherichia coli as well as enrichment of polystyrene beads were achieved in glass microfluidic devices as thin as 0.4 mm. Modifying the thickness of a microfluidic device can thus serve as a critical parameter for acoustic focusing when conventional parameters to achieve this effect are kept unchanged. Thus, our findings enable new approaches to the design of novel microfluidic devices.

A method of packaging molecule/cell-patterns in an open space into a glass microfluidic channel by combining pressure-based low/room temperature bonding and fluorosilane patterning

A fabrication method of a “post-molecule/cell patterned” glass microchip was developed by pressure-based bonding and patterning with a fluorosilane coupling reagent.

Dendritic cell CD83 homotypic interactions regulate inflammation and promote mucosal homeostasis

Dendritic cells (DCs) form an extensive network in the intestinal lamina propria, which orchestrates the mucosal immune response. Alterations in DC function can predispose to inflammatory bowel disease, although by unknown mechanisms. We show that CD83, a highly regulated DC cell surface protein, modulates the immune response to prevent colitis. Mice with a conditional knockout of CD83 in DCs develop exacerbated colitis following dextran sodium sulfate challenge, whereas mucosal overexpression of CD83 inhibits DC inflammatory response and protects against colitis. These CD83 perturbations can be modeled in vitro where we show that CD83 homotypic interaction occurs via cell-cell contact and inhibits pro-inflammatory responses. CD83 knockdown or cytoplasmic truncation abrogates the effects of homotypic binding. We demonstrate that CD83 homotypic interaction regulates DC activation via the mitogen-activated protein kinase pathway by inhibiting p38α phosphorylation. Our findings indicate that CD83 homotypic interactions regulate DC activation and promote mucosal homeostasis.

TNF superfamily member TL1A elicits type 2 innate lymphoid cells at mucosal barriers

Immune responses at mucosal barriers are regulated by innate type 2 lymphoid cells (ILC2s) that elaborate effector cytokines interleukins 5 and 13 (IL5 and IL13). IL25 and IL33 are key cytokines that support ILC2s; however, mice deficient in these pathways retain some functional ILC2s. Analysis of human and murine cells revealed that ILC2s highly express tumor necrosis factor (TNF)-receptor superfamily member DR3 (TNFRSF25). Engagement of DR3 with cognate ligand TL1A promoted ILC2 expansion, survival, and function. Exogenous protein or genetic overexpression of TL1A activated ILC2s independent of IL25 or IL33. Dr3(-/-) mice failed to control gut helminthic infections, and failed to mount ILC2 responses in the lung after nasal challenge with papain. Our data demonstrate a key role for TL1A in promoting ILC2s at mucosal barriers.

Estimation of regional cerebral blood flow using N-isopropyl-p-123I iodoamphetamine acquisition data from the lungs and brain

Aim: Previously, we devised a method for estimating 123I labeled N-isopropyl-p-iodo- amphetamine (123I IMP) arterial blood activity at 10 minutes after intravenous injection of 123I IMP (Ca10) without any blood sampling using 123I IMP autoradiography (ARG) acquisition data, and verified its usefulness for quantification of regional cerebral blood flow (rCBF). In this study, we attempted to develop an improved noninvasive method for estimating rCBF. Patients, methods: 123I IMP studies with 23 patients and 15O-H2O positron emission tomography (PET) ARG studies with 20 patients were evaluated. Multiple regression analysis was used to estimate an integral of the arterial blood counts during the time after injection of 123I (JCa) using parameters from the time series of the lung counts and brain counts as the explanatory variables and the fraction [brain single-photon emission computed tomography (SPECT) average count / the mean of rCBFs (mean CBF) measured by 15O-H2O PET ARG method] as the objective variable. Results: The regression equation was as follows: Estimated JCa = (7.09x10-3 · Cb12) - (1.57x10-4 · CbpreSPECT) + (9.48x10-5 · CbpostSPECT) + (1.35x10-4· L15) - (6.95x10-4· L33) + (7.61x10-4· L81) - (0.417), where Cb12: brain count at 12 minutes, Cbpre-SPECT: brain count before SPECT, Cbpost-SPECT: brain count after SPECT, L15, L33, and L81: lung count at 15, 33, and 81 seconds, respectively. The mean CBF values (ml/min/100g) calculated using the estimated JCa values more closely correlated with those measured by 15O-H2O PET ARG method (r = 0.833, p < 0.01) than those obtained by our previous method (r = 0.590, p < 0.01). Conclusion: The rCBFs obtained by this method approximated more accurately to the values measured by 15O-H2O PET ARG method than those obtained by our previous method.

Signals from the brainstem sleep/wake centers regulate behavioral timing via the circadian clock

Sleep-wake cycling is controlled by the complex interplay between two brain systems, one which controls vigilance state, regulating the transition between sleep and wake, and the other circadian, which communicates time-of-day. Together, they align sleep appropriately with energetic need and the day-night cycle. Neural circuits connect brain stem sites that regulate vigilance state with the suprachiasmatic nucleus (SCN), the master circadian clock, but the function of these connections has been unknown. Coupling discrete stimulation of pontine nuclei controlling vigilance state with analytical chemical measurements of intra-SCN microdialysates in mouse, we found significant neurotransmitter release at the SCN and, concomitantly, resetting of behavioral circadian rhythms. Depending upon stimulus conditions and time-of-day, SCN acetylcholine and/or glutamate levels were augmented and generated shifts of behavioral rhythms. These results establish modes of neurochemical communication from brain regions controlling vigilance state to the central circadian clock, with behavioral consequences. They suggest a basis for dynamic integration across brain systems that regulate vigilance states, and a potential vulnerability to altered communication in sleep disorders.

[Intraoperative direct angiography for the diminutive central pulmonary artery in a patient with major aortopulmonary collateral arteries]

It can be difficult to judge the degree of arborization of diminutive central pulmonary arteries (cPA) in patients with major aortopulmonary collateral arteries (MAPCA). Even through preoperative cardiac catheterization may not give adequate information. We introduce intra-operative direct angiography of diminutive cPA for patients with MAPCA. This would be one of the good options to judge the degree of arborization of the diminutive cPA, and to decide an initial surgical treatment. In this case, unifocalization of MAPCA without patch augmentation of pulmonary arteries, and an aortopulmonary shunt were performed at the 1st procedure. As enough growth of the cPA was obtained, this patient did not require additional patch augmentation of the pulmonary artery at the time of complete repair.

Opposing consequences of IL-23 signaling mediated by innate and adaptive cells in chemically induced colitis in mice

The interleukin-23 (IL-23) pathway has emerged as a promising therapeutic target for inflammatory bowel disease. Although the pathogenic role of IL-23 receptor (IL-23R) on T lymphocytes is well established, its function on innate immune cells has not been thoroughly examined. Here we investigate the consequence of IL-23R deletion in dextran sulfate sodium (DSS)-induced colitis. In IL23R(-/-) and IL23p19(-/-) mice, we observed decreased weight loss and reduced leukocyte infiltrate following DSS exposure. Surprisingly, when the IL-23R(-/-) allele was crossed into Rag2(-/-) mice, we observed exacerbated disease, increased epithelial damage, reduced pSTAT3 in the epithelium, and delayed recovery of IL23R(-/-)Rag2(-/-) mice. This phenotype was rescued with exogenous IL22-Fc, and epithelial pSTAT3 was restored. Depletion of Thy1(+) innate lymphoid cells eliminated the majority of IL-22 production in the colon lamina propria of DSS-treated Rag2(-/-) mice, suggesting that these are the major IL-23 responsive innate cells in this context. In summary, we provide evidence for opposing consequences of IL-23R on innate and adaptive lymphoid cells in murine colitis.

A novel pyridoxal 5'-phosphate-dependent amino acid racemase in the Aplysia californica central nervous system

D-aspartate (D-Asp) is found in specific neurons, transported to neuronal terminals and released in a stimulation-dependent manner. Because D-Asp formation is not well understood, determining its function has proved challenging. Significant levels of D-Asp are present in the cerebral ganglion of the F- and C-clusters of the invertebrate Aplysia californica, and D-Asp appears to be involved in cell-cell communication in this system. Here, we describe a novel protein, DAR1, from A. californica that can convert aspartate and serine to their other chiral form in a pyridoxal 5'-phosphate (PLP)-dependent manner. DAR1 has a predicted length of 325 amino acids and is 55% identical to the bivalve aspartate racemase, EC 5.1.1.13, and 41% identical to the mammalian serine racemase, EC 5.1.1.18. However, it is only 14% identical to the recently reported mammalian aspartate racemase, DR, which is closely related to glutamate-oxaloacetate transaminase, EC 2.6.1.1. Using whole-mount immunohistochemistry staining of the A. californica central nervous system, we localized DAR1-like immunoreactivity to the medial region of the cerebral ganglion where the F- and C-clusters are situated. The biochemical and functional similarities between DAR1 and other animal serine and aspartate racemases make it valuable for examining PLP-dependent racemases, promising to increase our knowledge of enzyme regulation and ultimately, D-serine and D-Asp signaling pathways.

Photochemical Smiles Rearrangement and Meisenheimer Complex Formation Catalyzed by Hydroxide Ion via Electron Hole Transfer Catalysis

[reactions: see text] Photochemical para-to-nitro Smiles rearrangement and para-to-nitro Meisenheimer complex formation occurs for nitrophenoxyethylamines with high concentrations of hydroxide ion in aqueous solution. Both photoreactions show first-order dependence on hydroxide ion concentration, but the mechanism involving hydroxide ion does not involve acid-base catalysis. The reactions take place from the triplet excited states of the nitrophenyl ethers. Analysis of quantum yields and kinetics is consistent with an electron hole transfer catalysis mechanism.

Estimation of plasma and saliva levels of coenzyme Q10 and influence of oral supplementation

Coenzyme Q10 (CoQ(10)) levels in human saliva were measured by HPLC with a highly sensitive electrochemical detector (ECD) and a special concentration column. This HPLC system showed satisfactory analytical results within the standard range of 0.78-50 ng/ml. We also found a significant correlation between CoQ(10) levels in plasma and in saliva from parotid glands, while this correlation was lacking between plasma CoQ10 and CoQ10 in whole saliva. Unlike in plasma, there are some fluctuations of saliva CoQ(10) levels throughout the day. A good correlation was obtained by collecting parotid gland saliva at times between meals. The mean saliva CoQ(10) level for 55 healthy volunteers was 17.0 ng/ml (S.D. 6.8 ng/ml); approximately one fiftieth of that in plasma. Regarding the influence of oral supplementation, CoQ(10) was analyzed in plasma and parotid gland saliva from 20 healthy volunteers supplemented daily with 100 mg of CoQ(10) for the first week and 200 mg for the second. The plasma CoQ(10) levels of all volunteers increased to different extents in accordance with the CoQ(10) daily intake and the corresponding change in saliva showed almost the same trend.

Differential expression of RANKL and osteoprotegerin in gingival crevicular fluid of patients with periodontitis

The receptor activator for NF-kappaB ligand (RANKL) plays an important role in osteoclast formation. A recent study with animal models suggests the involvement of RANKL in the pathogenesis of this periodontal disease. However, no one has examined the level of RANKL in the body fluid of human subjects. This communication reports on the in vivo concentrations of RANKL and the RANKL decoy receptor osteoprotegerin (OPG) in the gingival crevicular fluid (GCF) of periodontal subjects with severe, moderate, and mild forms of the disease. An increased concentration of RANKL and a decreased concentration of OPG were detected in GCF from patients with periodontitis (*p < 0.05 vs. control subjects). The ratio of the concentration of RANKL to that of OPG in the GCF was significantly higher for periodontal disease patients than for healthy subjects (*p < 0.01). Taken together, these data suggest that RANKL and OPG contribute to osteoclastic bone destruction in periodontal disease.

ABBREVIATIONS

GCF, gingival crevicular fluid; IL, interleukin; OPG, osteoprotegerin; RANKL, receptor activator for NF-kappaB ligand.

Association of a single nucleotide variant in the human tumour necrosis factor alpha promoter region with decreased bone mineral density

BACKGROUND

Tumour necrosis factor alpha (TNFalpha) has come to be regarded as a potential osteoporotic factor, because it has stimulatory effects on cells of the osteoclast lineage and has been implicated in the pathogenesis of bone loss associated with oestrogen deficiency. We recently described genetic linkage between the TNFalpha locus and human osteoporosis by sib-pair analysis. However, the molecular mechanism by which this locus regulates bone mineral density (BMD) remains unknown.

AIM

We investigated whether the observed linkage reflects a sequence variation which might affect expression of the TNFalpha gene or alter the function of TNFalpha protein.

SUBJECTS AND METHODS

We examined three single-nucleotide polymorphisms (SNPs) of the TNFalpha gene in a group of 390 postmenopausal Japanese women living in northern Japan. Minor-allele frequencies for the three SNPs (-1031C, -863A and -857T) in this population were 0.16, 0.13 and 0.20, respectively.

RESULTS

Among the three SNPs examined, we observed a significant correlation only between the presence of a T allele at nt -1031 and decreased BMD, by analysis of variance. Among the three genotypic groups at nt -1031, mean BMD values were significantly higher in the T-negative genotype (C/C homozygotes; mean SD = 0.342 +/- 0.052 g cm(-2)), compared with T-positive genotypes (T/T homozygotes, 0.309 +/- 0.062 g cm(-2); p = 0.0253 and T/C heterozygotes, 0.305 +/- 0.062 g cm(-2); p = 0.0164).

CONCLUSIONS

Given the lines of evidence from different genetic studies, we suggest that TNFalpha may play a role in pathogenesis of osteoporosis.

Molecular tryst peeping: detection of interactions between nonlabeled nucleic acids by fluorescence resonance energy transfer

We have developed a new method for monitoring the interactions between nonlabeled RNAs that involves detection of fluorescence resonance energy transfer (FRET) between two DNA probes with different fluorescent label. The sequences of the probes are complementary to those of the RNAs. In this study, we examined the interaction between a portion of the LTR RNA of HIV-1 and the corresponding antisense RNA. The antisense RNA was designed not to bind to the fluorescent DNA without prior hybridization to the target RNA. A mixture of RNAs and DNA probes with fluorescent labels was fractionated by electrophoresis on a nondenaturing polyacrylamide gel and then the gel was analyzed with a fluorescence imaging analyzer. FRET was observed only in the presence of target RNA, antisense RNA, and both of the fluorescent DNA probes. This strategy should be useful for the detection of interactions between nucleic acids that cannot be subjected to chemical modification, such as RNA transcripts inside cells.

Binding mode prediction for a flexible ligand in a flexible pocket using multi-conformation simulated annealing pseudo crystallographic refinement

We describe multi-conformation simulated annealing-pseudo-crystallographic refinement (MCSA-PCR), a technique developed for predicting the binding mode of a flexible ligand in a flexible binding pocket. To circumvent the local-minimum problem efficiently, this method performs multiple independent cycles of simulated annealing with explicit solvent, "growing" the ligand in the binding pocket each time. From the ensemble of structures, a pseudo-crystallographic electron density map is calculated, and then conventional crystallographic refinement methods are used to best fit a single, optimal structure into the density map. The advantage of the MCSA-PCR method is that it provides a direct means to evaluate the accuracy and uniqueness of the calculated solution, provides a measure of ligand and protein dynamics from the refined B-factors, and facilitates comparison with X-ray crystallographic data. Here, we show that our MCSA-PCR method succeeds in predicting the correct binding mode of the VSV8 peptide to the major histocompatibility complex (MHC) receptor. Importantly, there is a significant correlation between the experimentally determined crystallographic water molecules and water density observed in the pseudo map by MCSA-PCR. Furthermore, comparison of different approaches for extracting a single, most probable structure from the calculated ensemble reveals the power of the PCR method and provides insights into the nature of the energetic landscape.

A nucleotide variant in the promoter region of the interleukin-6 gene associated with decreased bone mineral density

Interleukin-6 (IL6) has come to be regarded as a potential osteoporotic factor because it has stimulatory effects on cells of the osteoclast lineage, and, thus, may play a role in the pathogenesis of bone loss associated with estrogen deficiency. We previously described association of the IL6 microsatellite with bone mineral density (BMD), as well as genetic linkage of the IL6 locus to human osteoporosis, by means of sib-pair analysis. However, the molecular mechanism by which this locus regulates BMD remains unknown. Accordingly, we searched for polymorphisms in the 5' and 3' flanking regions and in all five exons of the IL6 gene in a Japanese population sample. We identified three single-nucleotide sequence variations: a C/G substitution at nucleotide (nt) -634 in the promoter region, a G/A substitution at nt 4391 in the 3' noncoding region, and a variation in the AnTn tract around nt -447. The last of these had already been observed in Caucasians, as well as in Japanese. The single-nucleotide polymorphism at -634 created a restriction site for the BsrBI endonuclease, and the frequency of the minor (G) allele was 0.184. Five haplotypes were constructed among three variations examined in the population. Linkage disequilibrium was observed between the variation at -634 and the variation at 4391, as well as between the variation at -634 and the AnTn tract variation. We found a significant correlation, in 470 subjects, between the presence of the G allele and decreased BMD, by analysis of variance. When BMD values were compared among the three genotypic groups (G/G, G/C, C/C) at nt -634, BMD was lowest among the G/G homozygotes (mean +/- SD; 0.284 +/- 0.062g/cm2), highest among the C/C homozygotes (0.314 +/- 0.059g/cm2), and intermediate among the heterozygotes (0.303 +/- 0.066g/cm2; P < 0.05). Given the several lines of evidence from different genetic studies, we suggest that IL6 is, indeed, one of the genes affecting bone metabolism, in which variations can lead to osteoporosis.

Enzyme specificity under dynamic control II: Principal component analysis of alpha-lytic protease using global and local solvent boundary conditions

The contributions of conformational dynamics to substrate specificity have been examined by the application of principal component analysis to molecular dynamics trajectories of alpha-lytic protease. The wild-type alpha-lytic protease is highly specific for substrates with small hydrophobic side chains at the specificity pocket, while the Met190-->Ala binding pocket mutant has a much broader specificity, actively hydrolyzing substrates ranging from Ala to Phe. Based on a combination of multiconformation analysis of cryo-X-ray crystallographic data, solution nuclear magnetic resonance (NMR), and normal mode calculations, we had hypothesized that the large alteration in specificity of the mutant enzyme is mainly attributable to changes in the dynamic movement of the two walls of the specificity pocket. To test this hypothesis, we performed a principal component analysis using 1-nanosecond molecular dynamics simulations using either a global or local solvent boundary condition. The results of this analysis strongly support our hypothesis and verify the results previously obtained by in vacuo normal mode analysis. We found that the walls of the wild-type substrate binding pocket move in tandem with one another, causing the pocket size to remain fixed so that only small substrates are recognized. In contrast, the M190A mutant shows uncoupled movement of the binding pocket walls, allowing the pocket to sample both smaller and larger sizes, which appears to be the cause of the observed broad specificity. The results suggest that the protein dynamics of alpha-lytic protease may play a significant role in defining the patterns of substrate specificity. As shown here, concerted local movements within proteins can be efficiently analyzed through a combination of principal component analysis and molecular dynamics trajectories using a local solvent boundary condition to reduce computational time and matrix size.

Efficient transfer of intact oligonucleotides into the nucleus of ligament scar fibroblasts by HVJ-cationic liposomes is correlated with effective antisense gene inhibition

The efficacy of two different cationic liposomes, Lipofectin and hemagglutinating virus of Japan (HVJ)-cationic liposomes, on nuclear uptake of fluorescence-labeled phosphorothioate oligodeoxyribonucleotide (S-ODN) by ligament scar fibroblasts and suppression of decorin mRNA expression when antisense decorin S-ODN was transferred was investigated. There was no significant difference in nuclear uptake of fluorescent ODN between the two methods. However, only HVJ-cationic liposomes had a significant effect on suppression of decorin mRNA expression levels. To address the discrepancy, the molecular integrity of the transferred ODN in the cells was assessed by analysis of fluorescence resonance energy transfer (FRET) within double-fluorescence-labeled S-ODN. More than 70% of the ODN transfected by HVJ-cationic liposomes remained intact within the nucleus at 20 h after transfection, while the majority of the ODN transferred by Lipofectin was degraded at this point. These results suggest a strong relationship between the nuclear integrity of transfected antisense ODN and its suppression of target mRNA expression.

Five novel single-nucleotide polymorphisms of human interferon gamma identified by sequencing the entire gene

Interferon gamma (IFNG) plays important roles in the regulation of bone remodelling. We describe here six single-nucleotide polymorphisms (SNPs) in the IFNG gene, five of which are novel, and their allelic frequencies in the Japanese population, as determined by sequencing 48 alleles of the entire gene. Four of these polymorphisms were identified inside the third intron, at nucleotide (nt) positions 2459 (A/G), 2671 (T/C), 3177 (T/G), and 3273 (G/A). In exon 4, SNPs were identified at nt positions 5199 (A/T) and 5272 (A/G). These polymorphic sites will be useful for genetic studies of disorders that affect the inflammatory process or calcium metabolism.

Novel single nucleotide polymorphisms of the human colony-stimulating factor 2 (CSF2) gene identified by sequencing the entire gene

We describe three single nucleotide polymorphisms (SNPs) of the human colony-stimulating factor 2 (CSF2) gene and their allelic frequencies, as determined by direct sequencing of 48 alleles of the entire CSF2 gene. Three polymorphisms were identified, at nucleotide positions 1816 (T/C), 2284 (C/T), and 3079 (G/A). These polymorphisms will be useful in genetic studies not only of hematologic disorders but also of disorders of bone metabolism.

Novel single nucleotide polymorphisms of the human nuclear factor kappa-B 2 gene identified by sequencing the entire gene

The nuclear factor kappa-B 2 (NFKB2) gene is a member of the NFKB/Rel gene family, which is known to be a pivotal regulator of the acute phase of the inflammatory response and of immune responses. We identified three novel single nucleotide polymorphisms (SNPs) and determined their allelic frequencies, as determined by the sequencing of 48 alleles of the entire gene in a Japanese population sample. Two of the three polymorphisms were identified at nucleotide (nt) position 1837 (T/C) and nt position, 1867 (GG/G) in the upstream region of the gene. The other polymorphism was identified at nt position 2584 (G/T) within intron 1. These polymorphisms will be useful in genetic studies of the processes involved in inflammatory responses and in bone differentiation.

Linkage of human tumor necrosis factor-alpha to human osteoporosis by sib pair analysis

Osteoporosis as well as osteopenia are common human conditions considered to result from the interplay of multiple genetic and environmental factors. Twin and family studies have yielded strong correlation between measures of bone mass and a number of genetic factors. Certain genes (e.g., cytokines such as interleukin-1, interleukin-6, or tumor necrosis factor-alpha) are capable of regulating metabolism, formation, and resorption of bone; all processes that determine bone mass. We tested 192 sib-pairs of adult Japanese women from 136 families for genetic linkage between osteoporosis and osteopenia phenotypes and allelic variants at the tumor necrosis factor-alpha (TNFA) locus, using a dinucleotide-repeat polymorphism located near the gene. The TNFA locus showed evidence for linkage to osteoporosis, with mean allele sharing of 0.478 (P = 0.30) in discordant pairs and 0.637 (P = 0.001) in concordant affected pairs. Linkage with osteopenia was also significant in concordant affected pairs (P = 0.017). Analyses limited to the post-menopausal women in our cohort showed similar or even stronger linkage for both phenotypes. The results provide evidence that genetic variations within the TNFA locus or adjacent genes affect regulation of mineral metabolism in bone and some of them confer risk for osteoporosis in adult women.

Cloning and functional expression of rat kidney dipeptidyl peptidase II

Dipeptidyl peptidase II (DPP II; EC 3.4.14.2) from rat kidney was purified to a specific activity of 65.4 micromol/min per mg of protein for Lys-Ala-beta-naphthylamide. The N-terminal and partial amino acid sequences of the enzyme were determined. The peptide sequences were used to identify expressed sequence tag (EST) clones. By using the cDNA fragment of one of the EST clones as a probe, we isolated a cDNA clone with 1710 bp encoding DPP II from a rat kidney cDNA library. The cDNA of rat DPP II contained an open reading frame of 1500 bp, coding for a protein of 500 amino acids. The first 10 residues of the purified enzyme matched the deduced protein sequence starting with residue 37, suggesting the presence of a signal peptide. The mature enzyme (464 residues) had a calculated molecular mass of 51400 Da, which was lower than the value (about 60000 Da) determined by SDS/PAGE; and the deduced amino acid sequence showed six potential N-glycosylation sites. The deduced amino acid sequence of rat DPP II shared high similarity with quiescent-cell proline dipeptidase (78% identity) and prolyl carboxypeptidase (38% identity) and bore the putative catalytic triad (Ser, Asp, His) conserved in serine peptidase families. We transiently transfected COS-7 cells with pcDNA3.1 containing the cloned cDNA and obtained the overexpression of an immunoreactive protein (of about 60000 Da). The transfected cells showed Lys-Ala-methylcoumarinamide-hydrolysing activity that was 50 times higher than the control cells.

Association of bone mineral density with polymorphism of the human calcium-sensing receptor locus

A strong correlation between bone mass and genetic factors has been shown in twins and family studies. Some of the genes involved would regulate bone metabolism, bone formation, and resorption, all processes that determine bone mass. One candidate genes, calcium-sensing receptor (CASR) in the parathyroid gland, regulates calcium homeostasis by sensing decreases in extracellular calcium level and effecting an increase in secretion of parathyroid hormone (PTH) and calcium (Ca) reabsorption in the kidney. We have investigated a possible association between the CA-repeat polymorphism at the human CASR gene locus and the bone mineral density (BMD) of radial bone in 472 postmenopausal Japanese women. Genotypes were classified into nine groups according to the number of CA repeats present, from 20 to 12. BMD was expressed as the adjusted BMD, which was the body mass index (BMI), and age-adjusted average BMD. The 247 women who had an A3 allele (228 bp, containing 18 repeats of CA) had significantly lower adjusted BMD (mean +/- SD: 0.303 +/- 0.059 versus 0.316 +/- 0.063 g/cm(2); P = 0.0308) than the participants (n = 201) who did not carry an allele of that size. This result suggests that genetic variation at the CASR gene locus is associated with some determinants for BMD in postmenopausal women.

Non-Boltzmann thermodynamic integration (NBTI) for macromolecular systems: relative free energy of binding of trypsin to benzamidine and benzylamine

The relative free energies of binding of trypsin to two amine inhibitors, benzamidine (BZD) and benzylamine (BZA), were calculated using non-Boltzmann thermodynamic integration (NBTI). Comparison of the simulations with the crystal structures of both complexes, trypsin-BZD and trypsin-BZA, shows that NBTI simulations better sample conformational space relative to thermodynamic integration (TI) simulations. The relative binding free energy calculated using NBTI was much closer to the experimentally determined value than that obtained using TI. The error in the TI simulation was found to be primarily due to incorrect sampling of BZA's conformation in the binding pocket. In contrast, NBTI produces a smooth mutation from BZD to BZA using a surrogate potential, resulting in a much closer agreement between the inhibitors' conformations and the omit electron density maps. This superior agreement between experiment and simulation, of both relative binding free energy differences and conformational sampling, demonstrates NBTI's usefulness for free energy calculations in macromolecular simulations.

A common Ile796Val polymorphism of the human SREBP cleavage-activating protein (SCAP) gene

We identified a new common amino acid polymorphism of isoleucine/valine at codon 796 in exon 16 of the gene for human sterol regulatory element binding protein (SREBP) cleavage-activating protein (SCAP), a central regulator of lipid synthesis and metabolism in animal cells. It can be detected as an MslI restriction fragment length polymorphism. The allelic frequencies were: isoleucine (A) allele, 0.57 and valine (G) allele, 0.43. This polymorphism may be useful for genetic studies of disorders affecting intracellular lipid metabolism and hyperlipidemia.

Allelic variants in the interleukin-6 gene and essential hypertension in Japanese women

Genes that can be implicated in hypertension in experimental animals are plausible candidates in the pathogenesis of human hypertension. A recent genome-wide search for quantitative-trait loci (QTL) in hypertensive rats revealed a strong correlation between the interleukin-6 (IL-6) locus on rat chromosome 4 and systolic, diastolic, and mean arterial pressure in this mammalian species. To investigate a possible association between genetic variations of the IL-6 gene and hypertension in humans, we identified two novel single-nucleotide sequence variations, a C/G substitution at -634 in the promoter region and a G/A substitution at 4391 in a 3' non-coding portion of exon 5, and a previous reported sequence variant, an A/T variation in the composition of the AnTn tract around -447 in the promoter region (Fishman D et al. J Clin Invest 1998; 102: 1369-1376), within a test population of 96 Japanese subjects. Allelic associations involving these variations were analyzed in 150 hypertensive and 143 normotensive Japanese women. The distribution of alleles of the three polymorphisms, as well as a dinucleotide repeat present at the IL-6 locus, was similar in the two groups. Therefore, the IL-6 gene appears to play a minimal role in the genetic etiology of essential hypertension in Japanese women.