Protein structure prediction and structural genomics

Based on hydrogen and disulfide-mediated bonds, l-lysine and l-arginine enhanced the gel properties of low-salt mixed shrimp surimi (Antarctic krill and Pacific white shrimp)

This study delved into the effects of l-lysine (Lys) and l-arginine (Arg) on the gel properties and intermolecular interactions of low-salt (NaCl, 1 g/100 g) mixed shrimp surimi (Antarctic krill and Pacific white shrimp). The addition of Lys and Arg improved the gel strength and water holding capacity of low-salt gels, which were superior to the properties of STPP and high-salt (NaCl, 2.25 g/100 g) gels. These results can be attributed to the role of Lys and Arg in enhancing hydrogen and disulfide bonds within the low-salt gel system, promoting the solubilization of myofibrillar proteins (MP) and consequently increasing the number of MP molecules participating in gel formation. Antarctic krill MP did not show gel-forming ability and exerted a diluting effect on low-salt mixed shrimp surimi gels. Molecular docking analysis indicated the stable binding of Lys and Arg to myosin.

Effect of chlorpyrifos-exposure on the expression levels of CYP genes in Daphnia magna and examination of a possibility that an up-regulated clan 3 CYP, CYP360A8, reacts with pesticides

Daphnia magna is a test organism used for ecological risk assessments of pesticides, but little is known about the expression levels of cytochrome P450s (CYP)s and their changes after pesticide exposure in the less than 24-h-olds used for ecotoxicity tests. In this study, D. magna juveniles were exposed to 0.2 μg/L of chlorpyrifos under the conditions for acute immobilization test as specified by the OECD test guideline for 24 h, and then the gene expression was compared between the control and chlorpyrifos-exposure groups by RNA-sequencing analysis, with a focus on CYP genes. Among 38 CYP genes expressed in the control group, seven were significantly up-regulated while two were significantly down-regulated in the chlorpyrifos-exposure group. Although the sublethal concentration of chlorpyrifos did not change their expression levels so drastically (0.8 < fold change < 2.6), CY360A8 of D. magna (DmCYP360A8), which had been proposed to be responsible for metabolism of xenobiotics, was abundantly expressed in controls yet up-regulated by chlorpyrifos. Therefore, homology modeling of DmCYP360A8 was performed based on the amino acid sequence, and then molecular docking simulations with the insecticides that were indicated to be metabolized by CYPs in D. magna were conducted. The results indicated that DmCYP360A8 could contribute to the metabolism of diazinon and chlorfenapyr but not chlorpyrifos. These findings suggest that chlorpyrifos is probably detoxified by other CYP(s) including up-regulated and/or constitutively expressed one(s).

Site-specific mutagenesis screening in KRASG12D mutant library to uncover resistance mechanisms to KRASG12D inhibitors

KRAS plays a crucial role in regulating cell survival and proliferation and is one of the most commonly mutated oncogenes in human cancers. The novel KRASG12D inhibitor, MRTX1133, demonstrates promising antitumor efficacy in vitro and in vivo. However, the development of acquired resistance in treated patients presents a considerable challenge to sustained therapeutic effectiveness. In response to this challenge, we conducted site-specific mutagenesis screening to identify potential secondary mutations that could induce resistance to MRTX1133. We screened a range of KRASG12D variants harboring potential secondary mutations, and 44 representative variants were selected for in-depth validation of the pooled screening outcomes. We identified eight variants (G12D with V9E, V9W, V9Q, G13P, T58Y, R68G, Y96W, and Q99L) that exhibited substantial resistance, with V9W showing notable resistance, and downstream signaling analyses and structural modeling were conducted. We observed that secondary mutations in KRASG12D can lead to acquired resistance to MRTX1133 and BI-2865, a novel pan-KRAS inhibitor, in human cancer cell lines. This evidence is critical for devising new strategies to counteract resistance mechanisms and, ultimately, enhance treatment outcomes in patients with KRASG12D-mutant cancers.

Structure and function of engineered stromal cell-derived factor-1α

To design biologically active, collagen-based scaffolds for bone tissue engineering, we have synthesized chimeric proteins consisting of stromal cell-derived factor-1α (SDF) and the von Willebrand factor A3 collagen-binding domain (CBD). The chimeric proteins were used to evaluate the effect of domain linkage and its order on the structure and function of the SDF and CBD. The structure of the chimeric proteins was analyzed by circular dichroism spectroscopy, while functional analysis was performed by a cell migration assay for the SDF domain and a collagen-binding assay for the CBD domain. Furthermore, computational structural prediction was conducted for the chimeric proteins to examine the consistency with the results of structural and functional analyses. Our structural and functional analyses as well as structural prediction revealed that linking two domains can affect their functions. However, their order had minor effects on the three-dimensional structure of CBD and SDF in the chimeric proteins.

Computational Workflow for Refining AlphaFold Models in Drug Design Using Kinetic and Thermodynamic Binding Calculations: A Case Study for the Unresolved Inactive Human Adenosine A3 Receptor

A structure-based drug design pipeline that considers both thermodynamic and kinetic binding data of ligands against a receptor will enable the computational design of improved drug molecules. For unresolved GPCR-ligand complexes, a workflow that can apply both thermodynamic and kinetic binding data in combination with alpha-fold (AF)-derived or other homology models and experimentally resolved binding modes of relevant ligands in GPCR-homologs needs to be tested. Here, as test case, we studied a congeneric set of ligands that bind to a structurally unresolved G protein-coupled receptor (GPCR), the inactive human adenosine A3 receptor (hA3R). We tested three available homology models from which two have been generated from experimental structures of hA1R or hA2AR and one model was a multistate alphafold 2 (AF2)-derived model. We applied alchemical calculations with thermodynamic integration coupled with molecular dynamics (TI/MD) simulations to calculate the experimental relative binding free energies and residence time (τ)-random accelerated MD (τ-RAMD) simulations to calculate the relative residence times (RTs) for antagonists. While the TI/MD calculations produced, for the three homology models, good Pearson correlation coefficients, correspondingly, r = 0.74, 0.62, and 0.67 and mean unsigned error (mue) values of 0.94, 1.31, and 0.81 kcal mol-1, the τ-RAMD method showed r = 0.92 and 0.52 for the first two models but failed to produce accurate results for the multistate AF2-derived model. With subsequent optimization of the AF2-derived model by reorientation of the side chain of R1735.34 located in the extracellular loop 2 (EL2) that blocked ligand's unbinding, the computational model showed r = 0.84 for kinetic data and improved performance for thermodynamic data (r = 0.81, mue = 0.56 kcal mol-1). Overall, after refining the multistate AF2 model with physics-based tools, we were able to show a strong correlation between predicted and experimental ligand relative residence times and affinities, achieving a level of accuracy comparable to an experimental structure. The computational workflow used can be applied to other receptors, helping to rank candidate drugs in a congeneric series and enabling the prioritization of leads with stronger binding affinities and longer residence times.

Molecular Modeling Techniques and In-Silico Drug Discovery

Molecular modeling is the technique to determine the overall structure of an unknown molecule, be it a small one or a macromolecule. The technique encompasses the method of screening ligand libraries for the development of new candidate drug molecules. All these aspects have become an essential topic of research. This field is truly interdisciplinary and finds its applications in almost all fields of life science research. In this chapter, an overview of the protocol associated with molecular modeling techniques will be discussed.

