Biological Macromolecules Traditional molecular modeling is performed at atomic resolution, which relies on X-ray and NMR experiments to provide structural information. When dealing with biomolecular assemblies of millions of atoms, atomic description of molecular objects becomes very computational inefficient. We developed a method that uses map objects for molecular modeling to efficiently derive structural information from experimental maps, as well as conveniently manipulate map objects, perform conformational search directly using map objects. This development work has been implemented into CHARMM as the EMAP module. This implementation enables CHARMM to manipulate map objects, including map input, output, comparison, docking, etc. Other experiment such as transition metal ion FRET (tmFRET) is becoming a useful way to obtain protein structure information. A new focus of our research is to combine efficient simulation technique with structural information from experiment to assist high throughput protein structure determination. Structure determination from low resolution EM maps We developed a map-restrained self-guided Langevin dynamics (MapSGLD) simulation method for efficient targeted conformational search. The targeted conformational search represents simulations under restraints defined by experimental observations and/or by user specified structural requirements. Through map-restraints, this method provides an efficient way to maintain substructures and to set structure targets during conformational searching. With an enhanced conformational searching ability of self-guided Langevin dynamics, this approach is suitable for simulating large-scale conformational changes, such as the formation of macromolecular assemblies and transitions between different conformational states. A direct application of this method is to determine macromolecular structures by flexible fitting of atomic structures into density maps derived from cryo-electron microscopy. Analysis of waters positions identified by X-ray versus cyro-EM. Correlation between water molecules identified in atomic models of -galactosidase determined by cryo-EM and X-ray crystallography. This work is a collaboration with a cryo-EM laboratory at the National Cancer institute, NIH. It is based on their recently published 2.2- resolution solution structure of -galactosidase, that contains resolved water densities. This study analyses how the water positions determined by cryo-EM compare to conserved water across all crystallography-determined structures. Interactions between water molecules and amino acids occur both at the surface of protein and within the structure, including areas such as subunit interfaces, catalytic sites, and other cavities. Identifying protein solvation profiles are critical to understand both protein structure and function, as well as for the design of high affinity lead compounds for drug discovery. The global organization of protein binding sites. A global organization of protein binding sites is obtained by constructing a weighted network of all known binding sites based on their structural similarities and detecting communities of structurally similar binding sites based on the minimum description length principle. The analysis reveals that there are two central binding site communities that play the roles of the network hubs of smaller peripheral communities. The sizes of communities follow a power-law distribution, which indicates that the binding sites included in larger communities may be older and have been evolutionary structural scaffolds of more recent ones. Structurally similar binding sites in the same community bind to diverse ligands promiscuously and they are also embedded in diverse domain structures. Understanding the general principles of binding site interplay will pave the way for improved drug design and protein design. Structure mechanism of Glutamate receptor activation Ionotropic glutamate receptors are cation channels that mediate signal transmission by depolarizing the postsynapitic membrane in response to L-glutamate release from the presynaptic neuron. Within the iGluR family of receptors are a-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptor (AMPA), kainite(KA), and N-methyl-D-aspartate (NMDA) subtypes, receptors that are all activated by glutamate and related in amino acid sequence, yet distinct in overall architecture, pharmacology, and biophysical characteristics. AMPA receptors are tetrameric complexes composed of subunits with a modular domain arrangement, beginning with the amino-terminal domain(ATD), the ligand- or agonist-binding domain (LBD), and the pore-forming transmembrane domain (TMD). Because AMPA receptors undergo rapid and nearly complete desensitization in the continued presence of agonist, it has proven difficult to elucidate high-resolution structures of agonist-bound, activated states and to define mechanism by which the chemical potential of agonist binding is transduced into the mechanical force of ion channel gating. The map-restrained self-guided Langevin dynamics (MapSGLD) simulation method we developed previously can utilize structural information embedded in a force field to flexibly fit macromolecular systems into low resolution maps to obtain energetically favored atomic structures that satisfy the maps. We perform flexible fitting with MapSGLD to obtain atomic structures of the glutamate receptor from EM maps. The open state atomic structure of the glutamate receptor shows the LBD in the clamshell closed conformation that agrees with the LBD x-ray structure. In addition to structural determination, MapSGLD provides dynamic information about the transition between different states.

Project Start
Project End
Budget Start
Budget End
Support Year
35
Fiscal Year
2017
Total Cost
Indirect Cost
Name
U.S. National Heart Lung and Blood Inst
Department
Type
DUNS #
City
State
Country
Zip Code
Bartesaghi, Alberto; Aguerrebere, Cecilia; Falconieri, Veronica et al. (2018) Atomic Resolution Cryo-EM Structure of ?-Galactosidase. Structure 26:848-856.e3
Nagy, Gergely N; Suardíaz, Reynier; Lopata, Anna et al. (2016) Structural Characterization of Arginine Fingers: Identification of an Arginine Finger for the Pyrophosphatase dUTPases. J Am Chem Soc 138:15035-15045
Rinschen, Markus M; Bharill, Puneet; Wu, Xiongwu et al. (2016) The ubiquitin ligase Ubr4 controls stability of podocin/MEC-2 supercomplexes. Hum Mol Genet 25:1328-44
Bartesaghi, Alberto; Merk, Alan; Banerjee, Soojay et al. (2015) 2.2 Å resolution cryo-EM structure of ?-galactosidase in complex with a cell-permeant inhibitor. Science 348:1147-51
Frank, Gabriel A; Narayan, Kedar; Bess Jr, Julian W et al. (2015) Maturation of the HIV-1 core by a non-diffusional phase transition. Nat Commun 6:5854
Rinschen, Markus M; Pahmeyer, Caroline; Pisitkun, Trairak et al. (2015) Comparative phosphoproteomic analysis of mammalian glomeruli reveals conserved podocin C-terminal phosphorylation as a determinant of slit diaphragm complex architecture. Proteomics 15:1326-31
Lee, Juyong; Lee, Kiho; Joung, InSuk et al. (2015) Sigma-RF: prediction of the variability of spatial restraints in template-based modeling by random forest. BMC Bioinformatics 16:94
Konc, Janez; Miller, Benjamin T; Štular, Tanja et al. (2015) ProBiS-CHARMMing: Web Interface for Prediction and Optimization of Ligands in Protein Binding Sites. J Chem Inf Model 55:2308-14
Joo, Keehyoung; Joung, InSuk; Lee, Jinhyuk et al. (2015) Protein structure determination by conformational space annealing using NMR geometric restraints. Proteins 83:2251-62
Rinschen, Markus M; Wu, Xiongwu; König, Tim et al. (2014) Phosphoproteomic analysis reveals regulatory mechanisms at the kidney filtration barrier. J Am Soc Nephrol 25:1509-22

Showing the most recent 10 out of 21 publications