Journal
JOURNAL OF CHEMICAL THEORY AND COMPUTATION
Volume 14, Issue 2, Pages 784-798Publisher
AMER CHEMICAL SOC
DOI: 10.1021/acs.jctc.7b01236
Keywords
-
Funding
- NSF [MCB-1452464, OCI-1053575]
- NIH [R01-GM123169]
- U.S. Department of Energy (DOE) Office of Science, Biological and Environmental Research, Subsurface Biogeochemical Research (SBR) Program through the Mercury Science Focus Area Program (SFA) at Oak Ridge National Laboratory (ORNL)
- U.S. DOE [DE-AC05-00OR22725]
- Office of Science of the U.S. DOE [DE-AC02-05CH11231]
- Direct For Mathematical & Physical Scien
- Division Of Physics [1205878] Funding Source: National Science Foundation
- Div Of Molecular and Cellular Bioscience
- Direct For Biological Sciences [1452464] Funding Source: National Science Foundation
Ask authors/readers for more resources
Corrinoid cofactors such as cobalamin are used by many enzymes and are essential for most living organisms. Therefore, there is broad interest in investigating cobalamin protein-interactions with molecular dynamics simulations. Previously developed parameters for cobalamins are based mainly on crystal structure data. Here, we report CHARMM-compatible force field parameters for several corrinoids developed from quantum mechanical calculations. We provide parameters for corrinoids in three oxidation states, Co3+, Co2+, and Co1+ and with various axial ligands. Lennard-Jones parameters for the cobalt center in the Co(II) and Co(I) states were optimized using a helium atom probe, and partial atomic charges were obtained with a combination of natural population analysis (NPA) and restrained electrostatic potential (RESP) fitting approaches. The Force Field Toolkit was used to optimize all bonded terms. The resulting parameters, determined solely from calculations of cobalamin alone or in water, were then validated by assessing their agreement with density functional theory geometries and by analyzing molecular dynamics simulation trajectories of several corrinoid proteins for which X-ray crystal structures are available. In each case, we obtained excellent agreement with the reference data. In comparison to previous CHARMM-compatible parameters for cobalamin, we observe a better agreement for the fold angle and lower RMSD in the cobalamin binding site. The approach described here is readily adaptable for developing CHARMM-compatible force-field parameters for other corrinoids or large biomolecules.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available