Inhibition of dendritic Ca2+ spikes by GABAB receptors in cortical pyramidal neurons is mediated by a direct Gi/o- -subunit interaction with Cav1 channels

Voltage-dependent calcium channels (VDCCs) serve a wide range of physiological functions and their activity is modulated by different neurotransmitter systems. GABAergic inhibition of VDCCs in neurons has an important impact in controlling transmitter release, neuronal plasticity, gene expression and neuronal excitability. We investigated the molecular signalling mechanisms by which GABAB receptors inhibit calcium-mediated electrogenesis (Ca2+ spikes) in the distal apical dendrite of cortical layer 5 pyramidal neurons. Ca2+ spikes are the basis of coincidence detection and signal amplification of distal tuft synaptic inputs characteristic for the computational function of cortical pyramidal neurons. By combining dendritic whole-cell recordings with two-photon fluorescence Ca2+ imaging we found that all subtypes of VDCCs were present in the Ca2+ spike initiation zone, but that they contribute differently to the initiation and sustaining of dendritic Ca2+ spikes. Particularly, Cav1 VDCCs are the most abundant VDCC present in this dendritic compartment and they generated the sustained plateau potential characteristic for the Ca2+ spike. Activation of GABAB receptors specifically inhibited Cav1 channels. This inhibition of L-type Ca2+ currents was transiently relieved by strong depolarization but did not depend on protein kinase activity. Therefore, our findings suggest a novel membrane-delimited interaction of the Gi/o-βγ-subunit with Cav1 channels identifying this mechanism as the general pathway of GABAB receptor-mediated inhibition of VDCCs. Furthermore, the characterization of the contribution of the different VDCCs to the generation of the Ca2+ spike provides new insights into the molecular mechanism of dendritic computation.

Low-lying excited states and nonradiative processes of 9-methyl-2-aminopurine

he UV spectrum of the adenine analogue 9-methyl-2-aminopurine (9M-2AP) is investigated with one- and two-color resonant two-photon ionization spectroscopy at 0.3 and 0.05 cm−1 resolution in a supersonic jet. The electronic origin at 32 252 cm−1 exhibits methyl torsional subbands that originate from the 0A′′1 (l = 0) and 1E ″ (l = ±1) torsional levels. These and further torsional bands that appear up to 000+230 cm−1 allow to fit the threefold (V 3) barriers of the torsional potentials as ∣∣V′′3∣∣=50 cm−1 in the S 0 and ∣∣V′3∣∣=126 cm−1 in the S 1 state. Using the B3LYP density functional and correlated approximate second-order coupled cluster CC2 methods, the methyl orientation is calculated to be symmetric relative to the 2AP plane in both states, with barriers of V′′3=20 cm−1 and V′3=115 cm−1. The 000 rotational band contour is 75% in-plane (a/b) polarized, characteristic for a dominantly long-axis 1ππ* excitation. The residual 25% c-axis polarization may indicate coupling of the 1ππ* to the close-lying 1 nπ* state, calculated at 4.00 and 4.01 eV with the CC2 method. However, the CC2 calculated 1 nπ oscillator strength is only 6% of that of the 1ππ* transition. The 1ππ* vibronic spectrum is very complex, showing about 40 bands within the lowest 500 cm−1. The methyl torsion and the low-frequency out-of-plane ν′1 and ν′2 vibrations are strongly coupled in the 1ππ* state. This gives rise to many torsion-vibration combination bands built on out-of-plane fundamentals, which are without precedence in the 1ππ* spectrum of 9H-2-aminopurine [S. Lobsiger, R. K. Sinha, M. Trachsel, and S. Leutwyler, J. Chem. Phys.134, 114307 (2011)]. From the Lorentzian broadening needed to fit the 000 contour of 9M-2AP, the 1ππ* lifetime is τ ⩾ 120 ps, reflecting a rapid nonradiative transition.

Watson–Crick and Sugar-Edge Base Pairing of Cytosine in the Gas Phase: UV and Infrared Spectra of Cytosine·2-Pyridone

