Solvent Effects in Chemical Processes. Water-Assisted Proton Transfer Reaction of Pterin in Aqueous Environment

Pterins are members of a family of heterocyclic compounds present in a wide variety of biological systems and may exist in two forms, corresponding to an acid and a basic tautomer. In this work, the proton transfer reaction between these tautomeric forms was investigated in the gas phase and in aqueous solution. In gas phase, the intramolecular mechanism was carried out for die isolated pterin by quantum mechanical second-order Moller-Plesset Perturbation theory (MP2/aug-cc-pVDZ) calculations and it indicates that the acid form is more stable than the basic form by -1.4 kcal/mol with a barrier of 34.2 kcal/mol with respect to the basic form. In aqueous solution, the role of the water molecules in the proton transfer reaction was analyzed in two separated parts, the direct participation of one water molecule in the reaction path, called water-assisted mechanism, and the complementary participation of the aqueous solvation. The water-assisted mechanism was carried out for one pterin-water cluster by quantum mechanical calculations and it indicates that the acid form is still more stable by -3.3 kcal/mol with a drastic reduction of 70% of the barrier, The bulk solution effect on the intramolecular and water-assisted mechanisms was included by free energy perturbation implemented on Monte Carlo simulations. The bulk water effect is found to be substantial and decisive when the reaction path involves the water-assisted mechanism. In this case, the free energy barrier is only 6.7 kcal/mol and the calculated relative Gibbs free energy for the two tautomers is -11.2 kcal/mol. This value is used to calculate the pK(a) value of 8.2 +/- 0.6 that is in excellent agreement with the experimental result of 7.9.

Dipole polarizability and Rayleigh light scattering by the hydrated electron

Assuming the existence of a confined state of the electron in bulk water the polarizability of the hydrated electron is analyzed. Statistically uncorrelated supermolecular structures composed of seven water molecules (first solvation shell) with an extra electron were extracted from classical Monte Carlo simulation and used in quantum mechanical second-order Moller-Plesset calculations. It is found that the bound excess electron contributes with 274 a.u. to the total dipole polarizability of 345 a.u. for (H(2)O)(7)(-). From the calculated polarizabilities the Rayleigh elastic light scattering properties are inferred and found to considerably enhance activity and light depolarization. (C) 2009 Elsevier B.V. All rights reserved.

Compressible Sherrington-Kirkpatrick spin-glass model

We introduce a Sherrington-Kirkpatrick spin-glass model with the addition of elastic degrees of freedom. The problem is formulated in terms of an effective four-spin Hamiltonian in the pressure ensemble, which can be treated by the replica method. In the replica-symmetric approximation, we analyze the pressure-temperature phase diagram, and obtain expressions for the critical boundaries between the disordered and the ordered (spin-glass and ferromagnetic) phases. The second-order para-ferromagnetic border ends at a tricritical point, beyond which the transition becomes discontinuous. We use these results to make contact with the temperature-concentration phase diagrams of mixtures of hydrogen-bonded crystals.

Antiferromagnetic spherical spin-glass model

We study the thermodynamic properties and the phase diagrams of a multi-spin antiferromagnetic spherical spin-glass model using the replica method. It is a two-sublattice version of the ferromagnetic spherical p-spin glass model. We consider both the replica-symmetric and the one-step replica-symmetry-breaking solutions, the latter being the most general solution for this model. We find paramagnetic, spin-glass, antiferromagnetic and mixed or glassy antiferromagnetic phases. The phase transitions are always of second order in the thermodynamic sense, but the spin-glass order parameter may undergo a discontinuous change.

Photo-irradiation effects on the surface morphology of poly(p-phenylene vinylene) films

In this work, we have studied the surface morphology of photo-irradiated poly(p-phenylene vinylene) (PPV) thin films by using atomic force microscopy (AFM). We have analyzed the first-order statistical parameters, the height distribution and the distance between selected peaks. The second-order statistical analysis was introduced calculating the auto-covariance function to determine the correlation length between heights. We have observed that the photo-irradiation process produces a surface topology more homogeneous and isotropic such as a normal surface. In addition, the polymer surface irradiation can be used as a new methodology to obtain materials optically modified. (C) 2009 Elsevier B.V. All rights reserved.

Self-dual hopfions

We construct static and time-dependent exact soliton solutions with nontrivial Hopf topological charge for a field theory in 3 + 1 dimensions with the target space being the two dimensional sphere S(2). The model considered is a reduction of the so-called extended Skyrme-Faddeev theory by the removal of the quadratic term in derivatives of the fields. The solutions are constructed using an ansatz based on the conformal and target space symmetries. The solutions are said self-dual because they solve first order differential equations which together with some conditions on the coupling constants, imply the second order equations of motion. The solutions belong to a sub-sector of the theory with an infinite number of local conserved currents. The equation for the profile function of the ansatz corresponds to the Bogomolny equation for the sine-Gordon model.

