Molecular dissolved organic matter composition in lakes across Sweden as relative intensities of FT-ICR-MS peaks and PARAFAC components and optical indices

Whether intrinsic molecular properties or extrinsic factors such as environmental conditions control the decomposition of natural organic matter across soil, marine and freshwater systems has been subject to debate. Comprehensive evaluations of the controls that molecular structure exerts on organic matter's persistence in the environment have been precluded by organic matter's extreme complexity. Here we examine dissolved organic matter from 109 Swedish lakes using ultrahigh-resolution mass spectrometry and optical spectroscopy to investigate the constraints on its persistence in the environment. We find that degradation processes preferentially remove oxidized, aromatic compounds, whereas reduced, aliphatic and N-containing compounds are either resistant to degradation or tightly cycled and thus persist in aquatic systems. The patterns we observe for individual molecules are consistent with our measurements of emergent bulk characteristics of organic matter at wide geographic and temporal scales, as reflected by optical properties. We conclude that intrinsic molecular properties are an important control of overall organic matter reactivity.

Interaction of the Hsp70 molecular chaperone, DnaK, with its cochaperone DnaJ

Chaperones of the Hsp70 family bind to unfolded or partially folded polypeptides to facilitate many cellular processes. ATP hydrolysis and substrate binding, the two key molecular activities of this chaperone, are modulated by the cochaperone DnaJ. By using both genetic and biochemical approaches, we provide evidence that DnaJ binds to at least two sites on the Escherichia coli Hsp70 family member DnaK: under the ATPase domain in a cleft between its two subdomains and at or near the pocket of substrate binding. The lower cleft of the ATPase domain is defined as a binding pocket for the J-domain because (i) a DnaK mutation located in this cleft (R167H) is an allele-specific suppressor of the binding defect of the DnaJ mutation, D35N and (ii) alanine substitution of two residues close to R167 in the crystal structure, N170A and T173A, significantly decrease DnaJ binding. A second binding determinant is likely to be in the substrate-binding domain because some DnaK mutations in the vicinity of the substrate-binding pocket are defective in either the affinity (G400D, G539D) or rate (D526N) of both peptide and DnaJ binding to DnaK. Binding of DnaJ may propagate conformational changes to the nearby ATPase catalytic center and substrate-binding sites as well as facilitate communication between these two domains to alter the molecular properties of Hsp70.

X-ray studies on crystalline complexes involving amino acids and peptides. XXXIII. Crystal structures of L- and DL-arginine complexed with oxalic acid and a comparative study of amino acid oxalic acid complexes

The DL- and L-arginine complexes of oxalic acid are made up of zwitterionic positively charged amino acid molecules and semi-oxalate ions. The dissimilar molecules aggregate into separate alternating layers in the former. The basic unit in the arginine layer is a centrosymmetric dimer, while the semi-oxalate ions form hydrogen-bonded strings in their layer. In the L-arginine complex each semi-oxalate ion is surrounded by arginine molecules and the complex can be described as an inclusion compound. The oxalic acid complexes of basic amino acids exhibit a variety of ionization states and stoichiometry. They illustrate the effect of aggregation and chirality on ionization state and stoichiometry, and that of molecular properties on aggregation. The semi-oxalate/oxalate ions tend to be planar, but large departures from planarity are possible. The amino acid aggregation in the different oxalic acid complexes do not resemble one another significantly, but the aggregation of a particular amino acid in its oxalic acid complex tends to have similarities with its aggregation in other structures. Also, semi-oxalate ions aggregate into similar strings in four of the six oxalic acid complexes. Thus, the intrinsic aggregation propensities of individual molecules tend to be retained in the complexes.

Evolution of domain combinations in protein kinases and its implications for functional diversity

Protein kinases phosphorylating Ser/Thr/Tyr residues in several cellular proteins exert tight control over their biological functions. They constitute the largest protein family in most eukaryotic species. Protein kinases classified based on sequence similarity in their catalytic domains, cluster into subfamilies, which share gross functional properties. Many protein kinases are associated or tethered covalently to domains that serve as adapter or regulatory modules,naiding substrate recruitment, specificity, and also serve as scaffolds. Hence the modular organisation of the protein kinases serves as guidelines to their functional and molecular properties. Analysis of genomic repertoires of protein kinases in eukaryotes have revealed wide spectrum of domain organisation across various subfamilies of kinases. Occurrence of organism-specific novel domain combinations suggests functional diversity achieved by protein kinases in order to regulate variety of biological processes. In addition, domain architecture of protein kinases revealed existence of hybrid protein kinase subfamilies and their emerging roles in the signaling of eukaryotic organisms. In this review we discuss the repertoire of non-kinase domains tethered to multi-domain kinases in the metazoans. Similarities and differences in the domain architectures of protein kinases in these organisms indicate conserved and unique features that are critical to functional specialization. (C) 2009 Elsevier Ltd. All rights reserved.

