The dynamic action mechanism of small cationic antimicrobial peptides.

Antimicrobial peptides form part of the immune system as protection against the action of external pathogens. The differences that exist between mammalian and microbial cell membrane architectures are key aspects of the ability of these peptides to discriminate between pathogens and host cells. Given that the pathogen membrane is the non-specific target of these cationic peptides, different molecular mechanisms have been suggested to describe the rules that permit them to distinguish between pathogens and mammalian cells. In this context, and setting aside the old fashion idea that cationic peptides act through one mechanism alone, this work will provide insight into the molecular action mechanism of small antimicrobial peptides, based on molecular dynamics simulations of phospholipid bilayers that mimic different cell membrane architectures. After measuring different properties of these lipid bilayers, in the absence and presence of peptides, a four-step action mechanism was suggested on the basis of the formation of phospholipid rafts induced by the presence of these cationic peptides. Thus, this work shows how differences in the bending modulus (k(b)) of these lipid rafts and differences in the free energy profiles (ΔG(z)) associated with the insertion of these peptides into these lipid rafts are key aspects for explaining the action mechanism of these cationic peptides at the molecular level.

Production of gas diffusion layers with cotton fibers for their use in fuel cells.

The gas diffusion layer (GDL) is one of the most important parts of a proton exchange membrane fuel cell, that plays a key role transporting the current to the collector plates, distributing the reactant gases to the catalyst surface, and evacuating heat and water that is generated during the redox reactions inside the fuel cell. Speaking in terms of production cost, the GDL represents up to 45% of the total cost of the membrane electrode assembling (MEA). However, and despite its crucial role in a fuel cell, until recent years, the GDLs have not been studied with the same intensity as other MEA components, such as the catalyst or the proton exchange membrane. In this work, we present the production process, at laboratory scale, of a low cost GDL, using a non-woven paper-making process. A relevant aspect of this GDL is that up to 40% of their composition is natural cotton, despite which they present good electrical and thermal conductivity, high porosity, good pore morphology, high hydrophobicity as well as gas permeability. Furthermore, when the GDL with its optimum cotton content was tested in a single open cathode fuel cell, a good performance was obtained, which makes this GDL a promising candidate for its use in fuel cells.

Influence of Lipid Composition on the Insertion Process of Glyphosate into Membranes: A Thermodynamic Study.

In this work, molecular dynamics simulations were applied to investigate the influence of lipid composition of the model membrane on the insertion of glyphosate (in its charged state, GLYP2-). The profiles of free energy, entropy and enthalpy were obtained through umbrella sampling calculations, for lipid bilayers composed by only 1,2-dipalmitoyl-sn-glycerol-3-phosphocholine (DPPC), only 1,2-dipalmitoyl-sn-glycerol-3-phosphoserine (DPPS) or a symmetric binary mixture of DPPC and DPPS. In general, the location, the values of minima and maxima of the free energy, and the trend of free energy profiles are influenced by the lipid composition of the lipid bilayer. The driving force in the glyphosate insertion process depends on the lipid composition of the membrane model. If the lipid bilayer is composed solely of DPPS or DPPC, GLYP2- insertion is driven by a favorable enthalpic change. However, if the membrane is composed of a mixture of both lipids, this process is driven by a favorable entropic change. In the lipid bilayer containing DPPS, the glyphosate was found to penetrate hydrated and coordinated with Na+ ions, in contrast to the pure zwitterionic lipid bilayer which penetrated only hydrated. This effect is independent of the concentration of sodium ions present in the bulk solution.

Physicochemical study of the acetonitrile insertion into polypyrrole films.

