Macromolecular modeling and design in Rosetta: recent methods and frameworks

Growing Glycans in Rosetta: Accurate de novo glycan modeling, density fitting, and rational sequon design

Carbohydrates and glycoproteins modulate key biological functions. However, experimental structure determination of sugar polymers is notoriously difficult. Computational approaches can aid in carbohydrate structure prediction, structure determination, and design. In this work, we developed a glycan-modeling algorithm, GlycanTreeModeler, that computationally builds glycans layer-by-layer, using adaptive kernel density estimates (KDE) of common glycan conformations derived from data in the Protein Data Bank (PDB) and from quantum mechanics (QM) calculations. GlycanTreeModeler was benchmarked on a test set of glycan structures of varying lengths, or "trees". Structures predicted by GlycanTreeModeler agreed with native structures at high accuracy for both de novo modeling and experimental density-guided building. We employed these tools to design de novo glycan trees into a protein nanoparticle vaccine to shield regions of the scaffold from antibody recognition, and experimentally verified shielding. This work will inform glycoprotein model prediction, glycan masking, and further aid computational methods in experimental structure determination and refinement.

Antibody design using deep learning: from sequence and structure design to affinity maturation

Deep learning has achieved impressive results in various fields such as computer vision and natural language processing, making it a powerful tool in biology. Its applications now encompass cellular image classification, genomic studies and drug discovery. While drug development traditionally focused deep learning applications on small molecules, recent innovations have incorporated it in the discovery and development of biological molecules, particularly antibodies. Researchers have devised novel techniques to streamline antibody development, combining in vitro and in silico methods. In particular, computational power expedites lead candidate generation, scaling and potential antibody development against complex antigens. This survey highlights significant advancements in protein design and optimization, specifically focusing on antibodies. This includes various aspects such as design, folding, antibody-antigen interactions docking and affinity maturation.

A Perspective on Protein Structure Prediction Using Quantum Computers

Despite the recent advancements by deep learning methods such as AlphaFold2, in silico protein structure prediction remains a challenging problem in biomedical research. With the rapid evolution of quantum computing, it is natural to ask whether quantum computers can offer some meaningful benefits for approaching this problem. Yet, identifying specific problem instances amenable to quantum advantage and estimating the quantum resources required are equally challenging tasks. Here, we share our perspective on how to create a framework for systematically selecting protein structure prediction problems that are amenable for quantum advantage, and estimate quantum resources for such problems on a utility-scale quantum computer. As a proof-of-concept, we validate our problem selection framework by accurately predicting the structure of a catalytic loop of the Zika Virus NS3 Helicase, on quantum hardware.

Engineering and evaluation of FXa bypassing agents that restore hemostasis following Apixaban associated bleeding

Direct oral anticoagulants (DOACs) targeting activated factor Xa (FXa) are used to prevent or treat thromboembolic disorders. DOACs reversibly bind to FXa and inhibit its enzymatic activity. However, DOAC treatment carries the risk of anticoagulant-associated bleeding. Currently, only one specific agent, andexanet alfa, is approved to reverse the anticoagulant effects of FXa-targeting DOACs (FXaDOACs) and control life-threatening bleeding. However, because of its mechanism of action, andexanet alfa requires a cumbersome dosing schedule, and its use is associated with the risk of thrombosis. Here, we present the computational design, engineering, and evaluation of FXa-variants that exhibit anticoagulation reversal activity in the presence of FXaDOACs. Our designs demonstrate low DOAC binding affinity, retain FXa-enzymatic activity and reduce the DOAC-associated bleeding by restoring hemostasis in mice treated with apixaban. Importantly, the FXaDOACs reversal agents we designed, unlike andexanet alfa, do not inhibit TFPI, and consequently, may have a safer thrombogenic profile.

Can AlphaFold's breakthrough in protein structure help decode the fundamental principles of adaptive cellular immunity?

T cells are essential immune cells responsible for identifying and eliminating pathogens. Through interactions between their T-cell antigen receptors (TCRs) and antigens presented by major histocompatibility complex molecules (MHCs) or MHC-like molecules, T cells discriminate foreign and self peptides. Determining the fundamental principles that govern these interactions has important implications in numerous medical contexts. However, reconstructing a map between T cells and their antagonist antigens remains an open challenge for the field of immunology, and success of in silico reconstructions of this relationship has remained incremental. In this Perspective, we discuss the role that new state-of-the-art deep-learning models for predicting protein structure may play in resolving some of the unanswered questions the field faces linking TCR and peptide-MHC properties to T-cell specificity. We provide a comprehensive overview of structural databases and the evolution of predictive models, and highlight the breakthrough AlphaFold provided the field.

PPIscreenML: Structure-based screening for protein-protein interactions using AlphaFold

Protein-protein interactions underlie nearly all cellular processes. With the advent of protein structure prediction methods such as AlphaFold2 (AF2), models of specific protein pairs can be built extremely accurately in most cases. However, determining the relevance of a given protein pair remains an open question. It is presently unclear how to use best structure-based tools to infer whether a pair of candidate proteins indeed interact with one another: ideally, one might even use such information to screen amongst candidate pairings to build up protein interaction networks. Whereas methods for evaluating quality of modeled protein complexes have been co-opted for determining which pairings interact (e.g., pDockQ and iPTM), there have been no rigorously benchmarked methods for this task. Here we introduce PPIscreenML, a classification model trained to distinguish AF2 models of interacting protein pairs from AF2 models of compelling decoy pairings. We find that PPIscreenML out-performs methods such as pDockQ and iPTM for this task, and further that PPIscreenML exhibits impressive performance when identifying which ligand/receptor pairings engage one another across the structurally conserved tumor necrosis factor superfamily (TNFSF). Analysis of benchmark results using complexes not seen in PPIscreenML development strongly suggest that the model generalizes beyond training data, making it broadly applicable for identifying new protein complexes based on structural models built with AF2.

Mechanism of DNA unwinding by MCM8-9 in complex with HROB

HROB promotes the MCM8-9 helicase in DNA damage response. To understand how HROB activates MCM8-9, we defined their interaction interface. We showed that HROB makes important yet transient contacts with both MCM8 and MCM9, and binds the MCM8-9 heterodimer with the highest affinity. MCM8-9-HROB prefer branched DNA structures, and display low DNA unwinding processivity. MCM8-9 unwinds DNA as a hexamer that assembles from dimers on DNA in the presence of ATP. The hexamer involves two repeating protein-protein interfaces between the alternating MCM8 and MCM9 subunits. One of these interfaces is quite stable and forms an obligate heterodimer across which HROB binds. The other interface is labile and mediates hexamer assembly, independently of HROB. The ATPase site formed at the labile interface contributes disproportionally more to DNA unwinding than that at the stable interface. Here, we show that HROB promotes DNA unwinding downstream of MCM8-9 loading and ring formation on ssDNA.

The hinge-engineered IgG1-IgG3 hybrid subclass IgGh47 potently enhances Fc-mediated function of anti-streptococcal and SARS-CoV-2 antibodies