Application of vibrational microspectroscopy to biology and medicine

Vibrational microspectroscopic (Raman and infrared (IR)) techniques are rapidly emerging as effective tools to probe the basic processes of life. This review mainly focuses on the applications of Raman and IR microspectroscopy to biology and biomedicine, ranging from studies on cellular components in single cells to advancement in techniques for in vitro to in vivo applications. These techniques have proved to be instrumental in studying the biological specimen with minimum perturbation, i.e. without the use of dyes and contrast-inducing agents. These techniques probe the vibrational modes of the molecules and provide spectra that are specific to the molecular properties and chemical nature of the species.

Effect of alkyl substituents in BODIPYs: a comparative DFT computational investigation

Random changes in the alkyl substitution patterns of fluorescent dyes, e.g. BODIPYs, are often accompanied by significant changes in their photophysical properties. To understand such alterations in properties in closely related molecular systems, a comparative DFT (density functional theory) computational investigation was performed in order to comprehend the effects of alkyl substitution in controlling the structural and electronic nature of BODIPY dyes. In this context, a systematic strategy was utilized, considering all possible outcomes of constitutionally-isomeric molecules to understand the alkyl groups' effects on the BODIPY molecules. Four different computational methods {i.e. B3LYP/631G(d); B3LYP/6-311++ G(d,p); wb97xd/6-311++ G(d,p) and mpw1pw91/6-311++ G(d,p)} were employed to rationalize the agreement of the trends associated with the molecular properties. In line with experimental observations, it was found that alkyl substituents in BODIPY dyes situated at 3/5-positions effectively participate in stabilization as well as planarization of such molecules. Screening of all the possible isomeric molecular systems was used to understand the individual properties and overall effects of the typical alkyl substituents in controlling several basic properties of such BODIPY molecules.

Investigation of Fundamental Processes Governing Secondary Organic Aerosol Formation in Laboratory Chambers

Our understanding of the processes and mechanisms by which secondary organic aerosol (SOA) is formed is derived from laboratory chamber studies. In the atmosphere, SOA formation is primarily driven by progressive photooxidation of SOA precursors, coupled with their gas-particle partitioning. In the chamber environment, SOA-forming vapors undergo multiple chemical and physical processes that involve production and removal via gas-phase reactions; partitioning onto suspended particles vs. particles deposited on the chamber wall; and direct deposition on the chamber wall. The main focus of this dissertation is to characterize the interactions of organic vapors with suspended particles and the chamber wall and explore how these intertwined processes in laboratory chambers govern SOA formation and evolution.

A Functional Group Oxidation Model (FGOM) that represents SOA formation and evolution in terms of the competition between functionalization and fragmentation, the extent of oxygen atom addition, and the change of volatility, is developed. The FGOM contains a set of parameters that are to be determined by fitting of the model to laboratory chamber data. The sensitivity of the model prediction to variation of the adjustable parameters allows one to assess the relative importance of various pathways involved in SOA formation.

A critical aspect of the environmental chamber is the presence of the wall, which can induce deposition of SOA-forming vapors and promote heterogeneous reactions. An experimental protocol and model framework are first developed to constrain the vapor-wall interactions. By optimal fitting the model predictions to the observed wall-induced decay profiles of 25 oxidized organic compounds, the dominant parameter governing the extent of wall deposition of a compound is identified, i.e., wall accommodation coefficient. By correlating this parameter with the molecular properties of a compound via its volatility, the wall-induced deposition rate of an organic compound can be predicted based on its carbon and oxygen numbers in the molecule.

