Measurement and analysis of internal stress distributions created in gelatin simulated-brain tissue by a pulsed laser-induced liquid jet.

Transsphenoidal surgery is currently employed to treat complex lesions beyond the sella turcica; however, the procedure can be limited by difficulties encountered in dealing with small blood vessels, deep and narrow working spaces, and awkward working angles. To overcome these problems, we have developed a pulsed laser-induced liquid jet system that can dissect tumor tissue while preserving fine blood vessels within deep and narrow working spaces. We have previously evaluated the utility and safety of this procedure. However, the effects of the pulsejet after being injected into the brain are not yet well understood. Especially, the behavior of the stress distribution created by the jet is important because it has recently been reported that high acoustic pressures can affect the brain. In this study, we measured internal stress distributions in a gelatin simulated-brain using photoelasticity experiments. We used a high-speed camera with an image sensor on which an array of micropolarizers was attached to measure the stresses and the shear wave created when the pulsejet enters the simulated brain.

Analysis of liquid reservoir effect induced by pulsed laser liquid jet.

A pulsed-laser-induced liquid jet (LILJ) is a new device used in neurosurgery to simultaneously crush, incise, and aspirate tissues and tumors, preserving blood vessels and nerves. In addition, a feature of a pulsed LILJ is its ability to excavate tissue at constant depth while a liquid jet is being repeatedly focused at the same point. To clarify the mechanisms of constant depth of excavation, we employed a gelatin phantom and extracted brain tissue using a high-speed camera, and we then confirmed that the liquid-reservoir-induced LILJ played an important role in enabling the safe usage of an LILJ.

Evaluation of accuracy of non-linear finite element computations for surgical simulation: study using brain phantom.

In this paper, the accuracy of non-linear finite element computations in application to surgical simulation was evaluated by comparing the experiment and modelling of indentation of the human brain phantom. The evaluation was realised by comparing forces acting on the indenter and the deformation of the brain phantom. The deformation of the brain phantom was measured by tracking 3D motions of X-ray opaque markers, placed within the brain phantom using a custom-built bi-plane X-ray image intensifier system. The model was implemented using the ABAQUS(TM) finite element solver. Realistic geometry obtained from magnetic resonance images and specific constitutive properties determined through compression tests were used in the model. The model accurately predicted the indentation force-displacement relations and marker displacements. Good agreement between modelling and experimental results verifies the reliability of the finite element modelling techniques used in this study and confirms the predictive power of these techniques in surgical simulation.

Subject-specific non-linear biomechanical model of needle insertion into brain.

The previous models for predicting the forces acting on a needle during insertion into very soft organs (such as, e.g. brain) relied on oversimplifying assumptions of linear elasticity and specific experimentally derived functions for determining needle-tissue interactions. In this contribution, we propose a more general approach in which the needle forces are determined directly from the equations of continuum mechanics using fully non-linear finite element procedures that account for large deformations (geometric non-linearity) and non-linear stress-strain relationship (material non-linearity) of soft tissues. We applied these procedures to model needle insertion into a swine brain using the constitutive properties determined from the experiments on tissue samples obtained from the same brain (i.e. the subject-specific constitutive properties were used). We focused on the insertion phase preceding puncture of the brain meninges and obtained a very accurate prediction of the needle force. This demonstrates the utility of non-linear finite element procedures in patient-specific modelling of needle insertion into soft organs such as, e.g. brain.

Comparative study of translation termination sites and release factors (RF1 and RF2) in procaryotes.

