Base-catalyzed peptide hydrolysis is insensitive to mechanical stress

Oxadiazabicyclooctenone as a versatile monomer for the construction of pH sensitive functional polymers via ROMP

Herein, oxadiazabicyclooctenone is successfully developed as a versatile monomer for the construction of new pH-sensitive polymers by ROMP.

Structurally Linked Dynamics in Lactate Dehydrogenases of Evolutionarily Distinct Species

We present new findings about how primary and secondary structure affects the role of fast protein motions in the reaction coordinates of enzymatic reactions. Using transition path sampling and committor distribution analysis, we examined the difference in the role of these fast protein motions in the reaction coordinate of lactate dehydrogenases (LDHs) of Apicomplexa organisms Plasmodium falciparum and Cryptosporidium parvum. Having evolved separately from a common malate dehydrogenase ancestor, the two enzymes exhibit several important structural differences, notably a five-amino acid insertion in the active site loop of P. falciparum LDH. We find that these active site differences between the two organisms' LDHs likely cause a decrease in the contribution of the previously determined LDH rate-promoting vibration to the reaction coordinate of P. falciparum LDH compared to that of C. parvum LDH, specifically in the coupling of the rate-promoting vibration and the hydride transfer. This effect, while subtle, directly shows how changes in structure near the active site of LDH alter catalytically important motions. Insights provided by studying these alterations would prove to be useful in identifying LDH inhibitors that specifically target the isozymes of these parasitic organisms.

Ammonia as an In Situ Sanitizer: Influence of Virus Genome Type on Inactivation

Treatment of human excreta and animal manure (HEAM) is key in controlling the spread of persistent enteric pathogens, such as viruses. The extent of virus inactivation during HEAM storage and treatment appears to vary with virus genome type, although the reasons for this variability are not clear. Here, we investigated the inactivation of viruses of different genome types under conditions representative of HEAM storage or mesophilic digestion. The goals were to characterize the influence of HEAM solution conditions on inactivation and to determine the potential mechanisms involved. Specifically, eight viruses representing the four viral genome types (single-stranded RNA [ssRNA], double-stranded RNA [dsRNA], single-stranded DNA [ssDNA], and double-stranded DNA [dsDNA]) were exposed to synthetic solutions with well-controlled temperature (20 to 35°C), pH (8 to 9), and ammonia (NH3) concentrations (0 to 40 mmol liter(-1)). DNA and dsRNA viruses were considerably more resistant than ssRNA viruses, resulting in up to 1,000-fold-longer treatment times to reach a 4-log inactivation. The apparently slower inactivation of DNA viruses was rationalized by the higher stability of DNA than that of ssRNA in HEAM. Pushing the system toward harsher pH (>9) and temperature (>35°C) conditions, such as those encountered in thermophilic digestion and alkaline treatments, led to more consistent inactivation kinetics among ssRNA and other viruses. This suggests that the dependence of inactivation on genome type disappeared in favor of protein-mediated inactivation mechanisms common to all viruses. Finally, we recommend the use of MS2 as a conservative indicator to assess the inactivation of ssRNA viruses and the stable ΦX174 or dsDNA phages as indicators for persistent viruses.

IMPORTANCE

Viruses are among the most environmentally persistent pathogens. They can be present in high concentrations in human excreta and animal manure (HEAM). Therefore, appropriate treatment of HEAM is important prior to its reuse or discharge into the environment. Here, we investigated the factors that determine the persistence of viruses in HEAM, and we determined the main mechanisms that lead to their inactivation. Unlike other organisms, viruses can have four different genome types (double- or single-stranded RNA or DNA), and the viruses studied herein represent all four types. Genome type appeared to be the major determinant for persistence. Single-stranded RNA viruses are the most labile, because this genome type is susceptible to degradation in HEAM. In contrast, the other genome types are more stable; therefore, inactivation is slower and mainly driven by the degradation of viral proteins. Overall, this study allows us to better understand the behavior of viruses in HEAM.

Path Sampling Methods for Enzymatic Quantum Particle Transfer Reactions

The mechanisms of enzymatic reactions are studied via a host of computational techniques. While previous methods have been used successfully, many fail to incorporate the full dynamical properties of enzymatic systems. This can lead to misleading results in cases where enzyme motion plays a significant role in the reaction coordinate, which is especially relevant in particle transfer reactions where nuclear tunneling may occur. In this chapter, we outline previous methods, as well as discuss newly developed dynamical methods to interrogate mechanisms of enzymatic particle transfer reactions. These new methods allow for the calculation of free energy barriers and kinetic isotope effects (KIEs) with the incorporation of quantum effects through centroid molecular dynamics (CMD) and the full complement of enzyme dynamics through transition path sampling (TPS). Recent work, summarized in this chapter, applied the method for calculation of free energy barriers to reaction in lactate dehydrogenase (LDH) and yeast alcohol dehydrogenase (YADH). We found that tunneling plays an insignificant role in YADH but plays a more significant role in LDH, though not dominant over classical transfer. Additionally, we summarize the application of a TPS algorithm for the calculation of reaction rates in tandem with CMD to calculate the primary H/D KIE of YADH from first principles. We found that the computationally obtained KIE is within the margin of error of experimentally determined KIEs and corresponds to the KIE of particle transfer in the enzyme. These methods provide new ways to investigate enzyme mechanism with the inclusion of protein and quantum dynamics.

Enzymatic Kinetic Isotope Effects from First-Principles Path Sampling Calculations

