Linear response functions for a vibrational configuration interaction state

Linear response functions are implemented for a vibrational configuration interaction state allowing accurate analytical calculations of pure vibrational contributions to dynamical polarizabilities. Sample calculations are presented for the pure vibrational contributions to the polarizabilities of water and formaldehyde. We discuss the convergence of the results with respect to various details of the vibrational wave function description as well as the potential and property surfaces. We also analyze the frequency dependence of the linear response function and the effect of accounting phenomenologically for the finite lifetime of the excited vibrational states. Finally, we compare the analytical response approach to a sum-over-states approach

Basis set superposition error-counterpoise corrected potential energy surfaces : application to hydrogen peroxide ... X (X=F-, Cl-, Br-, Li+, Na+) complexes

Møller-Plesset (MP2) and Becke-3-Lee-Yang-Parr (B3LYP) calculations have been used to compare the geometrical parameters, hydrogen-bonding properties, vibrational frequencies and relative energies for several X- and X+ hydrogen peroxide complexes. The geometries and interaction energies were corrected for the basis set superposition error (BSSE) in all the complexes (1-5), using the full counterpoise method, yielding small BSSE values for the 6-311 + G(3df,2p) basis set used. The interaction energies calculated ranged from medium to strong hydrogen-bonding systems (1-3) and strong electrostatic interactions (4 and 5). The molecular interactions have been characterized using the atoms in molecules theory (AIM), and by the analysis of the vibrational frequencies. The minima on the BSSE-counterpoise corrected potential-energy surface (PES) have been determined as described by S. Simón, M. Duran, and J. J. Dannenberg, and the results were compared with the uncorrected PES

C-H···O H-bonded complexes: how does basis set superposition error change their potential-energy surfaces?

Geometries, vibrational frequencies, and interaction energies of the CNH⋯O3 and HCCH⋯O3 complexes are calculated in a counterpoise-corrected (CP-corrected) potential-energy surface (PES) that corrects for the basis set superposition error (BSSE). Ab initio calculations are performed at the Hartree-Fock (HF) and second-order Møller-Plesset (MP2) levels, using the 6-31G(d,p) and D95++(d,p) basis sets. Interaction energies are presented including corrections for zero-point vibrational energy (ZPVE) and thermal correction to enthalpy at 298 K. The CP-corrected and conventional PES are compared; the unconnected PES obtained using the larger basis set including diffuse functions exhibits a double well shape, whereas use of the 6-31G(d,p) basis set leads to a flat single-well profile. The CP-corrected PES has always a multiple-well shape. In particular, it is shown that the CP-corrected PES using the smaller basis set is qualitatively analogous to that obtained with the larger basis sets, so the CP method becomes useful to correctly describe large systems, where the use of small basis sets may be necessary

How does basis set superposition error change the potential surfaces for hydrogen-bonded dimers?

We describe a simple method to automate the geometric optimization of molecular orbital calculations of supermolecules on potential surfaces that are corrected for basis set superposition error using the counterpoise (CP) method. This method is applied to the H-bonding complexes HF/HCN, HF/H2O, and HCCH/H2O using the 6-31G(d,p) and D95 + + (d,p) basis sets at both the Hartree-Fock and second-order Møller-Plesset levels. We report the interaction energies, geometries, and vibrational frequencies of these complexes on the CP-optimized surfaces; and compare them with similar values calculated using traditional methods, including the (more traditional) single point CP correction. Upon optimization on the CP-corrected surface, the interaction energies become more negative (before vibrational corrections) and the H-bonding stretching vibrations decrease in all cases. The extent of the effects vary from extremely small to quite large depending on the complex and the calculational method. The relative magnitudes of the vibrational corrections cannot be predicted from the H-bond stretching frequencies alone

Producció i caracterització de variants de la ribonucleasa pancreàtica humana dissenyades per a adquirir propietats citotòxiques

