An ab initio two-center tight-binding approach to simulations of complex materials properties

Unraveling Bonding Mechanisms and Electronic Structure of Pyridine Oximes on Fe(110) Surface: Deeper Insights from DFT, Molecular Dynamics and SCC-DFT Tight Binding Simulations

The development of corrosion inhibitors with outstanding performance is a never-ending and complex process engaged in by researchers, engineers and practitioners. The computational assessment of organic corrosion inhibitors' performance is a crucial step towards the design of new task-specific materials. Herein, the electronic features, adsorption characteristics and bonding mechanisms of two pyridine oximes, namely 2-pyridylaldoxime (2POH) and 3-pyridylaldoxime (3POH), with the iron surface were investigated using molecular dynamics (MD), and self-consistent-charge density-functional tight-binding (SCC-DFTB) simulations. SCC-DFTB simulations revealed that the 3POH molecule can form covalent bonds with iron atoms in its neutral and protonated states, while the 2POH molecule can only bond with iron through its protonated form, resulting in interaction energies of -2.534, -2.007, -1.897, and -0.007 eV for 3POH, 3POH⁺, 2POH⁺, and 2POH, respectively. Projected density of states (PDOSs) analysis of pyridines-Fe(110) interactions indicated that pyridine molecules were chemically adsorbed on the iron surface. Quantum chemical calculations (QCCs) revealed that the energy gap and Hard and Soft Acids and Bases (HSAB) principles were efficient in predicting the bonding trend of the molecules investigated with an iron surface. 3POH had the lowest energy gap of 1.706 eV, followed by 3POH⁺ (2.806 eV), 2POH⁺ (3.121 eV), and 2POH (3.431 eV). In the presence of a simulated solution, MD simulation showed that the neutral and protonated forms of molecules exhibited a parallel adsorption mode on an iron surface. The excellent adsorption properties and corrosion inhibition performance of 3POH may be attributed to its low stability compared to 2POH molecules.

Theoretical study of mechanisms for the hydrolytic deamination of cytosine via steered molecular dynamic simulations

Gibbs free energy profiles of the cytosine deamination assisted by a water molecule in a discrete aqueous medium were obtained by the application of Steered Molecular Dynamic (SMD) simulations. Two pathways were considered to explain the mechanism of this process, where the water molecule attacks the C–N bond to give an intermediate (an amino–hydroxy–oxo structure in the A-path, and a hydroxy–oxo in the B-path) as the determinant step of reaction. Stationary structures along both energy profiles were analyzed at molecular dynamics level, obtaining states with higher free energies than those from electronic calculations in the gas phase and in solution described as a continuous medium. From the results obtained, the more complex A-pathway, with five steps, was kinetically the most favorable (with an endergonic reaction energy of 7.41 kcal mol⁻¹, a high barrier of 67.53 kcal mol⁻¹, and a small velocity constant k₂ = 1.80 × 10⁻³⁷ s⁻¹), concluding that the uracil base can participate in a spontaneous genetic mutation since the uracil–ammonia complex has a long lifetime of 6.10 × 10²⁷ s. This process turns out exergonic and faster when carried out in gas phase simulation or electronic calculation with a continuous medium, due to the disappearance of explicit water molecules that can compete with the assistant molecule.

Gibbs free energy profiles of the cytosine deamination assisted by a water molecule in a discrete aqueous medium were obtained by the application of Steered Molecular Dynamic (SMD) simulations.

Theoretical study of enzymatically catalyzed tautomerization of carbon acids in aqueous solution: quantum calculations and steered molecular dynamics simulations

The thermodynamics and kinetics of enzymatically assisted reactions of carbon acids were studied theoretically in this work. Quantum electronic (QE) structure calculations and steered molecular dynamics (SMD) simulations were carried out. Three 3-butenal tautomerization reactions that proceed from the β,γ-unsaturated reactant (R) to the α,β-unsaturated carbon acid product (P) and occur in two elementary steps through an intermediate (I) were studied, ignoring or including the surrounding aqueous medium in the calculations. The Gibbs free energies of activation of the R ⇆ I enolization and I ⇆ P ketonization steps were found to decrease considerably when residues simulating enzymes were introduced into these processes. Although the processes became slightly more favorable thermodynamically when the solution was included in the simulations, they became less favorable kinetically. The results from SMD simulations of these reactions were qualitatively consistent with the values we obtained using QE as well as those found by other authors in similar studies. Our simulations also allowed us to perform a detailed study of these reactions in solution.

Amide-imide tautomerism of acetohydroxamic acid in aqueous solution: quantum calculation and SMD simulations

The kinetics and the thermodynamics of the amide-imide tautomerizations of acetohydroxamic acid in aqueous solution are studied from a theoretical point of view.

Are current semiempirical methods better than force fields? A study from the thermodynamics perspective

The semiempirical Hamiltonians MNDO, AM1, PM3, RM1, PDDG/MNDO, PDDG/PM3, and SCC-DFTB, when used as part of a hybrid QM/MM scheme for the simulation of biological molecules, were compared on their abilities to reproduce experimental ensemble averages at or near room temperatures for the model system alanine dipeptide in water. Free energy surfaces in the (phi, psi) dihedral angle space, (3)J(H(N),H(alpha)) NMR dipolar coupling constants, basin populations, and peptide-water radial distribution functions (RDF) were calculated from replica exchange simulations and compared to both experiment and fully classical force field calculations using the Amber ff99SB force field. In contrast with the computational chemist's intuitive idea that the more expensive a method the better its accuracy, the ff99SB force field results were more accurate than most of the semiempirical methods, with the exception of RM1. None of the methods, however, was able to accurately reproduce the experimental data. Analysis of the results indicate that the specific QM/MM interactions have little influence on the sampling of free energy surfaces, and the differences are well explained simply by the intrinsic properties of the various QM methods.