Charge structure and counterion distribution in hexagonal DNA liquid crystal

A hexagonal liquid crystal of DNA fragments (double-stranded, 150 basepairs) with tetramethylammonium (TMA) counterions was investigated with small angle neutron scattering (SANS). We obtained the structure factors pertaining to the DNA and counterion density correlations with contrast matching in the water. Molecular dynamics (MD) computer simulation of a hexagonal assembly of nine DNA molecules showed that the inter-DNA distance fluctuates with a correlation time around 2 ns and a standard deviation of 8.5% of the interaxial spacing. The MD simulation also showed a minimal effect of the fluctuations in inter-DNA distance on the radial counterion density profile and significant penetration of the grooves by TMA. The radial density profile of the counterions was also obtained from a Monte Carlo (MC) computer simulation of a hexagonal array of charged rods with fixed interaxial spacing. Strong ordering of the counterions between the DNA molecules and the absence of charge fluctuations at longer wavelengths was shown by the SANS number and charge structure factors. The DNA-counterion and counterion structure factors are interpreted with the correlation functions derived from the Poisson-Boltzmann equation, MD, and MC simulation. Best agreement is observed between the experimental structure factors and the prediction based on the Poisson-Boltzmann equation and/or MC simulation. The SANS results show that TMA is too large to penetrate the grooves to a significant extent, in contrast to what is shown by MD simulation.     

Computer modeling demonstrates that electrostatic attraction of nucleosomal DNA is mediated by histone tails

We conducted molecular dynamics computer simulations of charged histone tail-DNA interactions in systems mimicking nucleosome core particles (NCP) . In a coarse-grained model, the NCP is modeled as a negatively charged spherical particle with flexible polycationic histone tails attached to it in a dielectric continuum with explicit mobile counterions and added salt. The size, charge, and distribution of the tails relative to the core were built to mimick real NCP. In this way, we incorporate attractive ion-ion correlation effects due to fluctuations in the ion cloud and the attractive entropic and energetic tail-bridging effects. In agreement with experimental data, increase of monovalent salt content from salt-free to physiological concentration leads to the formation of NCP aggregates; likewise, in the presence of MgCl₂, the NCPs form condensed systems via histone-tail bridging and accumulation of counterions. More detailed mechanisms of the histone tail-DNA interactions and dynamics have been obtained from all-atom molecular dynamics simulations (including water), comprising three DNA 22-mers and 14 short fragments of the H4 histone tail (amino acids 5–12) carrying three positive charges on lysine⁺ interacting with DNA. We found correlation of the DNA-DNA distance with the presence and association of the histone tail between the DNA molecules.     

Structural flexibility of the nucleosome core particle at atomic resolution studied by molecular dynamics simulation

Comparative explicit solvent molecular dynamics (MD) simulations have been performed on a complete nucleosome core particle with and without N-terminal histone tails for more than 20 ns. Main purpose of the simulations was to study the dynamics of mobile elements such as histone N-terminal tails and how packing and DNA-bending influences the fine structure and dynamics of DNA. Except for the tails, histone and DNA molecules stayed on average close to the crystallographic start structure supporting the quality of the current force field approach. Despite the packing strain, no increase of transitions to noncanonical nucleic acid backbone conformations compared to regular B-DNA was observed. The pattern of kinks and bends along the DNA remained close to the experiment overall. In addition to the local dynamics, the simulations allowed the analysis of the superhelical mobility indicating a limited relative mobility of DNA segments separated by one superhelical turn (mean relative displacement of approximately +/-0.2 nm, mainly along the superhelical axis). An even higher rigidity was found for relative motions (distance fluctuations) of segments separated by half a superhelical turn (approximately +/-0.1 nm). The N-terminal tails underwent dramatic conformational rearrangements on the nanosecond time scale toward partially and transiently wrapped states around the DNA. Many of the histone tail changes corresponded to coupled association and folding events from fully solvent-exposed states toward complexes with the major and minor grooves of DNA. The simulations indicate that the rapid conformational changes of the tails can modulate the DNA accessibility within a few nanoseconds.     

