Binding energies of protonated betaine complexes: a probe of zwitterion structure in the gas phase

The dissociation kinetics of proton-bound dimers of betaine with molecules of comparable gas-phase basicity were investigated using blackbody infrared radiative dissociation (BIRD). Threshold dissociation energies were obtained from these data using master equation modeling. For bases that have comparable or higher gas-phase basicity, the binding energy of the protonated base.betaine complex is approximately 1.4 eV. For molecules that are approximately 2 kcal/mol or more less basic, the dissociation energy of the complexes is approximately 1.2 eV. The higher binding energy of the former is attributed to an ion-zwitterion structure which has a much larger ion-dipole interaction. The lower binding energy for molecules that are approximately 2 kcal/mol or more less basic indicates that an ion-molecule structure is more favored. Semiempirical calculations at both the AM1 and PM3 levels indicate the most stable ion-molecule structure is one in which the base interacts with the charged quaternary ammonium end of betaine. These results indicate that the measurement of binding energies of neutral molecules to biological ions could provide a useful probe for the presence of zwitterions and salt bridges in the gas phase. From the BIRD data, the gas-phase basicity of betaine obtained from the kinetic method is found to be 239.2 +/- 1.0 kcal/mol. This value is in excellent agreement with the value of 239.3 kcal/mol (298 K) from ab initio calculations at the MP2/6-31+g** level. The measured value is slightly higher than those reported previously. This difference is attributed to entropy effects. The lower ion internal energy and longer time frame of BIRD experiments should provide values closer to those at standard temperature.

Unimolecular reactions of dihydrated alkaline earth metal dications M2+(H2O)2, M = Be, Mg, Ca, Sr, and Ba: salt-bridge mechanism in the proton-transfer reaction M2+(H2O)2 --> MOH+ + H3O

The unimolecular reactivity of M(2+)(H(2)O)(2), M = Be, Mg, Ca, Sr, and Ba, is investigated by density functional theory. Dissociation of the complex occurs either by proton transfer to form singly charged metal hydroxide, MOH(+), and protonated water, H(3)O(+), or by loss of water to form M(2+)(H(2)O) and H(2)O. Charge transfer from water to the metal forming H(2)O(+) and M(+)(H(2)O) is not favorable for any of the metal complexes. The relative energetics of these processes are dominated by the metal dication size. Formation of MOH(+) proceeds first by one water ligand moving to the second solvation shell followed by proton transfer to this second-shell water molecule and subsequent Coulomb explosion. These hydroxide formation reactions are exothermic with activation energies that are comparable to the water binding energy for the larger metals. This results in a competition between proton transfer and loss of a water molecule. The arrangement with one water ligand in the second solvation shell is a local minimum on the potential energy surface for all metals except Be. The two transition states separating this intermediate from the reactant and the products are identified. The second transition state determines the height of the activation barrier and corresponds to a M(2+)-OH(-)-H(3)O(+) "salt-bridge" structure. The computed B3LYP energy of this structure can be quantitatively reproduced by a simple ionic model in which Lewis charges are localized on individual atoms. This salt-bridge arrangement lowers the activation energy of the proton-transfer reaction by providing a loophole on the potential energy surface for the escape of H(3)O(+). Similar salt-bridge mechanisms may be involved in a number of proton-transfer reactions in small solvated metal ion complexes, as well as in other ionic reactions.

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.

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.

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.

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.

Hydration energies and structures of alkaline earth metal ions, M2+(H2O)n, n = 5-7, M = Mg, Ca, Sr, and Ba

The evaporation of water from hydrated alkaline earth metal ions, produced by electrospray ionization, was studied in a Fourier transform mass spectrometer. Zero-pressure-limit dissociation rate constants for loss of a single water molecule from the hydrated divalent metal ions, M(2+)(H(2)O)(n) (M = Mg, Ca, and Sr for n = 5-7, and M = Ba for n = 4-7), are measured as a function of temperature using blackbody infrared radiative dissociation. From these values, zero-pressure-limit Arrhenius parameters are obtained. By modeling the dissociation kinetics using a master equation formalism, threshold dissociation energies (E(o)) are determined. These reactions should have a negligible reverse activation barrier; therefore, E(o) values should be approximately equal to the binding energy or hydration enthalpy at 0 K. For the hepta- and hexahydrated ions at low temperature, binding energies follow the trend expected on the basis of ionic radii: Mg > Ca > Sr > Ba. For the hexahydrated ions at high temperature, binding energies follow the order Ca > Mg > Sr > Ba. The same order is observed for the pentahydrated ions. Collisional dissociation experiments on the tetrahydrated species result in relative dissociation rates that directly correlate with the size of the metals. These results indicate the presence of two isomers for hexahydrated magnesium ions: a low-temperature isomer in which the six water molecules are located in the first solvation shell, and a high-temperature isomer with the most likely structure corresponding to four water molecules in the inner shell and two water molecules in the second shell. These results also indicate that the pentahydrated magnesium ions have a structure with four water molecules in the first solvation shell and one in the outer shell. The dissociation kinetics for the hexa- and pentahydrated clusters of Ca(2+), Sr(2+), and Ba(2+) are consistent with structures in which all the water molecules are located in the first solvation shell.

