Multidimensional computational strategies enhance the thermostability of alpha-galactosidase

Alpha-Galactosidase has significant industrial application value in food processing, animal nutrition and medical applications. Microbial-derived α-galactosidases predominate industrial implementation due to high productivity, yet their inherent thermal instability necessitates systematic protein engineering. In this study, we established a dual-strategy protein engineering framework to enhance the thermostability of Aspergillus tubingensis α-galactosidase (AtWU_04653). Strategy I employed integrative computational design tools (ABACUS2/PROSS/DBD2) for mutational library construction, which yielded the dominant mutant A169P exhibiting remarkable performance: 78.52 % enhancement in thermal half-life at 55 °C (pH 4.0) and 52.04 % increase in catalytic efficiency (kcat /Km). Strategy II implemented a physics-based computational methodology combining GROMACS molecular dynamics simulations with Rosetta unfolding free energy calculations and SPIRED machine learning predictions, successfully deriving three stabilized variants (E429I, N380L, T64P) displaying 57.33 %, 67.17 %, and 41.34 % extended half-lives respectively. Notably, E429I and T64P demonstrated concurrent 85.25 % and 65.90 % catalytic activity augmentation (kcat /Km). Both strategies achieved substantial reduction in experimental screening workload while enabling synergistic thermostability-activity optimization. This study uses sequence conservation analysis, unfolding free energy calculation, molecular dynamics simulation, and innovative protein prediction models to establish multidimensional computational strategies for designing mutants, providing new and important technical references for computational design and functional optimization of enzymes.

Analytical control strategy for biologics. Part I: Foundations

Biologic therapeutics encompass different modalities with vastly different molecular profiles. Despite these differences, all products follow a similar approach to Pharmaceutical Development, which includes an integrated control strategy that relies on a clinical target product profile (TPP), a quality target product profile (QTPP), biophysical, biochemical and biological characterization, elucidation of critical quality attributes (CQAs), and development of an analytical control strategy. Technical and regulatory requirements for biologics development are established in numerous regulatory guidance documents issued by ICH, FDA, EMA, and other bodies. While there is substantial published knowledge on specific studies needed for development of a product, there is no specific guidance on establishing a comprehensive analytical control strategy as part of a modern integrated control strategy. This commentary is Part I of a two-part commentary series on analytical control strategy. In this part we present the foundations that are essential for developing an analytical control strategy to enable efficient lifecycle management across different biologic protein-based therapeutic modalities. In Part II, we will present a stage-appropriate roadmap to implementing an analytical control strategy from discovery research through the commercial life of the biologic.

Structure-based self-supervised learning enables ultrafast protein stability prediction upon mutation

Predicting free energy changes (ΔΔG) is essential for enhancing our understanding of protein evolution and plays a pivotal role in protein engineering and pharmaceutical development. While traditional methods offer valuable insights, they are often constrained by computational speed and reliance on biased training datasets. These constraints become particularly evident when aiming for accurate ΔΔG predictions across a diverse array of protein sequences. Herein, we introduce Pythia, a self-supervised graph neural network specifically designed for zero-shot ΔΔG predictions. Our comparative benchmarks demonstrate that Pythia outperforms other self-supervised pretraining models and force field-based approaches while also exhibiting competitive performance with fully supervised models. Notably, Pythia shows strong correlations and achieves a remarkable increase in computational speed of up to 105-fold. We further validated Pythia's performance in predicting the thermostabilizing mutations of limonene epoxide hydrolase, leading to higher experimental success rates. This exceptional efficiency has enabled us to explore 26 million high-quality protein structures, marking a significant advancement in our ability to navigate the protein sequence space and enhance our understanding of the relationships between protein genotype and phenotype. In addition, we established a web server at https://pythia.wulab.xyz to allow users to easily perform such predictions.

Biomimetic androgen receptor-based AIE biosensor for detecting bisphenol analogues: An integrating in silico topological analysis, molecular docking, and experimental validation study

