Molecular dynamics simulation of entropy driven ligand escape process in heme pocket

Generating High-Precision Force Fields for Molecular Dynamics Simulations to Study Chemical Reaction Mechanisms Using Molecular Configuration Transformer

Theoretical studies on chemical reaction mechanisms have been crucial in organic chemistry. Traditionally, calculating the manually constructed molecular conformations of transition states for chemical reactions using quantum chemical calculations is the most commonly used method. However, this way is heavily dependent on individual experience and chemical intuition. In our previous study, we proposed a research paradigm that used enhanced sampling in molecular dynamics simulations to study chemical reactions. This approach can directly simulate the entire process of a chemical reaction. However, the computational speed limited the use of high-precision potential energy functions for simulations. To address this issue, we presented a scheme for training high-precision force fields for molecular modeling using a previously developed graph-neural-network-based molecular model, molecular configuration transformer. This potential energy function allowed for highly accurate simulations at a low computational cost, leading to more precise calculations of the mechanism of chemical reactions. We applied this approach to study a Claisen rearrangement reaction and a carbonyl insertion reaction catalyzed by manganese.

Interfacial water effect on cooperativity and signal communication in Scapharca dimeric hemoglobin

Cooperativity is important in controlling the biological functions of allosteric proteins.

Internal water and microsecond dynamics in myoglobin

Myoglobin (Mb) binds diatomic ligands, like O₂, CO, and NO, in a cavity that is only transiently accessible. Crystallography and molecular simulations show that the ligands can migrate through an extensive network of transiently connected cavities but disagree on the locations and occupancy of internal hydration sites. Here, we use water ²H and ¹⁷O magnetic relaxation dispersion (MRD) to characterize the internal water molecules in Mb under physiological conditions. We find that equine carbonmonoxy Mb contains 4.5 ± 1.0 ordered internal water molecules with a mean survival time of 5.6 ± 0.5 μs at 25 °C. The likely locations of these water molecules are the four polar hydration sites, including one of the xenon-binding cavities, that are fully occupied in all high-resolution crystal structures of equine Mb. The finding that water escapes from these sites, located 17–31 Å apart in the protein, on the same μs time scale suggests a global exchange mechanism. We propose that this mechanism involves transient penetration of the protein by H-bonded water chains. Such a mechanism could play a functional role by eliminating trapped ligands. In addition, the MRD results indicate that 2 or 3 of the 11 histidine residues of equine Mb undergo intramolecular hydrogen exchange on a μs time scale.

Ligand migration in myoglobin: a combined study of computer simulation and x-ray crystallography

A ligand-migration mechanism of myoglobin was studied by a multidisciplinary approach that used x-ray crystallography and molecular dynamics simulation. The former revealed the structural changes of the protein along with the ligand migration, and the latter provided the statistical ensemble of protein conformations around the thermal average. We developed a novel computational method, homogeneous ensemble displacement, and generated the conformational ensemble of ligand-detached species from that of ligand-bound species. The thermally averaged ligand-protein interaction was illustrated in terms of the potential of mean force. Although the structural changes were small, the presence of the ligand molecule in the protein matrix significantly affected the 3D scalar field of the potential of mean force, in accordance with the self-opening model proposed in the previous x-ray study.

Force-Field Development and Molecular Dynamics Simulations of Ferrocene–Peptide Conjugates as a Scaffold for Hydrogenase Mimics

The increasing importance of hydrogenase enzymes in the new energy research field has led us to examine the structure and dynamics of potential hydrogenase mimics, based on a ferrocene-peptide scaffold, using molecular dynamics (MD) simulations. To enable this MD study, a molecular mechanics force field for ferrocene-bearing peptides was developed and implemented in the CHARMM simulation package, thus extending the usefulness of the package into peptide-bioorganometallic chemistry. Using the automated frequency-matching method (AFMM), optimized intramolecular force-field parameters were generated through quantum chemical reference normal modes. The partial charges for ferrocene were derived by fitting point charges to quantum-chemically computed electrostatic potentials. The force field was tested against experimental X-ray crystal structures of dipeptide derivatives of ferrocene-1,1'-dicarboxylic acid. The calculations reproduce accurately the molecular geometries, including the characteristic C2-symmetrical intramolecular hydrogen-bonding pattern, that were stable over 0.1 micros MD simulations. The crystal packing properties of ferrocene-1-(D)alanine-(D)proline-1'-(D)alanine-(D)proline were also accurately reproduced. The lattice parameters of this crystal were conserved during a 0.1 micros MD simulation and match the experimental values almost exactly. Simulations of the peptides in dichloromethane are also in good agreement with experimental NMR and circular dichroism (CD) data in solution. The developed force field was used to perform MD simulations on novel, as yet unsynthesized peptide fragments that surround the active site of [Ni-Fe] hydrogenase. The results of this simulation lead us to propose an improved design for synthetic peptide-based hydrogenase models. The presented MD simulation results of metallocenes thereby provide a convincing validation of our proposal to use ferrocene-peptides as minimal enzyme mimics.

Distal charge transport in peptides

Biological systems often transport charges and reactive processes over substantial distances. Traditional models of chemical kinetics generally do not describe such extreme distal processes. In this Review, an atomistic model for a distal transport of information, which was specifically developed for peptides, is considered. Chemical reactivity is taken as the result of distal effects based on two-step bifunctional kinetics involving unique, very rapid motional properties of peptides in the subpicosecond regime. The bifunctional model suggests highly efficient transport of charge and reactivity in an isolated peptide over a substantial distance; conversely, a very low efficiency in a water environment was found. The model suggests ultrafast transport of charge and reactivity over substantial molecular distances in a peptide environment. Many such domains can be active in a protein.

Selectivity principle of the ligand escape process from a two-gate tunnel in myoglobin: molecular dynamics simulation

We proposed a selectivity principle for the ligand escape process from two fluctuating bottlenecks in a cavity with a multigate inside a myoglobin pocket. Our previous analytical theory proposed a fluctuating bottleneck model for a Brownian particle passing through two gates on a cavity surface of an enzyme protein and has determined the escape rate in terms of the time-dependent gate function and the competition effect. It illustrated that with two (or more than two) gates on a cavity surface the gate modulation, which is controlled by protein fluctuation, dominates the ligand escape pathway. We have performed a molecular dynamics simulation to investigate the selectivity principle of the ligand escape process from two-gate tunnel in myoglobin. The simulation results confirm our theoretical conjecture. It indicates that the escape process is actually entropy driven, and the ligand escape pathway is chosen via the gate modulation. This suggests an interesting intrinsic property, that is, the oxymyoglobin tertiary structure is favorable to the departure of the ligand from one direction rather than through a biased random walk.

Modeling heme proteins using atomistic simulations

Heme proteins are found in all living organisms, and perform a wide variety of tasks ranging from electron transport, to the oxidation of organic compounds, to the sensing and transport of small molecules. In this work we review the application of classical and quantum-mechanical atomistic simulation tools to the investigation of several relevant issues in heme proteins chemistry: (i) conformational analysis, ligand migration, and solvation effects studied using classical molecular dynamics simulations; (ii) electronic structure and spin state energetics of the active sites explored using quantum-mechanics (QM) methods; (iii) the interaction of heme proteins with small ligands studied through hybrid quantum mechanics-molecular mechanics (QM-MM) techniques; (iv) and finally chemical reactivity and catalysis tackled by a combination of quantum and classical tools.