Non-linearity of the collagen triple helix in solution and implications for collagen function

Collagen Structured Hydration

Collagen is a triple-helical protein unique to the extracellular matrix, conferring rigidity and stability to tissues such as bone and tendon. For the [(PPG)10]3 collagen-mimetic peptide at room temperature, our molecular dynamics simulations show that these properties result in a remarkably ordered first hydration layer of water molecules hydrogen bonded to the backbone carbonyl (bb-CO) oxygen atoms. This originates from the following observations. The radius of gyration attests that the PPG triplets are organized along a straight line, so that all triplets (excepting the ends) are equivalent. The solvent-accessible surface area (SASA) for the bb-CO oxygens shows a repetitive regularity for every triplet. This leads to water occupancy of the bb-CO sites following a similar regularity. In the crystal-phase X-ray data, as well as in our 100 K simulations, we observe a 0-2-1 water occupancy in the P-P-G triplet. Surprisingly, a similar (0-1.7-1) regularity is maintained in the liquid phase, in spite of the sub-nsec water exchange rates, because the bb-CO sites rarely remain vacant. The manifested ordered first-shell water molecules are expected to produce a cylindrical electrostatic potential around the peptide, to be investigated in future work.

A solution structure analysis reveals a bent collagen triple helix in the complement activation recognition molecule mannan-binding lectin

Collagen triple helices are critical in the function of mannan-binding lectin (MBL), an oligomeric recognition molecule in complement activation. The MBL collagen regions form complexes with the serine proteases MASP-1 and MASP-2 in order to activate complement, and mutations lead to common immunodeficiencies. To evaluate their structure-function properties, we studied the solution structures of four MBL-like collagen peptides. The thermal stability of the MBL collagen region was much reduced by the presence of a GQG interruption in the typical (X-Y-Gly)n repeat compared to controls. Experimental solution structural data were collected using analytical ultracentrifugation and small angle X-ray and neutron scattering. As controls, we included two standard Pro-Hyp-Gly collagen peptides (POG)10-13, as well as three more peptides with diverse (X-Y-Gly)n sequences that represented other collagen features. These data were quantitatively compared with atomistic linear collagen models derived from crystal structures and 12,000 conformations obtained from molecular dynamics simulations. All four MBL peptides were bent to varying degrees up to 85o in the best-fit molecular dynamics models. The best-fit benchmark peptides (POG)n were more linear but exhibited a degree of conformational flexibility. The remaining three peptides showed mostly linear solution structures. In conclusion, the collagen helix is not strictly linear, the degree of flexibility in the triple helix depends on its sequence, and the triple helix with the GQG interruption showed a pronounced bend. The bend in MBL GQG peptides resembles the bend in the collagen of complement C1q and may be key for lectin pathway activation.

Discovering design principles of collagen molecular stability using a genetic algorithm, deep learning, and experimental validation

Collagen is the most abundant structural protein in humans, providing crucial mechanical properties, including high strength and toughness, in tissues. Collagen-based biomaterials are, therefore, used for tissue repair and regeneration. Utilizing collagen effectively during materials processing ex vivo and subsequent function in vivo requires stability over wide temperature ranges to avoid denaturation and loss of structure, measured as melting temperature (T_m). Although significant research has been conducted on understanding how collagen primary amino acid sequences correspond to T_m values, a robust framework to facilitate the design of collagen sequences with specific T_m remains a challenge. Here, we develop a general model using a genetic algorithm within a deep learning framework to design collagen sequences with specific T_m values. We report 1,000 de novo collagen sequences, and we show that we can efficiently use this model to generate collagen sequences and verify their T_m values using both experimental and computational methods. We find that the model accurately predicts T_m values within a few degrees centigrade. Further, using this model, we conduct a high-throughput study to identify the most frequently occurring collagen triplets that can be directly incorporated into collagen. We further discovered that the number of hydrogen bonds within collagen calculated with molecular dynamics (MD) is directly correlated to the experimental measurement of triple-helical quality. Ultimately, we see this work as a critical step to helping researchers develop collagen sequences with specific T_m values for intended materials manufacturing methods and biomedical applications, realizing a mechanistic materials by design paradigm.

Recent Advances in Collagen Mimetic Peptide Structure and Design