Molecular dynamics simulations of DNA and a protein-DNA complex including solvent

The results of a recent nanosecond (ns) molecular dynamics (MD) simulation of the d(CGCGAATTCGCG) double helix in water and a 100 ps MD study of the repressor-operator complex are described. The DNA simulations are analyzed in terms of the structural dynamics, fluctuations in the groove width and bending of the helical axis. The results indicate that the ns dynamical trajectory progresses through a series of three substates of B form DNA, with lifetimes of the order of hundreds of picoseconds (ps). An incipient dynamical equilibrium is evident. A comparison of the calculated axis bending with that observed in corresponding crystal structure data is presented. Simulation of the DNA in complex with the protein and that of the free DNA in solution, starting from the crystal conformation, reveal the dynamical changes that occur on complex formation.     

Modeling H3 histone N-terminal tail and linker DNA interactions

Molecular dynamics computer simulations were performed for the 25-residue N-terminal tail of the H3 histone protein in the proximity of a DNA segment of 10 base pairs (bp), representing a model for the linker DNA in chromatin. Several least biased configurations were used as initial configurations. The secondary structure content of the protein was increased by the presence of DNA close to it, but the locations of the secondary motifs were different for different initial orientations of the DNA grooves with respect to the protein. As a common feature to all simulations, the electrostatic attraction between negatively charged DNA and positively charged protein was screened by the water solvent and counterbalanced by the intrinsic compaction of the protein due to hydrophobic effects. The protein secondary structure limited the covering of DNA by the protein to 4-5 bp. The degree of compaction and charge density of the bound protein suggests a possible role of H3 tail in a nonspecific bending and plasticity of the linker DNA when the protein is located in the crowded dense chromatin.     

Molecular dynamics simulation study of interaction between a class IIa bacteriocin and its immunity protein

Molecular dynamics (MD) simulations of carnobacteriocin B2 (CbnB2), a structurally well-characterized class IIa bacteriocin, and its immunity protein (ImB2) in lipid bilayer environment have been conducted to explore the interaction between them. Six 30-ns simulations were conducted in DPPC or POPG bilayer systems. In these simulations, ImB2 was placed in the aqueous layer with different orientations facing CbnB2 to sample all the faces of ImB2. The MD results indicate that (i) while CbnB2 remained embedded in the bilayer, it tends to move toward the interface, and (ii) the presence of CbnB2 in the DPPC bilayer attracts ImB2 toward the bilayer. In one of the orientations in DPPC bilayer system (simulation 1), ImB2 penetrates the bilayer and interacts with CbnB2 by ion-pair interaction. At several instances toward the later half of the simulation (15-30 ns), ImB2 and CbnB2 were found to form salt-bridge between Arg95 of ImB2 and Glu24 of CbnB2. Simulation in POPG bilayer displayed strong interaction between the positively charged ImB2 and the negatively charged polar head groups of the POPG molecules at the lipid-water interface. However, ImB2 was not able to penetrate the bilayer thereby preventing any interaction between ImB2 and CbnB2.     

Dynamics of water and ions around DNA: What is so special about them?

Water around biomolecules is special for behaving strangely – both in terms of structure and dynamics, while ions are found to control various interactions in biomolecules such as DNA, proteins and lipids. The questions that how water and ions around these biomolecules behave in terms of their structure and dynamics, and how they affect the biomolecular functions have triggered tremendous research activities worldwide. Such activities not only unfolded important static and dynamic properties of water and ions around these biomolecules, but also provoked heated debate regarding their explanation and role in biological functions. DNA, being negatively charged, interacts strongly with the surrounding dipolar water and positively charged counterions, leading to complex electrostatic coupling of water and ions with the DNA. Recent time-resolved fluorescence Stokes shift experiments and related computer simulation studies from our and other laboratories have unfolded some unique dynamic characteristics of water and ions near different structures of DNA. These results are discussed here to showcase the specialty of water and ion dynamics around DNA.     

A structure-based simulation approach for electron paramagnetic resonance spectra using molecular and stochastic dynamics simulations