A study by molecular dynamics (MD) simulation of the acetonitrile diffusion into a polypyrrole film was carried out with atomic detail in a 0.1N lithium perchlorate solution. From the simulated trajectories, the acetonitrile behavior was estimated from bulk solution to the interior of the polypyrrole film, across the polypyrrole/solution interface, for a neutral (reduced) and charged (oxidized) state of the polymer. Among other properties, the translational diffusion coefficient and rotational relaxation time of the acetonitrile were calculated, where a diminution in the translational diffusion coefficient was measured in the interior of the polypyrrole matrix compared to bulk, independently of the oxidation state of the polymer, in contrast with the behavior of the rotational relaxation time that increases from bulk to the interior of the polymer for both oxidation states. In addition, the difference of free energy DeltaG associated to the acetonitrile penetration into the polymer was calculated. From the results, it was evidenced that the scarce affinity of acetonitrile to diffuse into the polymer in its reduced state is related with the positive uniform difference of free energy DeltaG approximately 20 kJ/mol, while in the oxidized state, an important free energy barrier of DeltaG approximately 10 kJ/mol has to pass trough for reaching stable sites inside the polymer with values of DeltaG up to -10 kJ/mol.

Antibacterial Effect of Chitosan-Gold Nanoparticles and Computational Modeling of the Interaction between Chitosan and a Lipid Bilayer Model.

Pathogenic bacteria have the ability to develop antibiotic resistance mechanisms. Their action consists mainly in the production of bacterial enzymes that inactivate antibiotics or the appearance of modifications that prevent the arrival of the drug at the target point or the alteration of the target point itself, becoming a growing problem for health systems. Chitosan-gold nanoparticles (Cs-AuNPs) have been shown as effective bactericidal materials avoiding damage to human cells. In this work, Cs-AuNPs were synthesized using chitosan as the reducing agent, and a systematic analysis of the influence of the synthesis parameters on the size and zeta potential of the Cs-AuNPs and their UV-vis spectra was carried out. We used a simulation model to characterize the interaction of chitosan with bacterial membranes, using a symmetric charged bilayer and two different chitosan models with different degrees of the chitosan amine protonation as a function of pH, with the aim to elucidate the antibacterial mechanism involving the cell wall disruption. The Cs-AuNP antibacterial activity was evaluated to check the simulation model.

Study of the benzocaine transfer from aqueous solution to the interior of a biological membrane.

The precise molecular mechanism of general anesthetics remains unknown. It is therefore important to understand where molecules with anesthetic properties localize within biological membranes. We have determined the free energy profile of a benzocaine molecule (BZC) across a biological membrane using molecular dynamics simulation. We use an asymmetric phospholipid bilayer with DPPS in one leaflet of a DPPC bilayer (Lopez Cascales et al. J. Phys. Chem. B 2006, 110, 2358-2363) to model a biological bilayer. From the free energy profile, we predict the zone of actuation of a benzocaine is located in the hydrocarbon region or at the end of the lipid head, depending of the presence of charged lipids (DPPS) in the leaflet. We observe a moderate increase in the disorder of the membrane and in particular an increase in the disorder of DPPS. Static and dynamic physicochemical properties of the benzocaine, such as its dipole orientation, translational diffusion coefficient, and rotational relaxation time were measured.

Molecular dynamics simulations of glyphosate in a DPPC lipid bilayer.

Extensive molecular dynamics simulations have been performed to study the effect of glyphosate (in their neutral and charged forms, GLYP and GLYP2-, respectively) on fully hydrated DiPalmitoylPhosphatidylCholine (DPPC) lipid bilayer. First, we calculated the free energy profile (using the Umbrella Sampling technique) for both states of charge of glyphosate. The minimum value for the free energy for GLYP is ∼-60 kJ mol-1 located at z = ±1.7 nm (from the lipid bilayer center), and there is almost no maximum at the center of the lipid bilayer. By contrast, the minimum for GLYP2- is ∼-35 kJ mol-1 located at z = ±â€¯1.4 nm (from the lipid bilayer center), and the maximum reaches ∼35 kJ mol-1 at the center of the lipid bilayer. Then, different lipid bilayer properties were analyzed for different glyphosate:lipid (G:L) ratios. The mean area per lipid was slightly affected, increasing only 5% (in the presence of glyphosate at high concentrations), which is in agreement with the slight decrease in deuterium order parameters. As for the thickness of the bilayer, it is observed that the state of charge produces opposite effects. On one hand, the neutral state produces an increase in the thickness of the lipid bilayer; on the other, the charged form produces a decrease in the thickness, which not depend linearly on the G:L ratios, either. The orientation of the DPPC head groups is practically unaffected throughout the range of the G:L ratios studied. Finally, the mobility of the lipids of the bilayer is strongly affected by the presence of glyphosate, considerably increasing its lateral diffusion coefficient noteworthy (one order of magnitude), with increasing G:L ratio.

