Molecular dynamics simulation of liquid trimethylphosphine

Structural and dynamical properties of liquid trimethylphosphine (TMP), (CH(3))(3)P, as a function of temperature is investigated by molecular dynamics (MD) simulations. The force field used in the MD simulations, which has been proposed from molecular mechanics and quantum chemistry calculations, is able to reproduce the experimental density of liquid TMP at room temperature. Equilibrium structure is investigated by the usual radial distribution function, g(r), and also in the reciprocal space by the static structure factor, S(k). On the basis of center of mass distances, liquid TMP behaves like a simple liquid of almost spherical particles, but orientational correlation due to dipole-dipole interactions is revealed at short-range distances. Single particle and collective dynamics are investigated by several time correlation functions. At high temperatures, diffusion and reorientation occur at the same time range as relaxation of the liquid structure. Decoupling of these dynamic properties starts below ca. 220 K, when rattling dynamics of a given TMP molecules due to the cage effect of neighbouring molecules becomes important. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3624408]

Cyclooligomerisation of azido-alkyne-functionalised sugars: synthesis of 1,6-linked cyclic pseudo-galactooligosaccharides and assessment of their sialylation by Trypanosoma cruzi trans-sialidase

Cyclic pseudo-galactooligosaccharides were synthesized by cyclooligomerisation of isomeric azido-alkyne derivatives of beta-D-galactopyranose under Cu(I)-catalysed azide-alkyne 1,3-dipolar cycloaddition reaction conditions. The principal products isolated were cyclic dimers and trimers, with lower amounts of cyclic tetramer and pentamer also evident in some cases. Molecular mechanics calculations suggest very compact but flexible structures for the cyclic trimers, with secondary OH groups exposed outside the macrocycle and available for enzymatic glycosylation. The cyclic dimers and trimers represent a new type of acceptor substrate for Trypanosoma cruzi trans-sialidase, giving rise to doubly and triply sialylated glycomacrocycles, respectively.

Detailed H-1 and C-13 NMR structural assignment and relative stereochemistry determination for three closely related sesquiterpene lactones

A complete analysis of H-1 and C-13 NMR spectra of the trypanocidal sesquiterpene lactone eremantholide C and two of its analogues is described. These structurally similar sesquiterpene lactones were submitted to H-1 NMR, C-13 (H-1) NMR, gCOSY, gHSQC, gHMBC, J-resolved and DPFGSE-NOE NMR techniques. The detailed analysis of those results, correlated to some computational calculations (molecular mechanics), led to the total and unequivocal assignment of all H-1 and C-13 NMR data. The determination of all H-1/H-1 coupling constants and all signal multiplicities, together with the elimination of previous ambiguities were also achieved. Copyright (C) 2008 John Wiley & Sons, Ltd.

Theoretical Analysis of Aggregation in Block-Copolymer Films: The Optical Signature

