Calculation of dose kernels for radiotherapy applications using the GEANT 4 Monte Carlo toolkit

Dose kernels may be used to calculate dose distributions in radiotherapy (as described by Ahnesjo et al., 1999). Their calculation requires use of Monte Carlo methods, usually by forcing interactions to occur at a point. The Geant4 Monte Carlo toolkit provides a capability to force interactions to occur in a particular volume. We have modified this capability and created a Geant4 application to calculate dose kernels in cartesian, cylindrical, and spherical scoring systems. The simulation considers monoenergetic photons incident at the origin of a 3 m x 3 x 9 3 m water volume. Photons interact via compton, photo-electric, pair production, and rayleigh scattering. By default, Geant4 models photon interactions by sampling a physical interaction length (PIL) for each process. The process returning the smallest PIL is then considered to occur. In order to force the interaction to occur within a given length, L_FIL, we scale each PIL according to the formula: PIL_forced = L_FIL 9 (1 - exp(-PIL/PILo)) where PILo is a constant. This ensures that the process occurs within L_FIL, whilst correctly modelling the relative probability of each process. Dose kernels were produced for an incident photon energy of 0.1, 1.0, and 10.0 MeV. In order to benchmark the code, dose kernels were also calculated using the EGSnrc Edknrc user code. Identical scoring systems were used; namely, the collapsed cone approach of the Edknrc code. Relative dose difference images were then produced. Preliminary results demonstrate the ability of the Geant4 application to reproduce the shape of the dose kernels; median relative dose differences of 12.6, 5.75, and 12.6 % were found for an incident photon energy of 0.1, 1.0, and 10.0 MeV respectively.

Recent emergence of dengue virus serotype 4 in French Polynesia results from multiple introductions from other south pacific islands

BACKGROUND: Infection by dengue virus (DENV) is a major public health concern in hundreds of tropical and subtropical countries. French Polynesia (FP) regularly experiences epidemics that initiate, or are consecutive to, DENV circulation in other South Pacific Island Countries (SPICs). In January 2009, after a decade of serotype 1 (DENV-1) circulation, the first cases of DENV-4 infection were reported in FP. Two months later a new epidemic emerged, occurring about 20 years after the previous circulation of DENV-4 in FP. In this study, we investigated the epidemiological and molecular characteristics of the introduction, spread and genetic microevolution of DENV-4 in FP. METHODOLOGY/PRINCIPAL FINDINGS: Epidemiological data suggested that recent transmission of DENV-4 in FP started in the Leeward Islands and this serotype quickly displaced DENV-1 throughout FP. Phylogenetic analyses of the nucleotide sequences of the envelope (E) gene of 64 DENV-4 strains collected in FP in the 1980s and in 2009-2010, and some additional strains from other SPICs showed that DENV-4 strains from the SPICs were distributed into genotypes IIa and IIb. Recent FP strains were distributed into two clusters, each comprising viruses from other but distinct SPICs, suggesting that emergence of DENV-4 in FP in 2009 resulted from multiple introductions. Otherwise, we observed that almost all strains collected in the SPICs in the 1980s exhibit an amino acid (aa) substitution V287I within domain I of the E protein, and all recent South Pacific strains exhibit a T365I substitution within domain III. CONCLUSIONS/SIGNIFICANCE: This study confirmed the cyclic re-emergence and displacement of DENV serotypes in FP. Otherwise, our results showed that specific aa substitutions on the E protein were present on all DENV-4 strains circulating in SPICs. These substitutions probably acquired and subsequently conserved could reflect a founder effect to be associated with epidemiological, geographical, eco-biological and social specificities in SPICs.

2-[(3,3-Dimethylindolin-2-ylidene)methyl]-4-[(3,3-dimethyl-3H-indol-1-ium-2-yl)methylidene]-3-oxocyclobut-1-en-1-olate chloroform disolvate

