Pressure-jump-induced kinetics reveals a hydration dependent folding/unfolding mechanism of ribonuclease A

Pressure-jump (p-jump)-induced relaxation kinetics was used to explore the energy landscape of protein folding/unfolding of Y115W, a fluorescent variant of ribonuclease A. Pressure-jumps of 40MPa amplitude (5ms dead-time) were conducted both to higher (unfolding) and to lower (folding) pressure, in the range from 100 to 500MPa, between 30 and 50°C. Significant deviations from the expected symmetrical protein relaxation kinetics were observed. Whereas downward p-jumps resulted always in single exponential kinetics, the kinetics induced by upward p-jumps were biphasic in the low pressure range and monophasic at higher pressures. The relative amplitude of the slow phase decreased as a function of both pressure and temperature. At 50°C, only the fast phase remained. These results can be interpreted within the framework of a two-dimensional energy surface containing a pressure- and temperature-dependent barrier between two unfolded states differing in the isomeric state of the Asn-113–Pro-114 bond. Analysis of the activation volume of the fast kinetic phase revealed a temperature-dependent shift of the unfolding transition state to a larger volume. The observed compensation of this effect by glycerol offers an explanation for its protein stabilizing effect

Paper del nucli hidrofòbic principal de la RNasa A en el plegament i desplegament induïts per pressió i temperatura

Utilitzant temperatura i pressió com a agents desnaturalitzants s'ha explorat la contribució a l'estabilitat de diferents residus del principal nucli hidrofòbic de la RNasa A. Aquests resutats suggereixen que el principal nucli hidrofòbic d'aquest enzim, està fortament empaquetat i ha revelat l'existència de reordenacions en l'interior de la proteïna. El mètode dels valors , han permès estudiar el paper de les interaccions hidrofòbiques establertes pels residus del principal nucli hidrofòbic de la RNasa A en el seu estat de transició induït per pressió. En conjunt, aquests resultats suggereixen que l'estat de transició de la RNasa A, s'assemblaria a un glòbul col·lapsat amb una cadena estructurada però amb un debilitat nucli hidrofòbic. S'ha explorat també, el paisatge energètic del plegament/desplegament proteic de la variant Y115W de la RNasa A. L'estat de transició sembla interaccionar fortament amb la capa d'hidratació d'aquest estat, tal i com indiquen els resultats en presència de glicerol.

Rhodium(I) catalyzed [2+2+2] cycloaddition reactions: experimental and theoretical studies

