Surviving an abdominal aortic aneurysm type II endoleak rupture without surgical intervention

We describe a patient who survived a ruptured abdominal aortic aneurysm without any surgical intervention. The patient had previously had endovascular repair of the aneurysm and surveillance of a stable persistent type II endoleak. This case highlights the difficulties surrounding type II endoleak, its natural history, and the ongoing controversies of its management.

All-atom empirical potential for molecular modeling and dynamics studies of proteins

New protein parameters are reported for the all-atom empirical energy function in the CHARMM program. The parameter evaluation was based on a self-consistent approach designed to achieve a balance between the internal (bonding) and interaction (nonbonding) terms of the force field and among the solvent-solvent, solvent-solute, and solute-solute interactions. Optimization of the internal parameters used experimental gas-phase geometries, vibrational spectra, and torsional energy surfaces supplemented with ab initio results. The peptide backbone bonding parameters were optimized with respect to data for N-methylacetamide and the alanine dipeptide. The interaction parameters, particularly the atomic charges, were determined by fitting ab initio interaction energies and geometries of complexes between water and model compounds that represented the backbone and the various side chains. In addition, dipole moments, experimental heats and free energies of vaporization, solvation and sublimation, molecular volumes, and crystal pressures and structures were used in the optimization. The resulting protein parameters were tested by applying them to noncyclic tripeptide crystals, cyclic peptide crystals, and the proteins crambin, bovine pancreatic trypsin inhibitor, and carbonmonoxy myoglobin in vacuo and in crystals. A detailed analysis of the relationship between the alanine dipeptide potential energy surface and calculated protein φ, χ angles was made and used in optimizing the peptide group torsional parameters. The results demonstrate that use of ab initio structural and energetic data by themselves are not sufficient to obtain an adequate backbone representation for peptides and proteins in solution and in crystals. Extensive comparisons between molecular dynamics simulations and experimental data for polypeptides and proteins were performed for both structural and dynamic properties. Energy minimization and dynamics simulations for crystals demonstrate that the latter are needed to obtain meaningful comparisons with experimental crystal structures. The presented parameters, in combination with the previously published CHARMM all-atom parameters for nucleic acids and lipids, provide a consistent set for condensed-phase simulations of a wide variety of molecules of biological interest.

Accelerating Virtual High-Throughput Ligand Docking: current technology and case study on a petascale supercomputer

In this paper we give the current state of high-throughput virtual screening. We describe a case study of using a task-parallel MPI (Message Passing Interface) version of Autodock4 [1], [2] to run a virtual high-throughput screen of one-million compounds on the Jaguar Cray XK6 Supercomputer at Oak Ridge National Laboratory. We include a description of scripts developed to increase the efficiency of the predocking file preparation and postdocking analysis. A detailed tutorial, scripts, and source code for this MPI version of Autodock4 are available online at http://www.bio.utk.edu/baudrylab/autodockmpi.htm.

Hydrolysis of DFP and the nerve agent (S)-sarin by DFPase proceeds along two different reaction pathways: implications for engineering bioscavengers

Organophosphorus (OP) nerve agents such as (S)-sarin are among the most highly toxic compounds that have been synthesized. Engineering enzymes that catalyze the hydrolysis of nerve agents (“bioscavengers”) is an emerging prophylactic approach to diminish their toxic effects. Although its native function is not known, diisopropyl fluorophosphatase (DFPase) from Loligo vulgaris catalyzes the hydrolysis of OP compounds. Here, we investigate the mechanisms of diisopropylfluorophosphate (DFP) and (S)-sarin hydrolysis by DFPase with quantum mechanical/molecular mechanical umbrella sampling simulations. We find that the mechanism for hydrolysis of DFP involves nucleophilic attack by Asp229 on phosphorus to form a pentavalent intermediate. P–F bond dissociation then yields a phosphoacyl enzyme intermediate in the rate-limiting step. The simulations suggest that a water molecule, coordinated to the catalytic Ca²⁺, donates a proton to Asp121 and then attacks the tetrahedral phosphoacyl intermediate to liberate the diisopropylphosphate product. In contrast, the calculated free energy barrier for hydrolysis of (S)-sarin by the same mechanism is highly unfavorable, primarily because of the instability of the pentavalent phosphoenzyme species. Instead, simulations suggest that hydrolysis of (S)-sarin proceeds by a mechanism in which Asp229 could activate an intervening water molecule for nucleophilic attack on the substrate. These findings may lead to improved strategies for engineering DFPase and related six-bladed β-propeller folds for more efficient degradation of OP compounds.