Enhanced positive cloud-to-ground lightning in thunderstorms ingesting smoke from fires

Smoke from forest fires in southern Mexico was advected into the U.S. southern plains from April to June 1998. Cloud-to-ground lightning (CG) flash data from the National Lightning Detection Network matched against satellite-mapped aerosol plumes imply that thunderstorms forming in smoke-contaminated air masses generated large amounts of lightning with positive polarity (+CGs). During 2 months, nearly half a million flashes in the southern plains exhibited +CG percentages that were triple the climatological norm. The peak currents in these +CGs were double the expected value. These thunderstorms also produced abnormally high numbers of mesospheric optical sprites.

Activation energies for dissociation of double strand oligonucleotide anions: evidence for watson-crick base pairing in vacuo

The dissociation kinetics of a series of complementary and noncomplementary DNA duplexes, (TGCA)(2) (3-), (CCGG)(2) (3-), (AATTAAT)(2) (3-), (CCGGCCG)(2) (3-), A(7)*T(7) (3-), A(7)*A(7) (3-), T(7)*T(7) (3-), and A(7)*C(7) (3-) were investigated using blackbody infrared radiative dissociation in a Fourier transform mass spectrometer. From the temperature dependence of the unimolecular dissociation rate constants, Arrhenius activation parameters in the zero-pressure limit are obtained. Activation energies range from 1.2 to 1.7 eV, and preexponential factors range from 10(13) to 10(19) s(-1). Dissociation of the duplexes results in cleavage of the noncovalent bonds and/or cleavage of covalent bonds leading to loss of a neutral nucleobase followed by backbone cleavage producing sequence-specific (a - base) and w ions. Four pieces of evidence are presented which indicate that Watson-Crick (WC) base pairing is preserved in complementary DNA duplexes in the gas phase: i. the activation energy for dissociation of the complementary dimer, A(7)*T(7) (3-), to the single strands is significantly higher than that for the related noncomplementary A(7)*A(7) (3-) and T(7)*T(7) (3-) dimers, indicating a stronger interaction between strands with a specific base sequence, ii. extensive loss of neutral adenine occurs for A(7)*A(7) (3-) and A(7)*C(7) (3-) but not for A(7)*T(7) (3-) consistent with this process being shut down by WC hydrogen bonding, iii. a correlation is observed between the measured activation energy for dissociation to single strands and the dimerization enthalpy (-DeltaH(d)) in solution, and iv. molecular dynamics carried out at 300 and 400 K indicate that WC base pairing is preserved for A(7)*T(7) (3-) duplex, although the helical structure is essentially lost. In combination, these results provide strong evidence that WC base pairing can exist in the complete absence of solvent.

Separation and characterization of amines from individual atrial gland vesicles of Aplysia californica

Several amine-containing components of individual vesicles from the atrial gland of Aplysia californica were identified with capillary electrophoresis (CE). On-line derivatization with naphthalene-2,3-dicarboxaldehyde was performed, and the derivatized amine-containing components were detected with laser-induced fluorescence (LIF). Amino acids, including taurine, that had not been determined previously in atrial gland vesicles were observed by using CE-LIF, and their identities were confirmed with CE, HPLC, NMR, and electrospray ionization mass spectrometry. The finding that taurine is packaged and stored into secretory vesicles supports the hypothesis that taurine may exhibit neuromodulatory activity. The bioactive peptides, well-known to be in atrial gland vesicles, were detected in lysed vesicle samples fractionated with HPLC and analyzed by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. These peptides were also observed in single-vesicle runs with CE-LIF. The atrial gland vesicles (ranging from 0.5 to 2 microns diameter and 65 aL to 4 fL volume, respectively) studied in this work represent the smallest biological entities to be analyzed chemically on an individual basis.

Analysis of isomeric mixtures using blackbody infrared radiative dissociation: determining isomeric purity and obtaining individual tandem mass spectra simultaneously

