Effect of salt on the interactions of antimicrobial peptides with zwitterionic lipid bilayers

The effect of salt on the binding of the antimicrobial peptide magainin to POPC lipid bilayers is studied by 40-50 ns molecular dynamics simulations of a POPC bilayer in the presence of different concentrations of Na+ and Cl− ions, corresponding to effective concentrations of 0, 100, 150, 200, 250 and 300 millimolar NaCl, with and without a single molecule of antimicrobial peptide magainin. Simulations without magainin showed that increasing salt concentration leads to the decrease in the area per lipid, a decrease in the head group tilt of the lipids, as well as increased order of lipid tails, in agreement with other recent simulations. Simulations with magainin show that peptide binding to the lipids is stronger at lower concentrations of salt. The peptides disorder the lipids in their immediate vicinity, but this effect is diminished as the salt concentration increases. Our studies indicate that while 50 ns simulations give information on peptide hydrogen bonding and lipid tail ordering that is insensitive to the initial peptide orientation, this run time is not sufficient to equilibrate the peptide position and orientation within the bilayer.     

pUC12-Wl: a plastnid vector for the study of cloned DNA conformation

pUC12-Wl is a new cloning vector for the study of torsion-induced structural transitions of insert DNA. It was derived from pUC12 by deleting three A + T-rich sequences which can undergo structural transitions when torsionally stressed. Transitions at these sites have low energy of activation and undefined structures. They complicate studies on transitions of DNA inserts by diverting torsional force and causing the vector to be undefined in helical and energetic terms. The new vector pUC12-Wl, from which these segments have been deleted, will facilitate studies of torsion-induced structural transitions of insert DNA.     

Probing the interactions of the solvated electron with DNA by molecular dynamics simulations: bromodeoxyuridine substituted DNA

Solvated electrons () are produced during water radiolysis and can interact with biological substrates, including DNA. To augment DNA damage, radiosensitizers such as bromo-deoxyuridine (BUdR), often referred to as an “electron affinic radiosensitizer”, are incorporated in place of isosteric thymidine. However, little is known about the primary interactions of with DNA. In the present study we addressed this problem by applying molecular modeling and molecular dynamics (MD) simulations to a system of normal (BUdR·A)-DNA and a hydrated electron, where the excess electron was modeled as a localized (H₂O)₆ anionic cluster. Our goals were to evaluate the suitability of the MD simulations for this application; to characterize the motion of around DNA (e.g., diffusion coefficients); to identify and describe configurational states of close localization to DNA; and to evaluate the structural dynamics of DNA in the presence of . The results indicate that has distinct space-preferences for forming close contacts with DNA and is more likely to interact directly with nucleotides other than BUdR. Several classes of DNA - contact sites, all within the major groove, were distinguished depending on the structure of the intermediate water layer H-bonding pattern (or its absence, i.e., a direct H-bonding of with DNA bases). Large-scale structural perturbations were identified during and after the approached the DNA from the major groove side, coupled with deeper penetration of sodium counterions in the minor groove. Figure A rare configuration showing direct interaction between the solvated electron and DNA, where (yellow) and N7(A16) are H-bonded. The close approach from the major groove side invokes deep Na⁺ (magenta) penetration into the minor DNA groove (Fig. 7a).     

A molecular-dynamics simulation study of diffusion of a single model carbonic chain on a graphite (001) surface

Molecular-dynamics simulations have been used to study the diffusion of a short single model carbonic chain on the graphite (001) surface. The calculated diffusion coefficient (D) first increases, then decreases with increasing chain length (N). This abnormal behavior is similar to polymer lateral diffusion at the solid–liquid interface. Furthermore, we have studied the relation between the mean-square gyration radius and N. Figure Log–log plot of the self-diffusion coefficient D versus the chain length N. The error bars are the standard deviation measured in three repeated simulations     

An elastic rod model to evaluate effects of ionic concentration on equilibrium configuration of DNA in salt solution

As a coarse-gained model, a super-thin elastic rod subjected to interfacial interactions is used to investigate the condensation of DNA in a multivalent salt solution. The interfacial traction between the rod and the solution environment is determined in terms of the Young–Laplace equation. Kirchhoff’s theory of elastic rod is used to analyze the equilibrium configuration of a DNA chain under the action of the interfacial traction. Two models are established to characterize the change of the interfacial traction and elastic modulus of DNA with the ionic concentration of the salt solution, respectively. The influences of the ionic concentration on the equilibrium configuration of DNA are discussed. The results show that the condensation of DNA is mainly determined by competition between the interfacial energy and elastic strain energy of the DNA itself, and the interfacial traction is one of forces that drive DNA condensation. With the change of concentration, the DNA segments will undergo a series of alteration from the original configuration to the condensed configuration, and the spiral-shape appearing in the condensed configuration of DNA is independent of the original configuration.     

