Molecular Basis of the Substrate Specificity of Phosphotriesterase from Pseudomonas diminuta: A Combined QM/MM MD and Electron Density Study.

The occurrence of organophosphorus compounds, pesticides, and flame-retardants in wastes is an emerging ecological problem. Bacterial phosphotriesterases are capable of hydrolyzing some of them. We utilize modern molecular modeling tools to study the hydrolysis mechanism of organophosphorus compounds with good and poor leaving groups by phosphotriesterase from Pseudomonas diminuta (Pd-PTE). We compute Gibbs energy profiles for enzymes with different cations in the active site: native Zn2+cations and Co2+cations, which increase the steady-state rate constant. Hydrolysis occurs in two elementary steps via an associative mechanism and formation of the pentacoordinated intermediate. The first step, a nucleophilic attack, occurs with a low energy barrier independently of the substrate. The second step has a low energy barrier and considerable stabilization of products for substrates with good leaving groups. For substrates with poor leaving groups, the reaction products are destabilized relative to the ES complex that suppresses the reaction. The reaction proceeds with low energy barriers for substrates with good leaving groups with both Zn2+and Co2+cations in the active site; thus, the product release is likely to be a limiting step. Electron density and geometry analysis of the QM/MM MD trajectories of the intermediate states with all considered compounds allow us to discriminate substrates by their ability to be hydrolyzed by the Pd-PTE. For hydrolyzable substrates, the cleaving bond between a phosphorus atom and a leaving group is elongated, and electron density depletion is observed on the Laplacian of electron density maps.

Dynamic features of virus protein 1 and substitutions in the 3-phenyl ring determine the potency and broad-spectrum activity of capsid-binding pyrazolo[3,4-d]pyrimidines against rhinoviruses.

Pyrazolo[3,4-d]pyrimidines represent one potent class of well tolerated and highly active rhinovirus (RV) inhibitors that act as capsid binders. The lead compound OBR-5-340 inhibits a broad-spectrum of RVs. Aiming to improve lead activity, we evaluated the impact of structural modifications in the 3-phenyl ring of OBR-5-340 on its potency and spectrum of anti-RV activity vitro. Our results demonstrate the crucial role of substitution at position 4 for strong, broad-spectrum anti-RV activity. The 4-methyl (RCB23137) and 4-chloro (RCB23138) derivatives outperformed OBR-5-340 in terms of potency and anti-RV activity spectrum. Based on these findings, the compounds were selected for computational binding studies. Molecular dynamic simulations with six RVs differing in OBR-5-340, RCB23137, and RCB23138 sensitivity proved the impact of dynamic features of two VP1 loops enveloping these inhibitors on antiviral potency.

Indole-core inhibitors of influenza a neuraminidase: iterative medicinal chemistry and molecular modeling.

Influenza viruses that cause seasonal and pandemic flu are a permanent health threat. The surface glycoprotein, neuraminidase, is crucial for the infectivity of the virus and therefore an attractive target for flu drug discovery campaigns. We have designed and synthesized more than 40 3-indolinone derivatives. We mainly investigated the role of substituents at the 2 position of the core as well as the introduction of substituents or a nitrogen atom in the fused phenyl ring of the core for inhibition of influenza virus neuraminidase activity and replication in vitro and in vivo. After evaluating the compounds for their ability to inhibit the viral neuraminidase, six potent inhibitors 3c, 3e, 7c, 12o, 12v, 18d were progressed to evaluate for cytotoxicity and inhibition of influenza virus A/PR/8/34 replication in in MDCK cells. Two hit compounds 3e and 12o were tested in an animal model of influenza virus infection. Molecular mechanism of the 3-indolinone derivatives interactions with the neuraminidase was revealed in molecular dynamic simulations. Proposed inhibitors bind to the 430-cavity that is different from the conventional binding site of commercial compounds. The most promising 3-indolinone inhibitors demonstrate stronger interactions with the neuraminidase in molecular models that supports proposed binding site.

Molecular mechanism of thiocyanate dehydrogenase at atomic resolution.

