Search

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.

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. This article is protected by copyright. All rights reserved.

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 A 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.

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

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.

Skp1, a subunit of E3 Skp1/Cullin-1/F-box protein ubiquitin ligases, is modified by a prolyl hydroxylase that mediates O(2) 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.

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.

Skp1, a subunit of E3 Skp1/Cullin-1/F-box protein ubiquitin ligases, is uniquely modified in protists by an O2-dependent prolyl hydroxylase that generates the attachment site for a defined pentasaccharide. Previous studies demonstrated the importance of the core glycan for growth of the parasite Toxoplasma gondii in fibroblasts, but the significance of the non-reducing terminal sugar was unknown. Here, we find that a homolog of glycogenin, an enzyme that can initiate and prime glycogen synthesis in yeast and animals, is required to catalyze the addition of an α-galactose in 3-linkage to the subterminal glucose to complete pentasaccharide assembly in cells. A strong selectivity of the enzyme (Gat1) for Skp1 in extracts is consistent with other evidence that Skp1 is the sole target of the glycosyltransferase pathway. gat1-disruption results in slow growth attesting to the importance of the terminal sugar. Molecular dynamics simulations provide an explanation for this finding and confirm the potential of the full glycan to control Skp1 organization as in the amoebozoan Dictyostelium despite the different terminal disaccharide assembled by different glycosyltransferases. Though Gat1 also exhibits low α-glucosyltransferase activity like glycogenin, autoglycosylation is not detected and gat1-disruption reveals no effect on starch accumulation. A crystal structure of the ortholog from the crop pathogen Pythium ultimum explains the distinct substrate preference and regiospecificity relative to glycogenin. A phylogenetic analysis suggests that Gat1 is related to the evolutionary progenitor of glycogenin, and acquired a role in glycogen formation following the ancestral disappearance of the underlying Skp1 glycosyltransferase prior to amoebozoan emergence.

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.

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.

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.

Background: Specificity and enzymology of glycan sialylation is poorly understood, despite its importance in biological recognition. Results: ST6GAL1 structure was determined, and substrate binding was modeled to probe active site specificity. Conclusion: The structure provides insights into the enzymatic basis of glycan sialylation. Significance: Knowledge of the enzyme structure can lead to broader understanding of enzymatic sialylation and selective inhibitor design. 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.

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 1C4-chair, 2SO-skew boat and less frequently, in 4C1-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 2SO-skew boat conformer. In contrast, the nonsulphated IdoA residue is pushed towards the 1C4-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.