Model of an asymmetric DPPC/DPPS membrane: effect of asymmetry on the lipid properties. A molecular dynamics simulation study.

The study of asymmetric lipid bilayers is of a crucial importance due to the great number of biological process in which they are involved such as exocytosis, intracellular fusion processes, phospholipid-protein interactions, and signal transduction pathway. In addition, the loss of this asymmetry is a hallmark of the early stages of apoptosis. In this regard, a model of an asymmetric lipid bilayer composed of DPPC and DPPS was simulated by molecular dynamics simulation. Thus, the asymmetric membrane was modeled by 264 lipids, of which 48 corresponded to DPPS- randomly distributed in the same leaflet with 96 DPPC. In the other leaflet, 120 DPPC were placed without DPPS-. Due to the presence of a net charge of -1 for the DPPS- in physiological conditions, 48 Na+ were introduced into the system to balance the charge. To ascertain whether the presence of the DPPS- in only one of the two leaflets perturbs the properties of the DPPC in the other leaflet composed only of DPPC, different properties were studied, such as the atomic density of the different components across the membrane, the electrostatic potential across the membrane, the translational diffusion of DPPC and DPPS, the deuterium order parameters, lipid hydration, and lipid-lipid charge bridges. Thus, we obtained that certain properties such as the surface area lipid molecule, lipid head orientation, order parameter, translational diffusion coefficient, or lipid hydration of DPPC in the leaflet without DPPS remain unperturbed by the presence of DPPS in the other leaflet, compared with a DPPC bilayer. On the other hand, in the leaflet containing DPPS, some of the DPPC properties were strongly affected by the presence of DPPS such as the order parameter or electrostatic potential.

Effect of lithium and sodium ions on a charged membrane of dipalmitoylphosphatidylserine: a study by molecular dynamics simulation.

We describe a series of molecular dynamics simulations performed on a model of charged lipid bilayer (dipalmitoylphosphatidylserine) and water, in presence of sodium and lithium ions, with an atomic detail. The structure of the lipid membranes was strongly affected by the presence of lithium, as manifested by the observation of a transition from a disordered to a gel state. Concerning the mechanism of such a transition, it was associated to the dehydration that we detected in the lipid-water interface in the presence of lithium. This dehydration introduced an increase in the lipid-lipid interactions, and as a consequence, a diminution of the disorder of the membrane. When both types of ions are present in the aqueous phase, lithium shown a special affinity for the lipid membrane displacing almost all the sodium ions toward the middle of the water layer. As a result, we observed remarkable differences in the atom and electric field distributions across the lipid membrane. Concerning the diffusion and orientation of water molecules across the lipid-water interface, we also observed a strong dependency of the type ion. On the other hand, the mobility and the hydration shell of lithium and sodium ions are strongly perturbed by the presence of the charged lipid bilayer. The lipid layer was responsible for a dehydration of the ions compared to bulk water. This dehydration was compensated by an increase of coordination number of the ions with the lipid oxygens. Also, the residence times of water in the first hydration shell of lithium and sodium ions are perturbed by the presence of the lipid membrane.

Perchlorate interchange during the redox process of PPy/PVS films in an acetonitrile medium. a voltammetric and EDX study.

