Second-order Møller-Plesset perturbation theory without basis set superposition error : II : open-shell systems

The basis set superposition error-free second-order MØller-Plesset perturbation theory of intermolecular interactions was studied. The difficulties of the counterpoise (CP) correction in open-shell systems were also discussed. The calculations were performed by a program which was used for testing the new variants of the theory. It was shown that the CP correction for the diabatic surfaces should be preferred to the adiabatic ones

Estimation of electronic coupling in π -stacked donor-bridge-acceptor systems: correction of the two-state model

Comparison of donor-acceptor electronic couplings calculated within two-state and three-state models suggests that the two-state treatment can provide unreliable estimates of Vda because of neglecting the multistate effects. We show that in most cases accurate values of the electronic coupling in a π stack, where donor and acceptor are separated by a bridging unit, can be obtained as Ṽ da = (E2 - E1) μ12 Rda + (2 E3 - E1 - E2) 2 μ13 μ23 Rda2, where E1, E2, and E3 are adiabatic energies of the ground, charge-transfer, and bridge states, respectively, μij is the transition dipole moments between the states i and j, and Rda is the distance between the planes of donor and acceptor. In this expression based on the generalized Mulliken-Hush approach, the first term corresponds to the coupling derived within a two-state model, whereas the second term is the superexchange correction accounting for the bridge effect. The formula is extended to bridges consisting of several subunits. The influence of the donor-acceptor energy mismatch on the excess charge distribution, adiabatic dipole and transition moments, and electronic couplings is examined. A diagnostic is developed to determine whether the two-state approach can be applied. Based on numerical results, we showed that the superexchange correction considerably improves estimates of the donor-acceptor coupling derived within a two-state approach. In most cases when the two-state scheme fails, the formula gives reliable results which are in good agreement (within 5%) with the data of the three-state generalized Mulliken-Hush model

Obtenció de nous anàlegs amb activitat brassinoesteroide mitjançant modelització molecular i síntesi

Els brassinoesteroides són productes naturals que actuen com a potents reguladors del creixement vegetal. Presenten aplicacions prometedores en l’agricultura degut a que, aplicats exògenament, augmenten la qualitat i la quantitat de les collites. Ara bé, el seu ús s’ha vist restringit degut a la seva costosa obtenció. Aquest fet ha motivat la recerca de nous compostos actius més assequibles. En aquest projecte es planteja el disseny i obtenció de nous anàlegs seguint diferents estratègies que impliquen tant l’ús de mètodes de modelització molecular com de síntesi orgànica. La primera d’aquestes estratègies consisteix en buscar compostos actius en bases de dades de compostos comercials a través de processos de Virtual Screening desenvolupats amb mètodes computacionals basats en Camps d’Interacció Molecular. Així, es van establir i interpretar models de Relacions Quantitatives Estructura-Activitat (QSAR) emprant descriptors independents de l’alineament (GRIND) i, amb col•laboració amb la Universitat de Perugia, aquest criteri de cerca es va ampliar amb l’aplicació de descriptors FLAP de nova generació. Una altra estratègia es va basar en intentar substituir l’esquelet esteroide dels brassinoesteroides per una estructura equivalent, fixant com a cadena lateral el grup (R)-hexahidromandelil. S’han aplicat dos criteris: mètodes computacionals basats en models QSAR establerts amb descriptors GRIND i també en la metodologia SHOP (scaffold hopping), i, per altra banda, anàlegs proposats racionalment a partir d’un estudi efectuat sobre disruptors endocrins no esteroïdals. Sobre les estructures trobades s’hi va unir la cadena lateral comercial esmentada per via sintètica, en la qual s’ha hagut de fer un èmfasi especial en grups protectors. En total, 49 estructures es proposen per a ser obtingudes sintèticament. També s’ha treballat en l’obtenció un agonista derivat de l’hipotètic antagonista KM-01. Totes les molècules candidates, ja siguin comercials o obtingudes sintèticament, estant sent avaluades en el test d’inclinació de la làmina d’arròs (RLIT).

Development and validation of AMANDA, a new algorithm for selecting highly relevant regions in Molecular Interaction Fields

Descriptors based on Molecular Interaction Fields (MIF) are highly suitable for drug discovery, but their size (thousands of variables) often limits their application in practice. Here we describe a simple and fast computational method that extracts from a MIF a handful of highly informative points (hot spots) which summarize the most relevant information. The method was specifically developed for drug discovery, is fast, and does not require human supervision, being suitable for its application on very large series of compounds. The quality of the results has been tested by running the method on the ligand structure of a large number of ligand-receptor complexes and then comparing the position of the selected hot spots with actual atoms of the receptor. As an additional test, the hot spots obtained with the novel method were used to obtain GRIND-like molecular descriptors which were compared with the original GRIND. In both cases the results show that the novel method is highly suitable for describing ligand-receptor interactions and compares favorably with other state-of-the-art methods.

