Mechanical aspects of nitrile hydratase enzymatic activity. Steered molecular dynamics simulations of Pseudonocardia thermophila JCM 3095

Structure-oriented engineering of nitrile hydratase: Reshaping of substrate access tunnel and binding pocket for efficient synthesis of cinnamamide

Cinnamamide and its derivatives are the most common and important building blocks widely present in natural products. Currently, nitrile hydratase (NHase, EC 4.2.1.84) has been widely used in large-scale industrial production of nicotinamide and acrylamide, while its catalytic activity is extremely low or inactive for bulky nitrile substrates such as cinnamonitrile. Therefore, beneficial variant βF37P/L48P/F51N were obtained from PtNHase of Pseudonocardia thermophila JCM3095 by reshaping of substrate access tunnel and binding pocket, which exhibited 14.88-fold improved catalytic efficiency compared to the wild-type PtNHase. Structure analysis, molecular dynamics simulations and dynamical cross-correlation matrix (DCCM) analysis revealed that the introduced mutations enlarged the substrate access tunnel and binding pocket, enhanced overall anti-correlated movements of enzymes, which would promote product release during the dynamic process of catalysis. In a hydration process, the complete conversion of 5 mM cinnamonitrile was achieved by βF37P/L48P/F51N in a 50 mL reaction, with cinnamamide yield of almost 100 % and productivity of 0.736 g L-1 h-1. The study demonstrates the co-evolution of substrate access tunnel and binding pocket is an effective strategy, and provides a valuable reference for future research. Furthermore, NHases have huge potential for catalyzing bulky nitriles to form corresponding amides in large-scale industrial production.

N-terminal loops at the tetramer interface of nitrile hydratase act as "hooks" determining resistance to high amide concentrations

Nitrile hydratase (NHase) has been extensively utilized in industrial acrylamide production. However, the vulnerability to high concentrations of acrylamide limits its further application. Herein, we redesigned the N-terminal loop at the tetramer interface of a thermophilic NHase from Pseudonocardia thermophila JCM3095 (PtNHase), and its catalytic activity, resistance to high acrylamide concentrations, and thermostability were improved. Amino acid residues located in the N-terminal loop of the tetramer interface that are responsible for enhancing the resistance to high acrylamide concentrations were identified via static structural analysis and molecular dynamics simulations. A variant library was used to fine-tune the tetramer interface. Variant αL6T exhibited 3.5-fold greater resistance to 50% (v/v) acrylamide, whereas its activity was 1.2-fold higher than that of the wild-type (WT) enzyme, revealing no activity-stability trade-off. Compared to the use of Escherichia coli harboring the WT enzyme, the use of E. coli harboring αL6T increased the acrylamide concentration from 398.1 g/L to 500 g/L. Crystal structure-guided analysis of αL6T and molecular dynamics simulations revealed that increased enzyme surface hydration and the introduction of positive cross-correlation into the N-terminal loop of the tetramer interface caused the two loop regions to hook to each other, thus improving the resistance to high acrylamide concentrations.

"Toolbox" construction of an extremophilic nitrile hydratase from Streptomyces thermoautotrophicus for the promising industrial production of various amides

Nitrile hydratase (NHase; EC 4.2.1.84) is widely used to synthesize the corresponding amides from nitriles, which is the most successful green biocatalyst. However, the limited acceptability of substrates and instability under harsh reaction conditions have hindered its widespread industrial application. Here, a gene encoding an extremophilic NHase from Streptomyces thermoautotrophicus (S.t NHase) was successfully overexpressed in Escherichia coli. The enzyme exhibited excellent thermostability, retaining >50 % of residual activity after heat treatment at 65 °C for 252 min. To further improve the catalytic performance of S.t NHase, semi-rational engineering of its substrate access tunnel was performed. A mutant βL48D showed a specific activity of 566.18 ± 18.86 U/mg towards 3-cyanopyridine, which was 7.7 times higher than its parent enzyme (73.80 ± 5.76 U/mg). Molecular dynamics simulation showed that the introduction of aspartic acid into βLeu48 resulted in a larger and more frequent opening of the substrate access tunnel entrance. On this basis, a "toolbox" containing various mutants on the substrate access tunnel was further established, whose catalytic activity towards various nitrile substrates was extensively improved, showing great potential for efficient synthesis of multiple high-value amides.