Characterization of terminal flowering cowpea (Vigna unguiculata (L.) Walp.) mutants obtained by induced mutagenesis digs out the loss-of-function of phosphatidylethanolamine-binding protein

Cowpea (Vigna unguiculata (L.) Walp) is one of the major food legume crops grown extensively in arid and semi-arid regions of the world. The determinate habit of cowpea has many advantages over the indeterminate and is well adapted to modern farming systems. Mutation breeding is an active research area to develop the determinate habit of cowpea. The present study aimed to develop new determinate habit mutants with terminal flowering (TFL) in locally well-adapted genetic backgrounds. Consequently, the seeds of popular cowpea cv P152 were irradiated with doses of gamma rays (200, 250, and, 300 Gy), and the M1 populations were grown. The M2 populations were produced from the M1 progenies and selected determinate mutants (TFLCM-1 and TFLCM-2) from the M2 generation (200 Gy) were forwarded up to the M5 generation to characterize the mutants and simultaneously they were crossed with P152 to develop a MutMap population. In the M5 generation, determinate mutants (80-81 days) were characterized by evaluating the TFL growth habit, longer peduncles (30.75-31.45 cm), erect pods (160°- 200°), number of pods per cluster (4-5 nos.), and early maturity. Further, sequencing analysis of the VuTFL1 gene in the determinate mutants and MutMap population revealed a single nucleotide transversion (A-T at 1196 bp) in the fourth exon and asparagine (N) to tyrosine (Y) amino acid change at the 143rd position of phosphatidylethanolamine-binding protein (PEBP). Notably, the loss of function PEPB with a higher confidence level modification of anti-parallel beta-sheets and destabilization of the protein secondary structure was observed in the mutant lines. Quantitative real-time PCR (qRT-PCR) analysis showed that the VuTFL1 gene was downregulated at the flowering stage in TFL mutants. Collectively, the insights garnered from this study affirm the effectiveness of induced mutation in modifying the plant's ideotype. The TFL mutants developed during this investigation have the potential to serve as a valuable resource for fostering determinate traits in future cowpea breeding programs and pave the way for mechanical harvesting.

Biochemical characterisation of a PL24 ulvan lyase from seaweed-associated Vibrio sp. FNV38

Ulvan is a green macroalgal cell wall polysaccharide that has tremendous potential for valorisation due to its unique composition of sulphated rhamnose, glucuronic acid, iduronic acid and xylose. Several potential applications such as production of biofuels, bioplastics and other value-added products necessitate the breakdown of the polysaccharide to oligomers or monomers. Research on ulvan saccharifying enzymes has been continually increasing over the last decade, with the increasing focus on valorisation of seaweed biomass for a biobased economy. Lyases are the first of several enzymes that are involved in saccharifying the polysaccharide and several ulvan lyases have been structurally and biochemically characterised to enable their effective use in the valorisation processes. This study investigates the whole genome of Vibrio sp. FNV38, an ulvan metabolising organism and biochemical characteristics of a PL24 ulvan lyase that it possesses. The genome of Vibrio sp. FNV38 has a diverse CAZy profile with several genes involved in the metabolism of ulvan, cellulose, agar, and alginate. The enzyme exhibits optimal activity at pH 8.5 in 100 mM Tris-HCl buffer and 30 °C. However, its thermal stability is poor with significant loss of activity after 2 h of incubation at temperatures above 25 °C. Breakdown product analysis reveals that the enzyme depolymerised the polysaccharide predominantly to disaccharides and tetrasaccharides.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10811-023-03136-3.

Computational approaches for modeling and structural design of biological systems: A comprehensive review

The convergence of biology and computational science has ushered in a revolutionary era, revolutionizing our understanding of biological systems and providing novel solutions to global problems. The field of genetic engineering has facilitated the manipulation of genetic codes, thus providing opportunities for the advancement of innovative disease therapies and environmental enhancements. The emergence of bio-molecular simulation represents a significant advancement in this particular field, as it offers the ability to gain microscopic insights into molecular-level biological processes over extended periods. Biomolecular simulation plays a crucial role in advancing our comprehension of organismal mechanisms by establishing connections between molecular structures, interactions, and biological functions. The field of computational biology has demonstrated its significance in deciphering intricate biological enigmas through the utilization of mathematical models and algorithms. The process of decoding the human genome has resulted in the advancement of therapies for a wide range of genetic disorders, while the simulation of biological systems contributes to the identification of novel pharmaceutical compounds. The potential of biomolecular simulation and computational biology is vast and limitless. As the exploration of the underlying principles that govern living organisms progresses, the potential impact of this understanding on cancer treatment, environmental restoration, and other domains is anticipated to be transformative. This review examines the notable advancements achieved in the field of computational biology, emphasizing its potential to revolutionize the comprehension and enhancement of biological systems.

Interaction of narcissoside with α-amylase from Bacillus subtilis and Porcine pancreatic by multi-spectral analysis and molecular dynamics simulation

In this work, interaction mechanism of narcissoside with two α-amylase from Bacillus subtilis (BSA) and Porcine pancreatic (PPA) are comparatively studied by multi-spectral analysis, molecular docking and molecular dynamics simulation. The results prove that narcissoside can statically quench fluorescence of BSA/PPA. Two complexes are mainly formed by hydrogen bond and van der Waals force. With the increase of temperature, the two complexes formed by narcissoside and two enzymes become unstable. At the same experimental temperature, the binding force of narcissoside to PPA is higher than that of BSA. The binding of narcissoside to PPA/BSA increases the hydrophobicity of microenvironment. Moreover, the secondary structure of PPA/BSA is mainly changed by decreasing the α-helix. The optimal binding modes of narcissoside with BSA/PPA are predicted by molecular docking, and the stability of the two complexes is evaluated by molecular dynamics simulations.

In Silico Analysis of the Dextransucrase Obtained From Leuconostoc mesenteroides Strain IBUN 91.2.98

The DSR-IBUN dextransucrase produced by Leuconostoc mesenteroides strain IBUN 91.2.98 has a short production time (4.5 hours), an enzymatic activity of 24.8 U/mL, and a specific activity of purified enzyme 2 times higher (331.6 U/mg) than that reported for similar enzymes. The aim of this study was to generate a structural model that, from an in silico approach, allows a better understanding, from the structural point of view, of the activity obtained by the enzyme of interest, which is key to continue with its study and industry application. For this, we translated the nucleotide sequence of the dsr_IBUN gene. With the primary structure of DSR-IBUN, the in silico prediction of physicochemical parameters, the possible subcellular localization, the presence of signal peptide, and the location of domains and functional and structural motifs of the protein were established. Subsequently, its secondary and tertiary structure were predicted and a homology model of the dextransucrase under study was constructed using Swiss-Model, performing careful template selection. The values obtained for the model, Global Model Quality Estimation (0.63), Quality Mean (-1.49), and root-mean-square deviation (0.09), allow us to affirm that the model for the enzyme dextransucrase DSR-IBUN is of adequate quality and can be used as a source of information for this protein.

Expression, Characterisation, Homology Modelling and Molecular Docking of a Novel M17 Family Leucyl-Aminopeptidase from Bacillus cereus CZ

