Near-attack Conformation Research Articles

Computation-Guided Rational Design of Cysteine-Less Protein Variants in Engineered hCGL.

The engineered human cystathionine-γ-lyase (hCGL) resulting in enhanced activity toward both cysteine and cystine unveils a potential robust antitumor activity. However, the presence of cysteine residues has the potential to induce oligomerization or incorrect disulfide bonding, which may decrease the bioavailability of biopharmaceuticals. Through a meticulous design process targeting the cysteine residues within engineered hCGL, a set of potential beneficial mutants were obtained by virtual screening employing Rosetta and ABACUS. Experimental measurements have revealed that most of the mutants showed increased activity toward both substrates l-Cys and CSSC. Furthermore, mutants C109V and C229D demonstrated Tm value increases of 8.2 and 1.8 °C, respectively. After an 80 min incubation at 60 °C, mutant C229D still maintained high residual activity. Unexpectedly, mutant C109V, displaying activity approximately 2-fold higher than the activity of wild type (WT) for both substrates, showed disappointing instability in plasma, which suggests that computational design still requires further consideration. Analysis of their structure and molecular dynamics (MD) simulation revealed the impact of hydrophobic interaction, hydrogen bonds, and near-attack conformation (NAC) stability on activity and stability. This study acquired information about mutants that exhibit enhanced activity or thermal resistance and serve as valuable guidance for subsequent specific cysteine modifications.

Modeling of a near-attack conformation of oxime in phosphorylated acetylcholinesterase via a reactivation product, a phosphorylated oxime

At the present, only four antidotes are in use in therapy for poisoning by organophosphorus compounds: 2-PAM, HI-6, obidoxime and trimedoxime. Numerous compounds have been designed and synthetized to be more effective reactivators than those currently in use. Many of those new compounds fail at the enzyme level because interactions formed within the AChE active site are not favourable ones that lead to a successful reactivation. The approach in which the modeling of a phosphorylated oxime (POX), a product of successful reactivation in the AChE active site, may be a way to better understand the role of active site residues during the process of formation of the Michaelis type of complex between an enzyme and oxime. After reactivation, a change in phosphorus stereochemistry occurs leading to a different spatial arrangement of attached substituents, now including an oxime. To study interactions between the AChE oxyanion hole and a phosphorylated oxime, an S203G mutant was used to avoid the steric hindrance caused by the catalytic serine. In this way, the POX could be positioned close to the oxyanion hole. In the final step, the oxime without a phosphoester moiety was transferred into the phosphorylated AChE and molecular dynamics was used to test the stability of the near-attack conformation of the oxime near the phosphorylated serine.

Molecular modeling of Mannich phenols as reactivators of human acetylcholinesterase inhibited by A-series nerve agents

The A-series is the most recent generation of chemical warfare nerve agents (CWA) which act directly on the inhibition of the human acetylcholinesterase (HssAChE) enzyme. These compounds lack accurate experimental data on their physicochemical properties, and there is no evidence that traditional antidotes effectively reactivate HssAChE inhibited by them. In the search for potential antidotes, we employed virtual screening, molecular docking, and molecular dynamics (MD) simulations for the theoretical assessment of the performance of a library of Mannich phenols as potential reactivators of HssAChE inhibited by the Novichok agents A-230, A-232, and A-234, in comparison with the commercial oximes pralidoxime (2-PAM), asoxime (HI-6), trimedoxime (TMB-4), and obidoxime. Following the near-attack conformation (NAC) approach, our results suggest that the compounds assessed would face difficulties in triggering the proposed nucleophilic in-line displacement mechanism. Despite this, it was observed that certain Mannich phenols presented similar or superior results to those obtained by reference oximes against A-232 and A-234 model, suggesting that these compounds can adopt more favourable conformations. Additional binding energy calculations confirmed the stability of the model/ligands complexes and the reactivating potential observed in the molecular docking and MD studies. Our findings indicate that the Mannich phenols could be alternative antidotes and that their efficacy should be evaluated experimentally against the A-series CWA.

Applications of the Near Attack Conformation (NAC) approach in the search for Acetylcholinesterase reactivators

