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.

SPH simulations of clumps formation by dissipative collisions of molecular clouds - II. Magnetic case

We performed computer simulations of interstellar cloud-cloud collisions using the three-dimensional smoothed particle magnetohydrodynamics method. In order to study the role of the magnetic field on the process of collision-triggered fragmentation, we focused our attention on head-on supersonic collisions between two identical spherical molecular-clouds. Two extreme configurations of the magnetic field were adopted: parallel and perpendicular to the initial clouds motion. The initial magnetic field strength was approximately 12.0 muG. In the parallel case, much more of the collision debris were retained in the shocking region than in the non-magnetic case where gas escaped freely throughout the symmetry plane. Differently from the non-magnetic case, eddy-like vortices were formed. The regions of highest vorticity and the the regions of highest density are offset. We found clumps formation only in the parallel case, however, they were larger, hotter and less dense than in the analogous non-magnetic case. In the perpendicular case, the compressed field works as a magnetic wall, preventing a stronger compression of the colliding clouds. This last effect inhibits direct contact of the two clouds. In both cases, we found that the field lines show a chaotic aspect in large scales. Also, the field magnitude is considerably amplified in the shock layer. However, the field distribution is almost coherent in the higher density regions.

SPH simulations of clumps formation by dissipative collision of molecular clouds - I. Non magnetic case

Computer experiments of interstellar cloud collisions were performed with a new smoothed-particle-hydrodynamics (SPH) code. The SPH quantities were calculated by using spatially adaptive smoothing lengths and the SPH fluid equations of motion were solved by means of a hierarchical multiple time-scale leapfrog. Such a combination of methods allows the code to deal with a large range of hydrodynamic quantities. A careful treatment of gas cooling by H, H(2), CO and H II, as well as a heating mechanism by cosmic rays and by H(2) production on grains surface, were also included in the code. The gas model reproduces approximately the typical environment of dark molecular clouds. The experiments were performed by impinging two dynamically identical spherical clouds onto each other with a relative velocity of 10 km s(-1) but with a different impact parameter for each case. Each object has an initial density profile obeying an r(-1)-law with a cutoff radius of 10 pc and with an initial temperature of 20 K. As a main result, cloud-cloud collision triggers fragmentation but in expense of a large amount of energy dissipated, which occurred in the head-on case only. Off-center collision did not allow remnants to fragment along the considered time (similar to 6 Myr). However, it dissipated a considerable amount of orbital energy. Structures as small as 0.1 pc, with densities of similar to 10(4) cm(-3), were observed in the more energetic collision.

Molecular dynamics simulations on devitrification: the PbF2 case

In this work molecular dynamics simulations were performed to reproduce the kinetic and thermodynamic transformations occurring during melt crystallization, vitrification, and glass crystallization (devitrification) of PbF2. Two potential parameters were analyzed in order to access the possibility of modeling these properties. These interionic potentials are models developed to describe specific characteristic of PbF2, and thermodynamic properties were well reproduced by one of them, while the other proved well adapted to simulate the crystalline structure of this fluoride. By a modeled nonisothermal heat treatment of the glass, it was shown that the devitrification of a cubic structure in which the Pb-Pb distances are in good agreement with theory and experiment. (C) 2002 American Institute of Physics.

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.

On the unzipping mechanisms of carbon nanotubes: Insights from reactive molecular dynamics simulations

Unzipping carbon nanotubes (CNTs) is considered one of the most promising approaches for the controlled and large-scale production of graphene nanoribbons (GNR). These structures are considered of great importance for the development of nanoelectronics because of its dimensions and intrinsic nonzero band gap value. Despite many years of investigations some details on the dynamics of the CNT fracture/unzipping processes remain unclear. In this work we have investigated some of these process through molecular dynamics simulations using reactive force fields (ReaxFF), as implemented in the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) code. We considered multi-walled CNTs of different dimensions and chiralities and under induced mechanical stretching. Our preliminary results show that the unzipping mechanisms are highly dependent on CNT chirality. Well-defined and distinct fracture patterns were observed for the different chiralities. Armchair CNTs favor the creation of GNRs with well-defined armchair edges, while zigzag and chiral ones produce GNRs with less defined and defective edges. © 2012 Materials Research Society.