Translation termination is catalyzed by release factors that recognize stop codons. However, previous works have shown that in some bacteria, the termination process also involves bases around stop codons. Recently, Ito et al. analyzed release factors and identified the amino acids therein that recognize stop codons. However, the amino acids that recognize bases around stop codons remain unclear. To identify the candidate amino acids that recognize the bases around stop codons, we aligned the protein sequences of the release factors of various bacteria and searched for amino acids that were conserved specifically in the sequence of bacteria that seemed to regulate translation termination by bases around stop codons. As a result, species having several highly conserved residues in RF1 and RF2 showed positive correlations between their codon usage bias and conservation of the bases around the stop codons. In addition, some of the residues were located very close to the SPF motif, which deciphers stop codons. These results suggest that these conserved amino acids enable the release factors to recognize the bases around the stop codons.

On biased distribution of introns in various eukaryotes.

We conducted comprehensive analyses on intron positions in the Mus musculus genome by comparing genomic sequences in the GenBank database and cDNA sequences in the mouse cDNA library recently developed by Riken Genomic Sciences Center. Our results confirm that introns have a tendency to be located toward the 5' end of the gene. The same type of analysis was conducted in the coding region of seven eukaryotes (Saccharomyces cerevisiae, Plasmodium falciparum, Caenorhabditis elegans, Drosophila melanogaster, M. musculus, Homo sapiens, Arabidopsis thaliana). Introns in genes with a single intron have a locational bias toward the 5' end in all species except A. thaliana. We also measured the distance from the start codon to the position of the intron, and found that single introns prefer the location immediately after the start codon in S. cerevisiae and P. falciparum. We discuss three possible explanations for these findings: (1) they are the consequence of intron loss by reverse-transcriptase; (2) they are necessary to accommodate the function; and (3) they are concerned with the mechanism of pre-mRNA splicing.

Comprehensive sequence analysis of translation termination sites in various eukaryotes.

Recent investigations into the translation termination sites of various organisms have revealed that not only stop codons but also sequences around stop codons have an effect on translation termination. To investigate the relationship between these sequence patterns and translation as well as its termination efficiency, we analysed the correlation between strength of consensus and translation efficiency, as predicted according to Codon Adaptation Index (CAI) value. We used RIKEN full-length mouse cDNA sequences and ten other eukaryotic UniGene datasets from NCBI for the analyses. First, we conducted sequence profile analyses following translation termination sites. We found base G and A at position +1 as a strong consensus for mouse cDNA. A similar consensus was found for other mammals, such as Homo sapiens, Rattus norvegicus and Bos taurus. However, some plants had different consensus sequences. We then analysed the correlation between the strength of consensus at each position and the codon biases of whole coding regions, using information content and CAI value. The results showed that in mouse cDNA, CAI value had a positive correlation with information content at positions +1. We also found that, for positions with strong consensus, the strength of the consensus is likely to have a positive correlation with CAI value in some other eukaryotes. Along with these observations, biological insights into the relationship between gene expression level, codon biases and consensus sequence around stop codons will be discussed.

Identification of human drug-metabolizing enzymes involved in the metabolism of SNI-2011.