In the title squaraine dye solvate, C26H24N2O2·2CHCl3, the dye molecule is essentially planar, except for the methyl groups, having a maximum deviation over the 26-membered delocalized bond system of 0.060 (2) Å. It possesses crystallographic twofold rotational symmetry with the indole ring systems adopting a syn conformation. The molecular structure features intramolecular N-HO hydrogen bonds enclosing conjoint S7 ring motifs about one of the dioxocyclobutene O atoms, while the two chloroform solvent molecules are linked to the second O atom through C-HO hydrogen bonds.

Hydrogen bonding in two ammonium salts of 5-sulfosalicylic acid : ammonium 3-carboxy-4-hydroxybenzenesulfonate monohydrate and triammonium 3-carboxy-4-hydroxybenzenesulfonate 3-carboxylato-4-hydroxybenzenesulfonate

The structures of two ammonium salts of 3-carboxy-4-hydroxybenzenesulfonic acid (5-sulfosalicylic acid, 5-SSA) have been determined at 200 K. In the 1:1 hydrated salt, ammonium 3-carboxy-4-hydroxybenzenesulfonate monohydrate, NH4+·C7H5O6S-·H2O, (I), the 5-SSA- monoanions give two types of head-to-tail laterally linked cyclic hydrogen-bonding associations, both with graph-set R44(20). The first involves both carboxylic acid O-HOwater and water O-HOsulfonate hydrogen bonds at one end, and ammonium N-HOsulfonate and N-HOcarboxy hydrogen bonds at the other. The second association is centrosymmetric, with end linkages through water O-HOsulfonate hydrogen bonds. These conjoined units form stacks down c and are extended into a three-dimensional framework structure through N-HO and water O-HO hydrogen bonds to sulfonate O-atom acceptors. Anhydrous triammonium 3-carboxy-4-hydroxybenzenesulfonate 3-carboxylato-4-hydroxybenzenesulfonate, 3NH4+·C7H4O6S2-·C7H5O6S-, (II), is unusual, having both dianionic 5-SSA2- and monoanionic 5-SSA- species. These are linked by a carboxylic acid O-HO hydrogen bond and, together with the three ammonium cations (two on general sites and the third comprising two independent half-cations lying on crystallographic twofold rotation axes), give a pseudo-centrosymmetric asymmetric unit. Cation-anion hydrogen bonding within this layered unit involves a cyclic R33(8) association which, together with extensive peripheral N-HO hydrogen bonding involving both sulfonate and carboxy/carboxylate acceptors, gives a three-dimensional framework structure. This work further demonstrates the utility of the 5-SSA- monoanion for the generation of stable hydrogen-bonded crystalline materials, and provides the structure of a dianionic 5-SSA2- species of which there are only a few examples in the crystallographic literature.

Proton-transfer compounds with 4-amino-N-(4,6-dimethylpyrimidin-2-yl)benzenesulfonamide (sulfamethazine): the structures and hydrogen-bonding in the salts with 5-nitrosalicylic acid and picric acid