Leucyl-aminopeptidase (LAP), an important metallopeptidase, hydrolyses amino acid residues from the N-terminus of polypeptides and proteins, acting preferentially on the peptide bond formed by N-terminus leucine. A new leucyl-aminopeptidase was found in Bacillus cereus CZ. Its gene (bclap) contained a 1485 bp ORF encoding 494 amino acids with a molecular weight of 54 kDa. The bcLAP protein was successfully expressed in E. coli BL21(DE3). Optimal activity is obtained at pH 9.0 and 58 °C. The bcLAP displays a moderate thermostability and an alkaline pH adaptation range. Enzymatic activity is dramatically enhanced by Ni2+. EDTA significantly inhibits the enzymatic activity, and bestatin and SDS also show strong inhibition. The three-dimensional model of bcLAP monomer and homohexamer is simulated byPHYRE2 server and SWISS-MODEL server. The docking of bestatin, Leu-Trp, Asp-Trp and Ala-Ala-Gly to bcLAP is performed using AutoDock4.2.5, respectively. Molecular docking results show that the residues Lys260, Asp265, Lys272, Asp283, Asp342, Glu344, Arg346, Gly372 and His437 are involved in the hydrogen bonding with the ligands and zinc ions. There may be two nucleophilic catalytic mechanisms in bcLAP, one involving His 437 or Arg346 and the other involving His437 and Arg346. The bcLAP can hydrolyse the peptide bonds in Leu-Trp, Asp-Trp and Ala-Ala-Gly.

Domain swap facilitates structural transitions of spider silk protein C-terminal domains

Domain swap is a mechanism of protein dimerization where the two interacting domains exchange parts of their structure. Web spiders make use of the process in the connection of C-terminal domains (CTDs) of spidroins, the soluble protein building blocks that form tough silk fibers. Besides providing connectivity and solubility, spidroin CTDs are responsible for inducing structural transitions during passage through an acidified assembly zone within spinning ducts. The underlying molecular mechanisms are elusive. Here, we studied the folding of five homologous spidroin CTDs from different spider species or glands. Four of these are domain-swapped dimers formed by five-helix bundles from spidroins of major and minor ampullate glands. The fifth is a dimer that lacks domain swap, formed by four-helix bundles from a spidroin of a flagelliform gland. Spidroins from this gland do not undergo structural transitions whereas the others do. We found a three-state mechanism of folding and dimerization that was conserved across homologues. In chemical denaturation experiments the native CTD dimer unfolded to a dimeric, partially structured intermediate, followed by full unfolding to denatured monomers. The energetics of the individual folding steps varied between homologues. Contrary to the common belief that domain swap stabilizes protein assemblies, the non-swapped homologue was most stable and folded four orders of magnitude faster than a swapped variant. Domain swap of spidroin CTDs induces an entropic penalty to the folding of peripheral helices, thus unfastening them for acid-induced unfolding within a spinning duct, which primes them for refolding into alternative structures during silk formation.

Systematic identification of conditionally folded intrinsically disordered regions by AlphaFold2

The AlphaFold Protein Structure Database contains predicted structures for millions of proteins. For the majority of human proteins that contain intrinsically disordered regions (IDRs), which do not adopt a stable structure, it is generally assumed that these regions have low AlphaFold2 confidence scores that reflect low-confidence structural predictions. Here, we show that AlphaFold2 assigns confident structures to nearly 15% of human IDRs. By comparison to experimental NMR data for a subset of IDRs that are known to conditionally fold (i.e., upon binding or under other specific conditions), we find that AlphaFold2 often predicts the structure of the conditionally folded state. Based on databases of IDRs that are known to conditionally fold, we estimate that AlphaFold2 can identify conditionally folding IDRs at a precision as high as 88% at a 10% false positive rate, which is remarkable considering that conditionally folded IDR structures were minimally represented in its training data. We find that human disease mutations are nearly fivefold enriched in conditionally folded IDRs over IDRs in general and that up to 80% of IDRs in prokaryotes are predicted to conditionally fold, compared to less than 20% of eukaryotic IDRs. These results indicate that a large majority of IDRs in the proteomes of human and other eukaryotes function in the absence of conditional folding, but the regions that do acquire folds are more sensitive to mutations. We emphasize that the AlphaFold2 predictions do not reveal functionally relevant structural plasticity within IDRs and cannot offer realistic ensemble representations of conditionally folded IDRs.

Structure-Based In Silico Screening of Marine Phlorotannins for Potential Walrus Calicivirus Inhibitor

A new calicivirus isolated from a walrus was reported in 2004. Since unknown marine mammalian zoonotic viruses could pose great risks to human health, this study aimed to develop therapeutic countermeasures to quell any potential outbreak of a pandemic caused by this virus. We first generated a 3D model of the walrus calicivirus capsid protein and identified compounds from marine natural products, especially phlorotannins, as potential walrus calicivirus inhibitors. A 3D model of the target protein was generated using homology modeling based on two publicly available template sequences. The sequence of the capsid protein exhibited 31.3% identity and 42.7% similarity with the reference templates. The accuracy and reliability of the predicted residues were validated via Ramachandran plotting. Molecular docking simulations were performed between the capsid protein 3D model and 17 phlorotannins. Among them, five phlorotannins demonstrated markedly stable docking profiles; in particular, 2,7-phloroglucinol-6,6-bieckol showed favorable structural integrity and stability during molecular dynamics simulations. The results indicate that the phlorotannins are promising walrus calicivirus inhibitors. Overall, the study findings showcase the rapid turnaround of in silico-based drug discovery approaches, providing useful insights for developing potential therapies against novel pathogenic viruses, especially when the 3D structures of the viruses remain experimentally unknown.

E2EDA: Protein Domain Assembly Based on End-to-End Deep Learning

With the development of deep learning, almost all single-domain proteins can be predicted at experimental resolution. However, the structure prediction of multi-domain proteins remains a challenge. Achieving end-to-end protein domain assembly and further improving the accuracy of the full-chain modeling by accurately predicting inter-domain orientation while improving the assembly efficiency will provide significant insights into structure-based drug discovery. In this work, we propose an End-to-End Domain Assembly method based on deep learning, named E2EDA. We first develop RMNet, an EfficientNetV2-based deep learning model that fuses multiple features using an attention mechanism to predict inter-domain rigid motion. Then, the predicted rigid motions are transformed into inter-domain spatial transformations to directly assemble the full-chain model. Finally, the scoring strategy RMscore is designed to select the best model from multiple assembled models. The experimental results show that the average TM-score of the model assembled by E2EDA on the benchmark set (282) is 0.827, which is better than those of other domain assembly methods SADA (0.792) and DEMO (0.730). Meanwhile, on our constructed multi-domain data set from AlphaFold DB, the model reassembled by E2EDA is 7.0% higher in TM-score compared to the full-chain model predicted by AlphaFold2, indicating that E2EDA can capture more accurate inter-domain orientations to improve the quality of the model predicted by AlphaFold2. Furthermore, compared to SADA and AlphaFold2, E2EDA reduced the average runtime on the benchmark by 64.7% and 19.2%, respectively, indicating that E2EDA can significantly improve assembly efficiency through an end-to-end approach. The online server is available at http://zhanglab-bioinf.com/E2EDA.

Functional annotation of haloacid dehalogenase superfamily structural genomics proteins

