Accurate ab initio potential energy surfaces for the (3)A '' and (3)A ' electronic states of the O(P-3) plus HBr system

In this work, we report the construction of potential energy surfaces for the (3)A '' and (3)A' states of the system O(P-3) + HBr. These surfaces are based on extensive ab initio calculations employing the MRCI+Q/CBS+SO level of theory. The complete basis set energies were estimated from extrapolation of MRCI+Q/aug-cc-VnZ(-PP) (n = Q, 5) results and corrections due to spin-orbit effects obtained at the CASSCF/aug-cc-pVTZ(-PP) level of theory. These energies, calculated over a region of the configuration space relevant to the study of the reaction O(P-3) + HBr -> OH + Br, were used to generate functions based on the many-body expansion. The three-body potentials were interpolated using the reproducing kernel Hilbert space method. The resulting surface for the (3)A '' electronic state contains van der Waals minima on the entrance and exit channels and a transition state 6.55 kcal/mol higher than the reactants. This barrier height was then scaled to reproduce the value of 5.01 kcal/mol, which was estimated from coupled cluster benchmark calculations performed to include high-order and core-valence correlation, as well as scalar relativistic effects. The (3)A' surface was also scaled, based on the fact that in the collinear saddle point geometry these two electronic states are degenerate. The vibrationally adiabatic barrier heights are 3.44 kcal/mol for the (3)A '' and 4.16 kcal/mol for the (3)A' state. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4705428]

Capillary bioreactors based on human purine nucleoside phosphorylase: A new approach for ligands identification and characterization

The enzyme purine nucleoside phosphorylase (PNP) is a target for the discovery of new lead compounds employed on the treatment severe T-cell mediated disorders. Within this context, the development of new, direct, and reliable methods for ligands screening is an important task. This paper describes the preparation of fused silica capillaries human PNP (HsPNP) immobilized enzyme reactor (IMER). The activity of the obtained IMER is monitored on line in a multidimensional liquid chromatography system, by the quantification of the product formed throughout the enzymatic reaction. The Km value for the immobilized enzyme was about twofold higher than that measured for the enzyme in solution (255 +/- 29.2 mu M and 133 +/- 114.9 mu M, respectively). A new fourth-generation immucillin derivative (DI4G: IC50 = 40.6 +/- 0.36 nM), previously identified and characterized in HsPNP free enzyme assays, was used to validate the IMER as a screening method for HsPNP ligands. The validated method was also used for mechanistic studies with this inhibitor. This new approach is a valuable tool to PNP ligand screening, since it directly measures the hypoxanthine released by inosine phosphorolysis, thus furnishing more reliable results than those one used in a coupled enzymatic spectrophotometric assay. (C) 2011 Elsevier B.V. All rights reserved.

Insights into Phosphate Cooperativity and Influence of Substrate Modifications on Binding and Catalysis of Hexameric Purine Nucleoside Phosphorylases

The hexameric purine nucleoside phosphorylase from Bacillus subtilis (BsPNP233) displays great potential to produce nucleoside analogues in industry and can be exploited in the development of new anti-tumor gene therapies. In order to provide structural basis for enzyme and substrates rational optimization, aiming at those applications, the present work shows a thorough and detailed structural description of the binding mode of substrates and nucleoside analogues to the active site of the hexameric BsPNP233. Here we report the crystal structure of BsPNP233 in the apo form and in complex with 11 ligands, including clinically relevant compounds. The crystal structure of six ligands (adenine, 2'deoxyguanosine, aciclovir, ganciclovir, 8-bromoguanosine, 6-chloroguanosine) in complex with a hexameric PNP are presented for the first time. Our data showed that free bases adopt alternative conformations in the BsPNP233 active site and indicated that binding of the co-substrate (2'deoxy) ribose 1-phosphate might contribute for stabilizing the bases in a favorable orientation for catalysis. The BsPNP233-adenosine complex revealed that a hydrogen bond between the 5' hydroxyl group of adenosine and Arg(43*) side chain contributes for the ribosyl radical to adopt an unusual C3'-endo conformation. The structures with 6-chloroguanosine and 8-bromoguanosine pointed out that the Cl-6 and Br-8 substrate modifications seem to be detrimental for catalysis and can be explored in the design of inhibitors for hexameric PNPs from pathogens. Our data also corroborated the competitive inhibition mechanism of hexameric PNPs by tubercidin and suggested that the acyclic nucleoside ganciclovir is a better inhibitor for hexameric PNPs than aciclovir. Furthermore, comparative structural analyses indicated that the replacement of Ser(90) by a threonine in the B. cereus hexameric adenosine phosphorylase (Thr(91)) is responsible for the lack of negative cooperativity of phosphate binding in this enzyme.