The [2+2+2] cycloaddition reaction involves the formation of three carbon-carbon bonds in one single step using alkynes, alkenes, nitriles, carbonyls and other unsaturated reagents as reactants. This is one of the most elegant methods for the construction of polycyclic aromatic compounds and heteroaromatic, which have important academic and industrial uses. The thesis is divided into ten chapters including six related publications. The first study based on the Wilkinson’s catalyst, RhCl(PPh3)3, compares the reaction mechanism of the [2+2+2] cycloaddition process of acetylene with the cycloaddition obtained for the model of the complex, RhCl(PH3)3. In an attempt to reduce computational costs in DFT studies, this research project aimed to substitute PPh3 ligands for PH3, despite the electronic and steric effects produced by PPh3 ligands being significantly different to those created by PH3 ones. In this first study, detailed theoretical calculations were performed to determine the reaction mechanism of the two complexes. Despite some differences being detected, it was found that modelling PPh3 by PH3 in the catalyst helps to reduce the computational cost significantly while at the same time providing qualitatively acceptable results. Taking into account the results obtained in this earlier study, the model of the Wilkinson’s catalyst, RhCl(PH3)3, was applied to study different [2+2+2] cycloaddition reactions with unsaturated systems conducted in the laboratory. Our research group found that in the case of totally closed systems, specifically 15- and 25-membered azamacrocycles can afford benzenic compounds, except in the case of 20-membered azamacrocycle (20-MAA) which was inactive with the Wilkinson’s catalyst. In this study, theoretical calculations allowed to determine the origin of the different reactivity of the 20-MAA, where it was found that the activation barrier of the oxidative addition of two alkynes is higher than those obtained for the 15- and 25-membered macrocycles. This barrier was attributed primarily to the interaction energy, which corresponds to the energy that is released when the two deformed reagents interact in the transition state. The main factor that helped to provide an explanation to the different reactivity observed was that the 20-MAA had a more stable and delocalized HOMO orbital in the oxidative addition step. Moreover, we observed that the formation of a strained ten-membered ring during the cycloaddition of 20-MAA presents significant steric hindrance. Furthermore, in Chapter 5, an electrochemical study is presented in collaboration with Prof. Anny Jutand from Paris. This work allowed studying the main steps of the catalytic cycle of the [2+2+2] cycloaddition reaction between diynes with a monoalkyne. First kinetic data were obtained of the [2+2+2] cycloaddition process catalyzed by the Wilkinson’s catalyst, where it was observed that the rate-determining step of the reaction can change depending on the structure of the starting reagents. In the case of the [2+2+2] cycloaddition reaction involving two alkynes and one alkene in the same molecule (enediynes), it is well known that the oxidative coupling may occur between two alkynes giving the corresponding metallacyclopentadiene, or between one alkyne and the alkene affording the metallacyclopentene complex. Wilkinson’s model was used in DFT calculations to analyze the different factors that may influence in the reaction mechanism. Here it was observed that the cyclic enediynes always prefer the oxidative coupling between two alkynes moieties, while the acyclic cases have different preferences depending on the linker and the substituents used in the alkynes. Moreover, the Wilkinson’s model was used to explain the experimental results achieved in Chapter 7 where the [2+2+2] cycloaddition reaction of enediynes is studied varying the position of the double bond in the starting reagent. It was observed that enediynes type yne-ene-yne preferred the standard [2+2+2] cycloaddition reaction, while enediynes type yne-yne-ene suffered β-hydride elimination followed a reductive elimination of Wilkinson’s catalyst giving cyclohexadiene compounds, which are isomers from those that would be obtained through standard [2+2+2] cycloaddition reactions. Finally, the last chapter of this thesis is based on the use of DFT calculations to determine the reaction mechanism when the macrocycles are treated with transition metals that are inactive to the [2+2+2] cycloaddition reaction, but which are thermally active leading to new polycyclic compounds. Thus, a domino process was described combining an ene reaction and a Diels-Alder cycloaddition.

An investigation of the germylene addition reaction, GeH2+C2H2: Time-resolved gas-phase kinetic studies and quantum chemical calculations of the reaction energy surface

