Voronoia 4-ever

MutationExplorer: a webserver for mutation of proteins and 3D visualization of energetic impacts

The possible effects of mutations on stability and function of a protein can only be understood in the context of protein 3D structure. The MutationExplorer webserver maps sequence changes onto protein structures and allows users to study variation by inputting sequence changes. As the user enters variants, the 3D model evolves, and estimated changes in energy are highlighted. In addition to a basic per-residue input format, MutationExplorer can also upload an entire replacement sequence. Previously the purview of desktop applications, such an upload can back-mutate PDB structures to wildtype sequence in a single step. Another supported variation source is human single nucelotide polymorphisms (SNPs), genomic coordinates input in VCF format. Structures are flexibly colorable, not only by energetic differences, but also by hydrophobicity, sequence conservation, or other biochemical profiling. Coloring by interface score reveals mutation impacts on binding surfaces. MutationExplorer strives for efficiency in user experience. For example, we have prepared 45 000 PDB depositions for instant retrieval and initial display. All modeling steps are performed by Rosetta. Visualizations leverage MDsrv/Mol*. MutationExplorer is available at: http://proteinformatics.org/mutation_explorer/.

Hydrophobic core evolution of major histocompatibility complex class I chain-related protein A for dramatic enhancing binding affinity

Interface residues at sites of protein-protein interaction (PPI) are the focus for affinity optimisation. However, protein hydrophobic cores (HCs) play critical roles and shape the protein surface. We hypothesise that manipulating protein HCs can enhance PPI interaction affinities. A cell stress molecule, major histocompatibility complex class I chain-related protein A (MICA), binds to the natural killer group 2D (NKG2D) homodimer to form three molecule interactions. MICA was used as a study subject to support our hypothesis. We redesigned MICA HCs by directed mutagenesis and isolated high-affinity variants through a newly designed partial-denature panning (PDP) method. A few mutations in MICA HCs increased the NKG2D-MICA interaction affinity by 325-5613-fold. Crystal structures of the NKG2D-MICA variant complexes indicated that mutagenesis of MICA HCs stabilised helical elements for decreasing intermolecular interactive free energy (ΔG) of the NKG2D-MICA heterotrimer. The repacking of MICA HC mutants maintained overall surface residues and the authentic binding specificity of MICA. In conclusion, this study provides a new method for MICA redesign and affinity optimisation through HC manipulation without mutating PPI interface residues. Our study introduces a novel approach to protein manipulation, potentially expanding the toolkit for protein affinity optimisation.

Structural Landscape of nsp Coding Genomic Regions of SARS-CoV-2-ssRNA Genome: A Structural Genomics Approach Toward Identification of Druggable Genome, Ligand-Binding Pockets, and Structure-Based Druggability

SARS-CoV-2 has a single-stranded RNA genome (+ssRNA), and synthesizes structural and non-structural proteins (nsps). All 16 nsp are synthesized from the ORF1a, and ORF1b regions associated with different life cycle preprocesses, including replication. The regions of ORF1a synthesizes nsp1 to 11, and ORF1b synthesizes nsp12 to 16. In this paper, we have predicted the secondary structure conformations, entropy & mountain plots, RNA secondary structure in a linear fashion, and 3D structure of nsp coding genes of the SARS-CoV-2 genome. We have also analyzed the A, T, G, C, A+T, and G+C contents, GC-profiling of these genes, showing the range of the GC content from 34.23 to 48.52%. We have observed that the GC-profile value of the nsp coding genomic regions was less (about 0.375) compared to the whole genome (about 0.38). Additionally, druggable pockets were identified from the secondary structure-guided 3D structural conformations. For secondary structure generation of all the nsp coding genes (nsp 1-16), we used a recent algorithm-based tool (deep learning-based) along with the conventional algorithms (centroid and MFE-based) to develop secondary structural conformations, and we found stem-loop, multi-branch loop, pseudoknot, and the bulge structural components, etc. The 3D model shows bound and unbound forms, branched structures, duplex structures, three-way junctions, four-way junctions, etc. Finally, we identified binding pockets of nsp coding genes which will help as a fundamental resource for future researchers to develop RNA-targeted therapeutics using the druggable genome.

Native architecture of a human GBP1 defense complex for cell-autonomous immunity to infection

All living organisms deploy cell-autonomous defenses to combat infection. In plants and animals, large supramolecular complexes often activate immune proteins for protection. In this work, we resolved the native structure of a massive host-defense complex that polymerizes 30,000 guanylate-binding proteins (GBPs) over the surface of gram-negative bacteria inside human cells. Construction of this giant nanomachine took several minutes and remained stable for hours, required guanosine triphosphate hydrolysis, and recruited four GBPs plus caspase-4 and Gasdermin D as a cytokine and cell death immune signaling platform. Cryo-electron tomography suggests that GBP1 can adopt an extended conformation for bacterial membrane insertion to establish this platform, triggering lipopolysaccharide release that activated coassembled caspase-4. Our "open conformer" model provides a dynamic view into how the human GBP1 defense complex mobilizes innate immunity to infection.

MUTATIONEXPLORER- A WEBSERVER FOR MUTATION OF PROTEINS AND 3D VISUALIZATION OF ENERGETIC IMPACTS