Some sulfur-oxidizing bacteria playing an important role in global geochemical cycles utilize thiocyanate as the sole source of energy and nitrogen. In these bacteria the process of thiocyanate into cyanate conversion is mediated by thiocyanate dehydrogenases - a recently discovered family of copper-containing enzymes with the three­copper active site unique among the other copper proteins. To get a deeper insight into the structure and molecular mechanism of action of thiocyanate dehydrogenases we isolated, purified, and comprehensively characterized an enzyme from the bacterium Pelomicrobium methylotrophicum. High-resolution crystal structures of the thiocyanate dehydrogenase in the free state and in the complexes with the transition state analog, thiourea, and the closest substrate analog, selenocyanate, unveiled the fine details of molecular events occurring at the enzyme active site. During the reaction thiocyanate dehydrogenase undergoes profound conformational change that affects the position of the constituent copper ions and results in the activation of the attacking water molecule. The structure of the enzyme complex with the selenium atom bridged in-between two copper ions was obtained representing an important transient intermediate. Structures of the complexes with inhibitors supplemented with quantum chemical calculations clarify the role of copper ions and refine molecular mechanism of catalysis by thiocyanate dehydrogenase.

Molecular Basis of Influence of A501X Mutations in Penicillin-Binding Protein 2 of Neisseria gonorrhoeae Strain 35/02 on Ceftriaxone Resistance.

The increase in the resistance of mutant strains of Neisseria gonorrhoeae to the antibiotic ceftriaxone is pronounced in the decrease in the second-order acylation rate constant, k2/KS, by penicillin-binding protein 2 (PBP2). These changes can be caused by both the decrease in the acylation rate constant, k2, and the weakening of the binding affinity, i.e., an increase in the substrate constant, KS. A501X mutations in PBP2 affect second-order acylation rate constants. The PBP2A501V variant exhibits a higher k2/KS value, whereas for PBP2A501R and PBP2A501P variants, these values are lower. We performed molecular dynamic simulations with both classical and QM/MM potentials to model both acylation energy profiles and conformational dynamics of four PBP2 variants to explain the origin of k2/KS changes. The acylation reaction occurs in two elementary steps, specifically, a nucleophilic attack by the oxygen atom of the Ser310 residue and C-N bond cleavage in the ß-lactam ring accompanied by the elimination of the leaving group of ceftriaxone. The energy barrier of the first step increases for PBP2 variants with a decrease in the observed k2/KS value. Submicrosecond classic molecular dynamic trajectories with subsequent cluster analysis reveal that the conformation of the ß3-ß4 loop switches from open to closed and its flexibility decreases for PBP2 variants with a lower k2/KS value. Thus, the experimentally observed decrease in the k2/KS in A501X variants of PBP2 occurs due to both the decrease in the acylation rate constant, k2, and the increase in KS.

Saturation Mutagenesis and Molecular Modeling: The Impact of Methionine 182 Substitutions on the Stability of ß-Lactamase TEM-1.

Serine ß-lactamase TEM-1 is the first ß-lactamase discovered and is still common in Gram-negative pathogens resistant to ß-lactam antibiotics. It hydrolyzes penicillins and cephalosporins of early generations. Some of the emerging TEM-1 variants with one or several amino acid substitutions have even broader substrate specificity and resistance to known covalent inhibitors. Key amino acid substitutions affect catalytic properties of the enzyme, and secondary mutations accompany them. The occurrence of the secondary mutation M182T, called a "global suppressor", has almost doubled over the last decade. Therefore, we performed saturating mutagenesis at position 182 of TEM-1 to determine the influence of this single amino acid substitution on the catalytic properties, thermal stability, and ability for thermoreactivation. Steady-state parameters for penicillin, cephalothin, and ceftazidime are similar for all TEM-1 M182X variants, whereas melting temperature and ability to reactivate after incubation at a higher temperature vary significantly. The effects are multidirectional and depend on the particular amino acid at position 182. The M182E variant of ß-lactamase TEM-1 demonstrates the highest residual enzymatic activity, which is 1.5 times higher than for the wild-type enzyme. The 3D structure of the side chain of residue 182 is of particular importance as observed from the comparison of the M182I and M182L variants of TEM-1. Both of these amino acid residues have hydrophobic side chains of similar size, but their residual activity differs by three-fold. Molecular dynamic simulations add a mechanistic explanation for this phenomenon. The important structural element is the V159-R65-E177 triad that exists due to both electrostatic and hydrophobic interactions. Amino acid substitutions that disturb this triad lead to a decrease in the ability of the ß-lactamase to be reactivated.

A highly specific aptamer for the SARS-CoV-2 spike protein from the authentic strain.