We focus this work on the theoretical investigation of the block-copolymer poly [oxyoctyleneoxy-(2,6-dimethoxy-1,4phenylene-1,2-ethinylene-phenanthrene-2,4diyl) named as LaPPS19, recently proposed for optoelectronic applications. We used for that a variety of methods, from molecular mechanics to quantum semiempirical techniques (AMI, ZINDO/S-CIS). Our results show that as expected isolated LaPPS19 chains present relevant electron localization over the phenanthrene group. We found, however, that LaPPS19 could assemble in a pi-stacked form, leading to impressive interchain interaction; the stacking induces electronic delocalization between neighbor chains and introduces new states below the phenanthrene-related absorption; these results allowed us to associate the red-shift of the absorption edge, seen in the experimental results, to spontaneous pi-stack aggregation of the chains. (C) 2009 Wiley Periodicals, Inc. Int J Quantum Chem 110: 885-892, 2010

Electronic properties of liquid ammonia: A sequential molecular dynamics/quantum mechanics approach

The electronic properties of liquid ammonia are investigated by a sequential molecular dynamics/quantum mechanics approach. Quantum mechanics calculations for the liquid phase are based on a reparametrized hybrid exchange-correlation functional that reproduces the electronic properties of ammonia clusters [(NH(3))(n); n=1-5]. For these small clusters, electron binding energies based on Green's function or electron propagator theory, coupled cluster with single, double, and perturbative triple excitations, and density functional theory (DFT) are compared. Reparametrized DFT results for the dipole moment, electron binding energies, and electronic density of states of liquid ammonia are reported. The calculated average dipole moment of liquid ammonia (2.05 +/- 0.09 D) corresponds to an increase of 27% compared to the gas phase value and it is 0.23 D above a prediction based on a polarizable model of liquid ammonia [Deng , J. Chem. Phys. 100, 7590 (1994)]. Our estimate for the ionization potential of liquid ammonia is 9.74 +/- 0.73 eV, which is approximately 1.0 eV below the gas phase value for the isolated molecule. The theoretical vertical electron affinity of liquid ammonia is predicted as 0.16 +/- 0.22 eV, in good agreement with the experimental result for the location of the bottom of the conduction band (-V(0)=0.2 eV). Vertical ionization potentials and electron affinities correlate with the total dipole moment of ammonia aggregates. (c) 2008 American Institute of Physics.

Combined Monte Carlo and quantum mechanics study of the solvatochromism of phenol in water. The origin of the blue shift of the lowest pi-pi* transition

A combined and sequential use of Monte Carlo simulations and quantum mechanical calculations is made to analyze the spectral shift of the lowest pi-pi* transition of phenol in water. The solute polarization is included using electrostatic embedded calculations at the MP2/aug-cc-pVDZ level giving a dipole moment of 2.25 D, corresponding to an increase of 76% compared to the calculated gas-phase value. Using statistically uncorrelated configurations sampled from the MC simulation,first-principle size-extensive calculations are performed to obtain the solvatochromic shift. Analysis is then made of the origin of the blue shift. Results both at the optimized geometry and in room-temperature liquid water show that hydrogen bonds of water with phenol promote a red shift when phenol is the proton-donor and a blue shift when phenol is the proton-acceptor. In the case of the optimized clusters the calculated shifts are in very good agreement with results obtained from mass-selected free jet expansion experiments. In the liquid case the contribution of the solute-solvent hydrogen bonds partially cancels and the total shift obtained is dominated by the contribution of the outer solvent water molecules. Our best result, including both inner and outer water molecules, is 570 +/- 35 cm(-1), in very good agreement with the small experimental shift of 460 cm(-1) for the absorption maximum.

The isotropic nuclear magnetic shielding constants of acetone in supercritical water: A sequential Monte Carlo/quantum mechanics study including solute polarization

The nuclear isotropic shielding constants sigma((17)O) and sigma((13)C) of the carbonyl bond of acetone in water at supercritical (P=340.2 atm and T=673 K) and normal water conditions have been studied theoretically using Monte Carlo simulation and quantum mechanics calculations based on the B3LYP/6-311++G(2d,2p) method. Statistically uncorrelated configurations have been obtained from Monte Carlo simulations with unpolarized and in-solution polarized solute. The results show that solvent effects on the shielding constants have a significant contribution of the electrostatic interactions and that quantitative estimates for solvent shifts of shielding constants can be obtained modeling the water molecules by point charges (electrostatic embedding). In supercritical water, there is a decrease in the magnitude of sigma((13)C) but a sizable increase in the magnitude of sigma((17)O) when compared with the results obtained in normal water. It is found that the influence of the solute polarization is mild in the supercritical regime but it is particularly important for sigma((17)O) in normal water and its shielding effect reflects the increase in the average number of hydrogen bonds between acetone and water. Changing the solvent environment from normal to supercritical water condition, the B3LYP/6-311++G(2d,2p) calculations on the statistically uncorrelated configurations sampled from the Monte Carlo simulation give a (13)C chemical shift of 11.7 +/- 0.6 ppm for polarized acetone in good agreement with the experimentally inferred result of 9-11 ppm. (C) 2008 American Institute of Physics.

Effects of chronic L-NAME treatment lung tissue mechanics, eosinophilic and extracellular matrix responses induced by chronic pulmonary inflammation

The importance of lung tissue in asthma pathophysiology has been recently recognized. Although nitric oxide mediates smooth muscle tonus control in airways, its effects on lung tissue responsiveness have not been investigated previously. We hypothesized that chronic nitric oxide synthase (NOS) inhibition by N-omega-nitro-L-arginine methyl ester (L-NAME) may modulate lung tissue mechanics and eosinophil and extracellular matrix remodeling in guinea pigs with chronic pulmonary inflammation. Animals were submitted to seven saline or ovalbumin exposures with increasing doses (1 similar to 5 mg/ml for 4 wk) and treated or not with L-NAME in drinking water. After the seventh inhalation (72 h), animals were anesthetized and exsanguinated, and oscillatory mechanics of lung tissue strips were performed in baseline condition and after ovalbumin challenge (0.1%). Using morphometry, we assessed the density of eosinophils, neuronal NOS (nNOS)- and inducible NOS (iNOS)-positive distal lung cells, smooth muscle cells, as well as collagen and elastic fibers in lung tissue. Ovalbumin-exposed animals had an increase in baseline and maximal tissue resistance and elastance, eosinophil density, nNOS- and iNOS-positive cells, the amount of collagen and elastic fibers, and isoprostane-8-PGF(2 alpha) expression in the alveolar septa compared with controls (P < 0.05). L-NAME treatment in ovalbumin-exposed animals attenuated lung tissue mechanical responses (P < 0.01), nNOS- and iNOS-positive cells, elastic fiber content (P < 0.001), and isoprostane-8-PGF(2 alpha) in the alveolar septa (P < 0.001). However, this treatment did not affect the total number of eosinophils and collagen deposition. These data suggest that NO contributes to distal lung parenchyma constriction and to elastic fiber deposition in this model. One possibility may be related to the effects of NO activating the oxidative stress pathway.

A Hopping Mechanism for Cargo Transport by Molecular Motors on Crowded Microtubules

Most models designed to study the bidirectional movement of cargos as they are driven by molecular motors rely on the idea that motors of different polarities can be coordinated by external agents if arranged into a motor-cargo complex to perform the necessary work Gross, Hither and yon: a review of bidirectional microtubule-based transport (Gross in Phys. Biol. 1:R1-R11, 2004). Although these models have provided us with important insights into these phenomena, there are still many unanswered questions regarding the mechanisms through which the movement of the complex takes place on crowded microtubules. For example (i) how does cargo-binding affect motor motility? and in connection with that-(ii) how does the presence of other motors (and also other cargos) on the microtubule affect the motility of the motor-cargo complex? We discuss these questions from a different perspective. The movement of a cargo is conceived here as a hopping process resulting from the transference of cargo between neighboring motors. In the light of this, we examine the conditions under which cargo might display bidirectional movement even if directed by motors of a single polarity. The global properties of the model in the long-time regime are obtained by mapping the dynamics of the collection of interacting motors and cargos into an asymmetric simple exclusion process (ASEP) which can be resolved using the matrix ansatz introduced by Derrida (Derrida and Evans in Nonequilibrium Statistical Mechanics in One Dimension, pp. 277-304, 1997; Derrida et al. in J. Phys. A 26: 1493-1517, 1993).

Electronic properties of liquid hydrogen fluoride: A sequential quantum mechanical/Born-Oppenheimer molecular dynamics approach

The electronic properties of liquid hydrogen fluoride (HF) were investigated by carrying out sequential quantum mechanics/Born-Oppenheimer molecular dynamics. The structure of the liquid is in good agreement with recent experimental information. Emphasis was placed on the analysis of polarisation effects, dynamic polarisability and electronic excitations in liquid HF. Our results indicate an increase in liquid phase of the dipole moment (similar to 0.5 D) and isotropic polarisability (5%) relative to their gas-phase values. Our best estimate for the first vertical excitation energy in liquid HF indicates a blue-shift of 0.4 +/- 0.2 eV relative to that of the gas-phase monomer (10.4 eV). (C) 2010 Elsevier B.V. All rights reserved.

The dynamics of cargo driven by molecular motors in the context of asymmetric simple exclusion processes

We consider the dynamics of cargo driven by a collection of interacting molecular motors in the context of ail asymmetric simple exclusion process (ASEP). The model is formulated to account for (i) excluded-volume interactions, (ii) the observed asymmetry of the stochastic movement of individual motors and (iii) interactions between motors and cargo. Items (i) and (ii) form the basis of ASEP models and have already been considered to study the behavior of motor density profile [A. Parmeggiani. T. Franosch, E. Frey, Phase Coexistence in driven one-dimensional transport, Phys. Rev. Lett. 90 (2003) 086601-1-086601-4]. Item (iii) is new. It is introduced here as an attempt to describe explicitly the dependence of cargo movement on the dynamics of motors in this context. The steady-state Solutions Of the model indicate that the system undergoes a phase transition of condensation type as the motor density varies. We study the consequences of this transition to the behavior of the average cargo velocity. (C) 2009 Elsevier B.V. All rights reserved.

Cellular automaton model for molecular traffic jams

We consider the time evolution of an exactly solvable cellular automaton with random initial conditions both in the large-scale hydrodynamic limit and on the microscopic level. This model is a version of the totally asymmetric simple exclusion process with sublattice parallel update and thus may serve as a model for studying traffic jams in systems of self-driven particles. We study the emergence of shocks from the microscopic dynamics of the model. In particular, we introduce shock measures whose time evolution we can compute explicitly, both in the thermodynamic limit and for open boundaries where a boundary-induced phase transition driven by the motion of a shock occurs. The motion of the shock, which results from the collective dynamics of the exclusion particles, is a random walk with an internal degree of freedom that determines the jump direction. This type of hopping dynamics is reminiscent of some transport phenomena in biological systems.

Molecular systematics of Thorea (Rhodophyta, Thoreales) species in Brazil

This study aimed to evaluate species level taxonomy and phylogenetic relationship among Thorea species in Brazil and other regions of the world using two molecular markers - RUBISCO large subunit plastid gene (rbcL) and nuclear small-subunit ribosomal DNA (SSU rDNA). Three samples of Thorea from Brazil (states of Mato Grosso do Sul and São Paulo) and one sample from Dominican Republic (DR) were sequenced. Analyses based on partial sequences of rbcL (1,282 bp) and complete sequences of SSU (1,752 bp) were essentially congruent and revealed that Thoreales formed a distinct monophyletic clade, which had two major branches with high support, representing the genera Thorea and Nemalionopsis. Thorea clade had four main branches with high support for all analyses, each one representing the species: 1) T. gaudichaudii C. Agardh from Asia (Japan and Philippines) - this clade occurred only in the rbcL analyses; 2) T. violacea Bory from Asia (Japan) and North America (U.S.A. and DR); 3) T. hispida (Thore) Desvaux from Europe (England) and Asia (Japan); 4) a distinct group with the three Brazilian samples (sequence identity: rbcL 97.2%, 1,246 bp; SSU 96.0-98.1%, 1,699-1,720 bp). The Brazilian samples clearly formed a monophyletic clade based on both molecular markers and was interpreted as a separate species, for which we resurrected the name T. bachmannii Pujals. Morphological and molecular evidences indicate that the Thoreales is well-resolved at ordinal and generic levels. In contrast, Thorea species recognized by molecular data require additional characters (e.g. reproductive and chromosome numbers) to allow consistent and reliable taxonomic circumscription aiming at a world revision based on molecular and morphological evidences.

