Integrating sequence information in microarray data analysis by free energy modeling

Microarray technology is a high-throughput method for genotyping and gene expression profiling. Limited sensitivity and specificity are one of the essential problems for this technology. Most of existing methods of microarray data analysis have an apparent limitation for they merely deal with the numerical part of microarray data and have made little use of gene sequence information. Because it's the gene sequences that precisely define the physical objects being measured by a microarray, it is natural to make the gene sequences an essential part of the data analysis. This dissertation focused on the development of free energy models to integrate sequence information in microarray data analysis. The models were used to characterize the mechanism of hybridization on microarrays and enhance sensitivity and specificity of microarray measurements. ^ Cross-hybridization is a major obstacle factor for the sensitivity and specificity of microarray measurements. In this dissertation, we evaluated the scope of cross-hybridization problem on short-oligo microarrays. The results showed that cross hybridization on arrays is mostly caused by oligo fragments with a run of 10 to 16 nucleotides complementary to the probes. Furthermore, a free-energy based model was proposed to quantify the amount of cross-hybridization signal on each probe. This model treats cross-hybridization as an integral effect of the interactions between a probe and various off-target oligo fragments. Using public spike-in datasets, the model showed high accuracy in predicting the cross-hybridization signals on those probes whose intended targets are absent in the sample. ^ Several prospective models were proposed to improve Positional Dependent Nearest-Neighbor (PDNN) model for better quantification of gene expression and cross-hybridization. ^ The problem addressed in this dissertation is fundamental to the microarray technology. We expect that this study will help us to understand the detailed mechanism that determines sensitivity and specificity on the microarrays. Consequently, this research will have a wide impact on how microarrays are designed and how the data are interpreted. ^

Quantum molecular dynamics and ab initio studies of the crystal structure of hydrogen and deuterium

Hydrogen isotopes play a critical role both in inertial and magnetic confinemen Nuclear Fusion. Since the preferent fuel needed for this technology is a mixture of deuterium and tritium. The study of these isotopes particularly at very low temperatures carries a technological interest in other applications. The present line promotes a deep study on the structural configuration that hydrogen and deuterium adopt at cryogenic temperatures and at high pressures. Typical conditions occurring in present Inertial Fusion target designs. Our approach is aims to determine the crystal structure characteristics, phase transitions and other parameters strongly correlated to variations of temperature and pressure.

On the formation of nitrogen oxides in the combustion processes of hydrogen in air

A research programme is being carried out at the Institute Nacional de Tecnica Aeroespacial of Spain, on several aspects of the formation of nitrogen oxides in continuous flow combustion systems, considering hydrogen and hydrocarbons as fuels. The research programme is fundamentally oriented on the basic aspects of the problem, although it also includes the study of the influence on the formation process of several operational and design variables of the combusters, such as type of fuels, fuel/air ratio, degree of mixing in premixed type flames, existence of droplets as compared with homogeneous combustion.This problem of nitrogen oxides formation is receiving lately great attention, specially in connection with automobile reciprocating engines and aircraft gas turbines. This is due to the fact of the increasing frequency and intensity of photochemical hazes or smog, typical of urban areas submitted to strong solar radiation, which are originated by the action on organic compounds of the oxidants resulting from the photochemical decomposition of nitrogen dioxide N02. In the combustion process almost all nitrogen oxides are in form of NO. This nitric oxide reacts with the oxygen of the air and forms N02, this reaction only taking place in or near the exhaust of tne motors, since the N0-02 reaction becomes frozen for the concentration existing in the atmosphere.

Soft solution fluorine-free synthesis of anatase nanoparticles with tailored morphology

TiO2 nanoparticles with tailored morphology have been synthesized under exceptionally soft conditions. The strategy is based on the use of a non-aqueous alcoholic reaction medium in which water traces, coming either from the air (atmospheric water) or from an ethanol–water azeotropic mixture (ethanol 96%), are incorporated in order to accelerate hydrolysis of the Ti–precursor. Moreover, organic surfactants have been used as capping agents so as to tailor crystal growth in certain preferential directions. Combinations of oleic acid and oleylamine, which lead to the formation of another surfactant, dioleamide, are employed instead of fluorine-based compounds, thus increasing the sustainability of the process. As a result, TiO2 nanostructured hierarchical microspheres and individual nanoparticles with exposed high-energy facets can be obtained at atmospheric pressure and temperatures as low as 78 °C.

