Computational Studies of Glycan Conformations in Glycoproteins

N-glycans refer to oligosaccharide chains covalently attached to the side chain of asparagine (Asn) residues, and the majority of proteins synthesized in the endoplasmic reticulum (ER) are N-glycosylated. N-glycans can modulate the structural properties of proteins due to their close proximity to their parent proteins and their interactions between the glycan and the protein surface residues. In addition, N-glycans provide specific regions of recognition for cellular and molecular recognition. Despite their biological importance, the structural understanding of glycans and the impact of glycosylation to glycan or protein structure are lacking. I have explored the conformational freedom of glycans and their conformational preferences in different environments using structural databases and computer simulations. First, I have developed an algorithm to reliably annotate a given atomic structure of glycans. This algorithm is important because many glycan molecules in the crystal structure database are misannotated or contain errors. Using the algorithm, a database of glycans found in the PDB is constructed and available to the public. Second, the impact of glycosylation on the glycan conformation has been examined. Contrary to the common belief that the glycan conformations are independent to the protein structure, it appears that the protein structure can significantly affect the glycan structure upon glycosylation. This observation is significant because it may provide insight into protein-glycan interaction and opens up the possibility of a template-based glycan modeling approach. Third, the differences in conformational preference between glycans in solution and in glycoproteins has been examined. Using molecular dynamics (MD) simulations, the conformational preference of N-glycan pentassacharide in solution is exhaustively studied. Surprisingly, the conformational distribution is dominated by a single major conformational state and several minor conformational states. The dominant conformational state adopts a more extended conformation, thus it appears that entropy plays an important role in determining the conformational state. On the other hand, in glycoproteins, glycans can interact with surrounding protein side chains and, as a result, several conformational states are more equally populated. Based on these observations, a protocol is proposed for modeling the glycan portion of a known protein structure. It is typically more managable to acquire an atomic resolution structure or aglycoprotein (glycoprotein without glycan). In addition, the glycoform and the glycosylation site can be identified independently by mass spectrometry or NMR. The proposed modeling protocol assumes the glycosylation site, glycoform, and aglycoprotein structure are already known, and builds glycan structure models on top of the known aglycoprotein structure. The performance of the modeling protocol is greatly improved by using appropriate template structures. This protocol can be used to generate the initial model for MD simulations or refinement of low resolution models from experiments (small angle X-ray scattering and electron microscopy)

Glycan fragment database: a database of PDB-based glycan 3D structures

The glycan fragment database (GFDB), freely available at http://www.glycanstructure.org, is a database of the glycosidic torsion angles derived from the glycan structures in the Protein Data Bank (PDB). Analogous to protein structure, the structure of an oligosaccharide chain in a glycoprotein, referred to as a glycan, can be characterized by the torsion angles of glycosidic linkages between relatively rigid carbohydrate monomeric units. Knowledge of accessible conformations of biologically relevant glycans is essential in understanding their biological roles. The GFDB provides an intuitive glycan sequence search tool that allows the user to search complex glycan structures. After a glycan search is complete, each glycosidic torsion angle distribution is displayed in terms of the exact match and the fragment match. The exact match results are from the PDB entries that contain the glycan sequence identical to the query sequence. The fragment match results are from the entries with the glycan sequence whose substructure (fragment) or entire sequence is matched to the query sequence, such that the fragment results implicitly include the influences from the nearby carbohydrate residues. In addition, clustering analysis based on the torsion angle distribution can be performed to obtain the representative structures among the searched glycan structures.University of Kansas General Research Fund [2301388-003]; Kanas-COBRE NIH P20 [GM103420]; NSF [MCB-0918374]; TeraGrid/XSEDE resources [TG-MCB070009]. Funding for open access charge: NSF [MCB-0918374]

Transmembrane Helix Orientation and Dynamics: Insights from Ensemble Dynamics with Solid-State NMR Observables

This is the publisher's version. Copyright 2011 by Elsevier.As the major component of membrane proteins, transmembrane helices embedded in anisotropic bilayer environments adopt preferential orientations that are characteristic or related to their functional states. Recent developments in solid-state nuclear magnetic resonance (SSNMR) spectroscopy have made it possible to measure NMR observables that can be used to determine such orientations in a native bilayer environment. A quasistatic single conformer model is frequently used to interpret the SSNMR observables, but important motional information can be missing or misinterpreted in the model. In this work, we have investigated the orientation of the single-pass transmembrane domain of viral protein ”u“ (VpuTM) from HIV-1 by determining an ensemble of structures using multiple conformer models based on the SSNMR ensemble dynamics technique. The resulting structure ensemble shows significantly larger orientational fluctuations while the ensemble-averaged orientation is compatible with the orientation based on the quasistatic model. This observation is further corroborated by comparison with the VpuTM orientation from comparative molecular dynamics simulations in explicit bilayer membranes. SSNMR ensemble dynamics not only reveals the importance of transmembrane helix dynamics in interpretation of SSNMR observables, but also provides a means to simultaneously extract both transmembrane helix orientation and dynamics information from the SSNMR measurements

