Use of a novel Hill-climbing genetic algorithm in protein folding simulations

We have developed a novel Hill-climbing genetic algorithm (GA) for simulation of protein folding. The program (written in C) builds a set of Cartesian points to represent an unfolded polypeptide's backbone. The dihedral angles determining the chain's configuration are stored in an array of chromosome structures that is copied and then mutated. The fitness of the mutated chain's configuration is determined by its radius of gyration. A four-helix bundle was used to optimise simulation conditions, and the program was compared with other, larger, genetic algorithms on a variety of structures. The program ran 50% faster than other GA programs. Overall, tests on 100 non-redundant structures gave comparable results to other genetic algorithms, with the Hill-climbing program running from between 20 and 50% faster. Examples including crambin, cytochrome c, cytochrome B and hemerythrin gave good secondary structure fits with overall alpha carbon atom rms deviations of between 5 and 5.6 Angstrom with an optimised hydrophobic term in the fitness function. (C) 2003 Elsevier Ltd. All rights reserved.

Advances in continuously profiling the thermodynamic state of the boundary layer: integration of measurements and methods

This paper describes advances in ground-based thermodynamic profiling of the lower troposphere through sensor synergy. The well-documented integrated profiling technique (IPT), which uses a microwave profiler, a cloud radar, and a ceilometer to simultaneously retrieve vertical profiles of temperature, humidity, and liquid water content (LWC) of nonprecipitating clouds, is further developed toward an enhanced performance in the boundary layer and lower troposphere. For a more accurate temperature profile, this is accomplished by including an elevation scanning measurement modus of the microwave profiler. Height-dependent RMS accuracies of temperature (humidity) ranging from 0.3 to 0.9 K (0.5–0.8 g m−3) in the boundary layer are derived from retrieval simulations and confirmed experimentally with measurements at distinct heights taken during the 2005 International Lindenberg Campaign for Assessment of Humidity and Cloud Profiling Systems and its Impact on High-Resolution Modeling (LAUNCH) of the German Weather Service. Temperature inversions, especially of the lower boundary layer, are captured in a very satisfactory way by using the elevation scanning mode. To improve the quality of liquid water content measurements in clouds the authors incorporate a sophisticated target classification scheme developed within the European cloud observing network CloudNet. It allows the detailed discrimination between different types of backscatterers detected by cloud radar and ceilometer. Finally, to allow IPT application also to drizzling cases, an LWC profiling method is integrated. This technique classifies the detected hydrometeors into three different size classes using certain thresholds determined by radar reflectivity and/or ceilometer extinction profiles. By inclusion into IPT, the retrieved profiles are made consistent with the measurements of the microwave profiler and an LWC a priori profile. Results of IPT application to 13 days of the LAUNCH campaign are analyzed, and the importance of integrated profiling for model evaluation is underlined.

Europe’s biocidal products directive: benefits and costs in urban pest management

The European Commission’s Biocidal Products Directive (Council Directive 98/8 EC), known as the BPD, is the largest regulatory exercise ever to affect the urban pest control industry. Although focussed in the European Union its impact is global because any company selling pest control products in the EU must follow its principles. All active substances, belonging to 23 different biocidal product types, come within the Directive’s scope of regulatory control. This will eventually involve re-registration of all existing products, as well as affecting any new product that comes to the market. Some active substances, such as the rodenticides and insecticides, are already highly regulated in Europe but others, such as embalming fluids, masonry preservatives, disinfectants and repellents/attractants will come under intensive regulatory scrutiny for the first time. One of the purposes of the Directive is to offer enhanced protection for human health and the environment. The potential benefit for suppliers of pest control products is mutual recognition of regulatory product dossiers across 25 Member States of the European Union. This process, requiring harmonisation of all regulatory decision-making processes, should reduce duplicated effort and, potentially, allow manufacturers speedier access to European markets. However, the cost to industry is enormous, both in terms of the regulatory resources required to assemble BPD dossiers and the development budgets required to conduct studies to meet its new standards. The cost to regulatory authorities is also tremendous, in terms of the need to upgrade staff capabilities to meet new challenges and the volume of the work expected by the Commission when they are appointed the Rapporteur Member State (RMS) for an active substance. Users of pest control products will pay a price too. The increased regulatory costs of maintaining products in the European market are likely to be passed on, at least in part, to users. Furthermore, where the costs of meeting new regulatory requirements cannot be recouped from product sales, many well-known products may leave the market. For example, it seems that in future few rodenticides that are not anticoagulants will be available within the EU. An understanding of the BPD is essential to those who intend to place urban pest control products on the European market and may be useful to those considering the harmonisation of regulatory processes elsewhere. This paper reviews the operation of the first stages of the BPD for rodenticides, examines the potential benefits and costs of the legislation to the urban pest control industry and looks forward to the next stages of implementation involving all insecticides used in urban pest management.

