SYNBIP: synthetic binding proteins for research, diagnosis and therapy

PROSCA: an online platform for humanized scaffold mining facilitating rational protein engineering

Protein scaffolds with small size, high stability and low immunogenicity show important applications in the field of protein engineering and design. However, no relevant computational platform has been reported yet to mining such scaffolds with the desired properties from massive protein structures in human body. Here, we developed PROSCA, a structure-based online platform dedicated to explore the space of the entire human proteome, and to discovery new privileged protein scaffolds with potential engineering value that have never been noticed. PROSCA accepts structure of protein as an input, which can be subsequently aligned with a certain class of protein structures (e.g. the human proteome either from experientially resolved or AlphaFold2 predicted structures, and the human proteins belonging to specific families or domains), and outputs humanized protein scaffolds which are structurally similar with the input protein as well as other related important information such as families, sequences, structures and expression level in human tissues. Through PROSCA, the user can also get excellent experience in visualizations of protein structures and expression overviews, and download the figures and tables of results which can be customized according to the user's needs. Along with the advanced protein engineering and selection technologies, PROSCA will facilitate the rational design of new functional proteins with privileged scaffolds. PROSCA is freely available at https://idrblab.org/prosca/.

Material-specific binding peptides empower sustainable innovations in plant health, biocatalysis, medicine and microplastic quantification

Material-binding peptides (MBPs) have emerged as a diverse and innovation-enabling class of peptides in applications such as plant-/human health, immobilization of catalysts, bioactive coatings, accelerated polymer degradation and analytics for micro-/nanoplastics quantification. Progress has been fuelled by recent advancements in protein engineering methodologies and advances in computational and analytical methodologies, which allow the design of, for instance, material-specific MBPs with fine-tuned binding strength for numerous demands in material science applications. A genetic or chemical conjugation of second (biological, chemical or physical property-changing) functionality to MBPs empowers the design of advanced (hybrid) materials, bioactive coatings and analytical tools. In this review, we provide a comprehensive overview comprising naturally occurring MBPs and their function in nature, binding properties of short man-made MBPs (<20 amino acids) mainly obtained from phage-display libraries, and medium-sized binding peptides (20-100 amino acids) that have been reported to bind to metals, polymers or other industrially produced materials. The goal of this review is to provide an in-depth understanding of molecular interactions between materials and material-specific binding peptides, and thereby empower the use of MBPs in material science applications. Protein engineering methodologies and selected examples to tailor MBPs toward applications in agriculture with a focus on plant health, biocatalysis, medicine and environmental monitoring serve as examples of the transformative power of MBPs for various industrial applications. An emphasis will be given to MBPs' role in detecting and quantifying microplastics in high throughput, distinguishing microplastics from other environmental particles, and thereby assisting to close an analytical gap in food safety and monitoring of environmental plastic pollution. In essence, this review aims to provide an overview among researchers from diverse disciplines in respect to material-(specific) binding of MBPs, protein engineering methodologies to tailor their properties to application demands, re-engineering for material science applications using MBPs, and thereby inspire researchers to employ MBPs in their research.

Computer-aided discovery of novel aryl hydrocarbon receptor ligands to regulate CYP1A1 expression in inflammatory macrophages

The environmental factor aryl hydrocarbon receptor (AhR), a key protein connecting the external environmental signals (e.g., environmental endocrine disruptor TCDD) to internal cellular processes, is involved in the activation of peripheral macrophages and inflammatory response in human body. Thus, there is widespread interest in finding compounds to anti-inflammatory response in macrophages by targeting human AhR. Here, ensemble docking based-virtual screening was first used to screen a library (~200,000 compounds) against human AhR ligand binding domain (LBD) and 25 compounds were identified as potential inhibitors. Then, 9 out of the 25 ligands were found to down-regulate the mRNA expression of CYP1A1 (a downstream gene of AhR signaling) in AhR overexpressing macrophages. The most potent compound AE-411/41415610 was selected for further study and found to reduce both mRNA and protein expressions level of CYP1A1 in mouse peritoneal macrophage. Moreover, protein chip signal pathway analysis indicated that AE-411/41415610 play a role in regulating JAK-STAT and AKT-mTOR pathways. In sum, the discovered hits with novel scaffolds provided a starting point for future design of more effective AhR-targeted lead compounds to regulate CYP1A1 expression of inflammatory peritoneal macrophages.

POLAT: Protein function prediction based on soft mask graph network and residue-Label ATtention

MOTIVATION

Elucidating protein function is a central problem in biochemistry, genetics, and molecular biology. Developing computational methods for protein function prediction is critical due to the significant gap between sequence and functional data. Recent advances in protein structure prediction, which strongly correlates with function, make it feasible to use structure to predict function. However, current structure-based methods overlook the fact that individual residues may contribute differently to the protein's function and do not take into account the correlation between protein residues and their functions. The challenge of effectively utilizing the relationship between protein residues and function-level information to predict protein function remains unsolved.

RESULT

We proposed a protein function prediction method based on Soft Mask Graph Networks and Residue-Label Attention (POLAT), which could combine sequence features, predicted structure features, and function-level information to get an accurate prediction. We use soft mask graph networks to adaptively extract the residues relevant to functions. A residue-label attention mechanism is adopted to obtain the protein-level encoded features of a protein, which are then concatenated with a protein-level embedding and fed into a dense classifier to determine the probabilities of each function. POLAT achieves 0.670, 0.515, 0.578 Fmax and 0.677, 0.409, 0.507 AUPR on the PDB cdhit test set for the MFO, BPO, and CCO domains, respectively, outperforming the existing structure-based SOTA method GAT-GO (Fmax 0.633, 0.492, 0.547; AUPR 0.660, 0.381, 0.479). POLAT is also competitive in extensive experiments among sequence-based and multimodal methods and achieves the SOTA performance in three out of six metrics.

Computational Chemistry in Structure-Based Solute Carrier Transporter Drug Design: Recent Advances and Future Perspectives

Solute carrier transporters (SLCs) are a class of important transmembrane proteins that are involved in the transportation of diverse solute ions and small molecules into cells. There are approximately 450 SLCs within the human body, and more than a quarter of them are emerging as attractive therapeutic targets for multiple complex diseases, e.g., depression, cancer, and diabetes. However, only 44 unique transporters (∼9.8% of the SLC superfamily) with 3D structures and specific binding sites have been reported. To design innovative and effective drugs targeting diverse SLCs, there are a number of obstacles that need to be overcome. However, computational chemistry, including physics-based molecular modeling and machine learning- and deep learning-based artificial intelligence (AI), provides an alternative and complementary way to the classical drug discovery approach. Here, we present a comprehensive overview on recent advances and existing challenges of the computational techniques in structure-based drug design of SLCs from three main aspects: (i) characterizing multiple conformations of the proteins during the functional process of transportation, (ii) identifying druggability sites especially the cryptic allosteric ones on the transporters for substrates and drugs binding, and (iii) discovering diverse small molecules or synthetic protein binders targeting the binding sites. This work is expected to provide guidelines for a deep understanding of the structure and function of the SLC superfamily to facilitate rational design of novel modulators of the transporters with the aid of state-of-the-art computational chemistry technologies including artificial intelligence.

Integrating Dynamic Network Analysis with AI for Enhanced Epitope Prediction in PD-L1:Affibody Interactions