Kinetics and product distribution of p-nitrophenyl phosphate dianion solvolysis in ternary DMSO/alcohol/water mixtures are compatible with metaphosphate formation

The rate of solvolysis of p-nitrophenyl phosphate (PNPP) dianion in DMSO/water strongly decreases by increasing water concentration. Addition of linear alcohols (methanol, propanol, butanol, pentanol, and hexanol) at constant DMSO/water molar ratio produced an even sharper rate decrease. Alkyl phosphate formation, resulting from PNPP solvolysis in ternary DMSO/water/alcohol mixtures, increased with alcohol concentration and was essentially temperature independent. Methanol and hexanol were the poorest nucleophiles under all conditions. Activation energies and enthalpies for solvolysis in ternary mixtures were similar and entropies varied with alcohol concentration. Taken together these results can be best interpreted in terms of a dissociative mechanism with the intervention of metaphosphate. Copyright (C) 2011 John Wiley & Sons, Ltd.

Electronic Properties of Water in Liquid Environment. A Sequential QM/MM Study Using the Free Energy Gradient Method

There is a continuous search for theoretical methods that are able to describe the effects of the liquid environment on molecular systems. Different methods emphasize different aspects, and the treatment of both the local and bulk properties is still a great challenge. In this work, the electronic properties of a water molecule in liquid environment is studied by performing a relaxation of the geometry and electronic distribution using the free energy gradient method. This is made using a series of steps in each of which we run a purely molecular mechanical (MM) Monte Carlo Metropolis simulation of liquid water and subsequently perform a quantum mechanical/molecular mechanical (QM/MM) calculation of the ensemble averages of the charge distribution, atomic forces, and second derivatives. The MP2/aug-cc-pV5Z level is used to describe the electronic properties of the QM water. B3LYP with specially designed basis functions are used for the magnetic properties. Very good agreement is found for the local properties of water, such as geometry, vibrational frequencies, dipole moment, dipole polarizability, chemical shift, and spin-spin coupling constants. The very good performance of the free energy method combined with a QM/MM approach along with the possible limitations are briefly discussed.

A simple analysis of the influence of the solvent-induced electronic polarization on the N-15 magnetic shielding of pyridine in water

Electronic polarization induced by the interaction of a reference molecule with a liquid environment is expected to affect the magnetic shielding constants. Understanding this effect using realistic theoretical models is important for proper use of nuclear magnetic resonance in molecular characterization. In this work, we consider the pyridine molecule in water as a model system to briefly investigate this aspect. Thus, Monte Carlo simulations and quantum mechanics calculations based on the B3LYP/6-311++G (d,p) are used to analyze different aspects of the solvent effects on the N-15 magnetic shielding constant of pyridine in water. This includes in special the geometry relaxation and the electronic polarization of the solute by the solvent. The polarization effect is found to be very important, but, as expected for pyridine, the geometry relaxation contribution is essentially negligible. Using an average electrostatic model of the solvent, the magnetic shielding constant is calculated as -58.7 ppm, in good agreement with the experimental value of -56.3 ppm. The explicit inclusion of hydrogen-bonded water molecules embedded in the electrostatic field of the remaining solvent molecules gives the value of -61.8 ppm.

Selectivity and Mechanisms Driven by Reaction Dynamics: The Case of the Gas-Phase OH-+CH3ONO2 Reaction

Well-established statistical approaches such as transition-state theory based on high-level calculated potential energy profiles are unable to account for the selectivity observed in the gas-phase OH- + CH3ONO2 reaction. This reaction can undergo bimolecular nucleophilic displacement at either the carbon center (S(N)2@C) or the nitrogen center (S(N)2@N) as well as a proton abstraction followed by dissociation (E(CO)2) pathway. Direct dynamics simulations yield an S(N)2:E(CO)2 product ratio in close agreement with experiment and show that the lack of reactivity at the nitrogen atom is due to the highly negative electrostatic potential generated by the oxygen atoms in the ONO2 group that scatters the incoming OH-. In addition to these dynamical effects, the nonstatistical behavior of these reactions is attributed to the absence of equilibrated reactant complexes and to the large number of recrossings, which might be present in several ion-molecule gas-phase reactions.

X-ray diffraction and theoretical investigation of the Gedunin crystal structure