Quantitative precipitation estimation over ocean using Bayesian approach from microwave observations during the typhoon season

We have developed a new Bayesian approach to retrieve oceanic rain rate from the Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI), with an emphasis on typhoon cases in the West Pacific. Retrieved rain rates are validated with measurements of rain gauges located on Japanese islands. To demonstrate improvement, retrievals are also compared with those from the TRMM/Precipitation Radar (PR), the Goddard Profiling Algorithm (GPROF), and a multi-channel linear regression statistical method (MLRS). We have found that qualitatively, all methods retrieved similar horizontal distributions in terms of locations of eyes and rain bands of typhoons. Quantitatively, our new Bayesian retrievals have the best linearity and the smallest root mean square (RMS) error against rain gauge data for 16 typhoon overpasses in 2004. The correlation coefficient and RMS of our retrievals are 0.95 and ~2 mm hr-1, respectively. In particular, at heavy rain rates, our Bayesian retrievals outperform those retrieved from GPROF and MLRS. Overall, the new Bayesian approach accurately retrieves surface rain rate for typhoon cases. Accurate rain rate estimates from this method can be assimilated in models to improve forecast and prevent potential damages in Taiwan during typhoon seasons.

Boundary-layer friction in midlatitude cyclone

Results from an idealized three-dimensional baroclinic life-cycle model are interpreted in a potential vorticity (PV) framework to identify the physical mechanisms by which frictional processes acting in the atmospheric boundary layer modify and reduce the baroclinic development of a midlatitude storm. Considering a life cycle where the only non-conservative process acting is boundary-layer friction, the rate of change of depth-averaged PV within the boundary layer is governed by frictional generation of PV and the flux of PV into the free troposphere. Frictional generation of PV has two contributions: Ekman generation, which is directly analogous to the well-known Ekman-pumping mechanism for barotropic vortices, and baroclinic generation, which depends on the turning of the wind in the boundary layer and low-level horizontal temperature gradients. It is usually assumed, at least implicitly, that an Ekman process of negative PV generation is the mechanism whereby friction reduces the strength and growth rates of baroclinic systems. Although there is evidence for this mechanism, it is shown that baroclinic generation of PV dominates, producing positive PV anomalies downstream of the low centre, close to developing warm and cold fronts. These PV anomalies are advected by the large-scale warm conveyor belt flow upwards and polewards, fluxed into the troposphere near the warm front, and then advected westwards relative to the system. The result is a thin band of positive PV in the lower troposphere above the surface low centre. This PV is shown to be associated with a positive static stability anomaly, which Rossby edge wave theory suggests reduces the strength of the coupling between the upper- and lower-level PV anomalies, thereby reducing the rate of baroclinic development. This mechanism, which is a result of the baroclinic dynamics in the frontal regions, is in marked contrast with simple barotropic spin-down ideas. Finally we note the implications of these frictionally generated PV anomalies for cyclone forecasting.

Assimilation of non-synoptic observations

A system for continuous data assimilation described recently (Bengtsson & Gustavsson, 1971) has been further developed and tested under more realistic conditions. A balanced barotropic model is used and the integration is performed over an octagon covering the area to the north of 20° N. Comparisons have been made between using data from the actual aerological network and data from a satellite in a polar orbit. The result of the analyses has been studied in different subregions situated in data sparse as well as in data dense areas. The errors of the analysis have also been studied in the wave spectrum domain. Updating is performed using data generated by the model but also by model-independent data. Rather great differences are obtained between the two experiments especially with respect to the ultra-long waves. The more realistic approach gives much larger analysis error. In general the satellite updating yields somewhat better result than the updating from the conventional aerological network especially in the data sparse areas over the oceans. Most of the experiments are performed by a satellite making 200 observations/track, a sidescan capability of 40° and with a RMS-error of 20 m. It is found that the effect of increasing the number of satellite observations from 100 to 200 per orbit is almost negligible. Similarly the effect is small of improving the observations by diminishing the RMS-error below a certain value. An observing system using two satellites 90° out of phase has also been investigated. This is found to imply a substantial improvement. Finally an experiment has been performed using actual SIRS-soundings from NIMBUS IV. With respect to the very small number of soundings at 500 mb, 142 during 48 hours, the result can be regarded as quite satisfactory.

An experiment in the Assimilation of Data in Dynamical Analysis

