A terminal α3-galactose modification regulates an E3 ubiquitin ligase subunit in Toxoplasma gondii.

Skp1, a subunit of E3 Skp1/Cullin-1/F-box protein ubiquitin ligases, is modified by a prolyl hydroxylase that mediates O2 regulation of the social amoeba Dictyostelium and the parasite Toxoplasma gondii The full effect of hydroxylation requires modification of the hydroxyproline by a pentasaccharide that, in Dictyostelium, influences Skp1 structure to favor assembly of Skp1/F-box protein subcomplexes. In Toxoplasma, the presence of a contrasting penultimate sugar assembled by a different glycosyltransferase enables testing of the conformational control model. To define the final sugar and its linkage, here we identified the glycosyltransferase that completes the glycan and found that it is closely related to glycogenin, an enzyme that may prime glycogen synthesis in yeast and animals. However, the Toxoplasma enzyme catalyzes formation of a Galα1,3Glcα linkage rather than the Glcα1,4Glcα linkage formed by glycogenin. Kinetic and crystallographic experiments showed that the glycosyltransferase Gat1 is specific for Skp1 in Toxoplasma and also in another protist, the crop pathogen Pythium ultimum The fifth sugar is important for glycan function as indicated by the slow-growth phenotype of gat1Δ parasites. Computational analyses indicated that, despite the sequence difference, the Toxoplasma glycan still assumes an ordered conformation that controls Skp1 structure and revealed the importance of nonpolar packing interactions of the fifth sugar. The substitution of glycosyltransferases in Toxoplasma and Pythium by an unrelated bifunctional enzyme that assembles a distinct but structurally compatible glycan in Dictyostelium is a remarkable case of convergent evolution, which emphasizes the importance of the terminal α-galactose and establishes the phylogenetic breadth of Skp1 glycoregulation.

Perturbing the energy landscape for improved packing during computational protein design.

The FastDesign protocol in the molecular modeling program Rosetta iterates between sequence optimization and structure refinement to stabilize de novo designed protein structures and complexes. FastDesign has been used previously to design novel protein folds and assemblies with important applications in research and medicine. To promote sampling of alternative conformations and sequences, FastDesign includes stages where the energy landscape is smoothened by reducing repulsive forces. Here, we discover that this process disfavors larger amino acids in the protein core because the protein compresses in the early stages of refinement. By testing alternative ramping strategies for the repulsive weight, we arrive at a scheme that produces lower energy designs with more native-like sequence composition in the protein core. We further validate the protocol by designing and experimentally characterizing over 4000 proteins and show that the new protocol produces higher stability proteins.

O2 sensing-associated glycosylation exposes the F-box-combining site of the Dictyostelium Skp1 subunit in E3 ubiquitin ligases.

Skp1 is a conserved protein linking cullin-1 to F-box proteins in SCF (Skp1/Cullin-1/F-box protein) E3 ubiquitin ligases, which modify protein substrates with polyubiquitin chains that typically target them for 26S proteasome-mediated degradation. In Dictyostelium (a social amoeba), Toxoplasma gondii (the agent for human toxoplasmosis), and other protists, Skp1 is regulated by a unique pentasaccharide attached to hydroxylated Pro-143 within its C-terminal F-box-binding domain. Prolyl hydroxylation of Skp1 contributes to O2-dependent Dictyostelium development, but full glycosylation at that position is required for optimal O2 sensing. Previous studies have shown that the glycan promotes organization of the F-box-binding region in Skp1 and aids in Skp1's association with F-box proteins. Here, NMR and MS approaches were used to determine the glycan structure, and then a combination of NMR and molecular dynamics simulations were employed to characterize the impact of the glycan on the conformation and motions of the intrinsically flexible F-box-binding domain of Skp1. Molecular dynamics trajectories of glycosylated Skp1 whose calculated monosaccharide relaxation kinetics and rotational correlation times agreed with the NMR data indicated that the glycan interacts with the loop connecting two α-helices of the F-box-combining site. In these trajectories, the helices separated from one another to create a more accessible and dynamic F-box interface. These results offer an unprecedented view of how a glycan modification influences a disordered region of a full-length protein. The increased sampling of an open Skp1 conformation can explain how glycosylation enhances interactions with F-box proteins in cells.

Applying Pose Clustering and MD Simulations To Eliminate False Positives in Molecular Docking.