Heterogeneous transformation of δ-hydroxycarbonyl, a major first-generation product from long-chain alkane photochemistry, is observed on the surface of particles and walls. The uniqueness of this reaction scheme is the production of substituted dihydrofuran, which is highly reactive towards ozone, OH, and NO3, thereby opening a reaction pathway that is not usually accessible to alkanes. A spectrum of highly-oxygenated products with carboxylic acid, ester, and ether functional groups is produced from the substituted dihydrofuran chemistry, thereby affecting the average oxidation state of the alkane-derived SOA.

The vapor wall loss correction is applied to several chamber-derived SOA systems generated from both anthropogenic and biogenic sources. Experimental and modeling approaches are employed to constrain the partitioning behavior of SOA-forming vapors onto suspended particles vs. chamber walls. It is demonstrated that deposition of SOA-forming vapors to the chamber wall during photooxidation experiments can lead to substantial and systematic underestimation of SOA. Therefore, it is likely that a lack of proper accounting for vapor wall losses that suppress chamber-derived SOA yields contribute substantially to the underprediction of ambient SOA concentrations in atmospheric models.

Carbon-13 magnetic resonance spectroscopy

High-resolution, natural-abundance ¹³C spectra have been obtained from a wide variety of organic compounds; ¹³C chemical shifts and coupling constants have been correlated with other molecular properties.

Geminal and vicinal, carbon-proton couplings in benzene and the five- and six-membered aromatic heterocycles have been related to the corresponding proton-proton couplings in substituted ethylenes. The carbon-proton coupling constants in benzene are J_CCH = + 1.0, J_CCCH = +7.4 and J_CCCH = -1.1 Hz. Extended Hückel wavefunctions are uniformly poor in explaining the long-range, carbon-proton couplings in aromatic systems.

Couplings between carbon and elements other than hydrogen have been observed in proton decoupled ¹³C spectra. All of the carbons in fluorobenzene and 1-fluoronaphthalene, but only six of the carbons in 2-fluoronaphthalene are coupled to the fluorine. One-bond, carbon-phosphorus coupling in trialkylphosphines is negative, while one-bond, carbon-phosphorus coupling in tetra-alkylphosphonium ions is positive. Atoms which do not use hybrid orbitals to form bonds to carbon (F, P(III), Se, Te) may have negative, one-bond coupling constants because of the failure of the average energy approximation. One-bond couplings between carbon and carbon, silicon, tin, lead and mercury appear to be explainable in terms of an effective nuclear charge and the s-bond order of the metal. Couplings between carbon and nitrogen and phosphorus (IV) have significant negative contributions to the Fermi contact coupling expression, though, within one series, correlations with s-bond order may be valid. Carbon-carbon coupling in cyclopropane derivatives (10-15 Hz) is consistent with a high degree of p character in the interior orbitals. Some two- and three-bond carbon-carbon coupling constants have also been observed.

Substituent effects of hydroxyl groups on the ¹³C chemical shifts of continuous-chain alkanes depend both on steric and electronic factors. The hydroxyl substituent effects in the long-chain, primary alcohols are α = -48.3, β = -10.2, and γ = +6.0 ppm. The upfield γ effect is attributed to steric crowding in the gauche conformations. Additivity of the hydroxyl and carbonyl and alkyl substituent effects in alkyl-substituted cyclohexanols and cyclohexanones has been demonstrated.

Chondrogenesis and mineralization during in vitro culture of human mesenchymal stem cells on three-dimensional woven scaffolds.

