The use of two-way linear mixed models in multitreatment meta-analysis.

Meta-analysis summarizes the results of a series of trials. When more than two treatments are included in the trials and when the set of treatments tested differs between trials, the combination of results across trials requires some care. Several methods have been proposed for this purpose, which feature under different labels, such as network meta-analysis or mixed treatment comparisons. Two types of linear mixed model can be used for meta-analysis. The one expresses the expected outcome of treatments as a contrast to a baseline treatment. The other uses a classical two-way linear predictor with main effects for treatment and trial. In this article, we compare both types of model and explore under which conditions they give equivalent results. We illustrate practical advantages of the two-way model using two published datasets. In particular, it is shown that between-trial heterogeneity as well as inconsistency between different types of trial is straightforward to account for.

Biomimetic reliability strategies for self-healing vascular networks in engineering materials.

Self-healing via a vascular network is an active research topic, with several recent publications reporting the application and optimization of these systems. This work represents the first consideration of the probable failure modes of a self-healing system as a driver for network design. The critical failure modes of a proposed self-healing system based on a vascular network were identified via a failure modes, effects and criticality analysis and compared to those of the human circulatory system. A range of engineering and biomimetic design concepts to address these critical failure modes is suggested with minimum system mass the overall design driver for high-performance systems. Plant vasculature has been mimicked to propose a segregated network to address the risk of fluid leakage. This approach could allow a network to be segregated into six separate paths with a system mass penalty of only approximately 25%. Fluid flow interconnections that mimic the anastomoses of animal vasculatures can be used within a segregated network to balance the risk of failure by leakage and blockage. These biomimetic approaches define a design space that considers the existing published literature in the context of system reliability.

A comparison of experimental designs for selection in breeding trials with nested treatment structure.

Plant breeders frequently evaluate large numbers of entries in field trials for selection. Generally, the tested entries are related by pedigree. The simplest case is a nested treatment structure, where entries fall into groups or families such that entries within groups are more closely related than between groups. We found that some plant breeders prefer to plant close relatives next to each other in the field. This contrasts with common experimental designs such as the alpha-design, where entries are fully randomized. A third design option is to randomize in such a way that entries of the same group are separated as much as possible. The present paper compares these design options by simulation. Another important consideration is the type of model used for analysis. Most of the common experimental designs were optimized assuming that the model used for analysis has fixed treatment effects. With many entries that are related by pedigree, analysis based on a model with random treatment effects becomes a competitive alternative. In simulations, we therefore study the properties of best linear unbiased predictions (BLUP) of genetic effects based on a nested treatment structure under these design options for a range of genetic parameters. It is concluded that BLUP provides efficient estimates of genetic effects and that resolvable incomplete block designs such as the alpha-design with restricted or unrestricted randomization can be recommended.

Construction of resolvable spatial row-column designs.

Resolvable row-column designs are widely used in field trials to control variation and improve the precision of treatment comparisons. Further gains can often be made by using a spatial model or a combination of spatial and incomplete blocking components. Martin, Eccleston, and Gleeson presented some general principles for the construction of robust spatial block designs which were addressed by spatial designs based on the linear variance (LV) model. In this article we define the two-dimensional form of the LV model and investigate extensions of the Martin et al. principles for the construction of resolvable spatial row-column designs. The computer construction of efficient spatial designs is discussed and some comparisons made with designs constructed assuming an autoregressive variance structure.

Identification of quantitative trait loci for wood and fibre properties in two full-sib properties of Eucalyptus globulus.

Regions of the genome influencing wood and fibre traits in Eucalyptus globulus Labill. have been identified in two full-sib pedigrees that share a common male parent. The first pedigree, cross A, contains 148 progeny, and the second pedigree, cross B, contains 135 progeny. Subsets of progeny of these two controlled crosses were planted at seven sites throughout Australia in 1990. Wood cores were taken at 0.9 m above ground in 1997, and wood and fibre traits were analysed for each individual. Three quantitative trait loci (QTL) affecting wood density, one QTL affecting pulp yield and one QTL affecting microfibril angle have been located in both pedigrees, using single-factor analysis of variance. Other QTLs affecting these traits, as well as fibre length and cellulose content, were located in cross A only.

Effect of metal ion and water coordination on the structure of a gas-phase amino acid.