Bisphenol analogues are the typical class of endocrine disrupting chemicals (EDCs) that interfere with binding of endogenous hormones to androgen receptor (AR). With the expansion of industrial activities and the intensification of environmental pollution, an increasing array of bisphenol analogues is being released into the environment and food chain. This highlights the urgency to develop sensitive methods for the detection of bisphenol analogues. Here, we propose a biomimetic AR-based biosensor platform for detecting bisphenol analogues (BPF, TBBPA, and TBBPS) by binding with Aggregation-Induced Emission (AIE) probes. Following a comparison of the PROSS and ABACUS methods, biomimetic AR was designed using the ABACUS approach and subsequently expressed in vitro via the E. coli expression system. Through molecular docking and the observation of fluorescence changes upon binding with biomimetic AR, BS-46006 was selected as the AIE probe for the biosensor. The biomimetic AR-based biosensor showed sensitive detections of BPF, TBBPA, and TBBPS within a range of 0-50 mM. To further elucidate the multi-residue recognition mechanism, molecular orbitals, Electron Localization Function (ELF), and Localized Orbital Locator (LOL) were systematically calculated in this study. Lowest unoccupied molecular orbital and highest occupied molecular orbital indicated the energy gap of BPF, TBBPA, and TBBPS, which correspond to 0.12812, 0.19689, and 0.18711 eV, respectively. ELF and LOL offered clearer perspective through heat maps to visually represent the electron delocalization in BPF, TBBPA, and TBBPS. The matrix effect analysis suggested that the responses of bisphenol analogues in soil matrices could be effectively mitigated through sample pretreatment. The analysis of spiked soil samples showed the acceptable recoveries ranged from 91 % to 105 %. Additionally, the biomimetic AR-based AIE biosensor, which combines multi-residue detection with Tolerable Daily Intakes, shows great promise for the risk assessment of bisphenol analogues. This research may present a viable approach for the analysis of environmental pollutants.

ProBID-Net: a deep learning model for protein-protein binding interface design

Protein-protein interactions are pivotal in numerous biological processes. The computational design of these interactions facilitates the creation of novel binding proteins, crucial for advancing biopharmaceutical products. With the evolution of artificial intelligence (AI), protein design tools have swiftly transitioned from scoring-function-based to AI-based models. However, many AI models for protein design are constrained by assuming complete unfamiliarity with the amino acid sequence of the input protein, a feature most suited for de novo design but posing challenges in designing protein-protein interactions when the receptor sequence is known. To bridge this gap in computational protein design, we introduce ProBID-Net. Trained using natural protein-protein complex structures and protein domain-domain interface structures, ProBID-Net can discern features from known target protein structures to design specific binding proteins based on their binding sites. In independent tests, ProBID-Net achieved interface sequence recovery rates of 52.7%, 43.9%, and 37.6%, surpassing or being on par with ProteinMPNN in binding protein design. Validated using AlphaFold-Multimer, the sequences designed by ProBID-Net demonstrated a close correspondence between the design target and the predicted structure. Moreover, the model's output can predict changes in binding affinity upon mutations in protein complexes, even in scenarios where no data on such mutations were provided during training (zero-shot prediction). In summary, the ProBID-Net model is poised to significantly advance the design of protein-protein interactions.

Microbial Cell Factories in the Bioeconomy Era: From Discovery to Creation

Microbial cell factories (MCFs) are extensively used to produce a wide array of bioproducts, such as bioenergy, biochemical, food, nutrients, and pharmaceuticals, and have been regarded as the "chips" of biomanufacturing that will fuel the emerging bioeconomy era. Biotechnology advances have led to the screening, investigation, and engineering of an increasing number of microorganisms as diverse MCFs, which are the workhorses of biomanufacturing and help develop the bioeconomy. This review briefly summarizes the progress and strategies in the development of robust and efficient MCFs for sustainable and economic biomanufacturing. First, a comprehensive understanding of microbial chassis cells, including accurate genome sequences and corresponding annotations; metabolic and regulatory networks governing substances, energy, physiology, and information; and their similarity and uniqueness compared with those of other microorganisms, is needed. Moreover, the development and application of effective and efficient tools is crucial for engineering both model and nonmodel microbial chassis cells into efficient MCFs, including the identification and characterization of biological parts, as well as the design, synthesis, assembly, editing, and regulation of genes, circuits, and pathways. This review also highlights the necessity of integrating automation and artificial intelligence (AI) with biotechnology to facilitate the development of future customized artificial synthetic MCFs to expedite the industrialization process of biomanufacturing and the bioeconomy.