The Near Attack Conformation (NAC) approach states that the efficiency of an enzyme-catalyzed reaction depends on the prior attainment of optimal conditions for substrate atom organization and positioning for bond formation. These conditions are prerequisites for the transition state (TS) in which the involved atoms are within the van der Waals range of contact and positioned at an angle similar to that achieved after bond formation. The successful application of this approach to investigate the reactivation mechanism of acetylcholinesterase inhibited by nerve agents has contributed to a better understanding of this mechanism and demonstrated consistent corroboration with experimental data. In this article, we summarize the accomplishments achieved thus far and outline future perspectives.

Rational Design of Cyclodextrin Glycosyltransferase with Improved Hesperidin Glycosylation Activity

Cyclodextrin glycosyltransferase (CGTase) can catalyze the glycosylation of hesperidin, resulting in α-glycosyl hesperidin with significantly improved water solubility. In this study, a rational design of CGTase to improve its hesperidin glycosylation activity was investigated. The strategy we employed involved docking hesperidin in its near-attack conformation and virtually mutating the surrounding residues, followed by calculating the changes in binding energy using Rosetta flex-ddG. The mutations with a stabilization effect were then subjected to an activity assay. Starting from CGTase-Y217F, we obtained three double-point mutants, Y217F/M351F, Y217F/M351L, and Y217F/D393H, with improved hesperidin glycosylation activities after screening twenty variants. The best variant, Y217F/D393H, exhibited a catalytic activity of 1305 U/g, and its kcat/KmA is 2.36 times higher compared to CGTase-Y217F and 15.14 times higher compared to the wild-type CGTase. Molecular dynamic simulations indicated that hesperidin was repulsed by CGTase-Y217F when bound in a near-attack conformation. However, by introducing a second-point mutation with a stabilization effect, the repulsion effect is weakened, resulting in a reduction in the distances between the bond-forming atoms and, thus, favoring the reaction.

Protoglobin-Catalyzed Formation of cis-Trifluoromethyl-Substituted Cyclopropanes by Carbene Transfer.

Trifluoromethyl-substituted cyclopropanes (CF3 -CPAs) constitute an important class of compounds for drug discovery. While several methods have been developed for synthesis of trans-CF3 -CPAs, stereoselective production of corresponding cis-diastereomers remains a formidable challenge. We report a biocatalyst for diastereo- and enantio-selective synthesis of cis-CF3 -CPAs with activity on a variety of alkenes. We found that an engineered protoglobin from Aeropyrnum pernix (ApePgb) can catalyze this unusual reaction at preparative scale with low-to-excellent yield (6-55 %) and enantioselectivity (17-99 % ee), depending on the substrate. Computational studies revealed that the steric environment in the active site of the protoglobin forced iron-carbenoid and substrates to adopt a pro-cis near-attack conformation. This work demonstrates the capability of enzyme catalysts to tackle challenging chemistry problems and provides a powerful means to expand the structural diversity of CF3 -CPAs for drug discovery.

Theoretical assessment of the performances of commercial oximes on the reactivation of acetylcholinesterase inhibited by the nerve agent A-242 (novichok)

The nerve agents of the A-series are relatively recent chemical weapons with no antidote available yet. Once inside the human body, those chemicals act similarly to the classic nerve agents, by binding to the catalytic residue Serine 203 (Ser203) of human acetylcholinesterase (HssAChE) and thus preventing the proper function of this enzyme. However, there is no experimental evidence yet if the current antidotes for intoxication by nerve agents are also capable of restoring AChE inhibited by the nerve agents of the A-series. In order to launch some light on this issue, we used computational techniques (molecular docking, molecular dynamics and MM-PBSA interaction energy calculations) to assess the performances of the four currently available commercial oximes (2-PAM, HI-6, obidoxime and trimedoxime) when in contact with HssAChE inhibited by the agent A-242. Based on the near-attack conformation (NAC) criterion, our results suggest that the commercial oximes would have limited efficacy to reactivate the enzyme since they are not able to properly approach the adduct Ser203-A-242. Among those oximes, trimedoxime seems to be the most promising, since it showed lower values of energy in the MM-PBSA calculations, a higher stability inside the catalytic anionic center (CAS) of HssAChE, and was able to adopt a position closer to the NAC that could enable the reactivation mechanism.

