How nature tunes sequences of disordered protein to yield the desired coiling properties is not yet well understood. To shed light on the relationship between protein sequence and elasticity, we here investigate four different natural disordered proteins with elastomeric function, namely: FG repeats in the nucleoporins; resilin in the wing tendon of dragonfly; PPAK in the muscle protein titin; and spider silk. We obtain force-extension curves for these proteins from extensive explicit solvent molecular dynamics simulations, which we compare to purely entropic coiling by modeling the four proteins as entropic chains. Although proline and glycine content are in general indicators for the entropic elasticity as expected, divergence from simple additivity is observed. Namely, coiling propensities correlate with polyproline II content more strongly than with proline content, and given a preponderance of glycines for sufficient backbone flexibility, nonlocal interactions such as electrostatic forces can result in strongly enhanced coiling, which results for the case of resilin in a distinct hump in the force-extension curve. Our results, which are directly testable by force spectroscopy experiments, shed light on how evolution has designed unfolded elastomeric proteins for different functions.

Folding pathways of the B domain of staphylococcal protein A have been sampled with a distributed computing approach. Starting from an extended structure, the method employs an index measuring topological similarity to the native structure to selectively sample trajectory branches leading to the native fold. Unperturbed and continuous folding trajectories are drawn on a physics-based atomic potential energy surface with an implicit solvent. The sampled folding trajectories demonstrate a similar sequence of events: the earlier stage involves a partial formation of helix 2 and to a less extent of helix 1 at their N terminals, followed by the hydrophobic collapse between residues F14, I17, and L18 on helix 1 and residues R28, F31, and I32 on helix 2, which results in the rigidification of the helix turn from R28 to I32. Helix 2 is then able to extend, allowing for the formation to turn 2. The above description explains one experimental result why a G30A mutant of the protein was observed to be the fastest folder among proteins of its size. And the ensemble of structures right before the final collapse is in good agreement with the transition state ensemble mapped by another recent experiment with Fersht Phi values. We emphasize that because the approach here does not provide quantifications of the free energy landscape, our model of the transition state ensemble emerges from comparisons of simulations and previous experimental results rather from the simulation results alone. On the other hand, as our approach does not rely on a low-dimensional free energy surface, it can complement methods based on the construction of free energy surfaces.

Acetylcholinesterase, with a deep, narrow active-site gorge, attracts enormous interest due to its particularly high catalytic efficiency and its inhibitors used for treatment of Alzheimer’s disease. To facilitate the massive pass-through of the substrate and inhibitors, “breathing” motions to modulate the size of the gorge are an important prerequisite. However, the molecular mechanism that governs such motions is not well explored. Here, to systematically investigate intrinsic motions of the enzyme, we performed microsecond molecular dynamics simulations on the monomer and dimer of Torpedo californica acetylcholinesterase (TcAChE) as well as the complex of TcAChE bound with the drug E2020. It has been revealed that protein-ligand interactions and dimerization both keep the gorge in bulk, and opening events of the gorge increase dramatically compared to the monomer. Dynamics of three subdomains, S3, S4 and the Ω-loop, are tightly associated with variations of the gorge size while the dynamics can be changed by ligand binding or protein dimerization. Moreover, high correlations among these subdomains provide a basis for remote residues allosterically modulating the gorge motions. These observations are propitious to expand our understanding of protein structure and function as well as providing clues for performing structure-based drug design.

Post-translational N-glycosylation of proteins is ubiquitous in eukaryotic cells, and has been shown to influence the thermodynamics of protein collapse and folding. However, the mechanism for this influence is not well understood. All-atom molecular dynamics simulations are carried out to study the collapse of a peptide linked to a single N-glycan. The glycan is shown to perturb the local water hydrogen-bonding network, rendering it less able to solvate the peptide and thus enhancing the hydrophobic contribution to the free energy of collapse. The enhancement of the hydrophobic collapse compensates for the weakened entropic coiling due to the bulky glycan chain and leads to a stronger burial of hydrophobic surface, presumably enhancing folding. This conclusion is reinforced by comparison with coarse-grained simulations, which contain no explicit solvent and correspondingly exhibit no significant thermodynamic changes on glycosylation.