Streptococcus pyogenes can cause invasive disease with high mortality despite adequate antibiotic treatments. To address this unmet need, we have previously generated an opsonic IgG1 monoclonal antibody, Ab25, targeting the bacterial M protein. Here, we engineer the IgG2-4 subclasses of Ab25. Despite having reduced binding, the IgG3 version promotes stronger phagocytosis of bacteria. Using atomic simulations, we show that IgG3's Fc tail has extensive movement in 3D space due to its extended hinge region, possibly facilitating interactions with immune cells. We replaced the hinge of IgG1 with four different IgG3-hinge segment subclasses, IgGhxx. Hinge-engineering does not diminish binding as with IgG3 but enhances opsonic function, where a 47 amino acid hinge is comparable to IgG3 in function. IgGh47 shows improved protection against S. pyogenes in a systemic infection mouse model, suggesting that IgGh47 has promise as a preclinical therapeutic candidate. Importantly, the enhanced opsonic function of IgGh47 is generalizable to diverse S. pyogenes strains from clinical isolates. We generated IgGh47 versions of anti-SARS-CoV-2 mAbs to broaden the biological applicability, and these also exhibit strongly enhanced opsonic function compared to the IgG1 subclass. The improved function of the IgGh47 subclass in two distant biological systems provides new insights into antibody function.

ALDELE: All-Purpose Deep Learning Toolkits for Predicting the Biocatalytic Activities of Enzymes

Rapidly predicting enzyme properties for catalyzing specific substrates is essential for identifying potential enzymes for industrial transformations. The demand for sustainable production of valuable industry chemicals utilizing biological resources raised a pressing need to speed up biocatalyst screening using machine learning techniques. In this research, we developed an all-purpose deep-learning-based multiple-toolkit (ALDELE) workflow for screening enzyme catalysts. ALDELE incorporates both structural and sequence representations of proteins, alongside representations of ligands by subgraphs and overall physicochemical properties. Comprehensive evaluation demonstrated that ALDELE can predict the catalytic activities of enzymes, and particularly, it identifies residue-based hotspots to guide enzyme engineering and generates substrate heat maps to explore the substrate scope for a given biocatalyst. Moreover, our models notably match empirical data, reinforcing the practicality and reliability of our approach through the alignment with confirmed mutation sites. ALDELE offers a facile and comprehensive solution by integrating different toolkits tailored for different purposes at affordable computational cost and therefore would be valuable to speed up the discovery of new functional enzymes for their exploitation by the industry.

High-Temperature Tolerance Protein Engineering through Deep Evolution

Protein engineering aimed at increasing temperature tolerance through iterative mutagenesis and high-throughput screening is often labor-intensive. Here, we developed a deep evolution (DeepEvo) strategy to engineer protein high-temperature tolerance by generating and selecting functional sequences using deep learning models. Drawing inspiration from the concept of evolution, we constructed a high-temperature tolerance selector based on a protein language model, acting as selective pressure in the high-dimensional latent spaces of protein sequences to enrich those with high-temperature tolerance. Simultaneously, we developed a variant generator using a generative adversarial network to produce protein sequence variants containing the desired function. Afterward, the iterative process involving the generator and selector was executed to accumulate high-temperature tolerance traits. We experimentally tested this approach on the model protein glyceraldehyde 3-phosphate dehydrogenase, obtaining 8 variants with high-temperature tolerance from just 30 generated sequences, achieving a success rate of over 26%, demonstrating the high efficiency of DeepEvo in engineering protein high-temperature tolerance.

Methyl rosmarinate is an allosteric inhibitor of SARS-CoV-2 3 CL protease as a potential candidate against SARS-cov-2 infection

The coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been ongoing for more than three years and urgently needs to be addressed. Traditional Chinese medicine (TCM) prescriptions have played an important role in the clinical treatment of patients with COVID-19 in China. However, it is difficult to uncover the potential molecular mechanisms of the active ingredients in these TCM prescriptions. In this paper, we developed a new approach by integrating the experimental assay, virtual screening, and the experimental verification, exploring the rapid discovery of active ingredients from TCM prescriptions. To achieve this goal, 4 TCM prescriptions in clinical use for different indications were selected to find the antiviral active ingredients in TCMs. The 3-chymotrypsin-like protease (3CLpro), an important target for fighting COVID-19, was utilized to determine the inhibitory activity of the TCM prescriptions and single herb. It was found that 10 single herbs had better inhibitory activity than other herbs by using a fluorescence resonance energy transfer (FRET) - based enzymatic assay of SARS-CoV-2 3CLpro. The ingredients contained in 10 herbs were thus virtually screened and the predicted active ingredients were experimentally validated. Thus, such a research strategy firstly removed many single herbs with no inhibitory activity against SARS-CoV-2 3CLpro at the very beginning by FRET-based assay, making our subsequent virtual screening more effective. Finally, 4 active components were found to have stronger inhibitory effects on SARS-CoV-2 3CLpro, and their inhibitory mechanism was subsequently investigated. Among of them, methyl rosmarinate as an allosteric inhibitor of SARS-CoV-2 3CLpro was confirmed and its ability to inhibit viral replication was demonstrated by the SARS-CoV-2 replicon system. To validate the binding mode via docking, the mutation experiment, circular dichroism (CD), enzymatic inhibition and surface plasmon resonance (SPR) assay were performed, demonstrating that methyl rosmarinate bound to the allosteric site of SARS-CoV-2 3CLpro. In conclusion, this paper provides the new ideas for the rapid discovery of active ingredients in TCM prescriptions based on a specific target, and methyl rosmarinate has the potential to be developed as an antiviral therapeutic candidate against SARS-CoV-2 infection.

The structural landscape of the immunoglobulin fold by large-scale de novo design

De novo designing immunoglobulin-like frameworks that allow for functional loop diversification shows great potential for crafting antibody-like scaffolds with fully customizable structures and functions. In this work, we combined de novo parametric design with deep-learning methods for protein structure prediction and design to explore the structural landscape of 7-stranded immunoglobulin domains. After screening folding of nearly 4 million designs, we have assembled a structurally diverse library of ~50,000 immunoglobulin domains with high-confidence AlphaFold2 predictions and structures diverging from naturally occurring ones. The designed dataset enabled us to identify structural requirements for the correct folding of immunoglobulin domains, shed light on β-sheet-β-sheet rotational preferences and how these are linked to functional properties. Our approach eliminates the need for preset loop conformations and opens the route to large-scale de novo design of immunoglobulin-like frameworks.

Targeting DCAF5 suppresses SMARCB1-mutant cancer by stabilizing SWI/SNF

Whereas oncogenes can potentially be inhibited with small molecules, the loss of tumour suppressors is more common and is problematic because the tumour-suppressor proteins are no longer present to be targeted. Notable examples include SMARCB1-mutant cancers, which are highly lethal malignancies driven by the inactivation of a subunit of SWI/SNF (also known as BAF) chromatin-remodelling complexes. Here, to generate mechanistic insights into the consequences of SMARCB1 mutation and to identify vulnerabilities, we contributed 14 SMARCB1-mutant cell lines to a near genome-wide CRISPR screen as part of the Cancer Dependency Map Project1-3. We report that the little-studied gene DDB1-CUL4-associated factor 5 (DCAF5) is required for the survival of SMARCB1-mutant cancers. We show that DCAF5 has a quality-control function for SWI/SNF complexes and promotes the degradation of incompletely assembled SWI/SNF complexes in the absence of SMARCB1. After depletion of DCAF5, SMARCB1-deficient SWI/SNF complexes reaccumulate, bind to target loci and restore SWI/SNF-mediated gene expression to levels that are sufficient to reverse the cancer state, including in vivo. Consequently, cancer results not from the loss of SMARCB1 function per se, but rather from DCAF5-mediated degradation of SWI/SNF complexes. These data indicate that therapeutic targeting of ubiquitin-mediated quality-control factors may effectively reverse the malignant state of some cancers driven by disruption of tumour suppressor complexes.