Critical radius for hot-jet ignition of hydrogen-air mixtures

This study addresses deflagration initiation of lean and stoichiometric hydrogen–air mixtures by the sudden discharge of a hot jet of their adiabatic combustion products. The objective is to compute the minimum jet radius required for ignition, a relevant quantity of interest for safety and technological applications. For sufficiently small discharge velocities, the numerical solution of the problem requires integration of the axisymmetric Navier–Stokes equations for chemically reacting ideal-gas mixtures, supplemented by standard descriptions of the molecular transport terms and a suitably reduced chemical-kinetic mechanism for the chemistry description. The computations provide the variation of the critical radius for hot-jet ignition with both the jet velocity and the equivalence ratio of the mixture, giving values that vary between a few tens microns to a few hundred microns in the range of conditions explored. For a given equivalence ratio, the critical radius is found to increase with increasing injection velocities, although the increase is only moderately large. On the other hand, for a given injection velocity, the smallest critical radius is found at stoichiometric conditions.

The dTTPase mechanism of T7 DNA helicase resembles the binding change mechanism of the F1-ATPase

Bacteriophage T7 DNA helicase is a ring-shaped hexamer that catalyzes duplex DNA unwinding using dTTP hydrolysis as an energy source. Of the six potential nucleotide binding sites on the hexamer, we have found that three are noncatalytic sites and three are catalytic sites. The noncatalytic sites bind nucleotides with a high affinity, but dTTPs bound to these sites do not dissociate or hydrolyze through many dTTPase turnovers at the catalytic sites. The catalytic sites show strong cooperativity which leads to sequential binding and hydrolysis of dTTP. The elucidated dTTPase mechanism of the catalytic sites of T7 helicase is remarkably similar to the binding change mechanism of the ATP synthase. Based on the similarity, a general mechanism for hexameric helicases is proposed. In this mechanism, an F1-ATPase-like rotational movement around the single-stranded DNA, which is bound through the central hole of the hexamer, is proposed to lead to unidirectional translocation along single-stranded DNA and duplex DNA unwinding.

A host–guest system to study structure–function relationships of membrane fusion peptides

We designed a host–guest fusion peptide system, which is completely soluble in water and has a high affinity for biological and lipid model membranes. The guest sequences are those of the fusion peptides of influenza hemagglutinin, which are solubilized by a highly charged unstructured C-terminal host sequence. These peptides partition to the surface of negatively charged liposomes or erythrocytes and elicit membrane fusion or hemolysis. They undergo a conformational change from random coil to an obliquely inserted (≈33° from the surface) α-helix on binding to model membranes. Partition coefficients for membrane insertion were measured for influenza fusion peptides of increasing lengths (n = 8, 13, 16, and 20). The hydrophobic contribution to the free energy of binding of the 20-residue fusion peptide at pH 5.0 is −7.6 kcal/mol (1 cal = 4.18 J). This energy is sufficient to stabilize a “stalk” intermediate if a typical number of fusion peptides assemble at the site of membrane fusion. The fusion activity of the fusion peptides increases with each increment in length, and this increase strictly correlates with the hydrophobic binding energy and the angle of insertion.

Recent improvements in prediction of protein structure by global optimization of a potential energy function