Human mesenchymal stem cells (hMSCs) and three-dimensional (3D) woven poly(ɛ-caprolactone) (PCL) scaffolds are promising tools for skeletal tissue engineering. We hypothesized that in vitro culture duration and medium additives can individually and interactively influence the structure, composition, mechanical, and molecular properties of engineered tissues based on hMSCs and 3D poly(ɛ-caprolactone). Bone marrow hMSCs were suspended in collagen gel, seeded on scaffolds, and cultured for 1, 21, or 45 days under chondrogenic and/or osteogenic conditions. Structure, composition, biomechanics, and gene expression were analyzed. In chondrogenic medium, cartilaginous tissue formed by day 21, and hypertrophic mineralization was observed in the newly formed extracellular matrix at the interface with underlying scaffold by day 45. Glycosaminoglycan, hydroxyproline, and calcium contents, and alkaline phosphatase activity depended on culture duration and medium additives, with significant interactive effects (all p < 0.0001). The 45-day constructs exhibited mechanical properties on the order of magnitude of native articular cartilage (aggregate, Young's, and shear moduli of 0.15, 0.12, and 0.033 MPa, respectively). Gene expression was characteristic of chondrogenesis and endochondral bone formation, with sequential regulation of Sox-9, collagen type II, aggrecan, core binding factor alpha 1 (Cbfα1)/Runx2, bone sialoprotein, bone morphogenetic protein-2, and osteocalcin. In contrast, osteogenic medium produced limited osteogenesis. Long-term culture of hMSC on 3D scaffolds resulted in chondrogenesis and regional mineralization at the interface between soft, newly formed engineered cartilage, and stiffer underlying scaffold. These findings merit consideration when developing grafts for osteochondral defect repair.

Volatile terpenoids and C13 norisoprenoids in wines: development of rapid methods of analysis and evaluation of the sesquiterpenoids biological potential

Vitis vinifera L., the most widely cultivated fruit crop in the world, was the starting point for the development of this PhD thesis. This subject was exploited following on two actual trends: i) the development of rapid, simple, and high sensitive methodologies with minimal sample handling; and ii) the valuation of natural products as a source of compounds with potential health benefits. The target group of compounds under study were the volatile terpenoids (mono and sesquiterpenoids) and C13 norisoprenoids, since they may present biological impact, either from the sensorial point of view, as regards to the wine aroma, or by the beneficial properties for the human health. Two novel methodologies for quantification of C13 norisoprenoids in wines were developed. The first methodology, a rapid method, was based on the headspace solid-phase microextraction combined with gas chromatography-quadrupole mass spectrometry operating at selected ion monitoring mode (HS-SPME/GC-qMS-SIM), using GC conditions that allowed obtaining a C13 norisoprenoid volatile signature. It does not require any pre-treatment of the sample, and the C13 norisoprenoid composition of the wine was evaluated based on the chromatographic profile and specific m/z fragments, without complete chromatographic separation of its components. The second methodology, used as reference method, was based on the HS-SPME/GC-qMS-SIM, allowing the GC conditions for an adequate chromatographic resolution of wine components. For quantification purposes, external calibration curves were constructed with β-ionone, with regression coefficient (r2) of 0.9968 (RSD 12.51 %) and 0.9940 (RSD of 1.08 %) for the rapid method and for the reference method, respectively. Low detection limits (1.57 and 1.10 μg L-1) were observed. These methodologies were applied to seventeen white and red table wines. Two vitispirane isomers (158-1529 L-1) and 1,1,6-trimethyl-1,2-dihydronaphthalene (TDN) (6.42-39.45 μg L-1) were quantified. The data obtained for vitispirane isomers and TDN using the two methods were highly correlated (r2 of 0.9756 and 0.9630, respectively). A rapid methodology for the establishment of the varietal volatile profile of Vitis vinifera L. cv. 'Fernão-Pires' (FP) white wines by headspace solid-phase microextraction combined with comprehensive two-dimensional gas chromatography with time-of-flight mass spectrometry (HS-SPME/GCxGC-TOFMS) was developed. Monovarietal wines from different harvests, Appellations, and producers were analysed. The study was focused on the volatiles that seem to be significant to the varietal character, such as mono and sesquiterpenic compounds, and C13 norisoprenoids. Two-dimensional chromatographic spaces containing the varietal compounds using the m/z fragments 93, 121, 161, 175 and 204 were established as follows: 1tR = 255-575 s, 2tR = 0,424-1,840 s, for monoterpenoids, 1tR = 555-685 s, 2tR = 0,528-0,856 s, for C13 norisoprenoids, and 1tR = 695-950 s, 2tR = 0,520-0,960 s, for sesquiterpenic compounds. For the three chemical groups under study, from a total of 170 compounds, 45 were determined in all wines, allowing defining the "varietal volatile profile" of FP wine. Among these compounds, 15 were detected for the first time in FP wines. This study proposes a HS-SPME/GCxGC-TOFMS based methodology combined with classification-reference sample to be used for rapid assessment of varietal volatile profile of wines. This approach is very useful to eliminate the majority of the non-terpenic and non-C13 norisoprenic compounds, allowing the definition of a two-dimensional chromatographic space containing these compounds, simplifying the data compared to the original data, and reducing the time of analysis. The presence of sesquiterpenic compounds in Vitis vinifera L. related products, to which are assigned several biological properties, prompted us to investigate the antioxidant, antiproliferative and hepatoprotective activities of some sesquiterpenic compounds. Firstly, the antiradical capacity of trans,trans-farnesol, cis-nerolidol, α-humulene and guaiazulene was evaluated using chemical (DPPH• and hydroxyl radicals) and biological (Caco-2 cells) models. Guaiazulene (IC50= 0.73 mM) was the sesquiterpene with higher scavenger capacity against DPPH•, while trans,trans-farnesol (IC50= 1.81 mM) and cis-nerolidol (IC50= 1.48 mM) were more active towards hydroxyl radicals. All compounds, with the exception of α-humulene, at non-cytotoxic levels (≤ 1 mM), were able to protect Caco-2 cells from oxidative stress induced by tert-butyl hydroperoxide. The activity of the compounds under study was also evaluated as antiproliferative agents. Guaiazulene and cis-nerolidol were able to more effectively arrest the cell cycle in the S-phase than trans,trans-farnesol and α-humulene, being the last almost inactive. The relative hepatoprotection effect of fifteen sesquiterpenic compounds, presenting different chemical structures and commonly found in plants and plant-derived foods and beverages, was assessed. Endogenous lipid peroxidation and induced lipid peroxidation with tert-butyl hydroperoxide were evaluated in liver homogenates from Wistar rats. With the exception of α-humulene, all the sesquiterpenic compounds under study (1 mM) were effective in reducing the malonaldehyde levels in both endogenous and induced lipid peroxidation up to 35% and 70%, respectively. The developed 3D-QSAR models, relating the hepatoprotection activity with molecular properties, showed good fit (R2LOO > 0.819) with good prediction power (Q2 > 0.950 and SDEP < 2%) for both models. A network of effects associated with structural and chemical features of sesquiterpenic compounds such as shape, branching, symmetry, and presence of electronegative fragments, can modulate the hepatoprotective activity observed for these compounds. In conclusion, this study allowed the development of rapid and in-depth methods for the assessment of varietal volatile compounds that might have a positive impact on sensorial and health attributes related to Vitis vinifera L. These approaches can be extended to the analysis of other related food matrices, including grapes and musts, among others. In addition, the results of in vitro assays open a perspective for the promising use of the sesquiterpenic compounds, with similar chemical structures such as those studied in the present work, as antioxidants, hepatoprotective and antiproliferative agents, which meets the current challenges related to diseases of modern civilization.