Strategies for designing biocatalysts with new functions

The engineering of natural enzymes has led to the availability of a broad range of biocatalysts that can be used for the sustainable manufacturing of a variety of chemicals and pharmaceuticals. However, for many important chemical transformations there are no known enzymes that can serve as starting templates for biocatalyst development. These limitations have fuelled efforts to build entirely new catalytic sites into proteins in order to generate enzymes with functions beyond those found in Nature. This bottom-up approach to enzyme development can also reveal new fundamental insights into the molecular origins of efficient protein catalysis. In this tutorial review, we will survey the different strategies that have been explored for designing new protein catalysts. These methods will be illustrated through key selected examples, which demonstrate how highly proficient and selective biocatalysts can be developed through experimental protein engineering and/or computational design. Given the rapid pace of development in the field, we are optimistic that designer enzymes will begin to play an increasingly prominent role as industrial biocatalysts in the coming years.

Continuous evolution of compact protein degradation tags regulated by selective molecular glues

Conditional protein degradation tags (degrons) are usually >100 amino acids long or are triggered by small molecules with substantial off-target effects, thwarting their use as specific modulators of endogenous protein levels. We developed a phage-assisted continuous evolution platform for molecular glue complexes (MG-PACE) and evolved a 36-amino acid zinc finger (ZF) degron (SD40) that binds the ubiquitin ligase substrate receptor cereblon in complex with PT-179, an orthogonal thalidomide derivative. Endogenous proteins tagged in-frame with SD40 using prime editing are degraded by otherwise inert PT-179. Cryo-electron microscopy structures of SD40 in complex with ligand-bound cereblon revealed mechanistic insights into the molecular basis of SD40's activity and specificity. Our efforts establish a system for continuous evolution of molecular glue complexes and provide ZF tags that overcome shortcomings associated with existing degrons.

When Data Are Lacking: Physics-Based Inverse Design of Biopolymers Interacting with Complex, Fluid Phases

Biomolecular research traditionally revolves around comprehending the mechanisms through which peptides or proteins facilitate specific functions, often driven by their relevance to clinical ailments. This conventional approach assumes that unraveling mechanisms is a prerequisite for wielding control over functionality, which stands as the ultimate research goal. However, an alternative perspective emerges from physics-based inverse design, shifting the focus from mechanisms to the direct acquisition of functional control strategies. By embracing this methodology, we can uncover solutions that might not have direct parallels in natural systems, yet yield crucial insights into the isolated molecular elements dictating functionality. This provides a distinctive comprehension of the underlying mechanisms.In this context, we elucidate how physics-based inverse design, facilitated by evolutionary algorithms and coarse-grained molecular simulations, charts a promising course for innovating the reverse engineering of biopolymers interacting with intricate fluid phases such as lipid membranes and liquid protein phases. We introduce evolutionary molecular dynamics (Evo-MD) simulations, an approach that merges evolutionary algorithms with the Martini coarse-grained force field. This method directs the evolutionary process from random amino acid sequences toward peptides interacting with complex fluid phases such as biological lipid membranes, offering significant promises in the development of peptide-based sensors and drugs. This approach can be tailored to recognize or selectively target specific attributes such as membrane curvature, lipid composition, membrane phase (e.g., lipid rafts), and protein fluid phases. Although the resulting optimal solutions may not perfectly align with biological norms, physics-based inverse design excels at isolating relevant physicochemical principles and thermodynamic driving forces governing optimal biopolymer interaction within complex fluidic environments. In addition, we expound upon how physics-based evolution using the Evo-MD approach can be harnessed to extract the evolutionary optimization fingerprints of protein-lipid interactions from native proteins. Finally, we outline how such an approach is uniquely able to generate strategic training data for predictive neural network models that cover the whole relevant physicochemical domain. Exploring challenges, we address key considerations such as choosing a fitting fitness function to delineate the desired functionality. Additionally, we scrutinize assumptions tied to system setup, the targeted protein structure, and limitations posed by the utilized (coarse-grained) force fields and explore potential strategies for guiding evolution with limited experimental data. This discourse encapsulates the potential and remaining obstacles of physics-based inverse design, paving the way for an exciting frontier in biomolecular research.

Context-dependent design of induced-fit enzymes using deep learning generates well-expressed, thermally stable and active enzymes

The potential of engineered enzymes in industrial applications is often limited by their expression levels, thermal stability, and catalytic diversity. De novo enzyme design faces challenges due to the complexity of enzymatic catalysis. An alternative approach involves expanding natural enzyme capabilities for new substrates and parameters. Here, we introduce CoSaNN (Conformation Sampling using Neural Network), an enzyme design strategy using deep learning for structure prediction and sequence optimization. CoSaNN controls enzyme conformations to expand chemical space beyond simple mutagenesis. It employs a context-dependent approach for generating enzyme designs, considering non-linear relationships in sequence and structure space. We also developed SolvIT, a graph NN predicting protein solubility in Escherichia coli, optimizing enzyme expression selection from larger design sets. Using this method, we engineered enzymes with superior expression levels, with 54% expressed in E. coli, and increased thermal stability, with over 30% having higher Tm than the template, with no high-throughput screening. Our research underscores AI's transformative role in protein design, capturing high-order interactions and preserving allosteric mechanisms in extensively modified enzymes, and notably enhancing expression success rates. This method's ease of use and efficiency streamlines enzyme design, opening broad avenues for biotechnological applications and broadening field accessibility.

Computational design of soluble functional analogues of integral membrane proteins

De novo design of complex protein folds using solely computational means remains a significant challenge. Here, we use a robust deep learning pipeline to design complex folds and soluble analogues of integral membrane proteins. Unique membrane topologies, such as those from GPCRs, are not found in the soluble proteome and we demonstrate that their structural features can be recapitulated in solution. Biophysical analyses reveal high thermal stability of the designs and experimental structures show remarkable design accuracy. The soluble analogues were functionalized with native structural motifs, standing as a proof-of-concept for bringing membrane protein functions to the soluble proteome, potentially enabling new approaches in drug discovery. In summary, we designed complex protein topologies and enriched them with functionalities from membrane proteins, with high experimental success rates, leading to a de facto expansion of the functional soluble fold space.

Unveiling the affinity-stability relationship in anti-measles virus antibodies: a computational approach for hotspots prediction

Recent years have seen an uptick in the use of computational applications in antibody engineering. These tools have enhanced our ability to predict interactions with antigens and immunogenicity, facilitate humanization, and serve other critical functions. However, several studies highlight the concern of potential trade-offs between antibody affinity and stability in antibody engineering. In this study, we analyzed anti-measles virus antibodies as a case study, to examine the relationship between binding affinity and stability, upon identifying the binding hotspots. We leverage in silico tools like Rosetta and FoldX, along with molecular dynamics (MD) simulations, offering a cost-effective alternative to traditional in vitro mutagenesis. We introduced a pattern in identifying key residues in pairs, shedding light on hotspots identification. Experimental physicochemical analysis validated the predicted key residues by confirming significant decrease in binding affinity for the high-affinity antibodies to measles virus hemagglutinin. Through the nature of the identified pairs, which represented the relative hydropathy of amino acid side chain, a connection was proposed between affinity and stability. The findings of the study enhance our understanding of the interactions between antibody and measles virus hemagglutinin. Moreover, the implications of the observed correlation between binding affinity and stability extend beyond the field of anti-measles virus antibodies, thereby opening doors for advancements in antibody research.