de Gennes narrowing describes the relative motion of protein domains

The relative motion of structural domains is essential for the biological function of many proteins. Here, by analyzing neutron scattering data and performing molecular dynamics simulations, we find that interdomain motion in several proteins obeys the principle of de Gennes narrowing, in which the wave vector dependence of the interdomain diffusion coefficient is inversely proportional to the interdomain structure factor. Thus, the rate of interdomain motion is inversely proportional to the probability distribution of the spatial configurations of domains.

Mercury methylation by HgcA: theory supports carbanion transfer to Hg(II)

Many proteins use corrinoid cofactors to facilitate methyl transfer reactions. Recently, a corrinoid protein, HgcA, has been shown to be required for the production of the neurotoxin methylmercury by anaerobic bacteria. A strictly conserved Cys residue in HgcA was predicted to be a lower-axial ligand to Co(III), which has never been observed in a corrinoid protein. Here, we use density functional theory to study homolytic and heterolytic Co-C bond dissociation and methyl transfer to Hg(II) substrates with model methylcobalamin complexes containing a lower-axial Cys or His ligand to cobalt, the latter of which is commonly found in other corrinoid proteins. We find that Cys thiolate coordination to Co facilitates both methyl radical and methyl carbanion transfer to Hg(II) substrates, but carbanion transfer is more favorable overall in the condensed phase. Thus, our findings are consistent with HgcA representing a new class of corrinoid protein capable of transferring methyl groups to electrophilic substrates.

Impact of resistance mutations on inhibitor binding to HIV-1 integrase

HIV-1 integrase (IN) is essential for HIV-1 replication, catalyzing two key reaction steps termed 3' processing and strand transfer. Therefore, IN has become an important target for antiviral drug discovery. However, mutants have emerged, such as E92Q/N155H and G140S/Q148H, which confer resistance to raltegravir (RAL), the first IN strand transfer inhibitor (INSTI) approved by the FDA, and to the recently approved elvitegravir (EVG). To gain insights into the molecular mechanisms of ligand binding and drug resistance, we performed molecular dynamics (MD) simulations of homology models of the HIV-1 IN and four relevant mutants complexed with viral DNA and RAL. The results show that the structure and dynamics of the 140s' loop, comprising residues 140 to 149, are strongly influenced by the IN mutations. In the simulation of the G140S/Q148H double mutant, we observe spontaneous dissociation of RAL from the active site, followed by an intrahelical swing-back of the 3'-OH group of nucleotide A17, consistent with the experimental observation that the G140S/Q148H mutant exhibits the highest resistance to RAL compared to other IN mutants. An important hydrogen bond between residues 145 and 148 is present in the wild-type IN but not in the G140S/Q148H mutant, accounting for the structural and dynamical differences of the 140s' loop and ultimately impairing RAL binding in the double mutant. End-point free energy calculations that broadly capture the experimentally known RAL binding profiles elucidate the contributions of the 140s' loop to RAL binding free energies and suggest possible approaches to overcoming drug resistance.

Long-time mean-square displacements in proteins