Characterisation of ZG16p, a unique mammalian lectin from pancreatic zymogen granules

The mechanisms of secretory granule biogenesis and regulated secretion of digestive enzymes in pancreatic acinar cells are still not well understood. To shed light on these processes, which are of biological and clinical importance (e.g., pancreatitis), a better molecular understanding of the components of the granule membrane, their functions and interactions is required. The application of proteomics has largely contributed to the identification of novel zymogen granule (ZG) proteins but was not yet accompanied by a better characterization of their functions. In this study we aimed at a) isolation and identification of novel membrane-associated ZG proteins; b) characterization of the biochemical properties and function of the secretory lectin ZG16p, a membrane-associated protein; c) exploring the potential of ZG16p as a new tool to label the endolysosomal compartment. First, we have performed a suborganellar proteomics approach by combining protein analysis by 2D-PAGE and identification by mass spectrometry, which has led to the identification of novel peripheral ZGM proteins with proteoglycan-binding properties (e.g., chymase, PpiB). Then, we have unveiled new molecular properties and (multiple) functions of the secretory lectin ZG16p. ZG16p is a unique mammalian lectin with glycan and proteoglycan binding properties. Here, I revealed for the first time that ZG16p is highly protease resistant by developing an enterokinase-digestion assay. In addition I revealed that ZG16p binds to a high molecular weight complex at the ZGM (which is also protease resistant) and forms highly stable dimers. In light of these findings I suggest that ZG16p is a key component of a predicted submembranous granule matrix attached to the luminal side of the ZGM that fulfils important functions during sorting and packaging of zymogens. ZG16p, may act as a linker between the matrix and aggregated zymogens due to dimer formation. Furthermore, ZG16p protease resistance might be of higher importance after secretion since it is known that ZG16p binds to pathogenic fungi in the gut. I have further investigated the role of ZG16p binding motifs in its targeting to ZG in AR42J cells, a pancreatic model system. Point mutations of the glycan and the proteoglycan binding motifs did not inhibit the targeting of ZG16p to ZG in AR42J cells. I have also demonstrated that when ZG16p is present in the cytoplasm it interacts with and modulates the endo-lysosomal compartment. Since it is known that impaired autophagy due to lysosomal malfunction is involved in the course of pancreatitis, a potential role of ZG16p in pancreatitis is discussed.