A green pathway for lignin valorization: Enzymatic lignin depolymerization in biocompatible ionic liquids and deep eutectic solvents

Lignin depolymerization, which enables the breakdown of a complex and heterogeneous aromatic polymer into relatively uniform derivatives, serves as a critical process in valorization of lignin. Enzymatic lignin depolymerization has become a promising biological strategy to overcome the heterogeneity of lignin, due to its mild reaction conditions and high specificity. However, the low solubility of lignin compounds in aqueous environments prevents efficient lignin depolymerization by lignin-degrading enzymes. The employment of biocompatible ionic liquids (ILs) and deep eutectic solvents (DESs) in lignin fractionation has created a promising pathway to enzymatically depolymerize lignin within these green solvents to increase lignin solubility. In this review, recent research progress on enzymatic lignin depolymerization, particularly in a consolidated process involving ILs/DESs is summarized. In addition, the interactions between lignin-degrading enzymes and solvent systems are explored, and potential protein engineering methodology to improve the performance of lignin-degrading enzymes is discussed. Consolidation of enzymatic lignin depolymerization and biocompatible ILs/DESs paves a sustainable, efficient, and synergistic way to convert lignin into value-added products.

Creating new-to-nature carbon fixation: A guide

Synthetic biology aims at designing new biological functions from first principles. These new designs allow to expand the natural solution space and overcome the limitations of naturally evolved systems. One example is synthetic CO2-fixation pathways that promise to provide more efficient ways for the capture and conversion of CO2 than natural pathways, such as the Calvin Benson Bassham (CBB) cycle of photosynthesis. In this review, we provide a practical guideline for the design and realization of such new-to-nature CO2-fixation pathways. We introduce the concept of "synthetic CO2-fixation", and give a general overview over the enzymology and topology of synthetic pathways, before we derive general principles for their design from their eight naturally evolved analogs. We provide a comprehensive summary of synthetic carbon-assimilation pathways and derive a step-by-step, practical guide from the theoretical design to their practical implementation, before ending with an outlook on new developments in the field.

Scaffold Matcher: A CMA-ES based algorithm for identifying hotspot aligned peptidomimetic scaffolds

The design of protein interaction inhibitors is a promising approach to address aberrant protein interactions that cause disease. One strategy in designing inhibitors is to use peptidomimetic scaffolds that mimic the natural interaction interface. A central challenge in using peptidomimetics as protein interaction inhibitors, however, is determining how best the molecular scaffold aligns to the residues of the interface it is attempting to mimic. Here we present the Scaffold Matcher algorithm that aligns a given molecular scaffold onto hotspot residues from a protein interaction interface. To optimize the degrees of freedom of the molecular scaffold we implement the covariance matrix adaptation evolution strategy (CMA-ES), a state-of-the-art derivative-free optimization algorithm in Rosetta. To evaluate the performance of the CMA-ES, we used 26 peptides from the FlexPepDock Benchmark and compared with three other algorithms in Rosetta, specifically, Rosetta's default minimizer, a Monte Carlo protocol of small backbone perturbations, and a Genetic algorithm. We test the algorithms' performance on their ability to align a molecular scaffold to a series of hotspot residues (i.e., constraints) along native peptides. Of the 4 methods, CMA-ES was able to find the lowest energy conformation for all 26 benchmark peptides. Additionally, as a proof of concept, we apply the Scaffold Match algorithm with CMA-ES to align a peptidomimetic oligooxopiperazine scaffold to the hotspot residues of the substrate of the main protease of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Our implementation of CMA-ES into Rosetta allows for an alternative optimization method to be used on macromolecular modeling problems with rough energy landscapes. Finally, our Scaffold Matcher algorithm allows for the identification of initial conformations of interaction inhibitors that can be further designed and optimized as high-affinity reagents.

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.

Opportunities and Challenges for Machine Learning-Assisted Enzyme Engineering

Enzymes can be engineered at the level of their amino acid sequences to optimize key properties such as expression, stability, substrate range, and catalytic efficiency-or even to unlock new catalytic activities not found in nature. Because the search space of possible proteins is vast, enzyme engineering usually involves discovering an enzyme starting point that has some level of the desired activity followed by directed evolution to improve its "fitness" for a desired application. Recently, machine learning (ML) has emerged as a powerful tool to complement this empirical process. ML models can contribute to (1) starting point discovery by functional annotation of known protein sequences or generating novel protein sequences with desired functions and (2) navigating protein fitness landscapes for fitness optimization by learning mappings between protein sequences and their associated fitness values. In this Outlook, we explain how ML complements enzyme engineering and discuss its future potential to unlock improved engineering outcomes.

Targeting high-risk multiple myeloma genotypes with optimized anti-CD70 CAR-T cells

Despite the success of BCMA-targeting CAR-Ts in multiple myeloma, patients with high-risk cytogenetic features still relapse most quickly and are in urgent need of additional therapeutic options. Here, we identify CD70, widely recognized as a favorable immunotherapy target in other cancers, as a specifically upregulated cell surface antigen in high risk myeloma tumors. We use a structure-guided design to define a CD27-based anti-CD70 CAR-T design that outperforms all tested scFv-based CARs, leading to >80-fold improved CAR-T expansion in vivo. Epigenetic analysis via machine learning predicts key transcription factors and transcriptional networks driving CD70 upregulation in high risk myeloma. Dual-targeting CAR-Ts against either CD70 or BCMA demonstrate a potential strategy to avoid antigen escape-mediated resistance. Together, these findings support the promise of targeting CD70 with optimized CAR-Ts in myeloma as well as future clinical translation of this approach.

Phage-assisted evolution of highly active cytosine base editors with enhanced selectivity and minimal sequence context preference

TadA-derived cytosine base editors (TadCBEs) enable programmable C•G-to-T•A editing while retaining the small size, high on-target activity, and low off-target activity of TadA deaminases. Existing TadCBEs, however, exhibit residual A•T-to-G•C editing at certain positions and lower editing efficiencies at some sequence contexts and with non-SpCas9 targeting domains. To address these limitations, we use phage-assisted evolution to evolve CBE6s from a TadA-mediated dual cytosine and adenine base editor, discovering mutations at N46 and Y73 in TadA that prevent A•T-to-G•C editing and improve C•G-to-T•A editing with expanded sequence-context compatibility, respectively. In E. coli, CBE6 variants offer high C•G-to-T•A editing and no detected A•T-to-G•C editing in any sequence context. In human cells, CBE6 variants exhibit broad Cas domain compatibility and retain low off-target editing despite exceeding BE4max and previous TadCBEs in on-target editing efficiency. Finally, we show that the high selectivity of CBE6 variants is well-suited for therapeutically relevant stop codon installation without creating unwanted missense mutations from residual A•T-to-G•C editing.

A novel major facilitator superfamily-type tripartite efflux system CprABC mediates resistance to polymyxins in Chryseobacterium sp. PL22-22A

Background