Assembling a xylanase-lichenase chimera through all-atom molecular dynamics simulations

Multifunctional enzyme engineering can improve enzyme cocktails for emerging biofuel technology. Molecular dynamics through structure-based models (SB) is an effective tool for assessing the tridimensional arrangement of chimeric enzymes as well as for inferring the functional practicability before experimental validation. This study describes the computational design of a bifunctional xylanase-lichenase chimera (XylLich) using the xynA and bglS genes from Bacillus subtilis. In silico analysis of the average solvent accessible surface area (SAS) and the root mean square fluctuation (RMSF) predicted a fully functional chimera, with minor fluctuations and variations along the polypeptide chains. Afterwards, the chimeric enzyme was built by fusing the xynA and bglS genes. XylLich was evaluated through small-angle X-ray scattering (SAXS) experiments, resulting in scattering curves with a very accurate fit to the theoretical protein model. The chimera preserved the biochemical characteristics of the parental enzymes, with the exception of a slight variation in the temperature of operation and the catalytic efficiency (k cat/Km). The absence of substantial shifts in the catalytic mode of operation was also verified. Furthermore, the production of chimeric enzymes could be more profitable than producing a single enzyme separately, based on comparing the recombinant protein production yield and the hydrolytic activity achieved for XylLich with that of the parental enzymes. © 2013 Elsevier B.V. All rights reserved.

Combining Experimental Evidence and Molecular Dynamic Simulations To Understand the Mechanism of Action of the Antimicrobial Octapeptide Jelleine-I

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

Role of peptide peptide interactions in aggregation: protonectins observed in equilibrium and replica exchange molecular dynamics simulations

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

Structural and functional insights into the conductive pili of Geobacter sulfurreducens revealed in molecular dynamics simulations

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

SMase II, a new sphingomyelinase D from Loxosceles laeta venom gland: Molecular cloning, expression, function and structural analysis

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

Quantum coherent tunneling between two atomic-molecular Bose-Einstein condensates

We study the quantum coherent tunneling dynamics of two weakly coupled atomic-molecular Bose-Einstein condensates (AMBEC). A weak link is supposed to be provided by a double-well trap. The regions of parameters where the macroscopic quantum localization of the relative atomic population occurs are revealed. The different dynamical regimes are found depending on the value of nonlinearity, namely, coupled oscillations of population imbalance of atomic and molecular condensate, including irregular oscillations regions, and macroscopic quantum self trapping regimes. Quantum means and quadrature variances are calculated for population of atomic and molecular condensates and the possibility of quadrature squeezing is shown via stochastic simulations within P-positive phase space representation method. Linear tunnel coupling between two AMBEC leads to correlations in quantum statistics.

Molecular Docking and Molecular Dynamic Studies of Semi-Synthetic Piperidine Alkaloids as Acetylcholinesterase Inhibitors

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

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.

Theoretical studies of liquids by computer simulations: the water-acetone mixture.

Monte Carlo simulation results for pure liquid acetone and water-acetone mixtures calculated in the isothermal and isobaric (NPT) ensemble at T=298K and p=1.0atm are presented. The TIP4P model was used for water and optimized potential for liquid simulation (OPLS) force field parameters used for acetone. The results obtained for the average configurational energy as a function of the mole fraction are in good accord with experimental data. Energy partitioning and co-ordination numbers results calculated for equimolar water-acetone solution are compared to similar data obtained for other water-organic liquid mixtures. These results show an increase in water-water interaction energy and co-ordination numbers when the interaction between water and organic liquid molecules decrease. Distribution functions for pure liquid acetone and water-acetone mixtures are presented. Dipole-dipole angular correlation functions obtained for pure liquid acetone show a predominance of dimers with parallel alignment of dipole moments. Radial distribution functions from water-acetone interaction show characteristic features of hydrogen bonded liquids. Radial and angular distribution functions for water-water correlation calculated in pure water and in equimolar water-acetone mixture are compared, showing very similar features in both systems. (C) 1999 Elsevier B.V. B.V. All rights reserved.