In this work, we developed a computational protocol that employs multiple molecular docking experiments, followed by pose clustering, molecular dynamic simulations (10 ns), and energy rescoring to produce reliable 3D models of antibody-carbohydrate complexes. The protocol was applied to 10 antibody-carbohydrate co-complexes and three unliganded (apo) antibodies. Pose clustering significantly reduced the number of potential poses. For each system, 15 or fewer clusters out of 100 initial poses were generated and chosen for further analysis. Molecular dynamics (MD) simulations allowed the docked poses to either converge or disperse, and rescoring increased the likelihood that the best-ranked pose was an acceptable pose. This approach is amenable to automation and can be a valuable aid in determining the structure of antibody-carbohydrate complexes provided there is no major side chain rearrangement or backbone conformational change in the H3 loop of the CDR regions. Further, the basic protocol of docking a small ligand to a known binding site, clustering the results, and performing MD with a suitable force field is applicable to any protein ligand system.

Synthesis of 3-O-Sulfated Oligosaccharides to Understand the Relationship between Structures and Functions of Heparan Sulfate.

The sulfation at the 3-OH position of glucosamine is an important modification in forming structural domains for heparan sulfate to enable its biological functions. Seven 3-O-sulfotransferase isoforms in the human genome are involved in the biosynthesis of 3-O-sulfated heparan sulfate. As a rare modification present in heparan sulfate, the availability of 3-O-sulfated oligosaccharides is very limited. Here, we report the use of a chemoenzymatic synthetic approach to synthesize six 3-O-sulfated oligosaccharides, including three hexasaccharides and three octasaccharides. The synthesis was achieved by rearranging the enzymatic modification sequence to accommodate the substrate specificity of 3-O-sulfotransferase 3. We studied the impact of 3-O-sulfation on the conformation of the pyranose ring of 2-O-sulfated iduronic acid using NMR, and on the correlation between ring conformation and anticoagulant activity. We identified a novel octasaccharide that interacts with antithrombin and displays anti factor Xa activity. Interestingly, the octasaccharide displays a faster clearance rate than fondaparinux, an FDA-approved pentasaccharide drug, in a rat model, making this octasaccharide a potential short-acting anticoagulant drug candidate that could reduce bleeding risk. Having access to a set of critically important 3-O-sulfated oligosaccharides offers the potential to develop new heparan sulfate-based therapeutics.

Enzymatic basis for N-glycan sialylation: structure of rat α2,6-sialyltransferase (ST6GAL1) reveals conserved and unique features for glycan sialylation.

Glycan structures on glycoproteins and glycolipids play critical roles in biological recognition, targeting, and modulation of functions in animal systems. Many classes of glycan structures are capped with terminal sialic acid residues, which contribute to biological functions by either forming or masking glycan recognition sites on the cell surface or secreted glycoconjugates. Sialylated glycans are synthesized in mammals by a single conserved family of sialyltransferases that have diverse linkage and acceptor specificities. We examined the enzymatic basis for glycan sialylation in animal systems by determining the crystal structures of rat ST6GAL1, an enzyme that creates terminal α2,6-sialic acid linkages on complex-type N-glycans, at 2.4 Å resolution. Crystals were obtained from enzyme preparations generated in mammalian cells. The resulting structure revealed an overall protein fold broadly resembling the previously determined structure of pig ST3GAL1, including a CMP-sialic acid-binding site assembled from conserved sialylmotif sequence elements. Significant differences in structure and disulfide bonding patterns were found outside the sialylmotif sequences, including differences in residues predicted to interact with the glycan acceptor. Computational substrate docking and molecular dynamics simulations were performed to predict and evaluate the CMP-sialic acid donor and glycan acceptor interactions, and the results were compared with kinetic analysis of active site mutants. Comparisons of the structure with pig ST3GAL1 and a bacterial sialyltransferase revealed a similar positioning of donor, acceptor, and catalytic residues that provide a common structural framework for catalysis by the mammalian and bacterial sialyltransferases.

High-resolution epitope mapping of commercial antibodies to ANCA antigens by yeast surface display.

Epitope mapping provides critical insight into antibody-antigen interactions. Epitope mapping of autoantibodies from patients with autoimmune diseases can help elucidate disease immunogenesis and guide the development of antigen-specific therapies. Similarly, epitope mapping of commercial antibodies targeting known autoantigens enables the use of those antibodies to test specific hypotheses. Anti-Neutrophil Cytoplasmic Autoantibody (ANCA) vasculitis results from the formation of autoantibodies to multiple autoantigens, including myeloperoxidase (MPO), proteinase-3 (PR3), plasminogen (PLG), and peroxidasin (PXDN). To perform high-resolution epitope mapping of commercial antibodies to these autoantigens, we developed a novel yeast surface display library based on a series of >5000 overlapping peptides derived from their protein sequences. Using both FACS and magnetic bead isolation of reactive yeast, we screened 19 commercially available antibodies to the ANCA autoantigens. This approach to epitope mapping resulted in highly specific, fine epitope mapping, down to single amino acid resolution in many cases. Our study also identified cross-reactivity between some commercial antibodies to MPO and PXDN, which suggests that patients with apparent autoantibodies to both proteins may be the result of cross-reactivity. Together, our data validate yeast surface display using maximally overlapping peptides as an excellent approach to linear epitope mapping.