Construction of highly active and stable recombinant nattokinase by engineered bacteria and computational design

Nattokinase (NK) is an enzyme that has been recognized as a new potential thrombolytic drug due to its strong thrombolytic activity. However, it is difficult to maintain the enzyme activity of NK during high temperature environment of industrial production. In this study, we constructed six NK mutants with potential for higher thermostability using a rational protein engineering strategy integrating free energy-based methods and molecular dynamics (MD) simulation. Then, wild-type NK and NK mutants were expressed in Escherichia coli (E. coli), and their thermostability and thrombolytic activity were tested. The results showed that, compared with wild-type NK, the mutants Y256P, Q206L and E156F all had improved thermostability. The optimal mutant Y256P showed a higher melting temperature (Tm) of 77.4 °C, an increase of 4 °C in maximum heat-resistant temperature and an increase of 51.8 % in activity at 37 °C compared with wild-type NK. Moreover, we also explored the mechanism of the increased thermostability of these mutants by analysing the MD trajectories under different simulation temperatures.

Using Protein Design and Directed Evolution to Monomerize a Bright Near-Infrared Fluorescent Protein

The small ultrared fluorescent protein (smURFP) is a bright near-infrared (NIR) fluorescent protein (FP) that forms a dimer and binds its fluorescence chromophore, biliverdin, at its dimer interface. To engineer a monomeric NIR FP based on smURFP potentially more suitable for bioimaging, we employed protein design to extend the protein backbone with a new segment of two helices that shield the original dimer interface while covering the biliverdin binding pocket in place of the second chain in the original dimer. We experimentally characterized 13 designs and obtained a monomeric protein with a weak fluorescence. We enhanced the fluorescence of this designed protein through two rounds of directed evolution and obtained designed monomeric smURFP (DMsmURFP), a bright, stable, and monomeric NIR FP with a molecular weight of 19.6 kDa. We determined the crystal structures of DMsmURFP both in the apo state and in complex with biliverdin, which confirmed the designed structure. The use of DMsmURFP in in vivo imaging of mammalian systems was demonstrated. The backbone design-based strategy used here can also be applied to monomerize other naturally multimeric proteins with intersubunit functional sites.

De novo design of cavity-containing proteins with a backbone-centered neural network energy function

The design of small-molecule-binding proteins requires protein backbones that contain cavities. Previous design efforts were based on naturally occurring cavity-containing backbone architectures. Here, we designed diverse cavity-containing backbones without predefined architectures by introducing tailored restraints into the backbone sampling driven by SCUBA (Side Chain-Unknown Backbone Arrangement), a neural network statistical energy function. For 521 out of 5816 designs, the root-mean-square deviations (RMSDs) of the Cα atoms for the AlphaFold2-predicted structures and our designed structures are within 2.0 Å. We experimentally tested 10 designed proteins and determined the crystal structures of two of them. One closely agrees with the designed model, while the other forms a domain-swapped dimer, where the partial structures are in agreement with the designed structures. Our results indicate that data-driven methods such as SCUBA hold great potential for designing de novo proteins with tailored small-molecule-binding function.

Sparks of function by de novo protein design

Information in proteins flows from sequence to structure to function, with each step causally driven by the preceding one. Protein design is founded on inverting this process: specify a desired function, design a structure executing this function, and find a sequence that folds into this structure. This 'central dogma' underlies nearly all de novo protein-design efforts. Our ability to accomplish these tasks depends on our understanding of protein folding and function and our ability to capture this understanding in computational methods. In recent years, deep learning-derived approaches for efficient and accurate structure modeling and enrichment of successful designs have enabled progression beyond the design of protein structures and towards the design of functional proteins. We examine these advances in the broader context of classical de novo protein design and consider implications for future challenges to come, including fundamental capabilities such as sequence and structure co-design and conformational control considering flexibility, and functional objectives such as antibody and enzyme design.

Protein sequence design on given backbones with deep learning