The structures of the anhydrous proton-transfer compounds of the sulfa drug sulfamethazine with 5-nitrosalicylic acid and picric acid, namely 2-(4-aminobenzenesulfonamido)-4,6-dimethylpyrimidinium 2-hydroxy-5-nitrobenzoate, C12H15N4O2S(+)·C7H4NO4(-), (I), and 2-(4-aminobenzenesulfonamido)-4,6-dimethylpyrimidinium 2,4,6-trinitrophenolate, C12H15N4O2S(+)·C6H2N3O7(-), (II), respectively, have been determined. In the asymmetric unit of (I), there are two independent but conformationally similar cation-anion heterodimer pairs which are formed through duplex intermolecular N(+)-H...Ocarboxylate and N-H...Ocarboxylate hydrogen-bond pairs, giving a cyclic motif [graph set R2(2)(8)]. These heterodimers form separate and different non-associated substructures through aniline N-H...O hydrogen bonds, one one-dimensional, involving carboxylate O-atom acceptors, the other two-dimensional, involving both carboxylate and hydroxy O-atom acceptors. The overall two-dimensional structure is stabilized by π-π interactions between the pyrimidinium ring and the 5-nitrosalicylate ring in both heterodimers [minimum ring-centroid separation = 3.4580 (8) Å]. For picrate (II), the cation-anion interaction involves a slightly asymmetric chelating N-H...O R2(1)(6) hydrogen-bonding association with the phenolate O atom, together with peripheral conjoint R1(2)(6) interactions between the same N-H groups and O atoms of the ortho-related nitro groups. An inter-unit amine N-H...Osulfone hydrogen bond gives one-dimensional chains which extend along a and inter-associate through π-π interactions between the pyrimidinium rings [centroid-centroid separation = 3.4752 (9) Å]. The two structures reported here now bring to a total of four the crystallographically characterized examples of proton-transfer salts of sulfamethazine with strong organic acids.

Stage-specific embryonic antigen-4 is not a marker for chondrogenic and osteogenic potential in cultured chondrocytes and mesenchymal progenitor cells

One important challenge for regenerative medicine is to produce a clinically relevant number of cells with consistent tissue-forming potential. Isolation and expansion of cells from skeletal tissues results in a heterogeneous population of cells with variable regenerative potential. A more consistent tissue formation could be achieved by identification and selection of potent progenitors based on cell surface molecules. In this study, we assessed the expression of stage-specific embryonic antigen-4 (SSEA-4), a classic marker of undifferentiated stem cells, and other surface markers in human articular chondrocytes (hACs), osteoblasts, and bone marrow-derived mesenchymal stromal cells (bmMSCs) and characterized their differentiation potential. Further, we sorted SSEA-4-expressing hACs and followed their potential to proliferate and to form cartilage in vitro. Cells isolated from cartilage and bone exhibited remarkably heterogeneous SSEA-4 expression profiles in expansion cultures. SSEA-4 expression levels increased up to approximately 5 population doublings, but decreased following further expansion and differentiation cultures; levels were not related to the proliferation state of the cells. Although SSEA-4-sorted chondrocytes showed a slightly better chondrogenic potential than their SSEA-4-negative counterparts, differences were insufficient to establish a link between SSEA-4 expression and chondrogenic potential. SSEA-4 levels in bmMSCs also did not correlate to the cells' chondrogenic and osteogenic potential in vitro. SSEA-4 is clearly expressed by subpopulations of proliferating somatic cells with a MSC-like phenotype. However, the predictive value of SSEA-4 as a specific marker of superior differentiation capacity in progenitor cell populations from adult human tissue and even its usefulness as a stem cell marker appears questionable.

Structural and electronic properties of diazonium functionalized (4, 4) single walled carbon nanotube : an ab initio study

In this work, ab initio density functional theory (DFT) calculations are performed to study the structural and electronic properties of diazonium reagent functionalized (4, 4) single-walled carbon nanotube (SWCNT). We find the aryl group covalently bonds with SWCNT and prefers to be perpendicular to the side wall of nanotube. It has a rotational barrier of 0.35 eV around the formed aryl-tube bond axis and should be thermodynamically stable at room temperature. Additionally, new peaks appeared around the Fermi energy in the density of state (DOS) due to the weak band dispersion. Increasing of the coverage of the functional group will result in significant upshift of the Fermi level.

Phantom Scatter and the A-SI EPID