Optimization of the valence energy variance of the CuH molecule

We developed the concept of split-'t to deal with the large molecules (in terms of the number of electrons and nuclear charge Z). This naturally leads to partitioning the local energy into components due to each electron shell. The minimization of the variation of the valence shell local energy is used to optimize a simple two parameter CuH wave function. Molecular properties (spectroscopic constants and the dipole moment) are calculated for the optimized and nearly optimized wave functions using the Variational Quantum Monte Carlo method. Our best results are comparable to those from the single and double configuration interaction (SDCI) method.

Desenvolupament, implementació i aplicació de noves metodologies per a l'estudi de la reactivitat química i de les propietats òptiques no lineals

El coneixement de la superfície d'energia potencial (PES) ha estat essencial en el món de la química teòrica per tal de discutir tant la reactivitat química com l'estructura i l'espectroscòpia molecular. En el camp de la reactivitat química es hem proposat continuar amb el desenvolupament de nova metodologia dins el marc de la teoria del funcional de la densitat conceptual. En particular aquesta tesis es centrarà en els següents punts: a) El nombre i la naturalesa dels seus punts estacionaris del PES poden sofrir canvis radicals modificant el nivell de càlcul utilitzats, de tal manera que per estar segurs de la seva naturalesa cal anar a nivells de càlcul molt elevats. La duresa és una mesura de la resistència d'un sistema químic a canviar la seva configuració electrònica, i segons el principi de màxima duresa on hi hagi un mínim o un màxim d'energia trobarem un màxim o un mínim de duresa, respectivament. A l'escollir tot un conjunt de reaccions problemàtiques des del punt de vista de presència de punts estacionaris erronis, hem observat que els perfils de duresa són més independents de la base i del mètode utilitzats, a més a més sempre presenten el perfil correcte. b) Hem desenvolupat noves expressions basades en les integracions dels kernels de duresa per tal de determinar la duresa global d'una molècula de manera més precisa que la utilitzada habitualment que està basada en el càlcul numèric de la derivada segona de l'energia respecte al número d'electrons. c) Hem estudiat la validesa del principis de màxima duresa i de mínima polaritzabiliat en les vibracions asimètriques en sistemes aromàtics. Hem trobat que per aquests sistemes alguns modes vibracionals incompleixen aquests principis i hem analitzat la relació d'aquest l'incompliment amb l'efecte de l'acoblament pseudo-Jahn-Teller. A més a més, hem postulat tot un conjunt de regles molt senzilles que ens permetien deduir si una molècula compliria o no aquests principis sense la realització de cap càlcul previ. Tota aquesta informació ha estat essencial per poder determinar exactament quines són les causes del compliment o l'incompliment del MHP i MPP. d) Finalment, hem realitzat una expansió de l'energia funcional en termes del nombre d'electrons i de les coordenades normals dintre del conjunt canònic. En la comparació d'aquesta expansió amb l'expansió de l'energia del nombre d'electrons i del potencial extern hem pogut recuperar d'una altra forma diferent tot un conjunt de relacions ja conegudes entre alguns coneguts descriptors de reactivitat del funcional de la densitat i en poden establir tot un conjunt de noves relacions i de nous descriptors. Dins del marc de les propietats moleculars es proposa generalitzar i millorar la metodologia pel càlcul de la contribució vibracional (Pvib) a les propietats òptiques no lineals (NLO). Tot i que la Pvib no s'ha tingut en compte en la majoria dels estudis teòrics publicats de les propietats NLO, recentment s'ha comprovat que la Pvib de diversos polímers orgànics amb altes propietats òptiques no lineals és fins i tot més gran que la contribució electrònica. Per tant, tenir en compte la Pvib és essencial en el disseny dels nous materials òptics no lineals utilitzats en el camp de la informàtica, les telecomunicacions i la tecnologia làser. Les principals línies d'aquesta tesis sobre aquest tema són: a) Hem calculat per primera vegada els termes d'alt ordre de Pvib de diversos polímers orgànics amb l'objectiu d'avaluar la seva importància i la convergència de les sèries de Taylor que defineixen aquestes contribucions vibracionals. b) Hem avaluat les contribucions electròniques i vibracionals per una sèrie de molècules orgàniques representatives utilitzant diferents metodologies, per tal de poder de determinar quina és la manera més senzilla per poder calcular les propietats NLO amb una precisió semiquantitativa.

