Water dynamics and salt-activation of enzymes in organic media: mechanistic implications revealed by NMR spectroscopy

The role of conformational dynamics on the activity of polymer-conjugated CalB in organic solvents

A perennial interest in enzyme catalysis has been expanding its applicability from aqueous phase where enzymes are naturally evolved to organic solvents in which the majority of industrial chemical synthesis...

Water Networks and Correlated Motions in Mutant Isocitrate Dehydrogenase 1 (IDH1) Are Critical for Allosteric Inhibitor Binding and Activity

Point mutations in human isocitrate dehydrogenase 1 (IDH1) can drive malignancies, including lower-grade gliomas and secondary glioblastomas, chondrosarcomas, and acute myeloid leukemias. These mutations, which usually affect residue R132, ablate the normal activity of catalyzing the NADP⁺-dependent oxidation of isocitrate to α-ketoglutarate (αKG) while also acquiring a neomorphic activity of reducing αKG to d-2-hydroxyglutarate (D2HG). Mutant IDH1 can be selectively therapeutically targeted due to structural differences that occur in the wild type (WT) versus mutant form of the enzyme, though the full mechanisms of this selectivity are still under investigation. Here we probe the mechanistic features of the neomorphic activity and selective small molecule inhibition through a new lens, employing WaterMap and molecular dynamics simulations. These tools identified a high-energy path of water molecules connecting the inhibitor binding site with the αKG and NADP⁺ binding sites in mutant IDH1. This water path aligns spatially with the α10 helix from WT IDH1 crystal structures. Mutating residues at the termini of this water path specifically disrupted inhibitor binding and/or D2HG production, revealing additional key residues to consider in optimizing druglike molecules against mutant IDH1. Taken together, our findings from molecular simulations and mutant enzyme kinetic assays provide insight into how disrupting water paths through enzyme active sites can impact not only inhibitor potency but also substrate recognition and activity.

Protein Solvent Shell Structure Provides Rapid Analysis of Hydration Dynamics

The solvation layer surrounding a protein is clearly an intrinsic part of protein structure-dynamics-function, and our understanding of how the hydration dynamics influences protein function is emerging. We have recently reported simulations indicating a correlation between regional hydration dynamics and the structure of the solvation layer around different regions of the enzyme Candida antarctica lipase B, wherein the radial distribution function (RDF) was used to calculate the pairwise entropy, providing a link between dynamics (diffusion) and thermodynamics (excess entropy) known as Rosenfeld scaling. Regions with higher RDF values/peaks in the hydration layer (the first peak, within 6 Å of the protein surface) have faster diffusion in the hydration layer. The finding thus hinted at a handle for rapid evaluation of hydration dynamics at different regions on the protein surface in molecular dynamics simulations. Such an approach may move the analysis of hydration dynamics from a specialized venture to routine analysis, enabling an informatics approach to evaluate the role of hydration dynamics in biomolecular function. This paper first confirms that the correlation between regional diffusive dynamics and hydration layer structure (via water center of mass around protein side-chain atom RDF) is observed as a general relationship across a set of proteins. Second, it seeks to devise an approach for rapid analysis of hydration dynamics, determining the minimum amount of information and computational effort required to get a reliable value of hydration dynamics from structural data in MD simulations based on the protein-water RDF. A linear regression model using the integral of the hydration layer in the water-protein RDF was found to provide statistically equivalent apparent diffusion coefficients at the 95% confidence level for a set of 92 regions within five different proteins. In summary, RDF analysis of 10 ns of data after simulation convergence is sufficient to accurately map regions of fast and slow hydration dynamics around a protein surface. Additionally, it is anticipated that a quick look at protein-water RDFs, comparing peak heights, will be useful to provide a qualitative ranking of regions of faster and slower hydration dynamics at the protein surface for rapid analysis when investigating the role of solvent dynamics in protein function.

Water structure and dynamics in the hydration layer of a type III anti-freeze protein