A system for continuous data assimilation is presented and discussed. To simulate the dynamical development a channel version of a balanced barotropic model is used and geopotential (height) data are assimilated into the models computations as data become available. In the first experiment the updating is performed every 24th, 12th and 6th hours with a given network. The stations are distributed at random in 4 groups in order to simulate 4 areas with different density of stations. Optimum interpolation is performed for the difference between the forecast and the valid observations. The RMS-error of the analyses is reduced in time, and the error being smaller the more frequent the updating is performed. The updating every 6th hour yields an error in the analysis less than the RMS-error of the observation. In a second experiment the updating is performed by data from a moving satellite with a side-scan capability of about 15°. If the satellite data are analysed at every time step before they are introduced into the system the error of the analysis is reduced to a value below the RMS-error of the observation already after 24 hours and yields as a whole a better result than updating from a fixed network. If the satellite data are introduced without any modification the error of the analysis is reduced much slower and it takes about 4 days to reach a comparable result to the one where the data have been analysed.

Estimation of sea surface temperature from the Spinning Enhanced Visible and Infra Red Imager, improved using numerical weather prediction

Most of the operational Sea Surface Temperature (SST) products derived from satellite infrared radiometry use multi-spectral algorithms. They show, in general, reasonable performances with root mean square (RMS) residuals around 0.5 K when validated against buoy measurements, but have limitations, particularly a component of the retrieval error that relates to such algorithms' limited ability to cope with the full variability of atmospheric absorption and emission. We propose to use forecast atmospheric profiles and a radiative transfer model to simulate the algorithmic errors of multi-spectral algorithms. In the practical case of SST derived from the Spinning Enhanced Visible and Infrared Imager (SEVIRI) onboard Meteosat Second Generation (MSG), we demonstrate that simulated algorithmic errors do explain a significant component of the actual errors observed for the non linear (NL) split window algorithm in operational use at the Centre de Météorologie Spatiale (CMS). The simulated errors, used as correction terms, reduce significantly the regional biases of the NL algorithm as well as the standard deviation of the differences with drifting buoy measurements. The availability of atmospheric profiles associated with observed satellite-buoy differences allows us to analyze the origins of the main algorithmic errors observed in the SEVIRI field of view: a negative bias in the inter-tropical zone, and a mid-latitude positive bias. We demonstrate how these errors are explained by the sensitivity of observed brightness temperatures to the vertical distribution of water vapour, propagated through the SST retrieval algorithm.

Infrared Observations of XTE J1118+480

