Insight into the structural stability of wild type and mutants of the tobacco etch virus protease with molecular dynamics simulations

Mechanism of Mutation-Induced Effects on the Catalytic Function of TEV Protease: A Molecular Dynamics Study

Tobacco etch virus protease (TEVp) is wildly exploited for various biotechnological applications. These applications take advantage of TEVp's ability to cleave specific substrate sequences to study protein function and interactions. A major limitation of this enzyme is its relatively slow catalytic rate. In this study, MD simulations were conducted on TEV enzymes and known highly active mutants (eTEV and uTEV3) to explore the relationship between mutation, conformation, and catalytic function. The results suggest that mutations distant from the active site can influence the substrate-binding pocket through interaction networks. MD analysis of eTEV demonstrates that, by stabilizing the orientation of the substrate at the catalytic site, mutations that appropriately enlarge the substrate-binding pocket will be beneficial for Kcat, enhancing the catalytic efficiency of the enzyme. On the contrary, mutations in uTEV3 reduced the flexibility of the active pocket and increased the hydrogen bonding between the substrate and enzyme, resulting in higher affinity. At the same time, the MD simulation demonstrates that mutations outside of the active site residues could affect the dynamic movement of the binding pocket by altering residue networks and communication pathways, thereby having a profound impact on reactivity. These findings not only provide a molecular mechanistic explanation for the excellent mutants, but also serve as a guiding framework for rational computational design.

Alchemical Binding Free Energy Calculations in AMBER20: Advances and Best Practices for Drug Discovery

Predicting protein-ligand binding affinities and the associated thermodynamics of biomolecular recognition is a primary objective of structure-based drug design. Alchemical free energy simulations offer a highly accurate and computationally efficient route to achieving this goal. While the AMBER molecular dynamics package has successfully been used for alchemical free energy simulations in academic research groups for decades, widespread impact in industrial drug discovery settings has been minimal because of the previous limitations within the AMBER alchemical code, coupled with challenges in system setup and postprocessing workflows. Through a close academia-industry collaboration we have addressed many of the previous limitations with an aim to improve accuracy, efficiency, and robustness of alchemical binding free energy simulations in industrial drug discovery applications. Here, we highlight some of the recent advances in AMBER20 with a focus on alchemical binding free energy (BFE) calculations, which are less computationally intensive than alternative binding free energy methods where full binding/unbinding paths are explored. In addition to scientific and technical advances in AMBER20, we also describe the essential practical aspects associated with running relative alchemical BFE calculations, along with recommendations for best practices, highlighting the importance not only of the alchemical simulation code but also the auxiliary functionalities and expertise required to obtain accurate and reliable results. This work is intended to provide a contemporary overview of the scientific, technical, and practical issues associated with running relative BFE simulations in AMBER20, with a focus on real-world drug discovery applications.

Rational design of a new mutant of tobacco etch virus protease in order to increase the in vitro solubility

Background and purpose:

Tobacco etch virus (TEV) protease is a protease with high sequence specificity which is useful for the cleavage of fusion proteins. A major limitation of this enzyme is its relatively poor solubility. This study aimed to investigate the effects of some suggested mutations by online tools and molecular dynamics simulation to improve the solubility of TEV protease in vitro.

Experimental approach:

We designed a rational multi-stage process to determine the solubilizing mutations of TEV protease. At the first stage, all the possible mutations were predicted using online tools such as PoPMuSiC and Eris servers, in which five mutations include N23F, N23L, Q74L, Q74V, and Q74I were suggested for further studies. In the next step, the three dimensional structure of the wild type (WT) and the best mutations were subjected to molecular dynamic simulations to evaluate the dynamic behaviour of the obtained structures. The selected mutation was introduced into the structure using site-directed mutagenesis and expressed in Escherichia coli BL21DE3. After purification, solubility and activity of the purified mutant and WT-TEV proteases were assayed.

Findings /Results:

By considering the analysis of various factors such as structural and solubility properties, one mutant, N23F, was selected for in vitro studies which led to a 1.5 times increase in the solubility compared to the WT while its activity was decreased somewhat.

Conclusion and implications:

We propose N23F mutation, according to computational and experimental analyses for TEV proteases which resulted in a 150% increase in solubility compared to the WT.

Tobacco etch virus (TEV) protease with multiple mutations to improve solubility and reduce self-cleavage exhibits enhanced enzymatic activity

Tobacco etch virus (TEV) protease is a 27‐kDa catalytic domain of the polyprotein nuclear inclusion a (NIa) in TEV, which recognizes the specific amino acid sequence ENLYFQG/S and cleaves between Q and G/S. Despite its substrate specificity, its use is limited by its autoinactivation through self‐cleavage and poor solubility during purification. It was previously reported that T17S/N68D/I77V mutations improve the solubility and yield of TEV protease and S219 mutations provide protection against self‐cleavage. In this study, we isolated TEV proteases with S219N and S219V mutations in the background of T17S, N68D, and I77V without the inclusion body, and measured their enzyme kinetics. The k_cat of two isolated S219N and S219V mutants in the background of T17S, N68D, and I77V mutations was highly increased compared to that of the control, and S219N was twofold faster than S219V without K_m change. This result indicates that combination of these mutations can further enhance TEV activity.

