Study of the antimicrobial peptide indolicidin and a mutant in micelle medium by molecular dynamics simulation

The antimicrobial peptide indolicidin (IND) and the mutant CP10A in hydrated micelles were studied using molecular dynamics simulations in order to observe whether the molecular dynamics and experimental data could be sufficiently correlated and a detailed description of the interaction of the antimicrobial peptides with a model of the membrane provided by a hydrated micelle system could be obtained. In agreement with the experiments, the simulations showed that the peptides are located near the surface of the micelles. Peptide insertions agree with available experimental data, showing deeper insertion of the mutant compared with the peptide IND. Major insertion into the hydrophobic core of the micelle by all tryptophan and mutated residues of CP10A in relation to IND was observed. The charged residues of the terminus regions of both peptides present similar behavior, indicating that the major differences in the interactions with the micelles of the peptides IND and CP10A occur in the case of the hydrophobic residues.

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]

The shoving model for the glass-former LiCl center dot 6H(2)O: A molecular dynamics simulation study

Molecular dynamics (MD) simulations of LiCl center dot 6H(2)O Showed that the diffusion coefficient D, and also I lie structural relaxation time , follow a power law at high temperatures, D(-1) proportional to (T - T(0))(-mu), with the same experimental parameters for viscosity (T(0) = 207 K, mu = 2.08). Decoupling between D and occurs at T(x) similar to 1.1 T(0). High frequency acoustic excitations for the LiCl center dot 6H(2)O model were obtained by the calculation of time correlation functions of mass current fluctuations. The temperature dependence of the instantaneous shear modulus, G,(T), was considered in the shoving model for supercooled liquids [J.C. Dyre, T. Christensen, N.B. Olsen, J. Non-Cryst. Solids 352 (2006) 4635] resulting in a linear relationship log (D(-1)) vs. G root T. (C) 2009 Elsevier B.V. All rights reserved.

Docking and molecular dynamics simulation of quinone compounds with trypanocidal activity

In this work, two different docking programs were used, AutoDock and FlexX, which use different types of scoring functions and searching methods. The docking poses of all quinone compounds studied stayed in the same region in the trypanothione reductase. This region is a hydrophobic pocket near to Phe396, Pro398 and Leu399 amino acid residues. The compounds studied displays a higher affinity in trypanothione reductase (TR) than glutathione reductase (GR), since only two out of 28 quinone compounds presented more favorable docking energy in the site of human enzyme. The interaction of quinone compounds with the TR enzyme is in agreement with other studies, which showed different binding sites from the ones formed by cysteines 52 and 58. To verify the results obtained by docking, we carried out a molecular dynamics simulation with the compounds that presented the highest and lowest docking energies. The results showed that the root mean square deviation (RMSD) between the initial and final pose were very small. In addition, the hydrogen bond pattern was conserved along the simulation. In the parasite enzyme, the amino acid residues Leu399, Met400 and Lys402 are replaced in the human enzyme by Met406, Tyr407 and Ala409, respectively. In view of the fact that Leu399 is an amino acid of the Z site, this difference could be explored to design selective inhibitors of TR.

Study of the antimicrobial peptide indolicidin and mutants in eukaryotic modelled membrane by molecular dynamics simulations

In this work the interaction of the antimicrobial peptide indolicidin (IND) and its mutants CP10A and CP11 with a eukaryotic membrane model was examined by molecular dynamics simulations. The aim was to analyse the behaviour of these antimicrobial peptides when they interact with a eukaryotic modelled membrane, thereby obtaining atomic detailed observations that are not experimentally available. In the simulations, the widely studied dipalmitoylphosphatidylcholine hydrated bilayer was used as a eukaryotic membrane model. In agreement with experimental observations, the peptides IND, CP10A, and CP11 insert into the bilayer differently; the peptides that insert more deeply present the major hemolytic activities. The hydrophobic residues are responsible for the insertion, but some Trp residues of the peptides remain at the bilayer/water interface because they interact with the bilayer choline groups by cation-pi interactions that should be important for recognition of eukaryotic membrane by the three studied peptides.

Molecular Basis of the Thermostability and Thermophilicity of Laminarinases: X-ray Structure of the Hyperthermostable Laminarinase from Rhodothermus marinus and Molecular Dynamics Simulations