The mode of metal ion and water binding to the amino acid valine is investigated using both theory and experiment. Computations indicate that without water, the structure of valine is nonzwitterionic. Both Li(+) and Na(+) are coordinated to the nitrogen and carbonyl oxygen (NO coordination), whereas K(+) coordinates to both oxygens (OO coordination) of nonzwitterionic valine. The addition of a single water molecule does not significantly affect the relative energies calculated for the cationized valine clusters. Experimentally, the rates of water evaporation from clusters of Val.M(+)(H(2)O)(1), M = Li, Na, and K, are measured using blackbody infrared radiative dissociation. The dissociation rate from the valine complex is compared to water evaporation rates from model complexes of known structure. These results indicate that the metal ion in the lithiated and the sodiated clusters is NO-coordinated to nonzwitterionic valine, while that in the potassiated cluster has OO coordination, in full agreement with theory. The zwitterionic vs nonzwitterionic character of valine in the potassiated cluster cannot be distinguished experimentally. Extensive modeling provides strong support for the validity of inferring structural information from the kinetic data.

Psychiatric status, somatisation, and health care utilization of frequent attenders at the emergency department: a comparison with routine attenders.

Seventy-seven frequent attenders at an emergency department (ED) in an inner-city hospital in the UK (defined as seven or more visits in the previous 12 months) were compared with 182 patients who were attending the same department on a routine basis. Patients completed the Schedules for Clinical Assessment in Neuropsychiatry (SCAN) and the Short Form (SF)-36. Information was obtained on 64% of the frequent attenders and 45% underwent a detailed psychiatric assessment. Of the frequent attenders, 45% had psychiatric disorder and 49% had some form of an alcohol-related disorder. Compared with routine attenders, frequent attenders reported lower health status, had more psychiatric disorder (odds ratio: OR=8.2, 95% confidence interval: CI=3.8--18.1), had more general hospital admissions (OR=19.9, 95% CI=8.3--47.8), more psychiatric admissions (OR=167.5, 95% CI=9.5--2959.0), and more GP visits (95% CI for difference=-10.2 to -5.7). There was no evidence that frequent attenders had more somatisation than routine attenders. Specific treatment and management strategies need to be developed for this group of patients, although a substantial proportion may be difficult to engage in the treatment process.

Gas-phase reactions of hydrated alkaline earth metal ions, M2+ (H2O)n (M = Mg, Ca, Sr, Ba and n = 4-7), with benzene.

Gas-phase reactions of hydrated divalent alkaline earth metal ions and benzene were investigated by electrospray ionization Fourier-transform mass spectrometry. Rate constants for solvent-exchange reactions were determined as a function of hydration extent for Mg2+, Ca2+, Sr2+, and Ba2+ clusters containing four to seven water molecules each. All of the strontium and barium clusters react quickly with benzene. Barium reacts slightly faster than the corresponding strontium cluster with the same number of water molecules attached. For calcium, clusters with four and five water molecules react quickly, whereas those with six and seven water molecules do not. Magnesium with four water molecules reacts quickly, but not when five through seven water molecules are attached. The slow reactivity observed for some of these clusters indicates that the cation-pi interaction between the metal ion and benzene is partially screened by the surrounding water molecules. The reactivity of magnesium with seven water molecules is intermediate that of the hexa- and pentahydrate and the tetrahydrate. This result is consistent with the seventh water molecule being in the outer shell and much more weakly bound. The unusual trend in reactivity observed for magnesium may be due to the presence of mixed shell structures observed previously. These results are the first to provide information about the relative importance of cation-pi interactions in divalent metal ions as a function of metal hydration extent. Such studies should also provide a model and some insight into the relative binding affinities of divalent metal ions to aromatic residues on peptides and proteins.

Supercharged protein and peptide ions formed by electrospray ionization.