Insight into the structural stability of wild-type and histidine mutants in Pin1 by experimental and computational methods

Pin1, a polypeptide proline isomerase parvulin, plays a key role in Alzheimer's disease (AD), common tumors and cancers. Two conservative histidine residues, His59 and His157, are important for maintaining the stability of the PPIase domain. Hence multiple spectral and computational techniques were performed to investigate the potential mechanism of two histidine residues. Thermal denaturation indicated that both residues His59 and His157 are not sensitive to the lower temperatures, while residue His59 is more sensitive to the higher temperatures than residue His157. Acidic denaturation suggested that influences of both residues His59 and His157 to acidic stability were the difference from Pin1-WT. ANS and RLS spectra hinted that there was no significant effect on hydrophobic change and aggregation by histidine mutations. The GndHCl-induced denaturation implied that residues His59 and His157 contributed the most to the chemical stability. MD simulations revealed that residues His59 and His157 mutations resulted in that the hydrogen bond network of the dual histidine motif was destroyed wholly. In summary, these histidine residues play an important role in maintaining the structural stability of the PPIase domain.

Molecular dynamics simulations provide insights into the origin of gleevec's selectivity toward human tyrosine kinases

Protein kinases are critical drug targets against cancer. Since the discovery of Gleevec, a specific inhibitor of Abl kinase, the capability of this drug to distinguish between Abl and other tyrosine kinases, such as Src, has been intensely investigated but the origin of Gleevec's selectivity to Abl against Src is less studied. Here, we performed molecular dynamics (MD) simulations, dynamical cross-correlation matrices (DCCM), dynamical network analysis, and binding free energy calculations to explore Gleevec's selectivity based on the crystal structures of Abl, Src, and their common ancestors (ANC-AS) and the two constructed mutation systems (AS→Abl and AS→Src). MD simulations revealed that the conformation of the phosphate-binding loop (P-loop) was altered significantly in the AS→Abl system. DCCM results unraveled that mutations increased anticorrelated motions in the AS→Abl system. Community network analysis suggested that the P-loop established special contacts in the AS→Abl system that are devoid in the AS→Src system. The binding free energy calculations unveiled that the affinity of Gleevec to AS→Abl increased to near the Abl level, whereas its affinity to AS→Src decreased to near the Src level. Analysis of individual residue contributions showed that the differences were located mainly at the P-loop. This study is valuable for understanding the sensitivity of Gleevec to human tyrosine kinases. Communicated by Ramaswamy H. Sarma.

Tobacco Etch Virus protease: A shortcut across biotechnologies

About thirty years ago, studies on the RNA genome of Tobacco Etch Virus revealed the presence of an efficient and specific protease, called Tobacco Etch Virus protease (TEVp), that was part of the Nuclear Inclusion a (NIa) enzyme. TEVp is an efficient and specific protease of 27kDa that has become a valuable biotechnological tool. Nowadays TEVp is a unique endopeptidase largely exploited in biotechnology from industrial applications to in vitro and in vivo cellular studies. A number of TEVp mutants with different rate of cleavage, stability and specificity have been reported. Similarly, a panel of different target cleavage sites, derived from the canonical ENLYFQ-G/S site, has been established. In this review we describe these aspects of TEVp and some of its multiple applications. A particular focus is on the use and molecular biology of TEVp in living cells and organisms.

Structural basis of valmerins as dual inhibitors of GSK3β/CDK5

Development of multi-target drugs is becoming increasingly attractive in the repertoire of protein kinase inhibitors discovery. In this study, we carried out molecular docking, molecular dynamics simulations, molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) binding free energy calculations, principal component analysis (PCA), and dynamical cross-correlation matrices (DCCM) to dissect the molecular mechanism for the valmerin-19 acting as a dual inhibitor for glycogen synthase kinase 3β (GSK3β) and cyclin-dependent kinase 5 (CDK5). Detailed MM-PBSA calculations revealed that the binding free energies of the valmerin-19 to GSK3β/CDK5 were calculated to be -12.60 ± 2.28 kcal mol(-1) and -11.85 ± 2.54 kcal mol(-1), respectively, indicating that valmerin-19 has the potential to act as a dual inhibitor of GSK3β/CDK5. The analyses of PCA and DCCM results unraveled that binding of the valmerin-19 reduced the conformational dynamics of GSK3β/CDK5 and the valmerin-19 bound to GSK3β/CDK5 might occur mostly through a conformational selection mechanism. This study may be helpful for the future design of novel and potent dual GSK3β/CDK5 inhibitors.