Mechanism of C-terminal intein cleavage in protein splicing from QM/MM molecular dynamics simulations

Molecular docking and MD simulations reveal protease inhibitors block the catalytic residues in Prp8 intein of Aspergillus fumigatus: a potential target for antimycotics

Resistance to azoles and amphotericin B especially in Aspergillus fumigatus is a growing concern towards the treatment of invasive fungal infection. At this critical juncture, intein splicing would be a productive, and innovative target to establish therapies against resistant strains. Intein splicing is the central event for the activation of host protein, essential for the growth and survival of various microorganisms including A. fumigatus. The splicing process is a four-step protease-like nucleophilic cascade. Thus, we hypothesise that protease inhibitors would successfully halt intein splicing and potentially restrict the growth of the aforementioned pathogen. Using Rosetta Fold and molecular dynamics simulations, we modelled Prp8 intein structure; resembling classic intein fold with horse shoe shaped splicing domain. To fully comprehend the active site of Afu Prp8 intein, C1, T62, H65, H818, N819 from intein sequences and S820, the first C-extein residue are selected. Molecular docking shows that two FDA-approved drugs, i.e. Lufotrelvir and Remdesivir triphosphate efficiently interact with Prp8 intein from the assortment of 212 protease inhibitors. MD simulation portrayed that Prp8 undergoes conformational change upon ligand binding, and inferred the molecular recognition and stability of the docked complexes. Per-residue decomposition analysis confirms the importance of F: block R802, V803, and Q807 binding pocket in intein splicing domain towards recognition of inhibitors, along with active site residues through strong hydrogen bonds and hydrophobic contacts. However, in vitro and in vivo assays are required to confirm the inhibitory action on Prp8 intein splicing; which may pave the way for the development of new antifungals for A. fumigatus.Communicated by Ramaswamy H. Sarma.

Extein residues regulate the catalytic function of Spl DnaX intein enzyme by restricting the near-attack conformations of the active-site residues

Intein enzymes catalyze the splicing of their flanking polypeptide chains and have found tremendous biotechnological applications. Their terminal residues form the catalytic core and participate in the splicing reaction. Hence, the neighboring N- and C-terminal extein residues influence the catalytic rate. As these extein residues vary depending on the substrate identity, we tested the influence of 20 amino acids at these sites in the Spl DnaX intein and observed significant variation of spliced product as well as N- and C-terminus cleavage product formation. We investigated the dependence of these reactions on the extein residues by molecular dynamics (MD) simulations on eight extein variants, and found that the conformational sampling of the active-site residues of the intein enzyme differed among these extein variants. We found that the extein variants that sample higher population of near-attack conformers (NACs) of the active-site residues undergo higher product formation in our activity assays. Ground state conformers that closely resemble the transition state are referred to as NACs. Very good correlation was observed between the NAC populations from the MD simulations of eight extein variants and the corresponding product formation from our activity assays. Furthermore, this molecular detail enabled us to elucidate the mechanistic roles of several conserved active-site residues in the splicing reaction. Overall, this study shows that the catalytic power of Spl DnaX intein enzyme, and most likely other inteins, depends on the efficiency of formation of NACs in the ground state, which is further modulated by the extein residues.

Mechanism of covalent binding of ibrutinib to Bruton's tyrosine kinase revealed by QM/MM calculations

Ibrutinib is the first covalent inhibitor of Bruton's tyrosine kinase (BTK) to be used in the treatment of B-cell cancers. Understanding the mechanism of covalent inhibition will aid in the design of safer and more selective covalent inhibitors that target BTK. The mechanism of covalent inhibition in BTK has been uncertain because there is no appropriate residue nearby that can act as a base to deprotonate the cysteine thiol prior to covalent bond formation. We investigate several mechanisms of covalent modification of C481 in BTK by ibrutinib using combined quantum mechanics/molecular mechanics (QM/MM) molecular dynamics reaction simulations. The lowest energy pathway involves direct proton transfer from C481 to the acrylamide warhead in ibrutinib, followed by covalent bond formation to form an enol intermediate. There is a subsequent rate-limiting keto-enol tautomerisation step (ΔG ^‡ = 10.5 kcal mol^-1) to reach the inactivated BTK/ibrutinib complex. Our results represent the first mechanistic study of BTK inactivation by ibrutinib to consider multiple mechanistic pathways. These findings should aid in the design of covalent drugs that target BTK and other similar targets.

A novel artificial intelligence protocol for finding potential inhibitors of acute myeloid leukemia