Recent improvements of a hierarchical ab initio or de novo approach for predicting both α and β structures of proteins are described. The united-residue energy function used in this procedure includes multibody interactions from a cumulant expansion of the free energy of polypeptide chains, with their relative weights determined by Z-score optimization. The critical initial stage of the hierarchical procedure involves a search of conformational space by the conformational space annealing (CSA) method, followed by optimization of an all-atom model. The procedure was assessed in a recent blind test of protein structure prediction (CASP4). The resulting lowest-energy structures of the target proteins (ranging in size from 70 to 244 residues) agreed with the experimental structures in many respects. The entire experimental structure of a cyclic α-helical protein of 70 residues was predicted to within 4.3 Å α-carbon (Cα) rms deviation (rmsd) whereas, for other α-helical proteins, fragments of roughly 60 residues were predicted to within 6.0 Å Cα rmsd. Whereas β structures can now be predicted with the new procedure, the success rate for α/β- and β-proteins is lower than that for α-proteins at present. For the β portions of α/β structures, the Cα rmsd's are less than 6.0 Å for contiguous fragments of 30–40 residues; for one target, three fragments (of length 10, 23, and 28 residues, respectively) formed a compact part of the tertiary structure with a Cα rmsd less than 6.0 Å. Overall, these results constitute an important step toward the ab initio prediction of protein structure solely from the amino acid sequence.

Free energy reconstruction from nonequilibrium single-molecule pulling experiments

Laser tweezers and atomic force microscopes are increasingly used to probe the interactions and mechanical properties of individual molecules. Unfortunately, using such time-dependent perturbations to force rare molecular events also drives the system away from equilibrium. Nevertheless, we show how equilibrium free energy profiles can be extracted rigorously from repeated nonequilibrium force measurements on the basis of an extension of Jarzynski's remarkable identity between free energies and the irreversible work.

Top-down free-energy minimization on protein potential energy landscapes

The hierarchical properties of potential energy landscapes have been used to gain insight into thermodynamic and kinetic properties of protein ensembles. It also may be possible to use them to direct computational searches for thermodynamically stable macroscopic states, i.e., computational protein folding. To this end, we have developed a top-down search procedure in which conformation space is recursively dissected according to the intrinsic hierarchical structure of a landscape's effective-energy barriers. This procedure generates an inverted tree similar to the disconnectivity graphs generated by local minima-clustering methods, but it fundamentally differs in the manner in which the portion of the tree that is to be computationally explored is selected. A key ingredient is a branch-selection algorithm that takes advantage of statistically predictive properties of the landscape to guide searches down the tree branches that are most likely to lead to the physically relevant macroscopic states. Using the computational folding of a β-hairpin-forming peptide as an example, we show that such predictive properties indeed exist and can be used for structure prediction by free-energy global minimization.

Interactions between tRNA identity nucleotides and their recognition sites in glutaminyl-tRNA synthetase determine the cognate amino acid affinity of the enzyme.

Sequence-specific interactions between aminoacyl-tRNA synthetases and their cognate tRNAs both ensure accurate RNA recognition and prevent the binding of noncognate substrates. Here we show for Escherichia coli glutaminyl-tRNA synthetase (GlnRS; EC 6.1.1.18) that the accuracy of tRNA recognition also determines the efficiency of cognate amino acid recognition. Steady-state kinetics revealed that interactions between tRNA identity nucleotides and their recognition sites in the enzyme modulate the amino acid affinity of GlnRS. Perturbation of any of the protein-RNA interactions through mutation of either component led to considerable changes in glutamine affinity with the most marked effects seen at the discriminator base, the 10:25 base pair, and the anticodon. Reexamination of the identity set of tRNA(Gln) in the light of these results indicates that its constituents can be differentiated based upon biochemical function and their contribution to the apparent Gibbs' free energy of tRNA binding. Interactions with the acceptor stem act as strong determinants of tRNA specificity, with the discriminator base positioning the 3' end. The 10:25 base pair and U35 are apparently the major binding sites to GlnRS, with G36 contributing both to binding and recognition. Furthermore, we show that E. coli tryptophanyl-tRNA synthetase also displays tRNA-dependent changes in tryptophan affinity when charging a noncognate tRNA. The ability of tRNA to optimize amino acid recognition reveals a novel mechanism for maintaining translational fidelity and also provides a strong basis for the coevolution of tRNAs and their cognate synthetases.

Crystal structure of the membrane-exposed domain from a respiratory quinol oxidase complex with an engineered dinuclear copper center.