While keto-amino cytosine is the dominant species in aqueous solution, spectroscopic studies in molecular beams and in noble gas matrices show that other cytosine tautomers prevail in apolar environments. Each of these offers two or three H-bonding sites (Watson–Crick, wobble, sugar-edge). The mass- and isomer-specific S1 ← S0 vibronic spectra of cytosine·2-pyridone (Cyt·2PY) and 1-methylcytosine·2PY are measured using UV laser resonant two-photon ionization (R2PI), UV/UV depletion, and IR depletion spectroscopy. The UV spectra of the Watson–Crick and sugar-edge isomers of Cyt·2PY are separated using UV/UV spectral hole-burning. Five different isomers of Cyt·2PY are observed in a supersonic beam. We show that the Watson–Crick and sugar-edge dimers of keto-amino cytosine with 2PY are the most abundant in the beam, although keto-amino-cytosine is only the third most abundant tautomer in the gas phase. We identify the different isomers by combining three different diagnostic tools: (1) methylation of the cytosine N1–H group prevents formation of both the sugar-edge and wobble isomers and gives the Watson–Crick isomer exclusively. (2) The calculated ground state binding and dissociation energies, relative gas-phase abundances, excitation and the ionization energies are in agreement with the assignment of the dominant Cyt·2PY isomers to the Watson–Crick and sugar-edge complexes of keto-amino cytosine. (3) The comparison of calculated ground state vibrational frequencies to the experimental IR spectra in the carbonyl stretch and NH/OH/CH stretch ranges strengthen this identification.

Naive B-cell trafficking is shaped by local chemokine availability and LFA-1-independent stromal interactions.

It is not known how naive B cells compute divergent chemoattractant signals of the T-cell area and B-cell follicles during in vivo migration. Here, we used two-photon microscopy of peripheral lymph nodes (PLNs) to analyze the prototype G-protein-coupled receptors (GPCRs) CXCR4, CXCR5, and CCR7 during B-cell migration, as well as the integrin LFA-1 for stromal guidance. CXCR4 and CCR7 did not influence parenchymal B-cell motility and distribution, despite their role during B-cell arrest in venules. In contrast, CXCR5 played a nonredundant role in B-cell motility in follicles and in the T-cell area. B-cell migration in the T-cell area followed a random guided walk model, arguing against directed migration in vivo. LFA-1, but not α4 integrins, contributed to B-cell motility in PLNs. However, stromal network guidance was LFA-1 independent, uncoupling integrin-dependent migration from stromal attachment. Finally, we observed that despite a 20-fold reduction of chemokine expression in virus-challenged PLNs, CXCR5 remained essential for B-cell screening of antigen-presenting cells. Our data provide an overview of the contribution of prototype GPCRs and integrins during naive B-cell migration and shed light on the local chemokine availability that these cells compute.

Development of light-responsive porous polycarbonate membranes for controlled caffeine delivery

For controlled caffeine release, light-responsive membranes were developed. It was possible to produce membranes that reduced their caffeine permeability resistance by about 97% when irradiated with UV-light compared to measurements at daylight. This was achieved by grafting polymers possessing photochromic units onto track-edged polycarbonate membranes. Covalently linked coatings on porous polycarbonate membranes were obtained by plasma activation of the membrane surface followed by plasma-induced graft polymerization. Copolymerization of spiro-compounds during the coating process as well as postmodification of preformed coatings with spiropyran resulted in photochromic membranes. For the copolymerization process, the synthesis of five photochromic methacrylic and acrylic spiropyrans and spirooxazines was successfully performed. Additionally, a spiropyran with carboxylic acid functionality was synthesized for the postmodification process. This enabled us to postmodify polymeric materials containing alcohol or amine groups to obtain photochromic materials. UV-irradiation of these light-responsive membranes resulted in a strong colouration of the membrane, in a reduction of surface tension, which resulted in a decreased caffeine permeability resistance. The membranes were characterized using XPS for the elemental composition of the coating, contact angle measurements for the surface tension, solid-state UV/VIS measurements for the determination of the kinetic and stability properties, and two-photon microscopy for the localisation of the photochromic substance in the porous membrane.

Excitonic Splitting, Delocalization, and Vibronic Quenching in the Benzonitrile Dimer