High-Level Expression of Nitrile Hydratase in Escherichia coli for 2-Amino-2,3-Dimethylbutyramide Synthesis

In the synthesis of imidazolinone herbicides, 2-Amino-2,3-dimethylbutyramide (ADBA) is an important intermedium. In this study, the recombinant production of nitrile hydratase (NHase) in Escherichia coli for ADBA synthesis was explored. A local library containing recombinant NHases from various sources was screened using a colorimetric method. NHase from Pseudonocardia thermophila JCM3095 was selected, fused with a His-tag and one-step purified. The enzymatic properties of recombinant NHase were studied and indicated robust thermal stability and inhibition of cyanide ions due to substrate degradation. After systematic optimization of fermentation conditions, the OD600 (optical density at 600 nm), enzyme activity and specific activity of recombinant strain E. coli BL21(DE3)/pET-28a+NHase reached 19.4, 3.72 U/mL and 1.04 U/mg protein at 42 h, representing 5.86-, 26.6- and 4-fold increases, respectively. These results offered an efficient recombinant whole-cell biocatalyst for ADBA synthesis.

Insight into the broadened substrate scope of nitrile hydratase by static and dynamic structure analysis

Mutations of two gating residues at the substrate access tunnel entrance direct the substrate scope of NHases.

Substrate access tunnel engineering for improving the catalytic activity of a thermophilic nitrile hydratase toward pyridine and pyrazine nitriles

Nitrile hydratase (NHase) is able to bio-transform nitriles into amides. As nitrile hydration being an exothermic reaction, a NHase with high activity and stability is needed for amide production. However, the widespread use of NHase for amide bio-production is limited by an activity-stability trade-off. In this study, through the combination of substrate access tunnel calculation, residue conservative analysis and site-saturation mutagenesis, a residue located at the substrate access tunnel entrance of the thermophilic NHase from extremophile Caldalkalibacillus thermarum TA2. A1, βLeu48, was semi-rationally identified as a potential gating residue that directs the enzymatic activity toward various pyridine and pyrazine nitriles. The specific activity of the corresponding mutant βL48H towards 3-cyanopyridine, 2-cyanopyridine and cyanopyrazine were 2.4-fold, 2.8-fold and 3.1-fold higher than that of its parent enzyme, showing a great potential in the industrial production of high-value pyridine and pyrazine carboxamides. Further structural analysis demonstrated that the βHis48 could form a long-lasting hydrogen bond with αGlu166, which contributes to the expansion of the entrance of substrate access tunnel and accelerate substrate migration.

Effect and mechanism analysis of different linkers on efficient catalysis of subunit-fused nitrile hydratase

Protein fusion using a linker plays an important role for protein evolution. However, designing suitable linkers for protein evolution is yet challenging and under-explored. To further clarify the regular pattern of suitable type of linker for fusion proteins, one nitrile hydratase (NHase) was used as a target protein and subunit fusion strategy was carried out to improve its efficient catalysis. Subunit-fused variants with three different types of linkers were constructed and characterized. All variants exhibited higher stability than that of the wild type. The longer the linker was, the higher stability NHase showed, however, too long linker affected NHase activity and expression. Among the three types of linkers, the α-helical linker seemed more suitable for NHase than flexible or rigid linkers. Though it is not clear how the linkers affecting the activity, structure analysis indicated that the stability improvement is dependent on the additional salt bridge, H-bond, and the subunit interface area increasing due to the linker insertion, among which the additional salt bridge and interface area were more important factors. The results described here may be useful for redesigning other enzymes through subunit fusion.

Pseudonocardia broussonetiae sp. nov., an endophytic actinomycete isolated from the roots of Broussonetia papyrifera