Understanding binding epitopes involved in protein-protein interactions and accurately determining their structure is a long standing goal with broad applicability in industry and biomedicine. Although various experimental methods for binding epitope determination exist, these approaches are typically low throughput and cost intensive. Computational methods have potential to accelerate epitope predictions, however, recently developed artificial intelligence (AI)-based methods frequently fail to predict epitopes of synthetic binding domains with few natural homologs. Here we have developed an integrated method employing generalized-correlation-based dynamic network analysis on multiple molecular dynamics (MD) trajectories, initiated from AlphaFold2 Multimer structures, to unravel the structure and binding epitope of the therapeutic PD-L1:Affibody complex. Both AlphaFold2 and conventional molecular dynamics trajectory analysis alone each proved ineffectual in differentiating between two putative binding models referred to as parallel and perpendicular. However, our integrated approach based on dynamic network analysis showed that the perpendicular mode was significantly more stable. These predictions were validated using a suite of experimental epitope mapping protocols including cross linking mass spectrometry and next-generation sequencing-based deep mutational scanning. Our research highlights the potential of deploying dynamic network analysis to refine AI-based structure predictions for precise predictions of protein-protein interaction interfaces.

Is fragment-based graph a better graph-based molecular representation for drug design? A comparison study of graph-based models

Graph Neural Networks (GNNs) have gained significant traction in various sectors of AI-driven drug design. Over recent years, the integration of fragmentation concepts into GNNs has emerged as a potent strategy to augment the efficacy of molecular generative models. Nonetheless, challenges such as symmetry breaking and potential misrepresentation of intricate cycles and undefined functional groups raise questions about the superiority of fragment-based graph representation over traditional methods. In our research, we undertook a rigorous evaluation, contrasting the predictive prowess of eight models-developed using deep learning algorithms-across 12 benchmark datasets that span a range of properties. These models encompass established methods like GCN, AttentiveFP, and D-MPNN, as well as innovative fragment-based representation techniques. Our results indicate that fragment-based methodologies, notably PharmHGT, significantly improve model performance and interpretability, particularly in scenarios characterized by limited data availability. However, in situations with extensive training, fragment-based molecular graph representations may not necessarily eclipse traditional methods. In summation, we posit that the integration of fragmentation, as an avant-garde technique in drug design, harbors considerable promise for the future of AI-enhanced drug design.

RVvictor: Virus RNA-directed molecular interactions for RNA virus infection

RNA viruses are major human pathogens that cause seasonal epidemics and occasional pandemic outbreaks. Due to the nature of their RNA genomes, it is anticipated that virus's RNA interacts with host protein (INTPRO), messenger RNA (INTmRNA), and non-coding RNA (INTncRNA) to perform their particular functions during their transcription and replication. In other words, thus, it is urgently needed to have such valuable data on virus RNA-directed molecular interactions (especially INTPROs), which are highly anticipated to attract broad research interests in the fields of RNA virus translation and replication. In this study, a new database was constructed to describe the virus RNA-directed interaction (INTPRO, INTmRNA, INTncRNA) for RNA virus (RVvictor). This database is unique in a) unambiguously characterizing the interactions between viruses RNAs and host proteins, b) providing, for the first time, the most systematic RNA-directed interaction data resources in providing clues to understand the molecular mechanisms of RNA viruses' translation, and replication, and c) in RVvictor, comprehensive enrichment analysis is conducted for each virus RNA based on its associated target genes/proteins, and the enrichment results were explicitly illustrated using various graphs. We found significant enrichment of a suite of pathways related to infection, translation, and replication, e.g., HIV infection, coronavirus disease, regulation of viral genome replication, and so on. Due to the devastating and persistent threat posed by the RNA virus, RVvictor constructed, for the first time, a possible network of cross-talk in RNA-directed interaction, which may ultimately explain the pathogenicity of RNA virus infection. The knowledge base might help develop new anti-viral therapeutic targets in the future. It's now free and publicly accessible at: https://idrblab.org/rvvictor/.

TTD: Therapeutic Target Database describing target druggability information

Target discovery is one of the essential steps in modern drug development, and the identification of promising targets is fundamental for developing first-in-class drug. A variety of methods have emerged for target assessment based on druggability analysis, which refers to the likelihood of a target being effectively modulated by drug-like agents. In the therapeutic target database (TTD), nine categories of established druggability characteristics were thus collected for 426 successful, 1014 clinical trial, 212 preclinical/patented, and 1479 literature-reported targets via systematic review. These characteristic categories were classified into three distinct perspectives: molecular interaction/regulation, human system profile and cell-based expression variation. With the rapid progression of technology and concerted effort in drug discovery, TTD and other databases were highly expected to facilitate the explorations of druggability characteristics for the discovery and validation of innovative drug target. TTD is now freely accessible at: https://idrblab.org/ttd/.

TheMarker: a comprehensive database of therapeutic biomarkers

Distinct from the traditional diagnostic/prognostic biomarker (adopted as the indicator of disease state/process), the therapeutic biomarker (ThMAR) has emerged to be very crucial in the clinical development and clinical practice of all therapies. There are five types of ThMAR that have been found to play indispensable roles in various stages of drug discovery, such as: Pharmacodynamic Biomarker essential for guaranteeing the pharmacological effects of a therapy, Safety Biomarker critical for assessing the extent or likelihood of therapy-induced toxicity, Monitoring Biomarker indispensable for guiding clinical management by serially measuring patients' status, Predictive Biomarker crucial for maximizing the clinical outcome of a therapy for specific individuals, and Surrogate Endpoint fundamental for accelerating the approval of a therapy. However, these data of ThMARs has not been comprehensively described by any of the existing databases. Herein, a database, named 'TheMarker', was therefore constructed to (a) systematically offer all five types of ThMAR used at different stages of drug development, (b) comprehensively describe ThMAR information for the largest number of drugs among available databases, (c) extensively cover the widest disease classes by not just focusing on anticancer therapies. These data in TheMarker are expected to have great implication and significant impact on drug discovery and clinical practice, and it is freely accessible without any login requirement at: https://idrblab.org/themarker.

ADCdb: the database of antibody-drug conjugates

Antibody-drug conjugates (ADCs) are a class of innovative biopharmaceutical drugs, which, via their antibody (mAb) component, deliver and release their potent warhead (a.k.a. payload) at the disease site, thereby simultaneously improving the efficacy of delivered therapy and reducing its off-target toxicity. To design ADCs of promising efficacy, it is crucial to have the critical data of pharma-information and biological activities for each ADC. However, no such database has been constructed yet. In this study, a database named ADCdb focusing on providing ADC information (especially its pharma-information and biological activities) from multiple perspectives was thus developed. Particularly, a total of 6572 ADCs (359 approved by FDA or in clinical trial pipeline, 501 in preclinical test, 819 with in-vivo testing data, 1868 with cell line/target testing data, 3025 without in-vivo/cell line/target testing data) together with their explicit pharma-information was collected and provided. Moreover, a total of 9171 literature-reported activities were discovered, which were identified from diverse clinical trial pipelines, model organisms, patient/cell-derived xenograft models, etc. Due to the significance of ADCs and their relevant data, this new database was expected to attract broad interests from diverse research fields of current biopharmaceutical drug discovery. The ADCdb is now publicly accessible at: https://idrblab.org/adcdb/.

Computational Protein Design - Where it goes?