Polymyxin B (PMB) and polymyxin E (colistin, CST) are polymyxin antibiotics, which are considered last-line therapeutic options against multidrug-resistant Gram-negative bacteria in serious infections. However, there is increasing risk of resistance to antimicrobial drugs. Effective efflux pump inhibitors (EPIs) should be developed to help combat efflux pump-mediated antibiotic resistance.

Methods

Chryseobacterium sp. PL22-22A was isolated from aquaculture sewage under selection with 8 mg/L PMB, and then its genome was sequenced using Oxford Nanopore and BGISEQ-500 platforms. Cpr (Chryseobacterium Polymyxins Resistance) genes encoding a major facilitator superfamily-type tripartite efflux system, were found in the genome. These genes, and the gene encoding a truncation mutant of CprB from which sequence called CprBc was deleted, were amplified and expressed/co-expressed in Escherichia coli DH5α. Minimum inhibitory concentrations (MICs) of polymyxins toward the various E. coli heterologous expression strains were tested in the presence of 2-128 mg/L PMB or CST. The pumping activity of CprABC was assessed via structural modeling using Discovery Studio 2.0 software. Moreover, the influence on MICs of baicalin, a novel MFS EPI, was determined, and the effect was analyzed based on homology modeling.

Results

Multidrug-resistant bacterial strain Chryseobacterium sp. PL22-22A was isolated in this work; it has notable resistance to polymyxin, with MICs for PMB and CST of 96 and 128 mg/L, respectively. A novel MFS-type tripartite efflux system, named CprABC, was identified in the genome of Chryseobacterium sp. PL22-22A. Heterologous expression and EPI assays indicated that the CprABC system is responsible for the polymyxin resistance of Chryseobacterium sp. PL22-22A. Structural modeling suggested that this efflux system provides a continuous conduit that runs from the CprB funnel through the CprC porin domain to pump polymyxins out of the cell. A specific C-terminal α-helix, CprBc, has an activation function on polymyxin excretion by CprB. The flavonoid compound baicalin was found to affect the allostery of CprB and/or obstruct the substrate conduit, and thus to inhibit extracellular polymyxin transport by CprABC.

Conclusion

Novel MFS-type tripartite efflux system CprABC in Chryseobacterium sp. PL22-22A mediates resistance to polymyxins, and baicalin is a promising EPI.

Disordered regions and folded modules in CAF-1 promote histone deposition in Schizosaccharomyces pombe

Genome and epigenome integrity in eukaryotes depends on the proper coupling of histone deposition with DNA synthesis. This process relies on the evolutionary conserved histone chaperone CAF-1 for which the links between structure and functions are still a puzzle. While studies of the Saccharomyces cerevisiae CAF-1 complex enabled to propose a model for the histone deposition mechanism, we still lack a framework to demonstrate its generality and in particular, how its interaction with the polymerase accessory factor PCNA is operating. Here, we reconstituted a complete SpCAF-1 from fission yeast. We characterized its dynamic structure using NMR, SAXS and molecular modeling together with in vitro and in vivo functional studies on rationally designed interaction mutants. Importantly, we identify the unfolded nature of the acidic domain which folds up when binding to histones. We also show how the long KER helix mediates DNA binding and stimulates SpCAF-1 association with PCNA. Our study highlights how the organization of CAF-1 comprising both disordered regions and folded modules enables the dynamics of multiple interactions to promote synthesis-coupled histone deposition essential for its DNA replication, heterochromatin maintenance, and genome stability functions.

Predictive modeling and cryo-EM: A synergistic approach to modeling macromolecular structure

Over the last 15 years, structural biology has seen unprecedented development and improvement in two areas: electron cryo-microscopy (cryo-EM) and predictive modeling. Once relegated to low resolutions, single-particle cryo-EM is now capable of achieving near-atomic resolutions of a wide variety of macromolecular complexes. Ushered in by AlphaFold, machine learning has powered the current generation of predictive modeling tools, which can accurately and reliably predict models for proteins and some complexes directly from the sequence alone. Although they offer new opportunities individually, there is an inherent synergy between these techniques, allowing for the construction of large, complex macromolecular models. Here, we give a brief overview of these approaches in addition to illustrating works that combine these techniques for model building. These examples provide insight into model building, assessment, and limitations when integrating predictive modeling with cryo-EM density maps. Together, these approaches offer the potential to greatly accelerate the generation of macromolecular structural insights, particularly when coupled with experimental data.

Computational Design of Peptides for Biomaterials Applications

Computer-aided molecular design and protein engineering emerge as promising and active subjects in bioengineering and biotechnological applications. On one hand, due to the advancing computing power in the past decade, modeling toolkits and force fields have been put to use for accurate multiscale modeling of biomolecules including lipid, protein, carbohydrate, and nucleic acids. On the other hand, machine learning emerges as a revolutionary data analysis tool that promises to leverage physicochemical properties and structural information obtained from modeling in order to build quantitative protein structure-function relationships. We review recent computational works that utilize state-of-the-art computational methods to engineer peptides and proteins for various emerging biomedical, antimicrobial, and antifreeze applications. We also discuss challenges and possible future directions toward developing a roadmap for efficient biomolecular design and engineering.

Leveraging Artificial Intelligence to Expedite Antibody Design and Enhance Antibody-Antigen Interactions

This perspective sheds light on the transformative impact of recent computational advancements in the field of protein therapeutics, with a particular focus on the design and development of antibodies. Cutting-edge computational methods have revolutionized our understanding of protein-protein interactions (PPIs), enhancing the efficacy of protein therapeutics in preclinical and clinical settings. Central to these advancements is the application of machine learning and deep learning, which offers unprecedented insights into the intricate mechanisms of PPIs and facilitates precise control over protein functions. Despite these advancements, the complex structural nuances of antibodies pose ongoing challenges in their design and optimization. Our review provides a comprehensive exploration of the latest deep learning approaches, including language models and diffusion techniques, and their role in surmounting these challenges. We also present a critical analysis of these methods, offering insights to drive further progress in this rapidly evolving field. The paper includes practical recommendations for the application of these computational techniques, supplemented with independent benchmark studies. These studies focus on key performance metrics such as accuracy and the ease of program execution, providing a valuable resource for researchers engaged in antibody design and development. Through this detailed perspective, we aim to contribute to the advancement of antibody design, equipping researchers with the tools and knowledge to navigate the complexities of this field.

DynamicBind: predicting ligand-specific protein-ligand complex structure with a deep equivariant generative model

While significant advances have been made in predicting static protein structures, the inherent dynamics of proteins, modulated by ligands, are crucial for understanding protein function and facilitating drug discovery. Traditional docking methods, frequently used in studying protein-ligand interactions, typically treat proteins as rigid. While molecular dynamics simulations can propose appropriate protein conformations, they're computationally demanding due to rare transitions between biologically relevant equilibrium states. In this study, we present DynamicBind, a deep learning method that employs equivariant geometric diffusion networks to construct a smooth energy landscape, promoting efficient transitions between different equilibrium states. DynamicBind accurately recovers ligand-specific conformations from unbound protein structures without the need for holo-structures or extensive sampling. Remarkably, it demonstrates state-of-the-art performance in docking and virtual screening benchmarks. Our experiments reveal that DynamicBind can accommodate a wide range of large protein conformational changes and identify cryptic pockets in unseen protein targets. As a result, DynamicBind shows potential in accelerating the development of small molecules for previously undruggable targets and expanding the horizons of computational drug discovery.

