Conformational and Thermodynamic Differences Underlying Wild-Type and Mutant Eleven-Nineteen-Leukemia YEATS Domain Specificity for Epigenetic Marks.

Histone post-translational modifications (PTMs) are interpreted by multiple reader domains and proteins to regulate gene expression. The eleven-nineteen-leukemia (ENL) YEATS domain is a prototypical PTM reader that recognizes multiple lysine acetylation marks on the histone H3 tails as a way of recruiting chromatin remodellers. Two ENL YEATS mutations have been identified which have been linked with leukemia, Wilms tumor, and other forms of cancer and result in either an insertion or deletion of residues in the loop connecting beta sheets distant from the protein active site. In vitro experiments have shown that these mutations modulate the selectivities of YEATS domains for various lysine acetylation marks, although different experiments have provided contrasting views on the abilities of the insertion and deletion mutants to discern specific PTMs. Here, we have performed multiple molecular dynamics simulations of wild-type and insertion and deletion mutant YEATS domains free from and in complex with two PTM peptides: one that is acetylated at K9 of H3 and the other that is acetylated at residue K27 of H3. Results show that these two peptides have distinct flexibilities and binding energetics when bound to YEATS domains and that these properties are affected by interactions with residues within and outside of the peptide consensus motif. Furthermore, these properties are modulated by the YEATS insertion and deletion mutants, which results in disparate binding effects in these systems. Together, these results suggest that only the partial exposure of histone tails is sufficient in the context of nucleosomes for YEATS-mediated recognition of acetylation marks on histone tails. They also caution against the overinterpretation of results obtained from experiments on reader domain-histone peptide binding in isolation and not in the full-length nucleosome context.

Peptide Dynamics and Metadynamics: Leveraging Enhanced Sampling Molecular Dynamics to Robustly Model Long-Timescale Transitions.

Molecular dynamics simulations can in theory reveal the thermodynamics and kinetics of peptide conformational transitions at atomic-level resolution. However, even with modern computing power, they are limited in the timescales they can sample, which is especially problematic for peptides that are fully or partially disordered. Here, we discuss how the enhanced sampling methods accelerated molecular dynamics (aMD) and metadynamics can be leveraged in a complementary fashion to quickly explore conformational space and then robustly quantify the underlying free energy landscape. We apply these methods to two peptides that have an intrinsically disordered nature, the histone H3 and H4 N-terminal tails, and use metadynamics to compute the free energy landscape along collective variables discerned from aMD simulations. Results show that these peptides are largely disordered, with a slight preference for α-helical structures.

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.

Carbamylation promotes amyloidogenesis and induces structural changes in Tau-core hexapeptide fibrils.

BACKGROUND: Carbamylation is a non-enzymatic post-translational modification (PTM), which involves the covalent modification of N-terminus of protein or ε-amino group of Lys. The role of carbamylation in several age-related disorders is well documented, however, the relationship between carbamylation and neurodegenerative disorders including Alzheimer's disease remains uncharted. METHODS: In the present study, using aggregation-prone tau-core hexapeptide fragments 306VQIVYK311 (PHF6) and 275VQIINK280 (PHF6*) as models, we have elucidated the effect of carbamylation on aggregation kinetics and the changes occurring in the 3-dimensional architecture of fibrils using biophysical assays and molecular dynamics simulations. RESULTS: We found that carbamylation aids in amyloid formation and can convert the unstructured off-pathway aggregates into robust amyloids, which were toxic to cells. Electron microscopy images and molecular dynamics simulations of PHF6 fibrils showed that carbamylated peptides can form excess hydrogen bonds and modulate the pitch length and twist of peptides fibrils. We have also compared N-terminal carbamylation to acetylation and further extended our finding to full length tau that exhibits aggregation upon carbamylation even in the absence of any external inducer. CONCLUSION: Our in vitro and in silico results together suggest that carbamylation can modulate the aggregation pathway of the amyloidegenic sequences and cause structural changes in fibril assemblies. GENERAL SIGNIFICANCE: Carbamylation acts as a switch, which triggers the aggregation in short amyloidogenic peptide fragments and modulate the structural changes in resulting amyloid fibrils.