There is currently no effective treatment for acute myeloid leukemia, and surgery is also ineffective as an important treatment for most tumors. Rapidly developing artificial intelligence technology can be applied to different aspects of drug development, and it plays a key role in drug discovery. Based on network pharmacology and virtual screening, candidates were selected from the molecular database. Nine artificial intelligence algorithm models were used to further verify the candidates' potential. The 350 training results of the deep learning model showed higher credibility, and the R-square of the training set and test set of the optimal model reached 0.89 and 0.84, respectively. The random forest model has an R-square of 0.91 and a mean square error of only 0.003. The R-square of the Adaptive Boosting model and the Bagging model reached 0.92 and 0.88, respectively. Molecular dynamics simulation evaluated the stability of the ligand-protein complex and achieved good results. Artificial intelligence models had unearthed the promising candidates for STAT3 inhibitors, and the good performance of most models showed that they still had practical value on small data sets.

Combined Quantum Mechanics and Molecular Mechanics Studies of Enzymatic Reaction Mechanisms

The combined quantum mechanics/molecular mechanics (QM/MM) methods have become a valuable tool in computational biochemistry and received versatile applications for studying the reaction mechanisms of enzymes. The approach combines the calculations of the electronic structure of the active site by QM, with modeling of the protein environment using MM force field, which allows the long-range electrostatics and steric effects on the enzyme reactivity to be accounted for. In this review, we review some key theoretical and computational aspects of the method and we also present some applications to particular enzymatic reactions such as tryptophan-7-halogenase, cyclooxygenase-1, and the epidermal growth factor receptor.

Coordination of the third step of protein splicing in two cyanobacterial inteins

The third step of protein splicing is cyclization of Asn coupled to peptide bond cleavage. In two related cyanobacterial inteins, this step is facilitated by Asn or Gln. For a Synechococcus sp. PCC7002 intein, the isolated third step of protein splicing is more efficient with its native Asn than with substitution to Gln. For a Trichodesmium erythraeum intein, its native Gln facilitates the third step as efficiently as with Asn. Despite these differences, the yield of splicing is not affected, suggesting that the third step is influenced by mechanism-linked conformational changes. A conserved catalytic His and the penultimate residue also play roles in promoting side-chain cyclization.

Quantum mechanics/molecular mechanics modeling of covalent addition between EGFR-cysteine 797 and N-(4-anilinoquinazolin-6-yl) acrylamide

Irreversible epidermal growth factor receptor (EGFR) inhibitors can circumvent resistance to first-generation ATP-competitive inhibitors in the treatment of nonsmall-cell lung cancer. They covalently bind a noncatalytic cysteine (Cys797) at the surface of EGFR active site by an acrylamide warhead. Herein, we used a hybrid quantum mechanics/molecular mechanics (QM/MM) potential in combination with umbrella sampling in the path-collective variable space to investigate the mechanism of alkylation of Cys797 by the prototypical covalent inhibitor N-(4-anilinoquinazolin-6-yl) acrylamide. Calculations show that Cys797 reacts with the acrylamide group of the inhibitor through a direct addition mechanism, with Asp800 acting as a general base/general acid in distinct steps of the reaction. The obtained reaction free energy is negative (ΔA = -12 kcal/mol) consistent with the spontaneous and irreversible alkylation of Cys797 by N-(4-anilinoquinazolin-6-yl) acrylamide. Our calculations identify desolvation of Cys797 thiolate anion as a key step of the alkylation process, indicating that changes in the intrinsic reactivity of the acrylamide would have only a minor impact on the inhibitor potency.

Dynamics differentiate between active and inactive inteins

The balance between stability and dynamics for active enzymes can be somewhat quantified by studies of intein splicing and cleaving reactions. Inteins catalyze the ligation of flanking host exteins while excising themselves. The potential for applications led to engineering of a mini-intein splicing domain, where the homing endonuclease domain of the Mycobacterium tuberculosis RecA (Mtu recA) intein was removed. The remaining domains were linked by several short peptides, but splicing activity in all was substantially lower than the full-length intein. Native splicing activity was restored in some cases by a V67L mutation. Using computations and experiments, we examine the impact of this mutation on the stability and conformational dynamics of the mini-intein splicing domain. Molecular dynamics simulations were used to delineate the factors that determine the active state, including the V67L mini-intein mutant, and peptide linker. We found that (1) the V67L mutation lowers the global fluctuations in all modeled mini-inteins, stabilizing the mini-intein constructs; (2) the connecting linker length affects intein dynamics; and (3) the flexibilities of the linker and intein core are higher in the active structure. We have observed that the interaction of the linker region and a turn region around residues 35-41 provides the pathway for the allostery interaction. Our experiments reveal that intein catalysis is characterized by non-linear Arrhenius plot, confirming the significant contribution of protein conformational dynamics to intein function. We conclude that while the V67L mutation stabilizes the global structure, cooperative dynamics of all intein regions appear more important for intein function than high stability. Our studies suggest that effectively quenching the conformational dynamics of an intein through engineered allosteric interactions could deactivate intein splicing or cleaving.