We propose a method for obtaining the intrinsic, long-time mean square displacement (MSD) of atoms and molecules in proteins from finite-time molecular dynamics (MD) simulations. Typical data from simulations are limited to times of 1 to 10 ns, and over this time period the calculated MSD continues to increase without a clear limiting value. The proposed method consists of fitting a model to MD simulation-derived values of the incoherent intermediate neutron scattering function, I(inc)(Q,t), for finite times. The infinite-time MSD, <r(2)>, appears as a parameter in the model and is determined by fits of the model to the finite-time I(inc)(Q,t). Specifically, the <r(2)> is defined in the usual way in terms of the Debye-Waller factor as I(Q,t=∞)=exp(-Q(2)<r(2)>/3). The method is illustrated by obtaining the intrinsic MSD <r(2)> of hydrated lysozyme powder (h=0.4 g water/g protein) over a wide temperature range. The intrinsic <r(2)> obtained from data out to 1 and to 10 ns is found to be the same. The intrinsic <r(2)> is approximately twice the value of the MSD that is reached in simulations after times of 1 ns which correspond to those observed using neutron instruments that have an energy resolution width of 1 μeV.

Soft collective fluctuations governing hydrophobic association

The interaction between two associating hydrophobic particles has traditionally been explained in terms of the release of entropically frustrated hydration shell water molecules. However, this picture cannot account for the kinetics of hydrophobic association and is therefore not capable of providing a microscopic description of the hydrophobic interaction (HI). Here, Monte Carlo simulations of a pair of molecular-scale apolar solutes in aqueous solution reveal the critical role of collective fluctuations in the hydrogen bond (HB) network for the microscopic picture of the HI. The main contribution to the HI is the relaxation of solute-water translational correlations. The existence of a heat capacity maximum at the desolvation barrier is shown to arise from softening of non-HB water fluctuations and the relaxation of many-body correlations in the labile HB network. The microscopic event governing the kinetics of hydrophobic association has turned out to be a relatively large critical collective fluctuation in hydration water displacing a substantial fraction of HB clusters from the inner to the outer region of the first hydration shell.

Zaccai neutron resilience and site-specific hydration dynamics in a globular protein

A discussion is presented of contributions of the Zaccai group to the understanding of flexibility in biological macromolecules using dynamic neutron scattering. The concept of resilience as introduced by Zaccai is discussed and investigated using molecular dynamics simulation on camphor-bound cytochrome P450. The resilience of hydrophilic residues is found to be more strongly affected by hydration than that of hydrophobic counterparts. The hydration-induced softening of protein propagates from the surface into the dry core. Moreover, buried hydrophilic residues behave more like those exposed on the protein surface, and are different from their hydrophobic counterparts.

GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit

MOTIVATION

Molecular simulation has historically been a low-throughput technique, but faster computers and increasing amounts of genomic and structural data are changing this by enabling large-scale automated simulation of, for instance, many conformers or mutants of biomolecules with or without a range of ligands. At the same time, advances in performance and scaling now make it possible to model complex biomolecular interaction and function in a manner directly testable by experiment. These applications share a need for fast and efficient software that can be deployed on massive scale in clusters, web servers, distributed computing or cloud resources.

RESULTS

Here, we present a range of new simulation algorithms and features developed during the past 4 years, leading up to the GROMACS 4.5 software package. The software now automatically handles wide classes of biomolecules, such as proteins, nucleic acids and lipids, and comes with all commonly used force fields for these molecules built-in. GROMACS supports several implicit solvent models, as well as new free-energy algorithms, and the software now uses multithreading for efficient parallelization even on low-end systems, including windows-based workstations. Together with hand-tuned assembly kernels and state-of-the-art parallelization, this provides extremely high performance and cost efficiency for high-throughput as well as massively parallel simulations.

AVAILABILITY

GROMACS is an open source and free software available from http://www.gromacs.org.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Three entropic classes of side chain in a globular protein

