Influence of ions on genome packaging and ejection: a molecular dynamics study

We, theoretically, investigate the effect of ions on the packing and ejection dynamics of flexible and semiflexible polymers from spherical viral capsids. We find that when the polymer charge is less screened, or the Debye length increases (corresponding to a buffer with low concentration of a monovalent salt, such as Na(+)), the packing becomes more difficult and it may stop midway. Ejection, instead, proceeds more easily if the electrostatic screening is small. On the other hand, more screening (corresponding, for example, to the addition of divalent ions such as Mg(2+)) results in easier packing and slower ejection. We interpret this as resulting from electrostatic forces among the various polymer sections, which can be tuned with the type of salt present in the solution. We also discuss how the DNA structure inside the capsid changes due to screened electrostatic interactions.     

Atomic detail investigation of the structure and dynamics of DNA.RNA hybrids: a molecular dynamics study

DNA.RNA hybrid duplexes are biologically important molecules and are shown to have potential therapeutic properties. To investigate the relationship between structures, energetics, solvation and RNase H activity of hybrid duplexes in comparison with pure DNA and RNA duplexes, a molecular dynamics study using the CHARMM27 force field was undertaken. The structural properties of all four nucleic acids considered are in very good agreement with the experimental data. The backbone dihedral angles and the puckering of the (deoxy)ribose indicate that the purine rich strands retain their A-/B-like properties but the pyrimidine rich DNA strand undergoes A-B conformational transitions. The minor groove widths of the hybrid structures are narrower than those in the RNA duplex, a requirement for RNase H binding. In addition, sampling of noncanonical phosphodiester backbone dihedrals by the DNA strands, differential solvation properties and helical properties, most notably rise, are suggested to contribute to hybrids being RNase H substrates. Differential RNase H activity toward hybrids containing purine versus pyrimidine rich RNA strands is suggested to be due to sampling of values of the phosphodiester backbone dihedrals in the DNA strands. Notably, the present results indicate that hybrids have decreased flexibility as compared to RNA, in contrast to previous reports.     

Structural properties of polystyrene oligomers in different environments: a molecular dynamics study

We have performed detailed molecular dynamics simulations to investigate the effects of solvation and confinement on the structure of polystyrene (PS) oligomers in four different environments, melt, concentrated solution, dilute solution and confined concentrated solution at 450 K and 1 bar, respectively. Local packing of the monomers and the solvent (toluene, good solvent) molecules were monitored by means of radial distribution functions (RDFs). We have also investigated bond, angle, and dihedral angle distributions of the monomers. End-to-end distances, radii of gyration and persistence lengths were calculated to characterize the static properties. The chain in the dilute solution was found to exhibit more stretched conformations. Dilution effect of the solvent was observed in the RDFs between the monomer centers. Only slight conformational changes in the polymers were observed by solvation. The effect of confinement was mainly seen in the density profiles, which showed an oscillatory behavior in the confined system.     

Effect of ions on a dipalmitoyl phosphatidylcholine bilayer. a molecular dynamics simulation study

The effect of physiological concentrations of different chlorides on the structure of a dipalmitoyl phosphatidylcholine (DPPC) bilayer has been investigated through atomistic molecular dynamics simulations. These calculations provide support to the concept that Li+, Na+, Ca2+, Mg2+, Sr2+, Ba2+, and Ac3+, but not K+, bind to the lipid-head oxygens. Ion binding exhibits an influence on lipid order, area per lipid, orientation of the lipid head dipole, the charge distribution in the system, and therefore the electrostatic potential across the head-group region of the bilayer. These structural effects are sensitive to the specific characteristics of each cation, i.e., radius, charge, and coordination properties. These results provide evidence aimed at shedding some light into the apparent contradictions among different studies reported recently regarding the ordering effect of ions on zwitterionic phosphatidylcholine lipid bilayers.     

Role of ions on structure and stability of a synthetic gramicidin ion channel in solution. A molecular dynamics study

