Principles of Protein Stability and Their Application in Computational Design.

Proteins are increasingly used in basic and applied biomedical research. Many proteins, however, are only marginally stable and can be expressed in limited amounts, thus hampering research and applications. Research has revealed the thermodynamic, cellular, and evolutionary principles and mechanisms that underlie marginal stability. With this growing understanding, computational stability design methods have advanced over the past two decades starting from methods that selectively addressed only some aspects of marginal stability. Current methods are more general and, by combining phylogenetic analysis with atomistic design, have shown drastic improvements in solubility, thermal stability, and aggregation resistance while maintaining the protein's primary molecular activity. Stability design is opening the way to rational engineering of improved enzymes, therapeutics, and vaccines and to the application of protein design methodology to large proteins and molecular activities that have proven challenging in the past.

Automated Design of Efficient and Functionally Diverse Enzyme Repertoires.

Substantial improvements in enzyme activity demand multiple mutations at spatially proximal positions in the active site. Such mutations, however, often exhibit unpredictable epistatic (non-additive) effects on activity. Here we describe FuncLib, an automated method for designing multipoint mutations at enzyme active sites using phylogenetic analysis and Rosetta design calculations. We applied FuncLib to two unrelated enzymes, a phosphotriesterase and an acetyl-CoA synthetase. All designs were active, and most showed activity profiles that significantly differed from the wild-type and from one another. Several dozen designs with only 3-6 active-site mutations exhibited 10- to 4,000-fold higher efficiencies with a range of alternative substrates, including hydrolysis of the toxic organophosphate nerve agents soman and cyclosarin and synthesis of butyryl-CoA. FuncLib is implemented as a web server (http://FuncLib.weizmann.ac.il); it circumvents iterative, high-throughput experimental screens and opens the way to designing highly efficient and diverse catalytic repertoires.

Integrating comparative modeling and accelerated simulations reveals conformational and energetic basis of actomyosin force generation.

Muscle contraction is performed by arrays of contractile proteins in the sarcomere. Serious heart diseases, such as cardiomyopathy, can often be results of mutations in myosin and actin. Direct characterization of how small changes in the myosin-actin complex impact its force production remains challenging. Molecular dynamics (MD) simulations, although capable of studying protein structure-function relationships, are limited owing to the slow timescale of the myosin cycle as well as a lack of various intermediate structures for the actomyosin complex. Here, employing comparative modeling and enhanced sampling MD simulations, we show how the human cardiac myosin generates force during the mechanochemical cycle. Initial conformational ensembles for different myosin-actin states are learned from multiple structural templates with Rosetta. This enables us to efficiently sample the energy landscape of the system using Gaussian accelerated MD. Key myosin loop residues, whose substitutions are related to cardiomyopathy, are identified to form stable or metastable interactions with the actin surface. We find that the actin-binding cleft closure is allosterically coupled to the myosin motor core transitions and ATP-hydrolysis product release from the active site. Furthermore, a gate between switch I and switch II is suggested to control phosphate release at the prepowerstroke state. Our approach demonstrates the ability to link sequence and structural information to motor functions.

Harnessing Deep Learning Methods for Voltage-Gated Ion Channel Drug Discovery.

Voltage-gated ion channels (VGICs) are pivotal in regulating electrical activity in excitable cells and are critical pharmaceutical targets for treating many diseases including cardiac arrhythmia and neuropathic pain. Despite their significance, challenges such as achieving target selectivity persist in VGIC drug development. Recent progress in deep learning, particularly diffusion models, has enabled the computational design of protein binders for any clinically relevant protein based solely on its structure. These developments coincide with a surge in experimental structural data for VGICs, providing a rich foundation for computational design efforts. This review explores the recent advancements in computational protein design using deep learning and diffusion methods, focusing on their application in designing protein binders to modulate VGIC activity. We discuss the potential use of these methods to computationally design protein binders targeting different regions of VGICs, including the pore domain, voltage-sensing domains, and interface with auxiliary subunits. We provide a comprehensive overview of the different design scenarios, discuss key structural considerations, and address the practical challenges in developing VGIC-targeting protein binders. By exploring these innovative computational methods, we aim to provide a framework for developing novel strategies that could significantly advance VGIC pharmacology and lead to the discovery of effective and safe therapeutics.

Massively parallel, computationally guided design of a proenzyme.

Confining the activity of a designed protein to a specific microenvironment would have broad-ranging applications, such as enabling cell type-specific therapeutic action by enzymes while avoiding off-target effects. While many natural enzymes are synthesized as inactive zymogens that can be activated by proteolysis, it has been challenging to redesign any chosen enzyme to be similarly stimulus responsive. Here, we develop a massively parallel computational design, screening, and next-generation sequencing-based approach for proenzyme design. For a model system, we employ carboxypeptidase G2 (CPG2), a clinically approved enzyme that has applications in both the treatment of cancer and controlling drug toxicity. Detailed kinetic characterization of the most effectively designed variants shows that they are inhibited by ∼80% compared to the unmodified protein, and their activity is fully restored following incubation with site-specific proteases. Introducing disulfide bonds between the pro- and catalytic domains based on the design models increases the degree of inhibition to 98% but decreases the degree of restoration of activity by proteolysis. A selected disulfide-containing proenzyme exhibits significantly lower activity relative to the fully activated enzyme when evaluated in cell culture. Structural and thermodynamic characterization provides detailed insights into the prodomain binding and inhibition mechanisms. The described methodology is general and could enable the design of a variety of proproteins with precise spatial regulation.

Accurate positioning of functional residues with robotics-inspired computational protein design.

SignificanceComputational protein design promises to advance applications in medicine and biotechnology by creating proteins with many new and useful functions. However, new functions require the design of specific and often irregular atom-level geometries, which remains a major challenge. Here, we develop computational methods that design and predict local protein geometries with greater accuracy than existing methods. Then, as a proof of concept, we leverage these methods to design new protein conformations in the enzyme ketosteroid isomerase that change the protein's preference for a key functional residue. Our computational methods are openly accessible and can be applied to the design of other intricate geometries customized for new user-defined protein functions.

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.

In Silico Screening, Molecular Dynamics Simulation and Binding Free Energy Identify Single-Point Mutations That Destabilize p53 and Reduce Binding to DNA.

Considering p53's pivotal role as a tumor suppressor protein, proactive identification and characterization of potentially harmful p53 mutations are crucial before they appear in the population. To address this, four computational prediction tools-SIFT, Polyphen-2, PhD-SNP, and MutPred2-utilizing sequence-based and machine-learning algorithms, were employed to identify potentially deleterious p53 nsSNPs (nonsynonymous single nucleotide polymorphisms) that may impact p53 structure, dynamics, and binding with DNA. These computational methods identified three variants, namely, C141Y, C238S, and L265P, as detrimental to p53 stability. Furthermore, molecular dynamics (MD) simulations revealed that all three variants exhibited heightened structural flexibility compared to the native protein, especially the C141Y and L265P mutations. Consequently, due to the altered structure of mutant p53, the DNA-binding affinity of all three variants decreased by approximately 1.8 to 9.7 times compared to wild-type p53 binding with DNA (14 µM). Notably, the L265P mutation exhibited an approximately ten-fold greater reduction in binding affinity. Consequently, the presence of the L265P mutation in p53 could pose a substantial risk to humans. Given that p53 regulates abnormal tumor growth, this research carries significant implications for surveillance efforts and the development of anticancer therapies.

Predicting success using response after lead implantation with sacral neuromodulation for urgency incontinence.

IMPORTANCE: Many women report inadequate symptom control after sacral neuromodulation (SNM), despite 50% reduction in urgency incontinence episodes (UUIE) after test stimulation. OBJECTIVE: To determine the ideal percent UUIE reduction after test stimulation that predicts 24-month success. STUDY DESIGN: Using data from a multicenter SNM trial, we constructed receiver operating characteristic curves to identify an ideal threshold of percent UUIE reduction after test stimulation. We defined 24-month success as Patient Global Impression of Improvement of "very much better" to "better." We compared predictive accuracy of two models predicting success: (1) percent UUIE reduction alone and (2) with baseline characteristics. RESULTS: Of 149 women (median [IQR] baseline daily UUIE 4.7 [3.7, 6.0]), the ideal threshold for 24-month success was 72% (95% confidence interval 64,76%) UUIE reduction with accuracy 0.54 (0.42, 0.66), sensitivity 0.71 (0.56, 0.86) and specificity 0.27 (0.05, 0.55). The accuracy of the 50% reduction threshold was 0.60 (0.49, 0.71), sensitivity 0.95 (0.88, 1.0) and specificity 0.04 (0.0, 0.12). Percent reduction in UUIE was not better than chance in predicting 24-month success (concordance index [c-index] 0.47 [0.46, 0.62]); adding age, body mass index, diabetes mellitus and visual or hearing impairment the c-index was 0.68 (0.61, 0.78). CONCLUSIONS: Among women who received an internal pulse generator (IPG) due to ≥50% UUIE reduction after test stimulation, we found no ideal threshold that better predicted 24-month success. Percent reduction in UUIE after test stimulation poorly predicts 24-month success with or without clinical factors. Given this, re-evaluating how we determine who should receive an IPG is needed.

LidPose: Real-Time 3D Human Pose Estimation in Sparse Lidar Point Clouds with Non-Repetitive Circular Scanning Pattern.

In this paper, we propose a novel, vision-transformer-based end-to-end pose estimation method, LidPose, for real-time human skeleton estimation in non-repetitive circular scanning (NRCS) lidar point clouds. Building on the ViTPose architecture, we introduce novel adaptations to address the unique properties of NRCS lidars, namely, the sparsity and unusual rosetta-like scanning pattern. The proposed method addresses a common issue of NRCS lidar-based perception, namely, the sparsity of the measurement, which needs balancing between the spatial and temporal resolution of the recorded data for efficient analysis of various phenomena. LidPose utilizes foreground and background segmentation techniques for the NRCS lidar sensor to select a region of interest (RoI), making LidPose a complete end-to-end approach to moving pedestrian detection and skeleton fitting from raw NRCS lidar measurement sequences captured by a static sensor for surveillance scenarios. To evaluate the method, we have created a novel, real-world, multi-modal dataset, containing camera images and lidar point clouds from a Livox Avia sensor, with annotated 2D and 3D human skeleton ground truth.

Cloning, Expression, Purification, and Characterization of Lactate Dehydrogenase from Plasmodium knowlesi: A Zoonotic Malaria Parasite.

Plasmodium knowlesi is the only Plasmodium that causes zoonotic disease among the Plasmodium that cause infection in humans. It is fatal due to its short asexual growth cycle within 24 h. Lactate dehydrogenase (LDH), an enzyme that catalyzes the final step of glycolysis, is a biomarker for diagnosing infection by Plasmodium spp. parasite. Therefore, this study aimed to efficiently produce the soluble form of P. knowlesi LDH (PkLDH) using a bacterial expression system for studying malaria caused by P. knowlesi. Recombinant pET-21a(+)-PkLDH plasmid was constructed by inserting the PkLDH gene into a pET-21a(+) expression vector. Subsequently, the recombinant plasmid was inserted into the protein-expressing Escherichia coli Rosetta(DE3) strain, and the optimal conditions for overexpression of the PkLDH protein were established using this strain. We obtained a yield of 52.0 mg/L PkLDH from the Rosetta(DE3) strain and confirmed an activity of 483.9 U/mg through experiments. This methodology for high-efficiency PkLDH production can be utilized for the development of diagnostic methods and drug candidates for distinguishing malaria caused by P. knowlesi.

A conserved set of mutations for stabilizing soluble envelope protein dimers from dengue and Zika viruses to advance the development of subunit vaccines.

Dengue viruses (DENV serotypes 1-4) and Zika virus (ZIKV) are related flaviviruses that continue to be a public health concern, infecting hundreds of millions of people annually. The traditional live-attenuated virus vaccine approach has been challenging for the four DENV serotypes because of the need to achieve balanced replication of four independent vaccine components. Subunit vaccines represent an alternative approach that may circumvent problems inherent with live-attenuated DENV vaccines. In mature virus particles, the envelope (E) protein forms a homodimer that covers the surface of the virus and is the major target of neutralizing antibodies. Many neutralizing antibodies bind to quaternary epitopes that span across both E proteins in the homodimer. For soluble E (sE) protein to be a viable subunit vaccine, the antigens should be easy to produce and retain quaternary epitopes recognized by neutralizing antibodies. However, WT sE proteins are primarily monomeric at conditions relevant for vaccination and exhibit low expression yields. Previously, we identified amino acid mutations that stabilize the sE homodimer from DENV2 and dramatically raise expression yields. Here, we tested whether these same mutations raise the stability of sE from other DENV serotypes and ZIKV. We show that the mutations raise thermostability for sE from all the viruses, increase production yields from 4-fold to 250-fold, stabilize the homodimer, and promote binding to dimer-specific neutralizing antibodies. Our findings suggest that these sE variants could be valuable resources in the efforts to develop effective subunit vaccines for DENV serotypes 1 to 4 and ZIKV.

Structural modeling of peptide toxin-ion channel interactions using RosettaDock.

Voltage-gated ion channels play essential physiological roles in action potential generation and propagation. Peptidic toxins from animal venoms target ion channels and provide useful scaffolds for the rational design of novel channel modulators with enhanced potency and subtype selectivity. Despite recent progress in obtaining experimental structures of peptide toxin-ion channel complexes, structural determination of peptide toxins bound to ion channels in physiologically important states remains challenging. Here we describe an application of RosettaDock approach to the structural modeling of peptide toxins interactions with ion channels. We tested this approach on 10 structures of peptide toxin-ion channel complexes and demonstrated that it can sample near-native structures in all tested cases. Our approach will be useful for improving the understanding of the molecular mechanism of natural peptide toxin modulation of ion channel gating and for the structural modeling of novel peptide-based ion channel modulators.

Sparse pseudocontact shift NMR data obtained from a non-canonical amino acid-linked lanthanide tag improves integral membrane protein structure prediction.

A single experimental method alone often fails to provide the resolution, accuracy, and coverage needed to model integral membrane proteins (IMPs). Integrating computation with experimental data is a powerful approach to supplement missing structural information with atomic detail. We combine RosettaNMR with experimentally-derived paramagnetic NMR restraints to guide membrane protein structure prediction. We demonstrate this approach using the disulfide bond formation protein B (DsbB), an α-helical IMP. Here, we attached a cyclen-based paramagnetic lanthanide tag to an engineered non-canonical amino acid (ncAA) using a copper-catalyzed azide-alkyne cycloaddition (CuAAC) click chemistry reaction. Using this tagging strategy, we collected 203 backbone HN pseudocontact shifts (PCSs) for three different labeling sites and used these as input to guide de novo membrane protein structure prediction protocols in Rosetta. We find that this sparse PCS dataset combined with 44 long-range NOEs as restraints in our calculations improves structure prediction of DsbB by enhancements in model accuracy, sampling, and scoring. The inclusion of this PCS dataset improved the Cα-RMSD transmembrane segment values of the best-scoring and best-RMSD models from 9.57 Å and 3.06 Å (no NMR data) to 5.73 Å and 2.18 Å, respectively.

Enhanced soluble expression of active recombinant human interleukin-29 using champion pET SUMO system.

Current research focuses on the soluble and high-level expression of biologically active recombinant human IL-29 protein in Escherichia coli. The codon-optimized IL-29 gene was cloned into the Champion™ pET SUMO expression system downstream of the SUMO tag under the influence of the T7 lac promoter. The expression of SUMO-fused IL-29 protein was compared in E. coli Rosetta 2(DE3), Rosetta 2(DE3) pLysS, and Rosetta-gami 2(DE3). The release of the SUMO fusion partner resulted in approximately 98 mg of native rhIL-29 protein with a purity of 99% from 1 l of fermentation culture. Purified rhIL-29 was found to be biologically active, as evaluated by its anti-proliferation assay. It was found that Champion™ pET SUMO expression system can be used to obtained high yield of biologically active soluble recombinant human protein compared to other expression vector.

Discovery of Marburg virus neutralizing antibodies from virus-naïve human antibody repertoires using large-scale structural predictions.

Marburg virus (MARV) disease is lethal, with fatality rates up to 90%. Neutralizing antibodies (Abs) are promising drug candidates to prevent or treat the disease. Current efforts are focused in part on vaccine development to induce such MARV-neutralizing Abs. We analyzed the antibody repertoire from healthy unexposed and previously MARV-infected individuals to assess if naïve repertoires contain suitable precursor antibodies that could become neutralizing with a limited set of somatic mutations. We computationally searched the human Ab variable gene repertoire for predicted structural homologs of the neutralizing Ab MR78 that is specific to the receptor binding site (RBS) of MARV glycoprotein (GP). Eight Ab heavy-chain complementarity determining region 3 (HCDR3) loops from MARV-naïve individuals and one from a previously MARV-infected individual were selected for testing as HCDR3 loop chimeras on the MR78 Ab framework. Three of these chimerized antibodies bound to MARV GP. We then tested a full-length native Ab heavy chain encoding the same 17-residue-long HCDR3 loop that bound to the MARV GP the best among the chimeric Abs tested. Despite only 57% amino acid sequence identity, the Ab from a MARV-naïve donor recognized MARV GP and possessed neutralizing activity against the virus. Crystallization of both chimeric and full-length native heavy chain-containing Abs provided structural insights into the mechanism of binding for these types of Abs. Our work suggests that the MARV GP RBS is a promising candidate for epitope-focused vaccine design to induce neutralizing Abs against MARV.

Expression, Purification, and Characterization of Plasmodium vivax Lactate Dehydrogenase from Bacteria without Codon Optimization.

Plasmodium vivax is the most widespread cause of malaria, especially in subtropical and temperate regions such as Asia-Pacific and America. P. vivax lactate dehydrogenase (PvLDH), an essential enzyme in the glycolytic pathway, is required for the development and reproduction of the parasite. Thus, LDH from these parasites has garnered attention as a diagnostic biomarker for malaria and as a potential molecular target for developing antimalarial drugs. In this study, we prepared a transformed Escherichia coli strain for the overexpression of PvLDH without codon optimization. We introduced this recombinant plasmid DNA prepared by insertion of the PvLDH gene in the pET-21a(+) expression vector, into the Rosetta(DE3), an E. coli strain suitable for eukaryotic protein expression. The time, temperature, and inducer concentration for PvLDH expression from this E. coli Rosetta(DE3), containing the original PvLDH gene, were optimized. We obtained PvLDH with a 31.0 mg/L yield and high purity (>95%) from this Rosetta(DE3) strain. The purified protein was characterized structurally and functionally. The PvLDH expressed and purified from transformed bacteria without codon optimization was successfully demonstrated to exhibit its potential tetramer structure and enzyme activity. These findings are expected to provide valuable insights for research on infectious diseases, metabolism, diagnostics, and therapeutics for malaria caused by P. vivax.

Improving the Modeling of Extracellular Ligand Binding Pockets in RosettaGPCR for Conformational Selection.

G protein-coupled receptors (GPCRs) are the largest class of drug targets and undergo substantial conformational changes in response to ligand binding. Despite recent progress in GPCR structure determination, static snapshots fail to reflect the conformational space of putative binding pocket geometries to which small molecule ligands can bind. In comparative modeling of GPCRs in the absence of a ligand, often a shrinking of the orthosteric binding pocket is observed. However, the exact prediction of the flexible orthosteric binding site is crucial for adequate structure-based drug discovery. In order to improve ligand docking and guide virtual screening experiments in computer-aided drug discovery, we developed RosettaGPCRPocketSize. The algorithm creates a conformational ensemble of biophysically realistic conformations of the GPCR binding pocket between the TM bundle, which is consistent with a knowledge base of expected pocket geometries. Specifically, tetrahedral volume restraints are defined based on information about critical residues in the orthosteric binding site and their experimentally observed range of Cα-Cα-distances. The output of RosettaGPCRPocketSize is an ensemble of binding pocket geometries that are filtered by energy to ensure biophysically probable arrangements, which can be used for docking simulations. In a benchmark set, pocket shrinkage observed in the default RosettaGPCR was reduced by up to 80% and the binding pocket volume range and geometric diversity were increased. Compared to models from four different GPCR homology model databases (RosettaGPCR, GPCR-Tasser, GPCR-SSFE, and GPCRdb), the here-created models showed more accurate volumes of the orthosteric pocket when evaluated with respect to the crystallographic reference structure. Furthermore, RosettaGPCRPocketSize was able to generate an improved realistic pocket distribution. However, while being superior to other homology models, the accuracy of generated model pockets was comparable to AlphaFold2 models. Furthermore, in a docking benchmark using small-molecule ligands with a higher molecular weight between 400 and 700 Da, a higher success rate in creating native-like binding poses was observed. In summary, RosettaGPCRPocketSize can generate GPCR models with realistic orthosteric pocket volumes, which are useful for structure-based drug discovery applications.

Ab Initio Modelling of the Structure of ToxA-like and MAX Fungal Effector Proteins.

Pathogenic fungal diseases in crops are mediated by the release of effector proteins that facilitate infection. Characterising the structure of these fungal effectors is vital to understanding their virulence mechanisms and interactions with their hosts, which is crucial in the breeding of plant cultivars for disease resistance. Several effectors have been identified and validated experimentally; however, their lack of sequence conservation often impedes the identification and prediction of their structure using sequence similarity approaches. Structural similarity has, nonetheless, been observed within fungal effector protein families, creating interest in validating the use of computational methods to predict their tertiary structure from their sequence. We used Rosetta ab initio modelling to predict the structures of members of the ToxA-like and MAX effector families for which experimental structures are known to validate this method. An optimised approach was then used to predict the structures of phenotypically validated effectors lacking known structures. Rosetta was found to successfully predict the structure of fungal effectors in the ToxA-like and MAX families, as well as phenotypically validated but structurally unconfirmed effector sequences. Interestingly, potential new effector structural families were identified on the basis of comparisons with structural homologues and the identification of associated protein domains.

Recent Advances in NMR Protein Structure Prediction with ROSETTA.

Nuclear magnetic resonance (NMR) spectroscopy is a powerful method for studying the structure and dynamics of proteins in their native state. For high-resolution NMR structure determination, the collection of a rich restraint dataset is necessary. This can be difficult to achieve for proteins with high molecular weight or a complex architecture. Computational modeling techniques can complement sparse NMR datasets (<1 restraint per residue) with additional structural information to elucidate protein structures in these difficult cases. The Rosetta software for protein structure modeling and design is used by structural biologists for structure determination tasks in which limited experimental data is available. This review gives an overview of the computational protocols available in the Rosetta framework for modeling protein structures from NMR data. We explain the computational algorithms used for the integration of different NMR data types in Rosetta. We also highlight new developments, including modeling tools for data from paramagnetic NMR and hydrogen-deuterium exchange, as well as chemical shifts in CS-Rosetta. Furthermore, strategies are discussed to complement and improve structure predictions made by the current state-of-the-art AlphaFold2 program using NMR-guided Rosetta modeling.