WRITING ELECTRONIC FERROMAGNETIC STATES IN A HIGH-TEMPERATURE PARAMAGNETIC NUCLEAR SPIN SYSTEM

In this paper, we use Nuclear Magnetic Resonance (NMR) to write electronic states of a ferromagnetic system into high-temperature paramagnetic nuclear spins. Through the control of phase and duration of radio frequency pulses, we set the NMR density matrix populations, and apply the technique of quantum state tomography to experimentally obtain the matrix elements of the system, from which we calculate the temperature dependence of magnetization for different magnetic fields. The effects of the variation of temperature and magnetic field over the populations can be mapped in the angles of spin rotations, carried out by the RF pulses. The experimental results are compared to the Brillouin functions of ferromagnetic ordered systems in the mean field approximation for two cases: the mean field is given by (i) B = B(0) + lambda M and (ii) B = B(0) + lambda M + lambda`M(3), where B(0) is the external magnetic field, and lambda, lambda` are mean field parameters. The first case exhibits second order transition, whereas the second case has first order transition with temperature hysteresis. The NMR simulations are in good agreement with the magnetic predictions.

A quasispecies approach to the evolution of sexual replication in unicellular organisms

This study develops a simplified model describing the evolutionary dynamics of a population composed of obligate sexually and asexually reproducing, unicellular organisms. The model assumes that the organisms have diploid genomes consisting of two chromosomes, and that the sexual organisms replicate by first dividing into haploid intermediates, which then combine with other haploids, followed by the normal mitotic division of the resulting diploid into two new daughter cells. We assume that the fitness landscape of the diploids is analogous to the single-fitness-peak approach often used in single-chromosome studies. That is, we assume a master chromosome that becomes defective with just one point mutation. The diploid fitness then depends on whether the genome has zero, one, or two copies of the master chromosome. We also assume that only pairs of haploids with a master chromosome are capable of combining so as to produce sexual diploid cells, and that this process is described by second-order kinetics. We find that, in a range of intermediate values of the replication fidelity, sexually reproducing cells can outcompete asexual ones, provided the initial abundance of sexual cells is above some threshold value. The range of values where sexual reproduction outcompetes asexual reproduction increases with decreasing replication rate and increasing population density. We critically evaluate a common approach, based on a group selection perspective, used to study the competition between populations and show its flaws in addressing the evolution of sex problem.

Genericity of Nondegeneracy for Light Rays in Stationary Spacetimes

Given a Lorentzian manifold (M, g), an event p and an observer U in M, then p and U are light conjugate if there exists a lightlike geodesic gamma : [0, 1] -> M joining p and U whose endpoints are conjugate along gamma. Using functional analytical techniques, we prove that if one fixes p and U in a differentiable manifold M, then the set of stationary Lorentzian metrics in M for which p and U are not light conjugate is generic in a strong sense. The result is obtained by reduction to a Finsler geodesic problem via a second order Fermat principle for light rays, and using a transversality argument in an infinite dimensional Banach manifold setup.

Minimal Lagrangian surfaces in the tangent bundle of a Riemannian surface

Given an oriented Riemannian surface (Sigma, g), its tangent bundle T Sigma enjoys a natural pseudo-Kahler structure, that is the combination of a complex structure 2, a pseudo-metric G with neutral signature and a symplectic structure Omega. We give a local classification of those surfaces of T Sigma which are both Lagrangian with respect to Omega and minimal with respect to G. We first show that if g is non-flat, the only such surfaces are affine normal bundles over geodesics. In the flat case there is, in contrast, a large set of Lagrangian minimal surfaces, which is described explicitly. As an application, we show that motions of surfaces in R(3) or R(1)(3) induce Hamiltonian motions of their normal congruences, which are Lagrangian surfaces in TS(2) or TH(2) respectively. We relate the area of the congruence to a second-order functional F = f root H(2) - K dA on the original surface. (C) 2010 Elsevier B.V. All rights reserved.

A matrix formula for the skewness of maximum likelihood estimators

We give a general matrix formula for computing the second-order skewness of maximum likelihood estimators. The formula was firstly presented in a tensorial version by Bowman and Shenton (1998). Our matrix formulation has numerical advantages, since it requires only simple operations on matrices and vectors. We apply the second-order skewness formula to a normal model with a generalized parametrization and to an ARMA model. (c) 2010 Elsevier B.V. All rights reserved.