A-kinase anchoring proteins are enriched in the central pair microtubules of motile cilia in Chlamydomonas reinhardtii

Cilia are microtubule-based sensory organelles present in a number of eukaryotic cells. Mutations in the genes encoding ciliary proteins cause ciliopathies in humans. A-kinase anchoring proteins (AKAPs) tether ciliary signaling proteins such as protein kinase A (PKA). The dimerization and docking domain (D/D) on the RIIα subunit of PKA interacts with AKAPs. Here, we show that AKAP240 from the central-pair microtubules of Chlamydomonas reinhardtii cilia uses two C-terminal amphipathic helices to bind to its partner FAP174, an RIIα-like protein with a D/D domain at the N-terminus. Co-immunoprecipitation using anti-FAP174 antibody with an enriched central-pair microtubule fraction isolated seven interactors whose mass spectrometry analysis revealed proteins from the C2a (FAP65, FAP70, and FAP147) and C1b (CPC1, HSP70A, and FAP42) microtubule projections and FAP75, a protein whose sub-ciliary localization is unknown. Using RII D/D and FAP174 as baits, we identified two additional AKAPs (CPC1 and FAP297) in the central-pair microtubules.

ProteinMPNN Recovers Complex Sequence Properties of Transmembrane β-barrels

Recent deep-learning (DL) protein design methods have been successfully applied to a range of protein design problems, including the de novo design of novel folds, protein binders, and enzymes. However, DL methods have yet to meet the challenge of de novo membrane protein (MP) and the design of complex β-sheet folds. We performed a comprehensive benchmark of one DL protein sequence design method, ProteinMPNN, using transmembrane and water-soluble β-barrel folds as a model, and compared the performance of ProteinMPNN to the new membrane-specific Rosetta Franklin2023 energy function. We tested the effect of input backbone refinement on ProteinMPNN performance and found that given refined and well-defined inputs, ProteinMPNN more accurately captures global sequence properties despite complex folding biophysics. It generates more diverse TMB sequences than Franklin2023 in pore-facing positions. In addition, ProteinMPNN generated TMB sequences that passed state-of-the-art in silico filters for experimental validation, suggesting that the model could be used in de novo design tasks of diverse nanopores for single-molecule sensing and sequencing. Lastly, our results indicate that the low success rate of ProteinMPNN for the design of β-sheet proteins stems from backbone input accuracy rather than software limitations.

De novo design of high-affinity binders of bioactive helical peptides

Many peptide hormones form an α-helix on binding their receptors1-4, and sensitive methods for their detection could contribute to better clinical management of disease5. De novo protein design can now generate binders with high affinity and specificity to structured proteins6,7. However, the design of interactions between proteins and short peptides with helical propensity is an unmet challenge. Here we describe parametric generation and deep learning-based methods for designing proteins to address this challenge. We show that by extending RFdiffusion8 to enable binder design to flexible targets, and to refining input structure models by successive noising and denoising (partial diffusion), picomolar-affinity binders can be generated to helical peptide targets by either refining designs generated with other methods, or completely de novo starting from random noise distributions without any subsequent experimental optimization. The RFdiffusion designs enable the enrichment and subsequent detection of parathyroid hormone and glucagon by mass spectrometry, and the construction of bioluminescence-based protein biosensors. The ability to design binders to conformationally variable targets, and to optimize by partial diffusion both natural and designed proteins, should be broadly useful.

An efflux pump in genomic island GI-M202a mediates the transfer of polymyxin B resistance in Pandoraea pnomenusa M202

BACKGROUND

Polymyxin B is considered a last-line therapeutic option against multidrug-resistant gram-negative bacteria, especially in COVID-19 coinfections or other serious infections. However, the risk of antimicrobial resistance and its spread to the environment should be brought to the forefront.

METHODS

Pandoraea pnomenusa M202 was isolated under selection with 8 mg/L polymyxin B from hospital sewage and then was sequenced by the PacBio RS II and Illumina HiSeq 4000 platforms. Mating experiments were performed to evaluate the transfer of the major facilitator superfamily (MFS) transporter in genomic islands (GIs) to Escherichia coli 25DN. The recombinant E. coli strain Mrc-3 harboring MFS transporter encoding gene FKQ53_RS21695 was also constructed. The influence of efflux pump inhibitors (EPIs) on MICs was determined. The mechanism of polymyxin B excretion mediated by FKQ53_RS21695 was investigated by Discovery Studio 2.0 based on homology modeling.

RESULTS

The MIC of polymyxin B for the multidrug-resistant bacterial strain P. pnomenusa M202, isolated from hospital sewage, was 96 mg/L. GI-M202a, harboring an MFS transporter-encoding gene and conjugative transfer protein-encoding genes of the type IV secretion system, was identified in P. pnomenusa M202. The mating experiment between M202 and E. coli 25DN reflected the transferability of polymyxin B resistance via GI-M202a. EPI and heterogeneous expression assays also suggested that the MFS transporter gene FKQ53_RS21695 in GI-M202a was responsible for polymyxin B resistance. Molecular docking revealed that the polymyxin B fatty acyl group inserts into the hydrophobic region of the transmembrane core with Pi-alkyl and unfavorable bump interactions, and then polymyxin B rotates around Tyr43 to externally display the peptide group during the efflux process, accompanied by an inward-to-outward conformational change in the MFS transporter. Additionally, verapamil and CCCP exhibited significant inhibition via competition for binding sites.

CONCLUSIONS

These findings demonstrated that GI-M202a along with the MFS transporter FKQ53_RS21695 in P. pnomenusa M202 could mediate the transmission of polymyxin B resistance.

TIMED-Design: flexible and accessible protein sequence design with convolutional neural networks

Sequence design is a crucial step in the process of designing or engineering proteins. Traditionally, physics-based methods have been used to solve for optimal sequences, with the main disadvantages being that they are computationally intensive for the end user. Deep learning-based methods offer an attractive alternative, outperforming physics-based methods at a significantly lower computational cost. In this paper, we explore the application of Convolutional Neural Networks (CNNs) for sequence design. We describe the development and benchmarking of a range of networks, as well as reimplementations of previously described CNNs. We demonstrate the flexibility of representing proteins in a three-dimensional voxel grid by encoding additional design constraints into the input data. Finally, we describe TIMED-Design, a web application and command line tool for exploring and applying the models described in this paper. The user interface will be available at the URL: https://pragmaticproteindesign.bio.ed.ac.uk/timed. The source code for TIMED-Design is available at https://github.com/wells-wood-research/timed-design.

Supercharged Phosphotriesterase for improved Paraoxon activity

Phosphotriesterases (PTEs) represent a class of enzymes capable of efficient neutralization of organophosphates (OPs), a dangerous class of neurotoxic chemicals. PTEs suffer from low catalytic activity, particularly at higher temperatures, due to low thermostability and low solubility. Supercharging, a protein engineering approach via selective mutation of surface residues to charged residues, has been successfully employed to generate proteins with increased solubility and thermostability by promoting charge-charge repulsion between proteins. We set out to overcome the challenges in improving PTE activity against OPs by employing a computational protein supercharging algorithm in Rosetta. Here, we discover two supercharged PTE variants, one negatively supercharged (with -14 net charge) and one positively supercharged (with +12 net charge) and characterize them for their thermodynamic stability and catalytic activity. We find that positively supercharged PTE possesses slight but significant losses in thermostability, which correlates to losses in catalytic efficiency at all temperatures, whereas negatively supercharged PTE possesses increased catalytic activity across 25°C-55°C while offering similar thermostability characteristic to the parent PTE. The impact of supercharging on catalytic efficiency will inform the design of shelf-stable PTE and criteria for enzyme engineering.