We report on a molecular dynamics study on the relation between the structure and the orientational (and hydrogen bond) dynamics of hydration water around the ocean pout AFP III anti-freeze protein. We find evidence for an increasing tetrahedral structure from the area opposite to the ice binding site (IBS) towards the protein IBS, with the strongest signal of tetrahedral structure around the THR-18 residue of the IBS. The tetrahedral structural parameter mostly positively correlates with increased reorientation decay times. Interestingly, for several key (polar) residues that are not part of the IBS but are in its vicinity, we observe a decrease of the reorientation time with increasing tetrahedral structure. A similar anti-correlation is observed for the hydrogen-bonded water molecules. These effects are enhanced at a lower temperature. We interpret these results in terms of the structure-making and structure-breaking residues. Moreover, we investigate the tetrahedral structure and dynamics of waters at a partially dehydrated IBS, and for the protein adsorbed at the air-water interface. We find that the mutation changes the preferred protein orientation upon adsorption at an air-water interface. These results are in agreement with the water-air Vibration Sum Frequency Generation spectroscopic experiments showing a strongly reduced tetrahedral signal upon mutation at the IBS.

How Does Solvation Layer Mobility Affect Protein Structural Dynamics?

Solvation is critical for protein structural dynamics. Spectroscopic studies have indicated relationships between protein and solvent dynamics, and rates of gas binding to heme proteins in aqueous solution were previously observed to depend inversely on solution viscosity. In this work, the solvent-compatible enzyme Candida antarctica lipase B, which functions in aqueous and organic solvents, was modeled using molecular dynamics simulations. Data was obtained for the enzyme in acetonitrile, cyclohexane, n-butanol, and tert-butanol, in addition to water. Protein dynamics and solvation shell dynamics are characterized regionally: for each α-helix, β-sheet, and loop or connector region. Correlations are seen between solvent mobility and protein flexibility. So, does local viscosity explain the relationship between protein structural dynamics and solvation layer dynamics? Halle and Davidovic presented a cogent analysis of data describing the global hydrodynamics of a protein (tumbling in solution) that fits a model in which the protein's interfacial viscosity is higher than that of bulk water's, due to retarded water dynamics in the hydration layer (measured in NMR τ2 reorientation times). Numerous experiments have shown coupling between protein and solvation layer dynamics in site-specific measurements. Our data provides spatially-resolved characterization of solvent shell dynamics, showing correlations between regional solvation layer dynamics and protein dynamics in both aqueous and organic solvents. Correlations between protein flexibility and inverse solvent viscosity (1/η) are considered across several protein regions and for a rather disparate collection of solvents. It is seen that the correlation is consistently higher when local solvent shell dynamics are considered, rather than bulk viscosity. Protein flexibility is seen to correlate best with either the local interfacial viscosity or the ratio of the mobility of an organic solvent in a regional solvation layer relative to hydration dynamics around the same region. Results provide insight into the function of aqueous proteins, while also suggesting a framework for interpreting and predicting enzyme structural dynamics in non-aqueous solvents, based on the mobility of solvents within the solvation layer. We suggest that Kramers' theory may be used in future work to model protein conformational transitions in different solvents by incorporating local viscosity effects.

Large scale analysis of protein conformational transitions from aqueous to non-aqueous media

Background

Biocatalysis in organic solvents is nowadays a common practice with a large potential in Biotechnology. Several studies report that proteins which are co-crystallized or soaked in organic solvents preserve their fold integrity showing almost identical arrangements when compared to their aqueous forms. However, it is well established that the catalytic activity of proteins in organic solvents is much lower than in water. In order to explain this diminished activity and to further characterize the behaviour of proteins in non-aqueous environments, we performed a large-scale analysis (1737 proteins) of the conformational diversity of proteins crystallized in aqueous and co-crystallized or soaked in non-aqueous media.

Results

Using proteins’ experimentally determined conformational diversity taken from CoDNaS database, we found that proteins in non-aqueous media display much lower conformational diversity when compared to the corresponding conformers obtained in water. When conformational diversity is compared between conformers obtained in different non-aqueous media, their structural differences are larger and mostly independent of the presence of cognate ligands. We also found that conformers corresponding to non-aqueous media have larger but less flexible cavities, lower number of disordered regions and lower active-site residue mobility.

Conclusions

Our results show that non-aqueous media conformers have specific structural features and that they do not adopt extreme conformations found in aqueous media. This makes them clearly different from their corresponding aqueous conformers.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2044-2) contains supplementary material, which is available to authorized users.

Coupled Valence-Bond State Molecular Dynamics Description of an Enzyme-Catalyzed Reaction in a Non-Aqueous Organic Solvent

