The stability and dynamics of computationally designed proteins

Modulation of Krüppel-like factors (KLFs) interaction with their binding partners in cancers through acetylation and phosphorylation

Post-translational modifications (PTMs) of transcription factors regulate transcriptional activity and play a key role in essentially all biological processes and generate indispensable insight towards biological function including activity state, subcellular localization, protein solubility, protein folding, substrate trafficking, and protein-protein interactions. Amino acids modified chemically via PTMs, function as molecular switches and affect the protein function and characterization and increase the proteome complexity. Krüppel-like transcription factors (KLFs) control essential cellular processes including proliferation, differentiation, migration, programmed cell death and various cancer-relevant processes. We investigated the interactions of KLF group-2 members with their binding partners to assess the role of acetylation and phosphorylation in KLFs on their binding affinity. It was observed that acetylation and phosphorylation at different positions in KLFs have a variable effect on binding with specific partners. KLF2-EP300, KLF4-SP1, KLF6-ATF3, KLF6-JUN, and KLF7-JUN show stabilization upon acetylation or phosphorylation at variable positions. On the other hand, KLF4-CBP, KLF4-EP300, KLF5-CBP, KLF5-WWP1, KLF6-SP1, and KLF7-ATF3 show stabilization or destabilization due to acetylation or phosphorylation at variable positions in KLFs. This provides a molecular explanation of the experimentally observed dual role of KLF group-2 members as a suppressor or activator of cancers in a PTM-dependent manner.

Computer-aided design to enhance the stability of aldo-keto reductase KdAKR

The aldo-keto reductase (AKR) KdAKR from Kluyvermyces dobzhanskii can reduce t-butyl 6-chloro-(5S)-hydroxy-3-oxohexanoate ((5S)-CHOH) to t-butyl 6-chloro-(3R,5S)-dihydroxyhexanoate ((3R,5S)-CDHH), which is the key chiral intermediate of rosuvastatin. Herein, a computer-aided design that combined the use of PROSS platform and consensus design was employed to improve the stability of a previously constructed mutant KdAKRM6 . Experimental verification revealed that S196C, T232A, V264I and V45L produced improved thermostability and activity. The "best" mutant KdAKRM10 (KdAKRM6 -S196C/T232A/V264I/V45L) was constructed by combining the four beneficial mutations, which displayed enhanced thermostability. Its T50 15 and Tm values were increased by 10.2 and 10.0°C, respectively, and half-life (t1/2 ) at 40°C was increased by 17.6 h. Additionally, KdAKRM10 demonstrated improved resistance to organic solvents compared to that of KdAKRM6 . Structural analysis revealed that the increased number of hydrogen bonds and stabilized hydrophobic core contributed to the rigidity of KdAKRM10 , thus improving its stability. The results validated the feasibility of the computer-aided design strategy in improving the stability of AKRs.

The Potential Mechanism of Liujunzi Decoction in the Treatment of Breast Cancer based on Network Pharmacology and Molecular Docking Technology

BACKGROUND

Liujunzi Decoction (LJZD) is a potential clinical treatment for Breast Cancer (BC), but the active ingredients and mechanisms underlying its effectiveness remain unclear.

OBJECTIVE

The study aimed to investigate the target gene of LJZD compatibility and the possible mechanism of action in the treatment of breast cancer by using network pharmacology and molecular docking.

METHODS

Based on TCMSP, ETCM, and BATMAN database searching and screening to obtain the ingredients of LJZD, the related targets were obtained. Breast cancer-related targets were collected through GEO, Geencards, OMIM, and other databases, and drug-disease Venn diagrams were drawn by R. The PPI network map was constructed by using Cytoscape. The intersecting targets were imported into the STRING database, and the core targets were analyzed and screened. The intersected targets were analyzed by the DAVID database for GO and KEGG enrichment. AutoDock Vina and Gromacs were used for molecular docking and simulation of the core targets and active ingredients.

RESULTS