Deep learning methods for protein sequence design focus on modeling and sampling the many- dimensional distribution of amino acid sequences conditioned on the backbone structure. To produce physically foldable sequences, inter-residue couplings need to be considered properly. These couplings are treated explicitly in iterative methods or autoregressive methods. Non-autoregressive models treating these couplings implicitly are computationally more efficient, but still await tests by wet experiment. Currently, sequence design methods are evaluated mainly using native sequence recovery rate and native sequence perplexity. These metrics can be complemented by sequence-structure compatibility metrics obtained from energy calculation or structure prediction. However, existing computational metrics have important limitations that may render the generalization of computational test results to performance in real applications unwarranted. Validation of design methods by wet experiments should be encouraged.

SPIN-CGNN: Improved fixed backbone protein design with contact map-based graph construction and contact graph neural network

Recent advances in deep learning have significantly improved the ability to infer protein sequences directly from protein structures for the fix-backbone design. The methods have evolved from the early use of multi-layer perceptrons to convolutional neural networks, transformers, and graph neural networks (GNN). However, the conventional approach of constructing K-nearest-neighbors (KNN) graph for GNN has limited the utilization of edge information, which plays a critical role in network performance. Here we introduced SPIN-CGNN based on protein contact maps for nearest neighbors. Together with auxiliary edge updates and selective kernels, we found that SPIN-CGNN provided a comparable performance in refolding ability by AlphaFold2 to the current state-of-the-art techniques but a significant improvement over them in term of sequence recovery, perplexity, deviation from amino-acid compositions of native sequences, conservation of hydrophobic positions, and low complexity regions, according to the test by unseen structures, "hallucinated" structures and diffusion models. Results suggest that low complexity regions in the sequences designed by deep learning, for generated structures in particular, remain to be improved, when compared to the native sequences.

Rational Redesign of Chitosanase to Enhance Thermostability and Catalytic Activity to Produce Chitooligosaccharides with a Relatively High Degree of Polymerization

Chitooligosaccharides (hdpCOS) with a high degree of polymerization (hdp, DP 4-10) generally have greater biological activities than those of low-DP (ldp, DP 2-3) COS. Chitosanase from Bacillus amyloliquefaciens KCP2 (Csn46) can degrade chitosan to more hdpCOS at high temperature (70 °C), but low thermal stability at this temperature makes it unsuitable for industrial application; the wild-type enzyme can only produce COS (DP 2-4) at lower temperatures. Several thermostable mutants were obtained by modifying chitosanase using a comprehensive strategy based on a computer-aided mutant design. A combination of four beneficial single-point mutations (A129L/T175 V/K70T/D34G) to Csn46 was selected to obtain a markedly improved mutant, Mut4, with a half-life at 60 °C extended from 34.31 to 690.80 min, and the specific activity increased from 1671.73 to 3528.77 U/mg. Mut4 produced COS with DPs of 2-4 and 2-7 at 60 and 70 °C, respectively. Therefore, Mut4 has the potential to be applied to the industrial-scale preparation of hdpCOS with high biological activity.

Exploring binding positions and backbone conformations of peptide ligands of proteins with a backbone-centred statistical energy function

When designing peptide ligands based on the structure of a protein receptor, it can be very useful to narrow down the possible binding positions and bound conformations of the ligand without the need to choose its amino acid sequence in advance. Here, we construct and benchmark a tool for this purpose based on a recently reported statistical energy model named SCUBA (Sidechain-Unknown Backbone Arrangement) for designing protein backbones without considering specific amino acid sequences. With this tool, backbone fragments of different local conformation types are generated and optimized with SCUBA-driven stochastic simulations and simulated annealing, and then ranked and clustered to obtain representative backbone fragment poses of strong SCUBA interaction energies with the receptor. We computationally benchmarked the tool on 111 known protein-peptide complex structures. When the bound ligands are in the strand conformation, the method is able to generate backbone fragments of both low SCUBA energies and low root mean square deviations from experimental structures of peptide ligands. When the bound ligands are helices or coils, low-energy backbone fragments with binding poses similar to experimental structures have been generated for approximately 50% of benchmark cases. We have examined a number of predicted ligand-receptor complexes by atomistic molecular dynamics simulations, in which the peptide ligands have been found to stay at the predicted binding sites and to maintain their local conformations. These results suggest that promising backbone structures of peptides bound to protein receptors can be designed by identifying outstanding minima on the SCUBA-modeled backbone energy landscape.

