A theoretical study of the thermodynamics and kinetics of small organosulfur compounds

Reliability of computed molecular structures

When the structures of 1342 molecules are optimized by 30 methods and 7 basis sets, there appear 289 (21.54%) problematic molecules and 112 (8.35%) failed ones. When 278 problematic molecules are compared, the best methods are BHandH and LC-wPBE, while B97D, BP86, HFS, VSXC, and HCTH are very unreliable. When 179 problematic molecules are computed with larger basis sets, the smallest mean absolute deviation (MAD) of bond angle (2.3°) is shown by QCISD(T)/cc-pVTZ, while the smallest MAD of bond length (0.021 Å), the best SUM1 (4.9 unit), and the best SUM2 (2.4 unit) are shown by DSDPBEP86(Full), DSDPBEP86, PBE1PBE-D3, MP2, and MP2(Full) in combination with aug-cc-pVQZ, cc-pVQZ, Def2QZVP, Def2TZVPP, and/or 6-311++G(3df,3pd). Very large basis sets, for example, larger than cc-pVTZ usually have to be used to obtain very good structures and the performances of many density-functional theory methods are encouraging. The best results may be the limit of modern computational chemistry.

Quantum Mechanical Methods Predict Accurate Thermodynamics of Biochemical Reactions

Thermodynamics plays a crucial role in regulating the metabolic processes in all living organisms. Accurate determination of biochemical and biophysical properties is important to understand, analyze, and synthetically design such metabolic processes for engineered systems. In this work, we extensively performed first-principles quantum mechanical calculations to assess its accuracy in estimating free energy of biochemical reactions and developed automated quantum-chemistry (QC) pipeline (https://appdev.kbase.us/narrative/45710) for the prediction of thermodynamics parameters of biochemical reactions. We benchmark the QC methods based on density functional theory (DFT) against different basis sets, solvation models, pH, and exchange-correlation functionals using the known thermodynamic properties from the NIST database. Our results show that QC calculations when combined with simple calibration yield a mean absolute error in the range of 1.60-2.27 kcal/mol for different exchange-correlation functionals, which is comparable to the error in the experimental measurements. This accuracy over a diverse set of metabolic reactions is unprecedented and near the benchmark chemical accuracy of 1 kcal/mol that is usually desired from DFT calculations.

Desulfination by 2'-hydroxybiphenyl-2-sulfinate desulfinase proceeds via electrophilic aromatic substitution by the cysteine-27 proton

Density functional theory shows that the rate-limiting desulfination step in biodesulfurization involves concerted electrophilic substitution with the Cys-27 proton.

Biodesulfurization is an attractive option for enzymatically removing sulfur from the recalcitrant thiophenic derivatives that comprise the majority of organosulfur compounds remaining in hydrotreated petroleum products. Desulfurization in the bacteria Rhodococcus erythropolis follows a four-step pathway culminating in C–S bond cleavage in the 2′-hydroxybiphenyl-2-sulfinate (HBPS) intermediate to yield 2-hydroxybiphenyl and bisulfite. The reaction, catalyzed by 2′-hydroxybiphenyl-2-sulfinate desulfinase (DszB), is the rate-limiting step and also the least understood, as experimental evidence points to a mechanism unlike that of other desulfinases. On the basis of structural and biochemical evidence, two possible mechanisms have been proposed: nucleophilic addition and electrophilic aromatic substitution. Density functional theory calculations showed that electrophilic substitution by a proton is the lower energy pathway and is consistent with previous kinetic and site-directed mutagenesis studies. C27 transfers its proton to HBPS, leading directly to the release of SO₂ without the formation of a carbocation intermediate. The H60–S25 dyad stabilizes the transition state by withdrawing the developing negative charge on cysteine. Establishing the desulfination mechanism and specific role of active site residues, accomplished in this study, is essential to protein engineering efforts to increase DszB catalytic activity, which is currently too low for industrial-scale application.

Reactivity of electrophilic chlorine atoms due to σ-holes: a mechanistic assessment of the chemical reduction of a trichloromethyl group by sulfur nucleophiles

σ-Holes are shown to promote the electrophilic behavior of chlorine atoms in a trichloromethyl group when bound to an electron-withdrawing moiety. A halogen bond-type non-covalent interaction between a chlorine atom and a negatively charged sulfur atom takes place, causing the abstraction of such a chlorine atom while leaving a carbanion, subsequently driving the chemical reduction of the trichloromethyl group to a sulfide in a stepwise process. The mechanism for the model reaction of trichloromethyl pyrimidine 1 with thiophenolate and thiophenol to yield phenylsulfide 4 was followed through ¹H-NMR and studied using DFT transition state calculations, and the energy profile for this transformation is fully discussed. MP2 calculations of the electrostatic potential were performed for a series of trichloromethyl compounds in order to assess the presence of σ-holes and quantify them by means of the maximum surface electrostatic potential. Such calculations showed that the chlorine atoms behave as electrophilic leaving groups toward a nucleophilic attack, opening a new possibility in the synthetic chemistry of the trichloromethyl group.

Automatic mechanism generation for pyrolysis of di-tert-butyl sulfide

The automated Reaction Mechanism Generator (RMG), using rate parameters derived from ab initio CCSD(T) calculations, is used to build reaction networks for the thermal decomposition of di-tert-butyl sulfide.

A kinetic and thermochemical database for organic sulfur and oxygen compounds

Potential energy surfaces and reaction kinetics were calculated for reactions involving sulfur and oxygen, which are potentially relevant in combustion and desulfurization chemistry.

Post-translational modification in the gas phase: mechanism of cysteine S-nitrosylation via ion-molecule reactions

The gas-phase mechanism of S-nitrosylation of thiols was studied in a quadrupole ion trap mass spectrometer. This was done via ion-molecule reactions of protonated cysteine and many of its derivatives and other thiol ions with neutral tert-butyl nitrite or nitrous acid. Our results showed that the presence of the carboxylic acid functional group, -COOH, in the vicinity of the thiol group is essential for the gas-phase nitrosylation of thiols. When the carboxyl proton is replaced by a methyl group (cysteine methyl ester) no nitrosylation was observed. Other thiols lacking a carboxylic acid functional group displayed no S-nitrosylation, strongly suggesting that the carboxyl hydrogen plays a key role in the nitrosylation process. These results are in excellent agreement with a solution-phase mechanism proposed by Stamler et al. (J. S. Stamler, E. J. Toone, S. A. Lipton, N. J. Sucher. Neuron 1997, 18, 691-696) who suggested a catalytic role for the carboxylic acid group adjacent to cysteine residues and with later additions by Ascenzi et al. (P. Ascenzi, M. Colasanti, T. Persichini, M. Muolo, F. Polticelli, G. Venturini, D. Bordo, M. Bolognesi. Biol. Chem. 2000, 381, 623-627) who postulated that the presence of the carboxyl in the cysteine microenvironment in proteins is crucial for S-nitrosylation. A concerted mechanism for the gas-phase S-nitrosylation was proposed based on our results and was further studied using theoretical calculations. Our calculations showed that this proposed pathway is exothermic by 44.0 kJ mol(-1). This is one of the few recent examples when a gas-phase mechanism matches one in solution.