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.

Abstract. Molecular interactions of optical clearing agents were investigated using a combination of molecular dynamics ͑MD͒ simulations and optical spectroscopy. For a series of sugar alcohols with low to high optical clearing potential, Raman spectroscopy and integrating sphere measurements were used to quantitatively characterize tissue water loss and reduction in light scattering following agent exposures. The rate of tissue water loss was found to correlate with agent optical clearing potential, but equivalent tissue optical clearing was measured in native and fixed tissue in vitro, given long-enough exposure times to the polyol series. MD simulations showed that the rate of tissue optical clearing correlated with the preferential formation of hydrogen bond bridges between agent and collagen. Hydrogen bond bridge formation disrupts the collagen hydration layer and facilitates replacement by a chemical agent to homogenize tissue refractive index. However, the reduction in tissue light scattering did not correlate with the agent index of refraction. Our results suggest that a necessary property of optical clearing agents is hyperosmolarity to tissue, but that the most effective agents with the highest rates of optical clearing are a subset with the highest collagen solubilities.

A generalized method, termed Fast-SAXS-pro, for computing small angle x-ray scattering (SAXS) profiles of proteins, nucleic acids, and their complexes is presented. First, effective coarse-grained structure factors of DNA nucleotides are derived using a simplified two-particle-per-nucleotide representation. Second, SAXS data of a 18-bp double-stranded DNA are measured and used for the calibration of the scattering contribution from excess electron density in the DNA solvation layer. Additional test on a 25-bp DNA duplex validates this SAXS computational method and suggests that DNA has a different contribution from its hydration surface to the total scattering compared to RNA and protein. To account for such a difference, a sigmoidal function is implemented for the treatment of non-uniform electron density across the surface of a protein/nucleic-acid complex. This treatment allows differential scattering from the solvation layer surrounding protein/nucleic-acid complexes. Finally, the applications of this Fast-SAXS-pro method are demonstrated for protein/DNA and protein/RNA complexes.

Structural determination of protein-protein complexes such as multidomain nuclear receptors has been challenging for high-resolution structural techniques. Here, we present a combined use of multiple biophysical methods, termed iSPOT, an integration of shape information from small-angle X-ray scattering (SAXS), protection factors probed by hydroxyl radical footprinting, and a large series of computationally docked conformations from rigid-body or molecular dynamics (MD) simulations. Specifically tested on two model systems, the power of iSPOT is demonstrated to accurately predict the structures of a large protein-protein complex (TGFβ-FKBP12) and a multidomain nuclear receptor homodimer (HNF-4α), based on the structures of individual components of the complexes. Although neither SAXS nor footprinting alone can yield an unambiguous picture for each complex, the combination of both, seamlessly integrated in iSPOT, narrows down the best-fit structures that are about 3.3 Å and 4.2 Å in RMSD from their corresponding crystal structures, respectively. Furthermore, this proof-of-principle study based on the data synthetically derived from available crystal structures shows that the iSPOT—using either rigid-body or MD-based flexible docking—is capable of overcoming the shortcomings of standalone computational methods, especially for HNF-4α. By taking advantage of the integration of SAXS-based shape information and footprinting-based protection/accessibility as well as computational docking, this iSPOT platform is set to be a powerful approach towards accurate integrated modeling of many challenging multiprotein complexes.