Suitability of GRIND-based principal properties for the description of molecular similarity and ligand-based virtual screening

The information provided by the alignment-independent GRid Independent Descriptors (GRIND) can be condensed by the application of principal component analysis, obtaining a small number of principal properties (GRIND-PP), which is more suitable for describing molecular similarity. The objective of the present study is to optimize diverse parameters involved in the obtention of the GRIND-PP and validate their suitability for applications, requiring a biologically relevant description of the molecular similarity. With this aim, GRIND-PP computed with a collection of diverse settings were used to carry out ligand-based virtual screening (LBVS) on standard conditions. The quality of the results obtained was remarkable and comparable with other LBVS methods, and their detailed statistical analysis allowed to identify the method settings more determinant for the quality of the results and their optimum. Remarkably, some of these optimum settings differ significantly from those used in previously published applications, revealing their unexplored potential. Their applicability in large compound database was also explored by comparing the equivalence of the results obtained using either computed or projected principal properties. In general, the results of the study confirm the suitability of the GRIND-PP for practical applications and provide useful hints about how they should be computed for obtaining optimum results.

Cannabinoid addiction: behavioral models and neural correlates

The use of cannabis sativa preparations as recreational drugs can be traced back to the earliest civilizations. However, animal models of cannabinoid addiction allowing the exploration of neural correlates of cannabinoid abuse have been developed only recently. We review these models and the role of the CB1 cannabinoid receptor, the main target of natural cannabinoids, and its interaction with opioid and dopamine transmission in reward circuits. Extensive reviews on the molecular basis of cannabinoid action are available elsewhere (Piomelli et al., 2000;Schlicker and Kathmann, 2001).

Matching stages of heavy-ion collision models

Heavy-ion reactions and other collective dynamical processes are frequently described by different theoretical approaches for the different stages of the process, like initial equilibration stage, intermediate locally equilibrated fluid dynamical stage, and final freeze-out stage. For the last stage, the best known is the Cooper-Frye description used to generate the phase space distribution of emitted, noninteracting particles from a fluid dynamical expansion or explosion, assuming a final ideal gas distribution, or (less frequently) an out-of-equilibrium distribution. In this work we do not want to replace the Cooper-Frye description, but rather clarify the ways of using it and how to choose the parameters of the distribution and, eventually, how to choose the form of the phase space distribution used in the Cooper-Frye formula. Moreover, the Cooper-Frye formula is used in connection with the freeze-out problem, while the discussion of transition between different stages of the collision is applicable to other transitions also. More recently, hadronization and molecular dynamics models have been matched to the end of a fluid dynamical stage to describe hadronization and freeze-out. The stages of the model description can be matched to each other on space-time hypersurfaces (just like through the frequently used freeze-out hypersurface). This work presents a generalized description of how to match the stages of the description of a reaction to each other, extending the methodology used at freeze-out, in simple covariant form which is easily applicable in its simplest version for most applications.

Molecular phylogenetics and temporal diversification in the genus Aeromonas based on the sequences of five housekeeping genes

Several approaches have been developed to estimate both the relative and absolute rates of speciation and extinction within clades based on molecular phylogenetic reconstructions of evolutionary relationships, according to an underlying model of diversification. However, the macroevolutionary models established for eukaryotes have scarcely been used with prokaryotes. We have investigated the rate and pattern of cladogenesis in the genus Aeromonas (γ-Proteobacteria, Proteobacteria, Bacteria) using the sequences of five housekeeping genes and an uncorrelated relaxed-clock approach. To our knowledge, until now this analysis has never been applied to all the species described in a bacterial genus and thus opens up the possibility of establishing models of speciation from sequence data commonly used in phylogenetic studies of prokaryotes. Our results suggest that the genus Aeromonas began to diverge between 248 and 266 million years ago, exhibiting a constant divergence rate through the Phanerozoic, which could be described as a pure birth process.

Definition of linear color models in the RGB vector color space to detect red peaches in orchard images taken under natural illumination