The Relationship between Enzyme Conformational Change, Proton Transfer, and Phosphoryl Transfer in β-Phosphoglucomutase

Molecular details for the timing and role of proton transfer in phosphoryl transfer reactions are poorly understood. Here, we have combined QM models, experimental NMR measurements, and X-ray structures to establish that the transition of an archetypal phosphoryl transfer enzyme, βPGM, from a very closed near-attack conformation to a fully closed transition state analogue (TSA) conformation triggers both partial proton transfer from the general acid–base residue to the leaving group oxygen and partial dissociation of the transferring phosphoryl group from the leaving group oxygen. Proton transfer continues but is not completed throughout the reaction path of the phosphoryl transfer with the enzyme in the TSA conformation. Moreover, using interacting quantum atoms (IQA) and relative energy gradient (REG) analysis approaches, we observed that the change in the position of the proton and the corresponding increased electrostatic repulsion between the proton and the phosphorus atom provide a stimulus for phosphoryl transfer in tandem with a reduction in the negative charge density on the leaving group oxygen atom. The agreement between solution-phase 19F NMR measurements and equivalent QM models of βPGMWT and βPGMD10N TSA complexes confirms the protonation state of G6P in the two variants, validating the employed QM models. Furthermore, QM model predictions of an AlF4 distortion in response to the proton position are confirmed using high resolution X-ray crystal structures, not only providing additional validation to the QM models but also further establishing metal fluorides as highly sensitive experimental predictors of active-site charge density distributions.

An Energy Optimization Strategy Based on the Perfect Conformation of Prolyl Endopeptidase for Improving Catalytic Efficiency.

Prolyl endopeptidases (PEPs) hydrolyze proteins to yield bioactive peptides and are effective in the treatment of celiac disease. However, the catalytic efficiency of PEPs still has the potential to be improved, which could further strengthen their industrial and therapeutic applications. Herein, a novel rational design strategy based on a "near-attack conformation" of the catalytic state of PEP was adopted. Constrained dynamic simulations were applied, followed by the virtual screening of potentially favorable mutants according to their binding free energy. We redesigned Sphaerobacter thermophiles PEP with high-temperature activity/stability, a wide range of pH stabilities, and high proline specificity. As a result, the kcat value of two PEP mutants (I462W and Q560Y) increased by 208.2 and 150.1%, respectively, and the kcat/KM increased by 32.7 and 6.3%, respectively. These data revealed that the PEP mutants had improved catalytic efficiency and that our strategy can be applied for enzyme engineering.

Lipase mediated enzymatic kinetic resolution of phenylethyl halohydrins acetates: A case of study and rationalization

Racemic phenylethyl halohydrins acetates containing several groups attached to the aromatic ring were resolved via hydrolysis reaction in the presence of lipase B from Candida antarctica (Novozym® 435). In all cases, the kinetic resolution was highly selective (E > 200) leading to the corresponding (S)-β-halohydrin with ee > 99 %. However, the time required for an ideal 50 % conversion ranged from 15 min for 2,4-dichlorophenyl chlorohydrin acetate to 216 h for 2-chlorophenyl bromohydrin acetate. Six chlorohydrins and five bromohydrins were evaluated, the latter being less reactive. For the β-brominated substrates, steric hindrance on the aromatic ring played a crucial role, which was not observed for the β-chlorinated derivatives. To shed light on the different reaction rates, docking studies were carried out with all the substrates using MD simulations. The computational data obtained for the β-brominated substrates, based on the parameters analysed such as NAC (near attack conformation), distance between Ser-O and carbonyl-C and oxyanion site stabilization were in agreement with the experimental results. On the other hand, the data obtained for β-chlorinated substrates suggested that physical aspects such as high hydrophobicity or induced change in the conformation of the enzymatic active site are more relevant aspects when compared to steric hindrance effects.

Reactivation of VX-Inhibited Human Acetylcholinesterase by Deprotonated Pralidoxime. A Complementary Quantum Mechanical Study.