A novel endophytic actinomycete, designated strain Gen 01^T, was isolated from the roots of Broussonetia papyrifera and characterized by using a polyphasic approach. The predominant cellular fatty acids were iso-C_16 : 0, summed feature 3, iso H-C_16 : 1, C_16 : 0 and iso-C_14 : 0. The polar lipid profile consisted of diphosphatidylglycerol, phosphatidylethanolamine, phosphatidylinositol mannosides, phospholipids of unknown structure containing glucosamine inositol, phosphatidylinositol and unidentified phospholipids. The major menaquinone was MK-8 (H₄). The DNA G+C content of the genome sequence, consisting of 7 177 725 bp, was 74.5 mol%. Phylogenetic analysis of the full-length 16S rRNA gene sequences showed that strain Gen 01^T belongs to the genus Pseudonocardia with the highest sequence similarity to Pseudonocardia petroleophila CGMCC 4.1532^T (98.9 %) and lower than 98.7 % similarity to other species of the genus Pseudonocardia with validly published names. The average nucleotide identity and digital DNA-DNAhybridization values between strain Gen 01^T and P. petroleophila CGMCC 4.1532^T were 84.6 and 30.9 %, respectively. Furthermore, the morphological, physiological and biochemical characteristics were sufficient to categorize strain Gen 01^T as being distinct from P. petroleophila CGMCC 4.1532^T. Consequently, based on phenotypic and genotypic characteristics, strain Gen 01^T represents a novel species of the genus Pseudonocardia, for which the name Pseudonocardia broussonetiae sp. nov. is proposed. The type strain is Gen 01^T (=CICC 24820^T=JCM 33840^T).

Computational Design of Nitrile Hydratase from Pseudonocardia thermophila JCM3095 for Improved Thermostability

High thermostability and catalytic activity are key properties for nitrile hydratase (NHase, EC 4.2.1.84) as a well-industrialized catalyst. In this study, rational design was applied to tailor the thermostability of NHase from Pseudonocardia thermophila JCM3095 (PtNHase) by combining FireProt server prediction and molecular dynamics (MD) simulation. Site-directed mutagenesis of non-catalytic residues provided by the rational design was subsequentially performed. The positive multiple-point mutant, namely, M10 (αI5P/αT18Y/αQ31L/αD92H/βA20P/βP38L/βF118W/βS130Y/βC189N/βC218V), was obtained and further analyzed. The Melting temperature (T_m) of the M10 mutant showed an increase by 3.2 °C and a substantial increase in residual activity of the enzyme at elevated temperatures was also observed. Moreover, the M10 mutant also showed a 2.1-fold increase in catalytic activity compared with the wild-type PtNHase. Molecular docking and MD simulations demonstrated better substrate affinity and improved thermostability for the mutant.

Recent Advances and Promises in Nitrile Hydratase: From Mechanism to Industrial Applications

Nitrile hydratase (NHase, EC 4.2.1.84) is one type of metalloenzyme participating in the biotransformation of nitriles into amides. Given its catalytic specificity in amide production and eco-friendliness, NHase has overwhelmed its chemical counterpart during the past few decades. However, unclear catalytic mechanism, low thermostablity, and narrow substrate specificity limit the further application of NHase. During the past few years, numerous studies on the theoretical and industrial aspects of NHase have advanced the development of this green catalyst. This review critically focuses on NHase research from recent years, including the natural distribution, gene types, posttranslational modifications, expression, proposed catalytic mechanism, biochemical properties, and potential applications of NHase. The developments of NHase described here are not only useful for further application of NHase, but also beneficial for the development of the fields of biocatalysis and biotransformation.

Substrate Profiling of the Cobalt Nitrile Hydratase from Rhodococcus rhodochrous ATCC BAA 870

The aromatic substrate profile of the cobalt nitrile hydratase from Rhodococcus rhodochrous ATCC BAA 870 was evaluated against a wide range of nitrile containing compounds (>60). To determine the substrate limits of this enzyme, compounds ranging in size from small (90 Da) to large (325 Da) were evaluated. Larger compounds included those with a bi-aryl axis, prepared by the Suzuki coupling reaction, Morita-Baylis-Hillman adducts, heteroatom-linked diarylpyridines prepared by Buchwald-Hartwig cross-coupling reactions and imidazo[1,2-a]pyridines prepared by the Groebke-Blackburn-Bienaymé multicomponent reaction. The enzyme active site was moderately accommodating, accepting almost all of the small aromatic nitriles, the diarylpyridines and most of the bi-aryl compounds and Morita-Baylis-Hillman products but not the Groebke-Blackburn-Bienaymé products. Nitrile conversion was influenced by steric hindrance around the cyano group, the presence of electron donating groups (e.g., methoxy) on the aromatic ring, and the overall size of the compound.