Electron paramagnetic resonance (EPR) spectroscopy using site-directed spin-labeling is an appropriate technique to analyze the structure and dynamics of flexible protein regions as well as protein-protein interactions under native conditions. The analysis of a set of protein mutants with consecutive spin-label positions leads to the identification of secondary and tertiary structure elements. In the first place, continuous-wave EPR spectra reflect the motional freedom of the spin-label specifically linked to a desired site within the protein. EPR spectra calculations based on molecular dynamics (MD) and stochastic dynamics simulations facilitate verification or refinement of predicted computer-aided models of local protein conformations. The presented spectra simulation algorithm implies a specialized in vacuo MD simulation at 600 K with additional restrictions to sample the entire accessible space of the bound spin-label without large temporal effort. It is shown that the distribution of spin-label orientations obtained from such MD simulations at 600 K agrees well with the extrapolated motion behavior during a long timescale MD at 300 K with explicit water. The following potential-dependent stochastic dynamics simulation combines the MD data about the site-specific orientation probabilities of the spin-label with a realistic rotational diffusion coefficient yielding a set of trajectories, each more than 700 ns long, essential to calculate the EPR spectrum. Analyses of a structural model of the loop between helices E and F of bacteriorhodopsin are illustrated to demonstrate the applicability and potentials of the reported simulation approach. Furthermore, effects on the motional freedom of bound spin-labels induced by solubilization of bacteriorhodopsin with Triton X-100 are examined.     

Homology modeling based solution structure of Hoxc8-DNA complex: role of context bases outside TAAT stretch

The 3D structure of neither Hoxc8 nor Hoxc8-DNA complex is known. The repressor protein Hoxc8 binds to the TAAT stretch of the promoter of the osteopontin gene and modulates its expression. Over expression of the osteopontin gene is related to diseases like osteoporosis, multiple sclerosis, cancer et cetera. In this paper we have proposed a 3D structure of Hoxc8-DNA complex obtained by Homology modeling and molecular dynamics (MD) simulation in explicit water. The crystal structure (9ant.pdb) of Antennapedia homeodomain in complex with its DNA sequence was chosen as the template based on (i) high sequence identity (85% for the protein and 60% for the DNA) and (ii) the presence of the TAAT stretch in interaction with the protein. The resulting model was refined by MD simulation for 2.0ns in explicit water. This refined model was then characterized in terms of the structural and the interactional features to improve our understanding of the mechanism of Hoxc8-DNA recognition. The interaction pattern shows that the residues Ile(195), Gln(198), and Asn(199), and the bases S2-(4)TAATG(8) are most important for recognition suggesting the stretch TAATG as the "true recognition element" in the present case. A strong and long-lived water bridge connecting Gln(198) and the base of S1-C(7) complementary to S2-G(8) was observed. Our predicted model of Hoxc8-DNA complex provides us with features that are consistent with the available experimental data on Hoxc8 and the general features of other homeodomain-DNA complexes. The predictions based on the model are also amenable to experimental verification.     

Simulation of the substrate cavity dynamics of quercetinase

Molecular dynamics (MD) simulations have been performed on quercetin 2,3 dioxygenase (2,3QD) to study the mobility and flexibility of the substrate cavity. 2,3QD is the only firmly established Cu-containing dioxygenase known so far. It catalyses the breakage of the O-heterocycle of flavonols. The substrates occupy a shallow and overall hydrophobic cavity proximal to the metal centre of the homo-dimeric enzyme. The linker connecting the C-terminal and N-terminal domains in the monomer is partly disordered in the crystal structure and part of it forms a flexible lid at the entrance of the substrate cavity. This loop has been tentatively assigned a role in the enzyme mechanism: it helps lock the substrate into place. The dynamics of this loop has been investigated by MD simulation. The initial coordinates were taken from the crystal structure of 2,3QD in the presence of the substrate kaempferol (KMP). After equilibration and simulation over 7.2ns the substrate was removed and another equilibration and simulation of 7.2ns was performed. The results show that the structures of the free enzyme as well as of the enzyme-substrate complex are stable in MD simulation. The linker shows strongly enhanced mobility in the loop region that is close to the entrance to the substrate cavity (residues 154-169). Movement of the loop takes place on a timescale of 5-10ns. To confirm the conclusions about the loop dynamics drawn from the 7.2ns simulation, the simulation was extended with another 8ns. When substrate binds into the cavity the loop orders remarkably, although mobility is retained by residues 155-158. Some regions of the loop (residues 154-160 and 164-176) move over a considerable distance and approach the substrate closely, reinforcing the idea that they lock the substrate in the substrate cavity. The enthalpic component of the interaction of the loop with the protein and the KMP appears to favour the locking of the substrate. Two water molecules were found immobilised in the cavity, one of which exhibited rotation on the picosecond timescale. When the substrate is removed, the empty cavity fills up with water within 200ps.     