Glycosyl hydrolases are enzymes capable of breaking the glycosidic linkage of polysaccharides and have considerable industrial and biotechnological applications. Driven by the later applications, it is frequently desirable that glycosyl hydrolases display stability and activity under extreme environment conditions, such as high temperatures and extreme pHs. Here, we present X-ray structure of the hyperthermophilic laminarinase from Rhodothermus marinus (RmLamR) determined at 1.95 angstrom resolution and molecular dynamics simulation studies aimed to comprehend the molecular basis, for the thermal stability of this class of enzymes. As most thermostable proteins, RmLamR contains a relatively large number of salt bridges, which are not randomly distributed on the structure. On the contrary, they form clusters interconnecting beta-sheets of the catalytic domain. Not all salt bridges, however, are beneficial for the protein thermostability: the existence of charge-charge interactions permeating the hydrophobic core of the enzymes actually contributes to destabilize the structure by facilitating water penetration into hydrophobic cavities, as can be seen in the case of mesophilic enzymes. Furthermore, we demonstrate that the mobility of the side-chains is perturbed differently in each class of enzymes. The side-chains of loop residues surrounding the catalytic cleft in the mesophilic laminarinase gain mobility and obstruct the active site at high temperature. By contrast, thermophilic laminarinases preserve their active site flexibility, and the active-site cleft remains accessible for recognition of polysaccharide substrates even at high temperatures. The present results provide structural insights into the role played by salt-bridges and active site flexibility on protein thermal stability and may be relevant for other classes of proteins, particularly glycosyl hydrolases.

Phase transitions and spatially ordered counterion association in ionic-lipid membranes: A statistical model

We propose a statistical model to account for the gel-fluid anomalous phase transitions in charged bilayer- or lamellae-forming ionic lipids. The model Hamiltonian comprises effective attractive interactions to describe neutral-lipid membranes as well as the effect of electrostatic repulsions of the discrete ionic charges on the lipid headgroups. The latter can be counterion dissociated (charged) or counterion associated (neutral), while the lipid acyl chains may be in gel (low-temperature or high-lateral-pressure) or fluid (high-temperature or low-lateral-pressure) states. The system is modeled as a lattice gas with two distinct particle types-each one associated, respectively, with the polar-headgroup and the acyl-chain states-which can be mapped onto an Ashkin-Teller model with the inclusion of cubic terms. The model displays a rich thermodynamic behavior in terms of the chemical potential of counterions (related to added salt concentration) and lateral pressure. In particular, we show the existence of semidissociated thermodynamic phases related to the onset of charge order in the system. This type of order stems from spatially ordered counterion association to the lipid headgroups, in which charged and neutral lipids alternate in a checkerboard-like order. Within the mean-field approximation, we predict that the acyl-chain order-disorder transition is discontinuous, with the first-order line ending at a critical point, as in the neutral case. Moreover, the charge order gives rise to continuous transitions, with the associated second-order lines joining the aforementioned first-order line at critical end points. We explore the thermodynamic behavior of some physical quantities, like the specific heat at constant lateral pressure and the degree of ionization, associated with the fraction of charged lipid headgroups.

Mechanism of point-defect diffusion in a two-dimensional colloidal crystal

The dynamics and mechanism of migration of a vacancy point defect in a two-dimensional (2D) colloidal crystal are studied using numerical simulations. We find that the migration of a vacancy is always realized by topology switching between its different configurations. From the temperature dependence of the topology switch frequencies, we obtain the activation energies for possible topology transitions associated with the vacancy diffusion in the 2D crystal. (C) 2011 American Institute of Physics. [doi:10.1063/1.3615287]

Density functional theory investigation of 3d, 4d, and 5d 13-atom metal clusters