Proteins have been playing a critical role in the regulation of diverse biological processes related to human life. With the increasing demand, functional proteins are sparse in this immense sequence space. Therefore, protein design has become an important task in various fields, including medicine, food, energy, materials, etc. Directed evolution has recently led to significant achievements. Molecular modification of proteins through directed evolution technology has significantly advanced the fields of enzyme engineering, metabolic engineering, medicine, and beyond. However, it is impossible to identify desirable sequences from a large number of synthetic sequences alone. As a result, computational methods, including data-driven machine learning and physics-based molecular modeling, have been introduced to protein engineering to produce more functional proteins. This review focuses on recent advances in computational protein design, highlighting the applicability of different approaches as well as their limitations.

Mutexa: A Computational Ecosystem for Intelligent Protein Engineering

Protein engineering holds immense promise in shaping the future of biomedicine and biotechnology. This Review focuses on our ongoing development of Mutexa, a computational ecosystem designed to enable "intelligent protein engineering". In this vision, researchers will seamlessly acquire sequences of protein variants with desired functions as biocatalysts, therapeutic peptides, and diagnostic proteins through a finely-tuned computational machine, akin to Amazon Alexa's role as a versatile virtual assistant. The technical foundation of Mutexa has been established through the development of a database that combines and relates enzyme structures and their respective functions (e.g., IntEnzyDB), workflow software packages that enable high-throughput protein modeling (e.g., EnzyHTP and LassoHTP), and scoring functions that map the sequence-structure-function relationship of proteins (e.g., EnzyKR and DeepLasso). We will showcase the applications of these tools in benchmarking the convergence conditions of enzyme functional descriptors across mutants, investigating protein electrostatics and cavity distributions in SAM-dependent methyltransferases, and understanding the role of nonelectrostatic dynamic effects in enzyme catalysis. Finally, we will conclude by addressing the future steps and fundamental challenges in our endeavor to develop new Mutexa applications that assist the identification of beneficial mutants in protein engineering.

Targeting the omicron variant of SARS-CoV-2 with phytochemicals from Saudi medicinal plants: molecular docking combined with molecular dynamics investigations

The new health crises caused by SARS-CoV-2 have resulted in millions of deaths worldwide. First discovered in November 2021, the omicron variant is more transmissible and is able to evade the immune system better than other previously identified SARS-CoV-2 variants, leading to a spike in cases. Great efforts have been made to discover inhibitors against SARS-CoV-2. Main protease (Mpro) inhibitors are considered promising anti-SARS-CoV-2 agents. The U.S. FDA has issued an Emergency Use Authorization for ritonavir-boosted nirmatrelvir. Nirmatrelvir is the first orally bioavailable inhibitor of SARS-CoV-2 Mpro. There is an urgent need to monitor the mutations and solve the problem of resistance, especially omicron Mpro, which contains one mutation - P132H. In the present study, 132,57 phytochemicals from 80 medicinal plants grown in Saudi Arabia were docked into the active site of Mpro omicron variant. Free binding energies were also calculated. This led to the discovery of five phytochemicals that showed better docking scores than the bound ligand nirmatrelvir. In addition, these molecules exhibited favorable free binding energies. The stability of compounds 1-5 with the protein was studied using molecular dynamics (MD) simulations. These compounds showed acceptable ADMET properties. The results were compared with the wild type. These candidates could be envisioned as new hits against SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

Molecular insights and optimization strategies for the competitive binding of engineered ACE2 proteins: a multiple replica molecular dynamics study

The COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) continues to spread globally, and rapid viral evolution and the emergence of new variants pose challenges to pandemic control. During infection, the spike protein of SARS-CoV-2 interacts with the human ACE2 protein via its receptor binding domain (RBD), and it is known that engineered forms of ACE2 can compete with wild-type (WT) ACE2 for binding to inhibit infection. Here, we conducted multiple replica molecular dynamics (MRMD) simulations to study the mechanisms of the engineered ACE2 variants 3N39 and 3N94 and provide directions for optimization. Our findings reveal that engineered ACE2 is notably more efficacious in systems that show weaker binding to WT ACE2 (i.e., WT and BA.1 RBD), but also faces immune escape as the virus evolves. Moreover, by modifying residue types near the binding interface, engineered ACE2 alters the electrostatic potential distribution and reconfigures the hydrogen bonding network, which results in modified binding to the RBD. However, this structural rearrangement does not occur in all RBD variants. In addition, we identified potentially engineerable beneficial residues and potentially engineerable detrimental residues in both ACE2 and RBD. Functional conservation can thus enable the optimization of these residues and improve the binding competitiveness of engineered ACE2, which therefore provides additional immune escape prevention. Finally, we conclude that these findings have implications for understanding the mechanisms responsible for engineered ACE2 and can help us to develop engineered ACE2 proteins that show superior performance.

The Important Role of Transporter Structures in Drug Disposition, Efficacy, and Toxicity

The ATP-binding cassette (ABC) and solute carrier (SLC) transporters are critical determinants of drug disposition, clinical efficacy, and toxicity as they specifically mediate the influx and efflux of various substrates and drugs. ABC transporters can modulate the pharmacokinetics of many drugs via mediating the translocation of drugs across biologic membranes. SLC transporters are important drug targets involved in the uptake of a broad range of compounds across the membrane. However, high-resolution experimental structures have been reported for a very limited number of transporters, which limits the study of their physiologic functions. In this review, we collected structural information on ABC and SLC transporters and described the application of computational methods in structure prediction. Taking P-glycoprotein (ABCB1) and serotonin transporter (SLC6A4) as examples, we assessed the pivotal role of structure in transport mechanisms, details of ligand-receptor interactions, drug selectivity, the molecular mechanisms of drug-drug interactions, and differences caused by genetic polymorphisms. The data collected contributes toward safer and more effective pharmacological treatments. SIGNIFICANCE STATEMENT: The experimental structure of ATP-binding cassette and solute carrier transporters was collected, and the application of computational methods in structure prediction was described. P-glycoprotein and serotonin transporter were used as examples to reveal the pivotal role of structure in transport mechanisms, drug selectivity, the molecular mechanisms of drug-drug interactions, and differences caused by genetic polymorphisms.

Mapping synthetic binding proteins epitopes on diverse protein targets by protein structure prediction and protein-protein docking

Synthetic binding proteins (SBPs) are a class of artificial proteins engineered from privileged protein scaffolds, which can form highly specific molecular recognition interfaces with a variety of targets. Due to the characteristics of small size, high stability, and good tissue permeability, SBPs have important applications in biomedical research, disease diagnosis and treatment. However, knowledge of SBPs epitopes on the structures of target proteins is still limited, which hinder the development of novel SBPs. In this study, based on the currently available information of SBPs and their targets, 96 pairs of interacting proteins referring to 96 representative SBPs and 80 different targets, were systemically investigated using the state-of-the-art computational modeling techniques including AlphaFold2 protein structure prediction and Rosetta protein-protein docking. As a result, 71 out of the 96 pairs were successfully docked, of which 18, 33, and 20 pairs were defined as models with high, medium, and acceptable quality, respectively. In addition, the interface information was analyzed to decipher the interaction types driven SBPs and targets recognition. Overall, this work not only provides important structural information for understanding the mechanism of action of other SBPs with same protein scaffold, but also for aiding the rational protein engineering and to design of novel SBPs with biomedical applications.

Cancer-inspired genomics mapper model for the generation of synthetic DNA sequences with desired genomics signatures