The relationship between the NMR methyl group axial order parameter and the side chain conformational entropy is investigated in inhibitor-bound and apo human HIV protease using molecular dynamics simulation. Three distinct entropic classes of methyl-bearing side chains, determined by the topological distance of the methyl group from the protein backbone (i.e., the number of χ-bonds between the Cα and the carbon of the CH3 group), are revealed by atomistic trajectory analyses performed in the local frame of reference of individual methyl probes. The results demonstrate that topologically equivalent methyl groups experience similar nonbonded microenvironments regardless of the type of residues to which they are attached. Similarly, methyl groups that belong to the same side chain but that are not topologically equivalent exhibit different thermodynamic and dynamic properties. The two-parameter classification (based upon entropy and methyl axial order parameter) of side chains described here permits improved estimates of the conformational entropies of proteins from NMR motional parameters.

Reorientation and dimerization of the membrane-bound antimicrobial peptide PGLa from microsecond all-atom MD simulations

The membrane-active antimicrobial peptide PGLa from Xenopus laevis is known from solid-state (2)H-, (15)N-, and (19)F-NMR spectroscopy to occupy two distinct α-helical surface adsorbed states in membranes: a surface-bound S-state with a tilt angle of ~95° at low peptide/lipid molar ratio (P/L = 1:200), and an obliquely tilted T-state with a tilt angle of 127° at higher peptide concentration (P/L = 1:50). Using a rapid molecular-dynamics insertion protocol in combination with microsecond-scale simulation, we have characterized the structure of both states in detail. As expected, the amphiphilic peptide resides horizontally on the membrane surface in a monomeric form at a low P/L, whereas the T-state is seen in the simulations to be a symmetric antiparallel dimer, with close contacts between small glycine and alanine residues at the interface. The computed tilt angles and azimuthal rotations, as well as the quadrupolar splittings predicted from the simulations agree with the experimental NMR data. The simulations reveal many structural details previously inaccessible, such as the immersion depth of the peptide in the membrane and the packing of the dimerization interface. The study highlights the ability and limitations of current state-of-the-art multimicrosecond all-atom simulations of membrane-active peptides to complement experimental data from solid-state NMR.

Elastic and conformational softness of a globular protein

Flexibility, or softness, is crucial for protein function and consists of a conformational component, involving jumps between potential wells, and an elastic component, involving fluctuations within the wells. Combining molecular dynamics simulation with incoherent neutron scattering and light scattering measurements on green fluorescent protein, we reveal a relationship between the intrawell fluctuations and elastic moduli of the protein. This finding leads to a simple means of experimentally separating the conformational from the elastic atomic displacements.

Time-dependent density functional theory assessment of UV absorption of benzoic acid derivatives

Benzoic acid (BA) derivatives of environmental relevance exhibit various photophysical and photochemical characteristics. Here, time-dependent density functional theory (TDDFT) is used to calculate photoexcitations of eight selected BAs and the results are compared with UV spectra determined experimentally. High-level gas-phase EOM-CCSD calculations and experimental aqueous-phase spectra were used as the references for the gas-phase and aqueous-phase TDDFT results, respectively. A cluster-continuum model was used in the aqueous-phase calculations. Among the 15 exchange-correlation (XC) functionals assessed, five functionals, including the meta-GGA hybrid M06-2X, double hybrid B2PLYPD, and range-separated functionals CAM-B3LYP, ωB97XD, and LC-ωPBE, were found to be in excellent agreement with the EOM-CCSD gas-phase calculations. These functionals furnished excitation energies consistent with the pH dependence of the experimental spectra with a standard deviation (STDEV) of ∼0.20 eV. A molecular orbital analysis revealed a πσ* feature of the low-lying transitions of the BAs. The CAM-B3LYP functional showed the best overall performance and therefore shows promise for TDDFT calculations of processes involving photoexcitations of benzoic acid derivatives.

Coupled flexibility change in cytochrome P450cam substrate binding determined by neutron scattering, NMR, and molecular dynamics simulation