The knowledge of the atomic structure of clusters composed by few atoms is a basic prerequisite to obtain insights into the mechanisms that determine their chemical and physical properties as a function of diameter, shape, surface termination, as well as to understand the mechanism of bulk formation. Due to the wide use of metal systems in our modern life, the accurate determination of the properties of 3d, 4d, and 5d metal clusters poses a huge problem for nanoscience. In this work, we report a density functional theory study of the atomic structure, binding energies, effective coordination numbers, average bond lengths, and magnetic properties of the 3d, 4d, and 5d metal (30 elements) clusters containing 13 atoms, M(13). First, a set of lowest-energy local minimum structures (as supported by vibrational analysis) were obtained by combining high-temperature first- principles molecular-dynamics simulation, structure crossover, and the selection of five well-known M(13) structures. Several new lower energy configurations were identified, e. g., Pd(13), W(13), Pt(13), etc., and previous known structures were confirmed by our calculations. Furthermore, the following trends were identified: (i) compact icosahedral-like forms at the beginning of each metal series, more opened structures such as hexagonal bilayerlike and double simple-cubic layers at the middle of each metal series, and structures with an increasing effective coordination number occur for large d states occupation. (ii) For Au(13), we found that spin-orbit coupling favors the three-dimensional (3D) structures, i.e., a 3D structure is about 0.10 eV lower in energy than the lowest energy known two-dimensional configuration. (iii) The magnetic exchange interactions play an important role for particular systems such as Fe, Cr, and Mn. (iv) The analysis of the binding energy and average bond lengths show a paraboliclike shape as a function of the occupation of the d states and hence, most of the properties can be explained by the chemistry picture of occupation of the bonding and antibonding states.

HFE gene mutations in patients with primary iron overload: Is there a significant improvement in molecular diagnosis yield with HFE sequencing?

