The self-diffusion and hydrogen bond interaction in neat liquid alkanols: a molecular dynamic simulation study

Self-diffusion of methanol, ethanol, 1-propanol and 2-propanol has been studied by molecular dynamics simulation in the temperature range between the melting pressure curve and 478 K at pressures up to 300 MPa. The simulation results on self-diffusion of methanol, ethanol and 2-propanol (for 2-propanol, at high temperatures) agree well with experiment, which suggests that the simulation method is a powerful tool to obtain self-diffusion coefficients over wide range of temperature and pressure, under which it is rather difficult for experiments. The local structures of methanol, ethanol and 2-propanol are investigated by calculating the radial distribution functions, H-bond numbers, coordination numbers and the ratios of H-bond number divided by coordination number. The correlation between self-diffusion and structural properties, and the influence of temperature and pressure on them are discussed. The degree of forming H-bond space network in methanol, ethanol and water is higher than that in 2-propanol, and they are all higher than those in ammonia and methylamine. The simulation results demonstrate that the effect of hydrogen bonding on the translational dynamics in methanol and ethanol is more pronounced than that in 2-propanol.     

Structural and dynamic differences of the estrogen receptor DNA-binding domain, binding as a dimer and as a monomer to DNA: molecular dynamics simulation studies

Molecular dynamics (MD) simulations of the estrogen receptor DNA-binding domain (ERDBD) as a dimer in complex with its DNA response element (ERE) show a significant difference in both structure and dynamics, compared to a MD simulation of monomeric ERDBD bound to its half-site response element (EREH). The C-terminal zinc binding domain (Zn_II), including a region (helix II) which is in a helical conformation in ERE-(ERDBD)₂, is considerably more flexible in EREH-ERDBD than in the dimeric complex. In EREH-ERDBD, all helical hydrogen bonds in helix II are broken and the entire Zn_II region is detached from a hydrogen bonding network that in ERE-(ERDBD)₂ connects to other parts of the protein as well as to the DNA. The regions that become flexible in EREH-ERDBD are identical to the regions where the NMR solution structure of free ERDBD is poorly ordered. This strongly suggests that dimerisation of ERDBD is required for ordering of the Zn_II region and that monomeric binding to DNA is not sufficient for the ordering. This contrasts to the glucocorticoid receptor DNA-binding domain (GRDBD) which has essentially the same mobility (uniform and limited), regardless of whether it is free as a monomer in solution, bound as a monomer to its half-site response element or in a dimeric complex with the full response element. The hydrogen bonding network that connects Zn_II with other parts of the protein and to DNA is almost identical in ERDBD and GRDBD. However, in GRDBD there is also a serine (in the N-terminal zinc coordinating region) with a central role in this network, connecting to the Zn_II region. This serine is replaced by a glycine in ERDBD and we suggest that this substitution is sufficient for destabilisation of the network, thus leading to a more flexible Zn_II region, which becomes ordered first upon forming a complex with another ERDBD and DNA. Received: 6 March 1998 / Revised version: 22 June 1998 / Accepted: 2 September 1998     

Experimental and computational studies on the effects of valganciclovir as an antiviral drug on calf thymus DNA

DNA-binding properties of an antiviral drug, valganciclovir (valcyte) was studied by using emission, absorption, circular dichroism, viscosity, differential pulse voltammetry, fluorescence techniques, and computational studies. The drug bound to calf thymus DNA (ct-DNA) in a groove-binding mode. The calculated binding constant of UV-vis, K_a, is comparable to groove-binding drugs. Competitive fluorimetric studies with Hoechst 33258 showed that valcyte could displace the DNA-bound Hoechst 33258. The drug could not displace intercalated methylene blue from DNA double helix. Furthermore, the induced detectable changes in the CD spectrum of ct-DNA as well as changes in its viscosity confirm the groove-binding mode. In addition, an integrated molecular docking was employed to further investigate the binding interactions between valcyte and calf thymus DNA.     

Triazolopyridyl ketones as a novel class of antileishmanial agents. DNA binding and BSA interaction