The excitonic S1/S2 state splitting and the localization/delocalization of the S1 and S2 electronic states are investigated in the benzonitrile dimer (BN)2 and its 13C and d5 isotopomers by mass-resolved two-color resonant two-photon ionization spectroscopy in a supersonic jet, complemented by calculations. The doubly hydrogen-bonded (BN-h5)2 and (BN-d5)2 dimers are C2h symmetric with equivalent BN moieties. Only the S0 → S2 electronic origin is observed, while the S0 → S1 excitonic component is electric-dipole forbidden. A single 12C/13C or 5-fold h5/d5 isotopic substitution reduce the dimer symmetry to Cs, so that the heteroisotopic dimers (BN)2-(h5 – h513C), (BN)2-(h5 – d5), and (BN)2-(h5 – h513C) exhibit both S0 → S1 and S0 → S2 origins. Isotope-dependent contributions Δiso to the excitonic splittings arise from the changes of the BN monomer zero-point vibrational energies; these range from Δiso(12C/13C) = 3.3 cm–1 to Δiso(h5/d5) = 155.6 cm–1. The analysis of the experimental S1/S2 splittings of six different isotopomeric dimers yields the S1/S2 exciton splitting Δexc = 2.1 ± 0.1 cm–1. Since Δiso(h5/d5) ≫ Δexc and Δiso(12C/13C) > Δexc, complete and near-complete exciton localization occurs upon 12C/13C and h5/d5 substitutions, respectively, as diagnosed by the relative S0 → S1 and S0 → S2 origin band intensities. The S1/S2 electronic energy gap of (BN)2 calculated by the spin-component scaled approximate second-order coupled-cluster (SCS-CC2) method is Δelcalc = 10 cm–1. This electronic splitting is reduced by the vibronic quenching factor Γ. The vibronically quenched exciton splitting Δelcalc·Γ = Δvibroncalc = 2.13 cm–1 is in excellent agreement with the observed splitting Δexc = 2.1 cm–1. The excitonic splittings can be converted to semiclassical exciton hopping times; the shortest hopping time is 8 ps for the homodimer (BN-h5)2, the longest is 600 ps for the (BN)2(h5 – d5) heterodimer.

Excited-State Structure, Vibrations, and Nonradiative Relaxation of Jet-Cooled 5-Fluorocytosine

The S0 → S1 vibronic spectrum and S1 state nonradiative relaxation of jet-cooled keto-amino 5-fluorocytosine (5FCyt) are investigated by two-color resonant two-photon ionization spectroscopy at 0.3 and 0.05 cm–1 resolution. The 000 rotational band contour is polarized in-plane, implying that the electronic transition is 1ππ*. The electronic transition dipole moment orientation and the changes of rotational constants agree closely with the SCS-CC2 calculated values for the 1ππ* (S1) transition of 5FCyt. The spectral region from 0 to 300 cm–1 is dominated by overtone and combination bands of the out-of-plane ν1′ (boat), ν2′ (butterfly), and ν3′ (HN–C6H twist) vibrations, implying that the pyrimidinone frame is distorted out-of-plane by the 1ππ* excitation, in agreement with SCS-CC2 calculations. The number of vibronic bands rises strongly around +350 cm–1; this is attributed to the 1ππ* state barrier to planarity that corresponds to the central maximum of the double-minimum out-of-plane vibrational potentials along the ν1′, ν2′, and ν3′ coordinates, which gives rise to a high density of vibronic excitations. At +1200 cm–1, rapid nonradiative relaxation (knr ≥ 1012 s–1) sets in, which we interpret as the height of the 1ππ* state barrier in front of the lowest S1/S0 conical intersection. This barrier in 5FCyt is 3 times higher than that in cytosine. The lifetimes of the ν′ = 0, 2ν1′, 2ν2′, 2ν1′ + 2ν2′, 4ν2′, and 2ν1′ + 4ν2′ levels are determined from Lorentzian widths fitted to the rotational band contours and are τ ≥ 75 ps for ν′ = 0, decreasing to τ ≥ 55 ps at the 2ν1′ + 4ν2′ level at +234 cm–1. These gas-phase lifetimes are twice those of S1 state cytosine and 10–100 times those of the other canonical nucleobases in the gas phase. On the other hand, the 5FCyt gas-phase lifetime is close to the 73 ps lifetime in room-temperature solvents. This lack of dependence on temperature and on the surrounding medium implies that the 5FCyt nonradiative relaxation from its S1 (1ππ*) state is essentially controlled by the same ∼1200 cm–1 barrier and conical intersection both in the gas phase and in solution.

Measurement of Higgs boson production in the diphoton decay channel in pp collisions at center-of-mass energies of 7 and 8 TeV with the ATLAS detector