Amb la finalitat d'aprofundir en les bases moleculars de la citotoxicitat de les ribonucleases pancreàtiques, es van construir variants derivades de l'HP-RNasa seguint dues estratègies. En la primera, es van generar variants de l'enzim resistents a l'acció de l'inhibidor proteic de les ribonucleases (hRI), substituint residus implicats en la interfície de contacte entre la ribonucleasa i l'hRI. En la segona, es va addicionar el motiu RGD en regions de superfície de la proteïna implicades en la formació del complex amb l'hRI, a fi de promoure la seva interacció amb la membrana plasmàtica de les cèl·lules i a la vegada disminuir l'afinitat de les variants per l'hRI. Es va comprovar que només les variants portadores de substitucions múltiples adquirien la capacitat de resistència a l'hRI. L'estudi del percentatge d'inhibició de la síntesi proteica en cèl·lules incubades amb cadascuna de les variants va mostrar que només dues de les variants construïdes havien adquirit propietats citotòxiques. La citotoxicitat més elevada la va presentar una variant que no era resistent a l'hRI, amb valors que eren només entre 5 i 15 vegades inferiors als de l'onconasa. Aquest resultat demostrà que la sensibilitat a l'hRI no és necessàriament un paràmetre limitant per a la citotoxicitat de les ribonucleases. Cap de les variants que incorporava un motiu RGD presentà citotoxicitat, evidenciant que aquest motiu no és efectiu a fi de dotar les ribonucleases pancreàtiques de propietats citotòxiques. Es van estudiar les bases moleculars de la citotoxicitat de la variant més citotòxica. En primer lloc, l'anàlisi de la internalització per marcatge radioactiu d'aquesta variant en relació amb l'onconasa i amb altres variants de l'HP-RNasa no citotòxiques, va posar en evidència que només l'onconasa era internalitzada eficientment. Es descartava així la possibilitat que l'acció citotòxica de l'enzim estudiat fos conseqüència d'una major eficiència d'endocitosi. També es va comprovar que l'addició del motiu RGD no era capaç de promoure la internalització de les proteïnes amb més eficàcia. Per microscòpia confocal de fluorescència, les variants humanes només es van començar a detectar a l'interior de la cèl·lula a partir de les 24 h d'incubació. Totes les variants generades van presentar una eficiència catalítica superior al 50 % de l'activitat de la seva proteïna parental, PM5, indicant que probablement l'estructura del centre actiu no havia estat afectada de manera dràstica per les substitucions introduïdes. No obstant, en tots els casos es va produir una disminució en la termoestabilitat respecte a PM5. Aquest resultat indicà que la correlació descrita a la bibliografia entre l'increment de termoestabilitat i l'increment de citotoxicitat per les ribonucleases no sempre es compleix. Per microscòpia confocal es va comprovar que tant la proteïna més citotòxica, com una variant no citotòxica resistent a l'hRI, així com la proteïna parental, seguien la via de degradació lisosomal. Aquesta ruta de trànsit no va ser afectada per l'addició de drogues que alteren les vies de trànsit retrògrad (monensina i brefeldina A), però sí per l'addició de la bafilomicina A1, una droga que neutralitza el pH endosomal i que va actuar alentint el trànsit de les proteïnes als lisosomes. D'acord amb aquests resultats, els valors de citotoxicitat de les variants es van incrementar de manera significativa només en presència de bafilomicina A1, suggerint que les ribonucleases transloquen al citoplasma a partir d'algun punt de la via de trànsit endosomal. Es va comprovar que l'acció de la variant més citotòxica era deguda a que l'addició d'un segon motiu de tres Arg en PE5 dota a aquesta proteïna amb un senyal de transport nuclear. La fracció d'enzim que aconsegueix translocar al citoplasma a partir d'algun punt de la via endosomal previ als lisosomes, és conduït ràpidament al nucli de la cèl·lula per mitjà del mecanisme clàssic de transport actiu. Per la seva afinitat amb l'rRNA, l'enzim es concentra en el nuclèol, on probablement duu a terme la seva activitat catalítica. La interacció d'aquesta variant amb els receptors nucleocitoplasmàtics, les importines, impediria per altra banda el bloqueig de l'enzim per part de l'hRI. Els resultats obtinguts presenten una nova estratègia de disseny de ribonucleases citotòxiques, basada en l'addició de segments NLS a fi de promoure el transport nuclear dels enzims. Aquesta estratègia podria permetre superar limitacions que fins al moment han estat descrites com a limitants de la citotoxicitat de les ribonucleases pancreàtiques, com la sensibilitat a l'hRI o la baixa eficiència d'internalització.

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.