Neutron scattering and nuclear magnetic resonance relaxation experiments are combined with molecular dynamics (MD) simulations in a novel, to our knowledge, approach to investigate the change in internal dynamics on substrate (camphor) binding to a protein (cytochrome P450cam). The MD simulations agree well with both the neutron scattering, which furnishes information on global flexibility, and the nuclear magnetic resonance data, which provides residue-specific order parameters. Decreased fluctuations are seen in the camphor-bound form using all three techniques, dominated by changes in specific regions of the protein. The combined experimental and simulation results permit a detailed description of the dynamical change, which involves modifications in the coupling between the dominant regions and concomitant substrate access channel closing, via specific salt-bridge, hydrogen-bonding, and hydrophobic interactions. The work demonstrates how the combination of complementary experimental spectroscopies with MD simulation can provide an in-depth description of functional dynamical protein changes.

Alteration of water structure by peptide clusters revealed by neutron scattering in the small-angle region (below 1 Å(-1))

Solution scattering of neutrons and x-rays can provide direct information on local interactions of importance for biomolecular folding and structure. Here, neutron scattering experiments are combined with molecular-dynamics simulation to interpret the scattering signal of a series of dipeptides with varying degrees of hydrophobicity (GlyAla, GlyPro, and AlaPro) in concentrated aqueous solution (1:20 solute/water ratio) in which the peptides form large segregates (up to 50-60 amino acids). Two main results are found: 1), the shift to lower Q of the so-called water-ring peak (Q ≈ 2 Å(-1)) arises mainly from an overlap of water-peptide and peptide-peptide correlations in the region of 1.3 <Q< 2 Å(-1), rather than from a shift of the water signal induced by the presence of the clusters; and 2), in the low-Q region (Q ≈ 0.6 Å(-1)) a positive peak is observed originating from both the solute-solute correlations and changes in the water structure induced by the formation of the clusters. In particular, the water molecules are found to be more connected than in the bulk with hydrogen-bonding directions tangential to the exposed hydrophobic surfaces, and this effect increases with increasing peptide hydrophobicity. This work demonstrates that important information on the (hydrophobic) hydration of biomolecules can be obtained in the very-small-angle region.

Dynamics of protein and its hydration water: neutron scattering studies on fully deuterated GFP

We present a detailed analysis of the picosecond-to-nanosecond motions of green fluorescent protein (GFP) and its hydration water using neutron scattering spectroscopy and hydrogen/deuterium contrast. The analysis reveals that hydration water suppresses protein motions at lower temperatures (<~ 200 K), and facilitates protein dynamics at high temperatures. Experimental data demonstrate that the hydration water is harmonic at temperatures <~ 180-190 K and is not affected by the proteins' methyl group rotations. The dynamics of the hydration water exhibits changes at ~ 180-190 K that we ascribe to the glass transition in the hydrated protein. Our results confirm significant differences in the dynamics of protein and its hydration water at high temperatures: on the picosecond-to-nanosecond timescale, the hydration water exhibits diffusive dynamics, while the protein motions are localized to <~3 Å. The diffusion of the GFP hydration water is similar to the behavior of hydration water previously observed for other proteins. Comparison with other globular proteins (e.g., lysozyme) reveals that on the timescale of 1 ns and at equivalent hydration level, GFP dynamics (mean-square displacements and quasielastic intensity) are of much smaller amplitude. Moreover, the suppression of the protein dynamics by the hydration water at low temperatures appears to be stronger in GFP than in other globular proteins. We ascribe this observation to the barrellike structure of GFP.

Down-regulation of the caffeic acid O-methyltransferase gene in switchgrass reveals a novel monolignol analog

BACKGROUND

