pDynamo3 Molecular Modeling and Simulation Program.

pDynamo3 is the first formal version of the Dynamo molecular modeling and simulation library that is written in Python 3. It follows from the previous pDynamo versions written in Python 2, the first of which was released in 2007. Both pDynamo and its predecessor, fDynamo, were designed with the aim of providing easy-to-use and flexible frameworks for performing molecular simulations at an atomistic level with a special emphasis on those employing hybrid quantum chemical and molecular mechanical potential methods. Although the use of pDynamo3 is quite similar to that of pDynamo2, it has added significant new capability and also undergone extensive restructuring that will make it much easier to extend with new functionality. The pDynamo3 code is issued under the GNU general public license at https://github.com/pdynamo/pDynamo3 with additional information on the pDynamo website https:www.pdynamo.org.

Density Fitting for QC/MM Interactions.

Density fitting is a standard technique in quantum chemistry as it helps to accelerate certain parts of a calculation, such as the computation of the electron repulsion energy, without significant loss of accuracy. This paper explores the effectiveness of this technique when it is extended to treat interactions with external electrostatic potentials, in particular those that arise from the environment in hybrid quantum chemical/molecular mechanical calculations. It is found that fitted densities are able to well reproduce the energies, forces, and properties obtained using unfitted densities, as long as a suitable operator is employed for the fitting. It is expected that this precision would be improved by the development of basis sets specifically designed to treat these types of interactions and that the use of this approximation could lead to substantial speed-ups in large hybrid potential simulations.

Theoretical characterization of the shikimate 5-dehydrogenase reaction from Mycobacterium tuberculosis by hybrid QC/MM simulations and quantum chemical descriptors.

In this study, we have investigated the enzyme shikimate 5-dehydrogenase from the causative agent of tuberculosis, Mycobacterium tuberculosis. We have employed a mixture of computational techniques, including molecular dynamics, hybrid quantum chemical/molecular mechanical potentials, relaxed surface scans, quantum chemical descriptors and free-energy simulations, to elucidate the enzyme's reaction pathway. Overall, we find a two-step mechanism, with a single transition state, that proceeds by an energetically uphill hydride transfer, followed by an energetically downhill proton transfer. Our mechanism and calculated free energy barrier for the reaction, 64.9 kJ mol- 1, are in good agreement with those predicted from experiment. An analysis of quantum chemical descriptors along the reaction pathway indicated a possibly important, yet currently unreported, role of the active site threonine residue, Thr65.

Modeling the Hydrolysis of Iron-Sulfur Clusters.

Iron-sulfur (FeS) clusters are essential metal cofactors involved in a wide variety of biological functions. Their catalytic efficiency, biosynthesis, and regulation depend on FeS stability in aqueous solution. Here, molecular modeling is used to investigate the hydrolysis of an oxidized (ferric) mononuclear FeS cluster by bare dissociation and water substitution mechanisms in neutral and acidic solution. First, approximate electronic structure descriptions of FeS reactions by density functional theory are validated against high-level wave function CCSD(T) calculations. Solvation contributions are included by an all-atom model with hybrid quantum chemical/molecular mechanical (QM/MM) potentials and enhanced sampling molecular dynamics simulations. The free energy profile obtained for FeS cluster hydrolysis indicates that the hybrid functional M06 together with an implicit solvent correction capture the most important aspects of FeS cluster reactivity in aqueous solution. Then, 20 reaction channels leading to two consecutive Fe-S bond ruptures were explored with this calibrated model. For all protonation states, nucleophilic substitution with concerted bond breaking and forming to iron is the preferred mechanism, both kinetic and thermodynamically. In neutral solution, proton transfer from water to the sulfur leaving group is also concerted. Dissociative reactions show higher barriers and will not be relevant for FeS reactivity when exposed to solvent. These hydrolysis mechanisms may help to explain the stability and catalytic mechanisms of FeS clusters of multiple sizes and proteins.

Electrocatalytic Hydrogen Evolution with a Cobalt Complex Bearing Pendant Proton Relays: Acid Strength and Applied Potential Govern Mechanism and Stability.