Collagen mimetic peptides (CMPs) fold into a polyproline type II triple helix, allowing the study of the structure and function (or misfunction) of the collagen family of proteins. This Perspective will focus on recent developments in the use of CMPs toward understanding the structure and controlling the stability of the triple helix. Triple helix assembly is influenced by various factors, including the single amino acid propensity for the triple helix fold, pairwise interactions between these amino acids, and long-range effects observed across the helix, such as bend, twist, and fraying. Important progress in creating a comprehensive and predictive understanding of these factors for peptides with exclusively natural amino acids has been made. In contrast, several groups have successfully developed unnatural amino acids that are engineered to stabilize the triple helical structure. A third approach to controlling the triple helical structure includes covalent cross-linking of the triple helix to stabilize the assembly, which eliminates the problematic equilibrium of unfolding into monomers and enforces compositional control. Advances in all these areas have resulted in significant improvements to our understanding and control of this important class of protein, allowing for the design and application of more chemically complex and well-controlled collagen mimetic biomaterials.

ColGen: An end-to-end deep learning model to predict thermal stability of de novo collagen sequences

Collagen is the most abundant structural protein in humans, with dozens of sequence variants accounting for over 30% of the protein in an animal body. The fibrillar and hierarchical arrangements of collagen are critical in providing mechanical properties with high strength and toughness. Due to this ubiquitous role in human tissues, collagen-based biomaterials are commonly used for tissue repairs and regeneration, requiring chemical and thermal stability over a range of temperatures during materials preparation ex vivo and subsequent utility in vivo. Collagen unfolds from a triple helix to a random coil structure during a temperature interval in which the midpoint or T_m is used as a measure to evaluate the thermal stability of the molecules. However, finding a robust framework to facilitate the design of a specific collagen sequence to yield a specific T_m remains a challenge, including using conventional molecular dynamics modeling. Here we propose a de novo framework to provide a model that outputs the T_m values of input collagen sequences by incorporating deep learning trained on a large data set of collagen sequences and corresponding T_m values. By using this framework, we are able to quickly evaluate how mutations and order in the primary sequence affect the stability of collagen triple helices. Specifically, we confirm that mutations to glycines, mutations in the middle of a sequence, and short sequence lengths cause the greatest drop in T_m values.

C1q binding to surface-bound IgG is stabilized by C1r2s2 proteases

Antibody-dependent complement activation plays a major role in various pathophysiological processes in our body, including infection, inflammation, autoimmunity, and transplant rejection. In order to activate complement, antibodies should bind to target cells and recruit complement component C1. C1 is a large, multimolecular complex that consists of the antibody recognition protein C1q and a heterotetramer of proteases (C1r₂s₂). Although it is believed that interactions between C1 and IgGs are solely mediated by C1q, we here show that C1r₂s₂ proteases affect the capacity of C1q to form an avid complex with surface-bound IgG molecules. Furthermore, we demonstrate that C1q–IgG stability is influenced by IgG oligomerization and that promoting IgG oligomerization improves phagocytosis of the pathogenic bacterium Staphylococcus aureus.

Complement is an important effector mechanism for antibody-mediated clearance of infections and tumor cells. Upon binding to target cells, the antibody’s constant (Fc) domain recruits complement component C1 to initiate a proteolytic cascade that generates lytic pores and stimulates phagocytosis. The C1 complex (C1qr₂s₂) consists of the large recognition protein C1q and a heterotetramer of proteases C1r and C1s (C1r₂s₂). While interactions between C1 and IgG-Fc are believed to be mediated by the globular heads of C1q, we here find that C1r₂s₂ proteases affect the capacity of C1q to form an avid complex with surface-bound IgG molecules (on various 2,4-dinitrophenol [DNP]-coated surfaces and pathogenic Staphylococcus aureus). The extent to which C1r₂s₂ contributes to C1q–IgG stability strongly differs between human IgG subclasses. Using antibody engineering of monoclonal IgG, we reveal that hexamer-enhancing mutations improve C1q–IgG stability, both in the absence and presence of C1r₂s₂. In addition, hexamer-enhanced IgGs targeting S. aureus mediate improved complement-dependent phagocytosis by human neutrophils. Altogether, these molecular insights into complement binding to surface-bound IgGs could be important for optimal design of antibody therapies.