Implementation and application of basis set superposition error-correction schemes to the theoretical modeling of weak intermolecular interactions

This thesis deals with the so-called Basis Set Superposition Error (BSSE) from both a methodological and a practical point of view. The purpose of the present thesis is twofold: (a) to contribute step ahead in the correct characterization of weakly bound complexes and, (b) to shed light the understanding of the actual implications of the basis set extension effects in the ab intio calculations and contribute to the BSSE debate. The existing BSSE-correction procedures are deeply analyzed, compared, validated and, if necessary, improved. A new interpretation of the counterpoise (CP) method is used in order to define counterpoise-corrected descriptions of the molecular complexes. This novel point of view allows for a study of the BSSE-effects not only in the interaction energy but also on the potential energy surface and, in general, in any property derived from the molecular energy and its derivatives A program has been developed for the calculation of CP-corrected geometry optimizations and vibrational frequencies, also using several counterpoise schemes for the case of molecular clusters. The method has also been implemented in Gaussian98 revA10 package. The Chemical Hamiltonian Approach (CHA) methodology has been also implemented at the RHF and UHF levels of theory for an arbitrary number interacting systems using an algorithm based on block-diagonal matrices. Along with the methodological development, the effects of the BSSE on the properties of molecular complexes have been discussed in detail. The CP and CHA methodologies are used for the determination of BSSE-corrected molecular complexes properties related to the Potential Energy Surfaces and molecular wavefunction, respectively. First, the behaviour of both BSSE-correction schemes are systematically compared at different levels of theory and basis sets for a number of hydrogen-bonded complexes. The Complete Basis Set (CBS) limit of both uncorrected and CP-corrected molecular properties like stabilization energies and intermolecular distances has also been determined, showing the capital importance of the BSSE correction. Several controversial topics of the BSSE correction are addressed as well. The application of the counterpoise method is applied to internal rotational barriers. The importance of the nuclear relaxation term is also pointed out. The viability of the CP method for dealing with charged complexes and the BSSE effects on the double-well PES blue-shifted hydrogen bonds is also studied in detail. In the case of the molecular clusters the effect of high-order BSSE effects introduced with the hierarchical counterpoise scheme is also determined. The effect of the BSSE on the electron density-related properties is also addressed. The first-order electron density obtained with the CHA/F and CHA/DFT methodologies was used to assess, both graphically and numerically, the redistribution of the charge density upon BSSE-correction. Several tools like the Atoms in Molecules topologycal analysis, density difference maps, Quantum Molecular Similarity, and Chemical Energy Component Analysis were used to deeply analyze, for the first time, the BSSE effects on the electron density of several hydrogen bonded complexes of increasing size. The indirect effect of the BSSE on intermolecular perturbation theory results is also pointed out It is shown that for a BSSE-free SAPT study of hydrogen fluoride clusters, the use of a counterpoise-corrected PES is essential in order to determine the proper molecular geometry to perform the SAPT analysis.

P-found: grid-enabling distributed repositories of protein folding and unfolding simulations for data mining

The P-found protein folding and unfolding simulation repository is designed to allow scientists to perform data mining and other analyses across large, distributed simulation data sets. There are two storage components in P-found: a primary repository of simulation data that is used to populate the second component, and a data warehouse that contains important molecular properties. These properties may be used for data mining studies. Here we demonstrate how grid technologies can support multiple, distributed P-found installations. In particular, we look at two aspects: firstly, how grid data management technologies can be used to access the distributed data warehouses; and secondly, how the grid can be used to transfer analysis programs to the primary repositories — this is an important and challenging aspect of P-found, due to the large data volumes involved and the desire of scientists to maintain control of their own data. The grid technologies we are developing with the P-found system will allow new large data sets of protein folding simulations to be accessed and analysed in novel ways, with significant potential for enabling scientific discovery.