Haloacid dehalogenases (HAD) are members of a large superfamily that includes many Structural Genomics proteins with poorly characterized functionality. This superfamily consists of multiple types of enzymes that can act as sugar phosphatases, haloacid dehalogenases, phosphonoacetaldehyde hydrolases, ATPases, or phosphate monoesterases. Here, we report on predicted functional annotations and experimental testing by direct biochemical assay for Structural Genomics proteins from the HAD superfamily. To characterize the functions of HAD superfamily members, nine representative HAD proteins and 21 structural genomics proteins are analyzed. Using techniques based on computed chemical and electrostatic properties of individual amino acids, the functions of five structural genomics proteins from the HAD superfamily are predicted and validated by biochemical assays. A dehalogenase-like hydrolase, RSc1362 (Uniprot Q8XZN3, PDB 3UMB) is predicted to be a dehalogenase and dehalogenase activity is confirmed experimentally. Four proteins predicted to be sugar phosphatases are characterized as follows: a sugar phosphatase from Thermophilus volcanium (Uniprot Q978Y6) with trehalose-6-phosphate phosphatase and fructose-6-phosphate phosphatase activity; haloacid dehalogenase-like hydrolase from Bacteroides thetaiotaomicron (Uniprot Q8A2F3; PDB 3NIW) with fructose-6-phosphate phosphatase and sucrose-6-phosphate phosphatase activity; putative phosphatase from Eubacterium rectale (Uniprot D0VWU2; PDB 3DAO) as a sucrose-6-phosphate phosphatase; and hypothetical protein from Geobacillus kaustophilus (Uniprot Q5L139; PDB 2PQ0) as a fructose-6-phosphate phosphatase. Most of these sugar phosphatases showed some substrate promiscuity.

AI/ML combined with next-generation sequencing of VHH immune repertoires enables the rapid identification of de novo humanized and sequence-optimized single domain antibodies: a prospective case study

Introduction: In this study, we demonstrate the feasibility of yeast surface display (YSD) and nextgeneration sequencing (NGS) in combination with artificial intelligence and machine learning methods (AI/ML) for the identification of de novo humanized single domain antibodies (sdAbs) with favorable early developability profiles. Methods: The display library was derived from a novel approach, in which VHH-based CDR3 regions obtained from a llama (Lama glama), immunized against NKp46, were grafted onto a humanized VHH backbone library that was diversified in CDR1 and CDR2. Following NGS analysis of sequence pools from two rounds of fluorescence-activated cell sorting we focused on four sequence clusters based on NGS frequency and enrichment analysis as well as in silico developability assessment. For each cluster, long short-term memory (LSTM) based deep generative models were trained and used for the in silico sampling of new sequences. Sequences were subjected to sequence- and structure-based in silico developability assessment to select a set of less than 10 sequences per cluster for production. Results: As demonstrated by binding kinetics and early developability assessment, this procedure represents a general strategy for the rapid and efficient design of potent and automatically humanized sdAb hits from screening selections with favorable early developability profiles.

Antimicrobial Peptides: The Production of Novel Peptide-Based Therapeutics in Plant Systems

The increased prevalence of antibiotic resistance is alarming and has a significant impact on the economies of emerging and underdeveloped nations. The redundancy of antibiotic discovery platforms (ADPs) and injudicious use of conventional antibiotics has severely impacted millions, across the globe. Potent antimicrobials from biological sources have been extensively explored as a ray of hope to counter the growing menace of antibiotic resistance in the population. Antimicrobial peptides (AMPs) are gaining momentum as powerful antimicrobial therapies to combat drug-resistant bacterial strains. The tremendous therapeutic potential of natural and synthesized AMPs as novel and potent antimicrobials is highlighted by their unique mode of action, as exemplified by multiple research initiatives. Recent advances and developments in antimicrobial discovery and research have increased our understanding of the structure, characteristics, and function of AMPs; nevertheless, knowledge gaps still need to be addressed before these therapeutic options can be fully exploited. This thematic article provides a comprehensive insight into the potential of AMPs as potent arsenals to counter drug-resistant pathogens, a historical overview and recent advances, and their efficient production in plants, defining novel upcoming trends in drug discovery and research. The advances in synthetic biology and plant-based expression systems for AMP production have defined new paradigms in the efficient production of potent antimicrobials in plant systems, a prospective approach to countering drug-resistant pathogens.

RNA contact prediction by data efficient deep learning

On the path to full understanding of the structure-function relationship or even design of RNA, structure prediction would offer an intriguing complement to experimental efforts. Any deep learning on RNA structure, however, is hampered by the sparsity of labeled training data. Utilizing the limited data available, we here focus on predicting spatial adjacencies ("contact maps") as a proxy for 3D structure. Our model, BARNACLE, combines the utilization of unlabeled data through self-supervised pre-training and efficient use of the sparse labeled data through an XGBoost classifier. BARNACLE shows a considerable improvement over both the established classical baseline and a deep neural network. In order to demonstrate that our approach can be applied to tasks with similar data constraints, we show that our findings generalize to the related setting of accessible surface area prediction.

Targeting SLC transporters: small molecules as modulators and therapeutic opportunities

Solute carrier (SLCs) transporters mediate the transport of a broad range of solutes across biological membranes. Dysregulation of SLCs has been associated with various pathologies, including metabolic and neurological disorders, as well as cancer and rare diseases. SLCs are therefore emerging as key targets for therapeutic intervention with several recently approved drugs targeting these proteins. Unlocking this large and complex group of proteins is essential to identifying unknown SLC targets and developing next-generation SLC therapeutics. Recent progress in experimental and computational techniques has significantly advanced SLC research, including drug discovery. Here, we review emerging topics in therapeutic discovery of SLCs, focusing on state-of-the-art approaches in structural, chemical, and computational biology, and discuss current challenges in transporter drug discovery.

Bioinformatic and literature assessment of toxicity and allergenicity of a CRISPR-Cas9 engineered gene drive to control Anopheles gambiae the mosquito vector of human malaria

BACKGROUND

Population suppression gene drive is currently being evaluated, including via environmental risk assessment (ERA), for malaria vector control. One such gene drive involves the dsxFCRISPRh transgene encoding (i) hCas9 endonuclease, (ii) T1 guide RNA (gRNA) targeting the doublesex locus, and (iii) DsRed fluorescent marker protein, in genetically-modified mosquitoes (GMMs). Problem formulation, the first stage of ERA, for environmental releases of dsxFCRISPRh previously identified nine potential harms to the environment or health that could occur, should expressed products of the transgene cause allergenicity or toxicity.

METHODS

Amino acid sequences of hCas9 and DsRed were interrogated against those of toxins or allergens from NCBI, UniProt, COMPARE and AllergenOnline bioinformatic databases and the gRNA was compared with microRNAs from the miRBase database for potential impacts on gene expression associated with toxicity or allergenicity. PubMed was also searched for any evidence of toxicity or allergenicity of Cas9 or DsRed, or of the donor organisms from which these products were originally derived.

RESULTS

While Cas9 nuclease activity can be toxic to some cell types in vitro and hCas9 was found to share homology with the prokaryotic toxin VapC, there was no evidence from previous studies of a risk of toxicity to humans and other animals from hCas9. Although hCas9 did contain an 8-mer epitope found in the latex allergen Hev b 9, the full amino acid sequence of hCas9 was not homologous to any known allergens. Combined with a lack of evidence in the literature of Cas9 allergenicity, this indicated negligible risk to humans of allergenicity from hCas9. No matches were found between the gRNA and microRNAs from either Anopheles or humans. Moreover, potential exposure to dsxFCRISPRh transgenic proteins from environmental releases was assessed as negligible.

