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. 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. EM study of maturation of the HIV-1 core The formation of the HIV-1 core is the final step in the viral maturation pathway, resulting in the formation of infectious virus. Most current models for HIV-1 core formation suggest that, upon proteolytic cleavage from the immature Gag, capsid (CA) dissociates into the viral interior before reforming into the core. Here we present evidence for an alternate view of core formation by taking advantage of our serendipitous observation of large membrane-enclosed structures in HIV-1 supernatants from infected cells. Cryo-electron tomographic studies show that these structures, which contain ordered arrays of what is likely the membrane-associated matrix protein, contain multiple cores that can be captured at different stages of maturation. Our studies suggest that HIV maturation involves a non-diffusional phase transition in which the detaching layer of the cleaved CA lattice is gradually converted into a roll that ultimately forms the surface of the mature conical core. Protein binding sites network The global organization of protein binding sites is analyzed by constructing a weighted network of 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.
Showing the most recent 10 out of 21 publications