Genome data are crucial in modern medicine, offering significant potential for diagnosis and treatment. Thanks to technological advancements, many millions of healthy and diseased genomes have already been sequenced; however, obtaining the most suitable data for a specific study, and specifically for validation studies, remains challenging with respect to scale and access. Therefore, in silico genomics sequence generators have been proposed as a possible solution. However, the current generators produce inferior data using mostly shallow (stochastic) connections, detected with limited computational complexity in the training data. This means they do not take the appropriate biological relations and constraints, that originally caused the observed connections, into consideration. To address this issue, we propose cancer-inspired genomics mapper model (CGMM), that combines genetic algorithm (GA) and deep learning (DL) methods to tackle this challenge. CGMM mimics processes that generate genetic variations and mutations to transform readily available control genomes into genomes with the desired phenotypes. We demonstrate that CGMM can generate synthetic genomes of selected phenotypes such as ancestry and cancer that are indistinguishable from real genomes of such phenotypes, based on unsupervised clustering. Our results show that CGMM outperforms four current state-of-the-art genomics generators on two different tasks, suggesting that CGMM will be suitable for a wide range of purposes in genomic medicine, especially for much-needed validation studies.

Rm-LR: A long-range-based deep learning model for predicting multiple types of RNA modifications

Recent research has highlighted the pivotal role of RNA post-transcriptional modifications in the regulation of RNA expression and function. Accurate identification of RNA modification sites is important for understanding RNA function. In this study, we propose a novel RNA modification prediction method, namely Rm-LR, which leverages a long-range-based deep learning approach to accurately predict multiple types of RNA modifications using RNA sequences only. Rm-LR incorporates two large-scale RNA language pre-trained models to capture discriminative sequential information and learn local important features, which are subsequently integrated through a bilinear attention network. Rm-LR supports a total of ten RNA modification types (m6A, m1A, m5C, m5U, m6Am, Ψ, Am, Cm, Gm, and Um) and significantly outperforms the state-of-the-art methods in terms of predictive capability on benchmark datasets. Experimental results show the effectiveness and superiority of Rm-LR in prediction of various RNA modifications, demonstrating the strong adaptability and robustness of our proposed model. We demonstrate that RNA language pretrained models enable to learn dense biological sequential representations from large-scale long-range RNA corpus, and meanwhile enhance the interpretability of the models. This work contributes to the development of accurate and reliable computational models for RNA modification prediction, providing insights into the complex landscape of RNA modifications.

Bioaffinity Nanoprobes for Foodborne Pathogen Sensing

Bioaffinity nanoprobes are a type of biosensor that utilize the specific binding properties of biological molecules, such as antibodies, enzymes, and nucleic acids, for the detection of foodborne pathogens. These probes serve as nanosensors and can provide highly specific and sensitive detection of pathogens in food samples, making them an attractive option for food safety testing. The advantages of bioaffinity nanoprobes include their ability to detect low levels of pathogens, rapid analysis time, and cost-effectiveness. However, limitations include the need for specialized equipment and the potential for cross-reactivity with other biological molecules. Current research efforts focus on optimizing the performance of bioaffinity probes and expanding their application in the food industry. This article discusses relevant analytical methods, such as surface plasmon resonance (SPR) analysis, Fluorescence Resonance Energy Transfer (FRET) measurements, circular dichroism, and flow cytometry, that are used to evaluate the efficacy of bioaffinity nanoprobes. Additionally, it discusses advances in the development and application of biosensors in monitoring foodborne pathogens.

Multi-state Model-Based Identification of Cryptic Allosteric Sites on Human Serotonin Transporter