We performed a molecular dynamics (MD) simulation to the investigate structure and stability of a synthetic gramicidin-like peptide in solution with and without ions. The starting structures of the MD simulations were taken from two recently solved NMR structures of this peptide in isotropic solution, which forms stable monomers or dimers in the presence or absence of ions, respectively. The monomeric structure is channel-like and is assumed to be stabilized by the presence of two Cs(+) ions bound in the channel, each one close to one channel entrance. In our MD simulations, we observed how the Cs(+) ions bind in the channel formed by the monomeric gramicidin-like peptide using implicit solvent and explicit ions with a concentration of 2 M. MD simulations were performed with and without explicit ions but with an implicit solvent model defined by the generalized Born approximation, which was used to mimic the dielectric properties of the solvent and to speed up the computations.     

Packing structures of single-stranded DNA and double-stranded DNA thiolates on Au(111): a molecular simulation study

The packing structures of single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA) thiolates on implicit gold surfaces were studied in explicit aqueous solutions of 1M NaCl using molecular dynamics simulations. The simulations were based on individual DNA chains placed in hexagonal simulation boxes of different sizes, representing various packing densities. The total potential energy per DNA chain was compared. The optimal packing structures were determined based on the minimal potential energy within the limits of the conditions that were evaluated in this study. The optimal packing density of ssDNA was found to be 0.19 DNA chains/nm(2), which is consistent with that determined experimentally. Furthermore, the optimal packing density of dsDNA was shown to be approximately 58% of the packing density for ssDNA, indicating that the packing of ssDNA should be approximately 58% of its optimal packing in order to achieve the best hybridization.     

Evidence of impurities in thiolated single-stranded DNA oligomers and their effect on DNA self-assembly on gold

The diversity of techniques used in the synthesis, treatment, and purification of the single-stranded DNA oligomers containing a thiol anchor group (SH-ssDNA) has led to a significant variation in the purity of commercially available SH-ssDNA. In this work, we use X-ray photoelectron spectroscopy (XPS) and time-of-flight secondary ion mass spectrometry (ToF-SIMS) to study how the impurities present in commercially synthesized SH-ssDNA oligomers affected the structure of the resulting DNA films on Au. XPS results indicate that two of the purchased SH-ssDNA oligomers contain excess carbon and sulfur. The molecular fragmentation patterns obtained with ToF-SIMS were used to determine the identity of several contaminants in the DNA films, including poly(dimethylsiloxane) (PDMS), lipid molecules, and sulfur-containing molecules. In particular, the ToF-SIMS results determined that the excess sulfur detected by XPS was due to the presence of dithiothreitol, a reductant often used to cleave disulfide precursors. Furthermore, we found that the SH-ssDNA self-assembly process is affected by the presence of these contaminants. When relatively pure SH-ssDNA is used to prepare the DNA films, the P, N, O, and C atomic percentages were observed by XPS to increase over a 24-h time period. In contrast, surfaces prepared using SH-ssDNA containing higher levels of contaminants did not follow this trend. XPS result indicates that, after the initial SH-ssDNA adsorption, the remaining material incorporated into these films was due to contamination.     

The structure of phosphine-functionalised silsesquioxane-based dendrimers: a molecular dynamics study

Molecular dynamics modelling has been used to simulate the structures of phopshine-functionalised, silsesquioxane-based dendrimers whose excellent catalytic properties have been previously demonstrated. The effect of changing the chemical composition of the dendrimer branches is simulated. The results indicate that adding a methylene unit to a branch increases the overall size of the dendrimer but replacing a methylene unit with an oxygen atom decreases the size of the dendrimer. The size and shape of the dendrimers have also been simulated on changing the temperature and polarity of the solvent. The distribution of phosphine groups on the exterior has also been modelled and this suggests that all are available for bonding to catalytic metals in all the compounds.     

Temperature effects on the structure and dynamics of liquid dimethyl sulfoxide: A molecular dynamics study