A new series of triazolopyridyl pyridyl ketones has been synthetized by regioselective lithiation of the corresponding [1,2,3]triazolo[1,5-a]pyridine at 7 position followed by reaction with different electrophiles. The in vitro antileishmanial activity of these compounds was evaluated against Leishmania infantum, Leishmania braziliensis, Leishmania guyanensis and Leishmania amazonensis. Compounds 6 and 7 were found to be the most active leishmanicidal agents. Both of them showed activities at micromolar concentration against cultured promastigotes of Leishmania spp. (IC₅₀ = 99.8–26.8 μM), without cytotoxicity on J774 macrophage cells. These two compounds were also tested in vivo in a murine model of acute infection by L. infantum. The triazolopyridine derivative 6 was effective against both spleen and liver parasites forms, while 7 was inactive against liver parasites. Mechanistic aspects of the antileishmanial activity were investigated by means of DNA binding studies (UV-titration and viscosimetry). Results have revealed that these active ligands are able to interact strongly with DNA [K_b = 1.14 × 10⁵ M⁻¹ (6) and 3.26 × 10⁵ M⁻¹ (7)]. Moreover, a DNA groove binding has been proposed for both 6 and 7. To provide more insight on the mode of action of compounds 6 and 7 under biological conditions, their interaction with bovine serum albumin (BSA) was monitored by fluorescence titrations and UV–visible spectroscopy. The quenching constants and binding parameters were determined. Triazolopyridine ketones 6 and 7 have exhibited significant affinity towards BSA [K_b = 2.5 × 10⁴ M⁻¹ (6) and 1.9 × 10⁴ M⁻¹ (7)]. Finally, to identify the binding location of compounds 6 and 7 on the BSA, competitive binding experiments were carried out, using warfarin, a characteristic marker for site I, and ibuprofen as one for site II. Results derived from these studies have indicated that both compounds interact at BSA site I and, to a lesser extent, at site II.     

Molecular dynamic simulations on an inhibitor of anti-apoptotic Bcl-2 proteins for insights into its interaction mechanism for anti-cancer activity

Inhibition of normal cellular apoptosis or programed cell death is the hallmark of all cancers. Apoptotic dysregulation can result in numerous pathological conditions, such as cancers, autoimmune disorders, and neurodegeneration. Members of the BCL-2 family of proteins regulate the process of apoptosis by its promotion or inhibition and overexpression of the pro-survival anti-apoptotic proteins (Bcl-2, Bcl-xL, and Mcl-1) has been associated with tumor maintenance, growth and progression Small molecules and peptides which bind the BH3 binding groove of these proteins have been explored in the recent times for their anticancer potential. The first anticancer agents targeting this family of proteins were aimed primarily toward inhibition of Bcl-2. An elevated level of Mcl-1, despite Bcl-2 inhibition, continues to be a cause for resistance in most cancers. However, in silico exploration of Mcl-1 specific drugs and their associated mechanisms have not been clearly elucidated. In order to understand the same, we have carried out docking and molecular dynamic simulations on ABT-263 (Navitoclax), an orally active inhibitor of Bcl-2, Bcl-xL, and Bcl-w proteins; Obatoclax, a pan-Bcl-2 inhibitor as well as Maritoclax, an Mcl-1 specific inhibitor. Docking studies revealed that binding to the hydrophobic grooves is a prerequisite for action on the BCL protein and the binding mechanism and chemical space utilization dictates stability as well as specificity of the inhibitor molecular dynamic simulations showed that on binding, the α-helices of these proteins exhibited less fluctuations than loop regions, also hydrophobic contacts and hydrogen bonding were observed to be the predominant interactions in the drug-receptor complexes.

Communicated by Ramaswamy H. Sarma     

Atomistic mechanism of polyphenol amyloid aggregation inhibitors: molecular dynamics study of Curcumin,Exifone, and Myricetin interaction with the segment of tau peptide oligomer