Serotonin transporter (SERT) plays a fundamental role in taking the synaptic cleft serotonin back to the presynaptic neuron. The discovery of allosteric SERT modulators represents the next-generation medication for psychiatric disorders such as depression. Here, based on the cryo-EM structures of ibogaine in complex with SERT in distinct conformations, the multiple functional structures of the transporter bound to serotonin, including outward-open (OO_holo), outward-occluded (OC_holo), and inward-open (IO_holo and IO_holo'), were carefully characterized by induced-fit docking Gaussian-accelerated molecular dynamics (IFD-GaMD) simulation and the free-energy landscape analysis. Further MM/GBSA binding free energy, per-residue contribution, and molecular interaction fingerprint calculations revealed the interaction variations of serotonin with SERT in functional structures, which confirmed the allostery of SERT during serotonin reuptake. Moreover, five unique cryptic allosteric sites, which are druggable and capable of targeting by small molecules, were identified on the characterized multistate structures. These results provide structural and energetic information for the molecular mechanism of serotonin reuptake and will provide opportunities for the development of novel therapeutics based on the identified new allosteric sites on SERT.

Is BF.7 more infectious than other Omicron subtypes: Insights from structural and simulation studies of BF.7 spike RBD variant

Fear of a fresh infection wave and a global health issue in the ongoing COVID-19 pandemic have been rekindled by the appearance of two new novel variants BF.7 and BA.4/5 of Omicron lineages. Predictions of increased antibody evasion capabilities and transmissibility have been recognised in addition to the existing lineages (BA.1.1, BA.2, BA.2.12.1 and BA.3) as cause for worry. In comparison to Omicron, BA.4 and BF.7 share nine mutations in the spike protein, Leu371Phe, Thr376Ala, Asp405Asn, Arg408Ser, Ser446Gly, Leu452Arg, Phe486Val, Arg493Gln, Ser496Gly, whereas BF.7 contains an additional mutation, Arg346Thr, in the receptor binding domain (RBD) region. Due to the critical need for analysis and data on the BA.4 and BF.7 variants, we have computationally analyzed the interaction pattern between the Omicron, BA.4 and BF.7 RBD and angiotensin-converting enzyme 2 (ACE2) to determine the influence of these unique mutations on the structures, functions, and binding affinity of RBD towards ACE2. These analyses also allow to compare molecular models to previously reported data to evaluate the robustness of our methods for quick prediction of emerging future variants. The docking results reveal that BA.4 and BF.7 have particularly strong interactions with ACE2 when compared to Omicron, as shown by several parameters such as salt bridge, hydrogen bond, and non-bonded interactions. In addition, the estimations of binding free energy corroborated the findings further. BA.4 and BF.7 were found to bind to ACE2 with similar affinities (-72.14 and - 71.54 kcal/mol, respectively) and slightly stronger than Omicron (-70.04 kcal/mol). The differences in the binding pattern between the Omicron, BA.4 and BF.7 variant complexes indicated that the BA.4 and BF.7 RBD substitutions Asp405Asn, Ser446Gly, Leu452Arg, Phe486Val and Arg493Gln caused additional interactions with ACE2. In addition, normal mode analyses also indicate more stable conformations of BA.4 and BF.7 RBDs against human ACE2. Based on these structural and simulation analyses, we hypothesized that these changes may affect the binding affinity of BA.4 and BF.7 with ACE2.

Unbiased molecular dynamics simulation of a first-in-class small molecule inhibitor binds to oncostatin M

Small molecule inhibitors (SMIs) targeting oncostatin M (OSM) signaling pathway represent new therapeutics to combat cancer, inflammatory bowel disease (IBD) and CNS disease. Recently, the first-in-class SMI named SMI-10B that target OSM and block its interaction with receptor (OSMR) were reported. However, the binding pocket and interaction mode of the compound on OSM remain poorly understood, which hampering the rational design of SMIs that target OSM. Here, using SMI-10B as a probe, the multiple pockets on OSM for small molecules binding were extensively explored by unbiased molecular dynamics (MD) simulations. Then, the near-native structure of the complex was identified by molecular mechanics generalized Born surface area (MM/GBSA) binding energy funnel. Moreover, the binding stabilities of the protein-ligand complexes in near- and non-native conformations were verified by additional independent MD runs and absolute free energy perturbation (FEP) calculation. In summary, the unique feature of SMI-10B spontaneously binds to OSM characterized here not only provide detailed information for understanding the molecular mechanism of SMI-10B binding to OSM, but also will facilitate the rational design of novel and more potent SMIs to block OSM signaling.

A systematic review of state-of-the-art strategies for machine learning-based protein function prediction

New drug discovery is inseparable from the discovery of drug targets, and the vast majority of the known targets are proteins. At the same time, proteins are essential structural and functional elements of living cells necessary for the maintenance of all forms of life. Therefore, protein functions have become the focus of many pharmacological and biological studies. Traditional experimental techniques are no longer adequate for rapidly growing annotation of protein sequences, and approaches to protein function prediction using computational methods have emerged and flourished. A significant trend has been to use machine learning to achieve this goal. In this review, approaches to protein function prediction based on the sequence, structure, protein-protein interaction (PPI) networks, and fusion of multi-information sources are discussed. The current status of research on protein function prediction using machine learning is considered, and existing challenges and prominent breakthroughs are discussed to provide ideas and methods for future studies.

A systematic review on the state-of-the-art strategies for protein representation

The study of drug-target protein interaction is a key step in drug research. In recent years, machine learning techniques have become attractive for research, including drug research, due to their automated nature, predictive power, and expected efficiency. Protein representation is a key step in the study of drug-target protein interaction by machine learning, which plays a fundamental role in the ultimate accomplishment of accurate research. With the progress of machine learning, protein representation methods have gradually attracted attention and have consequently developed rapidly. Therefore, in this review, we systematically classify current protein representation methods, comprehensively review them, and discuss the latest advances of interest. According to the information extraction methods and information sources, these representation methods are generally divided into structure and sequence-based representation methods. Each primary class can be further divided into specific subcategories. As for the particular representation methods involve both traditional and the latest approaches. This review contains a comprehensive assessment of the various methods which researchers can use as a reference for their specific protein-related research requirements, including drug research.

Development of a Functional Nanobody Targeting Programmed Cell Death Protein-1 as Immune Checkpoint Inhibitor

BACKGROUND

Programmed cell death protein 1 (PD-1) is a membrane receptor that is expressed on the surface of various immune cells, such as T cells, B cells, monocytes, natural killer T cells, and dendritic cells. In cancer, the interaction between PD-1 and its ligand PD-L1 suppresses the activation and function of T lymphocytes, leading to the impairment and apoptosis of tumor-specific T cells. This mechanism allows cancer cells to evade the immune response and promotes tumor progression.

METHODS

Recombinant PD-1 protein was produced and used to immunize a camel. A nanobody library was generated from the camel's peripheral blood lymphocytes and screened for PD-1 binding. A specific nanobody (3PD9) was selected and characterized by affinity measurement, western blotting, and flow cytometry analysis. The ability of the selected nanobody to block the inhibitory signal of PD-1 in peripheral blood mononuclear cells (PBMCs) was evaluated by measuring the level of interleukin-2 (IL-2).

RESULTS

The selected nanobody showed high specificity and affinity for human PD-1. Western blot and flow cytometry analysis confirmed that 3PD9 could recognize and bind to human PD-1 on the cell surface. It was demonstrated that the level of IL-2 was significantly increased in PBMCs treated with 3PD9 compared to the control group, indicating that the nanobody could enhance the T cell response by disrupting the PD-1/PD-L1 interaction.

CONCLUSION

The results suggested that the anti-PD-1 nanobody could be a promising candidate for cancer immunotherapy.

Identification of potential inhibitor targeting KRAS mutation in Papillary Thyroid Carcinoma through molecular docking and dynamic simulation analysis

Thyroid cancer is a predominant form of endocrine malignancy, which destabilizes the metabolic rate of the body. The rapid increase in the incidence rate of thyroid cancer in recent years has aroused great concern to be investigated and diagnosed at an early stage. This study aimed to analyze the pathogenic mutations in thyroid cancer to identify their potential inhibitors for therapeutic targets. RAS genes are the most common oncogenes, which encode proteins that play an essential role in cell signaling and have been frequently mutated in different cancer types. The mutation in these genes causes abnormal cell growth and fails to respond to death signals. In this study, we identified the most significant mutations in the RAS genes; thus, the highly pathogenic mutations were curated from thyroid cancer patients and analyzed for their pathogenicity effect. The physicochemical analysis predicted mutation in wild-type KRAS protein had adapted negative charge on single base substitution of G12D that may easily cause loss of interactions and result in critical differences in the structure and function of the protein. Furthermore, the native KRAS protein was mutated and screened against a library of druggable compounds from the ZINC drug repository. The molecular docking analysis revealed that G12D mutant KRAS protein form best-docked complex with Naldemedine with the highest binding affinity. The dynamic simulation results further justified the stability of Naldemedine as a potential inhibitor with high efficiency in MMPBSA value of -45.4867 kcal/mol of being treated as a potential drug for papillary thyroid carcinoma. Further in vivo and in vitro validation of Naldemedine and its efficiency as a drug for the targeted pathogenic KRAS mutation is required.

In silico analysis of a novel protein in CAR T- cell therapy for the treatment of hematologic cancer through molecular modelling, docking, and dynamics approach

Cellular therapy has emerged as a key tool in the treatment of hematological malignancies. An advanced cell therapy known as chimeric antigen receptor T cell therapy (CAR T-cell therapy) has been approved by the United States Food and Drug Administration (FDA) as KYMRIAH by Novartis and YESCARTA by Gilead/Kite pharma in the year 2017. A chimeric receptor is composed of an extracellular antigen recognition site along with some co-stimulating and signalling domains. On the whole, it turns out to be one of the most potent receptors on T cells targeting a specific type of cancer cell based on its antigenic marker. CD19 CAR T-cell therapy is the first clinically approved therapy for lymphoma with remarkable results in complete remission of B cell lymphoblastic leukemia up to 90%. The high rate of effectiveness of the CAR T-cell therapy against B-ALL justifies the investigation and application of this therapy for fatal diseases like all types of hematological malignancies. The most critical aspect of chimeric receptor therapy is designing and building an artificial receptor that is specific to a given type of cancer. For this reason, the in silico technique is an appropriate model to investigate the integrity and effectiveness of the engineered chimeric receptor prior to commencing in vitro experiments followed by clinical trials. This computerized experimental study aids in predicting the molecular mechanism of chimeric protein and how it interacts with both ligands. We have anticipated various features of the chimeric protein in terms of qualitative analysis (structure, protein modelling, physiological properties) and functional analysis (antigenicity, allergenicity, its receptor-ligand binding ability, involving signalling pathways). Furthermore, the reliability and validation of the binding mode of the chimeric protein against receptors were performed through a complex molecular dynamics simulation for a 100 ns timeframe in an aqueous environment. The obtained simulation study showed that CD30 was a better approachable marker as compared to CD20 due to its better binding energy score and also binding conformations stability.

DRESIS: the first comprehensive landscape of drug resistance information

Widespread drug resistance has become the key issue in global healthcare. Extensive efforts have been made to reveal not only diverse diseases experiencing drug resistance, but also the six distinct types of molecular mechanisms underlying this resistance. A database that describes a comprehensive list of diseases with drug resistance (not just cancers/infections) and all types of resistance mechanisms is now urgently needed. However, no such database has been available to date. In this study, a comprehensive database describing drug resistance information named 'DRESIS' was therefore developed. It was introduced to (i) systematically provide, for the first time, all existing types of molecular mechanisms underlying drug resistance, (ii) extensively cover the widest range of diseases among all existing databases and (iii) explicitly describe the clinically/experimentally verified resistance data for the largest number of drugs. Since drug resistance has become an ever-increasing clinical issue, DRESIS is expected to have great implications for future new drug discovery and clinical treatment optimization. It is now publicly accessible without any login requirement at: https://idrblab.org/dresis/.

DrugMAP: molecular atlas and pharma-information of all drugs

The efficacy and safety of drugs are widely known to be determined by their interactions with multiple molecules of pharmacological importance, and it is therefore essential to systematically depict the molecular atlas and pharma-information of studied drugs. However, our understanding of such information is neither comprehensive nor precise, which necessitates the construction of a new database providing a network containing a large number of drugs and their interacting molecules. Here, a new database describing the molecular atlas and pharma-information of drugs (DrugMAP) was therefore constructed. It provides a comprehensive list of interacting molecules for >30 000 drugs/drug candidates, gives the differential expression patterns for >5000 interacting molecules among different disease sites, ADME (absorption, distribution, metabolism and excretion)-relevant organs and physiological tissues, and weaves a comprehensive and precise network containing >200 000 interactions among drugs and molecules. With the great efforts made to clarify the complex mechanism underlying drug pharmacokinetics and pharmacodynamics and rapidly emerging interests in artificial intelligence (AI)-based network analyses, DrugMAP is expected to become an indispensable supplement to existing databases to facilitate drug discovery. It is now fully and freely accessible at: https://idrblab.org/drugmap/.

CovInter: interaction data between coronavirus RNAs and host proteins

Coronavirus has brought about three massive outbreaks in the past two decades. Each step of its life cycle invariably depends on the interactions among virus and host molecules. The interaction between virus RNA and host protein (IVRHP) is unique compared to other virus-host molecular interactions and represents not only an attempt by viruses to promote their translation/replication, but also the host's endeavor to combat viral pathogenicity. In other words, there is an urgent need to develop a database for providing such IVRHP data. In this study, a new database was therefore constructed to describe the interactions between coronavirus RNAs and host proteins (CovInter). This database is unique in (a) unambiguously characterizing the interactions between virus RNA and host protein, (b) comprehensively providing experimentally validated biological function for hundreds of host proteins key in viral infection and (c) systematically quantifying the differential expression patterns (before and after infection) of these key proteins. Given the devastating and persistent threat of coronaviruses, CovInter is highly expected to fill the gap in the whole process of the 'molecular arms race' between viruses and their hosts, which will then aid in the discovery of new antiviral therapies. It's now free and publicly accessible at: https://idrblab.org/covinter/.

Molecular dynamics simulations highlight the altered binding landscape at the spike-ACE2 interface between the Delta and Omicron variants compared to the SARS-CoV-2 original strain

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) B.1.1.529 variant (Omicron), represents a significant deviation in genetic makeup and function compared to previous variants. Following the BA.1 sublineage, the BA.2 and BA.3 Omicron subvariants became dominant, and currently the BA.4 and BA.5, which are quite distinct variants, have emerged. Using molecular dynamics simulations, we investigated the binding characteristics of the Delta and Omicron (BA.1) variants in comparison to wild-type (WT) at the interface of the spike protein receptor binding domain (RBD) and human angiotensin converting enzyme-2 (ACE2) ectodomain. The primary aim was to compare our molecular modelling systems with previously published observations, to determine the robustness of our approach for rapid prediction of emerging future variants. Delta and Omicron were found to bind to ACE2 with similar affinities (−39.4 and −43.3 kcal/mol, respectively) and stronger than WT (−33.5 kcal/mol). In line with previously published observations, the energy contributions of the non-mutated residues at the interface were largely retained between WT and the variants, with F456, F486, and Y489 having the strongest energy contributions to ACE2 binding. Further, residues N440K, Q498R, and N501Y were predicted to be energetically favourable in Omicron. In contrast to Omicron, which had the E484A and K417N mutations, intermolecular bonds were detected for the residue pairs E484:K31 and K417:D30 in WT and Delta, in accordance with previously published findings. Overall, our simplified molecular modelling approach represents a step towards predictive model systems for rapidly analysing arising variants of concern.