Enzymes are widely used in nonaqueous solvents to catalyze non-natural reactions. While experimental measurements showed that the solvent nature has a strong effect on the reaction kinetics, the molecular details of the catalytic mechanism in nonaqueous solvents have remained largely elusive. Here we study the transesterification reaction catalyzed by the paradigm subtilisin Carlsberg serine protease in an organic apolar solvent. The rate-limiting acylation step involves a proton transfer between active-site residues and the nucleophilic attack of the substrate to form a tetrahedral intermediate. We design the first coupled valence-bond state model that simultaneously describes both reactions in the enzymatic active site. We develop a new systematic procedure to parametrize this model on high-level ab initio QM/MM free energy calculations that account for the molecular details of the active site and for both substrate and protein conformational fluctuations. Our calculations show that the reaction energy barrier changes dramatically with the solvent and protein conformational fluctuations. We find that the mechanism of the tetrahedral intermediate formation during the acylation step is similar to that determined under aqueous conditions, and that the proton transfer and nucleophilic attack reactions occur concertedly. We identify the reaction coordinate to be mostly due to the rearrangement of some residual water molecules close to the active site.

Structural and functional stabilization of protein entities: state-of-the-art

Within the context of biomedicine and pharmaceutical sciences, the issue of (therapeutic) protein stabilization assumes particular relevance. Stabilization of protein and protein-like molecules translates into preservation of both structure and functionality during storage and/or targeting, and such stabilization is mostly attained through establishment of a thermodynamic equilibrium with the (micro)environment. The basic thermodynamic principles that govern protein structural transitions and the interactions of the protein molecule with its (micro)environment are, therefore, tackled in a systematic fashion. Highlights are given to the major classes of (bio)therapeutic molecules, viz. enzymes, recombinant proteins, (macro)peptides, (monoclonal) antibodies and bacteriophages. Modification of the microenvironment of the biomolecule via multipoint covalent attachment onto a solid surface followed by hydrophilic polymer co-immobilization, or physical containment within nanocarriers, are some of the (latest) strategies discussed aiming at full structural and functional stabilization of said biomolecules.

Mechanisms of acceleration and retardation of water dynamics by ions

There are fundamental and not yet fully resolved questions concerning the impact of solutes, ions in particular, on the structure and dynamics of water, which can be formulated as follows: Are the effects of ions local or long-ranged? Is the action of cations and anions on water cooperative or not? Here, we investigate how the reorientation and hydrogen-bond dynamics of water are affected by ions in dilute and concentrated aqueous salt solutions. By combining simulations and analytic modeling, we first show that ions have a short-ranged influence on the reorientation of individual water molecules and that depending on their interaction strength with water, they may accelerate or slow down water dynamics. A simple additive picture combining the effects of the cations and anions is found to provide a good description in dilute solutions. In concentrated solutions, we show that the average water reorientation time ceases to scale linearly with salt concentration due to overlapping hydration shells and structural rearrangements which reduce the translational displacements induced by hydrogen-bond switches and increase the solution viscosity. This effect is not ion-specific and explains why all concentrated salt solutions slow down water dynamics. Our picture, which is demonstrated to be robust vis-a-vis a change in the force-field, reconciles the seemingly contradictory experimental results obtained by ultrafast infrared and NMR spectroscopies, and suggests that there are no long-ranged cooperative ion effects on the dynamics of individual water molecules in dilute solutions.

Calcium-ion-induced stabilization of the protease from Bacillus cereus WQ9-2 in aqueous hydrophilic solvents: effect of calcium ion binding on the hydration shell and intramolecular interactions

The neutral protease WQ from Bacillus cereus is stable in various aqueous organic mixtures, with the exception of those containing acetonitrile (ACN) and dimethylformamide (DMF). The stability of the enzyme in aqueous hydrophilic solvents was dramatically enhanced with the addition of calcium ions, with the degree of improvement in the half-life relative to different solutions ranging from fourfold to more than 70-fold. Studies of the kinetic constants showed that calcium ions induced slight conformational changes in the active site of the enzyme in aqueous ACN. We investigated the molecular mechanisms underlying this stabilizing effect by employing a combination of biophysical techniques and molecular dynamics simulation. In aqueous ACN, the intrinsic fluorescence and circular dichroism analysis demonstrated that the addition of calcium ions induced a relatively compact conformation and maintained both the native-like microenvironment near the tryptophan residues and the secondary structure. Alternatively, homology modeling confirmed the location of four calcium-ion-binding sites in the enzyme, and molecular dynamics simulation revealed that three other calcium ions were bound to the surface of the enzyme. Calcium ions, known as a type of kosmotrope, can strongly bond with water molecules, thus aiding in the formation of the regional hydration shell required for the maintenance of enzyme activity. In addition, the introduction of calcium ions resulted in the formation of additional ionic interactions, providing propitious means for protein stabilization. Thus, the stronger intramolecular interactions were also expected to contribute partially to the enhanced stability of the enzyme in an aqueous organic solvent.