In vitro studies were conducted to identify human drug-metabolizing enzymes involved in the metabolism of SNI-2011 ((+/-)-cis-2-methylspiro [1,3-oxathiolane-5,3'-quinuclidine] monohydrochloride hemihydrate, cevimeline hydrochloride hydrate). When 14C-SNI-2011 was incubated with human liver microsomes, SNI-2011 trans-sulfoxide and cis-sulfoxide were detected as major metabolites. These oxidations required NADPH, and were markedly inhibited by SKF-525A, indicating that cytochrome P450 (CYP) was involved. In a chemical inhibition study, metabolism of SNI-2011 in liver microsomes was inhibited (35-65%) by CYP3A4 inhibitors (ketoconazole and troleandomycin) and CYP2D6 inhibitors (quinidine and chlorpromazine). Furthermore, using microsomes containing cDNA-expressed CYPs, it was found that high rates of sulfoxidation activities were observed with CYP2D6 and CYP3A4. On the other hand, when 14C-SNI-2011 was incubated with human kidney microsomes, SNI-2011 N-oxide was identified as a major metabolite. This N-oxidation required NADPH, and was completely inhibited by thiourea, indicating that flavin-containing monooxygenase (FMO) was involved. In addition, microsomes containing cDNA-expressed FMO1, a major isoform in human kidney, mainly catalyzed N-oxidation of SNI-2011, but microsomes containing FMO3, a major isoform in adult human liver, did not. These results suggest that SNI-2011 is mainly catalyzed to sulfoxides and N-oxide by CYP2D6/3A4 in liver and FMOI in kidney, respectively.

Correlation between sequence conservation of the 5' untranslated region and codon usage bias in Mus musculus genes.

The codon adaptation index (CAI) values of all protein-coding sequences of the full-length cDNA libraries of Mus musculus were computed based on the RIKEN mouse full-length cDNA library. We have also computed the extent of consensus in flanking sequences of the initiator ATG codon based on the 'relative entropy' values of respective nucleotide positions (from -20 to +12 bp relative to the initiator ATG codon) for each group of genes classified by CAI values. With regard to the two nucleotides positions (-3 and +4) known to be highly conserved in Kozak's consensus sequence, a clear correlation between CAI values and relative entropy values was observed at position -3 but this was not significant at position +4, although a significant correlation was found at position -1 of the consensus sequence. Further, although no correlation was observed at any additional positions, relative entropy values were very high at positions -4, -6, and -8 in genes with high CAI values. These findings suggest that the extent of conservation in the flanking sequence of the initiator ATG codon including Kozak's consensus sequence was an important factor in modulation of the translation efficiency as well as synonymous codon usage bias particularly in highly expressed genes.

Preferential usage of some minor codons in bacteria.

In many bacterial species, such as Deinococcus radiodurans, Haemophilus influenzae, and Methanobacterium thermoautotrophicum, some minor codons are preferentially used near the initiation codon. Among these codons, there are some minor codons that have strong preference for the initiation site in the high codon adaptation index (CAI) group (comprising of highly expressed genes) rather than in the low CAI genes group (comprising of low expressing genes). In the present study, codon usage in the initiation site and in the rest of the gene was systematically compared in the 27 complete bacterial genomes and Saccharomyces cerevisiae genome. Furthermore, we classified genes into two groups according to the CAI values and conducted the same analysis for each of the two groups. Our results suggest a role for some minor codons in the initiation site of the regulating translation system in many bacteria. We have summarized codons that are preferentially used in the initiation site and probably play a role in regulating genes expression in these organisms.

An experimental study on the interface strength between titanium mesh cage and vertebra in reference to vertebral bone mineral density.

STUDY DESIGN: Using human cadaver spines, the authors investigated mechanical properties of the interface between titanium mesh cage and vertebra in respect to vertebral bone mineral density. OBJECTIVES: The objective of this study is to examine the effects of the size of the mesh cage and an internal end ring system on interface mechanical properties in reference to vertebral bone mineral density. SUMMARY OF BACKGROUND DATA: A titanium mesh cage has recently been developed for anterior spinal reconstruction. The cage provides immediate postoperative stability and facilitates bony union with cancellous bone packed in the cage itself. Mechanical properties of the interface between the cage and vertebra, however, are yet to be clarified in osteoporotic spine. METHODS: Twenty-five lumbar vertebrae harvested from embalmed human cadavers (n = 20) were used. The vertebrae were divided into four experimental groups according to the applied cage conditions: Phi25 mm cage without internal end ring (L-), Phi19 mm cage without internal end ring (S-), Phi25 mm cage with internal end ring (L+), and Phi19 mm cage with internal end ring (S+). Bone mineral density of whole vertebral body was measured by dual energy radiograph absorptiometer (DXA). Peripheral quantitative computed tomography was used to determine local bone mineral density of subchondral cancellous bone of vertebral body. Each cage was compressed on vertebral endplate via a specially designed device connected to a material testing machine. Maximum load and stiffness of the interface between the cage and vertebra were measured from load-deformation data in quasi-static compression loading with a loading rate of 0.5 mm/min. Relationships between the mechanical properties and vertebral bone mineral densities were evaluated. In 11 specimens acoustic emission during compression loading was measured and simultaneously recorded in load-deformation data. After the mechanical test microradiograms of midsagittal sections of the vertebrae were taken to observe failure patterns of endplate or trabecular bone. RESULTS: Vertebrae compressed with large cages (group: L- or L+) showed greater maximum load than those compressed with small cages (group: S- or S+). The internal end ring contributed to higher maximumload. The size of the cage or the internal end ring, however, did not have any effect on stiffness. Maximum load and stiffness were positively correlated with whole vertebral bone mineral density measured by dual energy radiograph absorptiometer or local cancellous bone mineral density of subchondral bone measured by peripheral quantitative computed tomography. Correlation coefficient and P value were more significant in the association of the mechanical properties and subchondral bone mineral density measured by peripheral quantitative computed tomography than in the association of the parameters and whole vertebral bone mineral density measured by DXA. A load-deformation curve with an acoustic emission event count rate showed that significant acoustic emission signals were generated around maximum load. On microradiographic study most vertebrae compressed with the cage showed encroachment of the cage spikes into the endplate or trabecular structure, preserving structures of the most central portion of the vertebrae. CONCLUSION: A titanium mesh cage with larger diameter and/or augmentation of internal end ring produces a significant increase of the interface strength between the cage and the vertebra. A positive correlation between the interface strength and vertebral bone mineral density suggests that vertebral bone mineral density is an important parameter for successful spinal reconstruction, and also implies that in severe osteoporotic spine the stability of the cage is declined, and other instrumentation should be combined.

Correlation between Shine--Dalgarno sequence conservation and codon usage of bacterial genes.

In this study, we analyzed the correlation between codon usage bias and Shine--Dalgarno (SD) sequence conservation, using complete genome sequences of nine prokaryotes. For codon usage bias, we adopted the codon adaptation index (CAI), which is based on the codon usage preference of genes encoding ribosomal proteins, elongation factors, heat shock proteins, outer membrane proteins, and RNA polymerase subunit proteins. To compute SD sequence conservation, we used SD motif sequences predicted by Tompa and systematically aligned them with 5'UTR sequences. We found that there exists a clear correlation between the CAI values and SD sequence conservation in the genomes of Escherichia coli, Bacillus subtilis, Haemophilus influenzae, Archaeoglobus fulgidus, Methanobacterium thermoautotrophicum, and Methanococcus jannaschii, and no relationship is found in M. genitalium, M. pneumoniae, and Synechocystis. That is, genes with higher CAI values tend to have more conserved SD sequences than do genes with lower CAI values in these organisms. Some organisms, such as M. thermoautotrophicum, do not clearly show the correlation. The biological significance of these results is discussed in the context of the translation initiation process and translation efficiency.

Functional annotation of a full-length mouse cDNA collection.

The RIKEN Mouse Gene Encyclopaedia Project, a systematic approach to determining the full coding potential of the mouse genome, involves collection and sequencing of full-length complementary DNAs and physical mapping of the corresponding genes to the mouse genome. We organized an international functional annotation meeting (FANTOM) to annotate the first 21,076 cDNAs to be analysed in this project. Here we describe the first RIKEN clone collection, which is one of the largest described for any organism. Analysis of these cDNAs extends known gene families and identifies new ones.

Causal analysis of CpG suppression in the Mycoplasma genome.

Some bacterial genomes are known to have low CpG dinucleotide frequencies. While their causes are not clearly understood, the frequency of CpG is suppressed significantly in the genome of Mycoplasma genitalium, but not in that of Mycoplasma pneumoniae. We compared orthologous gene pairs of the two closely related species to analyze CpG substitution patterns between these two genomes. We also divided genome sequences into three regions: protein-coding, noncoding, and RNA-coding, and obtained the CpG frequencies for each region for each organism. It was found that the observed/expected ratio of CpG dinucleotides is low in both the protein-coding and noncoding regions; while that ratio is in the normal range in the RNA-coding region. Our results indicate that CpG suppression of the Mycoplasma genome is not caused by (1) biased usage amino acid; (2) biased usage of synonymous codon; or (3) methylation effects by the CpG methyltransferase in the genomes of their hosts. Instead, we consider it likely that a certain global pressure, such as genome-wide pressure for the advantages of DNA stability or replication, has the effect of decreasing CpG over the entire genome, which, in turn, resulted in the biased codon usage.

An experimental study of porcine lumbar segmental stiffness by the distraction-compression principle using a threaded interbody cage.

The objectives of this study were to quantify changes in stiffness and disk height of porcine functional spinal units (FSUs) by installation of a threaded interbody cage and those by gradual resection of the annulus fibrosus. Flexion, extension, bending, and torsion to the FSUs were performed in four sequential stages: stage I, intact FSU; stage II, the FSUs are fitted with a threaded fusion cage; stage III, the FSUs are fitted with a threaded fusion cage with the anterior one third of the annulus fibrosus excised, including excision of the anterior longitudinal ligament; and stage IV, in addition to stage III, the bilateral annulus fibrosus is excised. Segmental stiffness in each loading in the four stages and a change of disk height induced by the instrumentation were measured. After instrumentation, stiffnesses in all loading modes (p < 0.005) and disk height (p = 0.002) increased significantly. The stiffnesses of FSUs fixed by the cage decreased with gradual excision of the annulus fibrosus in flexion, extension, and bending. These results suggest that distraction of the annulus fibrosus and posterior ligamentous structures by installation of the cage increases the soft-tissue tension, resulting in compression to the cage and a stiffer motion segment.

Two-stage method for protein-ligand docking.

A two-stage method for the computational prediction of the structure of protein-ligand complexes is proposed. Given an experimentally determined structure of the protein, in the first stage a large number of plausible ligand conformations is generated using the fast docking algorithm FlexX. In the second stage these conformations are minimized and reranked using a method based on a classical force field. The two-stage method is tested for 10 different protein-ligand complexes. For 9 of them experimentally determined structures are known. It turns out that the two-stage method strongly improves the predictive power as compared to that of the fast docking stage alone. The tenth case is a bona fide prediction of a complex of thrombin with a new inhibitor for which no experimentally determined structure is available so far.

Comparative study of overlapping genes in the genomes of Mycoplasma genitalium and Mycoplasma pneumoniae.

Overlapping genes are defined, in this paper, as a pair of adjacent genes whose coding regions are partly overlapping. We systematically analyzed all overlapping genes in the genomes of two closely related species: Mycoplasma genitalium and Mycoplasma pneumoniae. Careful comparisons were made for homologous genes that are overlapped in one species but not in the other. This comparative analysis allows us to propose a model of how overlapping genes emerged in the course of evolution. It was found that overlapping genes were generated primarily due to the loss of a stop codon in either gene, in many cases, the absence of which resulted in elongation of the 3' end of the gene's coding region. More specifically, the loss of the stop codon took place as a result of the following events: deletion of the stop codon (64.4%), point mutation at the stop codon (4.4%), and frame shift at the end of the coding region (6.7%). Overlapping genes, in a sense, can be thought of as the results of evolutionary pressure to minimize genome size. However, our analysis indicates that many overlapping genes, at least in the genomes of M.genitalium and M.pneumoniae, are due to incidental elongation of the coding regions.

Analysis of complete genomes suggests that many prokaryotes do not rely on hairpin formation in transcription termination.

Free energy values of mRNA tertiary structures around stop codons were systematically calculated to surmise the hairpin-forming potential for all genes in each of the 16 complete prokaryote genomes. Instead of trying to detect each individual hairpin, we averaged the free energy values around the stop codons over the entire genome to predict how extensively the organism relies on hairpin formation in the process of transcription termination. The free energy values of Escherichia coli K-12 shows a sharp drop, as expected, at 30 bp downstream of the stop codon, presumably due to hairpin-forming sequences. Similar drops are observed for Haemophilus influenzae Rd, Bacillus subtilis and Chlamydia trachomatis, suggesting that these organisms also form hairpins at their transcription termination sites. On the other hand, 12 other prokaryotes- Mycoplasma genitalium, Mycoplasma pneumoniae, Synechocystis PCC6803, Helicobacter pylori, Borrelia burgdorferi, Methanococcus jannaschii, Archaeoglobus fulgidus, Methanobacterium thermoautotrophicum, Aquifex aeolicus, Pyrococcus horikoshii, Mycobacterium tuberculosis and Treponema pallidum -show no apparent decrease in free energy value at the corresponding regions. This result suggests that these prokaryotes, or at least some of them, may never form hairpins at their transcription termination sites.

Fracture risk during pedicle screw insertion in osteoporotic spine.

Although using a larger screw can enhance the stability of screw-bone interface in pedicular screw instrumentation, it may cause pedicle fracture during screw insertion in osteoporotic spine. We investigated structural changes of the pedicle with the advance of osteoporosis and its implication for the risk of pedicle fracture. Bone mineral density of trabecular, subcortical, and cortical bone of pedicle in an osteoporotic group, assessed by peripheral quantitative computed tomography, was significantly lower than that of those in the normal group. Cortical thickness also declined. There were no pedicle fractures in normal vertebrae, but seven (41.2%) fractures were observed among 17 pedicles in osteoporotic vertebrae in which bone mineral density measured with dual-energy x-ray absorptiometry was <0.7 g/cm2 and screw diameter >70% of the outer diameter of the pedicle. Our results suggest that screw diameter should not exceed 70% of the outer diameter of the pedicle in the case of osteoporosis in which vertebral bone mineral density is <0.7 g/cm2.

Structural characteristics of the pedicle and its role in screw stability.

STUDY DESIGN: Cross-sectional regional bone mineral density of the pedicle was measured by peripheral quantitative computed tomography. Biomechanical tests were performed to clarify the role of the pedicle in screw stability. OBJECTIVES: To identify the structural characteristics of the pedicle that supports pedicle screw stability and the differences in these characteristics between normal and osteoporotic vertebrae. SUMMARY OF BACKGROUND DATA: The pedicle screw is an essential component of many systems used to align the spine. The contribution of the pedicle to screw stability, however, has not been fully investigated. METHODS: Trabecular, subcortical, and cortical bone mineral density and the area of the pedicle were measured by peripheral quantitative computed tomography. Bone mineral density also was recalculated in four circumferential layers. These parameters were compared between normal and osteoporotic individuals. The relative contribution of the pedicle to screw stability was evaluated by caudocephalad and pull-out loading in a vertebra with or without its body. RESULTS: Inner trabecular, middle subcortical, and outer cortical bone mineral density and cortical bone area in the pedicle were significantly lower in osteoporotic vertebrae than those in normal vertebrae. In the pedicle, bone mineral density increased close to the outer layer. Bone mineral density not as thick even in the outer layer in osteoporotic subjects. Approximately 80% of the caudocephalad stiffness and 60% of the pullout strength of the pedicle screw depended on the pedicle rather than on the vertebral body. CONCLUSION: Screw stability depends on the structural characteristics of the pedicle. The pedicle was denser in the subcortical bone, in which the threads of the screw engage, than in trabecular bone. In osteoporosis, bone mineral density was not as dense even in the outer layer, and the cortex was thinner than normal. A larger screw would not enhance screw stability and may break the thin cortex in osteoporotic vertebrae.