Low Sequence Identity Research Articles

Identification and validation of fusidic acid and flufenamic acid as inhibitors of SARS-CoV-2 replication using DrugSolver CavitomiX

In this work, we present DrugSolver CavitomiX, a novel computational pipeline for drug repurposing and identifying ligands and inhibitors of target enzymes. The pipeline is based on cavity point clouds representing physico-chemical properties of the cavity induced solely by the protein. To test the pipeline’s ability to identify inhibitors, we chose enzymes essential for SARS-CoV-2 replication as a test system. The active-site cavities of the viral enzymes main protease (Mpro) and papain-like protease (Plpro), as well as of the human transmembrane serine protease 2 (TMPRSS2), were selected as target cavities. Using active-site point-cloud comparisons, it was possible to identify two compounds—flufenamic acid and fusidic acid—which show strong inhibition of viral replication. The complexes from which fusidic acid and flufenamic acid were derived would not have been identified using classical sequence- and structure-based methods as they show very little structural (TM-score: 0.1 and 0.09, respectively) and very low sequence (~ 5%) identity to Mpro and TMPRSS2, respectively. Furthermore, a cavity-based off-target screening was performed using acetylcholinesterase (AChE) as an example. Using cavity comparisons, the human carboxylesterase was successfully identified, which is a described off-target for AChE inhibitors.

Leveraging protein language models for accurate multiple sequence alignments.

Multiple sequence alignment (MSA) is a critical step in the study of protein sequence and function. Typically, MSA algorithms progressively align pairs of sequences and combine these alignments with the aid of a guide tree. These alignment algorithms use scoring systems based on substitution matrices to measure amino acid similarities. Although successful, standard methods struggle on sets of proteins with low sequence identity: the so-called twilight zone of protein alignment. For these difficult cases, another source of information is needed. Protein language models are a powerful new approach that leverages massive sequence data sets to produce high-dimensional contextual embeddings for each amino acid in a sequence. These embeddings have been shown to reflect physicochemical and higher-order structural and functional attributes of amino acids within proteins. Here, we present a novel approach to MSA, based on clustering and ordering amino acid contextual embeddings. Our method for aligning semantically consistent groups of proteins circumvents the need for many standard components of MSA algorithms, avoiding initial guide tree construction, intermediate pairwise alignments, gap penalties, and substitution matrices. The added information from contextual embeddings leads to higher accuracy alignments for structurally similar proteins with low amino-acid similarity. We anticipate that protein language models will become a fundamental component of the next generation of algorithms for generating MSAs.

Isolation of a novel feather-degrading Ectobacillus sp. JY-23 strain and characterization of a new keratinase in the M4 metalloprotease family

Microbial keratinases have prominent potential in biotransformation of recalcitrant keratin substrates to value-added products which has made keratinases a research focus in the past decades. In this study, an efficient feather-degrading bacterium was isolated and identified as a novel species in Ectobacillus genus and designated as Ectobacillus sp. JY-23. The degradation characteristics analysis revealed that Ectobacillus sp. JY-23 could utilize chicken feathers (0.4% w/v) as the sole nutrient source and degraded 92.95% of feathers in 72 h. A significant increase in sulfite and free sulfydryl group content detected in the feather hydrolysate (culture supernatant) indicated efficient reduction of disulfide bonds, which inferred that the degradation mechanism of isolated strain was a synergetic action of sulfitolysis and proteolysis. Moreover, abundant amino acids were also detected, among which proline and glycine were the predominant free amino acids. Then, the keratinase of Ectobacillus sp. JY-23 was mined and Y1_15990 was identified as the keratinase encoding gene of Ectobacillus sp. JY-23 and designated as kerJY-23. Escherichia coli strain overexpressing kerJY-23 degraded chicken feathers in 48 h. Finally, bioinformatics prediction of KerJY-23 demonstrated that it belonged to the M4 metalloprotease family, which was a third keratinase member in this family. KerJY-23 showed low sequence identity to the other two keratinase members, indicating the novelty of KerJY-23. Overall, this study presents a novel feather-degrading bacterium and a new keratinase in the M4 metalloprotease family with remarkable potential in feather keratin valorization.