Introduction: The use of amorphous-silicon electronic portal imaging devices (a-Si EPIDs) for dosimetry is complicated by the effects of scattered radiation. In photon radiotherapy, primary signal at the detector can be accompanied by photons scattered from linear accelerator components, detector materials, intervening air, treatment room surfaces (floor, walls, etc) and from the patient/phantom being irradiated. Consequently, EPID measurements which presume to take scatter into account are highly sensitive to the identification of these contributions. One example of this susceptibility is the process of calibrating an EPID for use as a gauge of (radiological) thickness, where specific allowance must be made for the effect of phantom-scatter on the intensity of radiation measured through different thicknesses of phantom. This is usually done via a theoretical calculation which assumes that phantom scatter is linearly related to thickness and field-size. We have, however, undertaken a more detailed study of the scattering effects of fields of different dimensions when applied to phantoms of various thicknesses in order to derive scattered-primary ratios (SPRs) directly from simulation results. This allows us to make a more-accurate calibration of the EPID, and to qualify the appositeness of the theoretical SPR calculations. Methods: This study uses a full MC model of the entire linac-phantom-detector system simulated using EGSnrc/BEAMnrc codes. The Elekta linac and EPID are modelled according to specifications from the manufacturer and the intervening phantoms are modelled as rectilinear blocks of water or plastic, with their densities set to a range of physically realistic and unrealistic values. Transmissions through these various phantoms are calculated using the dose detected in the model EPID and used in an evaluation of the field-size-dependence of SPR, in different media, applying a method suggested for experimental systems by Swindell and Evans [1]. These results are compared firstly with SPRs calculated using the theoretical, linear relationship between SPR and irradiated volume, and secondly with SPRs evaluated from our own experimental data. An alternate evaluation of the SPR in each simulated system is also made by modifying the BEAMnrc user code READPHSP, to identify and count those particles in a given plane of the system that have undergone a scattering event. In addition to these simulations, which are designed to closely replicate the experimental setup, we also used MC models to examine the effects of varying the setup in experimentally challenging ways (changing the size of the air gap between the phantom and the EPID, changing the longitudinal position of the EPID itself). Experimental measurements used in this study were made using an Elekta Precise linear accelerator, operating at 6MV, with an Elekta iView GT a-Si EPID. Results and Discussion: 1. Comparison with theory: With the Elekta iView EPID fixed at 160 cm from the photon source, the phantoms, when positioned isocentrically, are located 41 to 55 cm from the surface of the panel. At this geometry, a close but imperfect agreement (differing by up to 5%) can be identified between the results of the simulations and the theoretical calculations. However, this agreement can be totally disrupted by shifting the phantom out of the isocentric position. Evidently, the allowance made for source-phantom-detector geometry by the theoretical expression for SPR is inadequate to describe the effect that phantom proximity can have on measurements made using an (infamously low-energy sensitive) a-Si EPID. 2. Comparison with experiment: For various square field sizes and across the range of phantom thicknesses, there is good agreement between simulation data and experimental measurements of the transmissions and the derived values of the primary intensities. However, the values of SPR obtained through these simulations and measurements seem to be much more sensitive to slight differences between the simulated and real systems, leading to difficulties in producing a simulated system which adequately replicates the experimental data. (For instance, small changes to simulated phantom density make large differences to resulting SPR.) 3. Comparison with direct calculation: By developing a method for directly counting the number scattered particles reaching the detector after passing through the various isocentric phantom thicknesses, we show that the experimental method discussed above is providing a good measure of the actual degree of scattering produced by the phantom. This calculation also permits the analysis of the scattering sources/sinks within the linac and EPID, as well as the phantom and intervening air. Conclusions: This work challenges the assumption that scatter to and within an EPID can be accounted for using a simple, linear model. Simulations discussed here are intended to contribute to a fuller understanding of the contribution of scattered radiation to the EPID images that are used in dosimetry calculations. Acknowledgements: This work is funded by the NHMRC, through a project grant, and supported by the Queensland University of Technology (QUT) and the Royal Brisbane and Women's Hospital, Brisbane, Australia. The authors are also grateful to Elekta for the provision of manufacturing specifications which permitted the detailed simulation of their linear accelerators and amorphous-silicon electronic portal imaging devices. Computational resources and services used in this work were provided by the HPC and Research Support Group, QUT, Brisbane, Australia.

