Quantum Chemical Methods for Modeling Covalent Modification of Biological Thiols.

Targeted covalent inhibitor drugs require computational methods that go beyond simple molecular-mechanical force fields in order to model the chemical reactions that occur when they bind to their targets. Here, several semiempirical and density-functional theory (DFT) methods are assessed for their ability to describe the potential energy surface and reaction energies of the covalent modification of a thiol by an electrophile. Functionals such as PBE and B3LYP fail to predict a stable enolate intermediate. This is largely due to delocalization error, which spuriously stabilizes the prereaction complex, in which excess electron density is transferred from the thiolate to the electrophile. Functionals with a high-exact exchange component, range-separated DFT functionals, and variationally optimized exact exchange (i.e., the LC-B05minV functional) correct this issue to various degrees. The large gradient behavior of the exchange enhancement factor is also found to significantly affect the results, leading to the improved performance of PBE0. While ωB97X-D and M06-2X were reasonably accurate, no method provided quantitative accuracy for all three electrophiles, making this a very strenuous test of functional performance. Additionally, one drawback of M06-2X was that molecular dynamics (MD) simulations using this functional were only stable if a fine integration grid was used. The low-cost semiempirical methods, PM3, AM1, and PM7, provide a qualitatively correct description of the reaction mechanism, although the energetics is not quantitatively reliable. As a proof of concept, the potential of mean force for the addition of methylthiolate to methylvinyl ketone was calculated using quantum mechanical/molecular mechanical MD in an explicit polarizable aqueous solvent. © 2019 Wiley Periodicals, Inc.

Destruction of chemical warfare agents using metal-organic frameworks.

Chemical warfare agents containing phosphonate ester bonds are among the most toxic chemicals known to mankind. Recent global military events, such as the conflict and disarmament in Syria, have brought into focus the need to find effective strategies for the rapid destruction of these banned chemicals. Solutions are needed for immediate personal protection (for example, the filtration and catalytic destruction of airborne versions of agents), bulk destruction of chemical weapon stockpiles, protection (via coating) of clothing, equipment and buildings, and containment of agent spills. Solid heterogeneous materials such as modified activated carbon or metal oxides exhibit many desirable characteristics for the destruction of chemical warfare agents. However, low sorptive capacities, low effective active site loadings, deactivation of the active site, slow degradation kinetics, and/or a lack of tailorability offer significant room for improvement in these materials. Here, we report a carefully chosen metal-organic framework (MOF) material featuring high porosity and exceptional chemical stability that is extraordinarily effective for the degradation of nerve agents and their simulants. Experimental and computational evidence points to Lewis-acidic Zr(IV) ions as the active sites and to their superb accessibility as a defining element of their efficacy.

Prediction of (19)F NMR Chemical Shifts in Labeled Proteins: Computational Protocol and Case Study.

The structural analysis of ligand complexation in biomolecular systems is important in the design of new medicinal therapeutic agents; however, monitoring subtle structural changes in a protein's microenvironment is a challenging and complex problem. In this regard, the use of protein-based (19)F NMR for screening low-molecular-weight molecules (i.e., fragments) can be an especially powerful tool to aid in drug design. Resonance assignment of the protein's (19)F NMR spectrum is necessary for structural analysis. Here, a quantum chemical method has been developed as an initial approach to facilitate the assignment of a fluorinated protein's (19)F NMR spectrum. The epigenetic "reader" domain of protein Brd4 was taken as a case study to assess the strengths and limitations of the method. The overall modeling protocol predicts chemical shifts for residues in rigid proteins with good accuracy; proper accounting for explicit solvation of fluorinated residues by water is critical.

Empirical and theoretical insights into the structural features and host-guest chemistry of M8L4 tube architectures.

We demonstrate a general method for the construction of M8L4 tubular complexes via subcomponent self-assembly, starting from Cu(I) or Ag(I) precursors together with suitable elongated tetraamine and 2-formylpyridine subcomponents. The tubular architectures were often observed as equilibrium mixtures of diastereomers having two different point symmetries (D2d or D2 ⇄ D4) in solution. The equilibria between diastereomers were influenced through variation in ligand length, substituents, metal ion identity, counteranion, and temperature. In the presence of dicyanoaurate(I) and Au(I), the D4-symmetric hosts were able to bind linear Au(Au(CN)2)2(-) (with two different configurations) as the best-fitting guest. Substitution of dicyanoargentate(I) for dicyanoaurate(I) resulted in the formation of Ag(Au(CN)2)2(-) as the optimal guest through transmetalation. Density functional theory was employed to elucidate the host-guest chemistries of the tubes.

Design of a metal-organic framework with enhanced back bonding for separation of N2 and CH4.

Gas separations with porous materials are economically important and provide a unique challenge to fundamental materials design, as adsorbent properties can be altered to achieve selective gas adsorption. Metal-organic frameworks represent a rapidly expanding new class of porous adsorbents with a large range of possibilities for designing materials with desired functionalities. Given the large number of possible framework structures, quantum mechanical computations can provide useful guidance in prioritizing the synthesis of the most useful materials for a given application. Here, we show that such calculations can predict a new metal-organic framework of potential utility for separation of dinitrogen from methane, a particularly challenging separation of critical value for utilizing natural gas. An open V(II) site incorporated into a metal-organic framework can provide a material with a considerably higher enthalpy of adsorption for dinitrogen than for methane, based on strong selective back bonding with the former but not the latter.

Predicting paramagnetic 1H NMR chemical shifts and state-energy separations in spin-crossover host-guest systems.

The behaviour of metal-organic cages upon guest encapsulation can be difficult to elucidate in solution. Paramagnetic metal centres introduce additional dispersion of signals that is useful for characterisation of host-guest complexes in solution using nuclear magnetic resonance (NMR). However, paramagnetic centres also complicate spectral assignment due to line broadening, signal integration error, and large changes in chemical shifts, which can be difficult to assign even for known compounds. Quantum chemical predictions can provide information that greatly facilitates the assignment of NMR signals and identification of species present. Here we explore how the prediction of paramagnetic NMR spectra may be used to gain insight into the spin crossover (SCO) properties of iron(II)-based metal organic coordination cages, specifically examining how the structure of the local metal coordination environment affects SCO. To represent the tetrahedral metal-organic cage, a model system is generated by considering an isolated metal-ion vertex: fac-ML3(2+) (M = Fe(II), Co(II); L = N-phenyl-2-pyridinaldimine). The sensitivity of the (1)H paramagnetic chemical shifts to local coordination environments is assessed and utilised to shed light on spin crossover behaviour in iron complexes. Our data indicate that expansion of the metal coordination sphere must precede any thermal SCO. An attempt to correlate experimental enthalpies of SCO with static properties of bound guests shows that no simple relationship exists, and that effects are likely due to nuanced dynamic response to encapsulation.

Ab Initio Extension of the AMOEBA Polarizable Force Field to Fe(2.).

We extend the AMOEBA polarizable molecular mechanics force field to the Fe(2+) cation in its singlet, triplet, and quintet spin states. Required parameters are obtained either directly from first principles calculations or optimized so as to reproduce corresponding interaction energy components in a hexaaquo environment derived from quantum mechanical energy decomposition analyses. We assess the importance of the damping of point-dipole polarization at short distance as well as the influence of charge-transfer for metal-water interactions in hydrated Fe(2+); this analysis informs the selection of model systems employed for parametrization. We validate our final Fe(2+) model through comparison of molecular dynamics (MD) simulations to available experimental data for aqueous ferrous ion in its quintet electronic ground state.