Unfolding of α-helical 20-residue poly-glutamic acid analyzed by multiple runs of canonical molecular dynamics simulations

Genome-wide survey of D/E repeats in human proteins uncovers their instability and aids in identifying their role in the chromatin regulator ATAD2

D/E repeats are stretches of aspartic and/or glutamic acid residues found in over 150 human proteins. We examined genomic stability of D/E repeats and functional characteristics of D/E repeat-containing proteins vis-à-vis the proteins with poly-Q or poly-A repeats, which are known to undergo pathologic expansions. Mining of tumor sequencing data revealed that D/E repeat-coding regions are similar to those coding poly-Qs and poly-As in increased incidence of trinucleotide insertions/deletions but differ in types and incidence of substitutions. D/E repeat-containing proteins preferentially function in chromatin metabolism and are the more likely to be nuclear and interact with core histones, the longer their repeats are. One of the longest D/E repeats of unknown function is in ATAD2, a bromodomain family ATPase frequently overexpressed in tumors. We demonstrate that D/E repeat deletion in ATAD2 suppresses its binding to nascent and mature chromatin and to the constitutive pericentromeric heterochromatin, where ATAD2 represses satellite transcription.

•

Many human proteins contain runs of aspartic/glutamic acid residues (D/E repeats)

•

D/E repeats show increased incidence of in-frame insertions/deletions in tumors

•

Nuclear and histone-interacting proteins often have long D/E repeats

•

D/E repeat of the oncogene ATAD2 controls its binding to pericentric chromatin

Length-dependent collective vibrational dynamics in alpha-helices

Functions of protein molecules in nature are closely associated with their well-defined three-dimensional structures and dynamics in body fluid. So far, many efforts have been made to reveal the relation of protein structure, dynamics, and function, but the structural origin of protein dynamics, especially for secondary structures as building blocks of conformation transition, is still ambiguous. Here we theoretically uncover the collective vibrations of elastic poly-alanine α-helices and find vibration patterns that are distinctively different over residue numbers ranging from 20 to 80. Contrary to the decreasing vibration magnitude from ends to the middle region for short helices, the vibration magnitude for long helices takes the minimum at approximately 1/5 of helix length from ends but reaches a peak at the center. Further analysis indicates that major vibrational modes of helical structures strongly depend on their residue numbers, where the twist mode dominates in the vibrations of short helices while the bend mode dominates the long ones analogous to an elastic Euler beam. The helix-coil transition pathway is also affected by the alternation of the first-order mode in helices with different lengths. The dynamic properties of the helical polypeptides are promising to be harnessed for de novo design of protein-based materials and artificial biomolecules in clinical treatments.

pH-Induced Changes in Polypeptide Conformation: Force-Field Comparison with Experimental Validation

Microsecond-long all-atom molecular dynamics (MD) simulations, circular dichroism, laser Doppler velocimetry, and dynamic light-scattering techniques have been used to investigate pH-induced changes in the secondary structure, charge, and conformation of poly l-lysine (PLL) and poly l-glutamic acid (PGA). The employed combination of the experimental methods reveals for both PLL and PGA a narrow pH range at which they are charged enough to form stable colloidal suspensions, maintaining their α-helix content above 60%; an elevated charge state of the peptides required for colloidal stability promotes the peptide solvation as a random coil. To obtain a more microscopic view on the conformations and to verify the modeling performance, peptide secondary structure and conformations rising in MD simulations are also examined using three different force fields, i.e., OPLS-AA, CHARMM27, and AMBER99SB*-ILDNP. Ramachandran plots reveal that in the examined setup the α-helix content is systematically overestimated in CHARMM27, while OPLS-AA overestimates the β-sheet fraction at lower ionization degrees. At high ionization degrees, the OPLS-AA force-field-predicted secondary structure fractions match the experimentally measured distribution most closely. However, the pH-induced changes in PLL and PGA secondary structure are reasonably captured only by the AMBER99SB*-ILDNP force field, with the exception of the fully charged PGA in which the α-helix content is overestimated. The comparison to simulations results shows that the examined force fields involve significant deviations in their predictions for charged homopolypeptides. The detailed mapping of secondary structure dependency on pH for the polypeptides, especially finding the stable colloidal α-helical regime for both examined peptides, has significant potential for practical applications of the charged homopolypeptides. The findings raise attention especially to the pH fine tuning as an underappreciated control factor in surface modification and self-assembly.

Effects of number of parallel runs and frequency of bias-strength replacement in generalized ensemble molecular dynamics simulations

Background

The generalized ensemble approach with the molecular dynamics (MD) method has been widely utilized. This approach usually has two features. (i) A bias potential, whose strength is replaced during a simulation, is applied. (ii) Sampling can be performed by many parallel runs of simulations. Although the frequency of the bias-strength replacement and the number of parallel runs can be adjusted, the effects of these settings on the resultant ensemble remain unclear.

Method

In this study, we performed multicanonical MD simulations for a foldable mini-protein (Trp-cage) and two unstructured peptides (8- and 20-residue poly-glutamic acids) with various settings.

Results

As a result, running many short simulations yielded robust results for the Trp-cage model. Regarding the frequency of the bias-potential replacement, although using a high frequency enhanced the traversals in the potential energy space, it did not promote conformational changes in all the systems.