Comparison of metaheuristic strategies for peakbin selection in proteomic mass spectrometry data

Mass spectrometry (MS) data provide a promising strategy for biomarker discovery. For this purpose, the detection of relevant peakbins in MS data is currently under intense research. Data from mass spectrometry are challenging to analyze because of their high dimensionality and the generally low number of samples available. To tackle this problem, the scientific community is becoming increasingly interested in applying feature subset selection techniques based on specialized machine learning algorithms. In this paper, we present a performance comparison of some metaheuristics: best first (BF), genetic algorithm (GA), scatter search (SS) and variable neighborhood search (VNS). Up to now, all the algorithms, except for GA, have been first applied to detect relevant peakbins in MS data. All these metaheuristic searches are embedded in two different filter and wrapper schemes coupled with Naive Bayes and SVM classifiers.

Two-dimensional plasma expansion in a magnetic nozzle: Separation due to electron inertia

A previous axisymmetric model of the supersonic expansion of a collisionless, hot plasma in a divergent magnetic nozzle is extended here in order to include electron-inertia effects. Up to dominant order on all components of the electron velocity, electron momentum equations still reduce to three conservation laws. Electron inertia leads to outward electron separation from the magnetic streamtubes. The progressive plasma filling of the adjacent vacuum region is consistent with electron-inertia being part of finite electron Larmor radius effects, which increase downstream and eventually demagnetize the plasma. Current ambipolarity is not fulfilled and ion separation can be either outwards or inwards of magnetic streamtubes, depending on their magnetization. Electron separation penalizes slightly the plume efficiency and is larger for plasma beams injected with large pressure gradients. An alternative nonzero electron-inertia model [E. Hooper, J. Propul. Power 9, 757 (1993)] based on cold plasmas and current ambipolarity, which predicts inwards electron separation, is discussed critically. A possible competition of the gyroviscous force with electron-inertia effects is commented briefly.

Effect of impacting mass on firmness sensing of fruits.

There has been an increased interest in using impact techniques for sensing firmness of fruits and vegetables. When an impacter is used to impact a fruit, the impacting mass is an important parameter which affects both the impact signal and fruit damage. Results of theoretical analysis and tests conducted on two varieties of pears indicate that lowering the impacting mass results in amplifying the measured signal, reducing sensing errors, and minimizing damage to the fruit.

The p-Adaptive Biem Approach for Two-Dimensional Elasticity Analysis

This paper presents the development and application of the p-adaptive BIEM version in elastostatics. The basic concepts underlying the p-adaptive technique are summarized and discussed. Some Pascal pseudocodes which show the way how such a technique can be implemented easily in microcomputers are also provided. Both the applicability and the accuracy of the method proposed here are illustrated through a numerical example.

Two-dimensional mesh optimization in the finite element method

The solution to the problem of finding the optimum mesh design in the finite element method with the restriction of a given number of degrees of freedom, is an interesting problem, particularly in the applications method. At present, the usual procedures introduce new degrees of freedom (remeshing) in a given mesh in order to obtain a more adequate one, from the point of view of the calculation results (errors uniformity). However, from the solution of the optimum mesh problem with a specific number of degrees of freedom some useful recommendations and criteria for the mesh construction may be drawn. For 1-D problems, namely for the simple truss and beam elements, analytical solutions have been found and they are given in this paper. For the more complex 2-D problems (plane stress and plane strain) numerical methods to obtain the optimum mesh, based on optimization procedures have to be used. The objective function, used in the minimization process, has been the total potential energy. Some examples are presented. Finally some conclusions and hints about the possible new developments of these techniques are also given.

Hybridizing concentrated solar power (CSP) with biogas and biomethane as an alternative to natural gas: Analysis of environmental performance using LCA

Concentrating Solar Power (CSP) plants typically incorporate one or various auxiliary boilers operating in parallel to the solar field to facilitate start up operations, provide system stability, avoid freezing of heat transfer fluid (HTF) and increase generation capacity. The environmental performance of these plants is highly influenced by the energy input and the type of auxiliary fuel, which in most cases is natural gas (NG). Replacing the NG with biogas or biomethane (BM) in commercial CSP installations is being considered as a means to produce electricity that is fully renewable and free from fossil inputs. Despite their renewable nature, the use of these biofuels also generates environmental impacts that need to be adequately identified and quantified. This paper investigates the environmental performance of a commercial wet-cooled parabolic trough 50 MWe CSP plant in Spain operating according to two strategies: solar-only, with minimum technically viable energy non-solar contribution; and hybrid operation, where 12 % of the electricity derives from auxiliary fuels (as permitted by Spanish legislation). The analysis was based on standard Life Cycle Assessment (LCA) methodology (ISO 14040-14040). The technical viability and the environmental profile of operating the CSP plant with different auxiliary fuels was evaluated, including: NG; biogas from an adjacent plant; and BM withdrawn from the gas network. The effect of using different substrates (biowaste, sewage sludge, grass and a mix of biowaste with animal manure) for the production of the biofuels was also investigated. The results showed that NG is responsible for most of the environmental damage associated with the operation of the plant in hybrid mode. Replacing NG with biogas resulted in a significant improvement of the environmental performance of the installation, primarily due to reduced impact in the following categories: natural land transformation, depletion of fossil resources, and climate change. However, despite the renewable nature of the biofuels, other environmental categories like human toxicity, eutrophication, acidification and marine ecotoxicity scored higher when using biogas and BM.