Computational Design of Peptide Assemblies

With the ongoing development of peptide self-assembling materials, there is growing interest in exploring novel functional peptide sequences. From short peptides to long polypeptides, as the functionality increases, the sequence space is also expanding exponentially. Consequently, attempting to explore all functional sequences comprehensively through experience and experiments alone has become impractical. By utilizing computational methods, especially artificial intelligence enhanced molecular dynamics (MD) simulation and de novo peptide design, there has been a significant expansion in the exploration of sequence space. Through these methods, a variety of supramolecular functional materials, including fibers, two-dimensional arrays, nanocages, etc., have been designed by meticulously controlling the inter- and intramolecular interactions. In this review, we first provide a brief overview of the current main computational methods and then focus on the computational design methods for various self-assembled peptide materials. Additionally, we introduce some representative protein self-assemblies to offer guidance for the design of self-assembling peptides.

[Bioinformatics analysis and validation of the interaction between PML protein and TAB1 protein]

OBJECTIVE

To analyze the interaction between PML protein and TAB1 protein using bioinformatic approaches and experimentally verify the results.

METHODS

Using Rosetta software, a 3D model of TAB1 protein was constructed through a comparative modeling approach; the secondary structure of PML protein was retrieved in the PDB database and its crystal structure and 3D structure were resolved. Zdock 3.0.2 software was used to perform protein-protein docking of PML and TAB1, and the best conformation was extracted for molecular structure analysis of the docking model. The interaction between the two proteins was detected using immunoprecipitation in α-MMC-treated M1 inflammatory macrophages.

RESULTS

When 6IMQ of PML was used as the docking site, PML protein formed 3 salt bridges, 6 hydrogen bonds and 6 hydrophobic interactions with TAB1 proteins; when 5YUF of PML was used as the docking site, PML protein formed 1 hydrogen bond, 3 electrostatic interactions and 9 hydrophobic interactions with TAB1 proteins, and both of the docking modes formed good molecular docking and interactions. In the M1 inflammatory macrophages treated with α-MMC for 4 h, positive protein bands of PML and TAB1 were detected in the cell lysates in PML-IP group.

CONCLUSION

PML protein can interact strongly with TAB1 protein.

Ferric reduction by a CYBDOM protein counteracts increased iron availability in root meristems induced by phosphorus deficiency

To mobilize sparingly available phosphorus (P) in the rhizosphere, many plant species secrete malate to release P sorbed onto (hydr)oxides of aluminum and iron (Fe). In the presence of Fe, malate can provoke Fe over-accumulation in the root apoplast, triggering a series of events that inhibit root growth. Here, we identified HYPERSENSITIVE TO LOW P1 (HYP1), a CYBDOM protein constituted of a DOMON and a cytochrome b561 domain, as critical to maintain cell elongation and meristem integrity under low P. We demonstrate that HYP1 mediates ascorbate-dependent trans-plasma membrane electron transport and can reduce ferric and cupric substrates in Xenopus laevis oocytes and in planta. HYP1 expression is up-regulated in response to P deficiency in the proximal zone of the root apical meristem. Disruption of HYP1 leads to increased Fe and callose accumulation in the root meristem and causes significant transcriptional changes in roots. We further demonstrate that HYP1 activity overcomes malate-induced Fe accumulation, thereby preventing Fe-dependent root growth arrest in response to low P. Collectively, our results uncover an ascorbate-dependent metalloreductase that is critical to protect root meristems of P-deficient plants from increased Fe availability and provide insights into the physiological function of the yet poorly characterized but ubiquitous CYBDOM proteins.

A Perspective on the Prospective Use of AI in Protein Structure Prediction

AlphaFold2 (AF2) and RoseTTaFold (RF) have revolutionized structural biology, serving as highly reliable and effective methods for predicting protein structures. This article explores their impact and limitations, focusing on their integration into experimental pipelines and their application in diverse protein classes, including membrane proteins, intrinsically disordered proteins (IDPs), and oligomers. In experimental pipelines, AF2 models help X-ray crystallography in resolving the phase problem, while complementarity with mass spectrometry and NMR data enhances structure determination and protein flexibility prediction. Predicting the structure of membrane proteins remains challenging for both AF2 and RF due to difficulties in capturing conformational ensembles and interactions with the membrane. Improvements in incorporating membrane-specific features and predicting the structural effect of mutations are crucial. For intrinsically disordered proteins, AF2's confidence score (pLDDT) serves as a competitive disorder predictor, but integrative approaches including molecular dynamics (MD) simulations or hydrophobic cluster analyses are advocated for accurate dynamics representation. AF2 and RF show promising results for oligomeric models, outperforming traditional docking methods, with AlphaFold-Multimer showing improved performance. However, some caveats remain in particular for membrane proteins. Real-life examples demonstrate AF2's predictive capabilities in unknown protein structures, but models should be evaluated for their agreement with experimental data. Furthermore, AF2 models can be used complementarily with MD simulations. In this Perspective, we propose a "wish list" for improving deep-learning-based protein folding prediction models, including using experimental data as constraints and modifying models with binding partners or post-translational modifications. Additionally, a meta-tool for ranking and suggesting composite models is suggested, driving future advancements in this rapidly evolving field.

Targeting ion channels with ultra-large library screening for hit discovery

Ion channels play a crucial role in a variety of physiological and pathological processes, making them attractive targets for drug development in diseases such as diabetes, epilepsy, hypertension, cancer, and chronic pain. Despite the importance of ion channels in drug discovery, the vastness of chemical space and the complexity of ion channels pose significant challenges for identifying drug candidates. The use of in silico methods in drug discovery has dramatically reduced the time and cost of drug development and has the potential to revolutionize the field of medicine. Recent advances in computer hardware and software have enabled the screening of ultra-large compound libraries. Integration of different methods at various scales and dimensions is becoming an inevitable trend in drug development. In this review, we provide an overview of current state-of-the-art computational chemistry methodologies for ultra-large compound library screening and their application to ion channel drug discovery research. We discuss the advantages and limitations of various in silico techniques, including virtual screening, molecular mechanics/dynamics simulations, and machine learning-based approaches. We also highlight several successful applications of computational chemistry methodologies in ion channel drug discovery and provide insights into future directions and challenges in this field.

Programmable synthetic receptors: the next-generation of cell and gene therapies

Cell and gene therapies hold tremendous promise for treating a range of difficult-to-treat diseases. However, concerns over the safety and efficacy require to be further addressed in order to realize their full potential. Synthetic receptors, a synthetic biology tool that can precisely control the function of therapeutic cells and genetic modules, have been rapidly developed and applied as a powerful solution. Delicately designed and engineered, they can be applied to finetune the therapeutic activities, i.e., to regulate production of dosed, bioactive payloads by sensing and processing user-defined signals or biomarkers. This review provides an overview of diverse synthetic receptor systems being used to reprogram therapeutic cells and their wide applications in biomedical research. With a special focus on four synthetic receptor systems at the forefront, including chimeric antigen receptors (CARs) and synthetic Notch (synNotch) receptors, we address the generalized strategies to design, construct and improve synthetic receptors. Meanwhile, we also highlight the expanding landscape of therapeutic applications of the synthetic receptor systems as well as current challenges in their clinical translation.