Biosynthesis of meso-lanthionine in Fusobacterium nucleatum

The opportunistic oral pathogen, Fusobacterium nucleatum contains meso-lanthionine as the diaminodicarboxylic acid in the pentapeptide crosslink of the peptidoglycan layer. The diastereomer, l,l-lanthionine is formed by lanthionine synthase, a PLP-dependent enzyme that catalyzes the β-replacement of l-cysteine with a second equivalent of l-cysteine. In this study, we explored possible enzymatic mechanisms for the formation of meso-lanthionine. Our inhibition studies with lanthionine synthase, described herein, revealed that meso-diaminopimelate, a bioisostere of meso-lanthionine, is a more potent inhibitor of lanthionine synthase compared to the diastereomer, l,l-diaminopimelate. These results suggested that lanthionine synthase could also form meso-lanthionine by the β-replacement of l-cysteine with d-cysteine. Through steady-state and pre-steady state kinetic analysis, we confirm that d-cysteine reacts with the ⍺-aminoacylate intermediate with a kon that was 2-3-fold faster and Kd value that was 2–3fold lower compared to l-cysteine. However, given that intracellular levels of d-cysteine levels are assumed to be significantly lower than that of l-cysteine, we also determined if the gene product, FN1732, with low sequence identity to diaminopimelate epimerase could convert l,l-lanthionine to meso-lanthionine. Using diaminopimelate dehydrogenase in a coupled spectrophotometric assay, we show that FN1732 can convert l,l-lanthionine to meso-lanthionine with a kcat of 0.07 ± 0.001 s−1 and a KM of 1.9 ± 0.1 mM. In summary, our results provide two possible enzymatic mechanisms for the biosynthesis of meso-lanthionine in F. nucleatum.

Design and Synthesis of Novel Oxathiapiprolin Derivatives as Oxysterol Binding Protein Inhibitors and Their Application in Phytopathogenic Oomycetes.

Oomycetes, particularly those from the genus Phytophthora, are significant threats to global food security and natural ecosystems. Oxathiapiprolin (OXA) is an effective oomycete fungicide that targets an oxysterol binding protein (OSBP), while the binding mechanism of OXA is still unclear, which limits the pesticide design, induced by the low sequence identity of Phytophthora and template models. Herein, we generated the OSBP model of the well-reported Phytophthora capsici using AlphaFold 2 and studied the binding mechanism of OXA. Based on it, a series of OXA analogues were designed. Then, compound 2l, the most potent candidate, was successfully designed and synthesized, showing a control efficiency comparable to that of OXA. Moreover, field trial experiments showed that 2l exhibited nearly the same activity (72.4%) as OXA against cucumber downy mildew at 25 g/ha. The present work indicated that 2l could be used as a leading compound for the discovery of new OSBP fungicides.

The bacterial effector SidN/Lpg1083 promotes cell death by targeting Lamin-B2.

To facilitate survival, replication, and dissemination, the intracellular pathogen Legionella pneumophila relies on its unique type IVB secretion system (T4SS) to deliver over 330 effectors to hijack host cell pathways in a spatiotemporal manner. The effectors and their host targets are largely unexplored due to their low sequence identity to the known proteins and functional redundancy. The T4SS effector SidN (Lpg1083) is secreted into host cells during the late infection period. However, to the best of our knowledge, the molecular characterization of SidN has not been studied. Herein, we identified SidN as a nuclear envelope-localized effector. Its structure adopts a novel fold, and the N-terminal domain is crucial for its specific subcellular localization. Furthermore, we found that SidN is transported by eukaryotic karyopherin Importin-13 into the nucleus, where it attaches to the N-terminal region of Lamin-B2 to interfere with the integrity of the nuclear envelope, causing nuclear membrane disruption and eventually cell death. Our work provides new insights into the structure and function of an L. pneumophila effector protein, and suggests a potential strategy utilized by the pathogen to promote host cell death and then escape from the host for secondary infection.