Conformational sampling by NMR solution structures calculated with the program DIANA evaluated by comparison with long-time molecular dynamics calculations in explicit water

The NMR solution structure of bovine pancreatic trypsin inhibitor (BPTI) obtained by distance geometry calculations with the program DIANA is compared with groups of conformers generated by molecular dynamics (MD) simulations in explicit water at ambient temperature and pressure. The MD simulations started from a single conformer and were free or restrained either by the experimental NOE distance restraints or by time-averaged restraints; the groups of conformers were collected either in 10 ps intervals during 200 ps periods of simulation, or in 50 ps intervals during a 1 ns period of simulation. Overall, these comparisons show that the standard protein structure determination protocol with the program DIANA provides a picture of the protein structure that is in agreement with MD simulations using “realistic” potential functions over a nanosecond timescale. For well-constrained molecular regions there is a trend in the free MD simulation of duration 1 ns that the sampling of the conformation space is slightly increased relative to the DIANA calculations. In contrast, for surface-exposed side-chains that are less extensively constrained by the NMR data, the DIANA conformers tend to sample larger regions of conformational space than conformers selected from any of the MD trajectories. Additional insights into the behavior of surface side-chains come from comparison of the MD runs of 200 ps or 1 ns duration. In this time range the sampling of conformation space by the protein surface depends strongly on the length of the simulation, which indicates that significant side-chain transitions occur on the nanosecond timescale and that much longer simulations will be needed to obtain statistically significant data on side-chain dynamics.     

Molecular dynamics study of water penetration in staphylococcal nuclease

The ionization properties of Lys and Glu residues buried in the hydrophobic core of staphylococcal nuclease (SN) suggest that the interior of this protein behaves as a highly polarizable medium with an apparent dielectric constant near 10. This has been rationalized previously in terms of localized conformational relaxation concomitant with the ionization of the internal residue, and with contributions by internal water molecules. Paradoxically, the crystal structure of the SN V66E variant shows internal water molecules and the structure of the V66K variant does not. To assess the structural and dynamical character of interior water molecules in SN, a series of 10-ns-long molecular dynamics (MD) simulations was performed with wild-type SN, and with the V66E and V66K variants with Glu66 and Lys66 in the neutral form. Internal water molecules were identified based on their coordination state and characterized in terms of their residence times, average location, dipole moment fluctuations, hydrogen bonding interactions, and interaction energies. The locations of the water molecules that have residence times of several nanoseconds and display small mean-square displacements agree well with the locations of crystallographically observed water molecules. Additional, relatively disordered water molecules that are not observed crystallographically were found in internal hydrophobic locations. All of the interior water molecules that were analyzed in detail displayed a distribution of interaction energies with higher mean value and narrower width than a bulk water molecule. This underscores the importance of protein dynamics for hydration of the protein interior. Further analysis of the MD trajectories revealed that the fluctuations in the protein structure (especially the loop elements) can strongly influence protein hydration by changing the patterns or strengths of hydrogen bonding interactions between water molecules and the protein. To investigate the dynamical response of the protein to burial of charged groups in the protein interior, MD simulations were performed with Glu66 and Lys66 in the charged state. Overall, the MD simulations suggest that a conformational change rather than internal water molecules is the dominant determinant of the high apparent polarizability of the protein interior.     