A measurement of the production processes of the recently discovered Higgs boson is performed in the two-photon final state using 4.5  fb −1 of proton-proton collisions data at s √ =7  TeV and 20.3  fb −1 at s √ =8  TeV collected by the ATLAS detector at the Large Hadron Collider. The number of observed Higgs boson decays to diphotons divided by the corresponding Standard Model prediction, called the signal strength, is found to be μ=1.17±0.27 at the value of the Higgs boson mass measured by ATLAS, m H =125.4  GeV . The analysis is optimized to measure the signal strengths for individual Higgs boson production processes at this value of m H . They are found to be μ ggF =1.32±0.38 , μ VBF =0.8±0.7 , μ WH =1.0±1.6 , μ ZH =0.1 +3.7 −0.1 , and μ tt ¯ H =1.6 +2.7 −1.8 , for Higgs boson production through gluon fusion, vector-boson fusion, and in association with a W or Z boson or a top-quark pair, respectively. Compared with the previously published ATLAS analysis, the results reported here also benefit from a new energy calibration procedure for photons and the subsequent reduction of the systematic uncertainty on the diphoton mass resolution. No significant deviations from the predictions of the Standard Model are found.

Experimental and Calculated Spectra of π-Stacked Mild Charge-Transfer Complexes: Jet-Cooled Perylene·(Tetrachloroethene) n , n = 1,2

The S0 ↔ S1 spectra of the mild charge-transfer (CT) complexes perylene·tetrachloroethene (P·4ClE) and perylene·(tetrachloroethene)2 (P·(4ClE)2) are investigated by two-color resonant two-photon ionization (2C-R2PI) and dispersed fluorescence spectroscopy in supersonic jets. The S0 → S1 vibrationless transitions of P·4ClE and P·(4ClE)2 are shifted by δν = −451 and −858 cm–1 relative to perylene, translating to excited-state dissociation energy increases of 5.4 and 10.3 kJ/mol, respectively. The red shift is ∼30% larger than that of perylene·trans-1,2-dichloroethene; therefore, the increase in chlorination increases the excited-state stabilization and CT character of the interaction, but the electronic excitation remains largely confined to the perylene moiety. The 2C-R2PI and fluorescence spectra of P·4ClE exhibit strong progressions in the perylene intramolecular twist (1au) vibration (42 cm–1 in S0 and 55 cm–1 in S1), signaling that perylene deforms along its twist coordinate upon electronic excitation. The intermolecular stretching (Tz) and internal rotation (Rc) vibrations are weak; therefore, the P·4ClE intermolecular potential energy surface (IPES) changes little during the S0 ↔ S1 transition. The minimum-energy structures and inter- and intramolecular vibrational frequencies of P·4ClE and P·(4ClE)2 are calculated with the dispersion-corrected density functional theory (DFT) methods B97-D3, ωB97X-D, M06, and M06-2X and the spin-consistent-scaled (SCS) variant of the approximate second-order coupled-cluster method, SCS-CC2. All methods predict the global minima to be π-stacked centered coplanar structures with the long axis of tetrachloroethene rotated by τ ≈ 60° relative to the perylene long axis. The calculated binding energies are in the range of −D0 = 28–35 kJ/mol. A second minimum is predicted with τ ≈ 25°, with ∼1 kJ/mol smaller binding energy. Although both monomers are achiral, both the P·4ClE and P·(4ClE)2 complexes are chiral. The best agreement for adiabatic excitation energies and vibrational frequencies is observed for the ωB97X-D and M06-2X DFT methods.

The elusive S2 state, the S1/S2 splitting, and the excimer states of the benzene dimer