NtLTPI.38, a plasma membrane-localized protein, mediates lipid metabolism and salt tolerance in Nicotiana tabacum

Non-specific lipid transfer proteins (nsLTPs) typically have conserved structural resemblance, low sequence identity, and broad biological functions in plant growth and stress resistance. Here, a plasma membrane-localized nsLTP, NtLTPI.38, was identified in tobacco plants. Multi-omics integrated analysis revealed that NtLTPI.38 overexpression or knock out significantly changed glycerophospholipid and glycerolipid metabolism pathways. NtLTPI.38 overexpression remarkably increased phosphatidylcholine, phosphatidylethanolamine, triacylglycerol, and flavonoid levels, but decreased ceramides compared to wild type and mutant lines. Differentially expressed genes were associated with lipid metabolite and flavonoid synthesis. Many genes related to Ca2+ channels, abscisic acid (ABA) signal transduction, and ion transport pathways were upregulated in overexpressing plants. NtLTPI.38 overexpression in salt-stressed tobacco triggered a Ca2+ and K+ influx in leaves, increased the contents of chlorophyll, proline, flavonoids, and osmotic tolerance, and raised enzymatic antioxidant activities as well as the expression level of related genes. However, mutants accumulated more O2− and H2O2, exhibited ionic imbalance, gathered excess Na+, Cl−, and malondialdehyde, with more severe ion leakage. Therefore, NtLTPI.38 enhanced salt tolerance in tobacco by regulating lipid and flavonoid synthesis, antioxidant activity, ion homeostasis, and ABA signaling pathways.

Molecular characterisation of the current high prevalence of the new CPV-2c variants in the Southern Vietnamese dogs signifies a widespread in the worldwide dog population.

Canine parvovirus type 2 (CPV-2) is known as the primary etiological agent cause of acute gastroenteritis, myocarditis and death of canids worldwide. In Vietnam, although CPV-2 infection and its outbreaks are the most important risk factors of the canine's health concern, lack of available information about the molecular epidemiology of the CPV-2. In this study, the complete coding sequences of 10 CPV-2 strains collected from dogs vaccinated with CPV-2 vaccination were analysed to better understand the genomic characteristics of the current circulating CPV-2 in Vietnam. Ten CPV-specific PCR-positive rectal swab samples were collected from dogs with acute symptoms of haemorrhagic diarrhoea and vomiting in Vietnam in 2019. The complete coding sequences of these CPV strains were analysed to determine their phylogeny and genetic relationship with other available CPV strains globally. Analysis of the VP2 gene sequences demonstrated that the studied strains belonged to the new CPV-2c variants with the unique mutations at amino acids 5Ala-Gly and 447Iso-Met . Phylogenetic tree analysis indicated that the studied strains share a common evolutionary origin with the current CPV-2c strains circulating in dogs in Asia countries, including China, Thailand, Taiwan and Mongolia, in recent years. Low sequence identity between the studied strains and commercial vaccine strains was observed. This study provides deep insights into the molecular characteristics, genetic diversity, and evolution of circulating CPV-2 strains in Vietnam. We recommend more studies to estimate the effectiveness of the CPV vaccine and the need to continue developing other effective vaccination essential to better control the widespread of these new CPV-2 variants.

Taste 2 receptors in GtoPdb v.2023.1

Taste 2 receptors or Bitter taste receptors (TAS2Rs) are G protein-coupled receptors expressed in oral sensory cells and a variety of non-gustatory tissues. The ~25 human TAS2Rs share low amino acid sequence identities with other GPCR families and are classified as broadly tuned "generalist" receptors with numerous, chemically diverse bitter agonists, as narrowly tuned "specialist" receptors with very few activators, as intermediately tuned receptors with an average number of agonists, or receptors specialized to interact with chemically defined activators [32]. The number of functional bitter taste receptor genes varies among species and orthologues might not be functionally conserved. Due to their expression in various tissues, the signal transduction of TAS2Rs is complex. Some TAS2Rs interact with drugs such as analgesic, anti-inflammatory, and antibacterial compounds. The specialist database BitterDB contains additional information on bitter compounds and receptors [14].

