Insights into the structure and dynamics of lysyl oxidase propeptide, a flexible protein with numerous partners

Lysyl oxidase (LOX) catalyzes the oxidative deamination of lysine and hydroxylysine residues in collagens and elastin, which is the first step of the cross-linking of these extracellular matrix proteins. It is secreted as a proenzyme activated by bone morphogenetic protein-1, which releases the LOX catalytic domain and its bioactive N-terminal propeptide. We characterized the recombinant human propeptide by circular dichroism, dynamic light scattering, and small-angle X-ray scattering (SAXS), and showed that it is elongated, monomeric, disordered and flexible (Dmax: 11.7 nm, Rg: 3.7 nm). We generated 3D models of the propeptide by coarse-grained molecular dynamics simulations restrained by SAXS data, which were used for docking experiments. Furthermore, we have identified 17 new binding partners of the propeptide by label-free assays. They include four glycosaminoglycans (hyaluronan, chondroitin, dermatan and heparan sulfate), collagen I, cross-linking and proteolytic enzymes (lysyl oxidase-like 2, transglutaminase-2, matrix metalloproteinase-2), a proteoglycan (fibromodulin), one growth factor (Epidermal Growth Factor, EGF), and one membrane protein (tumor endothelial marker-8). This suggests new roles for the propeptide in EGF signaling pathway

A Simple Method for Discovering Druggable, Specific Glycosaminoglycan-Protein Systems. Elucidation of Key Principles from Heparin/Heparan Sulfate-Binding Proteins

Glycosaminoglycans (GAGs) affect human physiology and pathology by modulating more than 500 proteins. GAG-protein interactions are generally assumed to be ionic and nonspecific, but specific interactions do exist. Here, we present a simple method to identify the GAG-binding site (GBS) on proteins that in turn helps predict high specific GAG–protein systems. Contrary to contemporary thinking, we found that the electrostatic potential at basic arginine and lysine residues neither identifies the GBS consistently, nor its specificity. GBSs are better identified by considering the potential at neutral hydrogen bond donors such as asparagine or glutamine sidechains. Our studies also reveal that an unusual constellation of ionic and non-ionic residues in the binding site leads to specificity. Nature engineers the local environment of Asn45 of antithrombin, Gln255 of 3-O-sulfotransferase 3, Gln163 and Asn167 of 3-O-sulfotransferase 1 and Asn27 of basic fibroblast growth factor in the respective GBSs to induce specificity. Such residues are distinct from other uncharged residues on the same protein structure in possessing a significantly higher electrostatic potential, resultant from the local topology. In contrast, uncharged residues on nonspecific GBSs such as thrombin and serum albumin possess a diffuse spread of electrostatic potential. Our findings also contradict the paradigm that GAG-binding sites are simply a collection of contiguous Arg/Lys residues. Our work demonstrates the basis for discovering specifically interacting and druggable GAG-protein systems based on the structure of protein alone, without requiring access to any structure-function relationship data

Flexibility and Explicit Solvent in Molecular-Dynamics-Based Docking of Protein–Glycosaminoglycan Systems

We present Dynamic Molecular Docking (DMD), a novel targeted molecular dynamics-based protocol developed to address ligand and receptor flexibility as well as the inclusion of explicit solvent in local molecular docking. A class of ligands for which docking performance especially benefits from overcoming these challenges is the glycosaminoglycans (GAGs). GAGs are periodic, highly flexible, and negatively charged polysaccharides playing an important role in the extracellular matrix via interaction with proteins such as growth factors and chemokines. The goal of our work has been to develop a proof of concept for an MD-based docking approach and to analyze its applicability for protein–GAG systems. DMD exploits the electrostatics-driven attraction of a ligand to its receptor, treats both as entirely flexible, and considers solvent explicitly. We show that DMD has high predictive significance for systems dominated by electrostatic attraction and demonstrate its capability to reliably identify the receptor residues contributing most to binding

Toward Improving Understanding of the Structure and Biophysics of Glycosaminoglycans