The molecular dynamics (MD) simulation technique has been employed to investigate the thermodynamic properties and transport coefficients of the neat liquid dimethyl sulfoxide (DMSO). The fluid has been studied at temperatures in the range 298–353 K and at a pressure equal to 1 atm. The simulations employed a nine-site potential model, which is presented for the first time here, and all the available non-polarizable models. The performance of each model is tested using the same statistical mechanical ensemble and simulation method under the same conditions, revealing its weaknesses and strengths. Thermodynamic properties, microscopic structure and dynamic properties, such as transport coefficients, rotational and single-dipole correlation times have been calculated and compared with available experimental results. Estimations of transport coefficients from various theoretical and empirical models are tested against experimental and MD results. Translational and rotational dynamics suggest the existence of the cage effect and agree with the Stokes–Einstein–Debye relation. The dipole relaxation times calculated are discussed in terms of simple and useful approximations, such as the Glarum–Powles and Fatuzzo–Mason models.     

Proline-based chiral stationary phases: a molecular dynamics study of the interfacial structure

Proline chains have generated considerable interest as a possible basis for new selectors in chiral chromatography. In this article, we employ molecular dynamics simulations to examine the interfacial structure of two diproline chiral selectors, one with a terminal trimethylacetyl group and one with a terminal t-butyl carbamate group. The solvents consist of a relatively apolar n-hexane/2-propanol and a polar water/methanol mixture. We begin with electronic structure calculations for the two chiral selectors to assess the energetics of conformational changes, particularly along the backbone where the amide bonds can alternate between cis and trans conformations. Force fields have been developed for the two selectors, based on these ab initio calculations. Molecular dynamics simulations of the selective interfaces are performed to examine the preferred backbone conformations, as a function of end-group and solvent. The full chiral surface includes the diproline selectors, trimethylsilyl end-caps, and silanol groups. Connection is made with selectivity measurements on these interfaces, where significant differences are observed between these two very similar selectors.     

On the structure and desorption dynamics of DNA bases adsorbed on gold: a temperature-programmed study

The structure and desorption dynamics of mono- and multilayer samples of adenine, cytosine, guanine, and thymine on polycrystalline gold thin films are studied using temperature-programmed desorption-infrared reflection absorption spectroscopy (TPD-IRAS) and temperature-programmed desorption-mass spectroscopy (TPD-MS). It is shown that the pyrimidines, adenine and guanine, adsorb to gold in a complex manner and that both adhesive (adenine) and cohesive (guanine) interactions contribute the apparent binding energies to the substrate surface. Adenine displays at least two adsorption sites, including a high-energy site (210 degrees C, approximately 136 kJ/mol), wherein the molecule coordinates to the gold substrate via the NH2 group in an sp3-like, strongly perturbed, nonplanar configuration. The purines, cytosine and thymine, display a less complicated adsorption/desorption behavior. The desorption energy for cytosine (160 degrees C, approximately 122 kJ/mol) is similar to those obtained for adenine and guanine, but desorption occurs from a single site of dispersed, nonaggregated cytosine. Thymine desorbs also from a single site but at a significantly lower energy (100 degrees C, approximately 104 kJ/mol). Infrared data reveal that the monolayer architectures discussed herein are structurally very different from those observed for the bases in the bulk crystalline state. It is also evident that both pyrimidines and purines adsorb on gold with the plane of the molecule in a nonparallel orientation with respect to the substrate surface. The results of this work are discussed in the context of improving the understanding of the design of capturing oligonucleotides or DNA strands for bioanalytical applications, in particular, for gold nanoparticle-based assays.     

Quinoprotein methanol dehydrogenase: a molecular dynamics study and comparison with crystal structure