CONCLUSIONS

Bioinformatic and literature assessments found no convincing evidence to suggest that transgenic products expressed from dsxFCRISPRh were allergens or toxins, indicating that environmental releases of this population suppression gene drive for malaria vector control should not result in any increased allergenicity or toxicity in humans or animals. These results should also inform evaluations of other GMMs being developed for vector control and in vivo clinical applications of CRISPR-Cas9.

High-throughput screening and investigation of the inhibitory mechanism of α-glucosidase inhibitors in teas using an affinity selection-mass spectrometry method

An affinity selection-mass spectrometry method was applied for high-throughput screening of α-glucosidase (AGH) inhibitors from teas. Fourteen out of nineteen screened AGH inhibitor candidates were clustered as galloylated polyphenols (GPs). "AGH-GPs" interaction studies, including enzyme kinetics, fluorescence spectroscopy, circular dichroism, and molecular docking, jointly suggested that GPs noncompetitively inhibit AGH activity by interacting with amino acid residues near the active site of AGH and inducing changes in AGH secondary structure. Representative GPs and white tea extract (WTE) showed comparable AGH inhibition effects in Caco2 cells and postprandial hypoglycemic efficacy in diabetic mice as acarbose. The area under the curve of oral sucrose tolerance test was lower by 8.16%, 6.17%, and 7.37% than control group in 15 mg/kg EGCG, 15 mg/kg strictinin, and 150 mg/kg WTE group, respectively. Our study presents a high-efficiency approach to discover novel AGH inhibitors and elucidates a potential mechanism by which tea decreases diabetes risks.

An agnostic analysis of the human AlphaFold2 proteome using local protein conformations

Knowledge of the 3D structure of proteins is a valuable asset for understanding their precise biological mechanisms. However, the cost of production of 3D structures and experimental difficulties limit their obtaining. The proposal of 3D structural models is consequently an appealing alternative. The release of the AlphaFold Deep Learning approach has revolutionized the field. The recent near-complete human proteome proposal makes it possible to analyse large amounts of data and evaluate the results of the approach in greater depth. The 3D human proteome was thus analysed in light of the classic secondary structures, and many less-used protein local conformations (PolyProline II helices, type of γ-turns, of β-turns and of β-bulges, curvature of the helices, and a structural alphabet). Without questioning the global quality of the approach, this analysis highlights certain local conformations, which maybe poorly predicted and they could therefore be better addressed.

Protein Design Using Physics Informed Neural Networks

The inverse protein folding problem, also known as protein sequence design, seeks to predict an amino acid sequence that folds into a specific structure and performs a specific function. Recent advancements in machine learning techniques have been successful in generating functional sequences, outperforming previous energy function-based methods. However, these machine learning methods are limited in their interoperability and robustness, especially when designing proteins that must function under non-ambient conditions, such as high temperature, extreme pH, or in various ionic solvents. To address this issue, we propose a new Physics-Informed Neural Networks (PINNs)-based protein sequence design approach. Our approach combines all-atom molecular dynamics simulations, a PINNs MD surrogate model, and a relaxation of binary programming to solve the protein design task while optimizing both energy and the structural stability of proteins. We demonstrate the effectiveness of our design framework in designing proteins that can function under non-ambient conditions.

ModelCIF: An Extension of PDBx/mmCIF Data Representation for Computed Structure Models

ModelCIF (github.com/ihmwg/ModelCIF) is a data information framework developed for and by computational structural biologists to enable delivery of Findable, Accessible, Interoperable, and Reusable (FAIR) data to users worldwide. ModelCIF describes the specific set of attributes and metadata associated with macromolecular structures modeled by solely computational methods and provides an extensible data representation for deposition, archiving, and public dissemination of predicted three-dimensional (3D) models of macromolecules. It is an extension of the Protein Data Bank Exchange / macromolecular Crystallographic Information Framework (PDBx/mmCIF), which is the global data standard for representing experimentally-determined 3D structures of macromolecules and associated metadata. The PDBx/mmCIF framework and its extensions (e.g., ModelCIF) are managed by the Worldwide Protein Data Bank partnership (wwPDB, wwpdb.org) in collaboration with relevant community stakeholders such as the wwPDB ModelCIF Working Group (wwpdb.org/task/modelcif). This semantically rich and extensible data framework for representing computed structure models (CSMs) accelerates the pace of scientific discovery. Herein, we describe the architecture, contents, and governance of ModelCIF, and tools and processes for maintaining and extending the data standard. Community tools and software libraries that support ModelCIF are also described.

Application of computational methods for class A GPCR Ligand discovery

G protein-coupled receptors (GPCR) are integral membrane proteins of considerable interest as targets for drug development due to their role in transmitting cellular signals in a multitude of biological processes. Of the six classes categorizing GPCR (A, B, C, D, E, and F), class A contains the largest number of therapeutically relevant GPCR. Despite their importance as drug targets, many challenges exist for the discovery of novel class A GPCR ligands serving as drug precursors. Though knowledge of the structural and functional characteristics of GPCR has grown significantly over the past 20 years, a large portion of GPCR lack reported, experimentally determined structures. Furthermore, many GPCR have no known endogenous and/or synthetic ligands, limiting further exploration of their biochemical, cellular, and physiological roles. While many successes in GPCR ligand discovery have resulted from experimental high-throughput screening, computational methods have played an increasingly important role in GPCR ligand identification in the past decade. Here we discuss computational techniques applied to GPCR ligand discovery. This review summarizes class A GPCR structure/function and provides an overview of many obstacles currently faced in GPCR ligand discovery. Furthermore, we discuss applications and recent successes of computational techniques used to predict GPCR structure as well as present a summary of ligand- and structure-based methods used to identify potential GPCR ligands. Finally, we discuss computational hit list generation and refinement and provide comprehensive workflows for GPCR ligand identification.

Identification and Characterization of Three New Antimicrobial Peptides from the Marine Mollusk Nerita versicolor (Gmelin, 1791)

Mollusks have been widely investigated for antimicrobial peptides because their humoral defense against pathogens is mainly based on these small biomolecules. In this report, we describe the identification of three novel antimicrobial peptides from the marine mollusk Nerita versicolor. A pool of N. versicolor peptides was analyzed with nanoLC-ESI-MS-MS technology, and three potential antimicrobial peptides (Nv-p1, Nv-p2 and Nv-p3) were identified with bioinformatical predictions and selected for chemical synthesis and evaluation of their biological activity. Database searches showed that two of them show partial identity to histone H4 peptide fragments from other invertebrate species. Structural predictions revealed that they all adopt a random coil structure even when placed near a lipid bilayer patch. Nv-p1, Nv-p2 and Nv-p3 exhibited activity against Pseudomonas aeruginosa. The most active peptide was Nv-p3 with an inhibitory activity starting at 1.5 µg/mL in the radial diffusion assays. The peptides were ineffective against Klebsiella pneumoniae, Listeria monocytogenes and Mycobacterium tuberculosis. On the other hand, these peptides demonstrated effective antibiofilm action against Candida albicans, Candida parapsilosis and Candida auris but not against the planktonic cells. None of the peptides had significant toxicity on primary human macrophages and fetal lung fibroblasts at effective antimicrobial concentrations. Our results indicate that N. versicolor-derived peptides represent new AMP sequences and have the potential to be optimized and developed into antibiotic alternatives against bacterial and fungal infections.