DNA aptamers are oligonucleotides that specifically bind to target molecules, similar to how antibodies bind to antigens. We identified an aptamer named MEZ that is highly specific to the receptor-binding domain, RBD, of the SARS-CoV-2 spike protein from the Wuhan-Hu-1 strain. The SELEX procedure was utilized to enrich the initial 31-mer oligonucleotide library with the target aptamer. The aptamer identification was performed using the novel protocol based on nanopore sequencing developed in this study. The MEZ aptamer was chemically synthesized and tested for binding with the SARS-CoV-2 RBD of the spike protein from different strains. The Kd is 6.5 nM for the complex with the RBD from the Wuhan-Hu-1 strain, which is comparable with known aptamers; the advantage is that the MEZ aptamer is smaller than known analogs. The proposed aptamer is highly selective for the RBD protein from the Wuhan-Hu-1 strain and does not form complexes with the RBD from Beta, Delta and Omicron strains. Experimental and theoretical studies together revealed the molecular mechanism of aptamer binding. The aptamer occupies the same binding site as ACE2 when bound to the RBD. The 3'-end of the MEZ aptamer is important for complex formation and is responsible for the discrimination of the RBD protein from a specific strain. The 5'-end is responsible for the formation of a loop in the 3D structure of the aptamer, which is important for proper binding.

The role of the correlated motion(s) of the chromophore in photoswitching of green and red forms of the photoconvertible fluorescent protein mSAASoti.

Wild-type SAASoti and its monomeric variant mSAASoti can undergo phototransformations, including reversible photoswitching of the green form to a nonfluorescent state and irreversible green-to-red photoconversion. In this study, we extend the photochemistry of mSAASoti variants to enable reversible photoswitching of the red form. This result is achieved by rational and site-saturated mutagenesis of the M163 and F177 residues. In the case of mSAASoti it is M163T substitution that leads to the fastest switching and the most photostable variant, and reversible photoswitching can be observed for both green and red forms when expressed in eukaryotic cells. We obtained a 13-fold increase in the switching efficiency with the maximum switching contrast of the green form and the appearance of comparable switching of the red form for the C21N/M163T mSAASoti variant. The crystal structure of the C21N mSAASoti in its green on-state was obtained for the first time at 3.0 Å resolution, and it is in good agreement with previously calculated 3D-model. Dynamic network analysis reveals that efficient photoswitching occurs if motions of the 66H residue and phenyl fragment of chromophore are correlated and these moieties belong to the same community.

Complete genome sequence of new Microbacterium sp. strain ET2, isolated from roots of leafless orchid.

Whole-genome sequence of ET2 strain, isolated from the roots of leafless orchid, constitutes a single circular chromosome of 3,604,840 bp (69.44% G + C content). BLAST+-based average nucleotide identity (ANIb) and digital DNA-DNA hybridization values indicate that ET2 may be a novel Microbacterium species. Genes putatively involved in plant-microbial interactions were predicted.

Molecular mechanism of rhinovirus escape from the Pyrazolo[3,4-d]pyrimidine capsid-binding inhibitor OBR-5-340 via mutations distant from the binding pocket: Derivatives that brake resistance.

Rhinoviruses (RVs) cause the common cold. Attempts at discovering small molecule inhibitors have mainly concentrated on compounds supplanting the medium chain fatty acids residing in the sixty icosahedral symmetry-related hydrophobic pockets of the viral capsid of the Rhinovirus-A and -B species. High-affinity binding to these pockets stabilizes the capsid against structural changes necessary for the release of the ss(+) RNA genome into the cytosol of the host cell. However, single-point mutations may abolish this binding. RV-B5 is one of several RVs that are naturally resistant against the well-established antiviral agent pleconaril. However, RV-B5 is strongly inhibited by the pyrazolopyrimidine OBR-5-340. Here, we report on isolation and characterization of RV-B5 mutants escaping OBR-5-340 inhibition and show that substitution of amino acid residues not only within the binding pocket but also remote from the binding pocket hamper inhibition. Molecular dynamics network analysis revealed that strong inhibition occurs when an ensemble of several sequence stretches of the capsid proteins enveloping OBR-5-340 move together with OBR-5-340. Mutations abrogating this dynamic, regardless of whether being localized within the binding pocket or distant from it result in escape from inhibition. Pyrazolo [3,4-d]pyrimidine derivatives overcoming OBR-5-340 escape of various RV-B5 mutants were identified. Our work contributes to the understanding of the properties of capsid-binding inhibitors necessary for potent and broad-spectrum inhibition of RVs.