In the present work, we performed a complementary quantum mechanical (QM) study to describe the mechanism by which deprotonated pralidoxime (2-PAM) could reactivate human (Homo sapiens sapiens) acetylcholinesterase (HssAChE) inhibited by the nerve agent VX. Such a reaction is proposed to occur in subsequent addition–elimination steps, starting with a nucleophile bimolecular substitution (SN2) mechanism through the formation of a trigonal bipyramidal transition state (TS). A near attack conformation (NAC), obtained in a former study using molecular mechanics (MM) calculations, was taken as a starting point for this project, where we described the possible formation of the TS. Together, this combined QM/MM study on AChE reactivation shows the feasibility of the reactivation occurring via attack of the deprotonated form of 2-PAM against the Ser203-VX adduct of HssAChE.

Laser Ablation Assists Cyclization Reactions of Hydantoic Acid: A Proof for the Near-Attack Conformation Theory?

In the course of the investigation of the rotational spectrum of prebiotic hydantoic acid by Fourier transform microwave spectroscopy coupled to a laser ablation source in a supersonic expansion, rotational signatures of two cyclic molecules, hydantoin and 2,5-oxazolidinedione, have been unexpectedly observed along with the four most stable conformers of hydantoic acid. Interestingly, two of them presented folded geometric arrangements that might act as precursors in the cyclization reactions assisted by laser ablation. They could play the role of near-attack conformations (NACs) in the framework of the NAC theory for intramolecular reactions. A detailed analysis of the spectrum further revealed the simultaneous formation of other species in the jet, showing that the laser ablation of solid organic precursors constitutes an alternative tool in the generation of new chemical species.

Molecular modeling studies on the interactions of 7-methoxytacrine-4-pyridinealdoxime, 4-PA, 2-PAM, and obidoxime with VX-inhibited human acetylcholinesterase: a near attack conformation approach

7-methoxytacrine-4-pyridinealdoxime (7-MEOTA-4-PA, named hybrid 5C) is a compound formerly synthesized and evaluated in vitro, together with 4-pyridine aldoxime (4-PA) and commercial reactivators of acetylcholinesterase (AChE). This compound was designed with the purpose of being a prophylactic reactivator, capable of interacting with different subdomains of the active site of AChE. To investigate these interactions, theoretical results from docking were first compared with experimental data of hybrid 5C, 4-PA, and two commercial oximes, on the reactivation of human AChE (HssAChE) inhibited by VX. Then, further docking studies, molecular dynamics simulations, and molecular mechanics Poisson–Boltzmann surface area calculations, were carried out to investigate reactivation performances, considering the near attack conformation (NAC) approach, prior to the nucleophilic substitution mechanism. Our results helped to elucidate the interactions of such molecules with the different subdomains of the active site of HssAChE. Additionally, NAC poses of each oxime were suggested for further theoretical studies on the reactivation reaction.

Simple Design of an Enzyme-Inspired Supported Catalyst Based on a Catalytic Triad

Summary Enzyme active sites afford an intricate interplay of functional groups to mediate complex organic and inorganic reactions. Many hydrolytic enzymes use a catalytic triad comprising three different functional residues—(Ser(-OH), Hist(-imidazole), Asp(-CO 2 H))—that catalyze the hydrolysis of numerous unique substrates. Inspired by this design, we have developed a simple one-step synthesis for preparing a new supported catalytic system in which the three reactive groups of the catalytic triad (alcohol, imidazole, and carboxylate) are incorporated into a single functional unit. These artificial active sites can be coupled to a solid-phase support (Merrifield resin) by copper(I)-catalyzed azide-alkyne cycloaddition "click chemistry," and their effectiveness as esterolysis catalysts was demonstrated. Furthermore, tuning the local hydrophobicity of the resin particles with an approach analogous to the native enzyme hydrophobic pocket increased the catalytic efficiency. Quantum mechanics and molecular dynamics computational modeling were used to probe the catalytic effect and suggested a concerted two-step mechanism and hydrophobic nanoenvironment similar to that of hydrolytic enzymes.

Large-Scale Structure-Based Prediction and Identification of Novel Protease Substrates Using Computational Protein Design