Structural insights into the anti-cancer activity of quercetin on G-tetrad,mixed G-tetrad,and G-quadruplex DNA using quantum chemical and molecular dynamics simulations

Abstract

Human telomerase referred as ‘terminal transferase’ is a nucleoprotein enzyme which inhibits the disintegration of telomere length and act as a drug target for the anticancer therapy. The tandem repeating structure of telomere sequence forms the guanine-rich quadruplex structures that stabilize stacked tetrads. In our present work, we have investigated the interaction of quercetin with DNA tetrads using DFT. Geometrical analysis revealed that the influence of quercetin drug induces the structural changes into the DNA tetrads. Among DNA tetrads, the quercetin stacked with GCGC tetrad has the highest interaction energy of ?88.08?kcal/mol. The binding mode and the structural stability are verified by the absorption spectroscopy method. The longer wavelength was found at 380?nm and it exhibits bathochromic shift. The findings help us to understand the binding nature of quercetin drug with DNA tetrads and it also inhibits the telomerase activity. Further, the quercetin drug interacted with G-quadruplex DNA by using molecular dynamics (MD) simulation studies for 100?ns simulation at different temperatures and different pH levels (T?=?298 K, 320?K and pH = 7.4, 5.4). The structural stability of the quercetin with G-quadruplex structure is confirmed by RMSD. For the acidic condition (pH = 5.4), the binding affinity is higher toward G-quadruplex DNA, this result resembles that the quercetin drug is well interacted with G-quadruplex DNA at acidic condition (pH = 7.4) than the neutral condition. The obtained results show that quercetin drug stabilizes the G-quadruplex DNA, which regulates telomerase enzyme and it potentially acts as a novel anti-cancer agent.

Communicated by Ramaswamy H. Sarma     

On the interactions of histone H4 and H4 peptides with DNA. Electrooptical, hydrodynamic and electron microscopy studies

The interactions of DNA with histone H4 and with its fragments N-H4 (1-84) and C-H4 (85-102) have been studied by using electrooptical techniques, viscosity and electron microscopy. Electron microscopy reveals that histone H4 induces a large folding of DNA molecules : this is in agreement with electrooptical measurements which indicate that, with the increase of their ratio, H4/DNA complexes undergo a gradual process of condensation. Viscosity measurements show that complexes at ratios up to 0.20-0.25 become more rigid as compared to DNA. It appears that C-H4, and not the N-H4 fragment, causes a great distorsion to the structure of DNA, accompanied by an increase of rigidity at ratios up to 0.20-0.25, as occurs for H4/DNA complexes. Electrooptical studies of C-H4/DNA complexes show, along a range of histone/DNA ratios, an important permanent dipole component. These effects reveal a particular mode of interaction of C-H4 with DNA, indicating that some charged residues of the peptide are kept distant enough from the DNA backbone. As no dipole character, in addition to that shown for DNA, has been detected for H4/DNA complexes, it is concluded that the conformation of the H4 molecule modifies to some extent the interaction of the C-terminal region. Our results show that this histone, and particularly its C-terminal region, is important as a determinant factor in the folding of DNA within artificial complexes.     

Nucleosome composition regulates the histone H3 tail conformational ensemble and accessibility

Hexasomes and tetrasomes are intermediates in nucleosome assembly and disassembly. Their formation is promoted by histone chaperones, ATP-dependent remodelers, and RNA polymerase II. In addition, hexasomes are maintained in transcribed genes and could be an important regulatory factor. While nucleosome composition has been shown to affect the structure and accessibility of DNA, its influence on histone tails is largely unknown. Here, we investigate the conformational dynamics of the H3 tail in the hexasome and tetrasome. Using a combination of NMR spectroscopy, MD simulations, and trypsin proteolysis, we find that the conformational ensemble of the H3 tail is regulated by nucleosome composition. As has been found for the nucleosome, the H3 tails bind robustly to DNA within the hexasome and tetrasome, but upon loss of the H2A/H2B dimer, we determined that the adjacent H3 tail has an altered conformational ensemble, increase in dynamics, and increase in accessibility. Similar to observations of DNA dynamics, this is seen to be asymmetric in the hexasome. Our results indicate that nucleosome composition has the potential to regulate chromatin signaling and ultimately help shape the chromatin landscape.     