M6AREG: m6A-centered regulation of disease development and drug response

As the most prevalent internal modification in eukaryotic RNAs, N6-methyladenosine (m6A) has been discovered to play an essential role in cellular proliferation, metabolic homeostasis, embryonic development, etc. With the rapid accumulation of research interest in m6A, its crucial roles in the regulations of disease development and drug response are gaining more and more attention. Thus, a database offering such valuable data on m6A-centered regulation is greatly needed; however, no such database is as yet available. Herein, a new database named 'M6AREG' is developed to (i) systematically cover, for the first time, data on the effects of m6A-centered regulation on both disease development and drug response, (ii) explicitly describe the molecular mechanism underlying each type of regulation and (iii) fully reference the collected data by cross-linking to existing databases. Since the accumulated data are valuable for researchers in diverse disciplines (such as pathology and pathophysiology, clinical laboratory diagnostics, medicinal biochemistry and drug design), M6AREG is expected to have many implications for the future conduct of m6A-based regulation studies. It is currently accessible by all users at: https://idrblab.org/m6areg/.

Differential performance of RoseTTAFold in antibody modeling

Antibodies are essential to life, and knowing their structures can facilitate the understanding of antibody-antigen recognition mechanisms. Precise antibody structure prediction has been a core challenge for a prolonged period, especially the accuracy of H3 loop prediction. Despite recent progress, existing methods cannot achieve atomic accuracy, especially when the homologous structures required for these methods are not available. Recently, RoseTTAFold, a deep learning-based algorithm, has shown remarkable breakthroughs in predicting the 3D structures of proteins. To assess the antibody modeling ability of RoseTTAFold, we first retrieved the sequences of 30 antibodies as the test set and used RoseTTAFold to model their 3D structures. We then compared the models constructed by RoseTTAFold with those of SWISS-MODEL in a different way, in which we stratified Global Model Quality Estimate (GMQE) into three different ranges. The results indicated that RoseTTAFold could achieve results similar to SWISS-MODEL in modeling most CDR loops, especially the templates with a GMQE score under 0.8. In addition, we also compared the structures modeled by RoseTTAFold, SWISS-MODEL and ABodyBuilder. In brief, RoseTTAFold could accurately predict 3D structures of antibodies, but its accuracy was not as good as the other two methods. However, RoseTTAFold exhibited better accuracy for modeling H3 loop than ABodyBuilder and was comparable to SWISS-MODEL. Finally, we discussed the limitations and potential improvements of the current RoseTTAFold, which may help to further the accuracy of RoseTTAFold's antibody modeling.

Investigating the structure-activity relationship of marine polycyclic batzelladine alkaloids as promising inhibitors for SARS-CoV-2 main protease (Mpro)

