Glycan Reader: Automated Sugar Identification and Simulation Preparation for Carbohydrates and Glycoproteins

This is the peer reviewed version of the following article: Jo, S., Song, K. C., Desaire, H., MacKerell, A. D., & Im, W. (2011). Glycan Reader: Automated Sugar Identification and Simulation Preparation for Carbohydrates and Glycoproteins. Journal of Computational Chemistry, 32(14), 3135–3141. http://doi.org/10.1002/jcc.21886, which has been published in final form at http://doi.org/10.1002/jcc.21886. This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Self-Archiving.Understanding how glycosylation affects protein structure, dynamics, and function is an emerging and challenging problem in biology. As a first step toward glycan modeling in the context of structural glycobiology, we have developed Glycan Reader and integrated it into the CHARMMGUI, http://www.charmm-gui.org/input/glycan. Glycan Reader greatly simplifies the reading of PDB structure files containing glycans through (i) detection of carbohydrate molecules, (ii) automatic annotation of carbohydrates based on their three-dimensional structures, (iii) recognition of glycosidic linkages between carbohydrates as well as N-/O-glycosidic linkages to proteins, and (iv) generation of inputs for the biomolecular simulation program CHARMM with the proper glycosidic linkage setup. In addition, Glycan Reader is linked to other functional modules in CHARMM-GUI, allowing users to easily generate carbohydrate or glycoprotein molecular simulation systems in solution or membrane environments and visualize the electrostatic potential on glycoprotein surfaces. These tools are useful for studying the impact of glycosylation on protein structure and dynamics

Improving the force field description of tyrosine-choline cation-π interactions : QM investigation of phenol-N(Me)₄⁺ interactions

Cation-pi interactions between tyrosine amino acids and compounds containing N,N,N-trimethylethanolammonium (N(CH3)(3)) are involved in the recognition of histone tails by chromodomains and in the recognition of phosphatidylcholine (PC) phospholipids by membrane-binding proteins. Yet, the lack of explicit polarization or charge transfer effects in molecular mechanics force fields raises questions about the reliability of the representation of these interactions in biomolecular simulations. Here, we investigate the nature of phenol tetramethylammonium (TMA) interactions using quantum mechanical (QM) calculations, which we also use to evaluate the accuracy of the additive CHARIVIM36 and Drude polarizable force fields in modeling tyrosine-choline interactions. We show that the potential energy surface (PES) obtained using SAPT2+/aug-cc-pVDZ compares well with the large basis-set CCSD(T) PES when TMA approaches the phenol ring perpendicularly. Furthermore, the SAPT energy decomposition reveals comparable contributions from electrostatics and dispersion in phenol-TMA interactions. We then compared the SAPT2+/augcc-pVDZ PES obtained along various approach directions to the corresponding PES obtained with CHARMM, and we show that the force field accurately reproduces the minimum distances while the interaction energies are underestimated. The use of the Drude polarizable force field significantly improves the interaction energies but decreases the agreement on distances at energy minima. The best agreement between force field and QM PES is obtained by modifying the Lennard-Jones terms for atom pairs involved in the phenol-TMA cation-pi interactions. This is further shown to improve the correlation between the occupancy of tyrosine-choline cation-pi interactions obtained from molecular dynamics simulations of a bilayer-bound bacterial phospholipase and experimental affinity data of the wild-type protein and selected mutants

Identification of a novel transport system in Borrelia burgdorferi that links the inner and outer membranes

Borrelia burgdorferi, the spirochete that causes Lyme disease, is a diderm organism that is similar to Gram-negative organisms in that it contains both an inner and outer membrane. Unlike typical Gram-negative organisms, however, B. burgdorferi lacks lipopolysaccharide (LPS). Using computational genome analyses and structural modeling, we identified a transport system containing six proteins in B. burgdorferi that are all orthologs to proteins found in the lipopolysaccharide transport (LPT) system that links the inner and outer membranes of Gram-negative organisms and is responsible for placing LPS on the surface of these organisms. While B. burgdorferi does not contain LPS, it does encode over 100 different surface-exposed lipoproteins and several major glycolipids, which like LPS are also highly amphiphilic molecules, though no system to transport these molecules to the borrelial surface is known. Accordingly, experiments supplemented by molecular modeling were undertaken to determine whether the orthologous LPT system identified in B. burgdorferi could transport lipoproteins and/or glycolipids to the borrelial outer membrane. Our combined observations strongly suggest that the LPT transport system does not transport lipoproteins to the surface. Molecular dynamic modeling, however, suggests that the borrelial LPT system could transport borrelial glycolipids to the outer membrane

Unusual Sequence Effects on Nucleotide Excision Repair of Arylamine Lesions: DNA Bending/Distortion as a Primary Recognition Factor