Amyloid fibrils are highly ordered protein aggregates associated with many diseases affecting millions of people worldwide. Polyphenols such as Curcumin, Exifone, and Myricetin exhibit modest inhibition toward fibril formation of tau peptide which is associated with Alzheimer’s disease. However, the molecular mechanisms of this inhibition remain elusive. We investigated the binding of three polyphenol molecules to the protofibrils of an amyloidogenic fragment VQIVYK of tau peptide by molecular dynamics simulations in explicit solvent. We find that polyphenols induce conformational changes in the oligomer aggregate. These changes disrupt the amyloid H bonding, perturbing the aggregate. While the structural evolution of the control oligomer with no ligand is limited to the twisting of the β-sheets without their disassembly, the presence of polyphenol molecule pushes the β-sheets apart, and leads to a loosely packed structure where two of four β-sheets dissociate in each of the three cases considered here. The H-bonding capacity of polyphenols is responsible for the observed behavior. The calculated binding free energies and its individual components enabled better understanding of the binding. Results indicated that the contribution from Van der Waals interactions is more significant than electrostatic contribution to the binding. The findings from this study are expected to assist in the development of aggregation inhibitors. Significant binding between polyphenols and aggregate oligomer identified in our simulations confirms the previous experimental observations in which polyphenols refold the tau peptide without forming covalent bonds.     

Spectroscopic and molecular docking studies on the interaction of antiviral drug nevirapine with calf thymus DNA

The interaction of calf thymus DNA with nevirapine at physiological pH was studied by using absorption, circular dichroism, viscosity, differential pulse voltammetry, fluorescence techniques, salt effect studies and computational methods. The drug binds to ct-DNA in a groove binding mode, as shown by slight variation in the viscosity of ct-DNA. Furthermore, competitive fluorimetric studies with Hoechst 33258 indicate that nevirapine binds to DNA via groove binding. Moreover, the structure of nevirapine was optimized by DFT calculations and was used for the molecular docking calculations. The molecular docking results suggested that nevirapine prefers to bind on the minor groove of ct-DNA.     

Molecular Modeling Study on the Enantioselective Binding of S- and R-Ofloxacin to Various DNA Sequences

Abstract

The complex formation of S- and R-ofloxacin with the self-complementary oligonucleotides, namely d[ATAGCGCTAT]₂, d[GCGATATCGC]₂ and d[ATAICICTAT]₂, were investigated by the molecular dynamics (MD) simulation. Four starting positions, including two intercalation positions with different insertion directions and two minor groove binding positions, were considered. The total energy of both S- and R-ofloxacin-d[ATAGCGCTAT]₂ complex, in which ofloxacin binds in the minor groove of the oligonucleotide, were lower than any intercalation binding mode. For both enantiomers, formation of the complex with GC oligonucleotide is more favorable than AT and IC oligonucleotides. When S- and R-ofloxacin are compared, the S-enantiomer exhibits more favorable total energy and torsion angles in the complex formation. This result is in agreement with the experimental observation [Hwangbo et al., Eur J Pharm Sci 18, 197 (2003)]. In the complex, both enantiomers form two hydrogen bonds: one between the carbonyl group of ofloxacin and the amine group of G16 and the other between the fluorine group and the G6 amine for S-ofloxacin. However, only one hydrogen bond is formed between endocyclic hydrogen atom at the C2 position of adenine and inosine base and carbonyl group of ofloxacin, which may be the reason for the GC preferentiality of ofloxacin.     

Chrysophanol anthrone: a theoretical study on the potential energy surface