Characterizing the substrate specificity of protease enzymes is critical for illuminating the molecular basis of their diverse and complex roles in a wide array of biological processes. Rapid and accurate prediction of their extended substrate specificity would also aid in the design of custom proteases capable of selectively and controllably cleaving biotechnologically or therapeutically relevant targets. However, current in silico approaches for protease specificity prediction, rely on, and are therefore limited by, machine learning of sequence patterns in known experimental data. Here, we describe a general approach for predicting peptidase substrates de novo using protein structure modeling and biophysical evaluation of enzyme–substrate complexes. We construct atomic resolution models of thousands of candidate substrate–enzyme complexes for each of five model proteases belonging to the four major protease mechanistic classes—serine, cysteine, aspartyl, and metallo-proteases—and develop a discriminatory scoring function using enzyme design modules from Rosetta and AMBER's MMPBSA. We rank putative substrates based on calculated interaction energy with a modeled near-attack conformation of the enzyme active site. We show that the energetic patterns obtained from these simulations can be used to robustly rank and classify known cleaved and uncleaved peptides and that these structural-energetic patterns have greater discriminatory power compared to purely sequence-based statistical inference. Combining sequence and energetic patterns using machine-learning algorithms further improves classification performance, and analysis of structural models provides physical insight into the structural basis for the observed specificities. We further tested the predictive capability of the model by designing and experimentally characterizing the cleavage of four novel substrate motifs for the hepatitis C virus NS3/4 protease using an in vivo assay. The presented structure-based approach is generalizable to other protease enzymes with known or modeled structures, and complements existing experimental methods for specificity determination.

ChemInform Abstract: Theoretical Study on Highly Active Bifunctional Metalloporphyrin Catalysts for the Coupling Reaction of Epoxides with Carbon Dioxide

AbstractReview: [12 refs.

A Semiautomated Structure-Based Method To Predict Substrates of Enzymes via Molecular Docking: A Case Study with Candida antarctica Lipase B.

The discovery of unique substrates is important for developing potential applications of enzymes. However, the experimental procedures for substrate identification are laborious, time-consuming, and expensive. Although in silico structure-based approaches show great promise, recent extensive studies have shown that these approaches remain a formidable challenge for current biocomputational methodologies. Here we present an open-source, extensible, and flexible software platform for predicting enzyme substrates called THEMIS, which performs in silico virtual screening for potential catalytic targets of an enzyme on the basis of the enzyme's catalysis mechanism. On the basis of a generalized transition state theory of enzyme catalysis, we introduce a modified docking procedure called "mechanism-based restricted docking" (MBRD) for novel substrate recognition from molecular docking. Comprising a series of utilities written in C/Python, THEMIS automatically executes parallel-computing MBRD tasks and evaluates the results with various molecular mechanics (MM) criteria such as energy, distance, angle, and dihedral angle to help identify desired substrates. Exhaustive sampling and statistical measures were used to improve the robustness and reproducibility of the method. We used Candida antarctica lipase B (CALB) as a test system to demonstrate the effectiveness of our computational prediction of (non)substrates. A novel MM score function for CALB substrate identification derived from the near-attack conformation was used to evaluate the possibility of chemical transformation. A highly positive rate of 93.4% was achieved from a CALB substrate library with 61 known substrates and 35 nonsubstrates, and the screening rate has reached 103 compounds/day (96 CPU cores, 100 samples/compound). The performance shows that the present method is perhaps the first reported scheme to meet the requirement for practical applicability to enzyme studies. An additional study was performed to validate the universality of our method. In this verification we employed two distinct enzymes, nitrilase Nit6803 and SDR Gox2181, where the correct rates of both enzymes exceeded 90%. The source code used will be released under the GNU General Public License (GPLv3) and will be free to download. We believe that the present method will provide new insights into enzyme research and accelerate the development of novel enzyme applications.

Theoretical Study on Highly Active Bifunctional Metalloporphyrin Catalysts for the Coupling Reaction of Epoxides with Carbon Dioxide.

Highly active bifunctional metalloporphyrin catalysts were developed for the coupling reaction of epoxides with CO2 to produce cyclic carbonates. The bifunctional catalysts have both quaternary ammonium halide groups and a metal center. To elucidate the roles of these catalytic groups, DFT calculations were performed. Control reactions using tetrabutylammonium halide as a catalyst were also investigated for comparison. In the present article, the results of our computational studies are overviewed. The computational results are consistent with the experimental data and are useful for elucidating the structure-activity relationship. The key features responsible for the high catalytic activity of the bifunctional catalysts are as follows: 1) the cooperative action of the halide anion (nucleophile) and the metal center (Lewis acid); 2) the near-attack conformation, leading to the efficient opening of the epoxide ring in the rate-determining step; and 3) the conformational change of the quaternary ammonium cation to stabilize various anionic species generated during catalysis, in addition to the robustness (thermostability) of the catalysts.