Glycosaminoglycans (GAGs) are the linear carbohydrate components of proteoglycans (PGs) that mediate PG bioactivities, including signal transduction, tissue morphogenesis, and matrix assembly. To understand GAG function, it is important to understand GAG structure and biophysics at atomic resolution. This is a challenge for existing experimental and computational methods because GAGs are heterogeneous, conformationally complex, and polydisperse, containing up to 200 monosaccharides. Molecular dynamics (MD) simulations come close to overcoming this challenge but are only feasible for short GAG polymers. To address this problem, we developed an algorithm that applies conformations from unbiased all-atom explicit-solvent MD simulations of short GAG polymers to rapidly construct 3-D atomic-resolution models of GAGs of arbitrary length. MD simulations of GAG 10-mers (i.e., polymers containing 10 monosaccharides) and 20-mers were run and conformations of all monosaccharide rings and glycosidic linkages were analyzed and compared to existing experimental data. These analyses demonstrated that (1) MD-generated GAG conformations are in agreement with existing experimental data; (2) MD-generated GAG 10-mer ring and linkage conformations match those in corresponding GAG 20-mers, suggesting that these conformations are representative of those in longer GAG biopolymers; and (3) rings and linkages in GAG 10- and 20-mers behave randomly and independently in MD simulation. Together, these findings indicate that MD-generated GAG 20-mer ring and linkage conformations can be used to construct thermodynamically-correct models of GAG polymers. Indeed, our findings demonstrate that our algorithm constructs GAG 10- and 20-mer conformational ensembles that accurately represent the backbone flexibility seen in MD simulations. Furthermore, within a day, our algorithm constructs conformational ensembles of GAG 200-mers that we would reasonably expect from MD simulation, demonstrating the efficiency of the algorithm and reduction in its time and computational cost compared to simulation. While there are other programs that can quickly construct atomic-resolution models of GAGs, those programs use conformations from short GAG subunits in solid state. Our findings suggest that GAG 20-mers are more flexible than short GAG subunits, meaning our program constructs ensembles that more accurately represent GAG polymer backbone flexibility and provide valuable insights toward improving the understanding of the structure and biophysics of GAGs

NMR and in silico studies of fucosylated chondroitin sulfate (fCS) and its interactions with selectins

This thesis describes structural studies on the interactions between the fucosylated chondroitin sulfate (fCS) oligosaccharides and human proteins known as selectins. fCS is a carbohydrate obtained from sea cucumbers, that can be classified as a branched glycosaminoglycan (GAG). It has attracted much attention due to its anti-coagulant, anti-inflammatory, antimetastatic and anti-HIV properties and its structure was previously determined by NMR. Selectins constitute a family of proteins involved in cell adhesion processes, such as inflammation, attachment of viral particles and migration of tumour cells. fCS oligosaccharides have been shown to bind to selectins, which is likely a reason behind their biological activity. However, the mechanism of this interaction is currently unknown. The initial part of the thesis describes the experimental work on expression and purification of the recombinant L- and P-selectin constructs in Pichia pastoris, Escherichia coli and HEK 293 cells. The aim of these experiments was to produce two constructs for each selectin, a single domain construct, consisting of the C-type lectin domain only, and a double domain construct, consisting of both the C-type lectin and the EGF-like domains. The intention was that the recombinant proteins would be labelled with 13C and 15N to allow for the in-depth structural NMR studies on the fCS-selectin interaction. Various experimental approaches have been explored, including the use of different cell lines, modifications to construct design, as well as alterations to expression and purification conditions. Although it was not possible to produce soluble selectin constructs in either bacterial or yeast cells, protein expression tests in HEK293 cells, performed in collaboration with the Oxford Protein Production facility (OPPF), led to production of a soluble L-selectin construct, consisting of the L-selectin C-type lectin domain. The produced L-selectin construct, as well as two commercially available constructs of the Land P-selectin extracellular domains, were used in the Saturation Transfer Difference (STD) NMR experiments to provide new information about the nature of the fCS-selectin binding. The STD experiments allowed to identify the regions within the fCS oligosaccharides that are in direct contact with the protein and likely play an important role in this interaction. Experiments on different protein constructs allowed the comparison of fCS binding to P-selectin and to two different recombinant constructs of L-selectin. Results of these studies suggest that the binding occurs via a similar mechanism for both L- and P-selectins and that the fCS oligosaccharides bind to one-domain L-selectin construct with similar affinity as to a larger construct, consisting of the entire extracellular region of the protein. Alongside the experimental work, theoretical in silico studies on the fCS-selectin binding were undertaken as part of this project. The existing X-ray structures of selectin complexes were subjected to Molecular Dynamics (MD) simulations, which allowed to explore the dynamic behaviour of E-selectin upon binding to sialyl Lewis x (sLex). It was found that sLex forms a more favourable interaction with the extended conformation of E-selectin and that the protein in this conformation is characterised by a high degree of interdomain flexibility, with a new type of interdomain movement observed in the MD studies on this complex. In further in silico studies, the fCS oligosaccharides were docked to the existing P-selectin structures. The docking tests were performed on the computationally produced fCS trisaccharides with fucose branches either 2,4 or 3,4-sulfated. Results were evaluated with MD simulations and analysed in the light of current knowledge of selectin-ligand binding and the STD NMR experimental results. The in silico studies allowed to identify a subset of P-selectin residues that are likely involved in the interaction with fCS oligosaccharides in vivo. The conformational behaviour of P-selectin upon binding to fCS was also explored and it was found that the interdomain hinge is flexible during this interaction and allows transition from bent to extended conformational state. Finally, a new NMR method was developed to facilitate the studies of complex carbohydrates, incorporating the concepts of G-matrix Fourier Transform (GFT) NMR into 2D HSQC and 2D HSQC-TOCSY experiments. The method allows to separate peaks in the regions of high spectral overlap, providing information that can simplify the assignment process. The new experiments facilitated the structural evaluation of a sample containing a mixture of oligosaccharides resulting from the depolymerisation of fCS polysaccharide