Search

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.

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.

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.