Funcions densitat i semblança molecular quàntica: nous desenvolupaments i aplicacions

La present tesi, tot i que emmarcada dins de la teoria de les Mesures Semblança Molecular Quántica (MQSM), es deriva en tres àmbits clarament definits: - La creació de Contorns Moleculars de IsoDensitat Electrònica (MIDCOs, de l'anglès Molecular IsoDensity COntours) a partir de densitats electròniques ajustades. - El desenvolupament d'un mètode de sobreposició molecular, alternatiu a la regla de la màxima semblança. - Relacions Quantitatives Estructura-Activitat (QSAR, de l'anglès Quantitative Structure-Activity Relationships). L'objectiu en el camp dels MIDCOs és l'aplicació de funcions densitat ajustades, ideades inicialment per a abaratir els càlculs de MQSM, per a l'obtenció de MIDCOs. Així, es realitza un estudi gràfic comparatiu entre diferents funcions densitat ajustades a diferents bases amb densitats obtingudes de càlculs duts a terme a nivells ab initio. D'aquesta manera, l'analogia visual entre les funcions ajustades i les ab initio obtinguda en el ventall de representacions de densitat obtingudes, i juntament amb els valors de les mesures de semblança obtinguts prèviament, totalment comparables, fonamenta l'ús d'aquestes funcions ajustades. Més enllà del propòsit inicial, es van realitzar dos estudis complementaris a la simple representació de densitats, i són l'anàlisi de curvatura i l'extensió a macromolècules. La primera observació correspon a comprovar no només la semblança dels MIDCOs, sinó la coherència del seu comportament a nivell de curvatura, podent-se així observar punts d'inflexió en la representació de densitats i veure gràficament aquelles zones on la densitat és còncava o convexa. Aquest primer estudi revela que tant les densitats ajustades com les calculades a nivell ab initio es comporten de manera totalment anàloga. En la segona part d'aquest treball es va poder estendre el mètode a molècules més grans, de fins uns 2500 àtoms. Finalment, s'aplica part de la filosofia del MEDLA. Sabent que la densitat electrònica decau ràpidament al allunyar-se dels nuclis, el càlcul d'aquesta pot ser obviat a distàncies grans d'aquests. D'aquesta manera es va proposar particionar l'espai, i calcular tan sols les funcions ajustades de cada àtom tan sols en una regió petita, envoltant l'àtom en qüestió. Duent a terme aquest procés, es disminueix el temps de càlcul i el procés esdevé lineal amb nombre d'àtoms presents en la molècula tractada. En el tema dedicat a la sobreposició molecular es tracta la creació d'un algorisme, així com la seva implementació en forma de programa, batejat Topo-Geometrical Superposition Algorithm (TGSA), d'un mètode que proporcionés aquells alineaments que coincideixen amb la intuïció química. El resultat és un programa informàtic, codificat en Fortran 90, el qual alinea les molècules per parelles considerant tan sols nombres i distàncies atòmiques. La total absència de paràmetres teòrics permet desenvolupar un mètode de sobreposició molecular general, que proporcioni una sobreposició intuïtiva, i també de forma rellevant, de manera ràpida i amb poca intervenció de l'usuari. L'ús màxim del TGSA s'ha dedicat a calcular semblances per al seu ús posterior en QSAR, les quals majoritàriament no corresponen al valor que s'obtindria d'emprar la regla de la màxima semblança, sobretot si hi ha àtoms pesats en joc. Finalment, en l'últim tema, dedicat a la Semblança Quàntica en el marc del QSAR, es tracten tres aspectes diferents: - Ús de matrius de semblança. Aquí intervé l'anomenada matriu de semblança, calculada a partir de les semblances per parelles d'entre un conjunt de molècules. Aquesta matriu és emprada posteriorment, degudament tractada, com a font de descriptors moleculars per a estudis QSAR. Dins d'aquest àmbit s'han fet diversos estudis de correlació d'interès farmacològic, toxicològic, així com de diverses propietats físiques. - Aplicació de l'energia d'interacció electró-electró, assimilat com a una forma d'autosemblança. Aquesta modesta contribució consisteix breument en prendre el valor d'aquesta magnitud, i per analogia amb la notació de l'autosemblança molecular quàntica, assimilar-la com a cas particular de d'aquesta mesura. Aquesta energia d'interacció s'obté fàcilment a partir de programari mecanoquàntic, i esdevé ideal per a fer un primer estudi preliminar de correlació, on s'utilitza aquesta magnitud com a únic descriptor. - Càlcul d'autosemblances, on la densitat ha estat modificada per a augmentar el paper d'un substituent. Treballs previs amb densitats de fragments, tot i donar molt bons resultats, manquen de cert rigor conceptual en aïllar un fragment, suposadament responsable de l'activitat molecular, de la totalitat de l'estructura molecular, tot i que les densitats associades a aquest fragment ja difereixen degut a pertànyer a esquelets amb diferents substitucions. Un procediment per a omplir aquest buit que deixa la simple separació del fragment, considerant així la totalitat de la molècula (calcular-ne l'autosemblança), però evitant al mateix temps valors d'autosemblança no desitjats provocats per àtoms pesats, és l'ús de densitats de Forats de fermi, els quals es troben definits al voltant del fragment d'interès. Aquest procediment modifica la densitat de manera que es troba majoritàriament concentrada a la regió d'interès, però alhora permet obtenir una funció densitat, la qual es comporta matemàticament igual que la densitat electrònica regular, podent-se així incorporar dins del marc de la semblança molecular. Les autosemblances calculades amb aquesta metodologia han portat a bones correlacions amb àcids aromàtics substituïts, podent així donar una explicació al seu comportament. Des d'un altre punt de vista, també s'han fet contribucions conceptuals. S'ha implementat una nova mesura de semblança, la d'energia cinètica, la qual consisteix en prendre la recentment desenvolupada funció densitat d'energia cinètica, la qual al comportar-se matemàticament igual a les densitats electròniques regulars, s'ha incorporat en el marc de la semblança. A partir d'aquesta mesura s'han obtingut models QSAR satisfactoris per diferents conjunts moleculars. Dins de l'aspecte del tractament de les matrius de semblança s'ha implementat l'anomenada transformació estocàstica com a alternativa a l'ús de l'índex Carbó. Aquesta transformació de la matriu de semblança permet obtenir una nova matriu no simètrica, la qual pot ser posteriorment tractada per a construir models QSAR.

Anthropogenic snowfall events in the UK: examples of urban weather modification?

Snow in the UK is generally associated with synoptic or mesoscale weather systems, thus snowfall during quiescent anticyclonic conditions is surprising and might not even be forecast. Consequently it could present a hazard. Snowfall during anticyclonic freezing fog conditions at Didcot and Hereford in December 2006 is investigated here. These two snowfalls seem to present circumstances in which anthropogenically-produced aerosols could have provided ice nuclei within the freezing fog, and therefore might provide characteristic examples of Anthropogenic Snowfall Events (ASEs).

Evaluation of a large-eddy model simulation of a mixed-phase altocumulus cloud using microwave radiometer, lidar and Doppler radar data

Using the Met Office large-eddy model (LEM) we simulate a mixed-phase altocumulus cloud that was observed from Chilbolton in southern England by a 94 GHz Doppler radar, a 905 nm lidar, a dual-wavelength microwave radiometer and also by four radiosondes. It is important to test and evaluate such simulations with observations, since there are significant differences between results from different cloud-resolving models for ice clouds. Simulating the Doppler radar and lidar data within the LEM allows us to compare observed and modelled quantities directly, and allows us to explore the relationships between observed and unobserved variables. For general-circulation models, which currently tend to give poor representations of mixed-phase clouds, the case shows the importance of using: (i) separate prognostic ice and liquid water, (ii) a vertical resolution that captures the thin layers of liquid water, and (iii) an accurate representation the subgrid vertical velocities that allow liquid water to form. It is shown that large-scale ascents and descents are significant for this case, and so the horizontally averaged LEM profiles are relaxed towards observed profiles to account for these. The LEM simulation then gives a reasonable. cloud, with an ice-water path approximately two thirds of that observed, with liquid water at the cloud top, as observed. However, the liquid-water cells that form in the updraughts at cloud top in the LEM have liquid-water paths (LWPs) up to half those observed, and there are too few cells, giving a mean LWP five to ten times smaller than observed. In reality, ice nucleation and fallout may deplete ice-nuclei concentrations at the cloud top, allowing more liquid water to form there, but this process is not represented in the model. Decreasing the heterogeneous nucleation rate in the LEM increased the LWP, which supports this hypothesis. The LEM captures the increase in the standard deviation in Doppler velocities (and so vertical winds) with height, but values are 1.5 to 4 times smaller than observed (although values are larger in an unforced model run, this only increases the modelled LWP by a factor of approximately two). The LEM data show that, for values larger than approximately 12 cm s(-1), the standard deviation in Doppler velocities provides an almost unbiased estimate of the standard deviation in vertical winds, but provides an overestimate for smaller values. Time-smoothing the observed Doppler velocities and modelled mass-squared-weighted fallspeeds shows that observed fallspeeds are approximately two-thirds of the modelled values. Decreasing the modelled fallspeeds to those observed increases the modelled IWC, giving an IWP 1.6 times that observed.

Linear nucleation of polymers

The use of high-melting fibres as linear nuclei for quiescent polymeric melts is instrumental in providing the superior mechanical properties of polymeric self-composites. It also has inherent advantages in the elucidation of fundamental aspects of polymeric crystallization and self-organization, not least in allowing systematic microscopic studies of polymeric crystallization from nucleation through to the growth interface. This has demonstrated explicitly that lamellae develop in two distinct ways, for slower and faster growth, depending on whether fold packing has or has not time to order before the next molecular layer is added with only the former leading to banded growth in linear polyethylene. Other gains in understanding concern cellulation and morphological instability, internuclear interference, isothermal lamellar thickening and banded growth being a consequence of the partial relief of initial surface stress. (c) 2006 Elsevier Ltd. All rights reserved.

Wave functions for the methane molecule

If the potential field due to the nuclei in the methane molecule is expanded in terms of a set of spherical harmonics about the carbon nucleus, only the terms involving s, f, and higher harmonic functions differ from zero in the equilibrium configuration. Wave functions have been calculated for the equilibrium configuration, first including only the spherically symmetric s term in the potential, and secondly including both the s and the f terms. In the first calculation the complete Hartree-Fock S.C.F. wave functions were determined; in the second calculation a variation method was used to determine the best form of the wave function involving f harmonics. The resulting wave functions and electron density functions are presented and discussed

Vibrational intensities in methane

The absorption intensities of the two infra-red active vibrations in methane have been obtained from a perturbation calculation on the equilibrium wave functions derived in the preceding paper. The perturbation field is the change in the potential field due to the nuclei which results from moving the nuclei in the vibrational coordinate concerned, and a simplified form of second order perturbation theory, developed by Pople and Schofield, is used for the calculation. The main approximation involved is the neglect of f and higher harmonics in the spherical harmonic expansion of the nuclear field. The resulting dipole moment derivatives are approximately three times larger than the experimental values, but they show qualitative features and sign relationships which are significant.

Pollen ultrastructure in anther cultures of Datura innoxia. III. Incomplete microspore division

During the microspore division in Datura innoxia, the mitotic spindle is oriented in planes both perpendicular (PE) and oblique (OB) to the spore wall against which the nucleus is situated. However, irrespective of polarity, the usual type of hemispherical wall is laid down at cytokinesis and isolates the generative cell from the rest of the pollen grain (type A). In PE spores the vegetative nucleus initially occupies a central position in the pollen grain, whereas in OB spores the vegetative nucleus lies at the periphery of the grain close to the generative cell. In anther cultures initiated just before the microspore division is due to take place, no marked change can be observed in either orientation or symmetry of the mitotic spindle when the spores divide. In some, however, cytokinesis is disrupted and deposition of the hemispherical wall arrested. In the absence of a complete wall, differentiation of the generative cell cannot take place and binucleate pollen grains are formed having 2 vegetative-type nuclei (type B). The 2 nuclei in the B pollens are always situated against the pollen-grain wall, suggesting that the disruption phenomenon is related to the OB spores. The incomplete wall always makes contact with the intine on the intine-side of the spindle. Wall material may be represented merely as short stubs projecting out from the intine into the cytoplasm, in which event the 2 nuclei lie close to each other and are separated by only a narrow zone of cytoplasm. In other grains the wall is partially developed between the nuclei and terminates at varying distances from the tonoplast; in these, the nuclei are separated by a wider zone of cytoplasm. The significance of these binucleate grains in pollen embryogenesis is discussed.

Cell and nuclear degradation in sweetpotato [Ipomoea batatas (L.) Lam.] root meristems following exposure to salinity

Sodium chloride-induced cell and nuclear degradation in the root meristems of sweetpotato [Ipomoea batatas (L.) Lam.] were determined using fluorescent microscopy and flow cytometry analysis. Two sweetpotato cultivars were grown in liquid Murashige and Skoog medium and subjected to 0 mM and 500 mM NaCl, with or without 15 mM CaCl2, for periods up to 24 h. Changes to the nuclei of root meristematic cells showed a similar pattern of damage to the nuclei using both fluorescent microscopy and flow cytometry analysis. Damage occurring after only a few hours was followed by nuclear degradation at 24 h. Flow cytometry histograms showed a reduction in G1 and G2 nuclei and an increase in degraded nuclei in NaCl-stressed roots. Salinity-induced nuclear degradation was alleviated by the addition of CaCl2.

Microgemma vivaresi n. sp (Microsporidia, Tetramicridae), infecting liver and skeletal muscle of sea scorpions, Taurulus bubalis (Euphrasen 1786) (Osteichthyes, Cottidae), an inshore, littoral fish

The ultrastructure of a new microsporidian species Microgemmia vivaresi n. sp. causing liver cell xenoma formation in sea scorpions, Taurulus bubalis, is described. Stages of merogony, sporogony, and sporogenesis are mixed in the central cytoplasm of developing xenomas. All stages have unpaired nuclei. Uninucleate and multinucleate meronts lie within vacuoles formed from host endoplasmic reticulum and divide by binary or multiple fission. Sporonts, no longer in vacuoles, deposit plaques of surface coat on the plasma membrane that cause the surface to pucker. Division occurs at the Puckered stage into sporoblast mother cells, on which plaques join up to complete the surface coat. A final binary fission gives rise to sporoblasts. A dense globule, thought to be involved in polar tube synthesis, is gradually dispersed during spore maturation. Spores are broadly ovoid, have a large posterior vacuole, and measure 3.6 mu m x 2.1 pint (fresh). The polar tube has a short wide anterior section that constricts abruptly, then runs posteriad to coil about eight times around the posterior vacuole with granular contents. The polaroplast has up to 40 membranes arranged in pairs mostly attached to the wide region of the polar tube and directed posteriorty around a cytoplasm of a coarsely granular appearance. The species is placed alongside the type species Microgemmia hepaticus Ralphs and Matthews 1986 within the family Tetramicridae, which is transferred from the class Dihaplophasea to the class Haplophasea, as there is no evidence for the occurrence of a diplokaryotic phase.

Evidence for a reversal with age in the pattern of near-wall blood flow around aortic branches

The changes that occur with age in the distribution of atherosclerotic lesions around arterial branch points challenge accepted theories relating disease to haemodynamic stresses. We investigated whether flow near branch points changes with age in a way that can account for the different lesion distributions. Flow around 20 branches from immature and mature aortas was investigated by examining the length:width ratio and orientation of endothelial nuclei; these properties depend on the magnitude and direction of near-wall flows, respectively. There were significant changes in the pattern of nuclear shape with age, consistent with a reversal in the pattern of shear around branches. In control regions away from branches, there were no such changes. The role of haemodynamic stresses in atherogenesis may require re-evaluation in the light of these results. (C) 2003 Elsevier Ireland Ltd. All rights reserved.

Time course of gag protein assembly in HIV-1-infected cells: a study by immunoelectron microscopy

Recent biochemical studies have identified high molecular complexes of the HIV Gag precursor in the cytosol of infected cells. Using immunoelectron microscopy we studied the time course of the synthesis and assembly of a HIV Gag precursor protein (pr55gag) in Sf9 cells infected with recombinant baculovirus expressing the HIV gag gene. We also immunolabeled for pr55gag human T4 cells acutely or chronically infected with HIV-1. In Sf9 cells, the time course study showed that the first Gag protein appeared in the cytoplasm at 28-30 h p.i. and that budding started 6-8 h later. Colloidal gold particles, used to visualize the Gag protein, were first scattered randomly throughout the cytoplasm, but soon clusters representing 100 to 1000 copies of pr55gag were also observed. By contrast, in cells with budding or released virus-like particles the cytoplasm was virtually free of gold particles while the released virus-like particles were heavily labeled. Statistical analysis showed that between 80 and 90% of the gold particles in the cytoplasm were seen as singles, as doublets, or in small groups of up to five particles probably representing small oligomers. Clusters of gold particles were also observed in acutely infected lymphocytes as well as in multinuclear cells of chronically infected cultures of T4 cells. In a few cases small aggregates of gold particles were found in the nuclei of T4 lymphocytes. These observations suggest that the Gag polyprotein forms small oligomers in the cytoplasm of expressing cells but that assembly into multimeric complexes takes place predominantly at the plasma membrane. Large accumulations of Gag protein in the cytoplasm may represent misfolded molecules destined for degradation.

Evidence for a reversal with age in the pattern of near-wall blood flow around aortic branches

The changes that occur with age in the distribution of atherosclerotic lesions around arterial branch points challenge accepted theories relating disease to haemodynamic stresses. We investigated whether flow near branch points changes with age in a way that can account for the different lesion distributions. Flow around 20 branches from immature and mature aortas was investigated by examining the length:width ratio and orientation of endothelial nuclei; these properties depend on the magnitude and direction of near-wall flows, respectively. There were significant changes in the pattern of nuclear shape with age, consistent with a reversal in the pattern of shear around branches. In control regions away from branches, there were no such changes. The role of haemodynamic stresses in atherogenesis may require re-evaluation in the light of these results.

13C-2H correlation NMR spectroscopy studies of the in vivo transformations of natural products from Artemisia annua

13C-2H correlation NMR spectroscopy (13C-2H COSY) permits the identification of 13C and 2H nuclei which are connected to one another by a single chemical bond via the sizeable 1JCD coupling constant. The practical development of this technique is described using a 13C-2H COSY pulse sequence which is derived from the classical 13C-1H correlation experiment. An example is given of the application of 13C-2H COSY to the study of the biogenesis of natural products from the anti-malarial plant Artemisia annua, using a doubly-labelled precursor molecule. Although the biogenesis of artemisinin, the anti-malarial principle from this species, has been extensively studied over the past twenty years there is still no consensus as to the true biosynthetic route to this important natural product – indeed, some published experimental results are directly contradictory. One possible reason for this confusion may be the ease with which some of the metabolites from A. annua undergo spontaneous autoxidation, as exemplified by our recent in vitro studies of the spontaneous autoxidation of dihydroartemisinic acid, and the application of 13C-2H COSY to this biosynthetic problem has been important in helping to mitigate against such processes. In this in vivo application of 13C-2H COSY, [15-13C2H3]-dihydroartemisinic acid (the doubly-labelled analogue of the natural product from this species which was obtained through synthesis) was fed to A. annua plants and was shown to be converted into several natural products which have been described previously, including artemisinin. It is proposed that all of these transformations occurred via a tertiary hydroperoxide intermediate, which is derived from dihyroartemisinic acid. This intermediate was observed directly in this feeding experiment by the 13C-2H COSY technique; its observation by more traditional procedures (e.g., chromatographic separation, followed by spectroscopic analysis of the purified product) would have been difficult owing to the instability of the hydroperoxide group (as had been established previously by our in vitro studies of the spontaneous autoxidation of dihydroartemisinic acid). This same hydroperoxide has been reported as the initial product of the spontaneous autoxidation of dihydroartemisinic acid in our previous in vitro studies. Its observation in this feeding experiment by the 13C-2H COSY technique, a procedure which requires the minimum of sample manipulation in order to achieve a reliable identification of metabolites (based on both 13C and 2H chemical shifts at the 15-position), provides the best possible evidence for its status as a genuine biosynthetic intermediate, rather than merely as an artifact of the experimental procedure.

Oxidation of oleic acid at the air-water interface and its potential effects on cloud critical supersaturations

The oxidation of organic films on cloud condensation nuclei has the potential to affect climate and precipitation events. In this work we present a study of the oxidation of a monolayer of deuterated oleic acid (cis-9-octadecenoic acid) at the air-water interface by ozone to determine if oxidation removes the organic film or replaces it with a product film. A range of different aqueous sub-phases were studied. The surface excess of deuterated material was followed by neutron reflection whilst the surface pressure was followed using a Wilhelmy plate. The neutron reflection data reveal that approximately half the organic material remains at the air-water interface following the oxidation of oleic acid by ozone, thus cleavage of the double bond by ozone creates one surface active species and one species that partitions to the bulk (or gas) phase. The most probable products, produced with a yield of similar to(87 +/- 14)%, are nonanoic acid, which remains at the interface, and azelaic acid (nonanedioic acid), which dissolves into the bulk solution. We also report a surface bimolecular rate constant for the reaction between ozone and oleic acid of (7.3 +/- 0.9) x 10(-11) cm(2) molecule s(-1). The rate constant and product yield are not affected by the solution sub-phase. An uptake coefficient of ozone on the oleic acid monolayer of similar to 4 x 10(-6) is estimated from our results. A simple Kohler analysis demonstrates that the oxidation of oleic acid by ozone on an atmospheric aerosol will lower the critical supersaturation needed for cloud droplet formation. We calculate an atmospheric chemical lifetime of oleic acid of 1.3 hours, significantly longer than laboratory studies on pure oleic acid particles suggest, but more consistent with field studies reporting oleic acid present in aged atmospheric aerosol.

Potential energy surface and MULTIMODE vibrational analysis of C2H3+

A full dimensional, ab initio-based semiglobal potential energy surface for C2H3+ is reported. The ab initio electronic energies for this molecule are calculated using the spin-restricted, coupled cluster method restricted to single and double excitations with triples corrections [RCCSD(T)]. The RCCSD(T) method is used with the correlation-consistent polarized valence triple-zeta basis augmented with diffuse functions (aug-cc-pVTZ). The ab initio potential energy surface is represented by a many-body (cluster) expansion, each term of which uses functions that are fully invariant under permutations of like nuclei. The fitted potential energy surface is validated by comparing normal mode frequencies at the global minimum and secondary minimum with previous and new direct ab initio frequencies. The potential surface is used in vibrational analysis using the "single-reference" and "reaction-path" versions of the code MULTIMODE. (c) 2006 American Institute of Physics.