A new method that makes possible, for the first time, simultaneous acquisition of individual dissociation mass spectra of isomeric ions in mixtures is presented. This method exploits the exquisite sensitivity of blackbody infrared radiative dissociation kinetics to minor differences in ion structure. Instead of separating precursor ions based on mass (isomers have identical mass), fragment ions are related to their original precursor ions on the basis of rate constants for dissociation. Mixtures of the peptide isomers des-R1 and des-R9 bradykinin are dissociated simultaneously at several temperatures. By fitting the kinetic data to double-exponential functions, the dissociation rate constant and abundance of each isomer in the mixture are obtained. To overcome the difficulty of fitting double-exponential functions, a novel global analysis method is used in which several dissociation data sets collected at different temperatures are simultaneously fit. The kinetic data measured at multiple temperatures are modeled with the preexponentials (corresponding to the abundance of each isomer) as "global" parameters which are constant for all data sets and the exponentials (rate constants) as "local" variables which differ for each data set. The use of global parameters significantly improves the accuracy with which abundances and dissociation rate constants of each individual compound can be obtained from the mixture data. Fragment ions produced from a mixture of these two isomers are related back to their respective precursor ions from the kinetic data. Thus, not only can the composition of the isomeric mixture be determined but an individual tandem mass spectrum of each component in the mixture can be obtained.

Probing single secretory vesicles with capillary electrophoresis

Secretory vesicles obtained from the atrial gland of the gastropod mollusk Aplysia californica were chemically analyzed individually with a combination of optical trapping, capillary electrophoresis separation, and a laser-induced fluorescence detection. With the use of optical trapping, a single vesicle that had attoliters (10(-18) liters) of volume was introduced into the tapered inlet of a separation capillary. Once the vesicle was injected, it was lysed, and its components were fluorescently labeled with naphthalene-2, 3-dicarboxaldehyde before separation. The resultant electropherograms indicated distinct variations in the contents of single vesicles.

Activation of Peptide ions by blackbody radiation: factors that lead to dissociation kinetics in the rapid energy exchange limit

Unimolecular rate constants for blackbody infrared radiative dissociation (BIRD) were calculated for the model protonated peptide (AlaGly)(n) (n = 2-32) using a variety of dissociation parameters. Combinations of dissociation threshold energies ranging from 0.8 to 1.7 eV and transition entropies corresponding to Arrhenius preexponential factors ranging from very "tight" (A(infinity) = 10(9.9) s(-1)) to "loose" (A(infinity) = 10(16.8) s(-1)) were selected to represent dissociation parameters within the experimental temperature range (300-520 K) and kinetic window (k(uni) = 0.001-0.20 s(-1)) typically used in the BIRD experiment. Arrhenius parameters were determined from the temperature dependence of these values and compared to those in the rapid energy exchange (REX) limit. In this limit, the internal energy of a population of ions is given by a Boltzmann distribution, and kinetics are the same as those in the traditional high-pressure limit. For a dissociation process to be in this limit, the rate of photon exchange between an ion and the vacuum chamber walls must be significantly greater than the dissociation rate. Kinetics rapidly approach the REX limit either as the molecular size or threshold dissociation energy increases or as the transition-state entropy or experimental temperature decreases. Under typical experimental conditions, peptide ions larger than 1.6 kDa should be in the REX limit. Smaller ions may also be in the REX limit depending on the value of the threshold dissociation energy and transition-state entropy. Either modeling or information about the dissociation mechanism must be known in order to confirm REX limit kinetics for these smaller ions. Three principal factors that lead to the size dependence of REX limit kinetics are identified. With increasing molecular size, rates of radiative absorption and emission increase, internal energy distributions become relatively narrower, and the microcanonical dissociation rate constants increase more slowly over the energy distribution of ions. Guidelines established here should make BIRD an even more reliable method to obtain information about dissociation energetics and mechanisms for intermediate size molecules.

Dissociation energetics and mechanisms of leucine enkephalin (M + H)+ and (2M + X)+ ions (X = H, Li, Na, K, and Rb) measured by blackbody infrared radiative dissociation

The dissociation kinetics of protonated leucine enkephalin and its proton and alkali metal bound dimers were investigated by blackbody infrared radiative dissociation in a Fourier-transform mass spectrometer. From the temperature dependence of the unimolecular dissociation rate constants, Arrhenius activation parameters in the zero-pressure limit are obtained. Protonated leucine enkephalin dissociates to form b(4) and (M-H(2)O)(+) ions with an average activation energy (E(a)) of 1.1 eV and an A factor of 10(10.5) s(-1). The value of the A factor indicates that these dissociation processes are rearrangements. The b(4) ions subsequently dissociate to form a(4) ions via a process with a relatively high activation energy (1.3 eV), but one that is entropically favored. For the cationized dimers, the thermal stability decreases with increasing cation size, consistent with a simple electrostatic interaction in these noncovalent ion-molecule complexes. The E(a) and A factors are indistinguishable within experimental error with values of approximately 1.5 eV and 10(17) s(-1), respectively. Although not conclusive, results from master equation modeling indicate that all these BIRD processes, except for b(4) --> a(4), are in the rapid energy exchange limit. In this limit, the internal energy of the precursor ion population is given by a Boltzmann distribution and information about the energetics and dynamics of the reaction are obtained directly from the measured Arrhenius parameters.