Preferential binding of fullerene and fullerenol with the N-terminal and middle regions of amyloid beta peptide: an in silico investigation.

Amyloid beta (Aß) deposits are implicated in the pathogenesis of debilitating neurodegenerative disorders such as Alzheimer's disease. In the present study, the interactions of carbon-based nanoparticles (NPs) such as fullerene and fullerenol having different surface chemistry with Aß were investigated using molecular dynamics simulations and docking studies. A detailed analysis of docking results showed that in 68% of the Aß conformations, fullerene and fullerenol showed interactions with the N-terminal region of the peptide. However, the high-affinity binding site (E=-48.31 kJ/mol) of fullerene resides in the hydrophobic middle region of the peptide, whereas fullerenol interacts favorably with the charged N-terminal region with a binding energy of -50.42 kJ/mol. The above differences in binding could be attributed to the surface chemistry of fullerene and fullerenol. Moreover, the N-terminal and middle regions of Aß play an important role in Aß aggregation. Therefore, the binding of fullerene and fullerenol could inhibit amyloid aggregation. This information will be helpful in designing NPs for targeting amyloid-related disorders.

Pushing the Limits of Structure-Based Models: Prediction of Nonglobular Protein Folding and Fibrils Formation with Go-Model Simulations.

The development of computational efficient models is essential to obtain a detailed characterization of the mechanisms underlying the folding of proteins and the formation of amyloid fibrils. Structure-based computational models (Go-model) with Cα or all-atom resolutions have been able to successfully delineate the mechanisms of folding of several globular proteins and offer an interesting alternative to computationally intensive simulations with explicit solvent description. Here, we explore the limits of Go-model predictions by analyzing the folding of the nonglobular repeat domain proteins Notch Ankyrin and p16INK4 and the formation of human islet amyloid polypeptide (hIAPP) fibrils. Folding trajectories of the repeat domain proteins revealed that an all-atom resolution is required to capture the folding pathways and cooperativity reported in experimental studies. The all-atom Go-model was also successful in predicting the free-energy landscape of hIAPP fibrillation, suggesting a "dock and lock" mechanism of fibril elongation. We finally explored how mutations can affect the co-assembly of hIAPP fibrils by simulating a heterogeneous system composed of wild-type and mutated hIAPP peptides. Overall, this study shows that all-atom Go-model-based simulations have the potential of discerning the effects of mutations and post-translational modifications in protein folding and association and may help in resolving the dichotomy between experimental and theoretical studies on protein folding and amyloid fibrillation.

Impact of graphene-based nanomaterials (GBNMs) on the structural and functional conformations of hepcidin peptide.

Graphene-based nanomaterials (GBNMs) are widely used in various industrial and biomedical applications. GBNMs of different compositions, size and shapes are being introduced without thorough toxicity evaluation due to the unavailability of regulatory guidelines. Computational toxicity prediction methods are used by regulatory bodies to quickly assess health hazards caused by newer materials. Due to increasing demand of GBNMs in various size and functional groups in industrial and consumer based applications, rapid and reliable computational toxicity assessment methods are urgently needed. In the present work, we investigate the impact of graphene and graphene oxide nanomaterials on the structural conformations of small hepcidin peptide and compare the materials for their structural and conformational changes. Our molecular dynamics simulation studies revealed conformational changes in hepcidin due to its interaction with GBMNs, which results in a loss of its functional properties. Our results indicate that hepcidin peptide undergo severe structural deformations when superimposed on the graphene sheet in comparison to graphene oxide sheet. These observations suggest that graphene is more toxic than a graphene oxide nanosheet of similar area. Overall, this study indicates that computational methods based on structural deformation, using molecular dynamics (MD) simulations, can be used for the early evaluation of toxicity potential of novel nanomaterials.