On the importance of van der Waals interaction in the groove binding of DNA with ligands: restrained molecular dynamics study

Comparison of interaction energy between an oligonucleotide and a DNA-binding ligand in the minor and major groove modes was made by use of restrained molecular dynamics. Distortion in DNA was found for the major groove mode whereas less significant changes for both ligand and DNA were detected for the minor groove binding after molecular dynamics simulation. The conformation of the ligand obtained from the major groove mode resembles that computed with the ligand soaked in water. The van der Waals contact energy was found to be as significant as electrostatic energy and more important for difference in binding energy between these two binding modes. The importance of van der Waals force in groove binding was supported by computations on the complex formed by the repressor peptide fragment from the bacteriophage 434 and its operator oligonucleotide.     

Effect of external stresses on protein conformation: a computer modelling study

The increasing use of digital technologies such as mobile phones has led to major health concerns about the effects of non-ionizing pulsed radiation exposure. We believe that the health implications of exposure to radiation cannot be fully understood without establishing the molecular mechanisms of biological effects of pulsed microwaves. We aim to establish methods for studying the molecular mechanisms of protein structural and energetic changes occurring due to external stresses related to non-ionizing radiation by using a combination of experimental and theoretical approaches. In this paper, we present the results from our fully atomistic simulation study of chemical and thermal stress response of a prototype protein, insulin. We performed a series of molecular dynamics simulations of insulin in solution under equilibrium conditions, under chemical stress (imitated by reducing the disulfide bonds in the protein molecule), and under short-lived thermal stress (imitated by increasing simulation temperature for up to 2 ns). The resultant protein conformational behaviour was analysed for various properties with the aim of establishing analysis routines for classification of protein unfolding pathways and associated molecular mechanisms.Submitted as a record of the 2002 Australian Biophysical Society meeting     

The Interaction Between the Novel Intercalator Diethidium Cation and B-Form DNA: a Theoretical Study

Abstract

This research is an effort to further understand the physicochemical interaction between the novel drug molecule diethidium (2,7-diamino 9-[2,7 diamino 10-nN- phenanthridium] 10- nN- phenanthridium) and its biological receptor DNA. The ultimate goal is the elucidation of this novel class of drugs as potential pharmaceutical agents. Understanding the physico- chemical properties of this drug as well as the mechanism by which it interacts with DNA should ultimately allow the rational design of novel anti-cancer or anti-viral drugs.

A novel binding structure for the diethidium cation to B-form DNA is herein described. Molecular modeling on the complex formed between diethidium and a dodecamer of double-stranded B-form DNA, CGCGAATTCGCG, has shown that this complex is indeed fully capable of participating in the formation of a stable intercalation site. It was expected that diethidium would have a mechanism of intercalation significantly different from other classical intercalators because a) Its structure, that of two perpendicular planes, each known to have excellent intercalation properties, is novel b) The linker region length is zero c) The tilt between the two planes of the drug matches the geometry of the space available to this drug in the major groove.

We have studied the complex formed when diethidium enters the central site of the B-DNA dodecamer through the major groove. The complex forms several classes of intercalation structures, which are all stable and vary from “partially intercalate” to “fully intercalated”. Multiple minimizations show the drug to be very mobile within the intercalation site. Further, some structures show organization and concomitant stiffening of the DNA above the intercalation site, with a disorganization and disruption of the regular B-DNA structure immediately below the intercalation site. This particular phenomena may be expected to lead to significantly different physicochemical properties for the diethidium complex with respect to other known intercalators, because this sort of vectorial difference in structure above and below the site of intercalation is unknown in existing intercalators, as far as the authors are aware. In addition, we expect the mechanism of interaction between drug and DNA to be described by “direct ligand transfer”, wherein the drug is transferred from duplex DNA to duplex DNA without re-entering the solvent.¹

This work is the first instance known to the authors of a novel drug entity that was deduced solely by mathematical reasoning ² and described subsequently by computational methods. Evidence that diethidium should interact with its target site DNA differently from other known intercalators is strong.     

A comparative molecular dynamics study on thermostability of human and chicken prion proteins

To compare the thermostabilities of human and chicken normal cellular prion proteins (HuPrP(C) and CkPrP(C)), molecular dynamics (MD) simulations were performed for both proteins at an ensemble level (10 parallel simulations at 400 K and 5 parallel simulations at 300 K as a control). It is found that the thermostability of HuPrP(C) is comparable with that of CkPrP(C), which implicates that the non-occurrence of prion diseases in non-mammals cannot be completely attributed to the thermodynamic properties of non-mammalian PrP(C).     