[Co(bapbpy)Cl]+ (bapbpy: 6,6'-bis(2-aminopyridyl)-2,2'-bipyridine) is a polypyridyl cobalt(II) complex bearing both a redox-active bipyridine ligand and pendant proton relays. This compound catalyzes electro-assisted H2 evolution in DMF with distinct mechanisms depending on the strength of the acid used as the proton source (pKa values ranging from 3.4 to 13.5 in DMF) and the applied potential. Electrochemical studies combining cyclic voltammetry and bulk electrolysis measurements enabled one to bring out four distinct catalytic processes. Where applicable, relevant kinetic information were obtained using either foot-of-the-wave analysis (FOWA) or analytical treatment of bulk electrolysis experiments. Among the different catalytic pathways identified in this study, a clear relationship between the catalyst performances and stability was evidenced. These results draw attention to a number of interesting considerations and may help in the development of future adequately designed catalysts.

Negative Impact of Carbapenem Methylation on the Reactivity of ß-Lactams for Cysteine Acylation as Revealed by Quantum Calculations and Kinetic Analyses.

The Enterococcus faecium l,d-transpeptidase (Ldtfm) mediates resistance to most ß-lactam antibiotics in this bacterium by replacing classical peptidoglycan polymerases. The catalytic Cys of Ldtfm is rapidly acylated by ß-lactams belonging to the carbapenem class but not by penams or cephems. We previously reported quantum calculations and kinetic analyses for Ldtfm and showed that the inactivation profile is not determined by differences in drug binding (KD [equilibrium dissociation constant] values in the 50 to 80 mM range). In this study, we analyzed the reaction of a Cys sulfhydryl with various ß-lactams in the absence of the enzyme environment in order to compare the intrinsic reactivity of drugs belonging to the penam, cephem, and carbapenem classes. For this purpose, we synthesized cyclic Cys-Asn (cCys-Asn) to generate a soluble molecule with a sulfhydryl closely mimicking a cysteine in a polypeptide chain, thereby avoiding free reactive amino and carboxyl groups. Computational studies identified a thermodynamically favored pathway involving a concerted rupture of the ß-lactam amide bond and formation of an amine anion. Energy barriers indicated that the drug reactivity was the highest for nonmethylated carbapenems, intermediate for methylated carbapenems and cephems, and the lowest for penams. Electron-withdrawing groups were key reactivity determinants by enabling delocalization of the negative charge of the amine anion. Acylation rates of cCys-Asn determined by spectrophotometry revealed the same order in the reactivity of ß-lactams. We concluded that the rate of Ldtfm acylation is largely determined by the ß-lactam reactivity with one exception, as the enzyme catalytic pocket fully compensated for the detrimental effect of carbapenem methylation.

Proton reduction reaction catalyzed by homoleptic nickel bis-1,2-dithiolate complexes: Experimental and theoretical mechanistic investigations.

A series of homoleptic monoanionic nickel dithiolene complexes [Ni(bdt)2](NBu4), [Ni(tdt)2](NBu4), and [Ni(mnt)2](NBu4) containing the ligands benzene-1,2-dithiolate (bdt2-), toluene-3,4-dithiolate (tdt2-) and maleonitriledithiolate (mnt2-), respectively, have been employed as electrocatalysts in the hydrogen evolution reaction with trifluoroacetic acid as proton source in acetonitrile. All complexes were active catalysts with TONs reaching 113, 158 and 6 for [Ni(bdt)2](NBu4), [Ni(tdt)2](NBu4), and [Ni(mnt)2](NBu4), respectively. Faradaic yield for hydrogen evolution reaction reaches 88 % for 2- , which also displays the minimal overpotential requirement value (467 mV) within the series. Two pathways for H2 evolution can be hypothesized that differ on on the sequence of protonation and reduction steps. DFT calculations are in agreement with experimental data and indicate that protonation at sulfur follows reduction to the dianion. Hydrogen evolves from the direduced-diprotonated form via a highly distorted nickel hydride intermediate. The effects of acid strength and concentration in the hydrogen-evolving mechanism are also discussed.

Norovirus RNA-dependent RNA polymerase: A computational study of metal-binding preferences.