The Gedunin compound (C28H34O6) is a natural product extracted from Trichilia pallida that has shown a wide activity. The crystallographic structure shows two conformers in the asymmetric unit, which differ in a rotation of the furan group. To understand this molecular arrangement, the density functional calculations. Molecular Electrostatic Potential (MEP) and thermodynamic function calculation have been performed at the B3LYP/6-311++g(d,p) level. Both conformers were optimized and the agreement with the experimental structure was very good, making possible further theoretical analysis of the structure. The inter-conversion between two conformers depends on the energy barrier. This process is studied in the vacuum and shows two transition states with a low energetic barrier for a potential energy curve scanning rigid around furan group: 4.37 kcal/mol and 16.52 kcal/mol. As the first transition state has a notably lower energetic barrier, the preferred inter-conversion pathway between the conformers involves the first rather than the second transition state. Understanding this transition state in detail led us to perform its optimization, showing an energetic barrier around 3.66 kcal/mol. The negative free energy and low enthalpy confirm that the process is spontaneous and exothermic. The results show that this requirement makes the existence of the two conformers in the asymmetric unit possible. The structure of molecules in the asymmetric unit is better understood when the MEP is used on the interaction between molecules. For Gedunin, both molecules have shown MEP with well-defined regions, and this behavior contributes to the observed link between molecules and for the negative regions complementing positive regions of another molecule. (C) 2011 Elsevier B.V. All rights reserved.

Incorporation of PbF2 into Heavy Metal Oxide Borate Glasses. Structural Studies by Solid State NMR

A series of heavy metal oxide (HMO) glasses with composition 26.66B(2)O(3)-16GeO(2)-4 Bi2O3-(53.33-x)PbO-xPbF2 (0 <= x <= 40) were prepared and characterized with respect to their bulk (glass transition and crystallization temperatures, densities, molar volumes) and spectroscopic properties. Homogeneous glasses are formed up to x = 30, while crystallization of beta-PbF2 takes place at higher contents. Substitution of PbO by PbF2 shifts the optical band gap toward higher energies, thereby extending the UV transmission window significantly toward higher frequencies. Raman and infrared absorption spectra can be interpreted in conjunction with published reference data. Using B-11 and F-19 high-resolution solid state NMR as well as B-11/F-19 double resonance methodologies, we develop a quantitative structural description of this material. The fraction of four-coordinate boron is found to be moderately higher compared to that in glasses with the same PbO/B2O3 ratios, suggesting some participation of PbF2 in the network transformation process. This suggestion is confirmed by the F-19 NMR spectra. While the majority of the fluoride ions is present as ionic fluoride, similar to 20% of the fluorine inventory acts as a network modifier, resulting in the formation of four-coordinate BO3/2F- units. These units can be identified by F-19{B-11} rotational echo double resonance and B-11{F-19} cross-polarization magic angle spinning (CPMAS) data. These results provide the first unambiguous evidence of B-F bonding in a PbF2-modified glass system. The majority of the fluoride ions are found in a lead-dominated environment. F-19-F-19 homonuclear dipolar second moments measured by spin echo decay spectroscopy are quantitatively consistent with a model in which these ions are randomly distributed within the network modifier subdomain consisting of PbO, Bi2O3, and PbF2. This model, which implies both the features of atomic scale mixing with the network former borate species and some degree of fluoride ion clustering is consistent with all of the experimental data obtained on these glasses.

Interaction of the meso-tetrakis (4-N-methylpyridyl) porphyrin with gel and liquid state phospholipid vesicles

The interaction of the cationic meso-tetrakis 4-N-methylpyridyl porphyrin (TMPyP) with large unilamellar vesicles (LUVs) was investigated in the present study. LUVs were formed by mixtures of the zwitterionic 1,2-dipalmitoyl-sn-glycero-phosphatidylcholine (DPPC) and anionic 1,2-dipalmitoyl-sn-glycero-3-phosphoglycerol (DPPG) phospholipids, at different DPPG molar percentages. All investigations were carried out above (50 degrees C) and below (25 degrees C) the main phase transition temperature of the LUVs (similar to 41 degrees C). The binding constant values, K-b, estimated from the time-resolved fluorescence study, showed a significant increase of the porphyrin affinity at higher mol% DPPG. This affinity is markedly increased when the LUVs are in the liquid crystalline state. For both situations, the increase of the K-b value was also followed by a higher porphyrin fraction bound to the LUVs. The displacement of the vesicle-bound porphyrins toward the aqueous medium, upon titration with the salt potassium chloride (KCl), was also studied. Altogether, our steady-state and frequency-domain fluorescence quenching data results indicate that the TMPyP is preferentially located at the LUVs Stern layer. This is supported by the zeta potential studies, where a partial neutralization of the LUVs surface charge, upon porphyrin titration, was observed. Dynamic light scattering (DLS) results showed that, for some phospholipid systems, this partial neutralization leads to the LUVs flocculation. (C) 2012 Elsevier Inc. All rights reserved.