IgE cross-inhibition between Ara h 1 and Ara h 2 is explained by complex formation of both major peanut allergens

Surprisingly, IgE cross-reactivity between the major peanut allergens Ara h 1, 2, and 3 has been reported despite very low sequence identities. We investigated the unexpected cross-reactivity between peanut major allergens. Cross-contamination of purified natural Ara h 1, 2, 3, and 6 was assessed by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), Western blot test, liquid chromatography-tandem mass spectrometry (LC-MS/MS), and sandwich enzyme-linked immunosorbent assay (ELISA). IgE cross-reactivity was studied with sera of peanut-allergic patients (n= 43) by ELISA and ImmunoCAP inhibition using both intact natural and recombinant allergens and synthetic peptides representing postulated Ara h 1 and Ara h 2 cross-reactive epitopes. Both purified nAra h 1 and nAra h 3 were demonstrated to contain small but significant amounts of Ara h 2 and Ara h 6 (<1%) by sandwich ELISA, SDS-PAGE/Western blot analysis, and LC-MS/MS. IgE cross-inhibition between both 2S albumins and Ara h 1 and Ara h 3 was only observed when using natural purified allergens, not recombinant allergens or synthetic peptides. Apparent cross-reactivity was lost when purified nAra h 1 was pretreated under reducing conditions, suggesting that Ara h 2 and Ara h 6 contaminations may be covalently bound to Ara h 1 via disulfide interactions. True cross-reactivity of both peanut 2S albumins with Ara h 1 and Ara h 3 could not be demonstrated. Instead, cross-contamination with small quantities was shown to be sufficient to cause significant cross-inhibition that can be misinterpreted as molecular cross-reactivity. Diagnostic tests using purified nAra h 1 and nAra h 3 can overestimate their importance as major allergens as a result of the presence of contaminating 2S albumins, making recombinant Ara h 1 and Ara h 3 a preferred alternative.

Isolation and Genotypic Characterization of New Emerging Avian Reovirus Genetic Variants in Egypt

Avian reovirus (ARV) strains cause a variety of symptoms in chickens, including viral arthritis/tenosynovitis, a disease that has emerged as a significant cause of economic losses in commercial chicken flocks in recent years in various countries, including Egypt. Furthermore, ARV strains are frequently isolated from birds suffering from malabsorption. In the actual study, seventy-five samples were collected in 2021 and 2022 from broiler and vaccinated broiler breeder flocks at different farms in Giza Province, Egypt, with reovirus-like symptoms such as significant weight fluctuation and arthritis/malabsorption. ARV was screened using real-time PCR, and fifteen positive samples were detected (20%), which were then subjected to embryonated chicken egg (ECE) isolation and molecular characterization (11/15 sample) of a partial segment of the sigma (σ)C gene (S1-gene). Phylogenetically, nine strains were found to belong to genotypic cluster IV, with 82–89% identity with Israeli ARV 2018, and two strains belong to genotypic cluster V with a 78% nucleotide identity with Japan ARV 2021. No correlation between lesions and genotype was found. The strains under study had a low sequence identity (43–55%) when compared with various commercial vaccines belonging to genotypic cluster I (e.g., strain S1133). These findings imply that novel ARV genotypes representing clusters IV and V have recently been introduced to Egyptian poultry farms. A homologous vaccine is suggested; because this variation raises the possibility that commercial vaccines may not offer protection against circulating ARVs.

Identification and functional characterization of three cytochrome P450 genes for the abietane diterpenoid biosynthesis in Isodon lophanthoides