Recent trends in biocatalysis engineering

During the last 30 years the scope of biocatalysis has been expanding due to the advances in several technological fields. Diverse techniques as structural enzyme improvement (e.g. protein engineering, direct evolution), engineering approaches (e.g. ionic liquids, supercritical fluids) and physical stabilization (e.g. immobilization, CLEAS) have been developed, which in combination are powerful tools to improve biotransformation and to synthesize new products. In the present work, recent advances in biocatalysis are reviewed.

The dynamic determinants of reaction specificity in the IMPDH/GMPR family of (β/α)(8) barrel enzymes

The inosine monophosphate dehydrogenase (IMPDH)/guanosine monophosphate reductase (GMPR) family of (β/α)(8) enzymes presents an excellent opportunity to investigate how subtle changes in enzyme structure change reaction specificity. IMPDH and GMPR bind the same ligands with similar affinities and share a common set of catalytic residues. Both enzymes catalyze a hydride transfer reaction involving a nicotinamide cofactor hydride, and both reactions proceed via the same covalent intermediate. In the case of IMPDH, this intermediate reacts with water, while in GMPR it reacts with ammonia. In both cases, the two chemical transformations are separated by a conformational change. In IMPDH, the conformational change involves a mobile protein flap while in GMPR, the cofactor moves. Thus reaction specificity is controlled by differences in dynamics, which in turn are controlled by residues outside the active site. These findings have some intriguing implications for the evolution of the IMPDH/GMPR family.

Allosteric activation via kinetic control: potassium accelerates a conformational change in IMP dehydrogenase

Allosteric activators are generally believed to shift the equilibrium distribution of enzyme conformations to favor a catalytically productive structure; the kinetics of conformational exchange is seldom addressed. Several observations suggested that the usual allosteric mechanism might not apply to the activation of IMP dehydrogenase (IMPDH) by monovalent cations. Therefore, we investigated the mechanism of K(+) activation in IMPDH by delineating the kinetic mechanism in the absence of monovalent cations. Surprisingly, the K(+) dependence of k(cat) derives from the rate of flap closure, which increases by ≥65-fold in the presence of K(+). We performed both alchemical free energy simulations and potential of mean force calculations using the orthogonal space random walk strategy to computationally analyze how K(+) accelerates this conformational change. The simulations recapitulate the preference of IMPDH for K(+), validating the computational models. When K(+) is replaced with a dummy ion, the residues of the K(+) binding site relax into ordered secondary structure, creating a barrier to conformational exchange. K(+) mobilizes these residues by providing alternate interactions for the main chain carbonyls. Potential of mean force calculations indicate that K(+) changes the shape of the energy well, shrinking the reaction coordinate by shifting the closed conformation toward the open state. This work suggests that allosteric regulation can be under kinetic as well as thermodynamic control.

"Fluctuograms" reveal the intermittent intra-protein communication in subtilisin Carlsberg and correlate mechanical coupling with co-evolution

The mechanism of intra-protein communication and allosteric coupling is key to understanding the structure-property relationship of protein function. For subtilisin Carlsberg, the Ca²⁺-binding loop is distal to substrate-binding and active sites, yet the serine protease function depends on Ca²⁺ binding. The atomic molecular dynamics (MD) simulations of apo and Ca²⁺-bound subtilisin show similar structures and there is no direct evidence that subtilisin has alternative conformations. To model the intra-protein communication due to Ca²⁺ binding, we transform the sequential segments of an atomic MD trajectory into separate elastic network models to represent anharmonicity and nonlinearity effectively as the temporal and spatial variation of the mechanical coupling network. In analogy to the spectrogram of sound waves, this transformation is termed the “fluctuogram” of protein dynamics. We illustrate that the Ca²⁺-bound and apo states of subtilisin have different fluctuograms and that intra-protein communication proceeds intermittently both in space and in time. We found that residues with large mechanical coupling variation due to Ca²⁺ binding correlate with the reported mutation sites selected by directed evolution for improving the stability of subtilisin and its activity in a non-aqueous environment. Furthermore, we utilize the fluctuograms calculated from MD to capture the highly correlated residues in a multiple sequence alignment. We show that in addition to the magnitude, the variance of coupling strength is also an indicative property for the sequence correlation observed in a statistical coupling analysis. The results of this work illustrate that the mechanical coupling networks calculated from atomic details can be used to correlate with functionally important mutation sites and co-evolution.