The potential energy surface (PES) of chrysophanol anthrone, the active component of Goa Powder, was systematically explored and thoroughly scrutinised via density functional theory, in order to gain an understanding of its physicochemical properties. In particular, we focused on the rotations of the two hydroxyl-phenyl dihedral angles. A picture with a four stable rotamers emerged where only the most stable conformer has a planar structure and the less stable conformer has the maximum deviation from planarity. The computed PES shows that the energy barriers for the conformer interconversion are less than 15 kcal/mol. From the analysis of the calculated intramolecular hydrogen bond enthalpy, we conclude that the number of the intramolecular hydrogen bonds governs the conformer stability. Additionally, the conformational equilibrium was pursued by means of an analysis of the energy of OH internal rotation barriers. The total energy changes were decomposed in an electrostatic decomposition scheme in order to gain an insight into the effects governing the torsional barrier and preferred conformations. This analysis shows that the interplay between the repulsive and attractive potentials causes the conformer stability, where the attractive term dominates the conformer stabilisation.     

Electrochemical spectroscopic investigations on the interaction of an ytterbium complex with DNA and their analytical applications such as biosensor

Metal ion-DNA interactions are important in nature, often changing the genetic material's structure and function. A new Yb complex of YbCl₃ (tris(8-hydroxyquinoline-5-sulfonic acid) ytterbium) was synthesized and utilized as an electrochemical indicator for the detection of DNA oligonucleotide based on its interaction with Yb(QS)₃. Cyclic voltammetry (CV) and fluorescence spectroscopy were used to investigate the interaction of Yb(QS)₃ with ds-DNA. It was revealed that Yb(QS)₃ presented an excellent electrochemical activity on glassy carbon electrode (GCE) and could intercalate into the double helix of double-stranded DNA (ds-DNA). The binding mechanism of interaction was elucidated on glassy carbon electrode dipped in DNA solution and DNA modified carbon paste electrode by using differential pulse voltammetry and cyclic voltammetry. The binding ratio between this complex and ds-DNA was calculated to be 1:1. The extent of hybridization was evaluated on the basis of the difference between signals of Yb(QS)₃ with probe DNA before and after hybridization with complementary DNA. With this approach, this DNA could be quantified over the range from 1 × 10⁻⁸ to 1.1 × 10⁻⁷ M. The interaction mode between Yb(QS)₃ and DNA was found to be mainly intercalative interaction. These results were confirmed with fluorescence experiments.     

Molecular dynamics simulation study of the static and dynamic properties of a model colloidal suspension with explicit solvent

Molecular dynamics simulation was used to study a colloidal suspension with explicit solvent to determine how inclusion of the solvent affects the structure and dynamics of the system. The solute was modelled as a hard-core particle enclosed in a Weeks–Chandler–Andersen (WCA) potential shell, while the solvent was modelled as a simple WCA fluid. We found that when the solute–solvent interaction included a hard core equal to half of the solute hard-core diameter, large depletion effects arose, leading to an effective attraction and large deviations from hard-sphere structure for the colloidal component. It was found that these effects could be eliminated by reducing the hard-core distance parameter in the solute–solvent interaction, thus allowing the solvent to penetrate closer to the colloidal particles. Three different values for the solute–solvent hard-core parameter were systematically studied by comparing the static structure factor and radial distribution function to the predictions of the Percus–Yevick theory for hard spheres. When the optimal value of the solute–solvent hard-core interaction parameter was found, this model was then used to study the dynamical behaviour of the colloidal suspension. This was done by first measuring the velocity autocorrelation function (VACF) over a large range of packing fractions. We found that this model predicted the sign of the long-time tail in the VACF in agreement with experimental values, something that single component hard-sphere systems have failed to do. The intermediate scattering functions at low wavevector were briefly studied to determine their behaviour in a dilute system. It was found that they could be modelled using a simple diffusion equation with a wavevector independent diffusion coefficient, making this model an excellent analogue of experimentally studied hard-sphere colloids.     

A comparative study on the interaction of the putative anticancer alkaloids,sanguinarine and chelerythrine,with single- and double-stranded,and heat-denatured DNAs