Histone h3 lysine 4 methylation marks postreplicative human cytomegalovirus chromatin

In the nuclei of permissive cells, human cytomegalovirus genomes form nucleosomal structures initially resembling heterochromatin but gradually switching to a euchromatin-like state. This switch is characterized by a decrease in histone H3 K9 methylation and a marked increase in H3 tail acetylation and H3 K4 methylation across the viral genome. We used ganciclovir and a mutant virus encoding a reversibly destabilized DNA polymerase to examine the impact of DNA replication on histone modification dynamics at the viral chromatin. The changes in H3 tail acetylation and H3 K9 methylation proceeded in a DNA replication-independent fashion. In contrast, the increase in H3 K4 methylation proved to depend widely on viral DNA synthesis. Consistently, labeling of nascent DNA using "click chemistry" revealed preferential incorporation of methylated H3 K4 into viral (but not cellular) chromatin during or following DNA replication. This study demonstrates largely selective epigenetic tagging of postreplicative human cytomegalovirus chromatin.     

Conformational states of a TT mismatch from molecular dynamics simulation of duplex d (CGCGATTCGCG).

The TT mismatch region in duplex d (CGCGATTCGCG) was studied using a 500-ps molecular dynamics (MD) simulation in water, and a series of 1-ps MD simulations and energy minimizations in vacuum. The DNA maintained its duplex structure, although the mismatch region showed significantly higher flexibility than the GC regions. The predominant conformation in the 500-ps MD simulation involved an average -42 degrees propeller twist between T6 and T'6, and a -22 degree buckle between A5 and T'7. One hydrogen bond was formed between T6 and T'6, and another between T6 and the O2 of T'7, with both Watson-Crick hydrogen bonds between A5 and T'7 remaining intact. The minimizations resulted in conformations with the equivalent hydrogen-bonding pattern, as well as ones with "wobble pair" hydrogen bonds between T6 and T'6. However, the wobble pair conformation was found to be unstable in the water simulation.     

Molecular dynamics study of a gelsolin-derived peptide binding to a lipid bilayer containing phosphatidylinositol 4,5-bisphosphate

Gelsolin is an actin-severing protein whose action is initiated by Ca(2+) and inhibited by binding to phosphorylated inositol lipid or phosphoinositides. The regions of gelsolin responsible for phosphoinositide binding are comprised of residues 150-169 (G150-169) and 135-142 (G135-142). The corresponding peptides possess similar binding potency as native gelsolin. Their common feature is the presence of arginine and lysine residues that can bind to negatively charged phosphate groups of phosphoinositides. In this work the binding of the G150-169 peptide to a phosphatidylinositol 4,5-bisphosphate (PIP2) cluster in a lipid membrane model was investigated by molecular dynamics calculations (MD) with the AMBER 4.1 force field, taking into account explicit solvent molecules. Initially the structure of G150-169 was simulated by using the electrostatically driven Monte Carlo (EDMC) and MD methods, and the resulting structure agreed within 3.7 A backbone-atom root mean square deviation with the corresponding experimentally derived structure (PDB code: 1SOL). Using this model for the peptide, a subsequent MD simulation of G150-169 in a periodic box containing a model of dimyristoyl-phosphatidylcholine (DMPC) lipids with a cluster of four PIP2 molecules was carried out. During the simulation G150-169 interacted strongly with PIP2 molecules, initially by formation of salt bridges between its N-terminal basic groups and the phosphate groups of PIP2, followed by formation of hydrophobic bonds between the hydrophobic side chains of the peptide and the fatty acid tail of the lipid. As a result of the formation of hydrophobic bonds, the PIP2 molecules were pulled out from the lipid bilayer. This mode of binding differs from those of other PIP2-binding protein motifs such as PH domains that interact solely with the hydrophilic head group of PIP2. These results suggest that dissociation of gelsolin from actin by PIP2 lipids may involve entering of the PIP2 molecules to the gelsolin-actin interface, thereby weakening the interactions between these proteins.