Over a span of two years ago, since the emergence of the first case of the novel coronavirus (SARS-CoV-2) in China, the pandemic has crossed borders causing serious health emergencies, immense economic crisis and impacting the daily life worldwide. Despite the discovery of numerous forms of precautionary vaccines along with other recently approved orally available drugs, yet effective antiviral therapeutics are necessarily needed to hunt this virus and its variants. Historically, naturally occurring chemicals have always been considered the primary source of beneficial medications. Considering the SARS-CoV-2 main protease (M^pro) as the duplicate key element of the viral cycle and its main target, in this paper, an extensive virtual screening for a focused chemical library of 15 batzelladine marine alkaloids, was virtually examined against SARS-CoV-2 main protease (M^pro) using an integrated set of modern computational tools including molecular docking (MDock), molecule dynamic (MD) simulations and structure-activity relationships (SARs) as well. The molecular docking predictions had disclosed four promising compounds including batzelladines H–I (8–9) and batzelladines F-G (6–7), respectively according to their prominent ligand-protein energy scores and relevant binding affinities with the (M^pro) pocket residues. The best two chemical hits, batzelladines H–I (8–9) were further investigated thermodynamically though studying their MD simulations at 100 ns, where they showed excellent stability within the accommodated (M^pro) pocket. Moreover, SARs studies imply the crucial roles of the fused tricyclic guanidinic moieties, its degree of unsaturation, position of the N–OH functionality and the length of the side chain as a spacer linking between two active sites, which disclosed fundamental structural and pharmacophoric features for efficient protein-ligand interaction. Such interesting findings are greatly highlighting further in vitro/vivo examinations regarding those marine natural products (MNPs) and their synthetic equivalents as promising antivirals.

Exploration and validation of Taraxacum mongolicum anti-cancer effect

Taraxacum mongolicum gained a lot of concern and was applied in 93 formulas in China due to its fame as a traditional Chinese medicine. The earliest recorded application of Taraxacum mongolicum was traced back to the Han dynasty. Generations of doctors boosted the usage and enriched the pharmacological mechanism. Clinical application of the Taraxacum mongolicum is flourishing as it treats multiple diseases. This study aims to explore the anti-cancer effect, retrieve the active ingredients and screen the key targets of Taraxacum mongolicum in cancer therapy. We collected and evaluated 10 key active compounds to investigate the anti-cancer effect via 69 significant targets and a variety of biological processes and pathways. Gene Ontology (GO) enrichment analysis uncovered targets associated with protein phosphorylation, cell proliferation and apoptotic processes via regulation of kinases, ATP and enzyme binding activities. Half of the top 20 enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were directly involved in cancer. Based on standard selection criteria, seven hub targets were obtained. These targets functioned through distinct patterns and pathways in realizing the anti-cancer effect. Molecular docking was conducted to validate the potential combination between compounds and hub targets to explore the pharmacological mechanism of key compounds in Taraxacum mongolicum against cancer. In summary, our findings indicate that the famous and widely used Chinese herb, Taraxacum mongolicum, shows good anti-cancer effect through its active compounds, targeted genes, and multiple involved biological processes. The results may provide a theoretical basis for subsequent experimental validation and drug development of Taraxacum mongolicum extract against cancer.

ConSIG: consistent discovery of molecular signature from OMIC data

The discovery of proper molecular signature from OMIC data is indispensable for determining biological state, physiological condition, disease etiology, and therapeutic response. However, the identified signature is reported to be highly inconsistent, and there is little overlap among the signatures identified from different biological datasets. Such inconsistency raises doubts about the reliability of reported signatures and significantly hampers its biological and clinical applications. Herein, an online tool, ConSIG, was constructed to realize consistent discovery of gene/protein signature from any uploaded transcriptomic/proteomic data. This tool is unique in a) integrating a novel strategy capable of significantly enhancing the consistency of signature discovery, b) determining the optimal signature by collective assessment, and c) confirming the biological relevance by enriching the disease/gene ontology. With the increasingly accumulated concerns about signature consistency and biological relevance, this online tool is expected to be used as an essential complement to other existing tools for OMIC-based signature discovery. ConSIG is freely accessible to all users without login requirement at https://idrblab.org/consig/.

A multidrug efflux protein in Mycobacterium tuberculosis; tap as a potential drug target for drug repurposing

Tuberculosis (TB) is a serious communicative disease caused by Mycobacterium tuberculosis. Although there are vaccines and drugs available to treat the disease, they are not efficient, moreover, multidrug-resistant TB (MDR-TB) become a major hurdle in its therapy. These MDR strains utilize the multidrug efflux pump as a decisive weapon to fight against antitubercular drugs. Tap membrane protein was observed as a crucial multidrug efflux pump in M. tuberculosis and its critical implication in MDR-MTB development makes it an effective drug target. In the present study, we have utilized various in silico approaches to predict the applicability of FDA-approved ion channel inhibitors and blockers as therapeutic leads against Tuberculosis by targeting multidrug efflux protein; Tap in MTB. Tap protein structure is predicted by Phyre2 server followed by model refinement, validation, physio-chemical catheterization and target prediction. Further, the interaction between Tap protein and ligands were analysed by molecular docking and MD simulation run of 100 ns. Based on implication and compatibility, 18 FDA-approved ion channel inhibitors and blockers are selected as a ligand against the Tap protein and eventually observed five ligands; Glimepiride, Flecainide, Flupiritine, Nimodipine and Amlodipine as promising compounds which have exhibited the significant stable interaction with Tap protein and are proposed to modulate or interfere with its activity. These compounds illustrated the substantial docking score and total binding enthalpy more than -7 kcal/mol and -42 kcal/mol respectively which implies that the selected FDA-approved compounds can spontaneously interact with the Tap protein to modulate its function. This study proposed Tap protein as a prominent drug target in MTB and investigated compounds that show considerable interaction with the Tap protein as potential therapeutic molecules. These interactions may lead to modulating or inhibit the activity of drug efflux protein thereby making MTB susceptible to antitubercular drugs.

Biological activities of drug inactive ingredients

In a drug formulation (DFM), the major components by mass are not Active Pharmaceutical Ingredient (API) but rather Drug Inactive Ingredients (DIGs). DIGs can reach much higher concentrations than that achieved by API, which raises great concerns about their clinical toxicities. Therefore, the biological activities of DIG on physiologically relevant target are widely demanded by both clinical investigation and pharmaceutical industry. However, such activity data are not available in any existing pharmaceutical knowledge base, and their potentials in predicting the DIG-target interaction have not been evaluated yet. In this study, the comprehensive assessment and analysis on the biological activities of DIGs were therefore conducted. First, the largest number of DIGs and DFMs were systematically curated and confirmed based on all drugs approved by US Food and Drug Administration. Second, comprehensive activities for both DIGs and DFMs were provided for the first time to pharmaceutical community. Third, the biological targets of each DIG and formulation were fully referenced to available databases that described their pharmaceutical/biological characteristics. Finally, a variety of popular artificial intelligence techniques were used to assess the predictive potential of DIGs' activity data, which was the first evaluation on the possibility to predict DIG's activity. As the activities of DIGs are critical for current pharmaceutical studies, this work is expected to have significant implications for the future practice of drug discovery and precision medicine.

Identification of Anti-TNFα VNAR Single Domain Antibodies from Whitespotted Bambooshark (Chiloscyllium plagiosum)