The environmental arylamine mutagens are implicated in the etiology of various sporadic human cancers. Arylamine-modified dG lesions were studied in two fully paired 11-mer duplexes with a -G*C N- sequence context, in which G* is a C8-substituted dG adduct derived from fluorinated analogs of 4-aminobiphenyl (FABP), 2-aminofluorene (FAF) or 2-acetylaminofluorene (FAAF), and N is either dA or dT. The FABP and FAF lesions exist in a simple mixture of ‘stacked’ (S) and ‘B-type’ (B) conformers, whereas the N-acetylated FAAF also samples a ‘wedge’ (W) conformer. FAAF is repaired three to four times more efficiently than FABP and FAF. A simple A- to -T polarity swap in the G*C A/G*CT transition produced a dramatic increase in syn-conformation and resulted in 2- to 3-fold lower nucleotide excision repair (NER) efficiencies in Escherichia coli. These results indicate that lesion-induced DNA bending/thermodynamic destabilization is an important DNA damage recognition factor, more so than the local S/B-conformational heterogeneity that was observed previously for FAF and FAAF in certain sequence contexts. This work represents a novel 3\u27-next flanking sequence effect as a unique NER factor for bulky arylamine lesions in E. coli

Characterizing Structural Transitions Using Localized Free Energy Landscape Analysis

Structural changes in molecules are frequently observed during biological processes like replication, transcription and translation. These structural changes can usually be traced to specific distortions in the backbones of the macromolecules involved. Quantitative energetic characterization of such distortions can greatly advance the atomic-level understanding of the dynamic character of these biological processes.Molecular dynamics simulations combined with a variation of the Weighted Histogram Analysis Method for potential of mean force determination are applied to characterize localized structural changes for the test case of cytosine (underlined) base flipping in a GTCAGCGCATGG DNA duplex. Free energy landscapes for backbone torsion and sugar pucker degrees of freedom in the DNA are used to understand their behavior in response to the base flipping perturbation. By simplifying the base flipping structural change into a two-state model, a free energy difference of upto 14 kcal/mol can be attributed to the flipped state relative to the stacked Watson-Crick base paired state. This two-state classification allows precise evaluation of the effect of base flipping on local backbone degrees of freedom.The calculated free energy landscapes of individual backbone and sugar degrees of freedom expectedly show the greatest change in the vicinity of the flipping base itself, but specific delocalized effects can be discerned upto four nucleotide positions away in both 5' and 3' directions. Free energy landscape analysis thus provides a quantitative method to pinpoint the determinants of structural change on the atomic scale and also delineate the extent of propagation of the perturbation along the molecule. In addition to nucleic acids, this methodology is anticipated to be useful for studying conformational changes in all macromolecules, including carbohydrates, lipids, and proteins

A two-domain elevator mechanism for sodium/proton antiport

Sodium/proton (Na+/H+) antiporters, located at the plasma membrane in every cell, are vital for cell homeostasis1. In humans, their dysfunction has been linked to diseases, such as hypertension, heart failure and epilepsy, and they are well-established drug targets2. The best understood model system for Na+/H+ antiport is NhaA from Escherichia coli1, 3, for which both electron microscopy and crystal structures are available4, 5, 6. NhaA is made up of two distinct domains: a core domain and a dimerization domain. In the NhaA crystal structure a cavity is located between the two domains, providing access to the ion-binding site from the inward-facing surface of the protein1, 4. Like many Na+/H+ antiporters, the activity of NhaA is regulated by pH, only becoming active above pH 6.5, at which point a conformational change is thought to occur7. The only reported NhaA crystal structure so far is of the low pH inactivated form4. Here we describe the active-state structure of a Na+/H+ antiporter, NapA from Thermus thermophilus, at 3 Å resolution, solved from crystals grown at pH 7.8. In the NapA structure, the core and dimerization domains are in different positions to those seen in NhaA, and a negatively charged cavity has now opened to the outside. The extracellular cavity allows access to a strictly conserved aspartate residue thought to coordinate ion binding1, 8, 9 directly, a role supported here by molecular dynamics simulations. To alternate access to this ion-binding site, however, requires a surprisingly large rotation of the core domain, some 20° against the dimerization interface. We conclude that despite their fast transport rates of up to 1,500 ions per second3, Na+/H+ antiporters operate by a two-domain rocking bundle model, revealing themes relevant to secondary-active transporters in general

Computational Fragment-Based Binding Site Identification by Ligand Competitive Saturation

Fragment-based drug discovery using NMR and x-ray crystallographic methods has proven utility but also non-trivial time, materials, and labor costs. Current computational fragment-based approaches circumvent these issues but suffer from limited representations of protein flexibility and solvation effects, leading to difficulties with rigorous ranking of fragment affinities. To overcome these limitations we describe an explicit solvent all-atom molecular dynamics methodology (SILCS: Site Identification by Ligand Competitive Saturation) that uses small aliphatic and aromatic molecules plus water molecules to map the affinity pattern of a protein for hydrophobic groups, aromatic groups, hydrogen bond donors, and hydrogen bond acceptors. By simultaneously incorporating ligands representative of all these functionalities, the method is an in silico free energy-based competition assay that generates three-dimensional probability maps of fragment binding (FragMaps) indicating favorable fragment∶protein interactions. Applied to the two-fold symmetric oncoprotein BCL-6, the SILCS method yields two-fold symmetric FragMaps that recapitulate the crystallographic binding modes of the SMRT and BCOR peptides. These FragMaps account both for important sequence and structure differences in the C-terminal halves of the two peptides and also the high mobility of the BCL-6 His116 sidechain in the peptide-binding groove. Such SILCS FragMaps can be used to qualitatively inform the design of small-molecule inhibitors or as scoring grids for high-throughput in silico docking that incorporate both an atomic-level description of solvation and protein flexibility