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.

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.

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.

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

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.

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.

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.