Norovirus (NV) RNA-dependent RNA polymerase (RdRP) is essential for replicating the genome of the virus, which makes this enzyme a key target for the development of antiviral agents against NV gastroenteritis. In this work, a complex of NV RdRP bound to manganese ions and an RNA primer-template duplex was investigated using X-ray crystallography and hybrid quantum chemical/molecular mechanical simulations. Experimentally, the complex crystallized in a tetragonal crystal form. The nature of the primer/template duplex binding in the resulting structure indicates that the complex is a closed back-tracked state of the enzyme, in which the 3'-end of the primer occupies the position expected for the post-incorporated nucleotide before translocation. Computationally, it is found that the complex can accept a range of divalent metal cations without marked distortions in the active site structure. The highest binding energy is for copper, followed closely by manganese and iron, and then by zinc, nickel, and cobalt. Proteins 2017; 85:1435-1445. © 2017 Wiley Periodicals, Inc.

An Algorithm for Adaptive QC/MM Simulations.

An algorithm is proposed for the simulation of molecular systems with hybrid quantum chemical (QC) and molecular mechanical (MM) potentials that permits the adaptive partitioning of the atoms in the system between QC and MM regions. In contrast to existing methods, the algorithm requires only a single QC calculation of the QC/MM system per energy calculation and yet has consistent energy and forces, which makes it suitable for geometry optimizations and molecular dynamics calculations within the microcanonical ensemble, in addition to other types of simulation. This article describes the algorithm and its implementation, presents some simple test cases with both semiempirical and density functional theory QC/MM potentials, and discusses perspectives for future work.

Catalytic Mechanism of Peptidoglycan Deacetylase: A Computational Study.

Bacterial peptidoglycan deacetylase enzymes are potentially important targets for the design of new drugs. In pathogenic bacteria, they modify cell-wall peptidoglycan by removing the acetyl group, which makes the bacteria more resistant to the host's immune response and other forms of attack, such as degradation by lysozyme. In this study, we have investigated the mechanism of reaction of acetyl removal from a model substrate, the N-acetylglucosamine/N-acetylmuramic acid dimer, by peptidogylcan deacetylase from Helicobacter pylori. For this, we employed a range of computational approaches, including molecular docking, Poisson-Boltzmann electrostatic pKa calculations, molecular dynamics simulations, and hybrid quantum chemical/molecular mechanical potential calculations, in conjunction with reaction-path-finding algorithms. The active site of this enzyme is in a region of highly negative electrostatic potential and contains a zinc dication with a bound water molecule. In the docked enzyme-substrate complex, our pKa calculations indicate that in the most stable protonation states of the active site the zinc-bound water molecule is in its hydroxide form and that the adjacent histidine residue, His247, is doubly protonated. In addition, there are one or two excess protons, with the neighboring aspartate residues, Asp12 and/or Asp199, being protonated. Overall, we find five classes of feasible reaction mechanisms, with the favored mechanism depending heavily on the protonation state of the active site. In the major one-excess-proton form, the mechanism with the lowest barrier (84 kJ mol-1) involves an initial protonation of the substrate nitrogen, followed by nucleophilic attack of the zinc-bound hydroxide and rupture of the substrate's carbon-nitrogen bond. However, in the minor two-excess-proton form, four mechanisms are almost equienergetic (83-86 kJ mol-1), comprising both those that start with nitrogen protonation and those in which nucleophilic attack by hydroxide occurs first.

Hybrid Potential Simulation of the Acylation of Enterococcus faecium l,d-Transpeptidase by Carbapenems.

The l,d-transpeptidases, Ldts, catalyze peptidoglycan cross-linking in ß-lactam-resistant mutant strains of several bacteria, including Enterococcus faecium and Mycobacterium tuberculosis. Although unrelated to the essential d,d-transpeptidases, which are inactivated by the ß-lactam antibiotics, they are nevertheless inhibited by the carbapenem antibiotics, making them potentially useful targets in the treatment of some important diseases. In this work, we have investigated the acylation mechanism of the Ldt from E. faecium by the carbapenem, ertapenem, using computational techniques. We have employed molecular dynamics simulations in conjunction with QC/MM hybrid potential calculations to map out possible reaction paths. We have focused on determining the following: (i) the protonation state of the nucleophilic cysteine of the enzyme when it attacks; (ii) whether nucleophilic attack and ß-lactam ring-opening are concerted or stepwise, the latter occurring via an oxyanion intermediate; and (iii) the identities of the proton acceptors at the beginning and end of the reaction. Overall, we note that there is considerable plasticity in the mechanisms, owing to the significant flexibility of the enzyme, but find that the preferred pathways are ones in which nucleophilic attack of cysteine thiolate is concerted with ß-lactam ring-opening.