This work proposes the detection of red peaches in orchard images based on the definition of different linear color models in the RGB vector color space. The classification and segmentation of the pixels of the image is then performed by comparing the color distance from each pixel to the different previously defined linear color models. The methodology proposed has been tested with images obtained in a real orchard under natural light. The peach variety in the orchard was the paraguayo (Prunus persica var. platycarpa) peach with red skin. The segmentation results showed that the area of the red peaches in the images was detected with an average error of 11.6%; 19.7% in the case of bright illumination; 8.2% in the case of low illumination; 8.6% for occlusion up to 33%; 12.2% in the case of occlusion between 34 and 66%; and 23% for occlusion above 66%. Finally, a methodology was proposed to estimate the diameter of the fruits based on an ellipsoidal fitting. A first diameter was obtained by using all the contour pixels and a second diameter was obtained by rejecting some pixels of the contour. This approach enables a rough estimate of the fruit occlusion percentage range by comparing the two diameter estimates.

The embryonic blood-CSF barrier has molecular elements for specific glucose transport and for the general transport of molecules via transcellular routes.

In vertebrates, early brain development takes place at the expanded anterior end of the neural tube, which is filled with embryonic cerebrospinal fluid (E-CSF). We have recently identified a transient blood-CSF barrier that forms between embryonic days E3 and E4 in chick embryos and that is responsible for the transport of proteins and control of E-CSF homeostasis, including osmolarity. Here we examined the presence of glucose transporter GLUT-1 as well the presence of caveolae-structural protein Caveolin1 (CAV-1) in the embryonic blood-CSF barrier which may be involved in the transport of glucose and of proteins, water and ions respectively across the neuroectoderm. In this paper we demonstrate the presence of GLUT-1 and CAV-1 in endothelial cells of blood vessels as well as in adjacent neuroectodermal cells, located in the embryonic blood-CSF barrier. In blood vessels, these proteins were detected as early as E4 in chick embryos and E12.7 in rat embryos, i.e. the point at which the embryonic blood-CSF barrier acquires this function. In the neuroectoderm of the embryonic blood-CSF barrier, GLUT-1 was also detected at E4 and E12.7 respectively, and CAV-1 was detected shortly thereafter in both experimental models. These experiments contribute to delineating the extent to which the blood-CSF embryonic barrier controls E-CSF composition and homeostasis during early stages of brain development in avians and mammals. Our results suggest the regulation of glucose transport to the E-CSF by means of GLUT-1 and also suggest a mechanism by which proteins are transported via transcellular routes across the neuroectoderm, thus reinforcing the crucial role of E-CSF in brain development.

Molecular phylogenetics and temporal diversification in the genus Aeromonas based on the sequences of five housekeeping genes

Several approaches have been developed to estimate both the relative and absolute rates of speciation and extinction within clades based on molecular phylogenetic reconstructions of evolutionary relationships, according to an underlying model of diversification. However, the macroevolutionary models established for eukaryotes have scarcely been used with prokaryotes. We have investigated the rate and pattern of cladogenesis in the genus Aeromonas (γ-Proteobacteria, Proteobacteria, Bacteria) using the sequences of five housekeeping genes and an uncorrelated relaxed-clock approach. To our knowledge, until now this analysis has never been applied to all the species described in a bacterial genus and thus opens up the possibility of establishing models of speciation from sequence data commonly used in phylogenetic studies of prokaryotes. Our results suggest that the genus Aeromonas began to diverge between 248 and 266 million years ago, exhibiting a constant divergence rate through the Phanerozoic, which could be described as a pure birth process.

Molecular dynamics simulation of the spherical electrical double layer of a soft nanoparticle. Effect of the surface and counterion valence

Molecular dynamics simulations were performed to study the ion and water distribution around a spherical charged nanoparticle. A soft nanoparticle model was designed using a set of hydrophobic interaction sites distributed in six concentric spherical layers. In order to simulate the effect of charged functionalyzed groups on the nanoparticle surface, a set of charged sites were distributed in the outer layer. Four charged nanoparticle models, from a surface charge value of −0.035 Cm−2 to − 0.28 Cm−2, were studied in NaCl and CaCl2 salt solutions at 1 M and 0.1 M concentrations to evaluate the effect of the surface charge, counterion valence, and concentration of added salt. We obtain that Na + and Ca2 + ions enter inside the soft nanoparticle. Monovalent ions are more accumulated inside the nanoparticle surface, whereas divalent ions are more accumulated just in the plane of the nanoparticle surface sites. The increasing of the the salt concentration has little effect on the internalization of counterions, but significantly reduces the number of water molecules that enter inside the nanoparticle. The manner of distributing the surface charge in the nanoparticle (uniformly over all surface sites or discretely over a limited set of randomly selected sites) considerably affects the distribution of counterions in the proximities of the nanoparticle surface.