Time resolved studies of germylene, GeH2, generated by laser flash photolysis of 3,4-dimethylgermacyclopentene-3, have been carried out to obtain rate constants for its bimolecular reaction with acetylene, C2H2. The reaction was studied in the gas-phase over the pressure range 1-100 Tort, with SF6 as bath gas, at 5 temperatures in the range 297-553 K. The reaction showed a very slight pressure dependence at higher temperatures. The high pressure rate constants (obtained by extrapolation at the three higher temperatures) gave the Arrhenius equation: log(k(infinity)/cm(3) molecule(-1) s(-1)) (-10.94 +/- 0.05) + (6.10 +/- 0.36 kJ mol(-1))/RTln10. These Arrhenius parameters are consistent with a fast reaction occurring at approximately 30% of the collision rate at 298 K. Quantum chemical calculations (both DFT and ab initio G2//B3LYP and G2//QCISD) of the GeC2H4 potential energy surface (PES), show that GeH2 + C2H2 react initially to form germirene which can isomerise to vinylgermylene with a relatively low barrier. RRKM modelling, based on a loose association transition state, but assuming vinylgermylene is the end product (used in combination with a weak collisional deactivation model) predicts a strong pressure dependence using the calculated energies, in conflict with the experimental evidence. The detailed GeC2H4 PES shows considerable complexity with ten other accessible stable minima (B3LYP level), the three most stable of which are all germylenes. Routes through this complex surface were examined in detail. The only product combination which appears capable of satisfying the (P-3) + C2H4.C2H4 was confirmed as a product by GC observed lack of a strong pressure dependence is Ge(P-3) + C2H4. C2H4 was confirmed as a product by GC analysis. Although the formation of these products are shown to be possible by singlet-triplet curve crossing during dissociation of 1-germiranylidene (1-germacyclopropylidene), it seems more likely (on thermochernical grounds) that the triplet biradical, (GeCH2CH2.)-Ge-., is the immediate product precursor. Comparisons are made with the reaction of SiH2 with C2H2.

Investigation of the prototype silylene reaction, SiH2+H2O(and D2O): time-resolved gas-phase kinetic studies, isotope effects, RRKM calculations, and quantum chemical calculations of the reaction energy surface

Time-resolved kinetic studies of the reaction of silylene, SiH2, with H2O and with D2O have been carried out in the gas phase at 296 and at 339 K, using laser flash photolysis to generate and monitor SiH2. The reaction was studied over the pressure range 10-200 Torr with SF6 as bath gas. The second-order rate constants obtained were pressure dependent, indicating that the reaction is a third-body assisted association process. Rate constants at 339 K were about half those at 296 K. Isotope effects, k(H)/k(D), were small averaging 1.076 0.080, suggesting no involvement of H- (or D-) atom transfer in the rate determining step. RRKM modeling was undertaken based on a transition state appropriate to formation of the expected zwitterionic donoracceptor complex, H2Si...OH2. Because the reaction is close to the low pressure (third order) region, it is difficult to be definitive about the activated complex structure. Various structures were tried, both with and without the incorporation of rotational modes, leading to values for the high-pressure limiting (i.e., true secondorder) rate constant in the range 9.5 x 10(-11) to 5 x 10(-10) cm(3) molecule' s(-1). The RRKM modeling and mechanistic interpretation is supported by ab initio quantum calculations carried out at the G2 and G3 levels. The results are compared and contrasted with the previous studies.

endo-Selective (3+2) cycloaddition polymerizations of nitrone monomers with olefins utilising high pressure conditions

A high pressure mediated (3+2) cycloaddition polymerization strategy has been employed to afford linear poly(isoxazolidine) architectures. Under these high pressure conditions this cycloaddition process was found to afford primarily endoheterocycles which when translated to the polymerization should ultimately affect the tacticity and resultant properties of the polymer. The stereoselectivity occurred as a result of a lower volume of activation for the endo-transition state and the application of a 'type-I' regime (HOMODipole-LUMODipolarophile) cycloaddition process that features secondary orbital interactions within the extended molecular orbitals. A variety of linker segments were employed in an attempt to affect the physical properties of the polymeric cycloadducts such as T-g and solubility in order to tailor these materials for use in coating applications. (C) 2007 Elsevier Ltd. All rights reserved.

A DFT study of propylene polymerization using neutral salicyladiminato nickel(II) and palladium(II) as catalysts

Propylene polymerization using salicyladiminato metal catalalysts has been studied using density functional theory at the B3LYP/LANL2DZ level. In particular, the effects on the reaction mechanisms of changing the metal from Pd(II) to Ni(II) have been investigated. While the reaction mechanisms involving the salicyladiminato Ni(II) catalyst have been found to be similar to those established previously for the salicyladiminato Pd(II) catalyst, the nickel catalyst was found to differentiate the trans-O intermediate from the trans-.N intermediate with an energy difference of 46.63 U mol(-1) significantly more than the palladium catalyst for which the energy difference was calculated as 35.82 kJ mol(-1). The energy difference between the trans-O configuration and the trans-N configuration is decreased significantly when combining a molecule of propylene with the catalyst. For the Ni catalyst, the trans-O isomer is more stable than the trans-N isomer to a greater extent than for Pd, so that the insertion of propylene from 20 is relatively less favoured for Ni than for Pd. It is predicted that the mechanism of isomerization from 20 to 2N through a rotational transition state TS2O2N is more appropriate for the Ni catalyst system. The palladium system shows a larger preference for pi-coordination than its nickel counterpart, although the latter possesses a lower reaction barrier. It was found that the occupation of the trans-O position in the asymmetric salicyladiminato catalyst is also more favored by the alkene as it is by the alkyl so that insertion of the alkene may always start from a particular configuration so that specific products are obtained. (c) 2005 Elsevier B.V. All rights reserved.

2,2-Dialkyl-2H-benzimidazoles, the high energy tautomers of the corresponding 1,2-dialkyl-1H-benzimidazoles. Syntheses and their complexes with Cu(I) and Ag(I)

Reaction of Cu(1,2-phenylenediamine)(2)(ClO4)(2) with neat RR'=O (R = methyl and/or ethyl) (lives Cu(2,2-dialkyl-2H-benzimidazole)ClO4. demetallation of which by the action of aqueous ammonia yields Pure 2,2-dialkyl-2H-benzimidazoles. These are characterised by NMR. hi the X-ray crystal Structure, Ag(2,2-methyl-2H-benzimi-dazolc)NO3 is Found to be a spiral 1D coordination polymer where the 2H-benzimidazole acts as an N,N bridge between two Ag(I) centus. Although 2H-benzimidazoles are very unstable in the free state, they are quite stable in their Cu(I)(1) and Ag(I) complexes. The 1,2-tautomerisation in imidazole and benzimidazole have been Studied by means of transition state calculations at B3LYP/6-3 11 +G(2d,p)* level.

A low-temperature ammonium ylid rearrangement: enhanced reactivity engendered by rigidity

Sigmatropic rearrangement of tetrahydropyridine-derived ammonium is a valuable method for the preparation of substituted prolines. These reaction normally require elevated temperatures to proceed, but bicyclic tetrahydropyridine-like ylid I undergoes rearrangement at -15 degrees C; the extra rigidity of the azabicyclo[3.3.0]octene system preorganizes the transition state and lowers the activation energy for rearrangement.

Quantitative model for the kinetics of lyotropic phase transitions involving changes in monolayer curvature

We present a kinetic model for transformations between different self-assembled lipid structures. The model shows how data on the rates of phase transitions between mesophases of different geometries can be used to provide information on the mechanisms of the transformations and the transition states involved. This can be used, for example, to gain an insight into intermediate structures in cell membrane fission or fusion. In cases where the monolayer curvature changes on going from the initial to the final mesophase, we consider the phase transition to be driven primarily by the change in the relaxed curvature with pressure or temperature, which alters the relative curvature elastic energies of the two mesophase structures. Using this model, we have analyzed previously published kinetic data on the inter-conversion of inverse bicontinuous cubic phases in the 1-monoolein-30 wt% water system. The data are for a transition between QII(G) and QII(D) phases, and our analysis indicates that the transition state more closely resembles the QII(D) than the QII(G) phase. Using estimated values for the monolayer mean curvatures of the QII(G) and QII(D) phases of -0.123 nm(-1) and -0.133 nm(-1), respectively, gives values for the monolayer mean curvature of the transition state of between -0.131 nm(-1) and -0.132 nm(-1). Furthermore, we estimate that several thousand molecules undergo the phase transition cooperatively within one "cooperative unit", equivalent to 1-2 unit cells of QII(G) or 4-10 unit cells of QII(D).

A reaction surface Hamiltonian study of malonaldehyde

We report calculations using a reaction surface Hamiltonian for which the vibrations of a molecule are represented by 3N-8 normal coordinates, Q, and two large amplitude motions, s(1) and s(2). The exact form of the kinetic energy operator is derived in these coordinates. The potential surface is first represented as a quadratic in Q, the coefficients of which depend upon the values of s(1),s(2) and then extended to include up to Q(6) diagonal anharmonic terms. The vibrational energy levels are evaluated by solving the variational secular equations, using a basis of products of Hermite polynomials and appropriate functions of s(1),s(2). Our selected example is malonaldehyde (N=9) and we choose as surface parameters two OH distances of the migrating H in the internal hydrogen transfer. The reaction surface Hamiltonian is ideally suited to the study of the kind of tunneling dynamics present in malonaldehyde. Our results are in good agreement with previous calculations of the zero point tunneling splitting and in general agreement with observed data. Interpretation of our two-dimensional reaction surface states suggests that the OH stretching fundamental is incorrectly assigned in the infrared spectrum. This mode appears at a much lower frequency in our calculations due to substantial transition state character. (c) 2006 American Institute of Physics.

Evaluation of data for atmospheric models: master equation/RRKM calculations on the combination reaction, BrO+NO2 -> BrONO2, a conundrum

Experimental data for the title reaction were modeled using master equation (ME)/RRKM methods based on the Multiwell suite of programs. The starting point for the exercise was the empirical fitting provided by the NASA (Sander, S. P.; Finlayson-Pitts, B. J.; Friedl, R. R.; Golden, D. M.; Huie, R. E.; Kolb, C. E.; Kurylo, M. J.; Molina, M. J.; Moortgat, G. K.; Orkin, V. L.; Ravishankara, A. R. Chemical Kinetics and Photochemical Data for Use in Atmospheric Studies, Evaluation Number 15; Jet Propulsion Laboratory: Pasadena, California, 2006)(1) and IUPAC (Atkinson, R.; Baulch, D. L.; Cox, R. A.: R. F. Hampson, J.; Kerr, J. A.; Rossi, M. J.; Troe, J. J. Phys. Chem. Ref. Data. 2000, 29, 167) 2 data evaluation panels, which represents the data in the experimental pressure ranges rather well. Despite the availability of quite reliable parameters for these calculations (molecular vibrational frequencies (Parthiban, S.; Lee, T. J. J. Chem. Phys. 2000, 113, 145)3 and a. value (Orlando, J. J.; Tyndall, G. S. J. Phys. Chem. 1996, 100,. 19398)4 of the bond dissociation energy, D-298(BrO-NO2) = 118 kJ mol(-1), corresponding to Delta H-0(circle) = 114.3 kJ mol(-1) at 0 K) and the use of RRKM/ME methods, fitting calculations to the reported data or the empirical equations was anything but straightforward. Using these molecular parameters resulted in a discrepancy between the calculations and the database of rate constants of a factor of ca. 4 at, or close to, the low-pressure limit. Agreement between calculation and experiment could be achieved in two ways, either by increasing Delta H-0(circle) to an unrealistically high value (149.3 kJ mol(-1)) or by increasing , the average energy transferred in a downward collision, to an unusually large value (> 5000 cm(-1)). The discrepancy could also be reduced by making all overall rotations fully active. The system was relatively insensitive to changing the moments of inertia in the transition state to increase the centrifugal effect. The possibility of involvement of BrOONO was tested and cannot account for the difficulties of fitting the data.

Ferric and ferrous iron in nitroso-myoglobin: computer simulations of stable and metastable states and their infrared spectra

The binding of NO to iron is involved in the biological function of many heme proteins. Contrary to ligands like CO and O-2, which only bind to ferrous (Fe-II) iron, NO binds to both ferrous and ferric (Fe-II) iron. In a particular protein, the natural oxidation state can therefore be expected to be tailored to the required function. Herein, we present an ob initio potential-energy surface for ferric iron interacting with NO. This potential-energy surface exhibits three minima corresponding to eta'-NO coordination (the global minimum), eta(1)-ON coordination and eta(2) coordination. This contrasts with the potential-energy surface for Fe-II-NO, which ex- hibits only two minima (the eta(2) coordination mode for Fe-II is a transition state, not a minimum). In addition, the binding energies of NO are substantially larger for Fe-III than for Fe-II. We have performed molecular dynamics simulations for NO bound to ferric myoglobin (Mb(III)) and compare these with results obtained for Mb(II). Over the duration of our simulations (1.5 ns), all three binding modes are found to be stable at 200 K and transiently stable at 300 K, with eventual transformation to the eta(1)-NO global-minimum conformation. We discuss the implication of these results related to studies of rebinding processes in myoglobin.

Gas-phase kinetics of chlorosilylene reactions II. ClSiH + C2H4: absolute rate measurements and quantum chemical and RRKM calculations for the prototype pi addition reaction

Time-resolved studies of chlorosilylene, ClSiH, generated by the 193 nm laser flash photolysis of 1-chloro-1-silacyclopent-3-ene, are carried out to obtain rate constants for its bimolecular reaction with ethene, C2H4, in the gas-phase. The reaction is studied over the pressure range 0.13-13.3 kPa (with added SF6) at five temperatures in the range 296-562 K. The second order rate constants, obtained by extrapolation to the high pressure limits at each temperature, fitted the Arrhenius equation: log(k(infinity)/cm(3) molecule(-1) s(-1))=(-10.55 +/- 0.10) + (3.86 +/- 0.70) kJ mol(-1)/RT ln10. The Arrhenius parameters correspond to a loose transition state and the rate constant at room temperature is 43% of that for SiH2 + C2H4, showing that the deactivating effect of Cl-for-H substitution in the silylene is not large. Quantum chemical calculations of the potential energy surface for this reaction at the G3MP2//B3LYP level show that, as well as 1-chlorosilirane, ethylchlorosilylene is a viable product. The calculations reveal how the added effect of the Cl atom on the divalent state stabilisation of ClSiH influences the course of this reaction. RRKM calculations of the reaction pressure dependence suggest that ethylchlorosilylene should be the main product. The results are compared and contrasted with those of SiH2 and SiCl2 with C2H4.

Spectro-electrochemical and DFT studies of a planar Cu(II)–phenolate complex active in the aerobic oxidation of primary alcohols

A square-planar compound [Cu(pyrimol)Cl] (pyrimol = 4-methyl-2-N-(2-pyridylmethylene)aminophenolate) abbreviated as CuL–Cl) is described as a biomimetic model of the enzyme galactose oxidase (GOase). This copper(II) compound is capable of stoichiometric aerobic oxidation of activated primary alcohols in acetonitrile/water to the corresponding aldehydes. It can be obtained either from Hpyrimol (HL) or its reduced/hydrogenated form Hpyramol (4-methyl-2-N-(2-pyridylmethyl)aminophenol; H2L) readily converting to pyrimol (L-) on coordination to the copper(II) ion. Crystalline CuL–Cl and its bromide derivative exhibit a perfect square-planar geometry with Cu–O(phenolate) bond lengths of 1.944(2) and 1.938(2) Å. The cyclic voltammogram of CuL–Cl exhibits an irreversible anodic wave at +0.50 and +0.57 V versus ferrocene/ferrocenium (Fc/Fc+) in dry dichloromethane and acetonitrile, respectively, corresponding to oxidation of the phenolate ligand to the corresponding phenoxyl radical. In the strongly donating acetonitrile the oxidation path involves reversible solvent coordination at the Cu(II) centre. The presence of the dominant CuII–L. chromophore in the electrochemically and chemically oxidised species is evident from a new fairly intense electronic absorption at 400–480 nm ascribed to a several electronic transitions having a mixed pi-pi(L.) intraligand and Cu–Cl -> L. charge transfer character. The EPR signal of CuL–Cl disappears on oxidation due to strong intramolecular antiferromagnetic exchange coupling between the phenoxyl radical ligand (L.) and the copper(II) centre, giving rise to a singlet ground state (S = 0). The key step in the mechanism of the primary alcohol oxidation by CuL–Cl is probably the alpha-hydrogen abstraction from the equatorially bound alcoholate by the phenoxyl moiety in the oxidised pyrimol ligand, Cu–L., through a five-membered cyclic transition state.