Related Experiment Video
Updated: May 16, 2026

Synthesizing Amino Acids Modified with Reactive Carbonyls in Silico to Assess Structural Effects Using Molecular Dynamics Simulations
Published on: April 26, 2024
Charge parametrization of the DvH-c3 heme group: validation using constant-(pH,E) molecular dynamics simulations
João Henriques1, Paulo J Costa, Maria José Calhorda
1Centro de Química e Bioquímica and Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade de Lisboa, 1749-016 Lisboa, Portugal.
Abstract:
We studied the effect of using different heme group charge parametrization methods and schemes (Merz-Kollman, CHELPG, and single- and multiconformational RESP) on the quality of the results produced by the constant-(pH,E) MD method, applied to the redox titration of Desulfovibrio vulgaris Hildenborough cytochrome c(3). These new and more accurate charge sets enabled us to overcome the previously reported dependence of the method's performance on the dielectric constant, ε, assigned to the protein region. In particular, we found a systematic, clear shift of the E(mod) toward more negative values than those previously reported, in agreement with an electrostatics based reasoning. The simulations showed strong coupling between protonating/redox sites. We were also able to capture significant direct and, especially, indirect interactions between hemes, such as those mediated by histidine 67. Our results highlight the importance of having a good quantum description of the system before deriving atomic partial charges for classic force fields.
More Related Videos
08:54Vibrational Spectra of a N719-Chromophore/Titania Interface from Empirical-Potential Molecular-Dynamics Simulation, Solvated by a Room Temperature Ionic Liquid
Published on: January 25, 2020
08:0915N CPMG Relaxation Dispersion for the Investigation of Protein Conformational Dynamics on the µs-ms Timescale
Published on: April 19, 2021
Related Concept Videos
Spin–Spin Coupling: Three-Bond Coupling (Vicinal Coupling)
The extent of coupling depends on the C‑C bond length, the two H‑C‑C angles, any electron-withdrawing substituents, and the dihedral angle between the involved orbitals. The...
¹H NMR of Conformationally Flexible Molecules: Variable-Temperature NMR
Spin–Spin Coupling Constant: Overview
Qualitatively, any spin plus-half nucleus polarizes the spins of its electrons to the minus-half state. Consequently, the paired electron in the hydrogen–carbon bond must have a...
¹H NMR of Conformationally Flexible Molecules: Temporal Resolution
¹H NMR Chemical Shift Equivalence: Homotopic and Heterotopic Protons
IR Spectroscopy: Hooke's Law Approximation of Molecular Vibration
According to Hooke's law, the vibrational frequency is directly proportional to the...