Vibrational spectroscopic characterization of the phosphate mineral bermanite–Mn2+Mn23+(PO4)2(OH)2.4(H2O)

Bermanite Mn2þMn3þ2 ðPO4Þ2ðOHÞ2 � 4ðH2OÞ is a mixed valent hydrated hydroxy phosphate mineral. The mineral is reddish-brown and occurs in crystal aggregates and as lamellar masses. Bermanite is a common mineral in granitic pegmatites. The chemical composition of bermanite was obtained using EDS techniques. We have studied the molecular structure of bermanite using vibrational spectroscopy. The mineral is characterized by a Raman doublet at 991 and 999 cm-1 attributed to the phosphate stretching mode of two non-equivalent phosphate units. Raman bands at 1071, 1117 and 1142 cm-1 are assigned to the phosphate antisymmetric stretching modes. The hydroxyl stretching spectral region is complex with overlapping bands attributed to water and hydroxyl stretching vibrations. Vibrational spectroscopy proves most useful for the study of the mineral bermanite.

Vibrational spectroscopy of the borate mineral olshanskyite Ca3[B(OH)4]4(OH)2

The mineral olshanskyite is one of many calcium borate minerals which has never been studied using vibrational spectroscopy. The mineral is unstable and decomposes upon exposure to an electron beam. This makes the elemental analysis using EDX techniques difficult. Both the Raman and infrared spectra show complexity due to the complexity of the structure. Intense Raman bands are found at 989, 1,003, 1,025 and 1,069 cm-1 with a shoulder at 961 cm-1 and are assigned to trigonal borate units. The Raman bands at 1,141, 1,206 and 1,365 cm-1 are assigned to OH in-plane bending of BOH units. A series of Raman bands are observed in the 2,900–3,621cm-1 spectral range and are assigned to the stretching vibrations of OH and water. This complexity is also reflected in the infrared spectra. Vibrational spectroscopy enables aspects of the structure of olshanskyite to be elucidated.

Molecular cocrystals of 5-(4-bromophenyl)-1,3,4-thiadiazol-2-amine: hydrogen bonding in the structures of the 1:1 adduct with 2-(naphthalen-2-yloxy)acetic acid and the salt with 3,5-dinitrobenzoic acid

The structures of the anhydrous products from the interaction of 2-amino-5-(4-bromophenyl)-1,3,4-thiadiazole with (2-naphthoxy)acetic acid, the 1:1 adduct C8H6BrN3S . C12H10O3 (I) and 3,5-dinitrobenzoic acid, the salt C8H7BrN3S+ C7H3N2O6- (II) have been determined. In the adduct (I), a heterodimer is formed through a cyclic hydrogen-bonding motif [graph set R2/2(8)], involving carboxylic acid O-H...N(hetero)and amine N-H...O(carboxyl) interactions. The heterodimers are essentially planar with a thiadiazole to naphthyl ring dihedral angle of 15.9(2)deg. and the intramolecular thiadiazole to phenyl ring angle of 4.7(2)deg. An amine N-H...N(hetero) hydrogen bond between the heterodimers generates a one-dimensional chain structure extending down [001]. Also present are weak benzene-benzene and naphthalene-naphthalene pi-pi stacking interactions down the b axis [minimum ring centroid separation, 3.936(3) Ang.]. With the salt (II), the cation-anion association is also through a cyclic R2/2(8) motif but involving duplex N-H...O(carboxyl) hydrogen bonds, giving a heterodimer which is close to planar [dihedral angles between the thiadiazole ring and the two benzene rings, 5.00(16)deg. (intra) and 7.23(15)deg. (inter)]. A secondary centrosymmetric cyclic N-H...O(carboxyl) hydrogen-bonding association involving the second amino H-atom generates a heterotetramer. Also present in the crystal are weak pi-pi i-\p interactions between thiadiazolium rings [minimum ring centroid separation, 3.936(3)Ang.], as well as a short Br...O(nitro) interaction [3.314(4)Ang.]. The two structures reported here now provide a total of three crystallographically characterized examples of co-crystalline products from the interaction of 2-amino-5-(4-bromophenyl)-1,3,4-thiadiazole with carboxylic acids, of which only one involves proton-transfer.