Modulation of the pharmacokinetics of soluble ACE2 decoy receptors through glycosylation.

The Spike of SARS-CoV-2 recognizes a transmembrane protease, angiotensin-converting enzyme 2 (ACE2), on host cells to initiate infection. Soluble derivatives of ACE2, in which Spike affinity is enhanced and the protein is fused to Fc of an immunoglobulin, are potent decoy receptors that reduce disease in animal models of COVID-19. Mutations were introduced into an ACE2 decoy receptor, including adding custom N-glycosylation sites and a cavity-filling substitution together with Fc modifications, which increased the decoy's catalytic activity and provided small to moderate enhancements of pharmacokinetics following intravenous and subcutaneous administration in humanized FcRn mice. Most prominently, sialylation of native glycans increases exposures by orders of magnitude, and the optimized decoy is therapeutically efficacious in a mouse COVID-19 model. Ultimately, an engineered and highly sialylated decoy receptor produced using methods suitable for manufacture of representative drug substance has high exposure with a 5- to 9-day half-life. Finally, peptide epitopes at mutated sites in the decoys generally have low binding to common HLA class II alleles and the predicted immunogenicity risk is low. Overall, glycosylation is a critical molecular attribute of ACE2 decoy receptors and modifications that combine tighter blocking of Spike with enhanced pharmacokinetics elevate this class of molecules as viable drug candidates.

Stabilizing proteins, simplified: A Rosetta-based webtool for predicting favorable mutations.

Many proteins have low thermodynamic stability, which can lead to low expression yields and limit functionality in research, industrial and clinical settings. This article introduces two, web-based tools that use the high-resolution structure of a protein along with the Rosetta molecular modeling program to predict stabilizing mutations. The protocols were recently applied to three genetically and structurally distinct proteins and successfully predicted mutations that improved thermal stability and/or protein yield. In all three cases, combining the stabilizing mutations raised the protein unfolding temperatures by more than 20°C. The first protocol evaluates point mutations and can generate a site saturation mutagenesis heatmap. The second identifies mutation clusters around user-defined positions. Both applications only require a protein structure and are particularly valuable when a deep multiple sequence alignment is not available. These tools were created to simplify protein engineering and enable research that would otherwise be infeasible due to poor expression and stability of the native molecule.

Biophysical and Structural Characterization of Novel RAS-Binding Domains (RBDs) of PI3Kα and PI3Kγ.

Phosphatidylinositol-3-kinases (PI3Ks) are lipid kinases that phosphorylate phosphatidylinositol 4,5-bisphosphate to generate a key lipid second messenger, phosphatidylinositol 3,4,5-bisphosphate. PI3Kα and PI3Kγ require activation by RAS proteins to stimulate signaling pathways that control cellular growth, differentiation, motility and survival. Intriguingly, RAS binding to PI3K isoforms likely differ, as RAS mutations have been identified that discriminate between PI3Kα and PI3Kγ, consistent with low sequence homology (23%) between their RAS binding domains (RBDs). As disruption of the RAS/PI3Kα interaction reduces tumor growth in mice with RAS- and epidermal growth factor receptor driven skin and lung cancers, compounds that interfere with this key interaction may prove useful as anti-cancer agents. However, a structure of PI3Kα bound to RAS is lacking, limiting drug discovery efforts. Expression of full-length PI3K isoforms in insect cells has resulted in low yield and variable activity, limiting biophysical and structural studies of RAS/PI3K interactions. This led us to generate the first RBDs from PI3Kα and PI3Kγ that can be expressed at high yield in bacteria and bind to RAS with similar affinity to full-length PI3K. We also solved a 2.31 Å X-ray crystal structure of the PI3Kα-RBD, which aligns well to full-length PI3Kα. Structural differences between the PI3Kα and PI3Kγ RBDs are consistent with differences in thermal stability and may underly differential RAS recognition and RAS-mediated PI3K activation. These high expression, functional PI3K RBDs will aid in interrogating RAS interactions and could aid in identifying inhibitors of this key interaction.

Author Correction: Downstream Products are Potent Inhibitors of the Heparan Sulfate 2-O-Sulfotransferase.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Downstream Products are Potent Inhibitors of the Heparan Sulfate 2-O-Sulfotransferase.