Bifunctional Homodimeric Triokinase/FMN Cyclase: CONTRIBUTION OF PROTEIN DOMAINS TO THE ACTIVITIES OF THE HUMAN ENZYME AND MOLECULAR DYNAMICS SIMULATION OF DOMAIN MOVEMENTS

Mammalian triokinase, which phosphorylates exogenous dihydroxyacetone and fructose-derived glyceraldehyde, is neither molecularly identified nor firmly associated to an encoding gene. Human FMN cyclase, which splits FAD and other ribonucleoside diphosphate-X compounds to ribonucleoside monophosphate and cyclic X-phosphodiester, is identical to a DAK-encoded dihydroxyacetone kinase. This bifunctional protein was identified as triokinase. It was modeled as a homodimer of two-domain (K and L) subunits. Active centers lie between K1 and L2 or K2 and L1: dihydroxyacetone binds K and ATP binds L in different subunits too distant (≈ 14 Å) for phosphoryl transfer. FAD docked to the ATP site with ribityl 4'-OH in a possible near-attack conformation for cyclase activity. Reciprocal inhibition between kinase and cyclase reactants confirmed substrate site locations. The differential roles of protein domains were supported by their individual expression: K was inactive, and L displayed cyclase but not kinase activity. The importance of domain mobility for the kinase activity of dimeric triokinase was highlighted by molecular dynamics simulations: ATP approached dihydroxyacetone at distances below 5 Å in near-attack conformation. Based upon structure, docking, and molecular dynamics simulations, relevant residues were mutated to alanine, and kcat and Km were assayed whenever kinase and/or cyclase activity was conserved. The results supported the roles of Thr(112) (hydrogen bonding of ATP adenine to K in the closed active center), His(221) (covalent anchoring of dihydroxyacetone to K), Asp(401) and Asp(403) (metal coordination to L), and Asp(556) (hydrogen bonding of ATP or FAD ribose to L domain). Interestingly, the His(221) point mutant acted specifically as a cyclase without kinase activity.

Molecular modeling of different substrate-binding modes and their role in penicillin acylase catalysis

Molecular modeling was addressed to understand different substrate-binding modes and clarify the role of two positively charged residues of the penicillin G acylase active site - βR263 and αR145 - in binding of negatively charged substrates. Although the electrostatic contribution to productive substrate binding was dominated by βR263 rather than αR145, it was found that productive binding was not the only possible mode of substrate placement in the active site. Two extra binding modes - nonproductive and preproductive - were located by means of molecular docking and dynamics with binding affinities comparable with the productive one. A unique feature of nonproductive and preproductive complexes was that the substrate's acyl group did not penetrate the hydrophobic pocket, but occupied a patch on the protein interface spanning from βR263 to αR145. Nonproductive and preproductive complexes competed with each other and productive binding mode, giving rise to increased apparent substrate binding. Preproductive complex revealed an ability to switch to a productive one during molecular dynamics simulations, and conformational plasticity of the penicillin G acylase active site was shown to be crucial for that. Nonproductive binding observed at molecular modeling corresponded well with experimentally observed substrate inhibition in penicillin acylase catalysis. By combining estimated free energies of substrate binding in each mode, and accounting for two possible conformations of the penicillin G acylase active site (closed and open) quantitative agreement with experimentally measured K(M) values was achieved. Calculated near-attack conformation frequencies from corresponding molecular dynamics simulations were in a quantitative correlation with experimental k(cat) values and demonstrated adequate application of molecular modeling methods.