ZetaDesign: an end-to-end deep learning method for protein sequence design and side-chain packing

Computational protein design has been demonstrated to be the most powerful tool in the last few years among protein designing and repacking tasks. In practice, these two tasks are strongly related but often treated separately. Besides, state-of-the-art deep-learning-based methods cannot provide interpretability from an energy perspective, affecting the accuracy of the design. Here we propose a new systematic approach, including both a posterior probability and a joint probability parts, to solve the two essential questions once for all. This approach takes the physicochemical property of amino acids into consideration and uses the joint probability model to ensure the convergence between structure and amino acid type. Our results demonstrated that this method could generate feasible, high-confidence sequences with low-energy side conformations. The designed sequences can fold into target structures with high confidence and maintain relatively stable biochemical properties. The side chain conformation has a significantly lower energy landscape without delegating to a rotamer library or performing the expensive conformational searches. Overall, we propose an end-to-end method that combines the advantages of both deep learning and energy-based methods. The design results of this model demonstrate high efficiency, and precision, as well as a low energy state and good interpretability.

De Novo Computational Design of a Lipase with Hydrolysis Activity towards Middle-Chained Fatty Acid Esters

Innovations in biocatalysts provide great prospects for intolerant environments or novel reactions. Due to the limited catalytic capacity and the long-term and labor-intensive characteristics of mining enzymes with the desired functions, de novo enzyme design was developed to obtain industrial application candidates in a rapid and convenient way. Here, based on the catalytic mechanisms and the known structures of proteins, we proposed a computational protein design strategy combining de novo enzyme design and laboratory-directed evolution. Starting with the theozyme constructed using a quantum-mechanical approach, the theoretical enzyme-skeleton combinations were assembled and optimized via the Rosetta "inside-out" protocol. A small number of designed sequences were experimentally screened using SDS-PAGE, mass spectrometry and a qualitative activity assay in which the designed enzyme 1a8uD₁ exhibited a measurable hydrolysis activity of 24.25 ± 0.57 U/g towards p-nitrophenyl octanoate. To improve the activity of the designed enzyme, molecular dynamics simulations and the RosettaDesign application were utilized to further optimize the substrate binding mode and amino acid sequence, thus keeping the residues of theozyme intact. The redesigned lipase 1a8uD₁-M8 displayed enhanced hydrolysis activity towards p-nitrophenyl octanoate-3.34 times higher than that of 1a8uD₁. Meanwhile, the natural skeleton protein (PDB entry 1a8u) did not display any hydrolysis activity, confirming that the hydrolysis abilities of the designed 1a8uD₁ and the redesigned 1a8uD₁-M8 were devised from scratch. More importantly, the designed 1a8uD₁-M8 was also able to hydrolyze the natural middle-chained substrate (glycerol trioctanoate), for which the activity was 27.67 ± 0.69 U/g. This study indicates that the strategy employed here has great potential to generate novel enzymes exhibiting the desired reactions.

Exploiting Sequence-Dependent Rotamer Information in Global Optimization of Proteins

Rotamers, namely amino acid side chain conformations common to many different peptides, can be compiled into libraries. These rotamer libraries are used in protein modeling, where the limited conformational space occupied by amino acid side chains is exploited. Here, we construct a sequence-dependent rotamer library from simulations of all possible tripeptides, which provides rotameric states dependent on adjacent amino acids. We observe significant sensitivity of rotamer populations to sequence and find that the library is successful in locating side chain conformations present in crystal structures. The library is designed for applications with basin-hopping global optimization, where we use it to propose moves in conformational space. The addition of rotamer moves significantly increases the efficiency of protein structure prediction within this framework, and we determine parameters to optimize efficiency.

Rotamer-free protein sequence design based on deep learning and self-consistency

