MONTE-CARLO SIMULATION OF POLYAMPHOLYTE CHAINS

Polyampholyte copolymers containing both positive and negative monomers regularly dispersed along the chain were studied. The Monte Carlo method was used to simulate chains with charged monomers interacting by screened Coulomb potential. The neutral polyampholyte chains collapse due to the attractive electrostatic interactions. The nonneutral chains are in extended conformations due to the repulsive polyelectrolyte effects that dominate the attractive polyampholyte interactions. The results are in good agreement with experiment.

Modelling polymers with side chains: MEH-PPV and P3HT

Modelling polymers with side chains is always a challenge once the degrees of freedom are very high. In this study, we present a successful methodology to model poly[2-methoxy-5-(2′-ethyl-hexyloxy)-p-phenylenevinylene] (MEH-PPV) and poly[3-hexylthiophene] (P3HT) in solutions, taking into account the influence of side chains on the polymer conformation. Molecular dynamics and semi-empirical quantum mechanical methods were used for structure optimisation and evaluation of optical properties. The methodology allows to describe structural and optical characteristics of the polymers in a satisfactory way, as well as to evaluate some usual simplifications adopted for modelling these systems. Effective conjugation lengths of 8-14.6 and 21 monomers were obtained for MEH-PPV and P3HT, respectively, in accordance with experimental findings. In addition, anti/syn conformations of these polymers could be predicted based on intrinsic interactions of the lateral branches. © 2013 Copyright Taylor and Francis Group, LLC.

Quantum chemical topological analysis of hydrogen bonding in HX…HX and CH3X…HX dimers (X = Br, Cl, F)

We present a systematic investigation of the nature and strength of the hydrogen bonding in HX···HX and CH3X…HX (X = Br, Cl and F) dimers using ab initio MP2/aug-cc-pVTZ calculations in the framework of the quantum theory of atoms in molecules (QTAIM) and electron localisation functions (ELFs) methods. The electron density of the complexes has been characterised, and the hydrogen bonding energy, as well as the QTAIM and ELF parameters, is consistent, providing deep insight into the origin of the hydrogen bonding in these complexes. It was found that in both linear and angular HX…HX and CH3X…HX dimers, F atoms form stronger HB than Br and Cl, but they need short (∼2 Å) X…HX contacts.

Molecular dynamics simulation on devitrification: Isothermal devitrification and thermodynamics of PbF2 glasses

The vitrification and devitrification features of lead fluoride are investigated by means of molecular dynamic simulations. The influence of heating rate on the devitrification temperature as well as the dependence of the glass properties on its thermal history, i.e., the cooling rate employed, is identified. As expected, different glasses are obtained when the cooling rates differ. Diffusion coefficient analysis during heating of glass and crystal, indicates that the presence of defects on the glassy matrix favors the transition processes from the ionic to a superionic state, with high mobility of fluorine atoms, responsible for the high anionic conduction of lead fluoride. Nonisothermal and isothermal devitrification processes are simulated in glasses obtained at different cooling rates and structural organizations occurring during the heat treatments are clearly observed. When a fast cooling rate is employed during the glass formation, the devitrification of a single crystal (limited by the cell dimensions) is observed, while the glass obtained with slower cooling rate, allowing relaxations and organization of various regions on the glass bulk during the cooling process, devitrifies in more than one crystalline plane. (C) 2004 American Institute of Physics.

The use of molecular spectra simulation for diagnostics of reactive flows

The C 2 * radical is used as a system probe tool to the reactive flow diagnostic, and it was chosen due to its large occurrence in plasma and combustion in aeronautics and aerospace applications. The rotational temperatures of C 2 * species were determined by the comparison between experimental and theoretical data. The simulation code was developed by the authors, using C++ language and the object oriented paradigm, and it includes a set of new tools that increase the efficacy of the C 2 * probe to determine the rotational temperature of the system. A brute force approach for the determination of spectral parameters was adopted in this version of the computer code. The statistical parameter c 2 was used as an objective criterion to determine the better match of experimental and synthesized spectra. The results showed that the program works even with low-quality experimental data, typically collected from in situ airborne compact apparatus. The technique was applied to flames of a Bunsen burner, and the rotational temperature of ca. 2100 K was calculated.

Distribution of neutral prilocaine in a phospholipid bilayer: Insights from molecular dynamics simulations