Rare HFE variants have been shown to be associated with hereditary hemochromatosis (HH), an iron overload disease. The low frequency of the HFE p.C282Y mutation in HH-affected Brazilian patients may suggest that other HFE-related mutations may also be implicated in the pathogenesis of HH in this population. The main aim was to screen for new HFE mutations in Brazilian individuals with primary iron overload and to investigate their relationship with HH. Fifty Brazilian patients with primary iron overload (transferrin saturation >50% in females and 60% in males) were selected. Subsequent bidirectional sequencing for each HFE exon was performed. The effect of HFE mutations on protein structure were analyzed by molecular dynamics simulation and free binding energy calculations. p.C282Y in homozygosis or in heterozygosis with p.H63D were the most frequent genotypic combinations associated with HH in our sample population (present in 17 individuals, 34%). Thirty-six (72.0%) out of the 50 individuals presented at least one HFE mutation. The most frequent genotype associated with HH was the homozygous p.C282Y mutation (n = 11, 22.0%). One novel mutation (p.V256I) was indentified in heterozygosis with the p.H63D mutation. In silico modeling analysis of protein behavior indicated that the p.V256I mutation does not reduce the binding affinity between HFE and beta 2-microglobulin ((beta 2M) in the same way the p.C282Y mutation does compared with the native HFE protein. In conclusion, screening of HFE through direct sequencing, as compared to p.C282Y/p.H63D genotyping, was not able to increase the molecular diagnosis yield of HH. The novel p.V256I mutation could not be implicated in the molecular basis of the HH phenotype, although its role cannot be completely excluded in HH-phenotype development. Our molecular modeling analysis can help in the analysis of novel, previously undescribed, HFE mutations. (C) 2010 Elsevier Inc. All rights reserved.

Structural and functional properties of Bp-LAAO, a new L-amino acid oxidase isolated from Bothrops pauloensis snake venom

An L-amino acid oxidase (Bp-LAAO) from Bothrops pauloensis snake venom was highly purified using sequential chromatography steps on CM-Sepharose, Phenyl-Sepharose CL4B, Benzamidine Sepharose and C18 reverse-phase HPLC. Purified Bp-LAAO showed to be a homodimeric acidic glycoprotein with molecular weight around 65 kDa under reducing conditions in SDS-PAGE. The best substrates for Bp-LAAO were L-Met, L-Leu, L-Phe and L-Ile and the enzyme showed a strong reduction of its catalytic activity upon L-Met and L-Phe substrates at extreme temperatures. Bp-LAAO showed leishmanicidal, antitumoral and bactericidal activities dose dependently. Bp-LAAO induced platelet aggregation in platelet-rich plasma and this activity was inhibited by catalase. Bp-LAAC-cDNA of 1548 bp codified a mature protein with 516 amino acid residues corresponding to a theoretical isoelectric point and molecular weight of 6.3 and 58 kDa, respectively. Additionally, structural and phylogenetic studies identified residues under positive selection and their probable location in Elp-LAAO and other snake venom LAAOs (svLAAOs). Structural and functional investigations of these enzymes can contribute to the advancement of toxinology and to the elaboration of novel therapeutic agents. (C) 2009 Elsevier Masson SAS. All rights reserved.

Spin current in the Mobius cyclacene belts

We report on density functional theory studies of the electronic structure and magnetic properties of Mobius-[n]cyclacenes. The geometry of Mobius bands presents a modulation of bond lengths that is needed to accommodate the twist. This modulation takes the form of bond alternation defects analogous to those of solitons in polyacetylene. The ground state of all Mobius bands is a triplet, with a spin density distribution that follows the bond length modulation. A molecular dynamics simulation of the Mobius cyclacene at 300 K shows that the twist travels around the belt inducing a magnetic current. (C) 2009 Elsevier B. V. All rights reserved.

Ground-state energy of a classical artificial molecule

We study the ground-state energy of a classical artificial molecule formed by two-dimensional clusters (artificial atoms) of N/2 charged particles separated by a distance d. For the small molecules of N = 2 and 4, we obtain analytical expressions for this energy. For the larger ones, we calculate the ground-state energy using molecular dynamics simulation for N up to 128. From our numerical results, we are able to find out a function to approximate the ground-state energy of the molecules covering the range from atoms to molecules for any inter-atom distance d and for particle number from N = 8 to 128 within a difference less than one percent from the MD data.

Catalytic mechanism of inulinase from Arthrobacter sp S37

Detailed catalytic roles of the conserved Glu323, Asp460, and Glu519 of Arthrobacter sp. S37 inulinase (EnIA), a member of the glycoside hydrolase family 32, were investigated by site-directed mutagenesis and pH-dependence studies of the enzyme efficiency and homology modeling were carried out for EnIA and for D460E mutant. The enzyme efficiency (k(cat)/K-m) of the E323A and E519A mutants was significantly lower than that of the wild-type due to a substantial decrease in k(cat), but not due to variations in K-m, consistent with their putative roles as nucleophile and acid/base catalyst, respectively. The D460A mutant was totally inactive, whereas the D460E and D460N mutants were active to some extent, revealing Asp460 as a catalytic residue and demonstrating that the presence of a carboxylate group in this position is a prerequisite for catalysis. The pH-dependence studies indicated that the pK(a) of the acid/base catalyst decreased from 9.2 for the wild-type enzyme to 7.0 for the D460E mutant, implicating Asp460 as the residue that interacts with the acid/base catalyst Glu519 and elevates its pK(a). Homology modeling and molecular dynamics simulation of the wild-type enzyme and the D460E mutant shed light on the structural roles of Glu323, Asp460, and Glu519 in the catalytic activity of the enzyme. (C) 2008 Elsevier Inc. All rights reserved.

Effect of SO(2) on the Transport Properties of an Imidazolium Ionic Liquid and Its Lithium Solution

Transport coefficients have been measured as a function of the concentration of sulfur dioxide, SO(2), dissolved in 1-butyl-2,3-dimethylimidazolium bis(trifluoromethylsulfonyl)-imide, [BMMI][Tf(2)N], as well as in its lithium salt solution, Li[Tf(2)N]. The SO(2) reduces viscosity and density and increases conductivity and diffusion coefficients in both the neat [BMMI] [Tf(2)N] and the [BMMI][Tf(2)N]-Li[Tf(2)N] solution. The conductivity enhancement is not assigned to a simple viscosity effect; the weakening of ionic interactions upon SO(2) addition also plays a role. Microscopic details of the SO(2) effect were unraveled using Raman spectroscopy and molecular dynamics (MD) simulations. The Raman spectra suggest that the Li(+)-[Tf(2)N] interaction is barely affected by SO(2), and the SO(2)-[Tf(2)N] interaction is weaker than previously observed in an investigation of an ionic liquid containing the bromide anion. Transport coefficients calculated by MD simulations show the same trend as the experimental data with respect to SO(2) content. The MD simulations provide structural information on SO(2) molecules around [Tf(2)N], in particular the interaction of the sulfur atom of SO(2) with oxygen and fluorine atoms of the anion. The SO(2)-[BMMI] interaction is also important because the [BMMI] cations with above-average mobility have a larger number of nearest-neighbor SO(2) molecules.