New theoretical techniques are being developed and characterized. These efforts are usually coupled with software development, and involve the systematic testing and evaluation of new ideas. This development is driven by current needs and interests. Specific ongoing projects include: - Development of electric density map docking utility (EMAP) - The core-weighted fitting method to construct molecular assemblies from EM maps - Molecular modeling using low resolution maps - Thermodynamic properties in Self-Guided Molecular Dynamics (SGMD) simulations - Isotropic Periodic Sum method for the calculation of long range interactions - Anisotropic long range corrections for truncated van der Waal (LJ) interactions - Development and validation of methods for simulating interfacial systems in P21 boundary conditions (single surface) - Development of methods for examining reaction mechanism in complex systems - Unbiased forced sampling of complex conformational transitions and estimation of the potential of mean force along the reaction pathway - Development of the REPLICA/PATH method for determining reaction paths in complex systems using simulated annealing - Development of lineal scaling methods for ab initio combined Quantum Mechanical/Molecular Modeling - Enhancements of QM/MM potentials (using Gaussian delocalize MM charges, double link atom method) - GAMESS-UK and CHARMM integration for QM/MM applications - Density functional QM/MM using a double link atom interface - Development of accurate interaction energy calculations for macromolecules - Development of a rapid search strategy for docking two macromolecules Considerable efforts has been aimed at molecular modeling and structure determination using low resolution maps. This involved the development of electric density map docking utility (EMAP) module in CHARMM. we have developed a core-weighting approach to fit atomic structures into low resolution EM maps of biomolecular assembly with multiple components. The proposed core-weighted correlations have significantly improved sensitivity to distinguish the correct fit when compared with more traditional correlations. The construction of a molecular model for a complex macromolecular assembly is thus simplified from a many-body search problem to a series of single-body search problems, making the computational search for the correct fit much easier. Combined with the core-weighted correlation function, a grid-threading Monte Carlo (GTMC) approach is developed to search the best fit efficiently. Quantum mechanical/molecular mechanical (QM/MM) techniques are extremely useful in the theoretical examination of competing reaction pathways in enzyme mechanisms. GAMESS-UK has been tightly integrated into CHARMM to allow studies of catalytic paths in small molecules and enzyme complexes. This extends the QM/MM suite within CHARMM since GAMESS-UK provides DFT methods. Dr. Woodcock has primary been focused on developing and maintaining QM/MM interfaces as well as adding functionality to the existing QM/MM Replica/Path and Nudged Elastic Band (NEB) methods. A recent QM/MM enhancement added support for the ab initio software package Q-Chem to CHARMM. We have also developed a pathway sampling techniques to reproduce the potential of mean force (PMF) of complex chemical/biochemical reactions with reduced computational costs.

Agency
National Institute of Health (NIH)
Institute
National Heart, Lung, and Blood Institute (NHLBI)
Type
Intramural Research (Z01)
Project #
1Z01HL001051-07
Application #
6966892
Study Section
(LBC)
Project Start
Project End
Budget Start
Budget End
Support Year
7
Fiscal Year
2004
Total Cost
Indirect Cost
Name
U.S. National Heart Lung and Blood Inst
Department
Type
DUNS #
City
State
Country
United States
Zip Code
Lee, Juyong; Lee, In-Ho; Joung, InSuk et al. (2017) Finding multiple reaction pathways via global optimization of action. Nat Commun 8:15443
Lee, Juyong; Zhang, Zhong-Yuan; Lee, Jooyoung et al. (2017) Inverse Resolution Limit of Partition Density and Detecting Overlapping Communities by Link-Surprise. Sci Rep 7:12399
Huang, Jing; Simmonett, Andrew C; Pickard 4th, Frank C et al. (2017) Mapping the Drude polarizable force field onto a multipole and induced dipole model. J Chem Phys 147:161702
Tao, Peng; Hodoš?ek, Milan; Larkin, Joseph D et al. (2012) Comparison of Three Chain-of-States Methods: Nudged Elastic Band and Replica Path with Restraints or Constraints. J Chem Theory Comput 8:5035-5051
Ghysels, A; Van Neck, D; Brooks, B R et al. (2009) Normal modes for large molecules with arbitrary link constraints in the mobile block Hessian approach. J Chem Phys 130:084107
Woodcock, H Lee; Brooks, Bernard R; Pastor, Richard W (2008) Pathways and populations: stereoelectronic insights into the exocyclic torsion of 5-(hydroxymethyl)tetrahydropyran. J Am Chem Soc 130:6345-7
Miller, Benjamin T; Singh, Rishi P; Klauda, Jeffery B et al. (2008) CHARMMing: a new, flexible web portal for CHARMM. J Chem Inf Model 48:1920-9
Miller, Benjamin T; Zheng, Wenjun; Venable, Richard M et al. (2008) Langevin network model of myosin. J Phys Chem B 112:6274-81
O'Brien, Edward P; Ziv, Guy; Haran, Gilad et al. (2008) Effects of denaturants and osmolytes on proteins are accurately predicted by the molecular transfer model. Proc Natl Acad Sci U S A 105:13403-8
Zheng, Wenjun; Brooks, Bernard R; Hummer, Gerhard (2007) Protein conformational transitions explored by mixed elastic network models. Proteins 69:43-57

Showing the most recent 10 out of 37 publications