Infrared observations of the outbursting black hole XTE J1118+480 (ATEL #383) were performed using SQIID on the Kitt Peak National Observatory 2.1m telescope. Observations spanning 2005 January 15.42-15.58 found it somewhat fainter than the previous outburst (IAUC # 7394 , # 7407 ), at average brightness J=12.91+/-0.03, H=12.50+/-0.03, K=11.95+/-0.03. The colors again correspond to an approximately flat spectrum in F_nu. No orbital variation is apparent, but there is substantial unresolved rapid variability with rms amplitude 22% in K (between 2s exposures). Further observations are planned nightly until Jan 21.

Improvements in Aerosol Optical Depth Estimation Using Multiangle CHRIS/PROBA Images

A method has been developed to estimate aerosol optical depth (AOD) over land surfaces using high spatial resolution, hyperspectral, and multiangle Compact High Resolution Imaging Spectrometer (CHRIS)/Project for On Board Autonomy (PROBA) images. The CHRIS instrument is mounted aboard the PROBA satellite and provides up to 62 bands. The PROBA satellite allows pointing to obtain imagery from five different view angles within a short time interval. The method uses inversion of a coupled surface/atmosphere radiative transfer model and includes a general physical model of angular surface reflectance. An iterative process is used to determine the optimum value providing the best fit of the corrected reflectance values for a number of view angles and wavelengths with those provided by the physical model. This method has previously been demonstrated on data from the Advanced Along-Track Scanning Radiometer and is extended here to the spectral and angular sampling of CHRIS/PROBA. The values obtained from these observations are validated using ground-based sun-photometer measurements. Results from 22 image sets show an rms error of 0.11 in AOD at 550 nm, which is reduced to 0.06 after an automatic screening procedure.

Photometric redshifts and k-corrections for the Sloan Digital Sky Survey Data Release 7

We present a catalogue of galaxy photometric redshifts and k-corrections for the Sloan Digital Sky Survey Data Release 7 (SDSS-DR7), available on the World Wide Web. The photometric redshifts were estimated with an artificial neural network using five ugriz bands, concentration indices and Petrosian radii in the g and r bands. We have explored our redshift estimates with different training sets, thus concluding that the best choice for improving redshift accuracy comprises the main galaxy sample (MGS), the luminous red galaxies and the galaxies of active galactic nuclei covering the redshift range 0 < z < 0.3. For the MGS, the photometric redshift estimates agree with the spectroscopic values within rms = 0.0227. The distribution of photometric redshifts derived in the range 0 < z(phot) < 0.6 agrees well with the model predictions. k-corrections were derived by calibration of the k-correct_v4.2 code results for the MGS with the reference-frame (z = 0.1) (g - r) colours. We adopt a linear dependence of k-corrections on redshift and (g - r) colours that provide suitable distributions of luminosity and colours for galaxies up to redshift z(phot) = 0.6 comparable to the results in the literature. Thus, our k-correction estimate procedure is a powerful, low computational time algorithm capable of reproducing suitable results that can be used for testing galaxy properties at intermediate redshifts using the large SDSS data base.

AFM analysis of bleaching effects on dental enamel microtopography

The purpose of this in vitro study was to test a new methodology to evaluate the effects of 35% hydrogen peroxide agent on the microtopography of sound enamel using an atomic force microscope (AFM). The buccal sound surfaces of three extracted human lower incisors were used, without polishing the surfaces to maintain them with natural morphology. These unpolished surfaces were subjected to bleaching procedure with 35% hydrogen peroxide that consisted of 4 applications of the bleaching agent on enamel surfaces for 10 min each application. Surface images were obtained in a 15 mu m x 15 mu m area using an AFM. The roughness (Ra and RMS) and the power spectral density (PSD) were obtained before and after the bleaching treatment. As results we could inquire that the PSD analyses were very suitable to identifying the morphological changes on the surfaces, while the Ra and RMS parameters were insufficient to represent the morphological alterations promoted by bleaching procedure on enamel. The morphological wavelength in the range of visible light spectrum(380-750 nm) was analyzed, showing a considerable increase of the PSD with the bleaching treatment. (C) 2009 Elsevier B. V. All rights reserved.

Er : YAB nanoparticles and vitreous thin films by the polymeric precursor method

The synthesis of Y(0.9)Er(0.1)Al(3)(BO(3))(4) crystalline powders and vitreous thin films were studied. Precursor solutions were obtained using a modified polymeric precursor method using D-sorbitol as complexant agent. The chemical reactions were described. Y(0.)9Er(0.1)Al(3)(BO(3))(4) composition presents good thermal stability with regard to crystallization. The Y(0.9)Er(0.1)Al(3)(BO(3))(4) crystallized phase can be obtained at 1,150 degrees C, in agreement with other authors. Crack- and porosity-free films were obtained with very small grain size and low RMS roughness. The films thickness revealed to be linearly dependent on precursor solution viscosity, being the value of 25 mPa s useful to prepare high-quality amorphous multi-layers (up to similar to 800 nm) at 740 degrees C during 2 h onto silica substrates by spin coating with a gyrset technology.

Tetrachlorocarbonyliridates: Water-Soluble Carbon Monoxide Releasing Molecules Rate-Modulated by the Sixth Ligand

A new family of compounds is presented as potential carbon monoxide releasing molecules (CORMs). These compounds, based on tetrachlorocarbonyliridate(III) derivatives, were synthesized and fully characterized by X-ray diffraction, electrospray mass spectrometry, IR. NMR, and density functional theory calculations. The rate of CO release was studied via the myoglobin assay. The results showed that the rate depends on the nature of the sixth ligand, trans to CO, and that a significant modulation on the release rate can be produced by changing the ligand. The reported compounds are soluble in aqueous media, and the rates of CO release are comparable with those for known CORMs, releasing CO at a rate of 0.03-0.58 mu M min(-1) in a 10 mu M solution of myoglobin and 10 mu M of the complexes.

Characterization of anodic silicon oxide films grown in room temperature ionic liquids

The electroformation of silicon oxide was performed in two room temperature ionic liquids (RTIL), 1-butyl-3-methyl-imidazolium bis(trifluoromethane sulfonyl) imide (BMITFSI) and N-n-butyl-N-methylpiperidinium bis(trifluoromethane sulfonyl) imide (BMPTFSI). This phenomenon was studied by electrochemical techniques and it was observed that the oxide growth follows a high-field mechanism. X-ray Photoelectron Spectroscopy experiments have shown that a non-stoichiometric oxide film was formed, related to the low water content present in both RTILs (< 30 ppm). The roughness values obtained by using AFM technique of the silicon surface after etching with HF was 1.5 nm (RMS). The electrochemical impedance spectroscopy at low frequencies range was interpreted as a resistance in parallel with a CPE element, the capacitance obtained was associated with the dielectric nature of the oxide formed and the resistance was interpreted considering the chemical dissolution of the oxide by the presence of the TFSI anion. The CPE element was associated with the surface roughness and the very thin oxide film obtained. (C) 2007 Elsevier Ltd. All rights reserved.