Polypyrrole/poly(vinyl sulfonate) (PPy/PVS) films in acetonitrile containing 0.1 M LiClO4 were studied by cyclic voltammetry. Consecutive voltammograms pointed to a continuous increase in the charge involved in the process, suggesting a rise in the number of the electroactive participants involved in the redox process. However, voltammograms obtained for the PPy/ClO4 films in analogous conditions pointed to a steady-state behavior from the very early cycles. Theoretical studies based on the Nernst and Butler-Volmer equations indicated that perchlorate ions are involved during the oxidation/reduction process of the PPy/PVS films when the steady state is reached. This result was confirmed by "ex situ" energy-dispersive X-ray analysis of the films. In this regard, the electrochemical behavior of PPy/PVS polymers was similar to that of PPy/ClO4 films when a high number of cycles were carried out. The exchange of ClO4- during the redox reaction of the PPy/PVS films made it necessary to incorporate Li+ cations inside the polymer during the initial voltammetric cycles to compensate for the negative charges of PVS polyanions. Li+ cations are mainly stabilized inside the polymer by the ion pairs formed with the sulfonated groups of the PVS. An increase and shift of the voltammetric cycles indicated a restructuring of the polymeric chains with consecutive scans.

Molecular dynamics simulation of a dye molecule in the interior of a bilayer: 1,6-diphenyl-1,3,5-hexatriene in dipalmitoylphosphatidylcholine.

A molecular dynamics simulation was carried out for a dipalmitoylphosphatidylcholine (DPPC) membrane in its liquid crystalline state containing different concentrations of the dye molecule 1,6-diphenyl-1,3,5-hexatriene (DPH). From a numerical analysis of the trajectories, we obtained information concerning structural changes of the membrane due to the presence of the probe and some hydrodynamic information concerning the probe itself. The hydrodynamic properties regarding dye molecules that have been reported in this article are: rotational and translational diffusion coefficient and relaxation times. From this analysis, we estimated a range of values of 0.6-0.9 cP for the micro-viscosity in the mid-membrane. These simulations also afforded us some information regarding structural changes in the membrane as a consequence of the presence of the fluorescent dyes at different concentrations. Thus, the disorder inside the membrane, the surface area per lipid and thickness of the membrane were also investigated.

Bead-model calculation of scattering diagrams: Brownian dynamics study of flexibility in immunoglobulin IgG1.

We consider the calculation of the angular dependence of the scattering intensities for models composed of spherical elements of arbitrary size, in which some of the spheres may have a size close to that of the whole particle. This minimizes the number of spheres and makes it possible to use the same bead models for the prediction of scattering diagrams and hydrodynamics properties. A simple expression is employed for the scattering intensity. Particularly, we discuss the validity of some versions of the Debye scattering equation for bead models. The method could be used for any macromolecule either rigid or with conformational variability. The application of the method to a simple model for IgG1 shows that influence of flexibility in the scattering curve is stronger by the end of the first decade of decay and also shows that the range of linearity in the Guinier region is not a good method to characterize flexibility between arms, because it is very likely that experimental errors will hide the small differences between the extreme cases appreciated in our calculations.

Effect of the interfacial tension and ionic strength on the thermodynamic barrier associated to the benzocaine insertion into a cell membrane.

The insertion of local anaesthetics into a cell membrane is a key aspect for explaining their activity at a molecular level. It has been described how the potency and response time of local anaesthetics is improved (for clinical applications) when they are dissolved in a solution of sodium bicarbonate. With the aim of gaining insight into the physico-chemical principles that govern the action mechanism of these drugs at a molecular level, simulations of benzocaine in binary lipid bilayers formed by DPPC/DPPS were carried out for different ionic strengths of the aqueous solution. From these molecular dynamic simulations, we observed how the thermodynamic barrier associated with benzocaine insertion into the lipid bilayers diminished exponentially as the fraction of DPPS in the bilayer increased, especially when the ionic strength of the aqueous solution increased. In line with these results, we also observed how this thermodynamic barrier diminished exponentially with the phospholipid/water interfacial tension.

Binding of glutamate to the umami receptor.

The umami taste receptor is a heterodimer composed of two members of the T1R taste receptor family: T1R1 and T1R3. It detects glutamate in humans, and is a more general amino acid detector in other species. We have constructed homology models of the ligand binding domains of the human umami receptor (based on crystallographic structures of the metabotropic glutamate receptor of the central nervous system). We have carried out molecular dynamics simulations of the ligand binding domains, and we find that the likely conformation is that T1R1 receptor protein exists in the closed conformation, and T1R3 receptor in the open conformation in the heterodimer. Further, we have identified the important binding interactions and have made an estimate of the relative free energies associated with the two glutamate binding sites.