We identify two ferruginol synthases and a 11-hydroxyferruginol synthase from a traditional Chinese medicinal herb Isodon lophanthoides and propose their involvement in two independent abietane diterpenoids biosynthetic pathways. Isodon lophanthoides is a traditional Chinese medicinal herb rich in highly oxidized abietane-type diterpenoids. These compounds exhibit a wide range of pharmaceutical activities, yet the biosynthesis is barely known. Here, we describe the screening and functional characterization of P450s that oxidize the abietane skeleton abietatriene. We mainly focused on CYP76 family and identified 12 CYP76AHs by mining the RNA-seq data of I. lophanthoides. Among the 12 CYP76AHs, 6 exhibited similar transcriptional expression features as upstream diterpene synthases, including root or leaf-preferential expression pattern and highly MeJA inducibility. These six P450s were considered as first-tier candidates and functionally characterized in yeast and plant cells. In yeast assays showed that both CYP76AH42 and CYP76AH43 were ferruginol synthases hydroxylating the C12 position of abietatriene, whereas CYP76AH46 was characterized as a 11-hydroxyferruginol synthase whichcatalyzes two successive oxidations at C12 and C11 of abietatriene. Heterologous expression of three CYP76AHs in Nicotiana benthamiana resulted in the formation of ferruginol. qPCR analysis showed CYP76AH42 and CYP76AH43 were mainly expressed in the root, which was consistent with the distribution of ferruginol in the root periderms. CYP76AH46 was primarily expressed in the leaves where barely ferruginol or 11-hydroxyferruginol was detected. In addition to distinct organ-specific expression pattern, three CYP76AHs exhibited different genomic structures (w or w/o introns), low protein sequence identities (51-63%) and were placed in separate subclades in the phylogenetic tree. These results suggest that the identified CYP76AHs may be involved in at least two independent abietane biosynthetic pathways in the aerial and underground parts of I. lophanthoides.

Hic Sunt Dracones: Molecular Docking in Uncharted Territories with Structures from AlphaFold2 and RoseTTAfold.

AlphaFold2 and RoseTTAfold impress with their high accuracy in protein structure prediction. However, for structure-based virtual screenings, not only the overall structure but especially the binding sites need to be accurately predicted. In this work, the docking performance for 66 targets with known ligands but without experimental structures available in the protein data bank was elucidated. The results suggest that using an experimental surrogate-ligand complex is often superior over homology models, and only at low sequence identity to the closest homologue AlphaFold2 structures show an equal performance. The generally high fluctuation of receiver operating characteristic area under the curve values obtained for different homology models suggests that multiple combinations of docking programs and homology models should be tested prior to prospective virtual screenings, and in some cases post-processing of crude models might be necessary.

Experimental and computational approaches to characterize a novel amidase that initiates the biodegradation of the herbicide propanil in Bosea sp. P5

The herbicide propanil and its major metabolite 3,4-dichloroaniline (3,4-DCA) are difficult to biodegrade and pose great health and environmental risks. However, studies on the sole or synergistic mineralization of propanil by pure cultured strains are limited. A two-strain consortium (Comamonas sp. SWP-3 and Alicycliphilus sp. PH-34), obtained from a swep-mineralizing enrichment culture that can synergistically mineralize propanil, has been previously reported. Here, another propanil degradation strain, Bosea sp. P5, was successfully isolated from the same enrichment culture. A novel amidase, PsaA, responsible for initial propanil degradation, was identified from strain P5. PsaA shared low sequence identity (24.0–39.7 %) with other biochemically characterized amidases. PsaA exhibited optimal activity at 30 °C and pH 7.5 and had kcat and Km values of 5.7 s-1 and 125 μM, respectively. PsaA could convert the herbicide propanil to 3,4-DCA but exhibited no activity toward other herbicide structural analogs. This catalytic specificity was explained by using propanil and swep as substrates and then analyzed by molecular docking, molecular dynamics simulation and thermodynamic calculations, which revealed that Tyr138 is the key residue that affects the substrate spectrum of PsaA. This is the first propanil amidase with a narrow substrate spectrum identified, providing new insights into the catalytic mechanism of amidase in propanil hydrolysis

A novel thermostable and salt-tolerant carboxylesterase involved in the initial aerobic degradation pathway for pyrethroids in Glycomyces salinus