The possible effects of mutations on stability and function of a protein can only be understood in the context of protein 3D structure. The MutationExplorer webserver maps sequence changes onto protein structures and allows users to study variation by inputting sequence changes. As the user enters variants, the 3D model evolves, and estimated changes in energy are highlighted. In addition to a basic per-residue input format, MutationExplorer can also upload an entire replacement sequence. Previously the purview of desktop applications, such an upload can back-mutate PDB structures to wildtype sequence in a single step. Another supported variation source is human single nucelotide polymorphisms (SNPs), genomic coordinates input in VCF format.

Computer-Aided Reconstruction and Application of Bacillus halodurans S7 Xylanase with Heat and Alkali Resistance

β-1,4-Endoxylanase is the most critical hydrolase for xylan degradation during lignocellulosic biomass utilization. However, its poor stability and activity in hot and alkaline environments hinder its widespread application. In this study, BhS7Xyl from Bacillus halodurans S7 was improved using a computer-aided design through isothermal compressibility (βT) perturbation engineering and by combining three thermostability prediction algorithms (ICPE-TPA). The best variant with remarkable improvement in specific activity, heat resistance (70 °C), and alkaline resistance (both pH 9.0 and 70 °C), R69F/E137M/E145L, exhibited a 4.9-fold increase by wild-type in specific activity (1368.6 U/mg), a 39.4-fold increase in temperature half-life (458.1 min), and a 57.6-fold increase in pH half-life (383.1 min). Furthermore, R69F/E137M/E145L was applied to the hydrolysis of agricultural waste (corncob and hardwood pulp) to efficiently obtain a higher yield of high-value xylooligosaccharides. Overall, the ICPE-TPA strategy has the potential to improve the functional performance of enzymes under extreme conditions for the high-value utilization of lignocellulosic biomass.

Improving the Activity and Stability of Serine Protease ThAPT3 by Alleviating Self-Cleavage and Its Application in Deproteinization of Shrimp Shells

The self-cleavage properties of proteases result in low activity and instability, which limit their industrial application. In this study, the serine protease ThAPT3 from Torrubiella hemipterigena was successfully expressed in Komagataella phaffii. We investigated the self-degradation mechanism of ThAPT3 and presented a rational strategy to alleviate self-cleavage. A major self-degradation site (Leu²³⁸-Met²³⁹) and a primary autolysis region were identified. The autolysis regions (loop₁₈, α₈-helix, and loop₁₉) were redesigned and optimized using loop transplantation, energy calculations, surface cavity optimization, and loop anchoring. A triple-superposition mutant, ThAPT3-M9 (M²³⁹GKDGAVAAGLC²⁵⁰ → M²³⁹TLNRTTAANAC²⁵⁰/A251E/A254Q/R259L/A267E/S280N), was obtained. Compared to the wild type, the autolysis of M9 was significantly alleviated, and its half-life at 60 °C was increased approximately 39-fold (from 1.6 to 62.4 min). The optimal temperature and specific activity of M9 increased by 5 °C (from 60 to 65 °C) and 62% (4985 vs 3078 U/mg), respectively. M9 showed significant advantages in shrimp shell deproteinization.

Inside Out Computational Redesign of Cavities for Improving Thermostability and Catalytic Activity of Rhizomucor Miehei Lipase

Cavities are created by hydrophobic interactions between residue side chain atoms during the folding of enzymes. Redesigning cavities can improve the thermostability and catalytic activity of the enzyme; however, the synergistic effect of cavities remains unclear. In this study, Rhizomucor miehei lipase (RML) was used as a model to explore volume fluctuation and spatial distribution changes of the internal cavities, which could reveal the roles of internal cavities in the thermostability and catalytic activity. We present an inside out cavity engineering (CE) strategy based on computational techniques to explore how changes in the volumes and spatial distribution of cavities affect the thermostability and catalytic activity of the enzyme. We obtained 12 single-point mutants, among which the melting temperatures (Tm) of 8 mutants showed an increase of more than 2°C. Sixteen multipoint mutations were further designed by spatial distribution rearrangement of internal cavities. The Tm of the most stable triple variant, with mutations including T21V (a change of T to V at position 21), S27A, and T198L (T21V/S27A/T198L), was elevated by 11.0°C, together with a 28.7-fold increase in the half-life at 65°C and a specific activity increase of 9.9-fold (up to 5,828 U mg-1), one of the highest lipase activities reported. The possible mechanism of decreased volumes and spatial rearrangement of the internal cavities improved the stability of the enzyme, optimizing the outer substrate tunnel to improve the catalytic efficiency. Overall, the inside out computational redesign of cavities method could help to deeply understand the effect of cavities on enzymatic stability and activity, which would be beneficial for protein engineering efforts to optimize natural enzymes. IMPORTANCE In the present study, R. miehei lipase, which is widely used in various industries, provides an opportunity to explore the effects of internal cavities on the thermostability and catalytic activity of enzymes. Here, we execute high hydrostatic pressure molecular dynamics (HP-MD) simulations to screen the critical internal cavity and reshape the internal cavities through site-directed mutation. We show that as the global internal cavity volume decreases, cavity rearrangement can improve the stability of the protein while optimizing the substrate channel to improve the catalytic efficiency. Our results provide significant insights into understanding the mechanism of action of the internal cavity. Our strategy is expected to be applied to other enzymes to promote increases in thermostability and catalytic activity.