Assessment of the molecular structure of natrodufrénite – NaFe2+Fe3+5(PO4)4(OH)6•2(H2O), a secondary pegmatite phosphate mineral from Minas Gerais, Brazil

The mineral natrodufrénite a secondary pegmatite phosphate mineral from Minas Gerais, Brazil, has been studied by a combination of scanning electron microscopy and vibrational spectroscopic techniques. Electron probe analysis shows the formula of the studied mineral as (Na0.88Ca0.12)∑1.00(Mn0.11Mg0.08Ca0.04Zr0.01Cu0.01)∑0.97(Al0.02)∑4.91(PO4)3.96(OH6.15F0.07)6.22⋅2.05(H2O). Raman spectroscopy identifies an intense peak at 1003 cm−1 assigned to the ν1 symmetric stretching mode. Raman bands are observed at 1059 and 1118 cm−1 and are attributed to the ν3 antisymmetric stretching vibrations. A comparison is made with the spectral data of other hydrate hydroxy phosphate minerals including cyrilovite and wardite. Raman bands at 560, 582, 619 and 668 cm−1 are assigned to the ν4 bending modes and Raman bands at 425, 444, 477 and 507 cm−1 are due to the ν2 bending modes. Raman bands in the 2600–3800 cm−1 spectral range are attributed to water and OH stretching vibrations. Vibrational spectroscopy enables aspects of the molecular structure of natrodufrénite to be assessed.

Investigation between the S377G3 GATA-4 polymorphism and migraine

Migraine is a common and painful neurological disorder, with genetic and environmental components. Several conditions have been shown to be comorbid with migraine, notably a cardiac malformation affecting the interatrial septum and leading to patent foramen ovale (PFO). Mutations in the development regulatory gene GATA-4, located on human chromosome 8p23.1-p22, have been found to be responsible for some cases of congenital heart defects including PFO. To determine whether the GATA-4 gene is involved in migraine, the present study performed an association analysis of a common GATA-4 variant that results in a change of amino acid (S377G), in a large case/control population (275 unrelated Caucasian migraineurs versus 275 control individuals). The results showed that there was no significant association for this polymorphism between migraine and controls (χ² = 0.84, P = 0.66). Thus it appears that the GATA-4 (S377G) mutation does not play a significant role in common migraine susceptibility.

No role for estrogen receptor 1 gene intron 1 Pvu II and exon 4 C325G polymorphisms in migraine susceptibility

Background We have previously reported an association between the estrogen receptor 1 (ESR1) gene exon 8 G594A polymorphism and migraine susceptibility in two independent Australian cohorts. In this paper we report results of analysis of two further single nucleotide polymorphisms (SNPs) in the ESR1 gene in the same study group, the T/C Pvu II SNP in intron 1 and the C325G SNP in exon 4, as well as results of linkage disequilibrium (LD) analysis on these markers. Methods We investigated these variants by case-control association analysis in a cohort of 240 migraineurs and 240 matched controls. The SNPs were genotyped using specific restriction enzyme assays. Results were analysed using contingency table methods incorporating the chi-squared statistic. LD results are presented as D' statistics with associated P values. Results We found no evidence for association of the Pvu II T/C polymorphism and the C325G polymorphism and migraine susceptibility and no evidence for LD between these two SNPs and the previously implicated exon 8 G594A marker. Conclusion We have found no role for the polymorphisms in intron 1 and exon 4 with migraine susceptibility. To further investigate our previously implicated exon 8 marker, we suggest the need for studies with a high density of polymorphisms be undertaken, with particular focus on markers in LD with the exon 8 marker.