Effect of graphene oxide on the conformational transitions of amyloid beta peptide: A molecular dynamics simulation study.

The interactions between nanomaterials (NMs) and amyloid proteins are central to the nanotechnology-based diagnostics and therapy in neurodegenerative disorders such as Alzheimer's and Parkinson's. Graphene oxide (GO) and its derivatives have shown to modulate the aggregation pattern of disease causing amyloid beta (Aß) peptide. However, the mechanism is still not well understood. Using molecular dynamics simulations, the effect of graphene oxide (GO) and reduced graphene oxide (rGO) having carbon:oxygen ratio of 4:1 and 10:1, respectively, on the conformational transitions (alpha-helix to beta-sheet) and the dynamics of the peptide was investigated. GO and rGO decreased the beta-strand propensity of amino acid residues in Aß. The peptide displayed different modes of adsorption on GO and rGO. The adsorption on GO was dominated by electrostatic interactions, whereas on rGO, both van der Waals and electrostatic interactions contributed in the adsorption of the peptide. Our study revealed that the slight increase in the hydrophobic patches on rGO made it more effective inhibitor of conformational transitions in the peptide. Alpha helix-beta sheet transition in Aß peptide could be one of the plausible mechanism by which graphene oxide may inhibit amyloid fibrillation.

Hydration patterns of graphene-based nanomaterials (GBNMs) play a major role in the stability of a helical protein: a molecular dynamics simulation study.

Graphene-based nanomaterials (GBNMs) [graphene oxide (GO), reduced graphene oxide (rGO), and graphene] have been recognized as potential candidates for various biomedical applications ranging from biosensing platform to cellular delivery of proteins and peptides. However, GBNMs induced conformational changes in proteins are the major concerns in realizing their full potential in aforementioned applications. Despite several studies, the effect of GBNMs on the conformation of proteins is still not well understood. Therefore, an attempt was made to investigate the effect of GBNMs on the adsorption and conformation of positively charged cytoplasmic protein using molecular dynamics (MD) simulations. Our study showed that the adsorption of protein on GO was highly selective and mediated through electrostatic interactions (hydrogen bond/salt bridge interactions), whereas the van der Waals and π-π stacking interactions were the major driving forces for the adsorption of protein on rGO and graphene. The secondary structure analysis showed the conformational stability of the protein on GO may be attributed to the extensive hydration of GO surface and the absence of tyrosine residues in π-π stacking with π regions of GO. The GO surface acts as a hydrogen bond acceptor similar to the protein's natural receptor present in a physiological environment. This computational study has also explored the artificial protein receptor like potential of GO.

C60-fullerene binds with the ATP binding domain of human DNA topoiosmerase II alpha.

C60-fullerene has promising biological applications, such as drug delivery, biosensors, diagnosis and theraupetics. Despite of these applications, several in vitro studies have also reported the DNA damaging potential of this nanomaterial. Though, very little is known about the mechanism involved behind the fullerene mediated DNA damage. Our study was aimed at identifying the binding site of fullerene in the ATP binding domain of human topoisomerase II alpha, a major enzyme involved in maintaining DNA topology. In silico studies of fullerene with the enzyme demonstrated that it can interact with the active site residues of this enzyme through hydrophobic, pi-stacking and van der Waals interactions and could inhibit the activity of this enzyme.

Interaction of C60 fullerene with the proteins involved in DNA mismatch repair pathway.

Fullerenes are fascinating symmetric carbon nanostructures with a potential to bind DNA and proteins. Computational studies were performed to investigate the interaction of fullerenes with the proteins involved in DNA mismatch repair (MMR) pathway. Significant interactions of fullerene with PMS2, RFC3 and PCNA proteins were observed. These findings suggested that fullerene interferes in the human MMR process hence the subsequent disturbance may confer a large increase in spontaneous mutability and a strong predisposition to tumor development.