The multiple charging of large molecules in electrospray ionization provides key advantages for obtaining accurate molecular weights by mass spectrometry and for obtaining structural information by tandem mass spectrometry and MS(n) experiments. Addition of glycerol or m-nitrobenzyl alcohol into the electrospray solutions dramatically increases both the maximum observed charge state and the abundances of the high charge states of protein and peptide ions. Adding glycerol to acidified aqueous solutions of cytochrome c shifts the most abundant charge state from 17+ to 21+, shifts the maximum charge state from 20+ to 23+, and shifts the average charge state from 16.6+ to 20.9+. Much less m-nitrobenzyl alcohol (<1%) is required to produce similar results. With just 0.7% m-nitrobenzyl alcohol, even the 24+ charge state of cytochrome c is readily observed. Similar results are obtained with myoglobin and (Lys)4. For the latter molecule, the 5+ charge state is observed in the electrospray mass spectrum obtained from solutions containing 6.7% m-nitrobenzyl alcohol. This charge state corresponds to protonation of all basic sites in this peptide. Although the mechanism for enhanced charging is unclear, it does not appear to be a consequence of conformational changes of the analyte molecules. This method of producing highly charged protein ions should be useful for improving the performance of mass measurements on mass spectrometers with performances that decrease with increasing m/z. This should also be particularly useful for tandem mass spectrometry experiments, such as electron capture dissociation, for which highly charged ions are desired.

The Ampere and Electrical Standards.

This paper describes some of the major contributions to metrology and physics made by the NIST Electricity Division, which has existed since 1901. It was one of the six original divisions of the National Bureau of Standards. The Electricity Division provides dc and low-frequency calibrations for industrial, scientific, and research organizations, and conducts research on topics related to electrical metrology and fundamental constants. The early work of the Electricity Division staff included the development of precision standards, such as Rosa and Thomas standard resistors and the ac-dc thermal converter. Research contributions helped define the early international system of measurement units and bring about the transition to absolute units based on fundamental principles and physical and dimensional measurements. NIST research has helped to develop and refine electrical standards using the quantum Hall effect and the Josephson effect, which are both based on quantum physics. Four projects covering a number of voltage and impedance measurements are described in detail. Several other areas of current research at NIST are described, including the use of the Internet for international compatibility in metrology, determination of the fine-structure and Planck constants, and construction of the electronic kilogram.

Hysteresis and Related Error Mechanisms in the NIST Watt Balance Experiment.

The NIST watt balance experiment is being completely rebuilt after its 1998 determination of the Planck constant. That measurement yielded a result with an approximately 1×10(-7) relative standard uncertainty. Because the goal of the new incarnation of the experiment is a ten-fold decrease in uncertainty, it has been necessary to reexamine many sources of systematic error. Hysteresis effects account for a substantial portion of the projected uncertainty budget. They arise from mechanical, magnetic, and thermal sources. The new experiment incorporates several improvements in the apparatus to address these issues, including stiffer components for transferring the mass standard on and off the balance, better servo control of the balance, better pivot materials, and the incorporation of erasing techniques into the mass transfer servo system. We have carried out a series of tests of hysteresis sources on a separate system, and apply their results to the watt apparatus. The studies presented here suggest that our improvements can be expected to reduce hysteresis signals by at least a factor of 10-perhaps as much as a factor of 50-over the 1998 experiment.

Effects of solvent on the maximum charge state and charge state distribution of protein ions produced by electrospray ionization.

The effects of solvent composition on both the maximum charge states and charge state distributions of analyte ions formed by electrospray ionization were investigated using a quadrupole mass spectrometer. The charge state distributions of cytochrome c and myoglobin, formed from 47%/50%/3% water/solvent/acetic acid solutions, shift to lower charge (higher m/z) when the 50% solvent fraction is changed from water to methanol, to acetonitrile, to isopropanol. This is also the order of increasing gas-phase basicities of these solvents, although other physical properties of these solvents may also play a role. The effect is relatively small for these solvents, possibly due to their limited concentration inside the electrospray interface. In contrast, the addition of even small amounts of diethylamine (<0.4%) results in dramatic shifts to lower charge, presumably due to preferential proton transfer from the higher charge state ions to diethylamine. These results clearly show that the maximum charge states and charge state distributions of ions formed by electrospray ionization are influenced by solvents that are more volatile than water. Addition of even small amounts of two solvents that are less volatile than water, ethylene glycol and 2-methoxyethanol, also results in preferential deprotonation of higher charge state ions of small peptides, but these solvents actually produce an enhancement in the higher charge state ions for both cytochrome c and myoglobin. For instruments that have capabilities that improve with lower m/z, this effect could be taken advantage of to improve the performance of an analysis.