Molecular dynamics study of methane in water: diffusion and structure

In this study, theoretical analysis on the geometries and electronic properties of various conjugated oligomers based on thiophene (Th) or bicyclic non-classical Th units is reported. The dihedral angle, bond length, bond-length alternation, bond critical point (BCP) properties, nucleus-independent chemical shift (NICS) and Wiberg bond index (WBIs) are analysed and correlated with conduction properties. The changes of bond length, BCP properties, NICS and WBIs all show that the conjugational degree is increased systematically with main chain extension. As a result, the highest occupied molecular orbital–lowest unoccupied molecular orbital energy separation (E _g) is decreased upon chain elongation. The E _g of oligomers based on bicyclic non-classical Th unit is much lower than that of Th-based oligomers due to the narrower E _g of bicyclic non-classical Ths compared to Th, which indicates that the narrow E _g of the bicyclic non-classical Ths can be carried over to their polymers by using them as building blocks for the polymers. The band structures and density of states analysis show that the four polymers all have small band gap ( < 0.9 eV), wide highest occupied bandwidth and relatively small effective mass of hole, which indicate that those proposed polymers may be potential conductors.     

Wetting transition of nanodroplets of water on textured surfaces: a molecular dynamics study

The crossover behaviour of water droplet's state from the Wenzel state to the Cassie state with varying pillar height and surface fraction is examined critically using molecular dynamics. We report the effect of the system size on the wetting behaviour of water droplets by examining the contact angle for both regimes. We observe that when the droplet size is comparable to the pillar dimension, the contact angle of droplets fluctuates with increasing droplet size because of the contact line pinning, which is more pronounced in the Wenzel regime. We further demonstrate the phantom-wall method to evaluate free energy of intermediate wetting states.     

Investigations of peptide hydration using NMR and molecular dynamics simulations: A study of effects of water on the conformation and dynamics of antamanide

Summary The influence of water binding on the conformational dynamics of the cyclic decapeptide antamanide dissolved in the model lipophilic environment chloroform is investigated by NMR relaxation measurements. The water-peptide complex has a lifetime of 35 s at 250 K, which is longer than typical lifetimes of water-peptide complexes reported in aqueous solution. In addition, there is a rapid intracomplex mobility that probably involves librational motions of the bound water or water molecules hopping between different binding sites. Water binding restricts the flexibility of antamanide. The experimental findings are compared with GROMOS molecular dynamics simulations of antamanide with up to eight bound water molecules. Within the simulation time of 600 ps, no water molecule leaves the complex. Additionally, the simulations show a reduced flexibility for the complex in comparison with uncomplexed antamanide. Thus, there is a qualitative agreement between the experimental NMR results and the computer simulations.     

The behaviour of beta-carotene in the phosphatidylcholine bilayer as revealed by a molecular simulation study

A molecular dynamics (MD) simulation of the fully hydrated bilayer made of 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) and containing beta-carotene (beta-Car) molecules was carried out as a complementary approach to experimental techniques to investigate the orientation of beta-Car in the lipid membrane as well as its influence on the bilayer properties. The bilayer reached thermal equilibrium after 1200 ps of MD simulation and the productive run was carried out for 2800 ps. The results indicate that the carotene rings are located in the region occupied by the carbonyl groups of the POPC gamma-chain with no trace of penetration towards the centre of the bilayer. Carotene exhibits an ordering effect on both the beta- and the gamma-chain. While the fully saturated gamma-chain is affected evenly along, the order of the mono-unsaturated beta-chain is modified mainly below the double bond. In general, a high value of the order parameter and the chain tilt in the range from 11.4 degrees to 26.7 degrees were observed for the beta-Car molecules. However, for chain segment adjacent to methyl groups the value of the order parameter is low and the tilt angle is close to 75 degrees . Moreover, the probability of trans conformation being generally close to 1.0 along the beta-Car chain is reduced for these segments. Our MD simulation study suggests two pools of the preferential orientation of beta-Car: a slightly bent structure corresponding to a small chain tilt angle and a rather stretched structure that corresponds to a higher chain tilt. The results are discussed in the light of experimental findings.     

The glass transition temperatures of amorphous trehalose-water mixtures and the mobility of water: an experimental and in silico study