A QM/MM study of the catalytic mechanism of nicotinamidase

Nicotinamidase (Pnc1) is a member of Zn-dependent amidohydrolases that hydrolyzes nicotinamide (NAM) to nicotinic acid (NA), which is a key step in the salvage pathway of NAD(+) biosynthesis. In this paper, the catalytic mechanism of Pnc1 has been investigated by using a combined quantum-mechanical/molecular-mechanical (QM/MM) approach based on the recently obtained crystal structure of Pnc1. The reaction pathway, the detail of each elementary step, the energetics of the whole catalytic cycle, and the roles of key residues and Zn-binding site are illuminated. Our calculation results indicate that the catalytic water molecule comes from the bulk solvent, which is then deprotonated by residue D8. D8 functions as a proton transfer station between C167 and NAM, while the activated C167 serves as the nucleophile. The residue K122 only plays a role in stabilizing intermediates and transition states. The oxyanion hole formed by the amide backbone nitrogen atoms of A163 and C167 has the function to stabilize the hydroxyl anion of nicotinamide. The Zn-binding site rather than a single Zn(2+) ion acts as a Lewis acid to influence the reaction. Two elementary steps, the activation of C167 in the deamination process and the decomposition of catalytic water in the hydrolysis process, correspond to the large energy barriers of 25.7 and 28.1 kcal mol(-1), respectively, meaning that both of them contribute a lot to the overall reaction barrier. Our results may provide useful information for the design of novel and efficient Pnc1 inhibitors and related biocatalytic applications.

Conserved residues that modulate protein trans-splicing of Npu DnaE split intein

The first crystal trans-structure of a naturally occurring split intein has been determined for the Npu (Nostoc punctiforme PCC73102) DnaE split intein. Guided by this structure, the residues NArg50 and CSer35, well conserved in DnaE split inteins, are identified to be critical in the trans-splicing of Npu DnaE split intein. An in vitro splicing assay demonstrates that NArg50 and CSer35 play synergistic roles in modulating its intein activity. The C-terminal CAsn36 exhibits two orientations of its side chain and interacts with both NArg50 and CSer35 through hydrogen bonding. These interactions likely facilitate the cyclization of asparagine in the course of protein splicing. The mutation of either residue reduces intein activity, and correlates with the low activity of the Ssp (Cyanobacterium synechocystis sp. strain PCC6803) DnaE split intein. On the other hand, NArg50 also forms a hydrogen bond with the highly conserved F-block CAsp17, thus influencing the N-S acyl shift during N-terminal cleavage. Sequence alignments show that residues NArg50 and CSer35 are rather conserved in those split inteins that lack a penultimate histidine residue. The conserved non-catalytic residues of split inteins modulate the efficiency of protein trans-splicing by hydrogen-bond interactions with the catalytic residues at the splice junction.

QM/MM free-energy simulations of reaction in Serratia marcescens Chitinase B reveal the protonation state of Asp142 and the critical role of Tyr214

Serratia marcescens Chitinase B (ChiB), belonging to the glycosidase family 18 (GH18), catalyzes the hydrolysis of β-1,4-glycosidic bond, with retention of configuration, via an unusual substrate-assisted mechanism, in which the substrate itself acts as an intramolecular nucleophile. Here, both elementary steps (glycosylation and deglycosylation) of the ChiB-catalyzed reaction are investigated by means of combined quantum mechanics/molecular mechanics (QM/MM) umbrella sampling molecular dynamics (MD) simulations at the SCC-DFTB/CHARMM22 level of theory. We examine the influence of the Asp142 protonation state on the reaction and the role that this residue performs in the reaction. Our simulations show that reaction with a neutral Asp142 is preferred and demonstrate that this residue provides electrostatic stabilization of the oxazolinium ion intermediate formed in the reaction. Insight into the conformational itinerary ((1,4)B↔(4)H5↔(4)C1) adopted by the substrate (bound in subsite -1) along the preferred reaction pathway is also provided by the simulations. The relative energies of the stationary points found along the reaction pathway calculated with SCC-DFTB and B3LYP were compared. The results suggest that SCC-DFTB is an accurate method for estimating the relative barriers for both steps of the reaction; however, it was found to overestimate the relative energy of an intermediate formed in the reaction when compared with the higher level of theory. Glycosylation is suggested to be a rate-determining step in the reaction with calculated overall reaction free-energy barrier of 20.5 kcal/mol, in a reasonable agreement with the 16.1 kcal/mol barrier derived from the experiment. The role of Tyr214 in catalysis was also investigated with the results, indicating that the residue plays a critical role in the deglycosylation step of the reaction. Simulations of the enzyme-product complex were also performed with an unbinding event suggested to have been observed, affording potential new mechanistic insight into the release of the product of ChiB.