An Expanded Conformation of an Antibody Fab Region by X-Ray Scattering, Molecular Dynamics, and smFRET Identifies an Aggregation Mechanism

Protein aggregation is the underlying cause of many diseases, and also limits the usefulness of many natural and engineered proteins in biotechnology. Better mechanistic understanding and characterization of aggregation-prone states is needed to guide protein engineering, formulation, and drug-targeting strategies that prevent aggregation. While several final aggregated states-notably amyloids-have been characterized structurally, very little is known about the native structural conformers that initiate aggregation. We used a novel combination of small-angle x-ray scattering (SAXS), atomistic molecular dynamic simulations, single-molecule Förster resonance energy transfer, and aggregation-prone region predictions, to characterize structural changes in a native humanized Fab A33 antibody fragment, that correlated with the experimental aggregation kinetics. SAXS revealed increases in the native state radius of gyration, R_g, of 2.2% to 4.1%, at pH 5.5 and below, concomitant with accelerated aggregation. In a cutting-edge approach, we fitted the SAXS data to full MD simulations from the same conditions and located the conformational changes in the native state to the constant domain of the light chain (C_L). This C_L displacement was independently confirmed using single-molecule Förster resonance energy transfer measurements with two dual-labeled Fabs. These conformational changes were also found to increase the solvent exposure of a predicted APR, suggesting a likely mechanism through which they promote aggregation. Our findings provide a means by which aggregation-prone conformational states can be readily determined experimentally, and thus potentially used to guide protein engineering, or ligand binding strategies, with the aim of stabilizing the protein against aggregation.

Characterization by high-resolution crystal structure analysis of a triple-helix region of human collagen type III with potent cell adhesion activity

Collagen is one of the most abundant and important proteins in the human body. Human collagen type III (hCOL3A1) belongs to the fibril-forming collagens and is widely distributed in extensible connective tissue like skin, internal organs, or the vascular system. It plays key roles in wound healing, collagen fibrillogenesis, and normal cardiovascular development in human. The charged residues are considered to be an important characteristic of hCOL3A1, especially for collagen binding and recognition. Here we found that a triple helix fragment of hCOL3A1, Gly489-Gly510, contained multiple charged residues, as well as representative Glu-Lys-Gly and Glu-Arg-Gly charged triplets. We solved the crystal structure of this new fragment to a high-resolution of 1.50 Å and identified some important conformations of this new triple-helix region, including strong hydrogen bonds in interchain and interhelical interactions in addition to obvious flexible bending for the triple helix. We also found that the synthetic collagen peptides around this region exhibited potent activities through integrin-mediated peptide-membrane interaction. We then developed a method to produce a recombinant protein consisting of 16 tandem repeats of the triple-helix fragment of hCOL3A1 with strong activity without cytotoxicity. These results provide a strong base for further functional studies of human collagen type III and the method developed in this study can be applied to produce hCOL3A1-derived proteins or other tandem-repeat proteins with membrane adhesion activity.

Environmentally Controlled Curvature of Single Collagen Proteins

The predominant structural protein in vertebrates is collagen, which plays a key role in extracellular matrix and connective tissue mechanics. Despite its prevalence and physical importance in biology, the mechanical properties of molecular collagen are far from established. The flexibility of its triple helix is unresolved, with descriptions from different experimental techniques ranging from flexible to semirigid. Furthermore, it is unknown how collagen type (homo- versus heterotrimeric) and source (tissue derived versus recombinant) influence flexibility. Using SmarTrace, a chain-tracing algorithm we devised, we performed statistical analysis of collagen conformations collected with atomic force microscopy to determine the protein's mechanical properties. Our results show that types I, II, and III collagens-the key fibrillar varieties-exhibit similar molecular flexibilities. However, collagen conformations are strongly modulated by salt, transitioning from compact to extended as KCl concentration increases in both neutral and acidic pH. Although analysis with a standard worm-like chain model suggests that the persistence length of collagen can attain a wide range of values within the literature range, closer inspection reveals that this modulation of collagen's conformational behavior is not due to changes in flexibility but rather arises from the induction of curvature (either intrinsic or induced by interactions with the mica surface). By modifying standard polymer theory to include innate curvature, we show that collagen behaves as an equilibrated curved worm-like chain in two dimensions. Analysis within the curved worm-like chain model shows that collagen's curvature depends strongly on pH and salt, whereas its persistence length does not. Thus, we find that triple-helical collagen is well described as semiflexible irrespective of source, type, pH, and salt environment. These results demonstrate that collagen is more flexible than its conventional description as a rigid rod, which may have implications for its cellular processing and secretion.