One of the mechanisms proposed for methanol oxidation by methanol dehydrogenase (MDH) involves a hydride transfer to the quinone carbonyl carbon C5 of 2,7,9-tricarboxy-1H-pyrrolo[2,3-f]-quinoline-4,5-dione (PQQ), initiated by abstraction of a proton from the substrate methanol by a general base. Molecular dynamics studies are performed on MDH-bound to the C5 reduced intermediate (C5RI) of PQQ, for 3 ns. The structural features of the MD and X-ray structures are compared. An interesting feature observed during simulations is the strong hydrogen bond between oxyanion O5 of C5RI and Asp297-CO₂H in the active site. Asp297-CO₂⁻ is suggested to be the general-base catalyst for removing the hydroxyl proton of methanol in concert with the hydride ion transfer from the putative methoxide to C5 carbonyl of PQQ. The formed Asp297-CO₂H acts as the required proton source for the immediate product. Anticorreleated motions observed in the MD structure are not across the active site to influence the reaction mechanism of MDH.     

Internal structure of dendrimers in the melt under shear: a molecular dynamics study

The molecular structure of fluids composed of dendrimers of different generations is studied using nonequilibrium molecular dynamics (NEMD). NEMD results for dendrimer melts undergoing planar Couette flow are reported and analyzed with particular attention paid to the shear-induced changes in the internal structure of dendrimers. The radii of gyration, pair distribution functions and the fractal dimensionality of the dendrimers are determined at different strain rates. The location of the terminal groups is analyzed and found to be uniformly distributed throughout the space occupied by the molecules. The fractal dimension as a function of strain rate displays crossover behavior analogous to the Newtonian/non-Newtonian transition of shear viscosity.     

A molecular dynamics study of Lys-Trp-Lys: structure and dynamics in solution following photoexcitation

We report studies of the structure and dynamics of a tripeptide Lys-Trp-Lys (KWK) in aqueous solution following photoexcitation by molecular dynamics simulations. For ground-state KWK, we observe three stable conformations with free energy differences of less than 5.2 kJ/mol. Each conformer is stabilized by a pi-cation interaction between one of three protonated amino groups and the indole moiety. For the excited state of tryptophan in KWK, the simulated molecular dynamics of the three isomers are similar, all in good agreement with recent femtosecond experiments (J. Phys. Chem. B 2005, 109, 16901). Specifically, we observe: (1) the fluorescence anisotropy is dominated by a single-exponential component and decays in approximately 130 ps, (2) the total dynamic Stokes shift reaches approximately 2700 cm(-1), and (3) the excited state relaxation dynamics occurs on several time scales ranging from femtoseconds to tens of picoseconds. The relaxation dynamics involve rapid initial response of neighboring water, followed by local motions of flexible peptide chains. These processes drive global restructuring of the tripeptide on a rather flat energy surface, inducing slower dynamics evident in both the water and protein contributions to the stabilization energy of the photoexcited chromophore. The water and protein dynamics are strongly correlated. On a still longer time scale, we observe isomerization of two excited state conformers to the other most stable one, an analogue for evolution of trajectories along the funnel on the rugged free energy landscape to the final "native" state. Our studies suggest new experiments to detect this unique dynamics.     

The effect of pressure on the hydration structure around hydrophobic solute: a molecular dynamics simulation study

Molecular dynamics simulations are performed to study the effects of pressure on the hydrophobic interactions between neopentane molecules immersed in water. Simulations are carried out for five different pressure values ranging from 1 atm to 8000 atm. From potential of mean force calculations, we find that with enhancement of pressure, there is decrease in the well depth of contact minimum (CM) and the relative stability of solvent separated minimum over CM increases. Lower clustering of neopentane at high pressure is also observed in association constant and cluster-structure analysis. Selected site-site radial distribution functions suggest efficient packing of water molecules around neopentane molecules at elevated pressure. The orientational profile calculations of water molecules show that the orientation of water molecules in the vicinity of solute molecule is anisotropic and this distribution becomes flatter as we move away from the solute. Increasing pressure slightly changes the water distribution. Our hydrogen bond properties and dynamics calculations reveal pressure-induced formation of more and more number of water molecules with five and four hydrogen bond at the expense of breaking of two and three hydrogen bonded water molecules. We also find lowering of water-water continuous hydrogen bond lifetime on application of pressure. Implication of these results for relative dispersion of hydrophobic molecules at high pressure are discussed.     