Experimental and Theoretical Insight into Electrocatalytic Hydrogen Evolution with Nickel Bis(aryldithiolene) Complexes as Catalysts.

A series of neutral and monoanionic nickel dithiolene complexes with p-methoxyphenyl-substituted 1,2-dithiolene ligands have been prepared and characterized with physicochemical methods. Two of the complexes, the monoanion of the symmetric [Ni{S2C2(Ph-p-OCH3)2}2] (3(-)) with NBu4(+) as a counterion and the neutral asymmetric [Ni{S2C2(Ph)(Ph-p-OCH3)}2] (2), have been structurally characterized by single-crystal X-ray crystallography. All complexes have been employed as proton-reducing catalysts in N,N-dimethylformamide with trifluoroacetic acid as the proton source. The complexes are active catalysts with good faradaic yields, reaching 83% for 2 but relatively high overpotential requirements (0.91 and 1.55 V measured at the middle of the catalytic wave for two processes observed depending on the different routes of the mechanism). The similarity of the experimental data regardless of whether the neutral or anionic form of the complexes is used indicates that the neutral form acts as a precatalyst. On the basis of detailed density functional theory calculations, the proposed mechanism reveals two different main routes after protonation of the dianion of the catalyst in accordance with the experimental data, indicating the role of the concentration of the acid and the influence of the methoxy groups. Protonation at sulfur seems be more favorable than that at the metal, which is in marked contrast with the catalytic mechanism proposed for analogous cobalt dithiolene complexes.

Artificial Hydrogenases Based on Cobaloximes and Heme Oxygenase.

The insertion of cobaloxime catalysts in the heme-binding pocket of heme oxygenase (HO) yields artificial hydrogenases active for H2 evolution in neutral aqueous solutions. These novel biohybrids have been purified and characterized by using UV/visible and EPR spectroscopy. These analyses revealed the presence of two distinct binding conformations, thereby providing the cobaloxime with hydrophobic and hydrophilic environments, respectively. Quantum chemical/molecular mechanical docking calculations found open and closed conformations of the binding pocket owing to mobile amino acid residues. HO-based biohybrids incorporating a {Co(dmgH)2 } (dmgH2 =dimethylglyoxime) catalytic center displayed up to threefold increased turnover numbers with respect to the cobaloxime alone or to analogous sperm whale myoglobin adducts. This study thus provides a strong basis for further improvement of such biohybrids, using well-designed modifications of the second and outer coordination spheres, through site-directed mutagenesis of the host protein.

Ferric-Thiolate Bond Dissociation Studied with Electronic Structure Calculations.

The stability and reactivity of iron-sulfur clusters are fundamental properties for the biological function of these prosthetic groups. Here, we investigate the ferric-thiolate bond dissociation of model iron-sulfur tetrahedral complexes with high-level ab initio multiconfigurational electronic structure calculations. We find that the reaction mechanism is homolytic with a spin-crossing from the sextet state in the reactant to quartet state in the product. We also compare several density functionals and semiempirical configuration interaction with the high-level ab initio results to find an accurate but computationally more efficient method to describe the reaction. The functionals M06 and those based on the OPTX exchange functional show the best performance and may reasonably describe the various electron correlation effects involved in ferric-thiolate bond dissociation.

Pcetk: A pDynamo-based Toolkit for Protonation State Calculations in Proteins.

Pcetk (a pDynamo-based continuum electrostatic toolkit) is an open-source, object-oriented toolkit for the calculation of proton binding energetics in proteins. The toolkit is a module of the pDynamo software library, combining the versatility of the Python scripting language and the efficiency of the compiled languages, C and Cython. In the toolkit, we have connected pDynamo to the external Poisson-Boltzmann solver, extended-MEAD. Our goal was to provide a modern and extensible environment for the calculation of protonation states, electrostatic energies, titration curves, and other electrostatic-dependent properties of proteins. Pcetk is freely available under the CeCILL license, which is compatible with the GNU General Public License. The toolkit can be found on the Web at the address http://github.com/mfx9/pcetk. The calculation of protonation states in proteins requires a knowledge of pKa values of protonatable groups in aqueous solution. However, for some groups, such as protonatable ligands bound to protein, the pKa aq values are often difficult to obtain from experiment. As a complement to Pcetk, we revisit an earlier computational method for the estimation of pKa aq values that has an accuracy of ± 0.5 pKa-units or better. Finally, we verify the Pcetk module and the method for estimating pKa aq values with different model cases.