Unique Conformation in a Natural Interruption Sequence of Type XIX Collagen Revealed by Its High-Resolution Crystal Structure

Naturally occurring interruptions in nonfibrillar collagen play key roles in molecular flexibility, collagen degradation, and ligand binding. The structural feature of the interruption sequences and the molecular basis for their functions have not been well studied. Here, we focused on a G5G type natural interruption sequence G-POALO-G from human type XIX collagen, a homotrimer collagen, as this sequence possesses distinct properties compared with those of a pathological similar Gly mutation sequence in collagen mimic peptides. We determined the crystal structures of the host-guest peptide (GPO)₃-GPOALO-(GPO)₄ to 1.03 Å resolution in two crystal forms. In these structures, the interruption zone brings localized disruptions to the triple helix and introduces a light 6-8° bend with the same directional preference to the whole molecule, which may correspond structurally to the first physiological kink site in type XIX collagen. Furthermore, at the G5G interruption site, the presence of Ala and Leu residues, both with free N-H groups, allows the formation of more direct and water-mediated interchain hydrogen bonds than in the related Gly → Ala structure. These could partly explain the difference in thermal stability between the different interruptions. In addition, our structures provide a detailed view of the dynamic property of such an interrupted zone with respect to hydrogen bonding topology, torsion angles, and helical parameters. Our results, for the first time, also identified the binding of zinc to the end of the triple helix. These findings will shed light on how the interruption sequence influences the conformation of the collagen molecule and provide a structural basis for further functional studies.

Structure of the C1r-C1s interaction of the C1 complex of complement activation

C1 is a large complex that triggers the destruction of invading pathogens via lysis or by stimulation of innate and adaptive immune processes. It is composed of C1q, a protein with a bouquet-like architecture, together with a tetramer assembled from two copies each of the serine proteases C1r and C1s, which activate when C1q binds to a pathogen surface. Here we describe detailed structures that show how C1r and C1s interact via an extensive interface encompassing the N-terminal regions of both proteases. Our findings reveal how the protease tetramer is organized and suggest a mechanism for the assembly and activation of C1.

The multiprotein complex C1 initiates the classical pathway of complement activation on binding to antibody–antigen complexes, pathogen surfaces, apoptotic cells, and polyanionic structures. It is formed from the recognition subcomponent C1q and a tetramer of proteases C1r₂C1s₂ as a Ca²⁺-dependent complex. Here we have determined the structure of a complex between the CUB1-EGF-CUB2 fragments of C1r and C1s to reveal the C1r–C1s interaction that forms the core of C1. Both fragments are L-shaped and interlock to form a compact antiparallel heterodimer with a Ca²⁺ from each subcomponent at the interface. Contacts, involving all three domains of each protease, are more extensive than those of C1r or C1s homodimers, explaining why heterocomplexes form preferentially. The available structural and biophysical data support a model of C1r₂C1s₂ in which two C1r-C1s dimers are linked via the catalytic domains of C1r. They are incompatible with a recent model in which the N-terminal domains of C1r and C1s form a fixed tetramer. On binding to C1q, the proteases become more compact, with the C1r-C1s dimers at the center and the six collagenous stems of C1q arranged around the perimeter. Activation is likely driven by separation of the C1r-C1s dimer pairs when C1q binds to a surface. Considerable flexibility in C1s likely facilitates C1 complex formation, activation of C1s by C1r, and binding and activation of downstream substrates C4 and C4b-bound C2 to initiate the reaction cascade.

Proline hydroxylation at different sites in hypoxia-inducible factor 1α modulates its interactions with the von Hippel–Lindau tumor suppressor protein

Proline hydroxylation of HIF-1α affects the interaction affinity between pVHL and HIF-1α and allosterically induces the conformational change of pVHL.