Several previously proposed deep learning methods to design amino acid sequences that autonomously fold into a given protein backbone yielded promising results in computational tests but did not outperform conventional energy function-based methods in wet experiments. Here we present the ABACUS-R method, which uses an encoder-decoder network trained using a multitask learning strategy to predict the sidechain type of a central residue from its three-dimensional local environment, which includes, besides other features, the types but not the conformations of the surrounding sidechains. This eliminates the need to reconstruct and optimize sidechain structures, and drastically simplifies the sequence design process. Thus iteratively applying the encoder-decoder to different central residues is able to produce self-consistent overall sequences for a target backbone. Results of wet experiments, including five structures solved by X-ray crystallography, show that ABACUS-R outperforms state-of-the-art energy function-based methods in success rate and design precision.

DEPACT and PACMatch: A Workflow of Designing De Novo Protein Pockets to Bind Small Molecules

Engineering of new functional proteins such as enzymes and biosensors involves the design of new protein pockets for the specific binding of small molecules. Here, we report a workflow composed of two new computational methods to execute this task. The DEPACT (Design Pocket as a Cluster based on Templates) method is a data-driven approach to design and evaluate small-molecule-binding pockets as isolated clusters, while the PACMatch method is a computational approach to match pocket residues in a cluster model to positions on given protein scaffolds. Using DEPACT and its scoring function, pocket clusters of natural-pocket-like chemical compositions and protein-ligand interaction strength can be designed. DEPACT can design pocket clusters containing water- or metal-ion-mediated protein-ligand interactions. While being able to efficiently treat relatively large pocket cluster models (e.g., of around 10 pocket residues), PACMatch outperforms previous methods in test cases of recovering the native positions of pocket residues in natural enzyme-substrate complexes.

A backbone-centred energy function of neural networks for protein design

A protein backbone structure is designable if a substantial number of amino acid sequences exist that autonomously fold into it^1,2. It has been suggested that the designability of backbones is governed mainly by side chain-independent or side chain type-insensitive molecular interactions^3-5, indicating an approach for designing new backbones (ready for amino acid selection) based on continuous sampling and optimization of the backbone-centred energy surface. However, a sufficiently comprehensive and precise energy function has yet to be established for this purpose. Here we show that this goal is met by a statistical model named SCUBA (for Side Chain-Unknown Backbone Arrangement) that uses neural network-form energy terms. These terms are learned with a two-step approach that comprises kernel density estimation followed by neural network training and can analytically represent multidimensional, high-order correlations in known protein structures. We report the crystal structures of nine de novo proteins whose backbones were designed to high precision using SCUBA, four of which have novel, non-natural overall architectures. By eschewing use of fragments from existing protein structures, SCUBA-driven structure design facilitates far-reaching exploration of the designable backbone space, thus extending the novelty and diversity of the proteins amenable to de novo design.

De novo protein design by an energy function based on series expansion in distance and orientation dependence

MOTIVATION

Despite many successes, de novo protein design is not yet a solved problem as its success rate remains low. The low success rate is largely because we do not yet have an accurate energy function for describing the solvent-mediated interaction between amino acid residues in a protein chain. Previous studies showed that an energy function based on series expansions with its parameters optimized for side-chain and loop conformations can lead to one of the most accurate methods for side chain (OSCAR) and loop prediction (LEAP). Following the same strategy, we developed an energy function based on series expansions with the parameters optimized in four separate stages (recovering single-residue types without and with orientation dependence, selecting loop decoys and maintaining the composition of amino acids). We tested the energy function for de novo design by using Monte Carlo simulated annealing.

RESULTS

The method for protein design (OSCAR-Design) is found to be as accurate as OSCAR and LEAP for side-chain and loop prediction, respectively. In de novo design, it can recover native residue types ranging from 38% to 43% depending on test sets, conserve hydrophobic/hydrophilic residues at ∼75%, and yield the overall similarity in amino acid compositions at more than 90%. These performance measures are all statistically significantly better than several protein design programs compared. Moreover, the largest hydrophobic patch areas in designed proteins are near or smaller than those in native proteins. Thus, an energy function based on series expansion can be made useful for protein design.

AVAILABILITY AND IMPLEMENTATION