Protein splicing: how inteins escape from precursor proteins

Inteins are nature's escape artists; they facilitate their excision from flanking polypeptides (exteins) concomitant with extein ligation to produce a mature host protein. Splicing requires sequential nucleophilic displacement reactions catalyzed by strategies similar to proteases and asparagine lyases. Inteins require precise reaction coordination rather than rapid turnover or tight substrate binding because they are single turnover enzymes with covalently linked substrates. This has allowed inteins to explore alternative mechanisms with different steps or to use different methods for activation and coordination of the steps. Pressing issues include understanding the underlying details of catalysis and how the splicing steps are controlled.

Intein based bioprocess for production of a synthetic antimicrobial peptide: an alternative route to solid phase peptide synthesis

Intein based bioprocessing strategy for producing antimicrobial peptide candidate was found to be more sustainable compared to solid phase peptide synthesis strategy (SPPS).

Modeling catalytic promiscuity in the alkaline phosphatase superfamily

In recent years, it has become increasingly clear that promiscuity plays a key role in the evolution of new enzyme function. This finding has helped to elucidate fundamental aspects of molecular evolution. While there has been extensive experimental work on enzyme promiscuity, computational modeling of the chemical details of such promiscuity has traditionally fallen behind the advances in experimental studies, not least due to the nearly prohibitive computational cost involved in examining multiple substrates with multiple potential mechanisms and binding modes in atomic detail with a reasonable degree of accuracy. However, recent advances in both computational methodologies and power have allowed us to reach a stage in the field where we can start to overcome this problem, and molecular simulations can now provide accurate and efficient descriptions of complex biological systems with substantially less computational cost. This has led to significant advances in our understanding of enzyme function and evolution in a broader sense. Here, we will discuss currently available computational approaches that can allow us to probe the underlying molecular basis for enzyme specificity and selectivity, discussing the inherent strengths and weaknesses of each approach. As a case study, we will discuss recent computational work on different members of the alkaline phosphatase superfamily (AP) using a range of different approaches, showing the complementary insights they have provided. We have selected this particular superfamily, as it poses a number of significant challenges for theory, ranging from the complexity of the actual reaction mechanisms involved to the reliable modeling of the catalytic metal centers, as well as the very large system sizes. We will demonstrate that, through current advances in methodologies, computational tools can provide significant insight into the molecular basis for catalytic promiscuity, and, therefore, in turn, the mechanisms of protein functional evolution.

Quantum mechanics/molecular mechanics modeling of fatty acid amide hydrolase reactivation distinguishes substrate from irreversible covalent inhibitors

Carbamate and urea derivatives are important classes of fatty acid amide hydrolase (FAAH) inhibitors that carbamoylate the active-site nucleophile Ser241. In the present work, the reactivation mechanism of carbamoylated FAAH is investigated by means of a quantum mechanics/molecular mechanics (QM/MM) approach. The potential energy surfaces for decarbamoylation of FAAH covalent adducts, derived from the O-aryl carbamate URB597 and from the N-piperazinylurea JNJ1661610, were calculated and compared to that for deacylation of FAAH acylated by the substrate oleamide. Calculations show that a carbamic group bound to Ser241 prevents efficient stabilization of transition states of hydrolysis, leading to large increments in the activation barrier. Moreover, the energy barrier for the piperazine carboxylate was significantly lower than that for the cyclohexyl carbamate derived from URB597. This is consistent with experimental data showing slowly reversible FAAH inhibition for the N-piperazinylurea inhibitor and irreversible inhibition for URB597.

Modeling methods for studying post-translational and transcriptional modifying enzymes

Biological catalysis is a complex chemical process that involves not only electronic reorganization in the substrate but also the reorganization of the catalyst. This complexity is even larger in the case of post-transcriptional and post-translational modifications which may involve the interaction between two biomacromolecules. However, the development over the past decades of new computational methods and strategies is offering a detailed molecular picture of the catalytic event and a deep understanding of the mechanisms of chemical reactions in biological environments. Here we review the efforts made in the last years to model catalysis in post-transcriptional and post-translational processes. We stress on the advantages and problems of the different computational strategies and their applicability in different cases.