Two-dimensional isostatic meshes in the finite element method

In a Finite Element (FE) analysis of elastic solids several items are usually considered, namely, type and shape of the elements, number of nodes per element, node positions, FE mesh, total number of degrees of freedom (dot) among others. In this paper a method to improve a given FE mesh used for a particular analysis is described. For the improvement criterion different objective functions have been chosen (Total potential energy and Average quadratic error) and the number of nodes and dof's of the new mesh remain constant and equal to the initial FE mesh. In order to find the mesh producing the minimum of the selected objective function the steepest descent gradient technique has been applied as optimization algorithm. However this efficient technique has the drawback that demands a large computation power. Extensive application of this methodology to different 2-D elasticity problems leads to the conclusion that isometric isostatic meshes (ii-meshes) produce better results than the standard reasonably initial regular meshes used in practice. This conclusion seems to be independent on the objective function used for comparison. These ii-meshes are obtained by placing FE nodes along the isostatic lines, i.e. curves tangent at each point to the principal direction lines of the elastic problem to be solved and they should be regularly spaced in order to build regular elements. That means ii-meshes are usually obtained by iteration, i.e. with the initial FE mesh the elastic analysis is carried out. By using the obtained results of this analysis the net of isostatic lines can be drawn and in a first trial an ii-mesh can be built. This first ii-mesh can be improved, if it necessary, by analyzing again the problem and generate after the FE analysis the new and improved ii-mesh. Typically, after two first tentative ii-meshes it is sufficient to produce good FE results from the elastic analysis. Several example of this procedure are presented.

The structure of the two amino-terminal domains of human ICAM-1 suggests how it functions as a rhinovirus receptor and as an LFA-1 integrin ligand

The normal function of human intercellular adhesion molecule-1 (ICAM-1) is to provide adhesion between endothelial cells and leukocytes after injury or stress. ICAM-1 binds to leukocyte function-associated antigen (LFA-1) or macrophage-1 antigen (Mac-1). However, ICAM-1 is also used as a receptor by the major group of human rhinoviruses and is a catalyst for the subsequent viral uncoating during cell entry. The three-dimensional atomic structure of the two amino-terminal domains (D1 and D2) of ICAM-1 has been determined to 2.2-Å resolution and fitted into a cryoelectron microscopy reconstruction of a rhinovirus–ICAM-1 complex. Rhinovirus attachment is confined to the BC, CD, DE, and FG loops of the amino-terminal Ig-like domain (D1) at the end distal to the cellular membrane. The loops are considerably different in structure to those of human ICAM-2 or murine ICAM-1, which do not bind rhinoviruses. There are extensive charge interactions between ICAM-1 and human rhinoviruses, which are mostly conserved in both major and minor receptor groups of rhinoviruses. The interaction of ICAMs with LFA-1 is known to be mediated by a divalent cation bound to the insertion (I)-domain on the α chain of LFA-1 and the carboxyl group of a conserved glutamic acid residue on ICAMs. Domain D1 has been docked with the known structure of the I-domain. The resultant model is consistent with mutational data and provides a structural framework for the adhesion between these molecules.

Metabolism of an insect diuretic hormone by Malpighian tubules studied by liquid chromatography coupled with electrospray ionization mass spectrometry

The larger of two diuretic hormones of the tobacco hornworm, Manduca sexta, (Mas-DH) is a peptide of 41 residues. It is one of a family of seven currently known insect diuretic hormones that are similar to the corticotropin-releasing factor–urotensin–sauvagine family of peptides. We investigated the possible inactivation of Mas-DH by incubating it in vitro with larval Malpighian tubules (Mt), the target organ of the hormone. The medium was analyzed, and degradation products were identified, using on-line microbore reversed-phase liquid chromatography coupled to electrospray ionization mass spectrometry (RPLC-ESI-MS). This sensitive technique allows identification of metabolites of Mas-DH (present at an initial level of ≈1 μM). An accurate Mr value for a metabolite is usually sufficient for unambiguous identification. Mas-DH is cleaved by Mt proteases initially at L29–R30 and R30–A31 under our assay conditions; some Mas-DH is also oxidized, apparently at M2 and M11. The proteolysis can be inhibited by 5 mM EDTA, suggesting that divalent metals are needed for peptide cleavage. The oxidation of the hormone can be inhibited by catalase or 1 mM methionine, indicating that H2O2 or related reactive oxygen species are responsible for the oxidative degradation observed. RPLC-ESI-MS is shown here to be an elegant and efficient method for studying peptide hormone metabolism resulting from unknown proteases and pathways.