A hallmark of protein molecules is their machine-like behaviors while carrying out biological functions. At the molecular level, molecular signals such as binding a metal ion at an action site can cause long-range effects and alter protein function. Such phenomena are often referred to as intra-protein communication or allosteric coupling. Elucidating the underlying mechanisms could lead to novel discovery of molecular modulators to regulate protein function in a more specific and effective manner. A long-standing puzzle is the roles of the anharmonicity and nonlinearity in protein dynamics. To incorporate these characters in modeling intra-protein communication, we devise a “fluctuogram” analysis to record the choreography of allosteric coupling in an atomic molecular dynamics simulation. We show that fluctuogram analysis can bridge the results of physics-based simulation and sequence alignment in bioinformatics by capturing the residues that exhibit high correlation in a multiple sequence alignment. We also show that the fluctuograms calculated from atomic details have the potential to be applied as a tool to select mutation sites for modulating protein function.

On the roles of substrate binding and hinge unfolding in conformational changes of adenylate kinase

We characterized the conformational change of adenylate kinase (AK) between open and closed forms by conducting five all-atom molecular-dynamics simulations, each of 100 ns duration. Different initial structures and substrate binding configurations were used to probe the pathways of AK conformational change in explicit solvent, and no bias potential was applied. A complete closed-to-open and a partial open-to-closed transition were observed, demonstrating the direct impact of substrate-mediated interactions on shifting protein conformation. The sampled configurations suggest two possible pathways for connecting the open and closed structures of AK, affirming the prediction made based on available x-ray structures and earlier works of coarse-grained modeling. The trajectories of the all-atom molecular-dynamics simulations revealed the complexity of protein dynamics and the coupling between different domains during conformational change. Calculations of solvent density and density fluctuations surrounding AK did not show prominent variation during the transition between closed and open forms. Finally, we characterized the effects of local unfolding of an important hinge near Pro(177) on the closed-to-open transition of AK and identified a novel mechanism by which hinge unfolding modulates protein conformational change. The local unfolding of Pro(177) hinge induces alternative tertiary contacts that stabilize the closed structure and prevent the opening transition.

Enantioselective Transesterification Catalysis by Nanosized Serine Protease Subtilisin Carlsberg Particles in Tetrahydrofuran

Enzyme catalysis in organic solvents is a powerful tool for stereo-selective synthesis but the enantioselectivity is still hard to predict. To overcome this obstacle, we employed a nanoparticulate formulation of subtilisin Carlsberg (SC) and designed a series of 14 structurally related racemic alcohols. They were employed in the model transesterification reaction with vinyl butyrate and the enantioselectivities were determined. In general, short alcohol side chains led to low enantioselectivties, while larger and bulky side chains caused better discrimination of the enantiomers by the enzyme. With several bulky substrates high enantioselectivities with E>100 were obtained. Computational modeling highlighted that key to high enantioselectivity is the discrimination of the R and S substrates by the sole hydrophobic binding pocket based on their size and bulkiness. While bulky S enantiomer side chains could be accommodated within the binding pocket, bulky R enantiomer side chains could not. However, when also the S enantiomer side chain becomes too large and does not fit into the binding pocket anymore, enantioselectivity accordingly drops.

Biocatalyst activity in nonaqueous environments correlates with centisecond-range protein motions