Tumor necrosis factor α (TNFα), an important clinical testing factor and drug target, can trigger serious autoimmune diseases and inflammation. Thus, the TNFα antibodies have great potential application in diagnostics and therapy fields. The variable binding domain of IgNAR (VNAR), the shark single domain antibody, has some excellent advantages in terms of size, solubility, and thermal and chemical stability, making them an ideal alternative to conventional antibodies. This study aims to obtain VNARs that are specific for mouse TNF (mTNF) from whitespotted bamboosharks. After immunization of whitespotted bamboosharks, the peripheral blood leukocytes (PBLs) were isolated from the sharks, then the VNAR phage display library was constructed. Through phage display panning against mTNFα, positive clones were validated through ELISA assay. The affinity of the VNAR and mTNFα was measured using ELISA and Bio-Layer Interferometry. The binding affinity of 3B11 VNAR reached 16.7 nM. Interestingly, one new type of VNAR targeting mTNF was identified that does not belong to any known VNAR type. To understand the binding mechanism of VNARs to mTNFα, the models of VNARs-mTNFα complexes were predicted by computational modeling combining HawkDock and RosettaDock. Our results showed that four VNARs’ epitopes overlapped in part with that of mTNFR. Furthermore, the ELISA assay shows that the 3B11 potently inhibited mTNFα binding to mTNFR. This study may provide the basis for the TNFα blockers and diagnostics applications.

RGCNCDA: Relational graph convolutional network improves circRNA-disease association prediction by incorporating microRNAs

Recently, a large number of studies have indicated that circRNAs with covalently closed loops play important roles in biological processes and have potential as diagnostic biomarkers. Therefore, research on the circRNA-disease relationship is helpful in disease diagnosis and treatment. However, traditional biological verification methods require considerable labor and time costs. In this paper, we propose a new computational method (RGCNCDA) to predict circRNA-disease associations based on relational graph convolutional networks (R-GCNs). The method first integrates the circRNA similarity network, miRNA similarity network, disease similarity network and association networks among them to construct a global heterogeneous network. Then, it employs the random walk with restart (RWR) and principal component analysis (PCA) models to learn low-dimensional and high-order information from the global heterogeneous network as the topological features. Finally, a prediction model based on an R-GCN encoder and a DistMult decoder is built to predict the potential disease-associated circRNA. The predicted results demonstrate that RGCNCDA performs significantly better than the other six state-of-the-art methods in a 5-fold cross validation. Furthermore, the case study illustrates that RGCNCDA can effectively discover potential circRNA-disease associations.

PFmulDL: a novel strategy enabling multi-class and multi-label protein function annotation by integrating diverse deep learning methods

Bioinformatic annotation of protein function is essential but extremely sophisticated, which asks for extensive efforts to develop effective prediction method. However, the existing methods tend to amplify the representativeness of the families with large number of proteins by misclassifying the proteins in the families with small number of proteins. That is to say, the ability of the existing methods to annotate proteins in the 'rare classes' remains limited. Herein, a new protein function annotation strategy, PFmulDL, integrating multiple deep learning methods, was thus constructed. First, the recurrent neural network was integrated, for the first time, with the convolutional neural network to facilitate the function annotation. Second, a transfer learning method was introduced to the model construction for further improving the prediction performances. Third, based on the latest data of Gene Ontology, the newly constructed model could annotate the largest number of protein families comparing with the existing methods. Finally, this newly constructed model was found capable of significantly elevating the prediction performance for the 'rare classes' without sacrificing that for the 'major classes'. All in all, due to the emerging requirements on improving the prediction performance for the proteins in 'rare classes', this new strategy would become an essential complement to the existing methods for protein function prediction. All the models and source codes are freely available and open to all users at: https://github.com/idrblab/PFmulDL.

Computational studies on the interaction of SARS-CoV-2 Omicron SGp RBD with human receptor ACE2, limonin and glycyrrhizic acid

On November 24, 2021, the SARS-CoV-2 Omicron variant (B.1.1.529) was first identified in South Africa. The World Health Organization (WHO) declared the Omicron as a variant of concern (VoC) because of the unexpected and large numbers of mutations occurred in the genome, higher viral transmission and immune evasions. The present study was performed to explore the interactions of SARS-CoV-2 spike glycoprotein receptor-binding domain (SGp RBD) of the three variants (Omicron, Delta, and WT) with the receptor hACE2. The structural changes occurred in Omicron due to the mutations at key positions improved the ability to mediate SARS-CoV-2 viral infection compared to other VoCs. The phytochemicals limonin and glycyrrhizic acid were docked with the SGp RBD of the variants WT, Delta and Omicron. The computed dock score revealed that limonin and glycyrrhizic acid binds effectively at the SGp RBD of all three variants, and showed almost similar binding affinity at the binding interface of ACE2. Therefore, despite the multiple mutations occurred in Omicron and its viral transmission is comparatively high, the computed binding affinity of the phytochemicals limonin and glycyrrhizic acid supported that the traditional medicines can be useful in formulating adjuvant therapies to fight against the SARS-CoV-2 Omicron.

POSREG: proteomic signature discovered by simultaneously optimizing its reproducibility and generalizability

Mass spectrometry-based proteomic technique has become indispensable in current exploration of complex and dynamic biological processes. Instrument development has largely ensured the effective production of proteomic data, which necessitates commensurate advances in statistical framework to discover the optimal proteomic signature. Current framework mainly emphasizes the generalizability of the identified signature in predicting the independent data but neglects the reproducibility among signatures identified from independently repeated trials on different sub-dataset. These problems seriously restricted the wide application of the proteomic technique in molecular biology and other related directions. Thus, it is crucial to enable the generalizable and reproducible discovery of the proteomic signature with the subsequent indication of phenotype association. However, no such tool has been developed and available yet. Herein, an online tool, POSREG, was therefore constructed to identify the optimal signature for a set of proteomic data. It works by (i) identifying the proteomic signature of good reproducibility and aggregating them to ensemble feature ranking by ensemble learning, (ii) assessing the generalizability of ensemble feature ranking to acquire the optimal signature and (iii) indicating the phenotype association of discovered signature. POSREG is unique in its capacity of discovering the proteomic signature by simultaneously optimizing its reproducibility and generalizability. It is now accessible free of charge without any registration or login requirement at https://idrblab.org/posreg/.

Molecular Mechanism for the Allosteric Inhibition of the Human Serotonin Transporter by Antidepressant Escitalopram

Human serotine transporter (hSERT) is one of the most influential drug targets, and its allosteric modulators (e.g., escitalopram) have emerged to be the next-generation medication for psychiatric disorders. However, the molecular mechanism underlying the allosteric modulation of hSERT is still elusive. Here, the simulation strategies of conventional (cMD) and steered (SMD) molecular dynamics were applied to investigate this molecular mechanism from distinct perspectives. First, cMD simulations revealed that escitalopram's binding to hSERT's allosteric site simultaneously enhanced its binding to the orthosteric site. Then, SMD simulation identified that the occupation of hSERT's allosteric site by escitalopram could also block its dissociation from the orthosteric site. Finally, by comparing the simulated structures of two hSERT-escitalopram complexes with and without allosteric modulation, a new conformational coupling between an extracellular (Arg104-Glu494) and an intracellular (Lys490-Glu494) salt bridge was identified. In summary, this study explored the mechanism underlying the allosteric modulation of hSERT by collectively applying two MD simulation strategies, which could facilitate our understanding of the allosteric modulations of not only hSERT but also other clinically important therapeutic targets.

Recent advances in biological detection with rolling circle amplification: design strategy, biosensing mechanism, and practical applications

Rolling circle amplification (RCA) is a simple and isothermal DNA amplification technique that is used to generate thousands of repeating DNA sequences using circular templates under the catalysis of DNA polymerase.