126 active ingredients of LJZD were obtained; 241 targets related to breast cancer were sought after screening, and 180 intersection targets were identified through Venn diagram analysis. The core targets were FOS and ESR1. KEGG enrichment analysis mainly involved PI3K/Akt, MAPK, and other signaling pathways.

CONCLUSION

This study has explored the possible targets and signaling pathways of LJZD in treating breast cancer through network pharmacology and bioinformatics analysis. Molecular docking and simulation have further validated the potential mechanism of action of LJZD in breast cancer treatment, providing essential experimental data for future studies.

Characterizing the combined effects of cytochrome P450 missense variation within star allele definitions

Background: Cytochrome P450 (CYP) genetic variation largely impacts drug response. However, many CYP star alleles (haplotypes) lack functional annotation, impeding our understanding of drug metabolism mechanisms. We aimed to investigate the impact of missense variant combinations on CYP protein structures. Methods: Normal mode analysis was conducted on 261 missense variants within 91 CYP haplotypes. CYP2D6*2 and CYP2D6*17 were prioritized for molecular dynamics simulation. Results: Normal mode analysis and molecular dynamics highlight the effects of known CYP missense variants on protein stability and conformational dynamics. Missense variants within haplotypes may have intermodulating effects on protein structure and function. Conclusion: This study highlights the utility of multiscale modeling in interpreting CYP missense variants and particularly their combinations within various star alleles.

Simultaneous enhancement of thermostability and catalytic activity of κ-carrageenase from Pseudoalteromonas tetraodonis by rational design

κ-Carrageenase provides an attractive enzymatic approach to preparation of κ-carrageenan oligosaccharides. Pseudoalteromonas tetraodonis κ-carrageenase is active at the alkaline conditions but displays low thermostability. To further improve its enzymatic performance, two mutants of Q42V and I51H exhibiting both improved thermostability and enzyme activity were screened by the PoPMuSiC algorithm. Compared with the wild-type κ-carrageenase (WT), Q42V and I51H increased the enzyme activity by 20.9% and 25.4%, respectively. After treatment at 50 ℃ for 40 min, Q42V and I51H enhanced the residual activity by 31.1% and 25.9%, respectively. The T_m values of Q42V, I51H, and WT determined by differential scanning calorimetry were 58.2 ℃, 54.8 ℃, and 51.2 ℃, respectively. Compared with untreated and HCl-treated κ-carrageenans, Q42V-treated κ-carrageenan exhibited higher pancreatic lipase inhibitory activity. Molecular docking analysis indicated that the additional pi-sigma force and hydrophobic interaction in the enzyme-substrate complex could account for the increased catalytic activity of Q42V and I51H, respectively. Molecular dynamics analysis indicated that the improved thermostability of mutants Q42V and I51H could be attributed to the less structural deviation and the flexible changes of enzyme conformation at high temperature. This study provides new insight into κ-carrageenase performance improvement and identifies good candidates for their industrial applications.

Using High-Throughput Molecular Dynamics Simulation to Enhance the Computational Design of Kemp Elimination Enzymes

Computational enzyme design has been successfully applied to identify new alternatives to natural enzymes for the biosynthesis of important compounds. However, the moderate catalytic activities of de novo designed enzymes indicate that the modeling accuracy of current computational enzyme design methods should be improved. Here, high-throughput molecular dynamics simulations were used to enhance computational enzyme design, thus allowing the identification of variants with higher activities in silico. Different time schemes of high-throughput molecular dynamics simulations were tested to identify the catalytic features of evolved Kemp eliminases. The 20 × 1 ns molecular dynamics simulation scheme was sufficiently accurate and computationally viable to screen the computationally designed massive variants of Kemp elimination enzymes. The developed hybrid computational strategy was used to redesign the most active Kemp eliminase, HG3.17, and five variants were generated and experimentally confirmed to afford higher catalytic efficiencies than that of HG3.17, with one double variant (D52Q/A53S) exhibiting a 55% increase. The hybrid computational enzyme design strategy is general and computationally economical, with which we anticipate the efficient creation of practical enzymes for industrial biocatalysis.