Histidine-assisted reduction of arylnitrenes upon photo-activation of phenyl azide chromophores in GFP-like fluorescent proteins.

The photochemically active sites of the proteins sfGFP66azF and Venus66azF, members of the green fluorescent protein (GFP) family, contain a non-canonical amino acid residue p-azidophenylalanine (azF) instead of Tyr66. The light-induced decomposition of azF at these sites leads to the formation of reactive arylnitrene (nF) intermediates followed by the formation of phenylamine-containing chromophores. We report the first study of the reaction mechanism of the reduction of the arylnitrene intermediates in sfGFP66nF and Venus66nF using molecular modeling methods. The Gibbs energy profiles for the elementary steps of the chemical reaction in sfGFP66nF are computed using molecular dynamics simulations with quantum mechanics/molecular mechanics (QM/MM) potentials. Structures and energies along the reaction pathway in Venus66nF are evaluated using a QM/MM approach. According to the results of the simulations, arylnitrene reduction is coupled with oxidation of the histidine side chain on the His148 residue located near the chromophore.

Structural framework for the understanding spectroscopic and functional signatures of the cyanobacterial Orange Carotenoid Protein families.

The Orange Carotenoid Protein (OCP) is a unique photoreceptor crucial for cyanobacterial photoprotection. Best studied Synechocystis sp. PCC 6803 OCP belongs to the large OCP1 family. Downregulated by the Fluorescence Recovery Protein (FRP) in low-light, high-light-activated OCP1 binds to the phycobilisomes and performs non-photochemical quenching. Recently discovered families OCP2 and OCP3 remain structurally and functionally underexplored, and no systematic comparative studies have ever been conducted. Here we present two first crystal structures of OCP2 from morphoecophysiologically different cyanobacteria and provide their comprehensive structural, spectroscopic and functional comparison with OCP1, the recently described OCP3 and all-OCP ancestor. Structures enable correlation of spectroscopic signatures with the effective number of hydrogen and discovered here chalcogen bonds anchoring the ketocarotenoid in OCP, as well as with the rotation of the echinenone's ß-ionone ring in the CTD. Structural data also helped rationalize the observed differences in OCP/FRP and OCP/phycobilisome functional interactions. These data are expected to foster OCP research and applications in optogenetics, targeted carotenoid delivery and cyanobacterial biomass engineering.

KuINins as a New Class of HIV-1 Inhibitors That Block Post-Integration DNA Repair.

Integration of HIV-1 genomic cDNA results in the formation of single-strand breaks in cellular DNA, which must be repaired for efficient viral replication. Post-integration DNA repair mainly depends on the formation of the HIV-1 integrase complex with the Ku70 protein, which promotes DNA-PK assembly at sites of integration and its activation. Here, we have developed a first-class inhibitor of the integrase-Ku70 complex formation that inhibits HIV-1 replication in cell culture by acting at the stage of post-integration DNA repair. This inhibitor, named s17, does not affect the main cellular function of Ku70, namely its participation in the repair of double-strand DNA breaks through the non-homologous end-joining pathway. Using a molecular dynamics approach, we have constructed a model for the interaction of s17 with Ku70. According to this model, the interaction of two phenyl radicals of s17 with the L76 residue of Ku70 is important for this interaction. The requirement of two phenyl radicals in the structure of s17 for its inhibitory properties was confirmed using a set of s17 derivatives. We propose to stimulate compounds that inhibit post-integration repair by disrupting the integrase binding to Ku70 KuINins.

In search for structural targets for engineering d-amino acid transaminase: modulation of pH optimum and substrate specificity.

The development of biocatalysts requires reorganization of the enzyme's active site to facilitate the productive binding of the target substrate and improve turnover number at desired conditions. Pyridoxal-5'-phosphate (PLP) - dependent transaminases are highly efficient biocatalysts for asymmetric amination of ketones and keto acids. However, transaminases, being stereoselective enzymes, have a narrow substrate specificity due to the ordered structure of the active site and work only in neutral-alkaline media. Here, we investigated the d-amino acid transaminase from Aminobacterium colombiense, with the active site organized differently from that of the canonical d-amino acid transaminase from Bacillus sp. YM-1. Using a combination of site-directed mutagenesis, kinetic analysis, molecular modeling, and structural analysis we determined the active site residues responsible for substrate binding, substrate differentiation, thermostability of a functional dimer, and affecting the pH optimum. We demonstrated that the high specificity toward d-glutamate/α-ketoglutarate is due to the interactions of a γ-carboxylate group with K237 residue, while binding of other substrates stems from the effectiveness of their accommodation in the active site optimized for d-glutamate/α-ketoglutarate binding. Furthermore, we showed that the K237A substitution shifts the catalytic activity optimum to acidic pH. Our findings are useful for achieving target substrate specificity and demonstrate the potential for developing and optimizing transaminases for various applications.