Solid-State NMR Ensemble Dynamics as a Mediator between Experiment and Simulation

This is the publisher version. Copyright 2011 by Elsevier.Solid-state NMR (SSNMR) is a powerful technique to describe the orientations of membrane proteins and peptides in their native membrane bilayer environments. The deuterium (2H) quadrupolar splitting (DQS), one of the SSNMR observables, has been used to characterize the orientations of various single-pass transmembrane (TM) helices using a semistatic rigid-body model such as the geometric analysis of labeled alanine (GALA) method. However, dynamic information of these TM helices, which could be related to important biological function, can be missing or misinterpreted with the semistatic model. We have investigated the orientation of WALP23 in an implicit membrane of dimyristoylglycerophosphocholine by determining an ensemble of structures using multiple conformer models with a DQS restraint potential. When a single conformer is used, the resulting helix orientation (tilt angle (τ) of 5.6 ± 3.2° and rotation angle (ρ) of 141.8 ± 40.6°) is similar to that determined by the GALA method. However, as the number of conformers is increased, the tilt angles of WALP23 ensemble structures become larger (26.9 ± 6.7°), which agrees well with previous molecular dynamics simulation results. In addition, the ensemble structure distribution shows excellent agreement with the two-dimensional free energy surface as a function of WALP23's τ and ρ. These results demonstrate that SSNMR ensemble dynamics provides a means to extract orientational and dynamic information of TM helices from their SSNMR observables and to explain the discrepancy between molecular dynamics simulation and GALA-based interpretation of DQS data

CHARMM-GUI Membrane Builder for Mixed Bilayers and Its Application to Yeast Membranes

This is the publisher's version. Copyright 2009 by Elsevier.Most biological membranes are composed of many different kinds of lipid, and can be characterized by the composition of the lipids. Although more and more researchers have shown their interests in molecular dynamics simulation of lipid bilayer or protein-membrane complex system, the setup of such system remains quite challenging for even relatively experienced researchers. In the previous work [1], we have shown that the setup of molecular dynamics simulation for protein-membrane complex can be dramatically simplified by automating the process and providing intuitive and straightforward user interface. In this work, we have further elaborated the process to include 25 different kinds of lipid, which makes it possible to build more biologically relevant lipid bilayers, and we also added the facility to make a lipid bilayer system alone. The efficacy of the web interface at the CHARMM-GUI website [2] has been tested by building and simulating lipid bilayer systems that resemble yeast membrane, which is composed of cholesterol, DPPC, DOPC, POPE, POPA, and POPS. In this work, we will present the usages of the mixed bilayer generation in Membrane Builder and the simulation results of the yeast membrane systems

Restricted N-glycan Conformational Space in the PDB and Its Implication in Glycan Structure Modeling

A grant from the One-University Open Access Fund at the University of Kansas was used to defray the author’s publication fees in this Open Access journal. The Open Access Fund, administered by librarians from the KU, KU Law, and KUMC libraries, is made possible by contributions from the offices of KU Provost, KU Vice Chancellor for Research & Graduate Studies, and KUMC Vice Chancellor for Research. For more information about the Open Access Fund, please see http://library.kumc.edu/authors-fund.xml.Understanding glycan structure and dynamics is central to understanding protein-carbohydrate recognition and its role in protein-protein interactions. Given the difficulties in obtaining the glycan's crystal structure in glycoconjugates due to its flexibility and heterogeneity, computational modeling could play an important role in providing glycosylated protein structure models. To address if glycan structures available in the PDB can be used as templates or fragments for glycan modeling, we present a survey of the N-glycan structures of 35 different sequences in the PDB. Our statistical analysis shows that the N-glycan structures found on homologous glycoproteins are significantly conserved compared to the random background, suggesting that N-glycan chains can be confidently modeled with template glycan structures whose parent glycoproteins share sequence similarity. On the other hand, N-glycan structures found on non-homologous glycoproteins do not show significant global structural similarity. Nonetheless, the internal substructures of these N-glycans, particularly, the substructures that are closer to the protein, show significantly similar structures, suggesting that such substructures can be used as fragments in glycan modeling. Increased interactions with protein might be responsible for the restricted conformational space of N-glycan chains. Our results suggest that structure prediction/modeling of N-glycans of glycoconjugates using structure database could be effective and different modeling approaches would be needed depending on the availability of template structures