Role of Water Hydrogen Bonding on Transport of Small Molecules inside Hydrophobic Channels

We conducted a systematic analysis of water networking inside smooth hyperboloid hydrophobic structures (cylindrical, barrel, and hourglass shapes) to elucidate the role of water hydrogen bonding on the transport of small hydrophobic molecules (ligands). Through a series of molecular dynamics simulations, we established that a hydrogen-bonded network forming along the centerline resulted in a water exclusion zone adjacent to the walls. The size of the exclusion zone is a function of the geometry and the nonbonded interaction strength, defining the effective hydrophobicity of the structure. Exclusion of water molecules from this zone results in lower apparent viscosity, leading to acceleration of ligand transport up to 7 times faster than that measured in the bulk. Transport of ligands into and out of the hydrophobic structures was shown to be controlled by a single water molecule that capped the narrow regions in the structure. This mechanism provides physical insights into the behavior and role of water in the bottleneck regions of hydrophobic enzyme channels. These findings were then used in a sister publication [ Escalante , D. E. , Comput. Struct. Biotechnol. J. 2019 17 757 760 ] to develop a model that can accurately predict the transport of ligands along nanochannels of broad-substrate specificity enzymes.

Simulated and experimental force spectroscopy of lysozyme on silica

The force spectra of proteins detaching from oxide surfaces measured by atomic force microscopy (AFM) often present complex patterns of peaks, which are difficult to correlate with individual bond-breaking events at the atomic scale. In this work we rationalize experimental AFM force spectra of hen-egg-white lysozyme detaching from silica by means of all-atom steered molecular dynamics (SMD) simulations. In particular, we demonstrate that the native tertiary structure of lysozyme is preserved if, and only if, its four intramolecular disulfide bridges are intact. Otherwise, the protein pulled off the surface undergoes severe unfolding, which is well captured by SMD simulations in explicit solvent. Implicit solvent simulations, on the contrary, wrongly predict protein unfolding even in the presence of S-S bridges, due to the lack of additional structural stabilization provided by the water's hydrogen-bond network within and surrounding the protein. On the basis of our combined experimental and theoretical findings, we infer that the rugged force spectra characteristic of lysozyme/silica interfaces are not due to the successive breaking of internal disulfide bonds leading to partial unfolding events. Rather, they reflect the detachment of several molecules bound to the same AFM tip, each anchored to the surface via multiple hydrogen and ionic bonds.

Ligand diffusion in proteins via enhanced sampling in molecular dynamics

Computational simulations in biophysics describe the dynamics and functions of biological macromolecules at the atomic level. Among motions particularly important for life are the transport processes in heterogeneous media. The process of ligand diffusion inside proteins is an example of a complex rare event that can be modeled using molecular dynamics simulations. The study of physical interactions between a ligand and its biological target is of paramount importance for the design of novel drugs and enzymes. Unfortunately, the process of ligand diffusion is difficult to study experimentally. The need for identifying the ligand egress pathways and understanding how ligands migrate through protein tunnels has spurred the development of several methodological approaches to this problem. The complex topology of protein channels and the transient nature of the ligand passage pose difficulties in the modeling of the ligand entry/escape pathways by canonical molecular dynamics simulations. In this review, we report a methodology involving a reconstruction of the ligand diffusion reaction coordinates and the free-energy profiles along these reaction coordinates using enhanced sampling of conformational space. We illustrate the above methods on several ligand-protein systems, including cytochromes and G-protein-coupled receptors. The methods are general and may be adopted to other transport processes in living matter.

Memetic algorithms for ligand expulsion from protein cavities