Force-induced chemical reactions on the metal centre in a single metalloprotein molecule.

Metalloproteins play indispensable roles in biology owing to the versatile chemical reactivity of metal centres. However, studying their reactivity in many metalloproteins is challenging, as protein three-dimensional structure encloses labile metal centres, thus limiting their access to reactants and impeding direct measurements. Here we demonstrate the use of single-molecule atomic force microscopy to induce partial unfolding to expose metal centres in metalloproteins to aqueous solution, thus allowing for studying their chemical reactivity in aqueous solution for the first time. As a proof-of-principle, we demonstrate two chemical reactions for the FeS4 centre in rubredoxin: electrophilic protonation and nucleophilic ligand substitution. Our results show that protonation and ligand substitution result in mechanical destabilization of the FeS4 centre. Quantum chemical calculations corroborated experimental results and revealed detailed reaction mechanisms. We anticipate that this novel approach will provide insights into chemical reactivity of metal centres in metalloproteins under biologically more relevant conditions.

Dye-sensitized PS-b-P2VP-templated nickel oxide films for photoelectrochemical applications.

Moving from homogeneous water-splitting photocatalytic systems to photoelectrochemical devices requires the preparation and evaluation of novel p-type transparent conductive photoelectrode substrates. We report here on the sensitization of polystyrene-block-poly-(2-vinylpyridine) (PS-b-P2VP) diblock copolymer-templated NiO films with an organic push-pull dye. The potential of these new templated NiO film preparations for photoelectrochemical applications is compared with NiO material templated by F108 triblock copolymers. We conclude that NiO films are promising materials for the construction of dye-sensitized photocathodes to be inserted into photoelectrochemical (PEC) cells. However, a combined effort at the interface between materials science and molecular chemistry, ideally funded within a Global Artificial Photosynthesis Project, is still needed to improve the overall performance of the photoelectrodes and progress towards economically viable PEC devices.

Technical advances in molecular simulation since the 1980s.

This review describes how the theory and practice of molecular simulation have evolved since the beginning of the 1980s when the author started his career in this field. The account is of necessity brief and subjective and highlights the changes that the author considers have had significant impact on his research and mode of working.

Multiple proton relay routes in the reaction mechanism of RNAP II: assessing the effect of structural model.

RNA polymerase II catalyzes the nucleotidyl transfer reaction for messenger RNA synthesis in eukaryotes. Two crystal structures of this system have been resolved, each with its own defects in the coordination sphere of Mg(2+) (A) resulting from chemical modifications. We have used both structures and also a novel hybrid of the two that allows a better exploration of the parts of configuration space that reflect substrate-enzyme interactions. MD and QM/MM calculations have been performed, the latter with the semiempirical AM1/d-PhoT method, calibrated against density functional theory. Reaction path scans in 1-D provided insights about the role of Mg(2+) (A) which turns out to be more structural than catalytic. In contrast, 1-D scans of the incorporation of the nucleotidyl group yielded barriers that were much too high, necessitating the use of 2-D reaction coordinates. Three different proton acceptors for the initial reaction step were examined. For those models based on the two PDB structures the 2-D scans continued to yield very high barriers, indicating that the reaction is unlikely to proceed from these configurations. On the other hand, two hybrid models, chosen from the early and late parts of a 12 ns molecular dynamics simulation yielded greatly reduced barriers in the range of ∼ 17 to ∼ 27 kcal/mol for the three proton acceptors, as compared to the experimental estimate of 18 kcal/mol. The final step, release of pyrophosphate, was found to be facile. Our overall mechanism involves only active site residues or water without the need for external reactive agents such as the hydroxide ion previously proposed.

Theoretical modeling of low-energy electronic absorption bands in reduced cobaloximes.

The reduced Co(I) states of cobaloximes are powerful nucleophiles that play an important role in the hydrogen-evolving catalytic activity of these species. In this work we analyze the low-energy electronic absorption bands of two cobaloxime systems experimentally and use a variety of density functional theory and molecular orbital ab initio quantum chemical approaches. Overall we find a reasonable qualitative understanding of the electronic excitation spectra of these compounds but show that obtaining quantitative results remains a challenging task.