Cytochrome oxidase is a membrane protein complex that catalyzes reduction of molecular oxygen to water and utilizes the free energy of this reaction to generate a transmembrane proton gradient during respiration. The electron entry site in subunit II is a mixed-valence dinuclear copper center in enzymes that oxidize cytochrome c. This center has been lost during the evolution of the quinoloxidizing branch of cytochrome oxidases but can be restored by engineering. Herein we describe the crystal structures of the periplasmic fragment from the wild-type subunit II (CyoA) of Escherichia coli quinol oxidase at 2.5-A resolution and of the mutant with the engineered dinuclear copper center (purple CyoA) at 2.3-A resolution. CyoA is folded as an 11-stranded mostly antiparallel beta-sandwich followed by three alpha-helices. The dinuclear copper center is located at the loops between strands beta 5-beta 6 and beta 9-beta 10. The two coppers are at a 2.5-A distance and symmetrically coordinated to the main ligands that are two bridging cysteines and two terminal histidines. The residues that are distinct in cytochrome c and quinol oxidases are around the dinuclear copper center. Structural comparison suggests a common ancestry for subunit II of cytochrome oxidase and blue copper-binding proteins.

TRAP, the trp RNA-binding attenuation protein of Bacillus subtilis, is a toroid-shaped molecule that binds transcripts containing GAG or UAG repeats separated by two nucleotides.

The trp RNA-binding attenuation protein of Bacillus subtilis, TRAP, regulates both transcription and translation by binding to specific transcript sequences. The optimal transcript sequences required for TRAP binding were determined by measuring complex formation between purified TRAP protein and synthetic RNAs. RNAs were tested that contained repeats of different trinucleotide sequences, with differing spacing between the repeats. A transcript containing GAG repeats separated by two-nucleotide spacers was bound most tightly. In addition, transmission electron microscopy was used to examine the structure of TRAP and the TRAP-transcript complex. TRAP was observed to be a toroid-shaped oligomer when free or when bound to either a natural or a synthetic RNA.

Direct measurement of hydrogen bonding in DNA nucleotide bases by atomic force microscopy.

We have used self-assembled purines and pyrimidines on planar gold surfaces and on gold-coated atomic force microscope (AFM) tips to directly probe intermolecular hydrogen bonds. Electron spectroscopy for chemical analysis (ESCA) and thermal programmed desorption (TPD) measurements of the molecular layers suggested monolayer coverage and a desorption energy of about 25 kcal/mol. Experiments were performed under water, with all four DNA bases immobilized on AFM tips and flat surfaces. Directional hydrogen-bonding interaction between the tip molecules and the surface molecules could be measured only when opposite base-pair coatings were used. The directional interactions were inhibited by excess nucleotide base in solution. Nondirectional van der Waals forces were present in all other cases. Forces as low as two interacting base pairs have been measured. With coated AFM tips, surface chemistry-sensitive recognition atomic force microscopy can be performed.

The Policy Framework for the Promotion of Hydrogen and Fuel Cells in Europe: a Critical Assessment. Bruges European Economic Policy (BEEP) Briefing 19/2008

This paper reviews the current EU policy framework in view of its impact on hydrogen and fuel cell development. It screens EU energy policies, EU regulatory policies and EU spending policies. Key questions addressed are as follows: To what extent is the current policy framework conducive to hydrogen and fuel cell development? What barriers and inconsistencies can be identified? How can policies potentially promote hydrogen and fuel cells in Europe, taking into account the complex evolution of such a disruptive technology? How should the EU policy framework be reformed in view of a strengthened and more coherent approach? The paper concludes that the current EU policy framework does not hinder hydrogen development. Yet it does not constitute a strong push factor either. EU energy policies have the strongest impact on hydrogen and fuel cell development even though their potential is still underexploited. Regulatory policies have a weak but positive impact on hydrogen. EU spending policies show some inconsistencies. However, the large scale market development of hydrogen and fuel cells will require a new policy approach which comprises technology specific support as well as a supportive policy framework with a special regional dimension.