Heparan Sulfate (HS) is a cell signaling molecule linked to pathological processes ranging from cancer to viral entry, yet fundamental aspects of its biosynthesis remain incompletely understood. Here, the binding preferences of the uronyl 2-O-sulfotransferase (HS2ST) are examined with variably-sulfated hexasaccharides. Surprisingly, heavily sulfated oligosaccharides formed by later-acting sulfotransferases bind more tightly to HS2ST than those corresponding to its natural substrate or product. Inhibition assays also indicate that the IC50 values correlate simply with degree of oligosaccharide sulfation. Structural analysis predicts a mode of inhibition in which 6-O-sulfate groups located on glucosamine residues present in highly-sulfated oligosaccharides occupy the canonical binding site of the nucleotide cofactor. The unexpected finding that oligosaccharides associated with later stages in HS biosynthesis inhibit HS2ST indicates that the enzyme must be separated temporally and/or spatially from downstream products during biosynthesis in vivo, and highlights a challenge for the enzymatic synthesis of lengthy HS chains in vitro.

A Traveling Wave Ion Mobility Spectrometry (TWIMS) Study of the Robo1-Heparan Sulfate Interaction.

Roundabout 1 (Robo1) interacts with its receptor Slit to regulate axon guidance, axon branching, and dendritic development in the nervous system and to regulate morphogenesis and many cell functions in the nonneuronal tissues. This interaction is known to be critically regulated by heparan sulfate (HS). Previous studies suggest that HS is required to promote the binding of Robo1 to Slit to form the minimal signaling complex, but the molecular details and the structural requirements of HS for this interaction are still unclear. Here, we describe the application of traveling wave ion mobility spectrometry (TWIMS) to study the conformational details of the Robo1-HS interaction. The results suggest that Robo1 exists in two conformations that differ by their compactness and capability to interact with HS. The results also suggest that the highly flexible interdomain hinge region connecting the Ig1 and Ig2 domains of Robo1 plays an important functional role in promoting the Robo1-Slit interaction. Moreover, variations in the sulfation pattern and size of HS were found to affect its binding affinity and selectivity to interact with different conformations of Robo1. Both MS measurements and CIU experiments show that the Robo1-HS interaction requires the presence of a specific size and pattern of modification of HS. Furthermore, the effect of N-glycosylation on the conformation of Robo1 and its binding modes with HS is reported. Graphical Abstract ᅟ.

Vina-Carb: Improving Glycosidic Angles during Carbohydrate Docking.

Molecular docking programs are primarily designed to align rigid, drug-like fragments into the binding sites of macromolecules and frequently display poor performance when applied to flexible carbohydrate molecules. A critical source of flexibility within an oligosaccharide is the glycosidic linkages. Recently, Carbohydrate Intrinsic (CHI) energy functions were reported that attempt to quantify the glycosidic torsion angle preferences. In the present work, the CHI-energy functions have been incorporated into the AutoDock Vina (ADV) scoring function, subsequently termed Vina-Carb (VC). Two user-adjustable parameters have been introduced, namely, a CHI- energy weight term (chi_coeff) that affects the magnitude of the CHI-energy penalty and a CHI-cutoff term (chi_cutoff) that negates CHI-energy penalties below a specified value. A data set consisting of 101 protein-carbohydrate complexes and 29 apoprotein structures was used in the development and testing of VC, including antibodies, lectins, and carbohydrate binding modules. Accounting for the intramolecular energies of the glycosidic linkages in the oligosaccharides during docking led VC to produce acceptable structures within the top five ranked poses in 74% of the systems tested, compared to a success rate of 55% for ADV. An enzyme system was employed in order to illustrate the potential application of VC to proteins that may distort glycosidic linkages of carbohydrate ligands upon binding. VC represents a significant step toward accurately predicting the structures of protein-carbohydrate complexes. Furthermore, the described approach is conceptually applicable to any class of ligands that populate well-defined conformational states.

Uncovering the Relationship between Sulphation Patterns and Conformation of Iduronic Acid in Heparan Sulphate.

The L-iduronic acid (IdoA) residue is a critically important structural component in heparan sulphate polysaccharide for the biological functions. The pyranose ring of IdoA is present in (1)C4-chair, (2)SO-skew boat, and less frequently, in (4)C1-chair conformations. Here, we analyzed the conformation of IdoA residue in eight hexasaccharides by NMR. The data demonstrate a correlation between the conformation of IdoA and sulphations in the surrounding saccharide residues. For the 2-O-sulpho IdoA residue, a high degree of sulphation on neighboring residues drives ring dynamics towards the (2)SO-skew boat conformer. In contrast, the nonsulphated IdoA residue is pushed towards the (1)C4-chair conformer when the neighboring residues are highly sulphated. Our data suggest that the conformation of IdoA is regulated by the sulphation pattern of nearby saccharides that is genetically controlled by the heparan sulphate biosynthetic pathway.