The Linux executable version is freely available for academic users at http://zhouyq-lab.szbl.ac.cn/resources/.

De novo sequence redesign of a functional Ras-binding domain globally inverted the surface charge distribution and led to extreme thermostability

To acquire extremely thermostable proteins of given functions is challenging for conventional protein engineering. Here we applied ABACUS, a statistical energy function we developed for de novo amino acid sequence design, to globally redesign a Ras-binding domain (RBD), and obtained an extremely thermostable RBD that unfolds reversibly at above 110°C, the redesigned RBD experimentally confirmed to have expected structure and Ras-binding interface. Directed evolution of the redesigned RBD improved its Ras-binding affinity to the native protein level without excessive loss of thermostability. The designed amino acid substitutions were mostly at the protein surface. For many substitutions, strong epistasis or significantly differentiated effects on thermostability in the native sequence context relative to the redesigned sequence context were observed, suggesting the globally redesigned sequence to be unreachable through combining beneficial mutations of the native sequence. Further analyses revealed that by replacing 38 of a total of 48 non-interfacial surface residues at once, ABACUS redesign was able to globally "invert" the protein's charge distribution pattern in an optimized way. Our study demonstrates that computational protein design provides powerful new tools to solve challenging protein engineering problems.

FASPR: an open-source tool for fast and accurate protein side-chain packing

MOTIVATION

Protein structure and function are essentially determined by how the side-chain atoms interact with each other. Thus, accurate protein side-chain packing (PSCP) is a critical step toward protein structure prediction and protein design. Despite the importance of the problem, however, the accuracy and speed of current PSCP programs are still not satisfactory.

RESULTS

We present FASPR for fast and accurate PSCP by using an optimized scoring function in combination with a deterministic searching algorithm. The performance of FASPR was compared with four state-of-the-art PSCP methods (CISRR, RASP, SCATD and SCWRL4) on both native and non-native protein backbones. For the assessment on native backbones, FASPR achieved a good performance by correctly predicting 69.1% of all the side-chain dihedral angles using a stringent tolerance criterion of 20°, compared favorably with SCWRL4, CISRR, RASP and SCATD which successfully predicted 68.8%, 68.6%, 67.8% and 61.7%, respectively. Additionally, FASPR achieved the highest speed for packing the 379 test protein structures in only 34.3 s, which was significantly faster than the control methods. For the assessment on non-native backbones, FASPR showed an equivalent or better performance on I-TASSER predicted backbones and the backbones perturbed from experimental structures. Detailed analyses showed that the major advantage of FASPR lies in the optimal combination of the dead-end elimination and tree decomposition with a well optimized scoring function, which makes FASPR of practical use for both protein structure modeling and protein design studies.

AVAILABILITY AND IMPLEMENTATION

The web server, source code and datasets are freely available at https://zhanglab.ccmb.med.umich.edu/FASPR and https://github.com/tommyhuangthu/FASPR.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

To Improve Protein Sequence Profile Prediction through Image Captioning on Pairwise Residue Distance Map

Protein sequence profile prediction aims to generate multiple sequences from structural information to advance the protein design. Protein sequence profile can be computationally predicted by energy-based or fragment-based methods. By integrating these methods with neural networks, our previous method, SPIN2, has achieved a sequence recovery rate of 34%. However, SPIN2 employed only one-dimensional (1D) structural properties that are not sufficient to represent three-dimensional (3D) structures. In this study, we represented 3D structures by 2D maps of pairwise residue distances and developed a new method (SPROF) to predict protein sequence profiles based on an image captioning learning frame. To our best knowledge, this is the first method to employ a 2D distance map for predicting protein properties. SPROF achieved 39.8% in sequence recovery of residues on the independent test set, representing a 5.2% improvement over SPIN2. We also found the sequence recovery increased with the number of their neighbored residues in 3D structural space, indicating that our method can effectively learn long-range information from the 2D distance map. Thus, such network architecture using a 2D distance map is expected to be useful for other 3D structure-based applications, such as binding site prediction, protein function prediction, and protein interaction prediction. The online server and the source code is available at http://biomed.nscc-gz.cn and https://github.com/biomed-AI/SPROF , respectively.