Effect of temperature and pressure on the structure, dynamics, and hydrogen bond properties of liquid N-methylacetamide: a molecular dynamics study

The structure and dynamical properties of liquid N-methylacetamides (NMA) are calculated at five different temperatures and at four different pressures using classical molecular dynamics simulations. Our results are analyzed in terms of pressure-induced changes in structural properties by investigating the radial distribution functions of different atoms in NMA molecule. It is found that the first peak and also the second peak of C-O and N-H are well defined even at higher temperature and pressure. It is also observed that the number of hydrogen bonds increase with application of pressure at a given temperature. On the other hand, the calculated hydrogen bond energy (E(HB)) shows that the stability of hydrogen bond decreases with increasing of pressure and temperature. Various dynamical properties associated with translational and rotational motion of neat NMA are calculated and the self-diffusion coefficient of NMA is found to be in excellent agreement with the experiment and the behavior is non-Arrhenius at low temperatures with application of pressures. The single particle orientational relaxation time for dipole vector and N-C vector are also calculated and it is found that the orientational relaxation time follows Arrhenius behavior with a variation of temperature and pressure.     

Cooperativity in drug-DNA recognition: a molecular dynamics study.

NMR studies have shown that the minor groove-binding ligand Hoechst 33258 binds to the two T4/A4 tracts within the duplex d(CTTTTCGAAAAG)2 in a highly cooperative manner, such that in titration experiments no intermediate 1:1 complex can be detected. The NMR-derived structures of the free DNA and the 2:1 complex have been obtained, but can shed little light on what the origins of this cooperativity may be. Here we present the results of a series of molecular dynamics simulations on the free DNA, the 1:1 complex, and the 2:1 complex, which have been designed to enable us to calculate thermodynamic parameters associated with the molecular recognition events. The results of the molecular dynamics studies confirm that structural factors alone cannot explain the cooperativity observed, indeed when enthalpic and hydration factors are looked at in isolation, the recognition process is predicted to be slightly anticooperative. However, when changes in configurational entropy are taken into account as well, the overall free energy differences are such that the calculated cooperativity is in good agreement with that observed experimentally. The results indicate the power of molecular dynamics methods to provide reasonable explanations for phenomena that are difficult to explain on the basis of static models alone, and provide a nice example of the concept of "allostery without conformational change".     

Influence of a charged graphene surface on the orientation and conformation of covalently attached oligonucleotides: a molecular dynamics study

Molecular dynamics (MD) simulations of single-stranded (ss) and double-stranded (ds) oligonucleotides anchored via an aliphatic linker to a graphene surface were performed in order to investigate the role of the surface charge density in the structure and orientation of attached DNA. Two types of interactions of DNA with the surface are crucial for the stabilisation of the DNA-surface system. Whereas for a surface with a zero or low positive charge density the dispersion forces between the base(s) and the surface dominate, the higher charge densities applied on the surface lead to a strong electrostatic interaction between the phosphate groups of DNA, the surface and the ions. At high-charge densities, the interaction of the DNA with the surface is strongly affected by the formation of a low-mobility layer of counterions compensating for the charge of the surface. A considerable difference in the behaviour of the ds-DNA and ss-DNA anchored to the layer was observed. The ds-DNA interacts with the surface at low- and zero-charge densities exclusively by the nearest base pair. It keeps its geometry close to the canonical B-DNA form, even at surfaces with high-charge densities. The ss-DNA, owing to its much higher flexibility, has a tendency to maximise the attraction to the surface exploiting more bases for the interaction. The interaction of the polar amino group(s) of the base(s) of ss-DNA with a negatively charged surface also contributes significantly to the system stability.