A detailed investigation on the interaction of two benzophenanthridine alkaloids, sanguinarine (SGR) and chelerythrine (CHL), with the double-stranded (ds), heat-denatured (hd), and single-stranded (ss) DNA was performed by spectroscopy and calorimetry techniques. Binding to the three DNA conformations leads to quenching of fluorescence of SGR and enhancement in the fluorescence of CHL. The binding was cooperative for both of the alkaloids with all the three DNA conformations. The binding constant values of both alkaloids with the ds DNA were in the order of 10⁶ M^?1; binding was weak with hd and much weaker to the ss DNA. The fluorescence emission of the alkaloid molecules bound to the ds and hd DNAs was quenched much less compared to those bound to the ss DNA based on competition with the anionic quencher KI. For both double stranded and heat denatured structures the emission of the bound alkaloid molecules was polarized significantly and strong energy transfer from the DNA bases to the alkaloid molecules occurred. Intercalation of SGR and CHL to ds, hd, and ss DNA was proved from these fluorescence results. Calorimetric studies suggested that the binding to all DNA conformations was both enthalpy and entropy favored. Both the alkaloids preferred double-helical regions for binding, but SGR was a stronger binder than CHL to all the three DNA structures.     

Electronic and fluorescent studies on the interaction of DNA and BSA with a new ternary praseodymium complex containing 2,9-dimethyl 1,10-phenanthroline and antibacterial activities testing

In this study, fluorescence emission spectra, UV–vis absorption spectra, ethidium bromide (EB)-competition experiment, and iodide quenching experiment were used for the interaction study of the Fish salmon DNA (FS-DNA) with [Pr(dmp)2Cl3(OH2)] where dmp is 2,9-dimethyl 1,10-phenanthroline. The binding constant and the number of binding sites of the complex with FS-DNA were 6.09?±?0.04 M^?1 and 1.18, respectively. The free energy, enthalpy, and entropy changes (ΔG°, ΔH°, and ΔS°) in the binding process of the Pr(III) complex with FS-DNA were –8.02?kcal mol^?1, +39.44?kcal mol^?1, and +159.56?cal mol^?1 K^?1, respectively. Based on these results, the interaction process between FS-DNA with [Pr(dmp)2Cl3(OH2)] was spontaneous and the main binding interaction force was groove binding mode. Also, Fluorescence and electronic absorption spectroscopy were used in order to evaluate the binding characteristics, stoichiometry, and interaction mode of praseodymium(III) (Pr(III)) complex with bovine serum albumin (BSA). Title complex showed good binding propensity to BSA presenting moderately high Kb values. The fluorescence quenching of BSA by Pr(III) complex has been observed to be the static process. The positive ΔH° and ΔS° values showed that the hydrophobic interaction is the main force in the binding of Pr(III) complex and BSA. Eventually, the average aggregation number, <J>, of BSA potentially induced by title complex confirmed the 1:1 stoichiometry for title complex-BSA adducts. In vitro, antimicrobial activity of title complex was indicated that the complex is more active against both Escherichia coli and Enterococcus faecalis bacterial strains than Staphylococcus aureus, and Pseudomonas aeruginosa.

Communicated by Ramaswamy H. Sarma     

The study on the interaction between phytosterols and phospholipids in model membranes

Sterols are one of the major components of cellular membranes. Although in mammalian membranes cholesterol is a predominant sterol, in the human organism plant sterols (phytosterols) can also be found. Phytosterols, especially if present in concentrations higher than normal (phytosterolemia), may strongly affect membrane properties. In this work, we studied phytosterol-phospholipid interactions in mixed Langmuir monolayers serving as model membranes. Investigated were two phytosterols, beta-sitosterol and stigmasterol and a variety of phospholipids, both phosphatidylethanolamines and phosphatidylcholines. The phospholipids had different polar heads, different length and saturation of their hydrocarbon chains. The interactions between molecules in mixed sterol/phospholipid films were characterized with the mean area per molecule (A(12)) and the excess free energy of mixing (DeltaG(Exc)). The effect of the sterols on the molecular organization of the phospholipid monolayers was analyzed based on the compression modulus values. It was found that the incorporation of the phytosterols into the phospholipid monolayers increased their condensation. The plant sterols revealed higher affinity towards phosphatidylcholines as compared to phosphatidylethanolamines. The phytosterols interacted more strongly with phospholipids possessing longer and saturated chains. Moreover, both the length and the saturation of the phosphatidylcholines influenced the stoichiometry of the most stable complexes. Our results, compared with those presented previously for cholesterol/phospholipid monolayers, allowed us to draw a conclusion that the structure of sterol (cholesterol, beta-sitosterol, stigmasterol) does not affect the stoichiometry of the most stable complexes formed with particular phospholipids, but influences their stability. Namely, the strongest interactions were found for cholesterol/phospholipids mixtures, while the weakest for mixed systems containing stigmasterol.     