Study of the effect of Na+ and Ca2+ ion concentration on the structure of an asymmetric DPPC/DPPC + DPPS lipid bilayer by molecular dynamics simulation.

A molecular dynamics simulation study of the steady and dynamic properties of an asymmetric phospholipid bilayer was carried out in the presence of sodium or calcium ions. The asymmetric lipid bilayer was seen to resemble a cellular membrane of an eukaryotic cell, which was modeled by dipalmitoylphosphatidylcholine (DPPC) and dipalmitoylphosphatidylserine (DPPS), placing the DPPS in one of the two leaflets of the lipid bilayer. From a numerical analysis of the simulated trajectories, information was obtained with atomic resolution for both membrane leaflet concerning the effect of bilayer asymmetry on different properties of the lipid/water interface, such as the translational diffusion coefficient and rotational relaxation time of the water molecules, lipid hydration, and residence time of water around different lipid atoms. In addition, information related to lipid conformation, and lipid-lipid interactions was also analyzed.

Phase transition of a DPPC bilayer induced by an external surface pressure: from bilayer to monolayer behavior. a molecular dynamics simulation study.

Understanding the lipid phase transition of lipid bilayers is of great interest from biophysical, physicochemical, and technological points of view. With the aim of elucidating the structural changes that take place in a DPPC phospholipid bilayer induced by an external isotropic surface pressure, five computer simulations were carried out in a range from 0.1 to 40 mN/m. Molecular dynamics simulations provided insight into the structural changes that took place in the lipid structure. It was seen that low pressures ranging from 0.1 to 1 mN/m had hardly any effect on the structure, electrical properties, or hydration of the lipid bilayer. However, for pressures above 40 mN/m, there was a sharp change in the lipid-lipid interactions, hydrocarbon lipid fluidity, and electrostatic potential, corresponding to the mesomorphic transition from a liquid crystalline state (L(alpha)) to its gel state (P'(beta)). The head lipid orientation remained almost unaltered, parallel to the lipid layer, as the surface pressure was increased, although a noticeable change in its angular distribution function was evident with the phase transition.

Molecular dynamic simulation of the hydration and diffusion of chloride ions from bulk water to polypyrrole matrix.

A molecular dynamic simulation of wet polypyrrole film was carried out, in both oxidized and reduced state. The system was modeled by two layers of polypyrrole, water and chloride ions (as counterions required for charge balance in the oxidized state) in atomic detail to provide an insight into some dynamic and steady properties of the system. Our simulations pointed to a swelling of the polymer matrix after oxidation due to electrostatic repulsions between charged sites of the oxidized polypyrrole, followed by penetration of the polypyrrole by counterions to maintain the electroneutrality of the system. Associated with this penetration of counterions toward the core of the oxidized polypyrrole, dehydration of the counterions was observed. This dehydration was compensated (in part) by a strong coordination with the charged sites of the polymer. The remaining hydrophobicity inside the polymer also contributed to the dehydration of these counterions. The translational diffusion coefficient of chloride ions was also calculated at different positions of the polypyrrole/water interface, from bulk water to the inner polymer matrix. A value of 4.1 x 10(-5) cm(2) s(-1) was measured in the bulk water compared to 5 x 10(-7) cm(2) s(-1) inside the polymer, representing a diminution of two orders of magnitude for the translational diffusion coefficient from bulk water to the core of a oxidized polypyrrole matrix. These results were in good agreement with experimental data.

HYDRO: a computer program for the prediction of hydrodynamic properties of macromolecules.

HYDRO is a program for the calculation of sedimentation and diffusion coefficients, rotational relaxation times, and intrinsic viscosities of rigid macromolecules of arbitrary shape that are represented by bead models. Actually, HYDRO contains various FORTRAN callable subroutines that can be linked to the user's own programs to account for variability of shape or flexibility. Some hints are given for the use of HYDRO in various situations.