Measurements from hydroxyl radical footprinting (HRF) provide rich information about the solvent accessibility of amino acid side chains of a protein. Traditional HRF data analyses focus on comparing the difference in the modification/footprinting rate of a specific site to infer structural changes across two protein states, e.g., between a free and ligand-bound state. However, the rate information itself is not fully used for the purpose of comparing different protein sites within a protein on an absolute scale. To provide such a cross-site comparison, we present a new, to our knowledge, data analysis algorithm to convert the measured footprinting rate constant to a protection factor (PF) by taking into account the known intrinsic reactivity of amino acid side chain. To examine the extent to which PFs can be used for structural interpretation, this PF analysis is applied to three model systems where radiolytic footprinting data are reported in the literature. By visualizing structures colored with the PF values for individual peptides, a rational view of the structural features of various protein sites regarding their solvent accessibility is revealed, where high-PF regions are buried and low-PF regions are more exposed to the solvent. Furthermore, a detailed analysis correlating solvent accessibility and local structural contacts for gelsolin shows a statistically significant agreement between PF values and various structure measures, demonstrating that the PFs derived from this PF analysis readily explain fundamental HRF rate measurements. We also tested this PF analysis on alternative, chemical-based HRF data, showing improved correlations of structural properties of a model protein barstar compared to examining HRF rate data alone. Together, this PF analysis not only permits a novel, to our knowledge, approach of mapping protein structures by using footprinting data, but also elevates the use of HRF measurements from a qualitative, cross-state comparison to a quantitative, cross-site assessment of protein structures in the context of individual conformational states of interest.

Conformational stability of the collagen triple helix affects its turnover and determines tissue homeostasis. Although it is known that the presence of imino acids (prolines or hydroxyprolines) confer stability to the molecule, little is known regarding the stability of the imino-poor region lacking imino acids, which plays a key role in collagen cleavage. In particular, there have been continuing debates about the role of water in collagen stability. We addressed these issues using molecular dynamics simulations on 30-residue long collagen triple helices, including a structure that has a biologically relevant 9-residue imino-poor region from type III collagen (PDB ID: 1BKV). A torsional map approach was used to characterize the conformational motion of the molecule that differ between imino-rich and imino-poor regions. At temperatures 300 K and above, unwinding initiates at a common cleavage site, the glycine-isoleucine bond in the imino-poor region. This provides a linkage between previous observations that unwinding of the imino-poor region is a requirement for collagenase cleavage, and that isolated collagen molecules are unstable at body temperature. We found that unwinding of the imino-poor region is controlled by dynamic water bridges between backbone atoms with average lifetimes on the order of a few picoseconds, as the degree of unwinding strongly correlated with the loss of water bridges, and unwinding could be either prevented or enhanced, respectively by enforcing or forbidding water bridge formation. While individual water bridges were short-lived in the imino-poor region, the hydration shell surrounding the entire molecule was stable even at 330 K. The diameter of the hydrated collagen including the first hydration shell was about 14 A, in good agreement with the experimentally measured inter-collagen distances. These results elucidate the general role of water in collagen turnover: water not only affects collagen cleavage by controlling its torsional motion, but it also forms a larger-scale lubrication layer mediating collagen self-assembly.

Understanding protein-protein association is crucial in revealing the molecular basis of many biological processes. Here, we describe a theoretical simulation pipeline to study protein-protein association from an energy landscape perspective. First, a coarse-grained model is implemented and its applications are demonstrated via molecular dynamics simulations for several protein complexes. Second, an enhanced search method is used to efficiently sample a broad range of protein conformations. Third, multiple conformations are identified and clustered from simulation data and further projected on a three-dimensional globe specifying protein orientations and interacting energies. Results from several complexes indicate that the crystal-like conformation is favorable on the energy landscape even if the landscape is relatively rugged with metastable conformations. A closer examination on molecular forces shows that the formation of associated protein complexes can be primarily electrostatics-driven, hydrophobics-driven, or a combination of both in stabilizing specific binding interfaces. Taken together, these results suggest that the coarse-grained simulations and analyses provide an alternative toolset to study protein-protein association occurring in functional biomolecular complexes.

Nature provides numerous examples of self-assembly that can potentially be implemented for materials applications. Considerable attention has been given to one-dimensional cross-β or amyloid structures that can serve as templates for wire growth or strengthen materials such as glue or cement. Here, we demonstrate controlled amyloid self-assembly based on modifications of β-solenoid proteins. They occur naturally in several contexts (e.g., antifreeze proteins, drug resistance proteins) but do not aggregate in vivo due to capping structures or distortions at their ends. Removal of these capping structures and regularization of the ends of the spruce budworm and rye grass antifreeze proteins yield micron length amyloid fibrils with predictable heights, which can be a platform for biomaterial-based self-assembly. The design process, including all-atom molecular dynamics simulations, purification, and self-assembly procedures are described. Fibril formation with the predicted characteristics is supported by evidence from thioflavin-T fluorescence, circular dichroism, dynamic light scattering, and atomic force microscopy. Additionally, we find evidence for lateral assembly of the modified spruce budworm antifreeze fibrils with sufficient incubation time. The kinetics of polymerization are consistent with those for other amyloid formation reactions and are relatively fast due to the preformed nature of the polymerization nucleus.

Hydrotropes are a class of compounds that are widely used in chemical, cosmetic, and pharmaceutical industries. Though amphiphilic in character, they have short hydrophobic regions and thus differ from classical surfactants. Yet they display substantial (and often selective) ability to solubilize lipophilic compounds in water. Beyond the minimal hydrotropic concentration, the molecules self-aggregate to produce the operating assembly. We show in this paper, based on crystal structure analysis of several hydrotropes, that these compounds form open-layer assemblies, reminiscent of lamellar liquid crystals consisting of alternating hydrophobic clustering of the nonpolar regions adjacent to ionic or polar regions that are knitted together in a two-dimensional network. Stacking of aromatic rings is not seen. Two types of assemblies are seen, one with a more open and extended hydrophobic layer than the other. We suggest that the solubilizates enter the hydrophobic layers of microunits producing a cooperative and mutual stabilizing effect. The observed open layer structure of hydrotropes might also account for the occasional ability of these compounds to solubilize even better than micelles.