Down-regulation of the caffeic acid 3-O-methyltransferase EC 2.1.1.68 (COMT) gene in the lignin biosynthetic pathway of switchgrass (Panicum virgatum) resulted in cell walls of transgenic plants releasing more constituent sugars after pretreatment by dilute acid and treatment with glycosyl hydrolases from an added enzyme preparation and from Clostridium thermocellum. Fermentation of both wild-type and transgenic switchgrass after milder hot water pretreatment with no water washing showed that only the transgenic switchgrass inhibited C. thermocellum. Gas chromatography-mass spectrometry (GCMS)-based metabolomics were undertaken on cell wall aqueous extracts to determine the nature of the microbial inhibitors.

RESULTS

GCMS confirmed the increased concentration of a number of phenolic acids and aldehydes that are known inhibitors of microbial fermentation. Metabolomic analyses of the transgenic biomass additionally revealed the presence of a novel monolignol-like metabolite, identified as trans-3, 4-dimethoxy-5-hydroxycinnamyl alcohol (iso-sinapyl alcohol) in both non-pretreated, as well as hot water pretreated samples. iso-Sinapyl alcohol and its glucoside were subsequently generated by organic synthesis and the identity of natural and synthetic materials were confirmed by mass spectrometric and NMR analyses. The additional novel presence of iso-sinapic acid, iso-sinapyl aldehyde, and iso-syringin suggest the increased activity of a para-methyltransferase, concomitant with the reduced COMT activity, a strict meta-methyltransferase. Quantum chemical calculations were used to predict the most likely homodimeric lignans generated from dehydration reactions, but these products were not evident in plant samples.

CONCLUSIONS

Down-regulation of COMT activity in switchgrass resulted in the accumulation of previously undetected metabolites resembling sinapyl alcohol and its related metabolites, but that are derived from para-methylation of 5-hydroxyconiferyl alcohol, and related precursors and products; the accumulation of which suggests altered metabolism of 5-hydroxyconiferyl alcohol in switchgrass. Given that there was no indication that iso-sinapyl alcohol was integrated in cell walls, it is considered a monolignol analog. Diversion of substrates from sinapyl alcohol to free iso-sinapyl alcohol, its glucoside, and associated upstream lignin pathway changes, including increased phenolic aldehydes and acids, are together associated with more facile cell wall deconstruction, and to the observed inhibitory effect on microbial growth. However, iso-sinapyl alcohol and iso-sinapic acid, added separately to media, were not inhibitory to C. thermocellum cultures.

Surface hydration amplifies single-well protein atom diffusion propagating into the macromolecular core

The effect of surface hydration water on internal protein motion is of fundamental interest in molecular biophysics. Here, by decomposing the picosecond to nanosecond atomic motion in molecular dynamics simulations of lysozyme at different hydration levels into three components--localized single-well diffusion, methyl group rotation, and nonmethyl jumps--we show that the effect of surface hydration is mainly to increase the volume of the localized single-well diffusion. These diffusive motions are coupled in such a way that the hydration effect propagates from the protein surface into the dry core.

Derivation of mean-square displacements for protein dynamics from elastic incoherent neutron scattering

The derivation of mean-square displacements from elastic incoherent neutron scattering (EINS) of proteins is examined, with the aid of experiments on camphor-bound cytochrome P450cam and complementary molecular dynamics simulations. It is shown that a q(4) correction to the elastic incoherent structure factor (where q is the scattering vector) can be simply used to reliably estimate from the experiment both the average mean-square atomic displacement, <Δr(2)> of the nonexchanged hydrogen atoms in the protein and its variance, σ(2). The molecular dynamics simulation results are in broad agreement with the experimentally derived <Δr(2)> and σ(2) derived from EINS on instruments at two different energy resolutions, corresponding to dynamics on the ∼100 ps and ∼1 ns time scales. Significant dynamical heterogeneity is found to arise from methyl-group rotations. The easy-to-apply q(4) correction extends the information extracted from elastic incoherent neutron scattering experiments and should be of wide applicability.