Enhanced Thermostability and Molecular Insights for l-Asparaginase from Bacillus licheniformis via Structure- and Computation-Based Rational Design

l-Asparaginase has gained much attention for effectively treating acute lymphoblastic leukemia (ALL) and mitigating carcinogenic acrylamide in fried foods. Due to high-dose dependence for clinical treatment and low mitigation efficiency for thermal food processes caused by poor thermal stability, a method to achieve thermostable l-asparaginase has become a critical bottleneck. In this study, a rational design including free energy combined with structural and conservative analyses was applied to engineer the thermostability of l-asparaginase from Bacillus licheniformis (BlAsnase). Two enhanced thermostability mutants D172W and E207A were screened out by site-directed saturation mutagenesis. The double mutant D172W/E207A exhibited highly remarkable thermostability with a 65.8-fold longer half-life at 55 °C and 5 °C higher optimum reaction temperature and melting temperature (T_m) than those of wild-type BlAsnase. Further, secondary structure, sequence, molecular dynamics (MD), and 3D-structure analysis revealed that the excellent thermostability of the mutant D172W/E207A was on account of increased hydrophobicity and decreased flexibility, highly rigid structure, hydrophobic interactions, and favorable electrostatic potential. As the first report of rationally designing l-asparaginase with improved thermostability from B. licheniformis, this study offers a facile and efficient process to improve the thermostability of l-asparaginase for industrial applications.

Engineering and Design of Programmable Genome Editors

Programmable genome editors are enzymes that can be targeted to a specific location in the genome for making site-specific alterations or deletions. The engineering, design, and development of sequence-specific editors has resulted in a dramatic increase in the precision of editing for nucleotide sequences. These editors can target specific locations in a genome, in vivo. The genome editors are being deployed for the development of genetically modified organisms for agriculture and industry, and for gene therapy of inherited human genetic disorders, cancer, and immunotherapy. Experimental and computational studies of structure, binding, activity, dynamics, and folding, reviewed here, have provided valuable insights that have the potential for increasing the functional efficiency of these gene/genome editors. Biochemical and biophysical studies of the specificities of natural and engineered genome editors reveal that increased binding affinity can be detrimental because of the increase of off-target effects and that the engineering and design of genome editors with higher specificity may require modulation and control of the conformational dynamics.

Thermostability engineering of an inulin fructotransferase for the biosynthesis of difructose anhydride I

The thermostability of enzymes is an essential factor that performs a vital role during practical applications. Inulin fructotransferases can efficiently convert inulin into bio-functional difructose anhydrides (DFAs). The present study aimed to improve the thermostability of a previously reported inulin fructotransferase, SpIFTase, and apply it to the biosynthesis of DFA I. In silico rational design was used to predict mutation sites, based on sequential and structural information. Two triple-site mutants, Q69L/Q234L/K310G and E201I/Q234L/K310G, were characterized and exhibited enhanced thermostability with approximately 5 °C higher in melting temperature (T_m), respectively, and a 45-fold longer half-life (t_1/2) at 70 °C, compared to that of SpIFTase. Molecular dynamic simulations and elaborate structural analysis suggested that the combinations of hydrophobic interaction, electrostatic potential distribution, and decreased flexibility via stabilization of loops and α-helix improved the thermostability of SpIFTase. Additionally, the promising mutants exhibited great potential to the industrial production of DFA I.

Integrating dynamics into enzyme engineering

Enzyme engineering has become a widely adopted practice in research labs and industry. In parallel, the past decades have seen tremendous strides in characterizing the dynamics of proteins, using a growing array of methodologies. Importantly, links have been established between the dynamics of proteins and their function. Characterizing the dynamics of an enzyme prior to, and following, its engineering is beginning to inform on the potential of 'dynamic engineering', i.e. the rational modification of protein dynamics to alter enzyme function. Here we examine the state of knowledge at the intersection of enzyme engineering and protein dynamics, describe current challenges and highlight pioneering work in the nascent area of dynamic engineering.