MDASH: a multi-core-enabled program for structure solution from powder diffraction data

The simulated annealing approach to structure solution from powder diffraction data, as implemented in the DASH program, is easily amenable to parallelization at the individual run level. Modest increases in speed of execution can therefore be achieved by executing individual DASH runs on the individual cores of CPUs.

Boundary layer structure and decoupling from synoptic scale flow during NAMBLEX

This paper presents an overview of the meteorology and planetary boundary layer structure observed during the NAMBLEX field campaign to aid interpretation of the chemical and aerosol measurements. The campaign has been separated into five periods corresponding to the prevailing synoptic condition. Comparisons between meteorological measurements ( UHF wind profiler, Doppler sodar, sonic aneometers mounted on a tower at varying heights and a standard anemometer) and the ECMWF analysis at 10 m and 1100 m identified days when the internal boundary layer was decoupled from the synoptic flow aloft. Generally the agreement was remarkably good apart from during period one and on a few days during period four when the diurnal swing in wind direction implies a sea/land breeze circulation near the surface. During these periods the origin of air sampled at Mace Head would not be accurately represented by back trajectories following the winds resolved in ECMWF analyses. The wind profiler observations give a detailed record of boundary layer structure including an indication of its depth, average wind speed and direction. Turbulence statistics have been used to assess the height to which the developing internal boundary layer, caused by the increased surface drag at the coast, reaches the sampling location under a wide range of marine conditions. Sampling conducted below 10 m will be impacted by emission sources at the shoreline in all wind directions and tidal conditions, whereas sampling above 15 m is unlikely to be affected in any of the wind directions and tidal heights sampled during the experiment.

Can 4D-Var use dynamical information from targeted observations of a baroclinic structure?

Targeted observations are generally taken in regions of high baroclinicity, but often show little impact. One plausible explanation is that important dynamical information, such as upshear tilt, is not extracted from the targeted observations by the data assimilation scheme and used to correct initial condition error. This is investigated by generating pseudo targeted observations which contain a singular vector (SV) structure that is not present in the background field or routine observations, i.e. assuming that the background has an initial condition error with tilted growing structure. Experiments were performed for a single case-study with varying numbers of pseudo targeted observations. These were assimilated by the Met Office four-dimensional variational (4D-Var) data assimilation scheme, which uses a 6 h window for observations and background-error covariances calculated using the National Meteorological Centre (NMC) method. The forecasts were run using the operational Met Office Unified Model on a 24 km grid. The results presented clearly demonstrate that a 6 h window 4D-Var system is capable of extracting baroclinic information from a limited set of observations and using it to correct initial condition error. To capture the SV structure well (projection of 0.72 in total energy), 50 sondes over an area of 1×106 km2 were required. When the SV was represented by only eight sondes along an example targeting flight track covering a smaller area, the projection onto the SV structure was lower; the resulting forecast perturbations showed an SV structure with increased tilt and reduced initial energy. The total energy contained in the perturbations decreased as the SV structure was less well described by the set of observations (i.e. as fewer pseudo observations were assimilated). The assimilated perturbation had lower energy than the SV unless the pseudo observations were assimilated with the dropsonde observation errors halved from operational values. Copyright © 2010 Royal Meteorological Society

The effect of calcium removal from milk on casein micelle stability and structure

The total calcium level of raw skimmed milk was reduced by 10, 19, 29, 40 and 51% using Duolite® ion-exchange resin. The products were examined for concentrations of ionic calcium, sodium and potassium and the pH, ethanol stability, micelle diameter and ζ-potential were also measured. Ionic calcium decreased with removal of calcium and pH increased. Calcium removal resulted in an increase in the ethanol stability from 88% to above 100%. Casein micelle diameter increased as calcium was removed. The ζ-potential of the skimmed bulk milk was -24.4 mV, gradually becoming more negative with calcium removal to -30.6 mV after 51% calcium removal. The milk became more translucent as calcium was removed. To investigate the reversibility of this process, calcium chloride was added back to the depleted samples to restore their original total calcium content. At 51% removal, restoration of the total calcium level resulted in formation of clots. At levels of 10 and 19% calcium removal, the ethanol stability remained above 100%, but at higher levels of calcium removal the alcohol stability was adversely affected when the calcium was added back. Adding back calcium resulted in partial restoration of the original casein micelle diameter.

A benchmark method for global optimization problems in structure determination from powder diffraction data

Quasi-Newton-Raphson minimization and conjugate gradient minimization have been used to solve the crystal structures of famotidine form B and capsaicin from X-ray powder diffraction data and characterize the chi(2) agreement surfaces. One million quasi-Newton-Raphson minimizations found the famotidine global minimum with a frequency of ca 1 in 5000 and the capsaicin global minimum with a frequency of ca 1 in 10 000. These results, which are corroborated by conjugate gradient minimization, demonstrate the existence of numerous pathways from some of the highest points on these chi(2) agreement surfaces to the respective global minima, which are passable using only downhill moves. This important observation has significant ramifications for the development of improved structure determination algorithms.

Nuclear magnetic resonance structure of the N-terminal domain of nonstructural protein 3 from the severe acute respiratory syndrome coronavirus

This paper describes the structure determination of nsp3a, the N-terminal domain of the severe acute respiratory syndrome coronavirus (SARS-CoV) nonstructural protein 3. nsp3a exhibits a ubiquitin-like globular fold of residues 1 to 112 and a flexibly extended glutamic acid-rich domain of residues 113 to 183. In addition to the four beta-strands and two alpha-helices that are common to ubiquitin-like folds, the globular domain of nsp3a contains two short helices representing a feature that has not previously been observed in these proteins. Nuclear magnetic resonance chemical shift perturbations showed that these unique structural elements are involved in interactions with single-stranded RNA. Structural similarities with proteins involved in various cell-signaling pathways indicate possible roles of nsp3a in viral infection and persistence.