A comprehensive two-hybrid analysis to explore the yeast protein interactome

Protein–protein interactions play crucial roles in the execution of various biological functions. Accordingly, their comprehensive description would contribute considerably to the functional interpretation of fully sequenced genomes, which are flooded with novel genes of unpredictable functions. We previously developed a system to examine two-hybrid interactions in all possible combinations between the ≈6,000 proteins of the budding yeast Saccharomyces cerevisiae. Here we have completed the comprehensive analysis using this system to identify 4,549 two-hybrid interactions among 3,278 proteins. Unexpectedly, these data do not largely overlap with those obtained by the other project [Uetz, P., et al. (2000) Nature (London) 403, 623–627] and hence have substantially expanded our knowledge on the protein interaction space or interactome of the yeast. Cumulative connection of these binary interactions generates a single huge network linking the vast majority of the proteins. Bioinformatics-aided selection of biologically relevant interactions highlights various intriguing subnetworks. They include, for instance, the one that had successfully foreseen the involvement of a novel protein in spindle pole body function as well as the one that may uncover a hitherto unidentified multiprotein complex potentially participating in the process of vesicular transport. Our data would thus significantly expand and improve the protein interaction map for the exploration of genome functions that eventually leads to thorough understanding of the cell as a molecular system.

Observation via one-dimensional 13Calpha NMR of local conformational substates in thermal unfolding equilibria of a synthetic analog of the GCN4 leucine zipper.

Synthesis of a 33-residue, capped leucine zipper analogous to that in GCN4 is reported. Histidine and arginine residues are mutated to lysine to reduce the unfolding temperature. CD and ultracentrifugation studies indicate that the molecule is a two-stranded coiled coil under benign conditions. Versions of the same peptide are made with 99% 13Calpha at selected sites. One-dimensional 13C NMR spectra are assigned by inspection and used to study thermal unfolding equilibria over the entire transition from 8 to 73 degrees C. Spectra at the temperature extremes establish the approximate chemical shifts for folded and unfolded forms at each labeled site. Resonances for each amino acid appear at both locations at intermediate T, indicating that folded and unfolded forms interconvert slowly (> >2 ms) on the NMR time scale. Moreover, near room temperature, the structured form's resonance is double at several, but not all, sites, indicating at least two slowly interconverting, structured, local conformational substates. Analysis of the dynamics is possible. For example, near room temperature at the Val-9, Ala-24, and Gly-31 positions, the equilibrium constant for interconversion of the two structured forms is near unity and the time scale is > or= 10-20 ms.

Three-dimensional inverse modelling of magnetic anomaly sources based on a genetic algorithm.

We present a modelling method to estimate the 3-D geometry and location of homogeneously magnetized sources from magnetic anomaly data. As input information, the procedure needs the parameters defining the magnetization vector (intensity, inclination and declination) and the Earth's magnetic field direction. When these two vectors are expected to be different in direction, we propose to estimate the magnetization direction from the magnetic map. Then, using this information, we apply an inversion approach based on a genetic algorithm which finds the geometry of the sources by seeking the optimum solution from an initial population of models in successive iterations through an evolutionary process. The evolution consists of three genetic operators (selection, crossover and mutation), which act on each generation, and a smoothing operator, which looks for the best fit to the observed data and a solution consisting of plausible compact sources. The method allows the use of non-gridded, non-planar and inaccurate anomaly data and non-regular subsurface partitions. In addition, neither constraints for the depth to the top of the sources nor an initial model are necessary, although previous models can be incorporated into the process. We show the results of a test using two complex synthetic anomalies to demonstrate the efficiency of our inversion method. The application to real data is illustrated with aeromagnetic data of the volcanic island of Gran Canaria (Canary Islands).

Coherent electron dynamics in a two-dimensional random system with mobility edges

We study numerically the dynamics of a one-electron wavepacket in a two-dimensional random lattice with long-range correlated diagonal disorder in the presence of a uniform electric field. The time-dependent Schrodinger equation is used for this purpose. We find that the wavepacket displays Bloch-like oscillations associated with the appearance of a phase of delocalized states in the strong correlation regime. The amplitude of oscillations directly reflects the bandwidth of the phase and allows us to measure it. The oscillations reveal two main frequencies whose values are determined by the structure of the underlying potential in the vicinity of the wavepacket maximum.