Applications of Molecular Dynamics Simulations in Drug Discovery

In the current drug development process, molecular dynamics (MD) simulations have proven to be very useful. This chapter provides an overview of the current applications of MD simulations in drug discovery, from detecting protein druggable sites and validating drug docking outcomes to exploring protein conformations and investigating the influence of mutations on its structure and functions. In addition, this chapter emphasizes various strategies to improve the conformational sampling efficiency in molecular dynamics simulations. With a growing computer power and developments in the production of force fields and MD techniques, the importance of MD simulations in helping the drug development process is projected to rise significantly in the future.

Revolutionizing Peptide-Based Drug Discovery: Advances in the Post-AlphaFold Era

Peptide-based drugs offer high specificity, potency, and selectivity. However, their inherent flexibility and differences in conformational preferences between their free and bound states create unique challenges that have hindered progress in effective drug discovery pipelines. The emergence of AlphaFold (AF) and Artificial Intelligence (AI) presents new opportunities for enhancing peptide-based drug discovery. We explore recent advancements that facilitate a successful peptide drug discovery pipeline, considering peptides' attractive therapeutic properties and strategies to enhance their stability and bioavailability. AF enables efficient and accurate prediction of peptide-protein structures, addressing a critical requirement in computational drug discovery pipelines. In the post-AF era, we are witnessing rapid progress with the potential to revolutionize peptide-based drug discovery such as the ability to rank peptide binders or classify them as binders/non-binders and the ability to design novel peptide sequences. However, AI-based methods are struggling due to the lack of well-curated datasets, for example to accommodate modified amino acids or unconventional cyclization. Thus, physics-based methods, such as docking or molecular dynamics simulations, continue to hold a complementary role in peptide drug discovery pipelines. Moreover, MD-based tools offer valuable insights into binding mechanisms, as well as the thermodynamic and kinetic properties of complexes. As we navigate this evolving landscape, a synergistic integration of AI and physics-based methods holds the promise of reshaping the landscape of peptide-based drug discovery.

The PfRCR complex bridges malaria parasite and erythrocyte during invasion

The symptoms of malaria occur during the blood stage of infection, when parasites invade and replicate within human erythrocytes. The PfPCRCR complex1, containing PfRH5 (refs. 2,3), PfCyRPA, PfRIPR, PfCSS and PfPTRAMP, is essential for erythrocyte invasion by the deadliest human malaria parasite, Plasmodium falciparum. Invasion can be prevented by antibodies3-6 or nanobodies1 against each of these conserved proteins, making them the leading blood-stage malaria vaccine candidates. However, little is known about how PfPCRCR functions during invasion. Here we present the structure of the PfRCR complex7,8, containing PfRH5, PfCyRPA and PfRIPR, determined by cryogenic-electron microscopy. We test the hypothesis that PfRH5 opens to insert into the membrane9, instead showing that a rigid, disulfide-locked PfRH5 can mediate efficient erythrocyte invasion. We show, through modelling and an erythrocyte-binding assay, that PfCyRPA-binding antibodies5 neutralize invasion through a steric mechanism. We determine the structure of PfRIPR, showing that it consists of an ordered, multidomain core flexibly linked to an elongated tail. We also show that the elongated tail of PfRIPR, which is the target of growth-neutralizing antibodies6, binds to the PfCSS-PfPTRAMP complex on the parasite membrane. A modular PfRIPR is therefore linked to the merozoite membrane through an elongated tail, and its structured core presents PfCyRPA and PfRH5 to interact with erythrocyte receptors. This provides fresh insight into the molecular mechanism of erythrocyte invasion and opens the way to new approaches in rational vaccine design.

Engineering status of protein for improving microbial cell factories

With the development of metabolic engineering and synthetic biology, microbial cell factories (MCFs) have provided an efficient and sustainable method to synthesize a series of chemicals from renewable feedstocks. However, the efficiency of MCFs is usually limited by the inappropriate status of protein. Thus, engineering status of protein is essential to achieve efficient bioproduction with high titer, yield and productivity. In this review, we summarize the engineering strategies for metabolic protein status, including protein engineering for boosting microbial catalytic efficiency, protein modification for regulating microbial metabolic capacity, and protein assembly for enhancing microbial synthetic capacity. Finally, we highlight future challenges and prospects of improving microbial cell factories by engineering status of protein.

Evaluation of AlphaFold antibody-antigen modeling with implications for improving predictive accuracy

High resolution antibody-antigen structures provide critical insights into immune recognition and can inform therapeutic design. The challenges of experimental structural determination and the diversity of the immune repertoire underscore the necessity of accurate computational tools for modeling antibody-antigen complexes. Initial benchmarking showed that despite overall success in modeling protein-protein complexes, AlphaFold and AlphaFold-Multimer have limited success in modeling antibody-antigen interactions. In this study, we performed a thorough analysis of AlphaFold's antibody-antigen modeling performance on 427 nonredundant antibody-antigen complex structures, identifying useful confidence metrics for predicting model quality, and features of complexes associated with improved modeling success. Notably, we found that the latest version of AlphaFold improves near-native modeling success to over 30%, versus approximately 20% for a previous version, while increased AlphaFold sampling gives approximately 50% success. With this improved success, AlphaFold can generate accurate antibody-antigen models in many cases, while additional training or other optimization may further improve performance.

Immunoinformatics-based potential multi-peptide vaccine designing against Jamestown Canyon Virus (JCV) capable of eliciting cellular and humoral immune responses

Jamestown Canyon virus (JCV) is a deadly viral infection transmitted by various mosquito species. This mosquito-borne virus belongs to Bunyaviridae family, posing a high public health threat in the in tropical regions of the United States causing encephalitis in humans. Common symptoms of JCV include fever, headache, stiff neck, photophobia, nausea, vomiting, and seizures. Despite the availability of resources, there is currently no vaccine or drug available to combat JCV. The purpose of this study was to develop an epitope-based vaccine using immunoinformatics approaches. The vaccine aimed to be secure, efficient, bio-compatible, and capable of stimulating both innate and adaptive immune responses. In this study, the protein sequence of JCV was obtained from the NCBI database. Various bioinformatics methods, including toxicity evaluation, antigenicity testing, conservancy analysis, and allergenicity assessment were utilized to identify the most promising epitopes. Suitable linkers and adjuvant sequences were used in the design of vaccine construct. 50s ribosomal protein sequence was used as an adjuvant at the N-terminus of the construct. A total of 5 CTL, 5 HTL, and 5 linear B cell epitopes were selected based on non-allergenicity, immunological potential, and antigenicity scores to design a highly immunogenic multi-peptide vaccine construct. Strong interactions between the proposed vaccine and human immune receptors, i.e., TLR-2 and TLR-4, were revealed in a docking study using ClusPro software, suggesting their possible relevance in the immunological response to the vaccine. Immunological and physicochemical properties assessment ensured that the proposed vaccine demonstrated high immunogenicity, solubility and thermostability. Molecular dynamics simulations confirmed the strong binding affinities, as well as dynamic and structural stability of the proposed vaccine. Immune simulation suggest that the vaccine has the potential to effectively stimulate cellular and humoral immune responses to combat JCV infection. Experimental and clinical assays are required to validate the results of this study.