Medical clearance of psychiatric patients.

Patients with primarily psychiatric complaints commonly present to emergency departments. The challenges to emergency medicine physicians are to distinguish functional from organic disease and to address any underlying medical comorbidities. The key elements of this medical "clearance" evaluation of psychiatric patients are discussed. The tremendous importance of the history and physical examination are highlighted, as is judicious use of laboratory and radiographic testing. Toxicology testing is specifically addressed.

Energetics from slow infrared multiphoton dissociation of biomolecules.

Photodissociation kinetics of the protonated pentapeptide leucine enkephalin measured using a cw CO(2) laser and a Fourier-transform mass spectrometer are reported. A short induction period, corresponding to the time required to raise the internal energy of the ion population to a (dissociating) steady state, is observed. After this induction period, the dissociation data are accurately fit by first-order kinetics. A plot of the log of the unimolecular dissociation rate constant, k(uni), as a function of the log of laser power is linear at low laser powers (<9 W, k(uni) <0.05 s(-1)), but tapers off at high laser power (9-33 W, k(uni) = 0.05-7 s(-1)). The entire measured dissociation curve can be accurately fit by an exponential function plus a constant. The experiment is simulated using a master equation formalism. In the model, the laser radiation is described as an energetically flat-topped distribution which is spatially uniform. This description is consistent with experimental results which indicate that ion motion within the cell averages out spatial inhomogeneities in the laser light. The model has several adjustable parameters. The effect of varying these parameters on the calculated kinetics and power dependence curves is discussed. A procedure for determining a limited range of threshold dissociation energy, E(o), which fits both the measured induction period and power dependence curves, is presented. Using this procedure, E(o) of leucine enkephalin is determined to be 1.12-1.46 eV. This result is consistent with, although less precise than, values measured previously using blackbody infrared radiative dissociation. Although the blackbody dissociation results were used as a starting point to search for fits of the master equation model to experiment, these results demonstrate that it is, in principle, possible to determine a limited range of E(o) from slow infrared multiphoton dissociation data alone.

Structure of cationized glycine, gly.m (m = be, mg, ca, sr, ba), in the gas phase: intrinsic effect of cation size on zwitterion stability.

Interactions between divalent metal ions and biomolecules are common both in solution and in the gas phase. Here, the intrinsic effect of divalent alkaline earth metal ions (Be, Mg, Ca, Sr, Ba) on the structure of glycine in the absence of solvent is examined. Results from both density functional and Moller-Plesset theories indicate that for all metal ions except beryllium, the salt-bridge form of the ion, in which glycine is a zwitterion, is between 5 and 12 kcal/mol more stable than the charge-solvated structure in which glycine is in its neutral form. For beryllium, the charge-solvated structure is 5-8 kcal/mol more stable than the salt-bridge structure. Thus, there is a dramatic change in the structure of glycine with increased metal cation size. Using a Hartree-Fock-based partitioning method, the interaction between the metal ion and glycine is separated into electrostatic, charge transfer and deformation components. The charge transfer interactions are more important for stabilizing the charge-solvated structure of glycine with beryllium relative to magnesium. In contrast, the difference in stability between the charge-solvated and salt-bridge structure for magnesium is mostly due to electrostatic interactions that favor formation of the salt-bridge structure. These results indicate that divalent metal ions dramatically influence the structure of this simplest amino acid in the gas phase.

Effects of Gas-Phase Basicity on the Proton Transfer between Organic Bases and Trifluoroacetic Acid in the Gas Phase: Energetics of Charge Solvation and Salt Bridges.

The unimolecular dissociation pathways and kinetics of a series of protonated trimer ions consisting of two organic bases and trifluoroacetic acid were investigated using blackbody infrared radiative dissociation. Five bases with gas-phase basicities (GB) ranging from 238.4 to 246.2 kcal/mol were used. Both the dissociation pathways and the threshold dissociation energies depend on the GB of the base. Trimers consisting of the two most basic molecules dissociate to form protonated base monomers with an E(0) ~ 1.4 eV. Trimers consisting of the two least basic molecules dissociate to form protonated base dimers with an E(0) ~ 1.1-1.2 eV. These results indicate that the structures of the trimers change as a function of the GB of the basic molecule. The predominant structure of the protonated trimers consisting of the two most basic molecules is consistent with a salt bridge in which both of the basic molecules are protonated, and the trifluoroacetic acid molecule is deprotonated, whereas the predominant structure of the protonated trimers consisting of the two least basic molecules are consistent with charge-solvated complexes in which the proton is shared. The structure of the trimer consisting of the base of intermediate basicity is less clear; it dissociates to form primarily protonated base dimer, but has an E(0) ~ 1.2 eV. These results are consistent with the structure of this trimer as a salt bridge, but the resulting dissociation A(-). BH(+) product does not appear to be stable as an ion pair in the dissociative transition state.