Pyrrolo[2,3-e]indazole as a novel chemotype for both influenza A virus and pneumococcal neuraminidase inhibitors.

Influenza infections are often exacerbated by secondary bacterial infections, primarily caused by Streptococcus pneumoniae. Both respiratory pathogens have neuraminidases that support infection. Therefore, we hypothesized that dual inhibitors of viral and bacterial neuraminidases might be an advantageous strategy for treating seasonal and pandemic influenza pneumonia complicated by bacterial infections. By screening our in-house chemical library, we discovered a new chemotype that may be of interest for a further campaign to find small molecules against influenza. Our exploration of the pyrrolo[2,3-e]indazole space led to the identification of two hit compounds, 6h and 12. These molecules were well-tolerated by MDCK cells and inhibited the replication of H3N2 and H1N1 influenza A virus strains. Moreover, both compounds suppress viral and pneumococcal neuraminidases indicating their dual activity. Given its antiviral activity, pyrrolo[2,3-e]indazole has been identified as a promising scaffold for the development of novel neuraminidase inhibitors that are active against influenza A virus and S. pneumoniae.

Multiscale Simulations of the Covalent Inhibition of the SARS-CoV-2 Main Protease: Four Compounds and Three Reaction Mechanisms.

We report the results of computational modeling of the reactions of the SARS-CoV-2 main protease (MPro) with four potential covalent inhibitors. Two of them, carmofur and nirmatrelvir, have shown experimentally the ability to inhibit MPro. Two other compounds, X77A and X77C, were designed computationally in this work. They were derived from the structure of X77, a non-covalent inhibitor forming a tight surface complex with MPro. We modified the X77 structure by introducing warheads capable of reacting with the catalytic cysteine residue in the MPro active site. The reaction mechanisms of the four molecules with MPro were investigated by quantum mechanics/molecular mechanics (QM/MM) simulations. The results show that all four compounds form covalent adducts with the catalytic cysteine Cys 145 of MPro. From the chemical perspective, the reactions of these four molecules with MPro follow three distinct mechanisms. The reactions are initiated by a nucleophilic attack of the thiolate group of the deprotonated cysteine residue from the catalytic dyad Cys145-His41 of MPro. In the case of carmofur and X77A, the covalent binding of the thiolate to the ligand is accompanied by the formation of the fluoro-uracil leaving group. The reaction with X77C follows the nucleophilic aromatic substitution SNAr mechanism. The reaction of MPro with nirmatrelvir (which has a reactive nitrile group) leads to the formation of a covalent thioimidate adduct with the thiolate of the Cys145 residue in the enzyme active site. Our results contribute to the ongoing search for efficient inhibitors of the SARS-CoV-2 enzymes.

Combined Structural and Computational Study of the mRubyFT Fluorescent Timer Locked in Its Blue Form.

The mRubyFT is a monomeric genetically encoded fluorescent timer based on the mRuby2 fluorescent protein, which is characterized by the complete maturation of the blue form with the subsequent conversion to the red one. It has higher brightness in mammalian cells and higher photostability compared with other fluorescent timers. A high-resolution structure is a known characteristic of the mRubyFT with the red form chromophore, but structural details of its blue form remain obscure. In order to obtain insight into this, we obtained an S148I variant of the mRubyFT (mRubyFTS148I) with the blocked over time blue form of the chromophore. X-ray data at a 1.8 Å resolution allowed us to propose a chromophore conformation and its interactions with the neighboring residues. The imidazolidinone moiety of the chromophore is completely matured, being a conjugated π-system. The methine bridge is not oxidized in the blue form bringing flexibility to the phenolic moiety that manifests itself in poor electron density. Integration of these data with the results of molecular dynamic simulation disclosed that the OH group of the phenolic moiety forms a hydrogen bond with the side chain of the T163 residue. A detailed comparison of mRubyFTS148I with other available structures of the blue form of fluorescent proteins, Blue102 and mTagBFP, revealed a number of characteristic differences. Molecular dynamic simulations with the combined quantum mechanic/molecular mechanic potentials demonstrated that the blue form exists in two protonation states, anion and zwitterion, both sharing enolate tautomeric forms of the C=C-O- fragment. These two forms have similar excitation energies, as evaluated by calculations. Finally, excited state molecular dynamic simulations showed that excitation of the chromophore in both protonation states leads to the same anionic fluorescent state. The data obtained shed light on the structural features and spectral properties of the blue form of the mRubyFT timer.