In this work, we report a 20-ns constant pressure molecular dynamics simulation of prilocaine (PLC), in amine-amide local anesthetic, in a hydrated liquid crystal bilayer of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine. The partition of PLC induces the lateral expansion of the bilayer and a concomitant contraction in its thickness. PLC molecules are preferentially found in the hydrophobic acyl chains region, with a maximum probability at similar to 12 angstrom from the center of the bilayer (between the C(4) and C(5) methylene groups). A decrease in the acyl chain segmental order parameter, vertical bar S-CD vertical bar, compared to neat bilayers, is found, in good agreement with experimental H-2-NMR studies. The decrease in vertical bar S-CD vertical bar induced by PLC is attributed to a larger accessible volume per lipid in the acyl chain region. (C) 2008 Wiley Periodicals, Inc.

Preferential location of prilocaine and etidocaine in phospholipid bilayers: A molecular dynamics study

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Molecular characterization and phylogenetic analysis of JussuMP-I: A RGD-P-III class hemorrhagic metalloprotease from Bothrops jararacussu snake venom

Snake venom metalloproteases (SVMPs) embody zinc-dependent multidomain enzymes responsible for a relevant pathophysiology in envenomation. including local and systemic hemorrhage. The molecular features responsible for hemorrhagic potency of SVMPs have been associated with their multidomains structures which can target these proteins them to several receptors of different tissues and cellular types. BjussuMP-I. a SVMP isolated from the Bothrops jararacussu venom, has been characterized as a P-III hemorrhagic metalloprotease. The complete cDNA sequence of BjussuMP-I with 1641bp encodes open reading frames of 547 amino acid residues, which conserve the common domains of P-III high molecular weight hemorrhagic metalloproteases: (i) pre-pro-peptide, (ii) metalloprotease, (iii) disintegrin-like and (iv) rich cysteine domain. BjussuMP-I induced lyses in fibrin clots and inhibited collagen- and ADP-induced platelet aggregation. We are reporting, for the first time, the primary structure of an RGD-P-III class snake venom metalloprotease. A phylogenetic analysis of the BjussuMP-1 metalloprotease/catalytic domain was performed to get new insights into the molecular evolution of the metalloproteases. A theoretical molecular model of this domain was built through folding recognition (threading) techniques and refined by molecular dynamics simulation. Then, the final BjussuMP-I catalytic domain model was compared to other SVMPs and Reprolysin family proteins in order to identify eventual structural differences, which could help to understand the biochemical activities of these enzymes. The presence of large hydrophobic areas and some conserved surface charge-positive residues were identified as important features of the SVMPs and other metalloproteases. (C) 2006 Elsevier B.V. All rights reserved.

A deterministic simulation study of embryo marker-assisted selection for age at first calving in Nellore (Bos indicus) beef cattle

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Numerical simulation of magnetic-field-enhanced plasma immersion ion implantation in cylindrical geometry

Recent studies have demonstrated that the sheath dynamics in plasma immersion ion implantation (PIII) is significantly affected by an external magnetic field. In this paper, a two-dimensional computer simulation of a magnetic-field-enhanced PHI system is described. Negative bias voltage is applied to a cylindrical target located on the axis of a grounded vacuum chamber filled with uniform molecular nitrogen plasma. A static magnetic field is created by a small coil installed inside the target holder. The vacuum chamber is filled with background nitrogen gas to form a plasma in which collisions of electrons and neutrals are simulated by the Monte Carlo algorithm. It is found that a high-density plasma is formed around the target due to the intense background gas ionization by the magnetized electrons drifting in the crossed E x B fields. The effect of the magnetic field intensity, the target bias, and the gas pressure on the sheath dynamics and implantation current of the PHI system is investigated.

C-H center dot center dot center dot O and N-H center dot center dot center dot O hydrogen bonds in liquid amides investigated by Monte Carlo simulation

Monte Carlo simulations of liquid formamide, N-methylformamide (MF), and N,N-dimethylformamide (DMF) have been performed in the isothermal and isobaric ensemble at 298 K and 1 atm, aiming to investigate the C-H ... O and N-H ... O hydrogen bonds. The interaction energy was calculated using the classical 6-12 Lennard-Jones pairwise potential plus a Coulomb term on a rigid six-site molecular model with the potential parameters being optimized in this work. Theoretical values obtained for heat of vaporization and liquid densities are in good agreement with the experimental data. The radial distribution function [RDF, g(r)] obtained compare well with R-X diffraction data available. The RDF and molecular mechanics (MM2) minimization show that the C-H ... O interaction has a significant role in the structure of the three liquids. These results are supported by ab initio calculations. This Interaction is particularly important in the structure of MF. The intensity of the N-H ... O hydrogen bond is greater in the MF than formamide. This could explain some anomalous properties verified in MF. (C) 1997 John Wiley & Sons, Inc.