Molecular simulation investigation on the interaction between barrier-to-autointegration factor or its Gly25Glu mutant and DNA

Journal of Molecular Modeling - In this work, the capacity of three different imidazolium-based ionic liquids (ILs) for atmospheric mercury capture has been evaluated. Theoretical calculations...     

Molecular dynamics studies of a DNA-binding protein: 2. An evaluation of implicit and explicit solvent models for the molecular dynamics simulation of the Escherichia coli trp repressor.

Although aqueous simulations with periodic boundary conditions more accurately describe protein dynamics than in vacuo simulations, these are computationally intensive for most proteins. Trp repressor dynamic simulations with a small water shell surrounding the starting model yield protein trajectories that are markedly improved over gas phase, yet computationally efficient. Explicit water in molecular dynamics simulations maintains surface exposure of protein hydrophilic atoms and burial of hydrophobic atoms by opposing the otherwise asymmetric protein-protein forces. This properly orients protein surface side chains, reduces protein fluctuations, and lowers the overall root mean square deviation from the crystal structure. For simulations with crystallographic waters only, a linear or sigmoidal distance-dependent dielectric yields a much better trajectory than does a constant dielectric model. As more water is added to the starting model, the differences between using distance-dependent and constant dielectric models becomes smaller, although the linear distance-dependent dielectric yields an average structure closer to the crystal structure than does a constant dielectric model. Multiplicative constants greater than one, for the linear distance-dependent dielectric simulations, produced trajectories that are progressively worse in describing trp repressor dynamics. Simulations of bovine pancreatic trypsin were used to ensure that the trp repressor results were not protein dependent and to explore the effect of the nonbonded cutoff on the distance-dependent and constant dielectric simulation models. The nonbonded cutoff markedly affected the constant but not distance-dependent dielectric bovine pancreatic trypsin inhibitor simulations. As with trp repressor, the distance-dependent dielectric model with a shell of water surrounding the protein produced a trajectory in better agreement with the crystal structure than a constant dielectric model, and the physical properties of the trajectory average structure, both with and without a nonbonded cutoff, were comparable.     

Role of DNA conformation & energetic insights in Msx-1-DNA recognition as revealed by molecular dynamics studies on specific and nonspecific complexes

In most of homeodomain–DNA complexes, glutamine or lysine is present at 50th position and interacts with 5th and 6th nucleotide of core recognition region. Molecular dynamics simulations of Msx-1–DNA complex (Q50-TG) and its variant complexes, that is specific (Q50K-CC), nonspecific (Q50-CC) having mutation in DNA and (Q50K-TG) in protein, have been carried out. Analysis of protein–DNA interactions and structure of DNA in specific and nonspecific complexes show that amino acid residues use sequence-dependent shape of DNA to interact. The binding free energies of all four complexes were analysed to define role of amino acid residue at 50th position in terms of binding strength considering the variation in DNA on stability of protein–DNA complexes. The order of stability of protein–DNA complexes shows that specific complexes are more stable than nonspecific ones. Decomposition analysis shows that N-terminal amino acid residues have been found to contribute maximally in binding free energy of protein–DNA complexes. Among specific protein–DNA complexes, K50 contributes more as compared to Q50 towards binding free energy in respective complexes. The sequence dependence of local conformation of DNA enables Q50/Q50K to make hydrogen bond with nucleotide(s) of DNA. The changes in amino acid sequence of protein are accommodated and stabilized around TAAT core region of DNA having variation in nucleotides.