We observe the weak S 0 → S 2 transitions of the T-shaped benzene dimers (Bz)2 and (Bz-d 6)2 about 250 cm−1 and 220 cm−1 above their respective S 0 → S 1 electronic origins using two-color resonant two-photon ionization spectroscopy. Spin-component scaled (SCS) second-order approximate coupled-cluster (CC2) calculations predict that for the tipped T-shaped geometry, the S 0 → S 2 electronic oscillator strength f el (S 2) is ∼10 times smaller than f el (S 1) and the S 2 state lies ∼240 cm−1 above S 1, in excellent agreement with experiment. The S 0 → S 1 (ππ ∗) transition is mainly localized on the “stem” benzene, with a minor stem → cap charge-transfer contribution; the S 0 → S 2 transition is mainly localized on the “cap” benzene. The orbitals, electronic oscillator strengths f el (S 1) and f el (S 2), and transition frequencies depend strongly on the tipping angle ω between the two Bz moieties. The SCS-CC2 calculated S 1 and S 2 excitation energies at different T-shaped, stacked-parallel and parallel-displaced stationary points of the (Bz)2 ground-state surface allow to construct approximate S 1 and S 2 potential energy surfaces and reveal their relation to the “excimer” states at the stacked-parallel geometry. The f el (S 1) and f el (S 2) transition dipole moments at the C 2v -symmetric T-shape, parallel-displaced and stacked-parallel geometries are either zero or ∼10 times smaller than at the tipped T-shaped geometry. This unusual property of the S 0 → S 1 and S 0 → S 2 transition-dipole moment surfaces of (Bz)2 restricts its observation by electronic spectroscopy to the tipped and tilted T-shaped geometries; the other ground-state geometries are impossible or extremely difficult to observe. The S 0 → S 1/S 2 spectra of (Bz)2 are compared to those of imidazole ⋅ (Bz)2, which has a rigid triangular structure with a tilted (Bz)2 subunit. The S 0 → S 1/ S 2 transitions of imidazole-(benzene)2 lie at similar energies as those of (Bz)2, confirming our assignment of the (Bz)2 S 0 → S 2 transition.

Excitonic splittings in molecular dimers: why static ab initio calculations cannot match them

After decades of research on molecular excitons, only few molecular dimers are available on which exciton and vibronic coupling theories can be rigorously tested. In centrosymmetric H-bonded dimers consisting of identical (hetero)aromatic chromophores, the monomer electronic transition dipole moment vectors subtract or add, yielding S0 → S1 and S0 → S2 transitions that are symmetry-forbidden or -allowed, respectively. Symmetry breaking by 12C/13C or H/D isotopic substitution renders the forbidden transition weakly allowed. The excitonic coupling (Davydov splitting) can then be measured between the S0 → S1 and S0 → S2 vibrationless bands. We discuss the mass-specific excitonic spectra of five H-bonded dimers that are supersonically cooled to a few K and investigated using two-color resonant two-photon ionization spectroscopy. The excitonic splittings Δcalc predicted by ab initio methods are 5–25 times larger than the experimental excitonic splittings Δexp. The purely electronic ab initio splittings need to be reduced (“quenched”), reflecting the coupling of the electronic transition to the optically active vibrations of the monomers. The so-called quenching factors Γ < 1 can be determined from experiment (Γexp) and/or calculation (Γcalc). The vibronically quenched splittings Γ·Δcalc are found to nicely reproduce the experimental exciton splittings.

Intravital and whole-organ imaging reveals capture of melanoma-derived antigen by lymph node subcapsular macrophages leading to widespread deposition on follicular dendritic cells.

Aberrant antigens expressed by tumor cells, such as in melanoma, are often associated with humoral immune responses, which may in turn influence tumor progression. Despite recent data showing the central role of adaptive immune responses on cancer spread or control, it remains poorly understood where and how tumor-derived antigen (TDA) induces a humoral immune response in tumor-bearing hosts. Based on our observation of TDA accumulation in B cell areas of lymph nodes (LNs) from melanoma patients, we developed a pre-metastatic B16.F10 melanoma model expressing a fluorescent fusion protein, tandem dimer tomato, as a surrogate TDA. Using intravital two-photon microscopy (2PM) and whole-mount 3D LN imaging of tumor-draining LNs in immunocompetent mice, we report an unexpectedly widespread accumulation of TDA on follicular dendritic cells (FDCs), which were dynamically scanned by circulating B cells. Furthermore, 2PM imaging identified macrophages located in the subcapsular sinus of tumor-draining LNs to capture subcellular TDA-containing particles arriving in afferent lymph. As a consequence, depletion of macrophages or genetic ablation of B cells and FDCs resulted in dramatically reduced TDA capture in tumor-draining LNs. In sum, we identified a major pathway for the induction of humoral responses in a melanoma model, which may be exploitable to manipulate anti-TDA antibody production during cancer immunotherapy.

An efficient framework for photon Monte Carlo treatment planning

