Search

Shortly after the determination of the first protein x-ray crystal structures, researchers analyzed their cores and reported packing fractions $\phi \approx 0.75$, a value that is similar to close packing equal-sized spheres. A limitation of these analyses was the use of `extended atom' models, rather than the more physically accurate `explicit hydrogen' model. The validity of using the explicit hydrogen model is proved by its ability to predict the side chain dihedral angle distributions observed in proteins. We employ the explicit hydrogen model to calculate the packing fraction of the cores of over $200$ high resolution protein structures. We find that these protein cores have $\phi \approx 0.55$, which is comparable to random close-packing of non-spherical particles. This result provides a deeper understanding of the physical basis of protein structure that will enable predictions of the effects of amino acid mutations and design of new functional proteins., Comment: 5 pages, 4 figures

We perform extensive coarse-grained (CG) Langevin dynamics simulations of intrinsically disordered proteins (IDPs), which possess fluctuating conformational statistics between that for excluded volume random walks and collapsed globules. Our CG model includes repulsive steric, attractive hydrophobic, and electrostatic interactions between residues and is calibrated to a large collection of single-molecule fluorescence resonance energy transfer data on the inter-residue separations for 36 pairs of residues in five IDPs: $\alpha$-, $\beta$-, and $\gamma$-synuclein, the microtubule-associated protein $\tau$, and prothymosin $\alpha$. We find that our CG model is able to recapitulate the average inter-residue separations regardless of the choice of the hydrophobicity scale, which shows that our calibrated model can robustly capture the conformational dynamics of IDPs. We then employ our model to study the scaling of the radius of gyration with chemical distance in 11 known IDPs. We identify a strong correlation between the distance to the dividing line between folded proteins and IDPs in the mean charge and hydrophobicity space and the scaling exponent of the radius of gyration with chemical distance along the protein., Comment: 16 pages, 10 figures

We perform molecular dynamics simulations to compress binary hard spheres into jammed packings as a function of the compression rate $R$, size ratio $\alpha$, and number fraction $x_S$ of small particles to determine the connection between the glass-forming ability (GFA) and packing efficiency in bulk metallic glasses (BMGs). We define the GFA by measuring the critical compression rate $R_c$, below which jammed hard-sphere packings begin to form "random crystal" structures with defects. We find that for systems with $\alpha \gtrsim 0.8$ that do not de-mix, $R_c$ decreases strongly with $\Delta \phi_J$, as $R_c \sim \exp(-1/\Delta \phi_J^2)$, where $\Delta \phi_J$ is the difference between the average packing fraction of the amorphous packings and random crystal structures at $R_c$. Systems with $\alpha \lesssim 0.8$ partially de-mix, which promotes crystallization, but we still find a strong correlation between $R_c$ and $\Delta \phi_J$. We show that known metal-metal BMGs occur in the regions of the $\alpha$ and $x_S$ parameter space with the lowest values of $R_c$ for binary hard spheres. Our results emphasize that maximizing GFA in binary systems involves two competing effects: minimizing $\alpha$ to increase packing efficiency, while maximizing $\alpha$ to prevent de-mixing., Comment: 5 pages, 4 figures

Intrinsically disordered proteins (IDPs) do not possess well-defined three-dimensional structures in solution under physiological conditions. We develop all-atom, united-atom, and coarse-grained Langevin dynamics simulations for the IDP alpha-synuclein that include geometric, attractive hydrophobic, and screened electrostatic interactions and are calibrated to the inter-residue separations measured in recent smFRET experiments. We find that alpha-synuclein is disordered with conformational statistics that are intermediate between random walk and collapsed globule behavior. An advantage of calibrated molecular simulations over constraint methods is that physical forces act on all residues, not only on residue pairs that are monitored experimentally, and these simulations can be used to study oligomerization and aggregation of multiple alpha-synuclein proteins that may precede amyloid formation., Comment: 14 pages, 12 figures, 2 tables