Isothermal-isobaric molecular dynamics simulations are used to calculate the specific volume of models of trehalose and three amorphous trehalose-water mixtures (2.9%, 4.5% and 5.3% (w/w) water, respectively) as a function of temperature. Plots of specific volume versus temperature exhibit a characteristic change in slope when the amorphous systems change from the glassy to the rubbery state and the intersection of the two regression lines provides an estimate of the glass transition temperature T(g). A comparison of the calculated and experimental T(g) values, as obtained from differential scanning calorimetry, shows that despite the predicted values being systematically higher (about 21-26K), the trend and the incremental differences between the T(g) values have been computed correctly: T(g)(5.3%(w/w))相似文献   

Effect of gold nanoparticle on stability of the DNA molecule: A study of molecular dynamics simulation

An understanding of the mechanism of DNA interactions with gold nanoparticles is useful in today medicine applications. We have performed a molecular dynamics simulation on a B-DNA duplex (CCTCAGGCCTCC) in the vicinity of a gold nanoparticle with a truncated octahedron structure composed of 201 gold atoms (diameter ～1.8 nm) to investigate gold nanoparticle (GNP) effects on the stability of DNA. During simulation, the nanoparticle is closed to DNA and phosphate groups direct the particles into the major grooves of the DNA molecule. Because of peeling and untwisting states that are occur at end of DNA, the nucleotide base lies flat on the surface of GNP. The configuration entropy is estimated using the covariance matrix of atom-positional fluctuations for different bases. The results show that when a gold nanoparticle has interaction with DNA, entropy increases. The results of conformational energy and the hydrogen bond numbers for DNA indicated that DNA becomes unstable in the vicinity of a gold nanoparticle. The radial distribution function was calculated for water hydrogen–phosphate oxygen pairs. Almost for all nucleotide, the presence of a nanoparticle around DNA caused water molecules to be released from the DNA duplex and cations were close to the DNA.     

Molecular dynamics simulations of a cyclic-DP-240 amylose fragment in a periodic cell: Glass transition temperature and water diffusion

Molecular dynamics simulations using AMB06C, an in-house carbohydrate force field, (NPT ensembles, 1 atm) were carried out on a periodic cell that contained a cyclic 240 glucose residue amylose fragment (c-DP-240) and TIP3P water molecules. Molecular conformation and movement of the amylose fragment and water molecules at different temperatures were examined. The periodic cell volume, density, and potential energy were determined at temperatures above and below the glass transition temperature (T_g) in 25 K increments_. The amorphous cell is constructed through successive dynamic equilibration steps at temperatures above the assumed T_g value and the temperature successively lowered until several temperature points were obtained below T_g. Molecular dynamics simulations were continued for at least 500 ps or until the volume drift stopped and remained constant for several hundred picoseconds. The T_g values were found by noting the discontinuity in slope of the volume (V), potential energy (PE), or density (ρ) versus 1/T. The changes in flexibility and motion of the amylose chain as well as differences in self diffusion coefficients of water molecules are described. The final average T_g value found (316 K) is in agreement with experimental values, i.e. 320 K.     

Capecitabine as a minor groove binder of DNA: molecular docking,molecular dynamics,and multi-spectroscopic studies

The interaction mechanism and binding mode of capecitabine with ctDNA was extensively investigated using docking and molecular dynamics simulations, fluorescence and circular dichroism (CD) spectroscopy, DNA thermal denaturation studies, and viscosity measurements. The possible binding mode and acting forces on the combination between capecitabine and DNA had been predicted through molecular simulation. Results indicated that capecitabine could relatively locate stably in the G-C base-pairs-rich DNA minor groove by hydrogen bond and several weaker nonbonding forces. Fluorescence spectroscopy and fluorescence lifetime measurements confirmed that the quenching was static caused by ground state complex formation. This phenomenon indicated the formation of a complex between capecitabine and ctDNA. Fluorescence data showed that the binding constants of the complex were approximately 2 × 10⁴ M^?1. Calculated thermodynamic parameters suggested that hydrogen bond was the main force during binding, which were consistent with theoretical results. Moreover, CD spectroscopy, DNA melting studies, and viscosity measurements corroborated a groove binding mode of capecitabine with ctDNA. This binding had no effect on B-DNA conformation.     

Effect of a salt water intrusion on a freshwater chironomidae community: A paleolimnological study

Examination of two 90 cm sediment cores from Front Lake, a shallow (1 m) eutrophic habitat, revealed a chironomid community consisting of at least 24 species prior to the destruction of the fauna by an intrusion of salt water resulting from the high tides associated with the Saxby Gale of 1869. With the withdrawal of the salt water and a return to freshwater conditions, 23 species recolonized the lake. Of these, 22 were also found in the preintrusion community. The three species not common to both freshwater periods constituted only 2.3% of the head capsules recovered. Consequently, the chironomid community established in Front Lake both before and after the salt water intrusion had essentially identical species compositions, relative abundances of each species, and species diversity indices.