Average control of Markov decision processes with Feller transition probabilities and general action spaces

This paper studies the average control problem of discrete-time Markov Decision Processes (MDPs for short) with general state space, Feller transition probabilities, and possibly non-compact control constraint sets A(x). Two hypotheses are considered: either the cost function c is strictly unbounded or the multifunctions A(r)(x) = {a is an element of A(x) : c(x, a) <= r} are upper-semicontinuous and compact-valued for each real r. For these two cases we provide new results for the existence of a solution to the average-cost optimality equality and inequality using the vanishing discount approach. We also study the convergence of the policy iteration approach under these conditions. It should be pointed out that we do not make any assumptions regarding the convergence and the continuity of the limit function generated by the sequence of relative difference of the alpha-discounted value functions and the Poisson equations as often encountered in the literature. (C) 2012 Elsevier Inc. All rights reserved.

Analysis of confinement potential fluctuation and band-gap renormalization effects on excitonic transition in GaAs/AlGaAs multiquantum wells grown on (100) and (311)A GaAs surfaces

The competition between confinement potential fluctuations and band-gap renormalization (BGR) in GaAs/AlxGa1-xAs quantum wells grown on [1 0 0] and [3 1 1]A GaAs substrates is evaluated. The results clearly demonstrate the coexistence of the band-tail states filling related to potential fluctuations and the band-gap renormalization caused by an increase in the density of photogenerated carriers during the photoluminescence (PL) experiments. Both phenomena have strong influence on temperature dependence of the PL-peak energy (E-PL(T)). As the photon density increases, the E-PL can shift to either higher or lower energies, depending on the sample temperature. The temperature at which the displacement changes from a blueshift to a redshift is governed by the magnitude of the potential fluctuations and by the variation of BGR with excitation density. A simple band-tail model with a Gaussian-like distribution of the density of state was used to describe the competition between the band-tail filling and the BGR effects on E-PL(T). (C) 2012 Elsevier B.V. All rights reserved.

Complex mixed state of the Pauli-limited superconductor CeCoIn5

Magnetization measurements were performed on CeCoIn5 at temperatures down to 20 mK and magnetic fields up to 17 T applied along different crystallographic orientations. For field configurations nearly parallel to the ab plane (theta less than or similar to 40 degrees and T <= 50 mK), we have found an intriguing vortex dynamics regime revealed by a hysteretic and metastable anomalous peak effect (APE), which gives evidence of surface barrier effects enhanced by antiferromagnetic fluctuations in the mixed state of CeCoIn5. Furthermore, we have observed crossover features in the torque and magnetization traces at fields below H-c2, which are consistent with vortices lattice phase transitions and with the anomalies speculated to be the Fulde-Ferrel-Larkin-Ovchinnikov (FFLO) superconducting state in CeCoIn5. All of the above features were found to be dramatically perturbed in Ce0.98Gd0.02CoIn5.

SYNCHRONIZATION OF CHAOS AND THE TRANSITION TO WAVE TURBULENCE

We investigated the transition to wave turbulence in a spatially extended three-wave interacting model, where a spatially homogeneous state undergoing chaotic dynamics undergoes spatial mode excitation. The transition to this weakly turbulent state can be regarded as the loss of synchronization of chaos of mode oscillators describing the spatial dynamics.

STRUCTURE AND DYNAMICS: THE TRANSITION FROM NONEQUILIBRIUM TO EQUILIBRIUM IN INTEGRATE-AND-FIRE DYNAMICS

The transient and equilibrium properties of dynamics unfolding in complex systems can depend critically on specific topological features of the underlying interconnections. In this work, we investigate such a relationship with respect to the integrate-and-fire dynamics emanating from a source node and an extended network model that allows control of the small-world feature as well as the length of the long-range connections. A systematic approach to investigate the local and global correlations between structural and dynamical features of the networks was adopted that involved extensive simulations (one and a half million cases) so as to obtain two-dimensional correlation maps. Smooth, but diverse surfaces of correlation values were obtained in all cases. Regarding the global cases, it has been verified that the onset avalanche time (but not its intensity) can be accurately predicted from the structural features within specific regions of the map (i.e. networks with specific structural properties). The analysis at local level revealed that the dynamical features before the avalanches can also be accurately predicted from structural features. This is not possible for the dynamical features after the avalanches take place. This is so because the overall topology of the network predominates over the local topology around the source at the stationary state.