C - H ⋯ 0 and N - H ⋯ 0 hydrogen bonds in liquid amides investigated by monte carlo simulation

Monte Carlo simulations of liquid formamide, N-methylformamide (MF), and N, N-dimethytformamide (DMF) have been performed in the isothermal and isobaric ensemble at 298 K and 1 atm, aiming to investigate the C-H ⋯ O and N-H ⋯ O hydrogen bonds. The interaction energy was calculated using the classical 6-12 Lennard-Jones pairwise potential plus a Coulomb term on a rigid six-site molecular model with the potential parameters being optimized in this work. Theoretical values obtained for heat of vaporization and liquid densities are in good agreement with the experimental data. The radial distribution function [RDF, g(r)] obtained compare well with R-X diffraction data available. The RDF and molecular mechanics (MM2) minimization show that the C-H ⋯ O interaction has a significant role in the structure of the three liquids. These results are supported by ab initio calculations. This interaction is particularly important in the structure of MF. The intensity of the N - H ⋯ O hydrogen bond is greater in the MF than formamide. This could explain some anomalous properties verified in MF. © 1997 John Wiley & Sons, Inc.

Caracterização silvicultural, botânica e avaliação da variabilidade genética por meio do marcador molecular RAPD em um teste de progênies de Eucalyptus urophylla S. T. Blake

Molecular markers have recently been incorporated into genetic improvement programs. They are already considered as powerful tools with several different uses, for instance the monitoring of genetic variability in tree populations. The main objectives of this study were to evaluate genetic variability in Eucalyptus urophylla progenies and together with silvicultural and botanical information, provide assistance to the improvement program. The Eucalypts population is located at the Experimental Forestry Sciences Station, Anhembi, SP, which belongs to the College of Agriculture Luiz de Queiroz. Sixty-nine progenies were analysed representing one individual by family in open pollinated Eucalyptus urophylla trees. The RAPD technique allowed the identification of 72 loci that were analysed using Jaccard's Coefficient generating a genetic similarity matrix to permit estimation of genetic distances. The results obtained showed genetic distance between individuals of 0.40 with 12 groups of genetic variability using a standardised distance of 40%. The progenies showed different bark patterns, allowing the establishment of bark groups. The groups formed based on genetic distances obtained using DNA analysis did not correspond to those based on bark pattern. Genetic selection was simulated in which silvicultural and genetic variability data were linked, thus avoiding excessive variability losses. The simulation of controlled crossings allowed the maximum genetic difference to be obtained linked with height and individual bark roughness.

Molecular modeling databases: A new way in the search of protein targets for drug development

DBMODELING is a relational database of annotated comparative protein structure models and their metabolic, pathway characterization. It is focused on enzymes identified in the genomes of Mycobacterium tuberculosis and Xylella fastidiosa. The main goal of the present database is to provide structural models to be used in docking simulations and drug design. However, since the accuracy of structural models is highly dependent on sequence identity between template and target, it is necessary to make clear to the user that only models which show high structural quality should be used in such efforts. Molecular modeling of these genomes generated a database, in which all structural models were built using alignments presenting more than 30% of sequence identity, generating models with medium and high accuracy. All models in the database are publicly accessible at http://www.biocristalografia.df.ibilce.unesp.br/tools. DBMODELING user interface provides users friendly menus, so that all information can be printed in one stop from any web browser. Furthermore, DBMODELING also provides a docking interface, which allows the user to carry out geometric docking simulation, against the molecular models available in the database. There are three other important homology model databases: MODBASE, SWISSMODEL, and GTOP. The main applications of these databases are described in the present article. © 2007 Bentham Science Publishers Ltd.

Molecular Dynamics simulations of track formation at different ensembles

A series of Molecular Dynamics simulations of thermal spikes has been run in zircon. For two different ensembles: microcanonical one and a combination of microcanonical one acting on the simulation core with Langevin one on the side walls of simulation. Depending on the used ensemble, different track-formation threshold energies were found. When the combined ensemble is carried out, the total energy of the simulations varies with the temperature which can influence how annealing fission-track models should deal with the lattice recovery. A fission-track annealing model is tested with the simulation results. © 2012 Elsevier Ltd. All rights reserved.