To the Understanding of Catalysis by D-Amino Acid Transaminases: A Case Study of the Enzyme from Aminobacterium colombiense.

Pyridoxal-5'-phosphate (PLP)-dependent transaminases are highly efficient biocatalysts for stereoselective amination. D-amino acid transaminases can catalyze stereoselective transamination producing optically pure D-amino acids. The knowledge of substrate binding mode and substrate differentiation mechanism in D-amino acid transaminases comes down to the analysis of the transaminase from Bacillus subtilis. However, at least two groups of D-amino acid transaminases differing in the active site organization are known today. Here, we present a detailed study of D-amino acid transaminase from the gram-negative bacterium Aminobacterium colombiense with a substrate binding mode different from that for the transaminase from B. subtilis. We study the enzyme using kinetic analysis, molecular modeling, and structural analysis of holoenzyme and its complex with D-glutamate. We compare the multipoint binding of D-glutamate with the binding of other substrates, D-aspartate and D-ornithine. QM/MM MD simulation reveals that the substrate can act as a base and its proton can be transferred from the amino group to the α-carboxylate group. This process occurs simultaneously with the nucleophilic attack of the PLP carbon atom by the nitrogen atom of the substrate forming gem-diamine at the transimination step. This explains the absence of the catalytic activity toward (R)-amines that lack an α-carboxylate group. The obtained results clarify another substrate binding mode in D-amino acid transaminases and underpinned the substrate activation mechanism.

Extracellular Fe(III) reductase structure reveals a modular organization enabling S-layer insertion and electron transfer to insoluble substrates.

The thermophilic anaerobic Gram-positive bacterium Carboxydothermus ferrireducens utilizes insoluble Fe(III) oxides as electron acceptors in respiratory processes using an extracellular 11-heme cytochrome c OmhA as a terminal reductase. OmhA is able to transfer electrons to soluble and insoluble Fe(III) compounds, substrates of multiheme oxidoreductases, and soluble electron shuttles. The crystal structure of OmhA at 2.5 Å resolution shows that it consists of two functionally distinct parts: the cytochrome с electron transfer and the S-layer binding domains. Nonaheme C-terminal subdomain of the cytochrome с domain is structurally similar to the extracellular multiheme cytochrome OcwA from the metal-reducing Gram-positive bacterium "Thermincola potens." S-layer binding domain of OmhA is responsible for interaction with the S-layer that surrounds the Carboxydothermus ferrireducens cell envelope. The structural foundations enabling the embedding of extracellular multiheme cytochromes to the S-layer of a Gram-positive-type cell wall and putative electron transfer pathways to insoluble minerals are discussed.

Influence of the Active Site Flexibility on the Efficiency of Substrate Activation in the Active Sites of Bi-Zinc Metallo-ß-Lactamases.

The influence of the active site flexibility on the efficiency of catalytic reaction is studied by taking two members of metallo-ß-lactamases, L1 and NDM-1, with the same substrate, imipenem. Active sites of these proteins are covered by L10 loops, and differences in their amino acid compositions affect their rigidity. A more flexible loop in the NDM-1 brings additional flexibility to the active site in the ES complex. This is pronounced in wider distributions of key interatomic distances, such as the distance of the nucleophilic attack, coordination bond lengths, and covalent bond lengths in the substrate. Substrate activation, quantified by Fukui electrophilicity index of the carbonyl carbon atom of the substrate, is also sensitive to the active site flexibility. In the tighter and more rigid L1 enzyme-substrate complex, the substrate is activated more efficiently. In the NDM-1 containing system, only one third of the states are activated to the same extent. Other fractions demonstrate lower substrate activation. Efficiency of the substrate activation and rigidity of the ES complex influence the following chemical reaction. In the more rigid L1-containing system, the reaction barrier of the first step of the reaction is lower, and the first intermediate is more stabilized compared to the NDM-1 containing system.