Improved maximum-likelihood estimation in a regression model with general parametrization

We analyse the finite-sample behaviour of two second-order bias-corrected alternatives to the maximum-likelihood estimator of the parameters in a multivariate normal regression model with general parametrization proposed by Patriota and Lemonte [A. G. Patriota and A. J. Lemonte, Bias correction in a multivariate regression model with genereal parameterization, Stat. Prob. Lett. 79 (2009), pp. 1655-1662]. The two finite-sample corrections we consider are the conventional second-order bias-corrected estimator and the bootstrap bias correction. We present the numerical results comparing the performance of these estimators. Our results reveal that analytical bias correction outperforms numerical bias corrections obtained from bootstrapping schemes.

Bias correction in a multivariate normal regression model with general parameterization

This paper derives the second-order biases Of maximum likelihood estimates from a multivariate normal model where the mean vector and the covariance matrix have parameters in common. We show that the second order bias can always be obtained by means of ordinary weighted least-squares regressions. We conduct simulation studies which indicate that the bias correction scheme yields nearly unbiased estimators. (C) 2009 Elsevier B.V. All rights reserved.

INSECT CHYMOTRYPSINS: CHLOROMETHYL KETONE INACTIVATION AND SUBSTRATE SPECIFICITY RELATIVE TO POSSIBLE COEVOLUTIONAL ADAPTATION OF INSECTS AND PLANTS

Insect digestive chymotrypsins are present in a large variety of insect orders but their substrate specificity still remains unclear. Ewer insect chymotrypsins from 3 different insect orders (Dictyoptera, Coleoptera and two Lepidoptera) were isolated using affinity chromatography. Enzymes presented molecular masses in the range of 20 to 31 kDa and pH optima in the range of 7.5 to 10.0. Kinetic characterization. using different, colorimetric and fluorescent substrates indicated that insect chymotrypsins differ from, bovine chymotrypsin in their primary specificity toward small substrates (like N-benzoyl-L-Tyr p-nitroanilide) rather than on their preference for large substrates (exemplified by Succynil-Ala-Ala-Pro-Phe P-nitroanilide). Chloromethyl ketones (TPCK, N-alpha-tosyl-L-Phe chloromethyl ketone and Z-GGF-CK, N-carbobenzoxy-Gly-Gly-phe-CK) inactivated all chymotrypsins legated. Inactivation rates follow apparent first-order kinetics with variable second order rates (TPCK, 42 to 130 M(-1)s(-1); Z-GGF-CK, 150 to 450 M(-1)s(-1) that may be remarkably low for S. frugiperda chymotrypsin (TPCK, 6 M(-1)s(-1); Z-GGF-CK, 6.1 M(-1) s(-1)). Homology modelling and sequence alignment showed that. in lepidopteran chymotrypsins, differences in the amino acid residues in the neighborhood of the catalytic His 57 may affect its pKa, value. This is Proposed as the cause of the decrease in His 57 reactivity toward chloromethyl ketones. Such amino acid replacement in the active site is proposed. to be an adaptation to the presence of dietary ketones. (C) 2009 Wiley Periodicals, Inc.

Horseradish peroxidase compound I as a tool to investigate reactive protein-cysteine residues: from quantification to kinetics

Proteins containing reactive cysteine residues (protein-Cys) are receiving increased attention as mediators of hydrogen peroxide signaling. These proteins are mainly identified by mining the thiol proteomes of oxidized protein-Cys in cells and tissues. However, it is difficult to determine if oxidation occurs through a direct reaction with hydrogen peroxide or by thiol-disulfide exchange reactions. Kinetic studies with purified proteins provide invaluable information about the reactivity of protein-Cys residues with hydrogen peroxide. Previously, we showed that the characteristic UV-Vis spectrum of horseradish peroxidase compound I, produced from the oxidation of horseradish peroxidase by hydrogen peroxide, is a simple, reliable, and useful tool to determine the second-order rate constant of the reaction of reactive protein-Cys with hydrogen peroxide and peroxynitrite. Here, the method is fully described and extended to quantify reactive protein-Cys residues and micromolar concentrations of hydrogen peroxide. Members of the peroxiredoxin family were selected for the demonstration and validation of this methodology. In particular, we determined the pK(a) of the peroxidatic thiol of rPrx6 (5.2) and the second-order rate constant of its reactions with hydrogen peroxide ((3.4 +/- 0.2) x 10(7) M(-1) s(-1)) and peroxynitrite ((3.7 +/- 0.4) x 10(5) M(-1) s(-1)) at pH 7.4 and 25 degrees C. (C) 2011 Elsevier Inc. All rights reserved.