The long-term and excessive use of pyrethroid pesticides poses substantial health risks and ecosystem concerns. Several bacteria and fungi have been reported that could degrade pyrethroids. The ester-bond hydrolysis using hydrolases is the initial regulatory metabolic reaction of pyrethroids. However, the thoroughly biochemical characterization of hydrolases involved in this process is limited. Here, a novel carboxylesterase, designated as EstGS1 that could hydrolyze pyrethroid pesticides was characterized. EstGS1 showed low sequence identity (<27.03%) compared to other reported pyrethroid hydrolases and belonged to the hydroxynitrile lyase family that preferred short short-chain acyl esters (C2 to C8). EstGS1 displayed the maximal activity of 213.38 U/mg at 60 °C and pH 8.5 using pNPC2 as substrate, with Km and Vmax were 2.21 ± 0.72 mM and 212.90 ± 41.78 µM/min, respectively. EstGS1 is a halotolerant esterase and remains stable in 5.1 M NaCl. Based on molecular docking and mutational analysis, the catalytic triad of S74-D181-H212 and three other substrate-binding residues I108, S159, and G75 are critical for the enzymatic activity of EstGS1. Additionally, 61 and 40 mg/L of deltamethrin and λ-cyhalothrin were hydrolyzed by 20 U of EstGS1 in 4 h. This work presents the first report on a pyrethroid pesticide hydrolase characterized from a halophilic actinobacteria.

Molecular identification and functional analysis of X-linked inhibitor of apoptosis -associated factor-1 (XAF1) in grass carp, Ctenopharyngodon idella

X-linked inhibitor of apoptosis protein (XIAP) -associated factor 1 (XAF1) is an interferon-stimulated gene which exhibits pro-apoptosis effect. In this study, XAF1 was characterized from grass carp Ctenopharyngodon idella and its expression pattern and function were analyzed. The open reading frame (orf) of XAF1 is 789 nucleotides (nt) encoding 262 amino acids. SMART online search results showed that a C2H2-type and six C2HC-type zinc-fingers were found in XAF1, however, the XAF1 of grass carp showed high sequence identity to zebrafish (71%), low sequence identity to tetrapods (21–22%). Rt-qPCR results showed that XAF1 was constitutively expressed in all tested organs/tissues with highest expression in blood. An inductive expression of XAF1 at mRNA level was observed in peripheral blood leucocytes (PBLs) and C. idellus kidney cells (CIKs) after treatment with C. idellus recombinant interferon-γ (rIFNg). Overexpressing XAF1 in CIKs exhibited resistance against grass carp reovirus (GCRV) and more sensitivity to cisplatin. These results implied a functional homologue of XAF1 in evolution, however the mechanism may require further investigation.

Borate Transporters and SLC4 Bicarbonate Transporters Share Key Functional Properties.

Borate transporters are membrane transport proteins that regulate intracellular borate levels. In plants, borate is a micronutrient essential for growth but is toxic in excess, while in yeast, borate is unnecessary for growth and borate export confers tolerance. Borate transporters share structural homology with human bicarbonate transporters in the SLC4 family despite low sequence identity and differences in transported solutes. Here, we characterize the S. cerevisiae borate transporter Bor1p and examine whether key biochemical features of SLC4 transporters extend to borate transporters. We show that borate transporters and SLC4 transporters share multiple properties, including lipid-promoted dimerization, sensitivity to stilbene disulfonate-derived inhibitors, and a requirement for an acidic residue at the solute binding site. We also identify several amino acids critical for Bor1p function and show that disease-causing mutations in human SLC4A1 will eliminate in vivo function when their homologous mutations are introduced in Bor1p. Our data help elucidate mechanistic features of Bor1p and reveal significant functional properties shared between borate transporters and SLC4 transporters.

Sequence adaptation in various environments across organisms; similar yet different superoxide dismutase.