Structures of protonated arginine dimer and bradykinin investigated by density functional theory: further support for stable gas-phase salt bridges.

The gas-phase structures and energetics of both protonated arginine dimer and protonated bradykinin were investigated using a combination of molecular mechanics with conformational searching to identify candidate low-energy structures, and density functional theory for subsequent minimization and energy calculations. For protonated arginine dimer, a good correlation (R = 0.88) was obtained between the molecular mechanics and EDF1 6-31+G* energies, indicating that mechanics with MMFF is suitable for finding low-energy conformers. For this ion, the salt-bridge or ion-zwitterion form was found to be 5.7 and 7.2 kcal/mol more stable than the simple protonated or ion-molecule form at the EDF1 6-31++G** and B3LYP 6-311++G** levels. For bradykinin, the correlation between the molecular mechanics and DFT energies was poor (R = 0.28), indicating that many low-energy structures are likely passed over in the mechanics conformational searching. This result suggests that structures of this larger peptide ion obtained using mechanics calculations alone are not necessarily reliable. The lowest energy structure of the salt-bridge form of bradykinin is 10.6 kcal/mol lower in energy (EDF1) than the lowest energy simple protonated form at the 6-311G* level. Similarly, the average energy of all salt-bridge structures investigated is 13.6 kcal/mol lower than the average of all the protonated forms investigated. To the extent that a sufficient number of structures are investigated, these results provide some additional support for the salt-bridge form of bradykinin in the gas phase.

Dissociation energies of deoxyribose nucleotide dimer anions measured using blackbody infrared radiative dissociation.

The dissociation kinetics of deprotonated deoxyribose nucleotide dimers were measured using blackbody infrared radiative dissociation. Experiments were performed with noncovalently bound dimers of phosphate, adenosine (dAMP), cytosine (dCMP), guanosine (dGMP), thymidine (dTMP), and the mixed dimers dAMP.dTMP and dGMP.dCMP. The nucleotide dimers fragment through two parallel pathways, resulting in formation of the individual nucleotide or nucleotide + HPO3 ion. Master equation modeling of this kinetic data was used to determine threshold dissociation energies. The dissociation energy of (dGMP.dCMP-H)- is much higher than that for the other nucleotide dimers. This indicates that there is a strong interaction between the nucleobases in this dimer, consistent with the existence of Watson-Crick hydrogen bonding between the base pairs. Molecular mechanics simulations indicate that Watson-Crick hydrogen bonding occurs in the lowest energy structures of (dGMP.dCMP-H)-, but not in (dAMP.dTMP-H)-. The trend in gas phase dissociation energies is similar to the trend in binding energies measured in nonaqueous solutions within experimental error. Finally, the acidity ordering of the nucleotides is determined to be dTMP < dGMP < dCMP < dAMP, where dAMP has the highest acidity (largest delta Gacid).

Hydration of gas-phase ions formed by electrospray ionization.

The hydration of gas-phase ions produced by electrospray ionization was investigated. Evidence that the hydrated ions are formed by two mechanisms is presented. First, solvent condensation during the expansion inside the electrospray source clearly occurs. Second, some solvent evaporation from more extensively solvated ions or droplets is apparent. To the extent that these highly solvated ions have solution-phase structures, then the final isolated gas-phase structure of the ion will be determined by the solvent evaporation process. This process was investigated for hydrated gramicidin S in a Fourier-transform mass spectrometer. Unimolecular dissociation rate constants of isolated gramicidin S ions with between 2 and 14 associated water molecules were measured. These rate constants increased from 16 to 230 s-1 with increasing hydration, with smaller values corresponding to magic numbers.