Traditional molecular modeling is performed at atomic resolution, which relies on X-ray and NMR experiments to provide structural information. When deal 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. Particularly, we implemented the core-weighted correlation functions to effectively recognize correct fit of component maps in complex maps, and the grid-threading Monte Carlo search algorithm to efficiently construct complex structures from electron density maps. Using EMAP, we are conducting a series collaboration studies. Below is a list of the projects during the past year. Tubline formation Collabratoed with Prof. Ruxandra Dima at University of Cincinnati, the formation of a nanotube by tubline dimer was studied. Contacts between intra dimer and inter dimmers are identified, which are important for the geometry of the nanotube. Mutation effect on nanotube properties are studied. Conformational study of Thermosom Collaborated with Prof. George Stan at University of Cincinnati, the conformational states of Thermosome is studied with the map constrained simulation method. The open state conformation was obtained through self-guided molecular dynamics simulation combined with the map constrained simulation method. The simulation results provide insight to the functional pathway of theromsome. It also demonstrate the powerful capability of the map constrained simulation method in bridging the experimental map information to structural and dynamic studies. Conformational states of KIT extracellular domain in complex with stem cell factor Collaborated with Prof. Savvas Savvides and Dr. Jonathan Elegheert at Ghent University, we studied Conformational states of KIT extracellular domain in complex with stem cell factor based on cryo-EM and SAX images. Constrained molecular dynamic simulations are applied to identify conformation states and ensembles. Molecular modeling and simulation of the gp140/sugar system Colaborated with Dr. Sriram Subramaniam at NCI we performed molecular modeling and simulation study of gp140/sugar system. GP140 is homology modeled mainly based PDB structure, 3jwd. The v1v2 loop region was modeled based on a remote homologeous PDB structure 1ciy. V3 loop was modeled based on PDB structure, 2b4c. Glycan molecules are docked on the gp140 surface using the EMAP module1 of CHARMM2. The gp140-sugar system was fit into the EM map determined from their lab with the EMAP module1 of CHARMM2 to produce the trimer system. The N- and C- terminal motifs of the trimer are fixed by assuming they binding to gp41. The rest part of the trimer is simulated using the self-guided langevin dynamics (SGLD) simulation method to promote conformational changes. In a 1000 ps SGLD simulation, we observed the conformation changed from initial closed state to a open state that is similar to the structure, 3DNO.
Showing the most recent 10 out of 21 publications