Mycobacterium tuberculosis (M. tuberculosis) is a pathogenic bacterium that attack the host immune system and killed millions of people annually. The body first defense mechanism is through phagocytosis where an outburst of radical oxygen species (ROS) is released and potentially kill the bacterium. M. tuberculosis superoxide dismutase (MtSOD) survives the phagocytosis and suppresses the host immune system by converting reactive oxygen species to H2O2 and O2. SOD proteins are important for cell survival to prevent the buildup of reactive ROS in the cell that may cause damage to DNA, RNA, proteins, and even cell death. MtSOD (oxidized) is homologous to human SOD1 (hSOD1, reduced) which contains both zinc and copper metals along with the disulfide bond. Though, MtSOD lacks the zinc binding site which proposed to play a role in stability, dimerization, and function of SOD. Utilizing Circular Dichroism, an apparent correlation of MtSOD stabilization with the copper and disulfide bond as more stable secondary structure is observed for oxidized apoprotein compared to the reduced apo and oxidized holoprotein. To further investigate the novel role of zinc in SOD (ZnMtSOD), a mutation of zinc binding site consists of T119H, A123H, T192K (modeled from hSOD1) is introduced in the wild type MtSOD. With the low sequence identity of 19.7% but similar in protein structures, the thermodynamic and kinetics results are compared for hSOD1, ZnMtSOD, and wt MtSOD.

Target-template relationships in protein structure prediction and their effect on the accuracy of thermostability calculations.

Improving protein thermostability has been a labor- and time-consuming process in industrial applications of protein engineering. Advances in computational approaches have facilitated the development of more efficient strategies to allow the prioritization of stabilizing mutants. Among these is FEP+, a free energy perturbation implementation that uses a thoroughly tested physics-based method to achieve unparalleled accuracy in predicting changes in protein thermostability. To gauge the applicability of FEP+ to situations where crystal structures are unavailable, here we have applied the FEP+ approach to homology models of 12 different proteins covering 316 mutations. By comparing predictions obtained with homology models to those obtained using crystal structures, we have identified that local rather than global sequence conservation between target and template sequence is a determining factor in the accuracy of predictions. By excluding mutation sites with low local sequence identity (<40%) to a template structure, we have obtained predictions with comparable performance to crystal structures (R2 of 0.67 and 0.63 and an RMSE of 1.20 and 1.16 kcal/mol for crystal structure and homology model predictions, respectively) for identifying stabilizing mutations when incorporating residue scanning into a cascade screening strategy. Additionally, we identify and discuss inherent limitations in sequence alignments and homology modeling protocols that translate into the poor FEP+ performance of a few select examples. Overall, our retrospective study provides detailed guidelines for the application of the FEP+ approach using homology models for protein thermostability predictions, which will greatly extend this approach to studies that were previously limited by structure availability.

Modeling the Orthosteric Binding Site of the G Protein-Coupled Odorant Receptor OR5K1.

With approximately 400 encoding genes in humans, odorant receptors (ORs) are the largest subfamily of class A G protein-coupled receptors (GPCRs). Despite its high relevance and representation, the odorant-GPCRome is structurally poorly characterized: no experimental structures are available, and the low sequence identity of ORs to experimentally solved GPCRs is a significant challenge for their modeling. Moreover, the receptive range of most ORs is unknown. The odorant receptor OR5K1 was recently and comprehensively characterized in terms of cognate agonists. Here, we report two additional agonists and functional data of the most potent compound on two mutants, L1043.32 and L2556.51. Experimental data was used to guide the investigation of the binding modes of OR5K1 ligands into the orthosteric binding site using structural information from AI-driven modeling, as recently released in the AlphaFold Protein Structure Database, and from homology modeling. Induced-fit docking simulations were used to sample the binding site conformational space for ensemble docking. Mutagenesis data guided side chain residue sampling and model selection. We obtained models that could better rationalize the different activity of active (agonist) versus inactive molecules with respect to starting models and also capture differences in activity related to minor structural differences. Therefore, we provide a model refinement protocol that can be applied to model the orthosteric binding site of ORs as well as that of GPCRs with low sequence identity to available templates.