Thaptomys Thomas 1915 (Rodentia, Sigmodontinae, Akodontini) with karyotypes 2n = 50, FN = 48, and 2n = 52, FN = 52: two monophyletic lineages recovered by molecular phylogeny

A novel karyotype with 2n = 50, FN = 48, was described for specimens of Thaptomys collected at Una, State of Bahia, Brazil, which are morphologically indistinguishable from Thaptomys nigrita, 2n = 52, FN = 52, found in other localities. It was hence proposed that the 2n = 50 karyotype could belong to a distinct species, cryptic of Thaptomys nigrita, once chromosomal rearrangements observed, along with the geographic distance, might represent a reproductive barrier between both forms. Phylogenetic analyses using maximum parsimony and maximum likelihood based on partial cytochrome b sequences with 1077 bp were performed, attempting to establish the relationships among the individuals with distinct karyotypes along the geographic distribution of the genus; the sample comprised 18 karyotyped specimens of Thaptomys, encompassing 15 haplotypes, from eight different localities of the Atlantic Rainforest. The intra-generic relationships corroborated the distinct diploid numbers, once both phylogenetic reconstructions recovered two monophyletic lineages, a northeastern clade grouping the 2n = 50 and a southeastern clade with three subclades, grouping the 2n = 52 karyotype. The sequence divergence observed between their individuals ranged from 1.9% to 3.5%.