Molecular Dynamics Studies of Proton Transport in Hydrogenase and Hydrogenase Mimics.

There is extensive interest in hydrogenases based on their ability to rapidly and efficiently interconvert H2 with protons and electrons, and their (typically) superior function relative to molecular mimics. Understanding the function of enzymes is one approach to implementing design features to make better catalysts and is an approach we have implemented in our work. In this review, we will discuss our efforts to develop design principles from enzymes, with specific focus on proton transport. We will also present computational studies of the mimics we have investigated with similar methodologies. We will discuss the mechanisms used by small scaffolds on molecular mimics which in many cases are surprisingly similar to those used by nature, while in other cases, computational analysis allowed us to reveal an unexpected role. Computational methods provide one of the best ways, and in some cases, the only way, to gain insight into the mechanistic details of enzymes. In this review, we illustrate the general computational strategy we used to study the proton pathway of [FeFe]-hydrogenase, and the similar strategy to investigate small molecules. We present the main results we obtained and how our computational work stimulated or worked in concert with experimental investigations. We also focus on estimation of errors and convergence of properties in the simulations. These studies demonstrate the powerful results that can be obtained by the close pairing of experimental and theoretical approaches.

Mechanistic insights into hydride transfer for catalytic hydrogenation of CO(2) with cobalt complexes.

The catalytic hydrogenation of CO2 to formate by Co(dmpe)2H can proceed via direct hydride transfer or via CO2 coordination to Co followed by reductive elimination of formate. The different nature of the rate-determining step in the two mechanisms may provide new insights into designing catalysts with improved performance.

Structural aspects of the solvation shell of lysine and acetylated lysine: A Car-Parrinello and classical molecular dynamics investigation.

Lysine acetylation is a post-translational modification, which modulates the affinity of protein-protein and/or protein-DNA complexes. Its crucial role as a switch in signaling pathways highlights the relevance of charged chemical groups in determining the interactions between water and biomolecules. A great effort has been recently devoted to assess the reliability of classical molecular dynamics simulations in describing the solvation properties of charged moieties. In the spirit of these investigations, we performed classical and Car-Parrinello molecular dynamics simulations on lysine and acetylated-lysine in aqueous solution. A comparative analysis between the two computational schemes is presented with a focus on the first solvation shell of the charged groups. An accurate structural analysis unveils subtle, yet statistically significant, differences which are discussed in connection to the significant electronic density charge transfer occurring between the solute and the surrounding water molecules.

Large-scale motions and electrostatic properties of furin and HIV-1 protease.

We present a comparative study between two members of serine and aspartic proteases complexed with a peptide substrate. The same computational setup is used to characterize the structural, electrostatic, and electronic properties for the Michaelis complex of furin, a serine protease, and of the aspartic protease from HIV-1. In both cases plane-wave density functional theory (PW-DFT) and empirical force-field-based molecular dynamics calculations are used. For furin, calculations are extended to the complex with the intermediate of the first step of the reaction. Comparisons are also made with results from recent PW-DFT investigations on both families of enzymes and with the same chemical groups in an aqueous environment. It is found that the substrate carbonyl group is more polarized in the furin complex than in the HIV-1 protease one. A further difference regards the large-scale motions of the complexes as a whole and local conformational fluctuations at the active site. The global and local fluctuations are well coupled for HIV-1 protease but not for furin. Thus, despite some chemical analogies in the first step of the reaction mechanism, furin and HIV-1 protease complexes appear to be characterized by a different interplay of electrostatics and conformational fluctuations.