Currently photon Monte Carlo treatment planning (MCTP) for a patient stored in the patient database of a treatment planning system (TPS) can usually only be performed using a cumbersome multi-step procedure where many user interactions are needed. This means automation is needed for usage in clinical routine. In addition, because of the long computing time in MCTP, optimization of the MC calculations is essential. For these purposes a new graphical user interface (GUI)-based photon MC environment has been developed resulting in a very flexible framework. By this means appropriate MC transport methods are assigned to different geometric regions by still benefiting from the features included in the TPS. In order to provide a flexible MC environment, the MC particle transport has been divided into different parts: the source, beam modifiers and the patient. The source part includes the phase-space source, source models and full MC transport through the treatment head. The beam modifier part consists of one module for each beam modifier. To simulate the radiation transport through each individual beam modifier, one out of three full MC transport codes can be selected independently. Additionally, for each beam modifier a simple or an exact geometry can be chosen. Thereby, different complexity levels of radiation transport are applied during the simulation. For the patient dose calculation, two different MC codes are available. A special plug-in in Eclipse providing all necessary information by means of Dicom streams was used to start the developed MC GUI. The implementation of this framework separates the MC transport from the geometry and the modules pass the particles in memory; hence, no files are used as the interface. The implementation is realized for 6 and 15 MV beams of a Varian Clinac 2300 C/D. Several applications demonstrate the usefulness of the framework. Apart from applications dealing with the beam modifiers, two patient cases are shown. Thereby, comparisons are performed between MC calculated dose distributions and those calculated by a pencil beam or the AAA algorithm. Interfacing this flexible and efficient MC environment with Eclipse allows a widespread use for all kinds of investigations from timing and benchmarking studies to clinical patient studies. Additionally, it is possible to add modules keeping the system highly flexible and efficient.

Testing the GlaaS algorithm for dose measurements on low- and high-energy photon beams using an amorphous silicon portal imager

The GLAaS algorithm for pretreatment intensity modulation radiation therapy absolute dose verification based on the use of amorphous silicon detectors, as described in Nicolini et al. [G. Nicolini, A. Fogliata, E. Vanetti, A. Clivio, and L. Cozzi, Med. Phys. 33, 2839-2851 (2006)], was tested under a variety of experimental conditions to investigate its robustness, the possibility of using it in different clinics and its performance. GLAaS was therefore tested on a low-energy Varian Clinac (6 MV) equipped with an amorphous silicon Portal Vision PV-aS500 with electronic readout IAS2 and on a high-energy Clinac (6 and 15 MV) equipped with a PV-aS1000 and IAS3 electronics. Tests were performed for three calibration conditions: A: adding buildup on the top of the cassette such that SDD-SSD = d(max) and comparing measurements with corresponding doses computed at d(max), B: without adding any buildup on the top of the cassette and considering only the intrinsic water-equivalent thickness of the electronic portal imaging devices device (0.8 cm), and C: without adding any buildup on the top of the cassette but comparing measurements against doses computed at d(max). This procedure is similar to that usually applied when in vivo dosimetry is performed with solid state diodes without sufficient buildup material. Quantitatively, the gamma index (gamma), as described by Low et al. [D. A. Low, W. B. Harms, S. Mutic, and J. A. Purdy, Med. Phys. 25, 656-660 (1998)], was assessed. The gamma index was computed for a distance to agreement (DTA) of 3 mm. The dose difference deltaD was considered as 2%, 3%, and 4%. As a measure of the quality of results, the fraction of field area with gamma larger than 1 (%FA) was scored. Results over a set of 50 test samples (including fields from head and neck, breast, prostate, anal canal, and brain cases) and from the long-term routine usage, demonstrated the robustness and stability of GLAaS. In general, the mean values of %FA remain below 3% for deltaD equal or larger than 3%, while they are slightly larger for deltaD = 2% with %FA in the range from 3% to 8%. Since its introduction in routine practice, 1453 fields have been verified with GLAaS at the authors' institute (6 MV beam). Using a DTA of 3 mm and a deltaD of 4% the authors obtained %FA = 0.9 +/- 1.1 for the entire data set while, stratifying according to the dose calculation algorithm, they observed: %FA = 0.7 +/- 0.9 for fields computed with the analytical anisotropic algorithm and %FA = 2.4 +/- 1.3 for pencil-beam based fields with a statistically significant difference between the two groups. If data are stratified according to field splitting, they observed %FA = 0.8 +/- 1.0 for split fields and 1.0 +/- 1.2 for nonsplit fields without any significant difference.