Structural and functional characterisation of HepTH1-5 peptide as a potential hepcidin replacement

Hepcidin is a principal regulator of iron homeostasis and its dysregulation has been recognised as a causative factor in cancers and iron disorders. The strategy of manipulating the presence of hepcidin peptide has been used for cancer treatment. However, this has demonstrated poor efficiency and has been short-lived in patients. Many studies reported using minihepcidin therapy as an alternative way to treat hepcidin dysregulation, but this was only applied to non-cancer patients. Highly conserved fish hepcidin protein, HepTH1-5, was investigated to determine its potential use in developing a hepcidin replacement for human hepcidin (Hepc25) and as a therapeutic agent by targeting the tumour suppressor protein, p53, through structure-function analysis. The authors found that HepTH1-5 is stably bound to ferroportin, compared to Hepc25, by triggering the ferroportin internalisation via Lys42 and Lys270 ubiquitination, in a similar manner to the Hepc25 activity. Moreover, the residues Ile24 and Gly24, along with copper and zinc ligands, interacted with similar residues, Lys24 and Asp1 of Hepc25, respectively, showing that those molecules are crucial to the hepcidin replacement strategy. HepTH1-5 interacts with p53 and activates its function through phosphorylation. This finding shows that HepTH1-5 might be involved in the apoptosis signalling pathway upon a DNA damage response. This study will be very helpful for understanding the mechanism of the hepcidin replacement and providing insights into the HepTH1-5 peptide as a new target for hepcidin and cancer therapeutics.Communicated by Ramaswamy H. Sarma.

Electromagnetic Forces and Torques: From Dielectrophoresis to Optical Tweezers

Electromagnetic forces and torques enable many key technologies, including optical tweezers or dielectrophoresis. Interestingly, both techniques rely on the same physical process: the interaction of an oscillating electric field with a particle of matter. This work provides a unified framework to understand this interaction both when considering fields oscillating at low frequencies─dielectrophoresis─and high frequencies─optical tweezers. We draw useful parallels between these two techniques, discuss the different and often unstated assumptions they are based upon, and illustrate key applications in the fields of physical and analytical chemistry, biosensing, and colloidal science.

Contribution of κ-/ι-carrageenan on the gelling properties of shrimp myofibrillar protein and their interaction mechanism exploration

BACKGROUND

The contribution and mechanism of κ-/ι-carrageenan (CG) with different hydration characteristics on the gelling properties of shrimp myofibrillar protein (MP) gelation was studied.

RESULTS

The gel strength, water-holding capacity and viscoelastic properties of MP gels were significantly enhanced by 1.0% κ-/ι-CG (P < 0.05), but the microstructure showed that excessive carrageenan caused fragmentation of the gel network and a corresponding decrease in gel properties. Compared to MP-ιCG, MP-κCG showed larger breaking force and shorter breaking distance, thus enhancing the hardness and brittleness of the gel, which might be ascribed to a reinforced network skeleton and a tighter binding of κCG-myosin. However, MP-ιCG stabilized more moisture in the gel network, thereby improving the tenderness of the gel, which might be related to the electrostatic repulsion observed between the sulfate groups of ιCG and the myosin observed by molecular docking. In addition, the β-sheet content and intermolecular interactions might be positively correlated with gel properties.

CONCLUSION

In this study, a composite gel system was constructed based on the interaction of MP and CG. The quality differences of two kinds of CG-MP gels were clarified, which will provide guidance for the application of different kinds of carrageenan and the development of recombinant meat products with specific quality. © 2022 Society of Chemical Industry.

The In Silico Characterization of Monocotyledonous α-l-Arabinofuranosidases on the Example of Maize

Plant α-l-arabinofuranosidases remove terminal arabinose from arabinose-containing substrates such as plant cell wall polysaccharides, including arabinoxylans, arabinogalactans, and arabinans. In plants, de-arabinosylation of cell wall polysaccharides accompanies different physiological processes such as fruit ripening and elongation growth. In this report, we address the diversity of plant α-l-arabinofuranosidases of the glycoside hydrolase (GH) family 51 through their phylogenetic analysis as well as their structural features. The CBM4-like domain at N-terminus was found to exist only in GH51 family proteins and was detected in almost 90% of plant sequences. This domain is similar to bacterial CBM4, but due to substitutions of key amino acid residues, it does not appear to be able to bind carbohydrates. Despite isoenzymes of GH51 being abundant, in particular in cereals, almost half of the GH51 proteins in Poales have a mutation of the acid/base residue in the catalytic site, making them potentially inactive. Open-source data on the transcription and translation of GH51 isoforms in maize were analyzed to discuss possible functions of individual isoenzymes. The results of homology modeling and molecular docking showed that the substrate binding site can accurately accommodate terminal arabinofuranose and that arabinoxylan is a more favorable ligand for all maize GH51 enzymes than arabinan.

RCSB Protein Data Bank (RCSB.org): delivery of experimentally-determined PDB structures alongside one million computed structure models of proteins from artificial intelligence/machine learning

The Research Collaboratory for Structural Bioinformatics Protein Data Bank (RCSB PDB), founding member of the Worldwide Protein Data Bank (wwPDB), is the US data center for the open-access PDB archive. As wwPDB-designated Archive Keeper, RCSB PDB is also responsible for PDB data security. Annually, RCSB PDB serves >10 000 depositors of three-dimensional (3D) biostructures working on all permanently inhabited continents. RCSB PDB delivers data from its research-focused RCSB.org web portal to many millions of PDB data consumers based in virtually every United Nations-recognized country, territory, etc. This Database Issue contribution describes upgrades to the research-focused RCSB.org web portal that created a one-stop-shop for open access to ∼200 000 experimentally-determined PDB structures of biological macromolecules alongside >1 000 000 incorporated Computed Structure Models (CSMs) predicted using artificial intelligence/machine learning methods. RCSB.org is a 'living data resource.' Every PDB structure and CSM is integrated weekly with related functional annotations from external biodata resources, providing up-to-date information for the entire corpus of 3D biostructure data freely available from RCSB.org with no usage limitations. Within RCSB.org, PDB structures and the CSMs are clearly identified as to their provenance and reliability. Both are fully searchable, and can be analyzed and visualized using the full complement of RCSB.org web portal capabilities.

Contact-Assisted Threading in Low-Homology Protein Modeling

The ability to successfully predict the three-dimensional structure of a protein from its amino acid sequence has made considerable progress in the recent past. The progress is propelled by the improved accuracy of deep learning-based inter-residue contact map predictors coupled with the rising growth of protein sequence databases. Contact map encodes interatomic interaction information that can be exploited for highly accurate prediction of protein structures via contact map threading even for the query proteins that are not amenable to direct homology modeling. As such, contact-assisted threading has garnered considerable research effort. In this chapter, we provide an overview of existing contact-assisted threading methods while highlighting the recent advances and discussing some of the current limitations and future prospects in the application of contact-assisted threading for improving the accuracy of low-homology protein modeling.

Homology Modeling of Transporter Proteins