The structure and mechanism of serine acetyltransferase from Escherichia coli

Serine acetyltransferase (SAT) catalyzes the first step of cysteine synthesis in microorganisms and higher plants. Here we present the 2.2 Angstrom crystal structure of SAT from Escherichia coli, which is a dimer of trimers, in complex with cysteine. The SAT monomer consists of an amino-terminal alpha-helical domain and a carboxyl- terminal left-handed beta-helix. We identify His(158) and Asp(143) as essential residues that form a catalytic triad with the substrate for acetyl transfer. This structure shows the mechanism by which cysteine inhibits SAT activity and thus controls its own synthesis. Cysteine is found to bind at the serine substrate site and not the acetyl-CoA site that had been reported previously. On the basis of the geometry around the cysteine binding site, we are able to suggest a mechanism for the O-acetylation of serine by SAT. We also compare the structure of SAT with other left-handed beta-helical structures.

Structure and biogenesis of the capsular F1 antigen from Yersinia pestis: Preserved folding energy drives fiber formation

Most gram-negative pathogens express fibrous adhesive virulence organelles that mediate targeting to the sites of infection. The F1 capsular antigen from the plague pathogen Yersinia pestis consists of linear fibers of a single subunit (Caf1) and serves as a prototype for nonpilus organelles assembled via the chaperone/usher pathway. Genetic data together with high-resolution X-ray structures corresponding to snapshots of the assembly process reveal the structural basis of fiber formation. Comparison of chaperone bound Caf1 subunit with the subunit in the fiber reveals a novel type of conformational change involving the entire hydrophobic core of the protein. The observed conformational change suggests that the chaperone traps a high-energy folding intermediate of Caf1. A model is proposed in which release of the subunit allows folding to be completed, driving fiber formation.

Structure and rheology of aqueous micellar solutions and gels formed from an associative poly(oxybutylene)-poly(oxyethylene)-poly(oxybutylene) triblock copolymer

The structure and shear flow behaviour of aqueous micellar solutions and gels formed by an amphiphilic poly(oxybutylene)-poly(oxyethylene)-poly(oxybutylene) triblock copolymer with a lengthy hydrophilic poly(oxyethylene) block has been investigated by rheology, small angle neutron scattering (SANS) and small-angle X-ray scattering (SAXS). SANS revealed that bridging of chains between micelles introduces, in the micellar solution, an attractive long-range component which can be described through a potential of interaction corresponding to sticky soft spheres. The strength of the attractive interaction increases with increasing concentration. Rheology showed that the dependence of the storage modulus with temperature can be explained as a function of the micellar bridging, micellisation and phase morphology. SAXS studies showed that the orientation adopted by the system in the get phase under shear is similar to that previously observed by us for the gel phase of a poly(oxyethylene)-poly(oxybutylene) diblock copolymer with a long poly(oxyethylene) chain, suggesting that the micellar corona/core length ratio and not the architecture of the block copolymer influences the alignment of the gel phase under shear.

Crystal and molecular structure of the thioether-analogue of PEEK from X-ray powder data and diffraction-modelling

Analysis of X-ray powder data for the melt-crystallisable aromatic poly(thioether thioether ketone) [-S-Ar-S-Ar-CO-Ar](n), ('PTTK', Ar= 1,4-phenylene), reveals that it adopts a crystal structure very different from that established for its ether-analogue PEEK. Molecular modelling and diffraction-simulation studies of PTTK show that the structure of this polymer is analogous to that of melt-crystallised poly(thioetherketone) [-SAr-CO-Ar](n) in which the carbonyl linkages in symmetry-related chains are aligned anti-parallel to one another. and that these bridging units are crystallographically interchangeable. The final model for the crystal structure of PTTK is thus disordered, in the monoclinic space group 121a (two chains per unit cell), with cell dimensions a = 7.83, b = 6.06, c = 10.35 angstrom, beta = 93.47 degrees. (c) 2005 Elsevier Ltd. All rights reserved.

Structure of AuCN determined from total neutron diffraction

The structure of gold cyanide, AuCN, has been determined at 10 and 300 K using total neutron diffraction. The structure consists of infinite -Au-(CN)-Au-(CN)-Au-(CN)- linear chains, hexagonally packed, with the gold atoms in sheets. The Au-C and Au-N bond lengths are found to be identical, with d(Au-C/N) = 1.9703(5) Angstrom at 300 K. This work supersedes a previous study, by others, which used Rietveld analysis of neutron Bragg diffraction in isolation, and found these bonds to have significantly different lengths (Deltad = 0.24 Angstrom) at 300 K. The total correlation function, T(r), at 10 and 300 K, has been modeled using information derived from total diffraction. The broadening of inter- and intrachain correlations differs markedly due to random displacements of the chains in the direction of the chain axes. This is a consequence of the relatively weak bonding between the chains. An explanation for the negative thermal expansion in the c-direction, which occurs between 10 and 300 K, is presented.

A differential thermal expansion approach to crystal structure determination from powder diffraction data

Differential thermal expansion over the range 90-210 K has been applied successfully to determine the crystal structure of chlorothiazide from synchrotron powder diffraction data using direct methods. Key to the success of the approach is the use of a multi-data-set Pawley refinement to extract a set of reflection intensities that is more 'single-crystal-like' than those extracted from a single data set. The improvement in reflection intensity estimates is quantified by comparison with reference single-crystal intensities. (C) 2008 International Union of Crystallography Printed in Singapore - all rights reserved