Ligand diffusion through a protein interior is a fundamental process governing biological signaling and enzymatic catalysis. A complex topology of channels in proteins leads often to difficulties in modeling ligand escape pathways by classical molecular dynamics simulations. In this paper, two novel memetic methods for searching the exit paths and cavity space exploration are proposed: Memory Enhanced Random Acceleration (MERA) Molecular Dynamics (MD) and Immune Algorithm (IA). In MERA, a pheromone concept is introduced to optimize an expulsion force. In IA, hybrid learning protocols are exploited to predict ligand exit paths. They are tested on three protein channels with increasing complexity: M2 muscarinic G-protein-coupled receptor, enzyme nitrile hydratase, and heme-protein cytochrome P450cam. In these cases, the memetic methods outperform simulated annealing and random acceleration molecular dynamics. The proposed algorithms are general and appropriate in all problems where an accelerated transport of an object through a network of channels is studied.

A switch in a substrate tunnel for directing regioselectivity of nitrile hydratases towards α,ω-dinitriles

The β37 residue of nitrile hydratase (NHase) from Pseudomonas putida and NHase from Comamonas testosteroni played a critical role in directing enzyme regioselectivity.

Multiple reuses of Rhodococcus ruber TH3 free cells to produce acrylamide in a membrane dispersion microreactor

In this work, multiple reuses of Rhodococcus ruber TH3 free cells for the hydration of acrylonitrile to produce acrylamide in a membrane dispersion microreactor were carried out. Through using a centrifuge, the reactions reached 39.9, 39.5, 38.6 and 38.0wt% of the final acrylamide product concentration respectively within 35min in a four cycle reuse of free cells. In contrast, using a stirring tank, free cells could only be used once with the same addition speed of acrylonitrile with a microreactor. Through observing the dissolution behavior of acrylonitrile microdroplets in a free cell solution using a coaxial microfluidic device and microscope, it was found that the acrylonitrile microdroplets with a diameter of 75μm were rarely observed within a length of 2cm channel within 10s, which illustrated that the microreactor can intensify the reaction rate to reduce the inhibition of acrylonitrile and acrylamide.

Protein diffusion and long-term adsorption states at charged solid surfaces

The diffusion pathways of lysozyme adsorbed to a model charged ionic surface are studied using fully atomistic steered molecular dynamics simulation. The simulations start from existing protein adsorption trajectories, where it has been found that one particular residue, Arg128 at the N,C-terminal face, plays a crucial role in anchoring the lysozyme to the surface [Langmuir 2010 , 26 , 15954 - 15965]. We first investigate the desorption pathway for the protein by pulling the Arg128 side chain away from the surface in the normal direction, and its subsequent readsorption, before studying diffusion pathways by pulling the Arg128 side chain parallel to the surface. We find that the orientation of this side chain plays a decisive role in the diffusion process. Initially, it is oriented normal to the surface, aligning in the electrostatic field of the surface during the adsorption process, but after resorption it lies parallel to the surface, being unable to return to its original orientation due to geometric constraints arising from structured water layers at the surface. Diffusion from this alternative adsorption state has a lower energy barrier of ∼0.9 eV, associated with breaking hydrogen bonds along the pathway, in reasonable agreement with the barrier inferred from previous experimental observation of lysozyme surface clustering. These results show the importance of studying protein diffusion alongside adsorption to gain full insight into the formation of protein clusters and films, essential steps in the future development of functionalized surfaces.

Molecular jamming--the cystine slipknot mechanical clamp in all-atom simulations

A recent survey of 17 134 proteins has identified a new class of proteins which are expected to yield stretching induced force peaks in the range of 1 nN. Such high force peaks should be due to forcing of a slip-loop through a cystine ring, i.e., by generating a cystine slipknot. The survey has been performed in a simple coarse grained model. Here, we perform all-atom steered molecular dynamics simulations on 15 cystine knot proteins and determine their resistance to stretching. In agreement with previous studies within a coarse grained structure based model, the level of resistance is found to be substantially higher than in proteins in which the mechanical clamp operates through shear. The large stretching forces arise through formation of the cystine slipknot mechanical clamp and the resulting steric jamming. We elucidate the workings of such a clamp in an atomic detail. We also study the behavior of five top strength proteins with the shear-based mechanostability in which no jamming is involved. We show that in the atomic model, the jamming state is relieved by moving one amino acid at a time and there is a choice in the selection of the amino acid that advances the first. In contrast, the coarse grained model also allows for a simultaneous passage of two amino acids.