Reconstruction of protein side-chain conformational free energy surfaces from NMR-derived methyl axis order parameters

An analytical approach is developed for reconstructing site-specific methyl-bearing protein side-chain conformational energy surfaces from NMR methyl axis order parameters (O(axis)(2)). Application of an enhanced sampling algorithm (adaptive biasing force) to molecular dynamics simulation of a protein, calcium-bound calmodulin, reveals a nonlinear correlation between O(axis)(2) and the populations of rotamer states of protein side-chains, permitting the rotamer populations to be extracted directly from O(axis)(2). The analytical approach yields side-chain conformational distributions that are in excellent agreement with those obtained from the enhanced-sampling MD results.

Unwrapping of nucleosomal DNA ends: a multiscale molecular dynamics study

To permit access to DNA-binding proteins involved in the control and expression of the genome, the nucleosome undergoes structural remodeling including unwrapping of nucleosomal DNA segments from the nucleosome core. Here we examine the mechanism of DNA dissociation from the nucleosome using microsecond timescale coarse-grained molecular dynamics simulations. The simulations exhibit short-lived, reversible DNA detachments from the nucleosome and long-lived DNA detachments not reversible on the timescale of the simulation. During the short-lived DNA detachments, 9 bp dissociate at one extremity of the nucleosome core and the H3 tail occupies the space freed by the detached DNA. The long-lived DNA detachments are characterized by structural rearrangements of the H3 tail including the formation of a turn-like structure at the base of the tail that sterically impedes the rewrapping of DNA on the nucleosome surface. Removal of the H3 tails causes the long-lived detachments to disappear. The physical consistency of the CG long-lived open state was verified by mapping a CG structure representative of this state back to atomic resolution and performing molecular dynamics as well as by comparing conformation-dependent free energies. Our results suggest that the H3 tail may stabilize the nucleosome in the open state during the initial stages of the nucleosome remodeling process.

Increase of both order and disorder in the first hydration shell with increasing solute polarity

Monte Carlo simulations of a small model solute in an aqueous solution are used to examine the effects of solute polarity on hydration structure. A judicious definition of the orientational order parameter leads to reinterpretation of the conventional picture of hydration. As the solute varies from hydrophobic to hydrophilic the ordered first shell water simultaneously fractionates into a more highly ordered and a more disordered component. The hydrogen-bond network rearranges such that the more ordered component relaxes to configurations of optimal intermolecular angles, the other fraction being released from the network.

Role of histone tails in structural stability of the nucleosome

Histone tails play an important role in nucleosome structure and dynamics. Here we investigate the effect of truncation of histone tails H3, H4, H2A and H2B on nucleosome structure with 100 ns all-atom molecular dynamics simulations. Tail domains of H3 and H2B show propensity of -helics formation during the intact nucleosome simulation. On truncation of H4 or H2B tails no structural change occurs in histones. However, H3 or H2A tail truncation results in structural alterations in the histone core domain, and in both the cases the structural change occurs in the H2A3 domain. We also find that the contacts between the histone H2A C terminal docking domain and surrounding residues are destabilized upon H3 tail truncation. The relation between the present observations and corresponding experiments is discussed.

Histone tails are the most common sites of post-translational modifications. Tail modifications alter both inter and intra nucleosomal interactions to disrupt the condensed chromatin structure, thereby playing crucial role in gene access. Here we investigated histone tail functions on the stability of a single nucleosome in atomic detail by selectively truncating tail domains in molecular dynamics simulations. Our study revealed that truncation of H3 or H2A tail results in structural alterations in the nucleosome core whereas truncation of H4 or H2B tail does not. A potential role of H2A C terminal tail in regulating nucleosome stability is discussed. Finally, an -helical domain formation was observed in one of the H3 tails and, upon truncation of this tail, structural changes occurred in closely lying histone domains. The correlation between tail-truncation and structural changes likely sheds light on allosteric regulation of nucleosome stability.