Membrane transporter proteins are divided into channels/pores and carriers and constitute protein families of physiological and pharmacological importance. Several presently used therapeutic compounds elucidate their effects by targeting membrane transporter proteins, including anti-arrhythmic, anesthetic, antidepressant, anxiolytic and diuretic drugs. The lack of three-dimensional structures of human transporters hampers experimental studies and drug discovery. In this chapter, the use of homology modeling for generating structural models of membrane transporter proteins is reviewed. The increasing number of atomic resolution structures available as templates, together with improvements in methods and algorithms for sequence alignments, secondary structure predictions, and model generation, in addition to the increase in computational power have increased the applicability of homology modeling for generating structural models of transporter proteins. Different pitfalls and hints for template selection, multiple-sequence alignments, generation and optimization, validation of the models, and the use of transporter homology models for structure-based virtual ligand screening are discussed.

Protein-Protein Interaction Studies Using Molecular Dynamics Simulation

Protein-protein interaction (PPI) is a crucial event for many biological functions. Studying the molecular details of PPI requires structure determination using X-ray crystallography, nuclear magnetic resistance (NMR), and single particle Cryo-EM. However, sometimes it is not easy to solve the complex structure for various reasons. For example, complex may be unstable, not enough protein expression for structural studies, etc. Further, PPI are intricate processes, and its molecular details cannot be fully explained by experimental observations. Here, we describe a quick and simple method to study the PPI using the combinatorial approach of molecular dynamics simulation and biophysical methods.

Exploring microbial functional biodiversity at the protein family level-From metagenomic sequence reads to annotated protein clusters

Metagenomics has enabled accessing the genetic repertoire of natural microbial communities. Metagenome shotgun sequencing has become the method of choice for studying and classifying microorganisms from various environments. To this end, several methods have been developed to process and analyze the sequence data from raw reads to end-products such as predicted protein sequences or families. In this article, we provide a thorough review to simplify such processes and discuss the alternative methodologies that can be followed in order to explore biodiversity at the protein family level. We provide details for analysis tools and we comment on their scalability as well as their advantages and disadvantages. Finally, we report the available data repositories and recommend various approaches for protein family annotation related to phylogenetic distribution, structure prediction and metadata enrichment.

Structural basis for the binding of a selective inverse agonist AF64394 with the human G-protein coupled receptor 3 (GPR3)

The orphan class A G-protein coupled receptor 3 (GPR3) is highly expressed in brain and linked with various neuronal functions, and therefore, expected to play a vital role in the progression of Alzheimer's disease. In view of the lack of its experimental structure, we describe herein the three-dimensional structure and conformational dynamics of GPR3 complexed with the inverse agonist AF64394. The GPR3 model was predicted using the Iterative Threading ASSEmbly Refinement (I-TASSER) method. The Induced Fit Docking predicted two unique poses, Pose 1 and Pose 2, for AF64394, and then, molecular dynamics (MD) simulations followed by binding free-energy calculation revealed the Pose 1 as a very stable pose with the least fluctuation during the MD simulation while the Pose 2 underwent a significant fluctuation. The [1,2,4]triazolo[1,5-a]pyrimidine core was engaged in multiple hydrogen bonds (H-bonds), such as a water-mediated H-bond between the triazole nitrogen and T31, two direct H-bonds between the protonated triazole-ring nitrogen and V186 and T279, a direct H-bond between the secondary amine and V187. The phenyl substituent of AF64394 exhibited aromatic π-π stacking interactions with F97, F101, W43 and Y280. AF64394 showed a direct interaction with E28 and polar interactions with H96, T31 and T279. Throughout the MD simulation, the toggle switch residues, F120 and W260, remained in close contact, indicating that the GPR3 conformation represented an inactive state. The 4-(3-chloro-5-isopropoxyphenethyl) group resided near to the toggle switch residues. The insights gained here are expected to be useful in the structure-based design of new ligands targeting GPR3 modulation. Communicated by Ramaswamy H. Sarma.

Structure-Function Study of a Novel Inhibitor of Cyclin-Dependent Kinase C in Arabidopsis

The circadian clock, an internal time-keeping system with a period of about 24 h, coordinates many physiological processes with the day-night cycle. We previously demonstrated that BML-259 [N-(5-isopropyl-2-thiazolyl) phenylacetamide], a small molecule with mammal CYCLIN DEPENDENT KINASE 5 (CDK5)/CDK2 inhibition activity, lengthens Arabidopsis thaliana (Arabidopsis) circadian clock periods. BML-259 inhibits Arabidopsis CDKC kinase, which phosphorylates RNA polymerase II in the general transcriptional machinery. To accelerate our understanding of the inhibitory mechanism of BML-259 on CDKC, we performed structure-function studies of BML-259 using circadian period-lengthening activity as an estimation of CDKC inhibitor activity in vivo. The presence of a thiazole ring is essential for period-lengthening activity, whereas acetamide, isopropyl and phenyl groups can be modified without effect. BML-259 analog TT-539, a known mammal CDK5 inhibitor, did not lengthen the period nor did it inhibit Pol II phosphorylation. TT-361, an analog having a thiophenyl ring instead of a phenyl ring, possesses stronger period-lengthening activity and CDKC;2 inhibitory activity than BML-259. In silico ensemble docking calculations using Arabidopsis CDKC;2 obtained by a homology modeling indicated that the different binding conformations between these molecules and CDKC;2 explain the divergent activities of TT539 and TT361.

Research progress on unique paratope structure, antigen binding modes, and systematic mutagenesis strategies of single-domain antibodies

Single-domain antibodies (sdAbs) showed the incredible advantages of small molecular weight, excellent affinity, specificity, and stability compared with traditional IgG antibodies, so their potential in binding hidden antigen epitopes and hazard detection in food, agricultural and veterinary fields were gradually explored. Moreover, its low immunogenicity, easy-to-carry target drugs, and penetration of the blood-brain barrier have made sdAbs remarkable achievements in medical treatment, toxin neutralization, and medical imaging. With the continuous development and maturity of modern molecular biology, protein analysis software and database with different algorithms, and next-generation sequencing technology, the unique paratope structure and different antigen binding modes of sdAbs compared with traditional IgG antibodies have aroused the broad interests of researchers with the increased related studies. However, the corresponding related summaries are lacking and needed. Different antigens, especially hapten antigens, show distinct binding modes with sdAbs. So, in this paper, the unique paratope structure of sdAbs, different antigen binding cases, and the current maturation strategy of sdAbs were classified and summarized. We hope this review lays a theoretical foundation to elucidate the antigen-binding mechanism of sdAbs and broaden the further application of sdAbs.

Biology and medicine in the landscape of quantum advantages

Quantum computing holds substantial potential for applications in biology and medicine, spanning from the simulation of biomolecules to machine learning methods for subtyping cancers on the basis of clinical features. This potential is encapsulated by the concept of a quantum advantage, which is contingent on a reduction in the consumption of a computational resource, such as time, space or data. Here, we distill the concept of a quantum advantage into a simple framework to aid researchers in biology and medicine pursuing the development of quantum applications. We then apply this framework to a wide variety of computational problems relevant to these domains in an effort to (i) assess the potential of practical advantages in specific application areas and (ii) identify gaps that may be addressed with novel quantum approaches. In doing so, we provide an extensive survey of the intersection of biology and medicine with the current landscape of quantum algorithms and their potential advantages. While we endeavour to identify specific computational problems that may admit practical advantages throughout this work, the rapid pace of change in the fields of quantum computing, classical algorithms and biological research implies that this intersection will remain highly dynamic for the foreseeable future.