In this study, we develop and test a method to determine the rate of particle transfer and kinetic isotope effects in enzymatic reactions, specifically yeast alcohol dehydrogenase (YADH), from first-principles. Transition path sampling (TPS) and normal mode centroid dynamics (CMD) are used to simulate these enzymatic reactions without knowledge of their reaction coordinates and with the inclusion of quantum effects, such as zero-point energy and tunneling, on the transferring particle. Though previous studies have used TPS to calculate reaction rate constants in various model and real systems, it has not been applied to a system as large as YADH. The calculated primary H/D kinetic isotope effect agrees with previously reported experimental results, within experimental error. The kinetic isotope effects calculated with this method correspond to the kinetic isotope effect of the transfer event itself. The results reported here show that the kinetic isotope effects calculated from first-principles, purely for barrier passage, can be used to predict experimental kinetic isotope effects in enzymatic systems.

Force-Induced Reversal of β-Eliminations: Stressed Disulfide Bonds in Alkaline Solution

Understanding the impact of tensile forces on disulfide bond cleavage is not only crucial to the breaking of cross-linkers in vulcanized materials such as strained rubber, but also to the regulation of protein activity by disulfide switches. By using ab initio simulations in the condensed phase, we investigated the response of disulfide cleavage by β-elimination to mechanical stress. We reveal that the rate-determining first step of the thermal reaction, which is the abstraction of the β-proton, is insensitive to external forces. However, forces larger than about 1 nN were found to reshape the free-energy landscape of the reaction so dramatically that a second channel is created, where the order of the reaction steps is reversed, turning β-deprotonation into a barrier-free follow-up process to C-S cleavage. This transforms a slow and force-independent process with second-order kinetics into a unimolecular reaction that is greatly accelerated by mechanical forces.

Transition path sampling with quantum/classical mechanics for reaction rates

Predicting rates of biochemical reactions through molecular simulations poses a particular challenge for two reasons. First, the process involves bond formation and/or cleavage and thus requires a quantum mechanical (QM) treatment of the reaction center, which can be combined with a more efficient molecular mechanical (MM) description for the remainder of the system, resulting in a QM/MM approach. Second, reaction time scales are typically many orders of magnitude larger than the (sub-)nanosecond scale accessible by QM/MM simulations. Transition path sampling (TPS) allows to efficiently sample the space of dynamic trajectories from the reactant to the product state without an additional biasing potential. We outline here the application of TPS and QM/MM to calculate rates for biochemical reactions, by means of a simple toy system. In a step-by-step protocol, we specifically refer to our implementation within the MD suite Gromacs, which we have made available to the research community, and include practical advice on the choice of parameters.

A density functional theory model of mechanically activated silyl ester hydrolysis

To elucidate the mechanism of the mechanically activated dissociation of chemical bonds between carboxymethylated amylose (CMA) and silane functionalized silicon dioxide, we have investigated the dissociation kinetics of the bonds connecting CMA to silicon oxide surfaces with density functional calculations including the effects of force, solvent polarizability, and pH. We have determined the activation energies, the pre-exponential factors, and the reaction rate constants of candidate reactions. The weakest bond was found to be the silyl ester bond between the silicon and the alkoxy oxygen atom. Under acidic conditions, spontaneous proton addition occurs close to the silyl ester such that neutral reactions become insignificant. Upon proton addition at the most favored position, the activation energy for bond hydrolysis becomes 31 kJ mol(-1), which agrees very well with experimental observation. Heterolytic bond scission in the protonated molecule has a much higher activation energy. The experimentally observed bi-exponential rupture kinetics can be explained by different side groups attached to the silicon atom of the silyl ester. The fact that different side groups lead to different dissociation kinetics provides an opportunity to deliberately modify and tune the kinetic parameters of mechanically activated bond dissociation of silyl esters.

Molecular Mechanochemistry: Engineering and Implications of Inherently Strained Architectures

Mechanical activation of chemical bonds is usually achieved by applying external forces. However, nearly all molecules exhibit inherent strain of their chemical bonds and angles as a result of constraints imposed by covalent bonding and interactions with the surrounding environment. Particularly strong deformation of bonds and angles is observed in hyperbranched macromolecules caused by steric repulsion of densely grafted polymer branches. In addition to the tension amplification, macromolecular architecture allows for accurate control of strain distribution, which enables focusing of the internal mechanical tension to specific chemical bonds and angles. As such, chemically identical bonds in self-strained macromolecules become physically distinct because the difference in bond tension leads to the corresponding difference in the electronic structure and chemical reactivity of individual bonds within the same macromolecule. In this review, we outline different approaches to the design of strained macromolecules along with physical principles of tension management, including generation, amplification, and focusing of mechanical tension at specific chemical bonds.

Model studies of force-dependent kinetics of multi-barrier reactions

According to transition state theory, the rate of a reaction that traverses multiple energy barriers is determined by the least stable (rate-determining) transition state. The preceding ('inner') energy barriers are kinetically 'invisible' but mechanistically significant. Here we show experimentally and computationally that the reduction rate of organic disulphides by phosphines in water, which in the absence of force proceeds by an equilibrium formation of a thiophosphonium intermediate, measured as a function of force applied on the disulphide moiety yields a usefully accurate estimate of the height of the inner barrier. We apply varying stretching force to the disulphide by incorporating it into a series of increasingly strained macrocycles. This force accelerates the reduction, even though the strain-free rate-determining step is orthogonal to the pulling direction. The observed rate-force correlation is consistent with the simplest model of force-dependent kinetics of a multi-barrier reaction.