Recent studies exploring the relationship between enzymatic catalysis and protein dynamics in the aqueous phase have yielded evidence that dynamics and enzyme activity are strongly correlated. Given that protein dynamics are significantly attenuated in organic solvents and that proteins exhibit a wide range of motions depending on the specific solvent environment, the nonaqueous milieu provides a unique opportunity to examine the role of protein dynamics in enzyme activity. Variable-temperature kinetic measurements, X-band electron spin resonance spectroscopy, (1)H NMR relaxation, and (19)F NMR spectroscopy experiments were performed on subtilisin Carlsberg colyophilized with several inorganic salts and suspended in organic solvents. The results indicate that salt activation induces a greater degree of transition-state flexibility, reflected by a more positive DeltaDeltaS(dagger), for the more active biocatalyst preparations in organic solvents. In contrast, DeltaDeltaH(dagger) was negligible regardless of salt type or salt content. Electron spin resonance spectroscopy and (1)H NMR relaxation measurements, including spin-lattice relaxation, spin-lattice relaxation in the rotating frame, and longitudinal magnetization exchange, revealed that the enzyme's turnover number (k(cat)) was strongly correlated with protein motions in the centisecond time regime, weakly correlated with protein motions in the millisecond regime, and uncorrelated with protein motions on the piconanosecond timescale. In addition, (19)F chemical shift measurements and hyperfine tensor measurements of biocatalyst formulations inhibited with 4-fluorobenzenesulfonyl fluoride and 4-ethoxyfluorophosphinyl-oxy-TEMPO, respectively, suggest that enzyme activation was only weakly affected by changes in active-site polarity.

Activity of Subtilisin Carlsberg in macromolecular crowding

Enzymatic activity of a proteolytic enzyme Subtilisin Carlsberg (SC) in anionic sodium dodecyl sulfate (SDS) micellar medium has been explored and found to be retarded compared to that in bulk buffer. Circular dichroism (CD) study reveals that SDS, which is a potential protein denaturant, has an insignificant denaturation effect on SC. The structural integrity of the protein offers an opportunity to study the functionality of the enzyme SC in a macromolecular crowding of micelles. Dynamic light scattering (DLS) data indicates no sandwich-like micelle-SC complex formation ruling out the possibility of interaction of the enzyme with the hydrophobic core of the micelle. However, steady state and time resolved emission studies on specific and nonspecific fluorescent probes indicate the proximity effect at the surface of the enzyme due to macromolecular crowding of the micelles. The agreement of retarded enzymatic activity in the micellar crowd with a theoretical model ascribed to the facts that substrates are compartmentalized in the micelles and enzyme interacts with the micelle through stern layer.

Stability of biocatalysts

Despite their many favorable qualities, the marginal stability of biocatalysts in many types of reaction media often has prevented or delayed their implementation for industrial-scale synthesis of fine chemicals and pharmaceuticals. Consequently, there is great interest in understanding effects of solution conditions on protein stability, as well as in developing strategies to improve protein stability in desired reaction media. Recent methods include novel chemical modifications of protein, lyophilization in the presence of additives, and physical immobilization on novel supports. Rational and combinatorial protein engineering techniques have been used to yield unmodified proteins with exceptionally improved stability. Both have been aided by the development of computational tools and structure-guided heuristics aimed at reducing library sizes that must be generated and screened to identify improved mutants. The number of parameters used to indicate protein stability can complicate discussions and investigations, and care should be taken to identify whether thermodynamic or kinetic stability limits the observed stability of proteins. Although the useful lifetime of a biocatalyst is dictated by its kinetic stability, only 6% of protein stability studies use kinetic stability measures. Clearly, more effort is needed to study how solution conditions impact protein kinetic stability.

Mechanism for intein C-terminal cleavage: a proposal from quantum mechanical calculations

Inteins are autocatalytic protein cleavage and splicing elements. A cysteine to alanine mutation at the N-terminal of inteins inhibits splicing and isolates the C-terminal cleavage reaction. Experiments indicate an enhanced C-terminal cleavage reaction rate upon decreasing the solution pH for the cleavage mutant, which cannot be explained by the existing mechanistic framework. We use intein crystal structure data and the information about conserved amino acids to perform semiempirical PM3 calculations followed by high-level density functional theory calculations in both gas phase and implicit solvent environments. Based on these calculations, we propose a detailed "low pH" mechanism for intein C-terminal cleavage. Water plays an important role in the proposed reaction mechanism, acting as an acid as well as a base. The protonation of the scissile peptide bond nitrogen by a hydronium ion is an important first step in the reaction. That step is followed by the attack of the C-terminal asparagine side chain on its carbonyl carbon, causing succinimide formation and simultaneous peptide bond cleavage. The computed reaction energy barrier in the gas phase is approximately 33 kcal/mol and reduces to approximately 25 kcal/mol in solution, close to the 21 kcal/mol experimentally observed at pH 6.0. This mechanism is consistent with the observed increase in C-terminal cleavage activity at low pH for the cleavage mutant of the Mycobacterium tuberculosis RecA mini-intein.