Multi-state modeling of G-protein coupled receptors at experimental accuracy

The family of G-protein coupled receptors (GPCRs) is one of the largest protein families in the human genome. GPCRs transduct chemical signals from extracellular to intracellular regions via a conformational switch between active and inactive states upon ligand binding. While experimental structures of GPCRs remain limited, high-accuracy computational predictions are now possible with AlphaFold2. However, AlphaFold2 only predicts one state and is biased toward either the active or inactive conformation depending on the GPCR class. Here, a multi-state prediction protocol is introduced that extends AlphaFold2 to predict either active or inactive states at very high accuracy using state-annotated templated GPCR databases. The predicted models accurately capture the main structural changes upon activation of the GPCR at the atomic level. For most of the benchmarked GPCRs (10 out of 15), models in the active and inactive states were closer to their corresponding activation state structures. Median RMSDs of the transmembrane regions were 1.12 Å and 1.41 Å for the active and inactive state models, respectively. The models were more suitable for protein-ligand docking than the original AlphaFold2 models and template-based models. Finally, our prediction protocol predicted accurate GPCR structures and GPCR-peptide complex structures in GPCR Dock 2021, a blind GPCR-ligand complex modeling competition. We expect that high accuracy GPCR models in both activation states will promote understanding in GPCR activation mechanisms and drug discovery for GPCRs. At the time, the new protocol paves the way towards capturing the dynamics of proteins at high-accuracy via machine-learning methods.

Structure-Based Virtual Screening and De Novo Design to Identify Submicromolar Inhibitors of G2019S Mutant of Leucine-Rich Repeat Kinase 2

Missense mutations of leucine-rich repeat kinase 2 (LRRK2), including the G2019S mutant, are responsible for the pathogenesis of Parkinson's disease. In this work, structure-based virtual screening of a large chemical library was carried out to identify a number of novel inhibitors of the G2019S mutant of LRRK2, the biochemical potencies of which ranged from the low micromolar to the submicromolar level. The discovery of these potent inhibitors was made possible due to the modification of the original protein-ligand binding energy function in order to include an accurate ligand dehydration energy term. The results of extensive molecular docking simulations indicated that the newly identified inhibitors were bound to the ATP-binding site of the G2019S mutant of LRRK2 through the multiple hydrogen bonds with backbone amide groups in the hinge region as well as the hydrophobic interactions with the nonpolar residues in the P-loop, hinge region, and interdomain region. Among 18 inhibitors derived from virtual screening, 4-(2-amino-5-phenylpyrimidin-4-yl)benzene-1,3-diol (Inhibitor 2) is most likely to serve as a new molecular scaffold to optimize the biochemical potency, because it revealed submicromolar inhibitory activity in spite of its low molecular weight (279.3 amu). Indeed, a highly potent inhibitor (Inhibitor 2n) of the G2019S mutant was derived via the structure-based de novo design using the structure of Inhibitor 2 as the molecular core. The biochemical potency of Inhibitor 2n surged to the nanomolar level due to the strengthening of hydrophobic interactions in the ATP-binding site, which were presumably caused by the substitutions of small nonpolar moieties. Due to the high biochemical potency against the G2019S mutant of LRRK2 and the putatively good physicochemical properties, Inhibitor 2n is anticipated to serve as a new lead compound for the discovery of antiparkinsonian medicines.

BESFA: bioinformatics based evolutionary, structural & functional analysis of prostrate, Placenta, Ovary, Testis, and Embryo (POTE) paralogs

The POTE family comprises 14 paralogues and is primarily expressed in Prostrate, Placenta, Ovary, Testis, Embryo (POTE), and cancerous cells. The prospective function of the POTE protein family under physiological conditions is less understood. We systematically analyzed their cellular localization and molecular docking analysis to elucidate POTE proteins' structure, function, and Adaptive Divergence. Our results suggest that group three POTE paralogs (POTEE, POTEF, POTEI, POTEJ, and POTEKP (a pseudogene)) exhibits significant variation among other members could be because of their Adaptive Divergence. Furthermore, our molecular docking studies on POTE protein revealed the highest binding affinity with NCI-approved anticancer compounds. Additionally, POTEE, POTEF, POTEI, and POTEJ were subject to an explicit molecular dynamic simulation for 50ns. MM-GBSA and other essential electrostatics were calculated that showcased that only POTEE and POTEF have absolute binding affinities with minimum energy exploitation. Thus, this study’s outcomes are expected to drive cancer research to successful utilization of POTE genes family as a new biomarker, which could pave the way for the discovery of new therapies.

Identification of long-chain alkane-degrading (LadA) monooxygenases in Aspergillus flavus via in silico analysis

Efficient degradation of alkanes in crude oil by the isolated Aspergillus flavus MM1 alluded to the presence of highly active alkane-degrading enzymes in this fungus. A long-chain alkane-degrading, LadA-like enzyme family in A. flavus was identified, and possible substrate-binding modes were analyzed using a computational approach. By analyzing publicly available protein databases, we identified six uncharacterized proteins in A. flavus NRRL 3357, of which five were identified as class LadAα and one as class LadAβ, which are eukaryotic homologs of bacterial long-chain alkane monooxygenase (LadA). Computational models of A. flavus LadAα homologs (Af1-Af5) showed overall structural similarity to the bacterial LadA and the unique sequence and structural elements that bind the cofactor Flavin mononucleotide (FMN). A receptor-cofactor-substrate docking protocol was established and validated to demonstrate the substrate binding in the A. flavus LadAα homologs. The modeled Af1, Af3, Af4, and Af5 captured long-chain n-alkanes inside the active pocket, above the bound FMN. Isoalloxazine ring of reduced FMN formed a π–alkyl interaction with the terminal carbon atom of captured alkanes, C₁₆–C₃₀, in Af3–Af5 and C₁₆–C₂₄ in Af1. Our results confirmed the ability of identified A. flavus LadAα monooxygenases to bind long-chain alkanes inside the active pocket. Hence A. flavus LadAα monooxygenases potentially initiate the degradation of long-chain alkanes by oxidizing bound long-chain alkanes into their corresponding alcohol.

Characterization of a novel GH10 alkali-thermostable xylanase from a termite microbiome

The aim of the present study was to assess the biochemical and molecular structural characteristics of a novel alkali-thermostable GH10 xylanase (Xyl10B) identified in a termite gut microbiome by a shotgun metagenomic approach. This endoxylanase candidate was amplified, cloned, heterologously expressed in Escherichia coli and purified. The recombinant enzyme was active at a broad range of temperatures (37-60 ºC) and pH values (4-10), with optimal activity at 50 ºC and pH 9. Moreover, its activity remained at more than 80% of its maximum at 50 °C for 8 h. In addition, Xyl10B was found to be stable in the presence of salt and several ions and chemical reagents frequently used in the industry. These characteristics make this enzyme an interesting candidate for pulp and paper bleaching industries, since this process requires enzymes without cellulase activity and resistant to high temperatures and alkaline pH (thermo-alkaliphilic enzymes). The products of xylan hydrolysis by Xyl10B (short xylooligosaccharides, xylose and xylobiose) could be suitable for application as prebiotics and in the production of bioethanol.